Intra-Annual Cumulative Effects and Mechanisms of Climatic Factors on Global Vegetation Biomes’ Growth

: Previous studies have shown that climate change has signiﬁcant cumulative effects on vegetation growth. However, there remains a gap in understanding the characteristics of cumulative climatic effects on different vegetation types and the underlying driving mechanisms. In this study, using the normalized difference vegetation index data from 1982 to 2015, along with accumulated temperature, precipitation, and solar radiation data, we quantitatively investigated the intra-annual cumulative effects of climatic factors on global vegetation biomes across climatic zones. We also explored the underlying driving mechanisms. The results indicate that precipitation has a longer intra-annual cumulative effect on vegetation, with effects lasting up to 12 months for large percentages of most vegetation biomes. The cumulative effect of solar radiation is mostly concentrated within 0–6 months. Temperature has a shorter cumulative effect, with no signiﬁcant cumulative effect of temperature on large percentages of tree-type vegetation. Compared to other vegetation types, evergreen broadleaf forests, close shrublands, open shrublands, savannas, and woody savannas exhibit more complex cumulative climatic effects. Each vegetation type shows a weak-to-moderate correlation with accumulated precipitation while exhibiting strong-to-extremely-strong positive correlations with accumulated temperature and accumulated solar radiation. The climate-induced regulations of water, heat, and nutrient, as well as the intrinsic mechanisms of vegetation’s tolerance, resistance, and adaptation to climate change, account for the signiﬁcant heterogeneity of cumulative climatic effects across vegetation biomes in different climatic zones. This study contributes to enriching the theoretical understanding of the relationship between vegetation growth and climate change. It also offers crucial theoretical support for developing climate change adaptation strategies and improving future “vegetation-climate” models.


Introduction
Climate change poses one of the greatest challenges to humanity, with intricate mechanisms and far-reaching impacts.As a vital nexus between the soil, atmosphere, and water, green vegetation serves as a sensitive barometer of climate change's effects on terrestrial ecosystems.Through photosynthesis, green vegetation absorbs atmospheric carbon dioxide and sequesters it in biomass and the soil, while releasing oxygen.This process is critical for achieving carbon neutrality in land ecosystems and will be pivotal in humanity's efforts to mitigate global climate change.
Extensive research has examined climate change's impacts on vegetation growth.For example, studies have shown that warming induces earlier onset and extended duration of the growing season, altering vegetation's phenological traits.Additionally, warming enables the expansion of habitats suitable for vegetation growth towards higher latitudes and elevations [1][2][3].Warming can also boost vegetation's carbon sequestration capacity [4,5].However, increases in extreme weather like heatwaves and droughts can indirectly constrain vegetation growth, transforming some carbon sinks into sources [6][7][8].
Climate change is typified by fluctuations in factors like temperature, solar radiation, and precipitation-the three most widely used climatic indices [4,9,10].These factors shape plant growth and distribution by influencing processes including photosynthesis, respiration, and soil carbon decomposition.As key determinants of vegetation growth, temperature enhances productivity in colder regions [4,11], while solar radiation is the prime driver of tropical rainforest productivity [4,8,12].In arid areas, precipitation most strongly governs vegetation productivity changes [13,14].These factors are critical and frequently used to characterize vegetation growth shifts.Research suggests that over half of vegetation change can be attributed to variations in these three factors [4,9].
Climate change does not immediately impact vegetation.Instead, a noticeable time-lag exists.For instance, vegetation growth at middle and low latitudes exhibits a 1-3 month lag in responding to temperature changes, while middle-and high-latitude areas show a 1-3 month lag in responding to solar radiation [15].Precipitation also induces a 1 month lag in arid and semi-arid regions [15].Across China, a 0-1 month lag is typical in vegetation responses to temperature and precipitation changes [16].In the Amazon basin, lags of 0-6 months (temperature), 0-9 months (radiation), and 0-6 months (precipitation) have been observed [17].At the global scale, radiation, soil moisture, and temperature drive 0-3 month lags in vegetation growth [18].These studies indicate that historical climatic conditions often have the strongest impact on current growth.This lag provides vegetation time to respond to changes, preventing drastic effects from abrupt anomalies.
The lag duration is primarily determined by changes in vegetation water content [19].The timing of nutrient and water storage, transport, and excretion within the vegetation also play a role [20].In addition, the long-term balancing and synergistic mechanisms between vegetation growth and environmental factors like soil moisture and carbon/nitrogen [21][22][23], along with intrinsic mechanisms conferring climate resilience like drought tolerance [24] and the thermal adaptation of photosynthesis/respiration [25][26][27], are also important driving factors behind the lagged effects.
Furthermore, recent research shows that climate change has a long-term, cumulative impact on vegetation growth [28].This cumulative effect is widely observed across climates.For instance, temperature has a negligible cumulative effect in most of the Arctic and Northern Temperate Zones but a substantial 2-3 months effect in arid/semi-arid areas.Solar radiation induces a noticeable 3-month cumulative effect in the Arctic and Northern Temperate Zones, while it is shorter in other regions.Precipitation shows complex spatial cumulative effects.This research highlights the prolonged, compounding nature of climate change impacts on vegetation growth and indicates that considering the cumulative effect of climate, the relationship between "climate change-vegetation growth" can be better characterized.On a global scale, the average fitting capability is improved by 12.94%, with a maximum improvement of up to 76.85%.The fitting capability of climatic factors to vegetation growth can be effectively enhanced in approximately 77% of vegetation regions globally, especially in regions with strong interaction between vegetation and climatic factors, such as arid and semi-arid areas.These research findings have been widely confirmed by subsequent studies.
Based on the concept of cumulative effects [28], many scholars have further revealed that the intensity of cumulative climatic effects on vegetation growth varies over different growing seasons [29][30][31][32][33][34].Specifically, a study conducted in Siberia [29] found that climatic factors exhibit an average cumulative effect of 3.6 months on vegetation growth throughout the entire growing season, with more complex ecosystems showing longer cumulative effects.In a study conducted in Inner Mongolia [35], it was discovered that temperature and precipitation have a cumulative effect on vegetation growth for up to two months, but the contribution of climatic factors to the current vegetation growth varies depending on the month, with the highest contribution coming from the current month.It has been found that cumulative drought effects impact the photosynthesis process of 52.11% of global vegetation, with an average cumulative duration of 1-4 months [32].Furthermore, more than three-quarters (78.55%) of grasslands globally are affected by cumulative drought effects, lasting 8-10 months in arid regions of the northern hemisphere (40 • N-55 • N) [31] and 5-10 months in the Qinghai-Tibet Plateau region of China [34].Meanwhile, some researchers have attempted to couple cumulative climatic effects into the "climate changevegetation growth" relationship model in order to better characterize the response of vegetation growth to climate change [30,33,35].
Based on previous research, it can be summarized that climate change has significant cumulative effects on vegetation growth, exhibiting heterogeneity across different climate regions and vegetation types.However, the existing studies mentioned before did not systematically disentangled the effects of cumulative climatic effects on different vegetation types at global a scale, and the underlying driving mechanisms remain poorly understood.Specifically, as climatic conditions transition throughout the year, vegetation growth undergoes a distinct cyclic pattern of changes, these studies have not effectively captured the crucial intra-annual impacts of climatic variability on vegetation growth cycles, which hinders a comprehensive understanding of the relationship between climate change and vegetation growth.Therefore, this study utilized the satellite-derived Normalized Difference Vegetation Index (NDVI) and corresponding climate data (temperature, precipitation, solar radiation) from 1982 to 2015 to investigate the differential responses of major global vegetation types in different climatic zones to the cumulative effects of climatic factors.The present study also attempted to unravel the underlying driving mechanisms.Specifically, this study first extracted the spatial distribution of unchanged vegetation types during the study period.Then, through partial correlation and time-lag analyses, we examined the magnitude and duration of the cumulative effects of temperature, precipitation, and radiation on vegetation growth in these unchanged regions.We analyzed how these effects differ across vegetation types and climatic zones.Finally, we explored the potential mechanisms driving the observed cumulative climatic effects.Overall, this study provides new insights into differential vegetation responses to cumulative climatic effects over time.Our findings bridge knowledge gaps regarding vegetation-climate interactions and support climate-change-mitigation policies.

NDVI Data
This study utilized the Global Inventory Monitoring and Modelling Studies Group (GIMMS)'s third generation NDVI dataset (GIMMS NDVI3g v1.0) to characterize vegetation growth conditions.The NDVI3g dataset has a spatial resolution of 8 km, has a temporal resolution of 15 days, and spans July 1981 to December 2015 for global land areas outside of Antarctica.We processed the data into monthly values using the Maximum Value Composite (MVC) method to minimize contamination from non-vegetation signals.The NDVI3g dataset provides a long time series, extensive spatial coverage, and high accuracy for monitoring and evaluating large-scale vegetation dynamics [15,[36][37][38].Moreover, the statistical analyses in this study were based on the vegetation growing season, which was defined as those months with monthly average air temperature > 0 • C and NDVI > 0.2 [15,39].We used the subset of January 1982 to December 2015 in this study.

Climate Data
The climate data include monthly mean temperature (TEM), total monthly precipitation (PRE), and total monthly shortwave solar radiation (SOLAR) at 0.5 • spatial resolution from January 1982 to December 2015.TEM and PRE are from the Climatic Research Unit (CRU) Time-Series Version 3.24 dataset, which incorporates additional station data and corrects known errors relative to earlier versions [40].SOLAR is from the Princeton Global Meteorological Forcing Dataset version 2, which implements improved downscaling and bias correction methods [41].These datasets have been widely used in global vegetationclimate studies [15,28,36,37,[42][43][44]. The climatic datasets were interpolated into 8 km using kriging interpolation.

Vegetation Type Data
The vegetation classification was based on the 1 km University of Maryland (UMD) land cover product derived from 14 years (1981-1994) of Advanced Very High Resolution Radiometer (AVHRR) data.It provides 14 land use/cover types with a relatively high overall classification accuracy ranging from 81.4% to 90.3% [45].Furthermore, this study obtained the 500 m land use/cover product (MCD12Q1) from The Moderate Resolution Imaging Spectroradiometer (MODIS) satellite for the year of 2015.The UMD scheme with 11 vegetation types (Figure 1) was adopted, with non-vegetated areas masked out.Using the UMD and MCD12Q1 land use/cover data, the unchanged vegetation pixels for each vegetation biome were identified during the study period.Combined with the distribution of four climatic zones (Figure 1), we extracted the average values of NDVI and climatic factors for each vegetation biome in each climatic zone and then analyzed the climatic cumulative effects on different types of vegetation in different climatic zones.In addition, to conduct a more detailed analysis, we utilized the climatic zone data from IPCC (Figure A8) and the global distribution data of the four major crops (maize, rice, soybean, and wheat) provided by the Chinese Academy of Agricultural Sciences.
The climate data include monthly mean temperature (TEM), total monthly precipitation (PRE), and total monthly shortwave solar radiation (SOLAR) at 0.5° spatial resolution from January 1982 to December 2015.TEM and PRE are from the Climatic Research Unit (CRU) Time-Series Version 3.24 dataset, which incorporates additional station data and corrects known errors relative to earlier versions [40].SOLAR is from the Princeton Global Meteorological Forcing Dataset version 2, which implements improved downscaling and bias correction methods [41].These datasets have been widely used in global vegetation-climate studies [15,28,36,37,[42][43][44]. The climatic datasets were interpolated into 8 km using kriging interpolation.

Vegetation Type Data
The vegetation classification was based on the 1 km University of Maryland (UMD) land cover product derived from 14 years (1981-1994) of Advanced Very High Resolution Radiometer (AVHRR) data.It provides 14 land use/cover types with a relatively high overall classification accuracy ranging from 81.4% to 90.3% [45].Furthermore, this study obtained the 500 m land use/cover product (MCD12Q1) from The Moderate Resolution Imaging Spectroradiometer (MODIS) satellite for the year of 2015.The UMD scheme with 11 vegetation types (Figure 1) was adopted, with non-vegetated areas masked out.Using the UMD and MCD12Q1 land use/cover data, the unchanged vegetation pixels for each vegetation biome were identified during the study period.Combined with the distribution of four climatic zones (Figure 1), we extracted the average values of NDVI and climatic factors for each vegetation biome in each climatic zone and then analyzed the climatic cumulative effects on different types of vegetation in different climatic zones.In addition, to conduct a more detailed analysis, we utilized the climatic zone data from IPCC (Figure A8) and the global distribution data of the four major crops (maize, rice, soybean, and wheat) provided by the Chinese Academy of Agricultural Sciences.The data processing procedure and research methods are shown in Figure 2. In the figure, the blue boxes represent data processed using ArcGIS 10.8 software, and red box represents data processed using MATLAB.Please note that, to minimize the loss of data information, the post-processed datasets were resampled to 500 m to match the spatial resolution of MCD12Q1.Our analyses focused on areas where vegetation types remained unchanged during the study period.
The data processing procedure and research methods are shown in Figure 2. In the figure, the blue boxes represent data processed using ArcGIS 10.8 software, and red box represents data processed using MATLAB.Please note that, to minimize the loss of data information, the post-processed datasets were resampled to 500 m to match the spatial resolution of MCD12Q1.Our analyses focused on areas where vegetation types remained unchanged during the study period.

Accumulated Climatic Factors
We defined accumulated climatic factors, including accumulated temperature (ATEM), accumulated solar radiation (ASOLAR), and accumulated precipitation (APRE), to represent the historical magnitude of climate change.These cumulative factors were calculated by summing the monthly climate data exceeding specified thresholds over a given time period, as follows: where ACF is the accumulated climatic factors, namely, ATEM, ASOLAR, and APRE.CF is the time series data of the original climatic factors corresponding to ACF .
base CF is the threshold corresponding to the cumulative climate factors.Based on previous studies, the thresholds for ATEM, ASOLAR, and APRE were set at 0 °C, 0 W/m 2 , and 0 mm, respectively [28].m is the monthly time series, ranging from January 1982 to De- cember 2015.i is the cumulative time period.k is the difference between the current month and the cumulative months.When monthly data were below the thresholds, the cumulative values remained unchanged from the prior month.We set the maximum accumulation period ( n ) to 12 months, consistent with typical annual vegetation cycles [46,47].For example, the cumulative temperature data for December 2000 with an accumulative time of 3 months was obtained by summing the temperature values for October,

Accumulated Climatic Factors
We defined accumulated climatic factors, including accumulated temperature (ATEM), accumulated solar radiation (ASOLAR), and accumulated precipitation (APRE), to represent the historical magnitude of climate change.These cumulative factors were calculated by summing the monthly climate data exceeding specified thresholds over a given time period, as follows: where ACF is the accumulated climatic factors, namely, ATEM, ASOLAR, and APRE.CF is the time series data of the original climatic factors corresponding to ACF.CF base is the threshold corresponding to the cumulative climate factors.Based on previous studies, the thresholds for ATEM, ASOLAR, and APRE were set at 0 • C, 0 W/m 2 , and 0 mm, respectively [28].m is the monthly time series, ranging from January 1982 to December 2015.i is the cumulative time period.k is the difference between the current month and the cumulative months.When monthly data were below the thresholds, the cumulative values remained unchanged from the prior month.We set the maximum accumulation period (n) to 12 months, consistent with typical annual vegetation cycles [46,47].For example, the cumulative temperature data for December 2000 with an accumulative time of 3 months was obtained by summing the temperature values for October, November, and December of 2000.If the temperature for a month was below 0 • C, 0 was added instead.

Partial Correlation Analyses and Cumulative Effect Analyses
Partial correlation analysis was used to examine relationships between NDVI and one typical accumulated climatic factor while controlling others [28].The partial correlation coefficient takes values [−1, 1], representing complete negative correlation to complete positive correlation in statistics.The calculation formula is as follows: where X, Y and Z are the independent variable, dependent variable, and control variable, respectively.R X,Y|Z represents the partial correlation coefficient between X and Y after removing the influence of the Z variable.R XY , R XZ , and R YZ are the Pearson correlation coefficients between X and Y, X and Z, and Y and Z, respectively.When the absolute value of the partial correlation coefficient is between 0.8 and 1.0, it indicates a very strong correlation.A correlation between 0.6 and 0.8 is considered strong, 0.4 to 0.6 indicates a moderate degree of correlation, 0.2 to 0.4 suggests a weak correlation, and 0.0 to 0.2 indicates a very weak correlation or no correlation.
In this study, the accumulated duration ranges from 0 to 12 months.For each assumed accumulated duration, a partial correlation coefficient can be calculated.When a partial correlation coefficient reaches the maximum coefficient of determination, the corresponding accumulated duration is considered the optimal accumulated time (OAT).The partial correlation coefficient at this duration is also considered the optimal partial correlation coefficient (OPCC).The specific formula is as follows: where OPCC is the optimal partial correlation coefficient.i is the assumed accumulated duration, with values ranging from 0 to 12. k is the cumulative duration associated with the maximum deterministic coefficient.PCC k is the partial correlation coefficient at a certain assumed cumulative duration.PCC 2 k is the coefficient of determination for that correlation coefficient.OAT is the optimal accumulated time.
For example, we first calculate the partial correlation coefficients between the NDVI of a certain pixel and 13 cumulative temperature time series data from 0 to 12 months, with the corresponding cumulative precipitation and cumulative solar radiation data as control variables.The cumulative time corresponding to the largest squared value of these 13 partial correlation coefficients is the optimal cumulative time for the temperature of this pixel.

Global Distribution of Unchanged Vegetation Types
According to Figure 1, it can be observed that the unchanged vegetation types were widely distributed in the northern temperate (NTEM) and tropical (TRP) climatic zones during the study period.Most of the ENFs, DNFs, DBFs, MFs, WS, GLs, and CLs were located in NTEM, accounting for 99.20%, 99.13%, 88.65%, 99.10%, 72.40%, 80.32%, and 84.94%, respectively.EBFs were predominantly in the TRP at 95.24%.CS were mainly in the southern temperate (STEM) zone at 73.31%.OSs were more evenly distributed across climatic zones.Sava were mostly found in NTEM and TRP (85.94% combined).
Comparing with the global climate regions defined by the Intergovernmental Panel on Climate Change (IPCC), the following were found (Figure 3

Accumulated Climatic Effect Durations
The cumulative effects of precipitation on vegetation growth showed complex patterns across climatic zones and vegetation types (Figures 4 and 5).Different optimal accumulation durations were observed, varying significantly by region and vegetation type.Globally, ENFs exhibited cumulative precipitation effects mainly at 6-7 months (31.15%) and 12 months (32.79%).EBFs responded primarily at 0-3 months (40.49%).DNFs showed a major response at 2-3 months (58.07%).DBF effects accumulated over 12 months (43.33%) and 0-3 months (31.55%).MFs were mainly distributed in NTEM, and 59.00% of MFs in this climate zone exhibited a cumulative response to precipitation for up to 12 months, resulting in a 12-month cumulative effect on 58.59% of global MF growth.CSs in NTEM accumulated effects over 7-12 months (78.71%) versus 3-10 months in STEM (94.32%), showing a mixed global pattern.OSs, Sava, and WS demonstrated similar heterogeneity as CSs.However, it is worth noting that for Sava and WS in NTEM, 12-month effects were notable at 30.36% and 22.28%.GLs and CLs exhibited primary global responses at 12 months (38.62% and 40.52%) and 0-3 months (41.35% and 32.52%).

Correlations between Accumulated Climatic Factors and Vegetation Growth
Regarding APRE correlation (Figures 8 and A1), vegetation types generally exhibited a weak correlation.In ARC, on average, over 70% of pixels for each vegetation type showed a weak or very weak negative correlation, especially ENFs, DNFs, and GLs (81.41%, 93.02%, 88.52%).For MFs, 52.52% of the pixels showed a weak positive correlation.In NTEM and STEM, over 75% of the pixels for DNFs, DBFs, OSs, Sava, WS, GLs, and CLs were primarily weakly positively correlated.Notably, OSs and GLs had a significant proportion of pixels showing a strong positive correlation in STEM (averaging around 15%).ENFs was mainly weakly correlated in NTEM (43.71% positive, 54.11% negative), but weakly negatively correlated in STEM (76.84%).EBFs exhibited a weak or very weak negative correlation in NTEM (73.05%) but a weak or very weak positive correlation (72.94%) in STEM.The situation was reversed for MFs compared to EBFs.CSs was weakly negatively correlated in NTEM (75.90%) but moderately, strongly positively correlated in STEM (86.21%).In TRP, there was a relatively strong correlation between APRE and vegetation growth.Specifically, on average, around 65% of the pixels for each vegetation type exhibited a weak-to-moderate positive correlation.CSs, OSs, Sava, WS, GLs, and CLs had significant strong positive proportions (67.63%, 23.82%, 21.57%, 32.85%, 23.57%, respectively).

Correlations between Accumulated Climatic Factors and Vegetation Growth
Regarding APRE correlation (Figures A1 and 8), vegetation types generally exhibited a weak correlation.In ARC, on average, over 70% of pixels for each vegetation type showed a weak or very weak negative correlation, especially ENFs, DNFs, and GLs (81.41%, 93.02%, 88.52%).For MFs, 52.52% of the pixels showed a weak positive correlation.In NTEM and STEM, over 75% of the pixels for DNFs, DBFs, OSs, Sava, WS, GLs, and CLs were primarily weakly positively correlated.Notably, OSs and GLs had a significant proportion of pixels showing a strong positive correlation in STEM (averaging around 15%).ENFs was mainly weakly correlated in NTEM (43.71% positive, 54.11% negative), but weakly negatively correlated in STEM (76.84%).EBFs exhibited a weak or very weak negative correlation in NTEM (73.05%) but a weak or very weak positive correlation (72.94%) in STEM.The situation was reversed for MFs compared to EBFs.CSs was weakly negatively correlated in NTEM (75.90%) but moderately, strongly positively correlated in STEM (86.21%).In TRP, there was a relatively strong correlation between APRE and vegetation growth.Specifically, on average, around 65% of the pixels for each vegetation type exhibited a weak-to-moderate positive correlation.CSs, OSs, Sava, WS, GLs, and CLs had significant strong positive proportions (67.63%, 23.82%, 21.57%, 32.85%, 23.57%, respectively).For the ASOLAR correlation (Figures A2 and 8), vegetation types showed significant heterogeneity across climate zones.In ARC, all vegetation biomes exhibited strong to very strong positive correlations (>98% for ENFs, DNFs, MFs, Sava, WS).In NTEM, over 90% of DNFs, DBFs, MFs, WS, and CLs showed strong to very strong positive correlations.ENF, Sava, and GL showed moderate and above-positive correlations (~90%).A total of 86.82% of EBFs exhibited weak-to-strong negative correlations, while CSs correlated strongly positively (53.19%) and negatively (24.08%).In TRP, all vegetation types showed weak or lower correlations with roughly equal positive/negative proportions.In STEM, the correlation patterns were more complex.ENFs, DBFs, and MFs had weakly strongly positive correlations (77.79%, 90.48%, 87.51%), while over half of EBFs and WS had a weakly moderately positive correlation (58.42%, 55.42%).A total of 48.62% of CSs exhibited a weak-to-moderate negative correlation.A total of 81.63% of OSs displayed a weak or lower correlation (with positive correlation accounting for 49.69%).Sava and WS showed predominantly weak-to-moderate correlations (79.60% 78.85%) with positive proportions of 43.12% and 55.42%, respectively.A total of 64.61% of GLs exhibited moderate-to-strong positive correlations, while 51.89% of CLs showed strong and above-negative correlations.Lastly, 38.20% of CLs also exhibited moderate-to-strong positive correlations.
Regarding the correlation with ATEM (Figures A3 and 8), vegetation types differed significantly in strength.In ARC, most ENFs, DNFs, MFs, CSs, OSs, Sava, WS, and GLs exhibited strong-to-very strong positive correlations (>95% pixels).In NTEM, the majority of ENFs, DNFs, MFs, CSs, Sava, and WS showed strong-to-very strong positive correlations (>72%).In TRP, most vegetation biomes displayed weak or lower correlations.In STEM, 80.44% of ENFs showed negative weak-to-strong correlations, and 78.75% of EBF exhibited a weak or lower correlation.A total of 89.78% of DBF showed weak-to-moderate positive correlations.A total of 79.92% of MF showed moderate-to-strong negative correlations.A total of 70.04% of CS exhibited weak-to-moderate negative correlations, a total of 83.03% of OS showed moderate and below-negative correlations.Sava and WS exhibited weak-tomoderate correlations (80.65%, 78.55%) with positive and negative areas being roughly equal.A total of 78.98% of GL exhibited moderate-to-strong correlations (with negative correlations accounting for 47.06%), and 75.43% of CL showed negative moderate-tostrong correlations.

Discussion
This study found that a large proportion of ENFs, MFs, and DBFs had a cumulative effect on precipitation for up to 12 months (weakly positive correlation).These vegetation types are mainly distributed in the Cool Temperate Moist and Boreal Moist climate regions, accounting for over 85% of the total.This may be due to the deep root systems of these vegetation types, which allow them to effectively utilize moisture from the deeper soil layers to sustain growth while reducing dependency on current precipitation.However, the longterm accumulation of precipitation replenishes moisture in deeper soils, promoting the growth of these forests [48].Additionally, the complex ecosystem composition of MFs leads to niche partitioning and growth-promotion mechanisms.MFs exhibit strong adaptation mechanisms, including tolerance, avoidance, and recovery, to cope with drought [49][50][51][52].The higher evapotranspiration, lower ecological connectivity, and increased water demand during the growing season in DBFs likely contributes to its long-term cumulative precipitation effect [48,53,54].DNFs show a shorter cumulative effect of precipitation compared to DBFs, consistent with previous findings [15,48].The study also indicates that a significant proportion of DNF pixels exhibit a short response (2 months) to APRE, mainly in the Boreal Moist climate, suggesting a significant and urgent demand for precipitation by DNFs in this climate zone.EBFs primarily grow in hot and humid tropical and subtropical regions.This study found that EBFs in these two regions mainly exhibited weak or very weak negative correlations with APRE, potentially due to increased evapotranspiration, reduced radiation, and unfavorable temperatures for EBFs' growth [15].GANG C et al. [55] found that the NPP of GLs is positively correlated with annual mean precipitation, while its correlation with annual average temperature varies.This aligns with the current study's findings.The regions with a high GL distribution in NTEM, TRP, and STEM all showed a positive correlation of over 90% with APRE and a varying proportion of positive and negative correlations with ATEM.GLs are mainly distributed in the Cool Temperate Dry climate regions of North America and Eurasia in NTEM, and the Qinghai-Tibet Plateau (Cool Temperate Moist climate region).In the Cool Temperate Dry regions, limited water resources due to scarce precipitation and high evapotranspiration become the most important growth-limiting factors, resulting a moderate positive correlation between grassland and precipitation in NTEM as found by this study.The shallow root systems of GLs enable a quick response to alleviated water limitations from precipitation inputs.Due to the urgent precipitation demand, any rainfall during a growing season significantly promotes growth [56,57].Therefore, previous research [15] discovered a short lag effect (1 month) of precipitation on GLs, while this study found a longer cumulative effect of precipitation (12 months).On the other hand, GLs in the Qinghai-Tibet Plateau (Cool Temperate Moist climate region) can mitigate water limitations through snow-cover melting.As a result, the demand for precipitation by GLs is not strong, leading to a significantly smaller cumulative effect of precipitation (3 months).CLs are primarily distributed in the Cool Temperate Dry, Cool Temperate Moist, and Warm Temperate Dry climate regions.In these areas, approximately half of the CLs exhibited a cumulative effect of precipitation lasting up to 12 months, indicating a long-term influence of precipitation.
Vegetation growing in high-latitude and high-altitude regions is primarily temperaturedependent [15].Temperature influences vegetation growth by increasing photosynthetic enzyme activity [58] and reducing autumn foliage degradation rates [59], promoting vegetation growth.This rapid response [15,48] has led to urgent temperature response mechanisms in these vegetation types.This may be an important reason for the strong positive temperature correlation and short cumulative temperature effects (0-1 month) for Cool Temperate climate vegetation types such as ENFs, DNFs, DBFs, MFs, WS, Sava, and, specifically, GLs on the Qinghai-Tibet Plateau.Furthermore, we found that the GLs growing in the North American Great Plains and the Eurasian continent west of Lake Balkhash to the Caspian Sea were significantly negatively correlated with temperature, and the cumulative response to temperature was relatively long (more than 4 months), which may be because the two regions belong to a temperate humid climate, and the increase in temperature significantly reduces soil moisture, which is unfavorable for local vegetation growth.In contrast, vegetation in arid regions has evolved mechanisms to withstand drought, allowing them to tolerate prolonged environmental anomalies [60].The heterogeneous response durations of EBFs, OSs, CSs, Sava, WS, GLs, and CLs suggest complex temperature response mechanisms.
Increasing growing-season solar radiation enhances photosynthetic enzyme activity [58], providing more photosynthetically active radiation (PAR) to improve plant photosynthetic capacity and carbon sequestration [61,62].Sunlight also delays leaf abscisic acid accumulation, slowing down aging and extending the growing season [63,64].Wen et al. (2019) [28] documented that solar radiation promotes snow melting and indirectly increases surrounding temperatures, thus facilitating vegetation growth.However, this energy conversion process requires a certain accumulation period.It is likely why this study found that ENFs, DNFs, MFs, Sava, and WS, which are widely distributed in ARC and NTEM, had a strong or abovepositive correlation with ATEM and ASOLAR, and they mainly exhibited a 0-1 month cumulative effect of temperature, while the cumulative effect of solar radiation mainly occurred within 2-3 months.This result further validates previous findings [28,48] that vegetation growth in ARC and NTEM can respond rapidly to temperature changes but requires a certain amount of accumulated solar radiation.This study found varying cumulative effects of solar radiation on vegetation growth in low-latitude arid and semi-arid regions, confirming previous research findings [15,48].This relationship is closely related to the complex vegetation composition in these regions.The growth of vegetation in these regions shows a negative correlation with ASOLAR.This could be because an increase in long-term solar radiation promotes evapotranspiration and leads to water loss, which is unfavorable for vegetation growth [24].However, in the EBF of the Amazon rainforest, there is a high positive correlation with solar radiation, consistent with previous study [65], which indicates that solar radiation is an important driving factor for leaf shedding and new leaf growth in the Amazon rainforest.Additionally, this study's results suggest that the cumulative response of grasslands in the North American Great Plains to temperature and solar radiation is longer than that of grasslands in the Eurasian continent west of Lake Balkhash to the Caspian Sea.This indicates that grasslands in the Great Plains have developed stronger mechanisms to tolerate drought and light conditions.We believe this may be related to the widespread presence of C4 grasslands in North America, which have strong characteristics of adaptation to high temperatures, intense light, and arid environments [60].
This study found that, globally, the average OAT for vegetation types is generally ranked as precipitation > solar radiation > temperature, aligning with the research conducted by DING et al. [48].Previous studies have indicated that green vegetation can adapt to seasonal variations in water, heat, and solar radiation resources through leaf shedding [66].For tropical rain forests, the peak of leaf litter occurs during the dry spring or winter seasons and is predominantly regulated by precipitation and radiation, while for boreal forests, the peak of leaf litter occurs in different seasons and is mostly regulated by radiation and temperature [66].Furthermore, studies have shown that temperature drives the distribution of deciduous and evergreen forests through its impact on nitrogen mineralization rates [67].The effect of temperature on aboveground biomass accumulation in DBFs transitions from initial inhibition to subsequent promotion during the maturation stage [67].Additionally, precipitation contributes to biomass accumulation in maturing ENFs but inhibits biomass accumulation in maturing EBFs, DNFs, and DBFs [67].This indicates that different vegetation types exhibit complex and not absolute response mechanisms to climatic factors, which is also a major cause of spatiotemporal heterogeneity in cumulative climatic effects on global vegetation.
This study found stronger cumulative solar radiation effects, followed by temperature, and then precipitation, for widely distributed Arctic (ARC) vegetation types, including OSs, Sava, WS, and GLs, congruent with the results of SHI et al. [29].DBFs, Sava, and WS exhibited prolonged intra-annual response durations, which is consistent with the findings of DING et al. [48].This study indicates that the global distributions of the cumulative effects of climate on GLs and CLs, Sava and WS, and CSs and OSs share a similarity pattern, aligning with previous research [15,48].This may be related to the similarity in community types, intrinsic characteristics, and growth environments shared by these vegetation types [10,68].For Sava and WS, the cumulative effects of each climate factor in NTEM and STEM are more similar to those of OSs, but in TRP, they are more similar to those of GLs.This may be because in TRP, the NDVI information captured by satellites is mainly derived from grasslands, while in NTEM and STEM, it is more derived from trees [69].CSs are predominantly distributed in STEM, while OSs are ubiquitous across climatic zones.This study reveals that there is no clear absolute advantage in the cumulative effects of precipitation on OSs and CSs, which is true also for solar radiation on CSs, indicating complex water and radiation demand for OSs' and CSs' growth [15,48].Additionally, this study reveals that solar radiation mainly exhibits cumulative effects on OSs within 0-3 months, while temperature exhibits cumulative effects of 0 months on OSs and CSs, illustrating the urgent demand of different shrub categories for light and heat energy [15,48].However, the specific reasons behind these effects need further investigation.
Combining the global distribution of four major crops as mapped by the Chinese Academy of Agricultural Sciences, it can be observed that (Figures A5-A7) different crops exhibit varying cumulative effects on temperature, precipitation, and solar radiation.For example, a larger proportion of corn, soybean, and wheat show cumulative effects of precipitation lasting up to 12 months, while rice primarily shows 2-3 months.Rice exhibits cumulative effects of solar radiation mainly within 4-8 months, while corn and soybean show mainly within 4-6 months, and wheat exhibits shorter times (around 3-4 months).Additionally, four major crops exhibit cumulative effects on temperature mainly within 0 and around 6 months.It is worth noting that CLs are a vegetation type that is heavily influenced by human activities [70].People cultivate different crops based on factors such as seasons, climate, soil conditions, and the water requirements of these crops, which can be met through irrigation [71].The growth of CLs is also influenced by local weather conditions, pests and diseases, fertilizer application, and other factors [72][73][74].These factors contribute to the complexity of the duration and relationships of climatic factors' cumulative effects on CLs, which may be more intricate than what is revealed in this study.
In summary, the intrinsic tolerance, resistance, and adaptation mechanisms evolved in diverse vegetation types over time, alongside shifts in water, thermal, and nutrient limitations imposed by changing climatic factors, represent critical determinants of the spatiotemporal heterogeneities in cumulative climatic effects across vegetation biomes.The present findings enrich the theoretical comprehension of vegetation-climate interactions, elucidating the intricate processes governing vegetative responses to climate change.These insights contribute to formulating strategies for addressing and increasing adaptation to climate change.However, it is important to note the complex, multifaceted nature of climate change impacts on vegetation.This study relied solely on linear modeling to characterize vegetation-climate linkages, potentially oversimplifying true interactions.Furthermore, quantitative assessments of anthropogenic effects, extreme events, and other disturbances were not conducted here.Future work could explore nonlinear techniques to further disentangle cumulative climatic effect traits across vegetation types.Examining additional climatic and non-climatic parameters (e.g., human activity, natural disturbances) is also needed to elucidate specific factors eliciting pronounced heterogeneous responses among tropical/Australian vegetation biomes, EBFs, OSs, CSs, Sava, WS, and CLs.This study spans over 30 years of data, and through comparative analyses of land use at the beginning and end of the study period, efforts were made to eliminate regions with significant changes in vegetation types during the research period.However, few pixels might still have exhibited land use changes that were not captured, which could have introduced a certain degree of error into the research results.Future research could attempt to further eliminate these pixels.In addition, global climate change is a complex phenomenon, with cumulative effects influenced by both natural and anthropogenic factors.This study only employed three main climatic factors to investigate the cumulative effects of climate change on vegetation growth, and, thus, our conclusions explain the main cumulative effects reflected by these climatic factors.Future research could incorporate more influencing factors, thereby more comprehensively characterizing this phenomenon.Furthermore, consistent with previous studies [28], heterogeneities exist in the desert of central Australia, as well as the African and the Amazon rainforests.The reason might be that NDVI is sensitive to soil background, which makes it difficult for the vegetation index based on red and near-infrared reflectance to distinguish perennial vegetation in arid and semi-arid ecosystems from the red soil, as these plants often lack contrast between red and near-infrared reflectance [75,76].NDVI tends to saturate over regions like the Amazon rainforest, making it more susceptible to residual cloud or aerosol noise than other vegetation indices [77,78].We argue and strongly recommend that field-based studies be conducted to reduce this confusion.

Conclusions
This study utilized accumulated temperature, precipitation, and solar radiation data, as well as GIMMS NDVI3g data, to examine the intra-annual cumulative climatic effects on global vegetation biomes from 1982 to 2015 quantitatively.The main conclusions are as follows.
Precipitation exerted the longest cumulative effects on vegetation growth worldwide with significant intra-annual heterogeneity.Specifically, precipitation primarily displayed 6-7 and 12-month cumulative effects on ENFs, 0-3-and 12-month effects on GLs and CLs, 2-3-month effects on DNFs, and up to 12-month effects on over 40% of DBFs and MFs.Solar radiation's cumulative effects on vegetation growth were concentrated within 0-6 months, with 3-5-month effects dominating in ENFs and DBFs and 3-month effects in DNFs and MFs.Temperature exhibited the shortest cumulative effects, with no significant impacts observed for many tree types.EBFs, CSs, OSs, Sava, and WS showed the most complex cumulative effects without a clear duration dominance.Vegetation types demonstrated weak-to-moderate correlations with APRE, highlighting the complex role of precipitation.Apart from EBFs, CSs, OSs, and Sava, other vegetation types showed strong-to-extremelystrong positive ATEM and ASOLAR correlations, indicating temperature and radiation dependency during growth cycles.This research indicates that the intrinsic tolerance, resistance, and adaptation mechanisms evolved in vegetation types over time, alongside shifting water, thermal, and nutrient limitations imposed by climate changes, are critical determinants of the spatiotemporal heterogeneities in intra-annual cumulative climatic effects.The results of this study provide key theoretical underpinnings to construct more accurate vegetation-climate relational models moving forward.

8 of 27 28 Figure 4 .
Figure 4.The distribution of accumulation effects durations of climatic factors (a) APRE, (b) ATEM, and (c) ASOLAR on vegetation growth from 1982 to 2015.The white areas represent regions with a p-value greater than 0.05, areas with vegetation type transitions, or non-vegetated areas.

Figure 4 .
Figure 4.The distribution of accumulation effects durations of climatic factors (a) APRE, (b) ATEM, and (c) ASOLAR on vegetation growth from 1982 to 2015.The white areas represent regions with a p-value greater than 0.05, areas with vegetation type transitions, or non-vegetated areas.

Figure 8 .
Figure 8.The distribution of OPCC between global vegetation growth and accumulated climate factors (a) APRE, (b) ATEM, and (c) ASOLAR from 1982 to 2015.The white areas represent regions with a p-value greater than 0.05, areas with vegetation type transitions, or bare areas.

Figure 8 .
Figure 8.The distribution of OPCC between global vegetation growth and accumulated climate factors (a) APRE, (b) ATEM, and (c) ASOLAR from 1982 to 2015.The white areas represent regions with a p-value greater than 0.05, areas with vegetation type transitions, or bare areas.

Figure A2 .
Figure A2.The rectangular tree map illustrates the proportion of the partial correlation strength between the growth of different vegetation types in various climatic zones and ASOLAR.Different color schemes represent different climatic zones, and each vegetation type is labeled in the top left corner of the corresponding box.The map displays the large proportion of partial correlations and merges the smaller proportions into the category "REST".

Figure A2 .
Figure A2.The rectangular tree map illustrates the proportion of the partial correlation strength between the growth of different vegetation types in various climatic zones and ASOLAR.Different color schemes represent different climatic zones, and each vegetation type is labeled in the top left corner of the corresponding box.The map displays the large proportion of partial correlations and merges the smaller proportions into the category "REST".

Figure A3 .
Figure A3.The rectangular tree map illustrates the proportion of the partial correlation strength between the growth of different vegetation types in various climatic zones and ATEM.Different color schemes represent different climatic zones, and each vegetation type is labeled in the top left corner of the corresponding box.The map displays the large proportion of partial correlations and merges the smaller proportions into the category "REST".

Figure A4 .
Figure A4.The climatic conditions of vegetation growth regions.

Figure A3 . 28 Figure A3 .
Figure A3.The rectangular tree map illustrates the proportion of the partial correlation strength between the growth of different vegetation types in various climatic zones and ATEM.Different color schemes represent different climatic zones, and each vegetation type is labeled in the top left corner of the corresponding box.The map displays the large proportion of partial correlations and merges the smaller proportions into the category "REST".

Figure A4 .
Figure A4.The climatic conditions of vegetation growth regions.Figure A4.The climatic conditions of vegetation growth regions.

Figure A4 .
Figure A4.The climatic conditions of vegetation growth regions.Figure A4.The climatic conditions of vegetation growth regions.

Figure A5 .
Figure A5.The percentage of the OAT of accumulation effects of precipitation on four major crops globally (maize, rice, soybean, and wheat).

Figure A6 .
Figure A6.The percentage of the OAT of accumulation effects of solar radiation on four major crops globally (maize, rice, soybean, and wheat).

Figure A7 .
Figure A7.The percentage of the OAT of accumulation effects of temperature on four major crops globally (maize, rice, soybean, and wheat).

Figure A5 . 28 Figure A5 .
Figure A5.The percentage of the OAT of accumulation effects of precipitation on four major crops globally (maize, rice, soybean, and wheat).

Figure A6 .
Figure A6.The percentage of the OAT of accumulation effects of solar radiation on four major crops globally (maize, rice, soybean, and wheat).

Figure A7 .
Figure A7.The percentage of the OAT of accumulation effects of temperature on four major crops globally (maize, rice, soybean, and wheat).

Figure A6 . 28 Figure A5 .
Figure A6.The percentage of the OAT of accumulation effects of solar radiation on four major crops globally (maize, rice, soybean, and wheat).

Figure A6 .
Figure A6.The percentage of the OAT of accumulation effects of solar radiation on four major crops globally (maize, rice, soybean, and wheat).

Figure A7 .
Figure A7.The percentage of the OAT of accumulation effects of temperature on four major crops globally (maize, rice, soybean, and wheat).

Figure A7 .
Figure A7.The percentage of the OAT of accumulation effects of temperature on four major crops globally (maize, rice, soybean, and wheat).

Table A1 .
Percentage of unchanged vegetation in different climate regions.

Table A2 .
The proportion of the top ten optimal accumulation durations for APRE, ATEM, and ASOLAR across all climatic zones for different vegetation types.

Table A1 .
Percentage of unchanged vegetation in different climate regions.

Table A2 .
The proportion of the top ten optimal accumulation durations for APRE, ATEM, and ASOLAR across all climatic zones for different vegetation types.