Next Article in Journal
Ecological Risk Assessment of Saltwater Intrusion and Urban Ecosystem Management in Shenzhen City
Next Article in Special Issue
Spatial–Temporal Pattern Evolution and Differentiation Mechanism of Urban Dual Innovation: A Case Study of China’s Three Major Urban Agglomerations
Previous Article in Journal
Ecosystem Health Assessment of the Manas River Basin: Application of the CC-PSR Model Improved by Coupling Coordination Degree
Previous Article in Special Issue
Mapping the Functional Structure of Urban Agglomerations at the Block Level: A New Spatial Classification That Goes beyond Land Use
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Land
Volume 13
Issue 9
10.3390/land13091337
land-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Vegetation Dynamics and Driving Mechanisms Considering Time-Lag and Accumulation Effects: A Case Study of Hubao–Egyu Urban Agglomeration

by
Xi Liu
1,
Guoming Du
1,*,
Xiaodie Zhang
2,
Xing Li
1,*,
Shining Lv
3 and
Yinghao He
1
1
School of Geography and Planning, Sun Yat-sen University, Guangzhou 510006, China
2
Beidou Research Institute, South China Normal University, Foshan 528225, China
3
School of Traffic and Transportation Engineering, Changsha University of Science and Technology, Changsha 410114, China
*
Authors to whom correspondence should be addressed.
Land 2024, 13(9), 1337; https://doi.org/10.3390/land13091337
Submission received: 28 May 2024 / Revised: 17 August 2024 / Accepted: 21 August 2024 / Published: 23 August 2024
(This article belongs to the Special Issue Celebrating the 100th Anniversary of Sun Yat-Sen University in Urban Planning and Land)
Browse Figures
Review Reports Versions Notes

Abstract

:
The Hubao–Egyu Urban Agglomeration (HBEY) was a crucial ecological barrier in northern China. To accurately assess the impact of climate change on vegetation growth, it is essential to consider the effects of time lag and accumulation. In this study, we used a newly proposed kernel Normalized Difference Vegetation Index (kNDVI) as the metric for vegetation condition, and employed partial correlation analysis to ascertain the lag and accumulation period of vegetation response to climate by considering different scenarios (No/Lag/Acc/LagAcc) and various combinations. Moreover, we further modified the traditional residual analysis model. The results are as follows: (1) From 2000 to 2022, the HBEY experienced extensive and persistent greening, with a kNDVI slope of 0.0163/decade. Precipitation was identified as the dominant climatic factor influencing vegetation dynamics. (2) In HBEY, the lag effect of temperature was most distinct, particularly affecting the vegetation in cropland and grassland. The accumulation effect of precipitation was pronounced in grassland. (3) Incorporating lag and accumulation effects into models increases the explanatory power of climate impacts on vegetation dynamics by 6.95% compared to traditional residual models. Our findings hold essential implications for regional ecological regulation and climate change response research.

1. Introduction

Vegetation, a significant element of terrestrial ecosystems, not only provides food and habitat for numerous species through its regulation of material cycles and energy exchanges but also plays an indispensable role in stabilizing climate by sequestering atmospheric carbon dioxide through photosynthesis [1]. However, in recent years, vegetation ecosystems have faced unprecedented challenges due to intensified pressures [2,3]. Global warming and frequent extreme weather events directly impact the growth cycles and rhythm of vegetation [4]. Since entering the “Anthropocene”, the impact of human activities on some regional ecosystems has surpassed that of climate change [5]. Particularly, land-use changes are continually encroaching on vegetation habitats, with this effect still accelerating [6]. Under the combined influence of these factors, the issues are becoming increasingly prominent, posing severe threats to human survival and development [7]. The Hubao–Egyu Urban Agglomeration sits at the junction of the Ordos Plateau and the Loess Plateau, making its vegetation dynamics highly sensitive to climate change. Furthermore, this region has undergone rapid urbanization with frequent land use changes, and it serves as a crucial site for national ecological projects [8]. This makes it an ideal case for examining the human effects on vegetation.
The Normalized Difference Vegetation Index (NDVI) was extensively used for large-scale research due to its excellent performance, long time series, and easy accessibility [9]. However, Huete et al. [10] noted that the NDVI exhibited non-linearities and saturation effects, which could limit its monitoring capabilities in high vegetation coverage areas. Ryu et al. [1] found that NDVI and other Vegetation Indices (VIs) did not reflect photosynthesis itself but only indicate the presence of leaves, introducing some uncertainty in assessing vegetation physiological conditions. Based on these findings, Camps-Valls et al. [11] applied the kernel method to NDVI, creating a more stable and complex phenology-adaptable vegetation index—kNDVI—which had shown excellent performance across various climate zones and latitude zones [12]. The kNDVI has demonstrated superior performance in ecological assessments, including more accurate estimation of vegetation traits, as well as enhanced crop monitoring and forecasting capabilities [13,14,15,16]. Wang et al. has employed kNDVI as an evaluation metric for assessing the dynamics and resilience of vegetation in the Loess Plateau [17]. Additionally, its application in monitoring the arid region’s ecological changes indicated that kNDVI is instrumental in understanding the climate’s effects on vegetation [18,19].
Vegetation change is a complex and prolonged process. The continuous updates of remote sensing data and the extension of meteorological monitoring data have made vegetation–climate response mechanisms a hot topic in ecological remote sensing [20]. To quantitatively analyze these driving mechanisms, researchers often use correlation analysis methods such as Pearson and Spearman methods to quantify the relationships between NDVI and drivers [21]. To further eliminate the mutual influences of meteorological factors, partial correlation coefficients can be calculated [22]. Additionally, a linear regression model [23], geodetector [24], and structural equation model [25] also reveal the pathways and intensities of the effects on vegetation. Currently, the impacts of climate warming, drought, and the fertilization effect of carbon dioxide are constantly being confirmed [26]. Furthermore, studies indicated that plants in the high latitudes and high-altitude mountain areas are more significantly affected by global warming [27]. A study by Nemani et al. [28] on the driving mechanisms in global arid and semi-arid areas indicates that temperature, Vapor Pressure Deficit (VPD), solar radiation, and precipitation all have varying impacts on vegetation growth. Zhang et al. [29] found that potential evapotranspiration (PET) plays a pivotal role in shaping the carbon sequestration capacity of grassland ecosystems in arid regions. It is noteworthy that these studies all indicate that water supply is the primary limiter for vegetation growth in arid areas, especially in grassland ecosystems.
With further research, studies indicated that the response of vegetation to climate change has a certain temporal effect which included a lag and accumulation effect [30]. Due to the lag in vegetation’s response to climate change, previous climate conditions may have a more powerful impact, known as the lag effect [31]. Similarly, the accumulation effect refers to the substantial impact of prolonged and continuous climatic factors on current vegetation [32]. Researchers often quantify the influence of climatic factors on vegetation in different temporal scenarios using correlation or the coefficient of determination from multiple linear regression models [33,34]. Zhao et al. [35] studied the vegetation–climate response of the Loess Plateau, showing that there is a one-month lag between NDVI and precipitation, but no lag effect with temperature. Zhang et al. [36] found that extreme precipitation has a significant impact on at least 98.5% of the Qilian Mountains, with two months’ lag. In addition, some researchers have found the existence of an accumulative effect [37]. Wu et al. [31] believe that incorporating the temporal effect can better predict the real vegetation response and reduce uncertainty. Under rapidly changing climate conditions, the temporal effect is a key factor to be considered. However, the impacts of lag and accumulation effects have always been discussed separately, and the bidirectional coupling of the two has had paid less attention paid to it. At the same time, many researchers have tried to separate human activities from climate change and to judge the respective contributions of vegetation dynamics. At present, researchers mostly use residual analysis to isolate the effects of climatic change on vegetation [38]. However, the attribution usually only considers the climate variables of the corresponding period, and it is unable to deal with the effects of lag and accumulative effects on the model.
This study focuses on the Hubao–Egyu Urban Agglomeration, a typical ecologically fragile area, which experiences large human activity impacts. The kNDVI is employed as a monitoring indicator for vegetation dynamics. The study is structured around three clear objectives: (1) Characterize Vegetation Cover: delineate the distribution patterns and change characteristics of vegetation cover in the study area using trend analysis in the Google Earth Engine; (2) Understand kNDVI Responses: we seek to explore how kNDVI responds to climatic factors by establishing four temporal effects—no temporal effect (No), lag effect (Lag), accumulation effect (Acc), and bidirectional coupling of lag and accumulation (LagAcc). This analysis evaluates the explanatory power of climate change on vegetation dynamics under various scenarios; (3) Refine Model Uncertainties: We seek to refine the uncertainties of the residual model by incorporating temporal effects. This enhancement will allow us to better quantify the impacts of human activity (HA) and climate change (CC) on vegetation dynamics, thereby providing a more robust foundation for our research approach. The findings are crucial for understanding the mechanisms of regional vegetation behavior and maintaining ecosystem balance.

2. Materials and Methods

2.1. Study Area

The Hubao–Egyu Urban Agglomeration (HBEY), situated in Northwest China, comprises the cities of Hohhot, Baotou, Ordos, and Yulin, encompassing 175,000 km2. The respective cities contribute 9.8%, 15.9%, 49.7%, and 24.6% to this total (Figure 1). Characterized by an average elevation of around 1300 m, the terrain predominantly slopes towards the southeast. The region experiences an arid to semi-arid climate, leading to uneven precipitation and a relative scarcity of water resources, which contribute to a sensitive and fragile ecological environment [39]. The heart of the region is the expansive Ordos Plateau, with its northwestern frontier abutting the Kubuqi Desert, China’s seventh-largest desert, often referred to as “China’s last ecological barrier”. This area is subject to intense wind erosion, posing significant challenges to regional ecological conservation [40]. Mu Us Sandy Land has evolved into a “desert oasis” following more than half a century of ecological rehabilitation efforts [41]. In recent years, leveraging its abundant coal resources, the urban agglomeration has undergone rapid economic development. However, this growth has brought with it substantial ecological and environmental pressures, now recognized as a critical strategic challenge for the area.

2.2. Datasets

We used multiple datasets, including kNDVI, the China land cover dataset (CLCD), DEM, meteorological, and administrative division data (Table S1), all of which were converted to the same projection coordinates by ArcGIS 10.8.

2.2.1. kNDVI

kNDVI is a state-of-the-art index produced by nonlinear processing of NDVI using the kernel method. Specifically, we defined NDVI in Hilbert space using a Radial Basis Function (RBF) reproducing kernel, which generated the kernel NDVI. The NDVIs were sourced from the MOD13Q1, with 16 days and 250 m resolution from 2000 to 2022 (accessed on 13 January 2024). Due to its high resolution and lower error incidence, it is widely used in mid- to high-latitude areas [42]. To ensure the reliability and mitigate the effects of clouds, aerosols, and solar zenith angles, we used a maximum value composite method to preprocess the raw data and generate a monthly NDVI dataset. Furthermore, we used Savitzky–Golay filter to smooth the data, effectively reducing noise [43]. The kNDVI have the ability to adapt effectively to areas with either sparse or dense vegetation coverage, and have exhibited robust performance in other applications [44].

2.2.2. Meteorological Data

Temperature (T), precipitation (P), and potential evapotranspiration (PET) data were sourced from the National Tibetan Plateau Data Center (TPDC, accessed on 27 January 2024), encompassing monthly records from 2000 to 2022 [45]. The potential evapotranspiration of this product was calculated using the Hargreaves formula based on monthly average, minimum, and maximum temperature data. This product has been extensively validated by various studies, and its credibility has been confirmed for assessing climate change impacts across different regions in China [46,47].

2.2.3. LULC Data

The land use data were sourced from the 30 m resolution annual China Land Cover Dataset (CLCD, accessed on 27 January 2024) [48]. The dataset includes nine land use types. In the Hubao–Egyu region, eight categories are presented. Given the negligible pixel number of Shrub and Wetland categories, these were overlooked, to streamline the data interpretation.

2.3. Methods

Our research methodology is illustrated in Figure 2. Specifically, we first constructed a dataset of kNDVI and driving factors on the GEE [49]. We then analyzed the trend and significance of kNDVI using Theil–Sen analysis and Mann–Kendall test. Next, we calculate the second-order partial correlation coefficients between kNDVI and climatic factors under different temporal scenarios to determine the lag and accumulation periods. Finally, we refine the traditional residual model by considering the temporal effects, and analyze the contributions under various temporal scenarios.

2.3.1. Theil–Sen Trend and Mann–Kendall Test

This study used the Theil–Sen trend analysis and the Mann–Kendall test to analyze the change characteristics of kNDVI at the pixel level. The Theil–Sen trend was used to calculate the trends of the kNDVI [50]. The Mann–Kendall test can detect abrupt changes while assessing the significance of trends and does not require data to adhere to a specific distribution, such as the normal distribution [51]. This test utilized the statistic Z to assess the significance [52]. It is frequently used in conjunction with the Theil–Sen trend as a basis for analyzing long-term series trends (Table S2).

2.3.2. Evaluating Time-Lag and Accumulation Effects

This study incorporated T, P, and PET as drivers to quantify their effects on kNDVI. We used partial correlation analysis to accurately assess the impact of each variable on kNDVI while eliminating interference from others [53]. Additionally, it was used to evaluate the influence of time-lag and accumulation effects. Based on this, our study calculated the second-order partial correlation coefficients between kNDVI and meteorological factors under different temporal effects. We then compared a total of ten sets of partial correlation coefficients to identify the combination with the strongest correlation, which represents the lag-accumulation period (Figure 3).

2.3.3. Modified Residual Model

We propose a modified model that allows for the examination of the impacts of climate change and human activity on vegetation under varying temporal effects. To verify the performance of our modified residual model, we conducted a systematic assessment using three key performance metrics: R2, Root Mean Square Error (RMSE), and Mean Absolute Error (MAE) [54]. R2 represents the proportion of the dependent variable that is accounted for by the model, RMSE gives an estimate of the error magnitude, and MAE offers an average of the absolute errors. To further evaluate the impacts of CC and HA on vegetation, we delineate interaction scenarios and calculate contribution rates under four temporal effects according to the rules specified in Table S3 [55].

3. Results

3.1. Spatiotemporal Variations of Vegetation Dynamic

The kNDVI in HBEY was lower in the west and higher in the east. High coverage was primarily found in the eastern part of Yulin and the northern part of Hohhot, aligning with the geomorphological distribution of the western highland barren areas and the eastern hilly grassland regions of the study area (Supplementary Figure S1). HBEY’s average kNDVI ranged between 0.030 and 0.072, with the minimum and maximum values occurring in 2001 and 2018, respectively. Typically, the values clustered near the lower quartile, suggesting that the vegetation coverage consistently remained at a relatively low level (Supplementary Figure S2). Throughout the study period, kNDVI exhibited an upward trend (slope = 0.0163/decades, p < 0.005, R2 = 0.85), with minor interannual fluctuations. It is worth mentioning that after 2011, fluctuations in kNDVI noticeably increased, and the improvement trend slowed, signaling a need for enhanced continuous monitoring of the restoration areas to consolidate ecological recovery gains.
The vegetation dynamic in HBEY was illustrated in Figure 4. The results indicate that 70.96% of the area exhibited significant change, of which nearly all (99.49%) were positive, demonstrating overall greening. It revealed distinct differences in trend across the cities, with the order of improvement areas (SLI and SII) in descending order as follows: Yulin (97.44%) > Ordos (83.63%) > HBEY (78.87%) > Hohhot (71.80%) > Baotou (39.31%). Yulin showed an extremely high proportion of significant improvements (95.76%), with almost no degradation. The Mu Us Sandy Land has also continued greening over the past two decades, reflecting the substantial efforts in ecological protection in this area. Baotou and Hohhot had relatively high proportions of stable areas (59.12% and 26.73%), with scattered regions of degradation sporadically in the central, northern, and southern fringes. Although these proportions are small, they still warrant attention. Furthermore, there were areas where kNDVI showed a declining trend (Slope < 0), although only 10.24%, often concentrated in urban centers, along the Yellow River, and in some heavy industrial areas. This indicates that excessive land development and urban expansion are still significant constraints on vegetation growth.

3.2. Time-Lag and Accumulation Effects of Climate Factors on kNDVI

The results indicated that kNDVI was generally positively correlated with temperature, precipitation, and PET across the study area (Figure 5a–c). Notably, the partial correlation coefficient with precipitation was obviously higher than with the other factors, suggesting that precipitation was the dominant climatic driver of kNDVI changes in HBEY. Specifically, the central and eastern parts of HBEY showed strong positive correlations with precipitation, indicating a substantial influence of water availability on vegetation growth. Conversely, kNDVI exhibited a relatively weak correlation with temperature overall, with a noticeable negative correlation particularly concentrated in the northern Baotou and the western Ordos Plateau. In these areas, a negative correlation between kNDVI and PET was also observed. Additionally, we examined the vegetation response under nine different temporal effect scenarios (Supplementary Figures S3–S5). When incorporating the temporal effect, the correlation coefficients of kNDVI and climate factors changed to different degrees. Specifically, the correlation changes between PET and kNDVI were minimal, whereas the correlation with precipitation showed a marked decrease, and the correlation with temperature exhibited an increase. We also found that when the lag and accumulation time were sufficiently extended, the correlation may reverse, often occurring in areas with weak initial correlations. For example, in the southeastern region of the study area, what was initially a positive correlation between precipitation and kNDVI shifted to a negative correlation in scenarios with extended lags. Similarly, in eastern Yulin, an initially negative correlation with temperature transitioned to a positive correlation when the accumulation period was prolonged.
The strength and direction of the climatic factors driving vegetation change—including both direct and lag/accumulation effects—vary greatly across regions. Generally, the lag effects of temperature and PET were more pronounced, while the accumulation effects of precipitation were stronger. Specifically, for lag and accumulation effects, temperature, PET, and precipitation exhibited lag periods of 1.708 ± 0.98, 1.16 ± 0.79, and 0.155 ± 0.33 months, and accumulation periods of 0.565 ± 0.63, 0.282 ± 0.31, and 0.462 ± 0.42 months, respectively. According to Figure 5d–f, the temporal effects of temperature on kNDVI were predominantly characterized by the scenario of Lag2 and Acc0 (41.13%), Lag2 and Acc1 (30.64%), Lag1 and Acc0 (13.2%), and Lag1 and Acc2 (12.54%). The response of vegetation to temperature was primarily dominated by the lag effects and the bidirectional coupling effect (LagAcc), accounting for 54.79% and 43.33% (Supplementary Table S4), respectively. The bidirectional coupling effects mainly affected the eastern part, while the lag effects were dominant across most of the Ordos Plateau and Baotou. Among the different land use types (Supplementary Figure S6b), almost none show accumulation effects solely. Nearly half of the forest areas exhibited no temporal effects. Cropland exhibited the most pronounced LagAcc effects, at 70.65%.
In terms of precipitation and PET, the predominant temporal effects were: Lag0 and Acc0 (53.55%), Lag0 and Acc1 (38.85%), and Lag2 and Acc1 (6.53%) for precipitation, and Lag1 and Acc0 (49.76%), Lag1 and Acc1 (23.7%), and Lag2 and Acc0 (18.03%) for PET. In terms of precipitation, no temporal effect (53.55%) and acc effect (38.87%) had the greatest impact, predominantly affecting the central and peripheral areas. From Supplementary Figure S6, the vast majority of forest vegetation showed no temporal response to precipitation (88.03%), while grasslands exhibited a relatively pronounced accumulation effect (42.89%). Conversely, for PET, the lag effect (69.65%) and LagAcc effect (24.15%) were the dominant temporal effects. The accumulation effect of forest was particularly prominent (46.14%), while cropland and grassland showed the most pronounced lag effects (Supplementary Figure S6 and Supplementary Table S5). This phenomenon highlights the complex interactions and adaptive mechanisms between vegetation and different climatic factors, underscoring adaptability as a key to vegetation survival in varied climatic environments.

3.3. Contributions of Human Activity and Climate Change to Vegetation Dynamics

According to Figure 6, the vegetation changes in HBEY were mainly influenced by both HA and CC (greening: 68.3%, browning: 9.21%). Among these were regions where vegetation browning was primarily distributed in urban built-up areas and along the Yellow River, while the most extensive greening driven by the two is concentrated in the southeastern and northeastern parts of HBEY. Vegetation greening driven by human activity accounts for 17.32%, mainly distributed in the Ordos Plateau and the Kubuqi Desert. In contrast, greening driven by climate change accounts for only 4.15% and is sporadically distributed in the northern and southwestern corners of the study area. Comparatively, human activity exerted a stronger influence on vegetation than climate change. The extensive greening experienced by the urban agglomeration is mainly influenced by both human activity and climate change. As shown in Figure 6, in built-up urban areas (A), the significant traces of urban expansion reveal the direct and drastic impact of HA on vegetation degradation. In the Jizi Bay (B), floodplain sedimentation caused by CC is the main reason for vegetation browning [56]. The case of area (C) further demonstrates this complexity. On the one hand, humans actively improve vegetation by reclaiming barren land ③; on the other hand, climate change affects vegetation through natural processes such as grassland succession ④. In area (D), the joint action of human activity and climate change is fully reflected in sand control, demonstrating their synergistic effects on vegetation improvement [57]. In addition, the vegetation degradation caused by the development of artificial lakes ⑥ and the vegetation improvement brought by industrial park ecological management ⑦ in area (E) coexist, highlighting the duality of human activities and climate change in influencing vegetation. More importantly, our findings also demonstrate that the modified residual model can maintain objectively true results within individual pixels, indirectly demonstrating the importance of time effects in attributing vegetation changes, and of adaptability as a key to vegetation survival in varied climatic environments.
Under the consideration of four temporal effects—No, Lag, Acc, and LagAcc—the explanatory power of CC on vegetation shows an increasing trend (Figure 7), with values of 40.90% ± 21.32%, 46.91% ± 17.56%, 43.03% ± 21.03%, and 47.85% ± 18.58%, respectively. It confirmed that considering the lag effect rather than the accumulation effect can more strongly reveal the contribution of climate change. Compared to the No temporal effect, the LagAcc effect shows a relative increase of 6.95% (R2 increase of 0.08), particularly noticeable within the Ordos Plateau. This finding implied that the contribution of CC is somewhat obscured by temporal effects. By accounting for lag and accumulation effects, we can effectively restore the true contribution of CC and reduce the uncertainty in the climate–vegetation response. Under the LagAcc effect, the regions where climate change dominate were mainly located in the northern part of Baotou, Hohhot City, and the eastern area of Yulin. According to the results presented in Figure 7e,f, human activity consistently plays a dominant role on vegetation dynamics in HBEY. In regions with harsh natural environments, as well as in highly urbanized areas, a significant human contribution is observed. However, it is important to note that these contributions are not always positive, underscoring the need for a thorough understanding of the mechanisms driving the impact of CC and HA on vegetation dynamics. This finding is important for making informed, sustainable decisions within this complex ecosystem.
To further assess our modified model’s effectiveness, we calculated its performance metrics (R2, MAE, RMSE) using ten-fold cross-validation under four scenarios [30]. The results indicated that incorporating lag and accumulation effects enhanced model performance (Figure 8). Specifically, the introduction of lag effects, by capturing the delayed response of vegetation to climate change, greatly improved the model’s explanatory power and predictive accuracy. The inclusion of accumulation effects further refined the model’s understanding of vegetation dynamics. Additionally, the model that considered both lag and accumulation effects (LagAcc) exceled across all performance metrics. This suggests that it can most effectively capture the complexity of vegetation dynamics, fully demonstrating the effectiveness of the modified residual model.

4. Discussion

4.1. The Temporal Effect Mechanism of Climate Factors on Vegetation

The existence of temporal effects was attributed to the specific rhythm of vegetation growth and development, which requires a certain cycle to occur [31]. When climate changes, vegetation needs time to adapt and exhibit a corresponding response. The adaptive mechanisms of vegetation physiology can also mitigate the impacts of temporary climate changes [58,59]. Wen’s research indicated that, on a global scale, the accumulation effects of vegetation and temperature were dominant in low-latitude regions, while in arid and semi-arid areas, the lag effects of precipitation prevail [60]. Our findings showed that the temporal effect of temperature in HBEY is predominantly characterized by lag effects, followed by accumulation effects. This seems contradictory to the conclusion of Wen [60] and Ding [61], who claimed that vegetation response to temperature was mostly accumulative or instantaneous in Central and North Asia. However, it was worth noting that Liu’s research suggested that temporal effects vary with different climates and vegetation types. Furthermore, Kong found that in arid regions, grasslands exhibited a more pronounced lag effect to temperature compared to other vegetation types [62]. This is in line with our findings, given that HBEY primarily consists of grassland and cropland, both of which exhibited longer lag periods (Tables S5 and S6). The root system of grassland can store nutrients to regulate growth imbalances caused by temperature fluctuations, which aligns with previous research [63]. Ma [30] found that vegetation in high-altitude areas required a longer time to accumulate heat to support its growth and development in relatively cold environments, which also corroborated the accumulation effects observed on the Ordos Plateau. Additionally, Wan [64] found that in the northern temperate regions, the large seasonal temperature variations lead to stronger resistance of vegetation to temperature fluctuations. Therefore, combined with previous studies, we believe that in temperature-limited areas, such as high-altitude and high-latitude regions, vegetation responses to temperature tends to have relatively longer lag times.
The temporal effect of PET is generally consistent with temperature but relatively weaker in magnitude, which was basically consistent with the studies carried on the arid areas [60]. Zhao’s study on the temporal effects of arid grassland vegetation also found similar phenomena [35]. In semi-arid areas, an increase in PET implied more severe soil moisture depletion. Ershadi [65] suggests that higher PET requires vegetation to accumulate more precipitation to maintain the water necessary for photosynthesis. Without additional water supply, this could further exacerbate water stress, requiring time for vegetation to adjust its physiological strategies to cope with the increased challenges [66]. Precipitation’s temporal effect on vegetation was primarily characterized by an instantaneous response, followed by accumulation effect. Anderegg [67] found that vegetation in semi-arid region exhibited heightened sensitivity to precipitation changes, coupled with a comparatively weak resilience to water scarcity. From this perspective, measures such as artificial irrigation and improving water availability were effective interventions to influence vegetation growth conditions. Additionally, a previous study by Paudel and Andersen [68] suggested that vegetation does not directly respond to precipitation but rather to actual soil moisture. The soil’s memory of water was long-term and accumulative: when rainfall is sufficient, the soil can store excess water to meet the constant growth requirements. Moreover, Braswell [69] and Zhang [70] argued that the process of soil moisture being transferred to the surface also requires time. In comparison, root systems in grasslands can effectively absorb soil moisture, but they also require more accumulated water to bring about significant changes in vegetation [71].

4.2. The Necessity and Effectiveness of the Modified Residual Model

Relying solely on immediate climate data may not fully explain the growth patterns of vegetation in arid areas [61]. However, traditional residual models mainly consider the influence of independent variables within the same period. This approach does not fully capture the lag and accumulation effects of vegetation’s response to climate change, potentially leading to errors in vegetation attribution analysis [31,72,73]. Ding’s study indicated that accounting for lag effects increases the explanatory power of residual models by 9% during the growing season [61]. Similarly, Ma [34] discovered that extreme climate factors also exhibited temporal effects, and by establishing linear models between vegetation indices (VIs) and climate variables under different lag and accumulation scenarios, the explanatory power of residual models improved by approximately 8%. However, traditional linear relationships may involve interactions between seasonal variations and climate change, leading to uncertainty [74]. Our proposed model addresses this by using partial correlation analysis to eliminate the interdependence among climate variables, thereby enhancing reliability. The results confirmed that scenarios incorporating both lag and accumulation effects are the most effective, enhancing the explanatory power of climate factors on vegetation dynamics, particularly in barren and grassland areas (Supplementary Table S4). Moreover, Wu [31] found that when considering lag effects, the contribution of climate change to global vegetation dynamics increased from 58% to 64%, and when considering cumulative effects, it increased to 68%. Ma [30] integrated both lag and accumulation effects and discovered that the contribution of temperature and precipitation on vegetation in the Qilian Mountains rose from 35.6% to 54.1%. Those were largely consistent with our findings. Therefore, we claimed that climate variables that consider temporal effects more accurately reflect the actual contribution of climate change, effectively reducing the interference and errors caused by non-climatic factors. This allowed for a more precise isolation of the independent contribution of climate change to vegetation dynamics.

4.3. Uncertainty and Future Outlook

Considering temporal effects, we propose a second-order lag accumulation partial correlation analysis and modify the residual analysis model. It quantifies vegetation response characteristics in four scenarios and ten sets, providing new insights into the interaction among climate, vegetation, and human activity. However, there are limitations. Firstly, the driving factors of vegetation change are complex and difficult to fully enumerate. This study attributes vegetation changes in the HBEY to three influential driving factors in arid to semi-arid regions: temperature, precipitation, and potential evapotranspiration. The effects of solar radiation, CO2 fertilization, nitrogen deposition, and extreme weather were not considered [75,76,77]. Moreover, human activities such as ecological engineering and land abandonment, as well as vegetation pests cannot be quantified at the pixel scale [78,79]. Additionally, the driving mechanisms are extremely complex and cannot be fully reflected by simple statistical models. Future research may consider using the SEM model to identify the interaction effects [25]. This could further deepen our understanding of vegetation driving mechanisms.

5. Final Remarks

Our study explored the spatiotemporal dynamics of kNDVI in the Hubao–Egyu Urban Agglomeration from 2000 to 2022, employing second-order partial correlation analysis and a modified residual model to more accurately quantify the contributions of climate change and human activity to vegetation dynamics. Over the study period, notable vegetation increases were observed in HBEY, particularly in Yulin City. Precipitation emerged as the primary climatic factor driving these changes, with early climate impacts proving more significance due to lag effects. Considering temporal effects, the modified residual model restored the contribution of climate change to vegetation dynamics (increasing explanatory power by 6.95%), demonstrating the necessity and effectiveness of modifying temporal effects in residual models. Vegetation dynamics in HBEY are primarily influenced by both human activity and climate change, with human activities playing a dominant role in greening, especially in urban areas. This approach not only offers new insights into the interplay among climate, vegetation, and human activity in arid to semi-arid regions, but also holds profound implications for ecological conservation and resource management in these areas. Future research will further explore the underlying mechanisms of vegetation dynamics to promote the sustainable development of ecosystems.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/land13091337/s1, Figures S1–S7, Tables S1–S6.

Author Contributions

Conceptualization, G.D.; Data curation, S.L. and Y.H.; Formal analysis, X.L. (Xi Liu); Funding acquisition, G.D. and X.L. (Xing Li); Investigation, Y.H.; Methodology, X.L. (Xi Liu); Software, X.L. (Xi Liu), X.Z. and S.L.; Supervision, G.D. and X.L. (Xing Li); Visualization, X.Z.; Writing—original draft, X.L. (Xi Liu); Writing—review and editing, X.L. (Xi Liu), X.Z. and X.L. (Xing Li). All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Program on Key Research Projects of China [grant numbers 2017YFC1502706] and the starting grant for introduced talents from Sun Yat-sen University (to X.L., 37000-12240012).

Data Availability Statement

The data sources for this paper can be found in Table S1. Other data that support the findings of this study are available from the author, upon reasonable request.

Acknowledgments

We appreciate anonymous reviewers for their constructive comments on the manuscript.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Ryu, Y.; Berry, J.A.; Baldocchi, D.D. What is global photosynthesis? History, uncertainties and opportunities. Remote Sens. Environ. 2019, 223, 95–114. [Google Scholar] [CrossRef]
  2. Seddon, A.W.R.; Macias-Fauria, M.; Long, P.R.; Benz, D.; Willis, K.J. Sensitivity of global terrestrial ecosystems to climate variability. Nature 2016, 531, 229–232. [Google Scholar] [CrossRef]
  3. Li, X.; Ryu, Y.; Xiao, J.; Dechant, B.; Liu, J.; Li, B.; Jeong, S.; Gentine, P. New-generation geostationary satellite reveals widespread midday depression in dryland photosynthesis during 2020 western U.S. heatwave. Sci. Adv. 2023, 9, eadi0775. [Google Scholar] [CrossRef] [PubMed]
  4. Bonan, G.B. Forests and climate change: Forcings, feedbacks, and the climate benefits of forests. Science 2008, 320, 1444–1449. [Google Scholar] [CrossRef] [PubMed]
  5. MacDougall, A.S.; McCann, K.S.; Gellner, G.; Turkington, R. Diversity loss with persistent human disturbance increases vulnerability to ecosystem collapse. Nature 2013, 494, 86–89. [Google Scholar] [CrossRef]
  6. Bala, G.; Caldeira, K.; Wickett, M.; Phillips, T.J.; Lobell, D.B.; Delire, C.; Mirin, A. Combined climate and carbon-cycle effects of large-scale deforestation. Proc. Natl. Acad. Sci. USA 2007, 104, 6550–6555, Corrected in Proc. Natl. Acad. Sci. USA 2007, 104, 9911. [Google Scholar] [CrossRef]
  7. Chen, Y.; Duo, L.; Zhao, D.; Zeng, Y.; Guo, X. The response of ecosystem vulnerability to climate change and human activities in the Poyang lake city group, China. Environ. Res. 2023, 233, 116473. [Google Scholar] [CrossRef]
  8. Qi, T.; Ren, Q.; Zhang, D.; Lu, W.; He, C. Impacts of urban expansion on vegetation in drylands: A multiscale analysis based on the vegetation disturbance index. Ecol. Indic. 2023, 147, 109984. [Google Scholar] [CrossRef]
  9. Hilker, T.; Lyapustin, A.I.; Hall, F.G.; Myneni, R.; Knyazikhin, Y.; Wang, Y.; Tucker, C.J.; Sellers, P.J. On the measurability of change in Amazon vegetation from MODIS. Remote Sens. Environ. 2015, 166, 233–242. [Google Scholar] [CrossRef]
  10. Huete, A.; Didan, K.; Miura, T.; Rodriguez, E.P.; Gao, X.; Ferreira, L.G. Overview of the radiometric and biophysical performance of the MODIS vegetation indices. Remote Sens. Environ. 2002, 83, 195–213. [Google Scholar] [CrossRef]
  11. Camps-Valls, G.; Campos-Taberner, M.; Moreno-Martínez, A.; Walther, S.; Duveiller, G.; Cescatti, A.; Mahecha, M.D.; Muñoz-Marí, J.; García-Haro, F.J.; Guanter, L.; et al. A unified vegetation index for quantifying the terrestrial biosphere. Sci. Adv. 2021, 7, eabc7447. [Google Scholar] [CrossRef]
  12. Wang, Q.; Moreno-Martínez, Á.; Muñoz-Marí, J.; Campos-Taberner, M.; Camps-Valls, G. Estimation of vegetation traits with kernel NDVI. ISPRS J. Photogramm. Remote Sens. 2023, 195, 408–417. [Google Scholar] [CrossRef]
  13. Gu, Z.; Chen, X.; Ruan, W.; Zheng, M.; Gen, K.; Li, X.; Deng, H.; Chen, Y.; Liu, M. Quantifying the direct and indirect effects of terrain, climate and human activity on the spatial pattern of kNDVI-based vegetation growth: A case study from the Minjiang River Basin, Southeast China. Ecol. Inform. 2024, 80, 102493. [Google Scholar] [CrossRef]
  14. Yu, H.; Yang, Q.; Jiang, S.; Zhan, B.; Zhan, C. Detection and Attribution of Vegetation Dynamics in the Yellow River Basin Based on Long-Term Kernel NDVI Data. Remote Sens. 2024, 16, 1280. [Google Scholar] [CrossRef]
  15. Amin, E.; Pipia, L.; Belda, S.; Perich, G.; Graf, L.V.; Aasen, H.; Van Wittenberghe, S.; Moreno, J.; Verrelst, J. In-season forecasting of within-field grain yield from Sentinel-2 time series data. Int. J. Appl. Earth Obs. Geoinf. 2024, 126, 103636. [Google Scholar] [CrossRef]
  16. Feng, X.; Tian, J.; Wang, Y.; Wu, J.; Liu, J.; Ya, Q.; Li, Z. Spatio-Temporal Variation and Climatic Driving Factors of Vegetation Coverage in the Yellow River Basin from 2001 to 2020 Based on kNDVI. Forests 2023, 14, 620. [Google Scholar] [CrossRef]
  17. Wang, Z.; Fu, B.; Wu, X.; Li, Y.; Feng, Y.; Wang, S.; Wei, F.; Zhang, L. Vegetation resilience does not increase consistently with greening in China’s Loess Plateau. Commun. Earth Environ. 2023, 4, 1–10. [Google Scholar] [CrossRef]
  18. Wu, J.; Sun, Z.; Yao, Y.; Liu, Y. Trends of Grassland Resilience under Climate Change and Human Activities on the Mongolian Plateau. Remote Sens. 2023, 15, 2984. [Google Scholar] [CrossRef]
  19. Chen, Z.; Zhang, X.; Jiao, Y.; Cheng, Y.; Zhu, Z.; Wang, S.; Zhang, H. Investigating the spatio-temporal pattern evolution characteristics of vegetation change in Shendong coal mining area based on kNDVI and intensity analysis. Front. Ecol. Evol. 2023, 11, 1344664. [Google Scholar] [CrossRef]
  20. Sun, Z.; Wang, X.; Yamamoto, H.; Tani, H.; Zhong, G.; Yin, S.; Guo, E. Spatial pattern of GPP variations in terrestrial ecosystems and its drivers: Climatic factors, CO2 concentration and land-cover change, 1982–2015. Ecol. Inform. 2018, 46, 156–165. [Google Scholar] [CrossRef]
  21. Wang, J.; Price, K.P.; Rich, P.M. Spatial patterns of NDVI in response to precipitation and temperature in the central Great Plains. Int. J. Remote Sens. 2001, 22, 3827–3844. [Google Scholar] [CrossRef]
  22. Xu, H.-J.; Wang, X.-P.; Yang, T.-B. Trend shifts in satellite-derived vegetation growth in Central Eurasia, 1982–2013. Sci. Total Environ. 2017, 579, 1658–1674. [Google Scholar] [CrossRef] [PubMed]
  23. Chen, B.; Zhang, X.; Tao, J.; Wu, J.; Wang, J.; Shi, P.; Zhang, Y.-J.; Yu, C. The impact of climate change and anthropogenic activities on alpine grassland over the Qinghai-Tibet Plateau. Agric. For. Meteorol. 2014, 189–190, 11–18. [Google Scholar] [CrossRef]
  24. Peng, W.; Kuang, T.; Tao, S. Quantifying influences of natural factors on vegetation NDVI changes based on geographical detector in Sichuan, western China. J. Clean. Prod. 2019, 233, 353–367. [Google Scholar] [CrossRef]
  25. Yang, L.; Shen, F.; Zhang, L.; Cai, Y.; Yi, F.; Zhou, C. Quantifying influences of natural and anthropogenic factors on vegetation changes using structural equation modeling: A case study in Jiangsu Province, China. J. Clean. Prod. 2021, 280, 124330. [Google Scholar] [CrossRef]
  26. Yan, Y.; Piao, S.; Hammond, W.M.; Chen, A.; Hong, S.; Xu, H.; Munson, S.M.; Myneni, R.B.; Allen, C.D. Climate-induced tree-mortality pulses are obscured by broad-scale and long-term greening. Nat. Ecol. Evol. 2024, 8, 912–923. [Google Scholar] [CrossRef]
  27. A Walker, D.; E Epstein, H.; Raynolds, M.K.; Kuss, P.; A Kopecky, M.; Frost, G.V.; A Daniëls, F.J.; O Leibman, M.; Moskalenko, N.G.; Matyshak, G.V.; et al. Environment, vegetation and greenness (NDVI) along the North America and Eurasia Arctic transects. Environ. Res. Lett. 2012, 7, 015504. [Google Scholar] [CrossRef]
  28. Nemani, R.R.; Keeling, C.D.; Hashimoto, H.; Jolly, W.M.; Piper, S.C.; Tucker, C.J.; Myneni, R.B.; Running, S.W. Climate-Driven Increases in Global Terrestrial Net Primary Production from 1982 to 1999. Science 2003, 300, 1560–1563. [Google Scholar] [CrossRef]
  29. Zhang, X.; Chen, Y.; Zhang, Q.; Xia, Z.; Hao, H.; Xia, Q. Potential evapotranspiration determines changes in the carbon sequestration capacity of forest and grass ecosystems in Xinjiang, Northwest China. Glob. Ecol. Conserv. 2023, 48, e02737. [Google Scholar] [CrossRef]
  30. Ma, Y.; Guan, Q.; Sun, Y.; Zhang, J.; Yang, L.; Yang, E.; Li, H.; Du, Q. Three-dimensional dynamic characteristics of vegetation and its response to climatic factors in the Qilian Mountains. Catena 2022, 208, 105694. [Google Scholar] [CrossRef]
  31. Wu, D.H.; Zhao, X.; Liang, S.L.; Zhou, T.; Huang, K.C.; Tang, B.J.; Zhao, W.Q. Time-lag effects of global vegetation responses to climate change. Glob. Chang. Biol. 2015, 21, 3520–3531. [Google Scholar] [CrossRef] [PubMed]
  32. Guo, L.; Wang, J.; Li, M.; Liu, L.; Xu, J.; Cheng, J.; Gang, C.; Yu, Q.; Chen, J.; Peng, C.; et al. Distribution margins as natural laboratories to infer species’ flowering responses to climate warming and implications for frost risk. Agric. For. Meteorol. 2019, 268, 299–307. [Google Scholar] [CrossRef]
  33. Ivits, E.; Horion, S.; Fensholt, R.; Cherlet, M. Drought footprint on European ecosystems between 1999 and 2010 assessed by remotely sensed vegetation phenology and productivity. Glob. Chang. Biol. 2014, 20, 581–593. [Google Scholar] [CrossRef] [PubMed]
  34. Ma, M.; Wang, Q.; Liu, R.; Zhao, Y.; Zhang, D. Effects of climate change and human activities on vegetation coverage change in northern China considering extreme climate and time-lag and -accumulation effects. Sci. Total Environ. 2023, 860, 160527. [Google Scholar] [CrossRef]
  35. Zhao, J.; Huang, S.; Huang, Q.; Wang, H.; Leng, G.; Fang, W. Time-lagged response of vegetation dynamics to climatic and teleconnection factors. Catena 2020, 189, 104474. [Google Scholar] [CrossRef]
  36. Zhang, S.; Zhang, J.; Liang, S.; Liu, S.; Zhou, Y. A perception of the nexus “resistance, recovery, resilience” of vegetations responded to extreme precipitation pulses in arid and semi-arid regions: A case study of the Qilian Mountains Nature Reserve, China. Sci. Total Environ. 2022, 843, 157105. [Google Scholar] [CrossRef] [PubMed]
  37. Liu, L.; Guan, J.; Zheng, J.; Wang, Y.; Han, W.; Liu, Y. Cumulative effects of drought have an impact on net primary productivity stability in Central Asian grasslands. J. Environ. Manag. 2023, 344, 118734. [Google Scholar] [CrossRef]
  38. Wang, J.; Wang, K.; Zhang, M.; Zhang, C. Impacts of climate change and human activities on vegetation cover in hilly southern China. Ecol. Eng. 2015, 81, 451–461. [Google Scholar] [CrossRef]
  39. Zhang, L.; Fang, C.; Zhao, R.; Zhu, C.; Guan, J. Spatial–temporal evolution and driving force analysis of eco-quality in urban agglomerations in China. Sci. Total Environ. 2023, 866, 161465. [Google Scholar] [CrossRef]
  40. Ren, M.; Chen, W.; Wang, H. Ecological Policies Dominated the Ecological Restoration over the Core Regions of Kubuqi Desert in Recent Decades. Remote Sens. 2022, 14, 5243. [Google Scholar] [CrossRef]
  41. Gao, W.; Zheng, C.; Liu, X.; Lu, Y.; Chen, Y.; Wei, Y.; Ma, Y. NDVI-based vegetation dynamics and their responses to climate change and human activities from 1982 to 2020: A case study in the Mu Us Sandy Land, China. Ecol. Indic. 2022, 137, 108745. [Google Scholar] [CrossRef]
  42. Cunha, A.; Alvalá, R.; Nobre, C.; Carvalho, M. Monitoring vegetative drought dynamics in the Brazilian semiarid region. Agric. For. Meteorol. 2015, 214–215, 494–505. [Google Scholar] [CrossRef]
  43. Chen, Y.; Cao, R.; Chen, J.; Liu, L.; Matsushita, B. A practical approach to reconstruct high-quality Landsat NDVI time-series data by gap filling and the Savitzky–Golay filter. ISPRS J. Photogramm. Remote Sens. 2021, 180, 174–190. [Google Scholar] [CrossRef]
  44. Forzieri, G.; Dakos, V.; McDowell, N.G.; Ramdane, A.; Cescatti, A. Emerging signals of declining forest resilience under climate change. Nature 2022, 608, 534–539. [Google Scholar] [CrossRef]
  45. Peng, S.; Ding, Y.; Liu, W.; Li, Z. 1 km monthly temperature and precipitation dataset for China from 1901 to 2017. Earth Syst. Sci. Data 2019, 11, 1931–1946. [Google Scholar] [CrossRef]
  46. Cheng, W.; Feng, Q.; Xi, H.; Yin, X.; Cheng, L.; Sindikubwabo, C.; Zhang, B.; Chen, Y.; Zhao, X. Modeling and assessing the impacts of climate change on groundwater recharge in endorheic basins of Northwest China. Sci. Total Environ. 2024, 918, 170829. [Google Scholar] [CrossRef]
  47. Bai, Y.; Zhu, Y.; Liu, Y.; Wang, S. Vegetation Greening and Its Response to a Warmer and Wetter Climate in the Yellow River Basin from 2000 to 2020. Remote Sens. 2024, 16, 790. [Google Scholar] [CrossRef]
  48. Yang, J.; Huang, X. The 30 m annual land cover dataset and its dynamics in China from 1990 to 2019. Earth Syst. Sci. Data 2021, 13, 3907–3925. [Google Scholar] [CrossRef]
  49. Velastegui-Montoya, A.; Montalván-Burbano, N.; Carrión-Mero, P.; Rivera-Torres, H.; Sadeck, L.; Adami, M. Google Earth Engine: A Global Analysis and Future Trends. Remote Sens. 2023, 15, 3675. [Google Scholar] [CrossRef]
  50. Zhang, Y.; Zhang, C.; Wang, Z.; Chen, Y.; Gang, C.; An, R.; Li, J. Vegetation dynamics and its driving forces from climate change and human activities in the Three-River Source Region, China from 1982 to 2012. Sci. Total Environ. 2016, 563–564, 210–220. [Google Scholar] [CrossRef]
  51. Ren, X.; Xu, E.; Smith, C.K.; Vrahnakis, M.; Jing, W.; Zhao, W.; Wang, R.; Jia, X.; Yan, C.; Liu, R. Changes in Surface Runoff and Temporal Dispersion in a Restored Montane Watershed on the Qinghai–Tibetan Plateau. Land 2024, 13, 583. [Google Scholar] [CrossRef]
  52. Cao, S.; Zhang, L.; He, Y.; Zhang, Y.; Chen, Y.; Yao, S.; Yang, W.; Sun, Q. Effects and contributions of meteorological drought on agricultural drought under different climatic zones and vegetation types in Northwest China. Sci. Total Environ. 2022, 821, 153270. [Google Scholar] [CrossRef]
  53. Li, Q.; Bao, Y.; Wang, Z.; Chen, X.; Lin, X. Trade-offs and synergies of ecosystem services in karst multi-mountainous cities. Ecol. Indic. 2024, 159, 111637. [Google Scholar] [CrossRef]
  54. Mehdizadeh, S. Estimation of daily reference evapotranspiration (ETo) using artificial intelligence methods: Offering a new approach for lagged ETo data-based modeling. J. Hydrol. 2018, 559, 794–812. [Google Scholar] [CrossRef]
  55. Zhang, F.; Zeng, B.; Cao, Y.; Li, F.; Tang, Z.; Qi, J. Human activities have markedly altered the pattern and trend of net primary production in the Ili River basin of Northwest China under current climate change. Land Degrad. Dev. 2022, 33, 2585–2595. [Google Scholar] [CrossRef]
  56. Li, H.; He, Y.; Zhang, L.; Cao, S.; Sun, Q. Spatiotemporal changes of Gross Primary Production in the Yellow River Basin of China under the influence of climate-driven and human-activity. Glob. Ecol. Conserv. 2023, 46, e02550. [Google Scholar] [CrossRef]
  57. Dang, X.; Na, Y.; Chi, W.; Zhao, J.; Zhao, Y.; Wang, Y.; Wu, X.; Wang, Y. Spatio-Temporal Evolution of Sandy Land and its Impact on Soil Wind Erosion in the Kubuqi Desert in Recent 30 Years. Front. Environ. Sci. 2022, 10, 950196. [Google Scholar] [CrossRef]
  58. Gunderson, C.A.; Norby, R.J.; Wullschleger, S.D. Acclimation of photosynthesis and respiration to simulated climatic warming in northern and southern populations of Acer saccharum: Laboratory and field evidence. Tree Physiol. 2000, 20, 87–96. [Google Scholar] [CrossRef] [PubMed]
  59. Yamori, W.; Hikosaka, K.; Way, D.A. Temperature response of photosynthesis in C3, C4, and CAM plants: Temperature acclimation and temperature adaptation. Photosynth. Res. 2014, 119, 101–117. [Google Scholar] [CrossRef]
  60. Wen, Y.; Liu, X.; Pei, F.; Li, X.; Du, G. Non-uniform time-lag effects of terrestrial vegetation responses to asymmetric warming. Agric. For. Meteorol. 2018, 252, 130–143. [Google Scholar] [CrossRef]
  61. Ding, Y.; Li, Z.; Peng, S. Global analysis of time-lag and -accumulation effects of climate on vegetation growth. Int. J. Appl. Earth Obs. Geoinf. 2020, 92, 102179. [Google Scholar] [CrossRef]
  62. Kong, D.; Miao, C.; Wu, J.; Zheng, H.; Wu, S. Time lag of vegetation growth on the Loess Plateau in response to climate factors: Estimation, distribution, and influence. Sci. Total Environ. 2020, 744, 140726. [Google Scholar] [CrossRef] [PubMed]
  63. Kuzyakov, Y.; Gavrichkova, O. REVIEW: Time lag between photosynthesis and carbon dioxide efflux from soil: A review of mechanisms and controls. Glob. Chang. Biol. 2010, 16, 3386–3406. [Google Scholar] [CrossRef]
  64. Wan, S.; Xia, J.; Liu, W.; Niu, S. Photosynthetic overcompensation under nocturnal warming enhances grassland carbon sequestration. Ecology 2009, 90, 2700–2710. [Google Scholar] [CrossRef] [PubMed]
  65. Ershadi, A.; McCabe, M.; Evans, J.; Wood, E. Impact of model structure and parameterization on Penman–Monteith type evaporation models. J. Hydrol. 2015, 525, 521–535. [Google Scholar] [CrossRef]
  66. Kong, D.; Miao, C.; Borthwick, A.G.L.; Lei, X.; Li, H. Spatiotemporal variations in vegetation cover on the Loess Plateau, China, between 1982 and 2013: Possible causes and potential impacts. Environ. Sci. Pollut. Res. 2018, 25, 13633–13644. [Google Scholar] [CrossRef]
  67. Anderegg, W.R.L.; Schwalm, C.R.; Biondi, F.; Camarero, J.J.; Koch, G.W.; Litvak, M.; Ogle, K.; Shaw, J.D.; Shevliakova, E.; Williams, A.P.; et al. Pervasive drought legacies in forest ecosystems and their implications for carbon cycle models. Science 2015, 349, 528–532. [Google Scholar] [CrossRef]
  68. Paudel, K.P.; Andersen, P. Response of rangeland vegetation to snow cover dynamics in Nepal Trans Himalaya. Clim. Chang. 2013, 117, 149–162. [Google Scholar] [CrossRef]
  69. Braswell, B.H.; Schimel, D.S.; Linder, E.; Moore, B., III. The response of global terrestrial ecosystems to interannual temperature variability. Science 1997, 278, 870–873. [Google Scholar] [CrossRef]
  70. Zhang, T.; Yang, S.; Guo, R.; Guo, J. Correction: Warming and Nitrogen Addition Alter Photosynthetic Pigments, Sugars and Nutrients in a Temperate Meadow Ecosystem. PLoS ONE 2016, 11, e0158249. [Google Scholar] [CrossRef]
  71. Lehmann, C.E.R.; Anderson, T.M.; Sankaran, M.; Higgins, S.I.; Archibald, S.; Hoffmann, W.A.; Hanan, N.P.; Williams, R.J.; Fensham, R.J.; Felfili, J.; et al. Savanna Vegetation-Fire-Climate Relationships Differ Among Continents. Science 2014, 343, 548–552. [Google Scholar] [CrossRef]
  72. Zhao, X.; Hu, H.; Shen, H.; Zhou, D.; Zhou, L.; Myneni, R.B.; Fang, J. Satellite-indicated long-term vegetation changes and their drivers on the Mongolian Plateau. Landsc. Ecol. 2015, 30, 1599–1611. [Google Scholar] [CrossRef]
  73. Zhou, X.; Yamaguchi, Y.; Arjasakusuma, S. Distinguishing the vegetation dynamics induced by anthropogenic factors using vegetation optical depth and AVHRR NDVI: A cross-border study on the Mongolian Plateau. Sci. Total Environ. 2018, 616–617, 730–743. [Google Scholar] [CrossRef] [PubMed]
  74. Peng, J.; Jiang, H.; Liu, Q.; Green, S.M.; Quine, T.A.; Liu, H.; Qiu, S.; Liu, Y.; Meersmans, J. Human activity vs. climate change: Distinguishing dominant drivers on LAI dynamics in karst region of southwest China. Sci. Total Environ. 2021, 769, 144297. [Google Scholar] [CrossRef] [PubMed]
  75. Liu, J.; Wei, L.; Zheng, Z.; Du, J. Vegetation cover change and its response to climate extremes in the Yellow River Basin. Sci. Total Environ. 2023, 905, 167366. [Google Scholar] [CrossRef]
  76. Wang, Y.; Xiao, C.; Liu, C.; He, N. The response of grassland productivity to atmospheric nitrogen deposition in northern China. Agric. Ecosyst. Environ. 2024, 359, 108764. [Google Scholar] [CrossRef]
  77. Zhang, C.; Long, D.; Liu, T.; Slater, L.J.; Wang, G.; Zuo, D.; Duan, L.; Cui, Y.; Cui, Y. Grassland Greening and Water Resource Availability May Coexist in a Warming Climate in Northern China and the Tibetan Plateau. Earth’s Futur. 2023, 11, e2023EF004037. [Google Scholar] [CrossRef]
  78. Zhang, Y.; Liang, S. Changes in forest biomass and linkage to climate and forest disturbances over Northeastern China. Glob. Chang. Biol. 2014, 20, 2596–2606. [Google Scholar] [CrossRef]
  79. Cai, D.; Ge, Q.; Wang, X.; Liu, B.; Goudie, A.S.; Hu, S. Contributions of ecological programs to vegetation restoration in arid and semiarid China. Environ. Res. Lett. 2020, 15, 114046. [Google Scholar] [CrossRef]
Land 13 01337 g001
Figure 1. Hubao–Egyu Urban Agglomeration ((a), Location; (b), Administration division; (c), Elevation).
Figure 1. Hubao–Egyu Urban Agglomeration ((a), Location; (b), Administration division; (c), Elevation).
Land 13 01337 g001
Land 13 01337 g002
Figure 2. Flow chart for the data analysis.
Figure 2. Flow chart for the data analysis.
Land 13 01337 g002
Land 13 01337 g003
Figure 3. Temporal effects of vegetation response to climate factors. Taking kNDVI in May as an example (pointer), the loop from inside to outside represents (i) no time effect (Lag0 and Acc0), with kNDVI and climate factors in May; (ii) lag effect (Lag1 and Acc0), with kNDVI in May and climate factors with a lag of one month (i.e., April); (iii) accumulation effect (Lag0 and Acc2), with kNDVI in May and climatic factors accumulated for two months (i.e., March + April + May); (iv) bidirectional coupling effect (Lag1 and Acc2), with kNDVI in May is correlated with climatic factors with a lag of one month and an accumulation of two months (i.e., March + April).
Figure 3. Temporal effects of vegetation response to climate factors. Taking kNDVI in May as an example (pointer), the loop from inside to outside represents (i) no time effect (Lag0 and Acc0), with kNDVI and climate factors in May; (ii) lag effect (Lag1 and Acc0), with kNDVI in May and climate factors with a lag of one month (i.e., April); (iii) accumulation effect (Lag0 and Acc2), with kNDVI in May and climatic factors accumulated for two months (i.e., March + April + May); (iv) bidirectional coupling effect (Lag1 and Acc2), with kNDVI in May is correlated with climatic factors with a lag of one month and an accumulation of two months (i.e., March + April).
Land 13 01337 g003
Land 13 01337 g004
Figure 4. Trends in each region of kNDVI from 2000 to 2022. Full meaning of the abbreviation in Legend: ESD: Extremely Significant Degradation; SID: Significant Degradation; STA: Stable; SII: Significant Improvement; ESI: Extremely Significant Improvement.
Figure 4. Trends in each region of kNDVI from 2000 to 2022. Full meaning of the abbreviation in Legend: ESD: Extremely Significant Degradation; SID: Significant Degradation; STA: Stable; SII: Significant Improvement; ESI: Extremely Significant Improvement.
Land 13 01337 g004
Land 13 01337 g005
Figure 5. Partial correlation and the lag and accumulation period of vegetation kNDVI with precipitation (a,d), temperature (b,e), and potential evapotranspiration (c,f) in the Hubao–Egyu Urban Agglomeration. Here, Lag i and Acc j means i months lag and j months accumulation.
Figure 5. Partial correlation and the lag and accumulation period of vegetation kNDVI with precipitation (a,d), temperature (b,e), and potential evapotranspiration (c,f) in the Hubao–Egyu Urban Agglomeration. Here, Lag i and Acc j means i months lag and j months accumulation.
Land 13 01337 g005
Land 13 01337 g006
Figure 6. The driving mechanism of CC and HA on vegetation dynamics. The meaning of abbreviations in the legend is shown in Table S3. Columns 1 to 4 in the submap are the land use maps and Google Earth Maps in 2000 and 2022, respectively, which reflect the land use change and real images. The colors of the land use types are consistent with those in Supplementary Figure S7.
Figure 6. The driving mechanism of CC and HA on vegetation dynamics. The meaning of abbreviations in the legend is shown in Table S3. Columns 1 to 4 in the submap are the land use maps and Google Earth Maps in 2000 and 2022, respectively, which reflect the land use change and real images. The colors of the land use types are consistent with those in Supplementary Figure S7.
Land 13 01337 g006
Land 13 01337 g007
Figure 7. The contribution of CC (a–d) and HA (e–h) to kNDVI in HBEY from 2000 to 2022, based on a modified model. Columns 1 to 4 represent the contribution rate under four scenarios: no temporal effect (No), lag effect (Lag), accumulation effect (Acc) and bidirectional coupling effects (LagAcc).
Figure 7. The contribution of CC (a–d) and HA (e–h) to kNDVI in HBEY from 2000 to 2022, based on a modified model. Columns 1 to 4 represent the contribution rate under four scenarios: no temporal effect (No), lag effect (Lag), accumulation effect (Acc) and bidirectional coupling effects (LagAcc).
Land 13 01337 g007
Land 13 01337 g008
Figure 8. Performance of residual models between kNDVI and climate variables under four scenarios. The diagram depicted the violin diagram of the three evaluation indicators: (a) R2, (b) MAE, (c) RMSE.
Figure 8. Performance of residual models between kNDVI and climate variables under four scenarios. The diagram depicted the violin diagram of the three evaluation indicators: (a) R2, (b) MAE, (c) RMSE.
Land 13 01337 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, X.; Du, G.; Zhang, X.; Li, X.; Lv, S.; He, Y. Vegetation Dynamics and Driving Mechanisms Considering Time-Lag and Accumulation Effects: A Case Study of Hubao–Egyu Urban Agglomeration. Land 2024, 13, 1337. https://doi.org/10.3390/land13091337

AMA Style

Liu X, Du G, Zhang X, Li X, Lv S, He Y. Vegetation Dynamics and Driving Mechanisms Considering Time-Lag and Accumulation Effects: A Case Study of Hubao–Egyu Urban Agglomeration. Land. 2024; 13(9):1337. https://doi.org/10.3390/land13091337

Chicago/Turabian Style

Liu, Xi, Guoming Du, Xiaodie Zhang, Xing Li, Shining Lv, and Yinghao He. 2024. "Vegetation Dynamics and Driving Mechanisms Considering Time-Lag and Accumulation Effects: A Case Study of Hubao–Egyu Urban Agglomeration" Land 13, no. 9: 1337. https://doi.org/10.3390/land13091337

APA Style

Liu, X., Du, G., Zhang, X., Li, X., Lv, S., & He, Y. (2024). Vegetation Dynamics and Driving Mechanisms Considering Time-Lag and Accumulation Effects: A Case Study of Hubao–Egyu Urban Agglomeration. Land, 13(9), 1337. https://doi.org/10.3390/land13091337

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop