Next Article in Journal
The Impact of Nanomaterials on Photosynthesis and Antioxidant Mechanisms in Gramineae Plants: Research Progress and Future Prospects
Next Article in Special Issue
Experimental Induction of Extreme Indented Growth Rings (Hazel Wood) in Pinus halepensis Miller by Wide and Long Parallel Bark and Vascular Cambium Woundings
Previous Article in Journal
Transcriptome-Wide N6-Methyladenosine (m6A) Methylation Analyses in a Compatible Wheat–Puccinia striiformis f. sp. tritici Interaction
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 13
Issue 7
10.3390/plants13070980
plants-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

Process-Based Modeling of Phenology and Radial Growth in Pinus tabuliformis in Response to Climate Factors over a Cold and Semi-Arid Region

by
Zihong Man
1,2,
Junzhou Zhang
1,3,*,
Junjun Liu
1,3,
Li Liu
1,3,
Jiqin Yang
1,2,3 and
Zongying Cao
1,3
1
Gansu Liancheng Forest Ecosystem Field Observation and Research Station, Lanzhou University, Lanzhou 730333, China
2
Liancheng National Nature Reserve in Gansu, Lanzhou 730300, China
3
Key Laboratory of Western China’s Environmental Systems (Ministry of Education), College of Earth and Environmental Sciences, Lanzhou University, Lanzhou 730000, China
*
Author to whom correspondence should be addressed.
Plants 2024, 13(7), 980; https://doi.org/10.3390/plants13070980
Submission received: 18 February 2024 / Revised: 18 March 2024 / Accepted: 28 March 2024 / Published: 29 March 2024
(This article belongs to the Special Issue Dendrochronology and Wood Anatomy: Revealing How Plants Respond to Changing Environment)
Browse Figures
Versions Notes

Abstract

:
(1) Background: Climate change significantly impacts the phenology and dynamics of radial tree growth in alpine dryland forests. However, there remains a scarcity of reliable information on the physiological processes of tree growth and cambial phenology in response to long-term climate change in cold and semi-arid regions. (2) Methods: We employed the process-based Vaganov–Shashkin (VS) model to simulate the phenology and growth patterns of Chinese pine (Pinus tabuliformis) in the eastern Qilian Mountains, northeastern Tibetan Plateau. The model was informed by observed temperature and precipitation data to elucidate the relationships between climate factors and tree growth. (3) Results: The simulated tree-ring index closely aligned with the observed tree-ring chronology, validating the VS model’s effectiveness in capturing the climatic influences on radial growth and cambial phenology of P. tabuliformis. The model outputs revealed that the average growing season spanned from mid-April to mid-October and experienced an extension post-1978 due to ongoing warming trends. However, it is important to note that an increase in the duration of the growing season did not necessarily result in a higher level of radial growth. (4) Conclusions: While the duration of the growing season was primarily determined by temperature, the growth rate was predominantly influenced by water conditions during the growing season, making it the most significant factor contributing to ring formation. Our study provides valuable insights into the potential mechanisms underlying tree growth responses to climate change in cold and semi-arid regions.

1. Introduction

Radial tree growth (wood formation), which is the product of cambial cell division and differentiation [1], represents the most important way for long-term carbon sequestration in forests [2,3]. However, the increasing air temperatures and shifting precipitation patterns present a challenge to this carbon sequestration capacity by affecting various physiological processes [4], including decreasing growth rates [5,6] and shifts in the growth phenology of both leaves and cambia [1,7]. Given the projected climate changes expected in the coming decades, it is essential to assess long-term growth phenology and dynamics in response to climate to predict and quantify the future carbon sequestration potential.
To assess growth phenology and seasonal radial growth of trees, researchers employ methods such as monitoring of microcoring and simulation of physiological tree growth models. Microcoring provided a continuous and direct record of phenology and growth dynamics, facilitating the elucidation of uncertainties in radial tree growth at the cellular level, which is crucial for understanding the response mechanism between climate and tree-ring formation [8,9]. However, the labor-intensive nature of direct monitoring constrains its utility for long-term phenological and growth trend analysis [9]. Therefore, there has been a shift towards the development and application of process-based physiological tree growth models to address these challenges [10,11].
The Vaganov–Shashkin (VS) model is one of the most efficient and widely used physiological models for verifying the relationships between tree growth and climate factors [12]. It outputs growth parameters such as cambial phenology and growth rates, which are instrumental in assessing long-term changes over the long-term study period, typically defined by the span of available climate data or tree-ring indices [11,12]. The VS model has proven its precision in simulating tree growth processes across various regions and has provided reliable growth parameters, forming an important foundation for explicating the physiological mechanisms underlying regional tree growth in response to climate changes [13,14,15,16,17,18,19,20]. Its reliability in simulating and evaluating climate–growth relationships across diverse regions and environmental conditions is well documented [17,21,22].
Semi-arid alpine forests, as indicators of global change, are pivotal for soil and water conservation and carbon sequestration in semi-arid regions [23]. Their sensitivity and vulnerability to climate change necessitate a focus on the growth phenology and radial growth of trees in these forests under the background of climate change, as well as an evaluation of their growth responses and adaptations [24]. The Qilian Mountains are a critical area for examining the effects of climate change on semi-arid alpine forests. This region has experienced continuous warming and a notable increase in the frequency of extreme climatic events [25]. Chinese pine (Pinus tabuliformis), the dominant coniferous species of northern China, reaches its westernmost distribution in the eastern Qilian Mountains. Previous studies have shown that the radial growth of P. tabuliformis is highly sensitive to climate change and has been widely used to reconstruct precipitation, temperature, and drought indices over the past several hundred years [26,27]. However, although several studies have focused on the phenology and radial growth dynamics of this species using microcoring monitoring [15,28,29], there remains a lack of reliable information on the long-term phenology and growth dynamics of this species in response to climate change.
In this study, we aim to explore the efficacy of the VS model in replicating physiological parameters such as growth phenology and variability patterns of P. tabuliformis in the eastern Qilian Mountains over recent decades. Based on previous monitoring studies in this region [29,30], we hypothesize that (i) the onset of radial tree growth is primarily determined by the temperature, while water availability is crucial for the end of tree growth, and (ii) water availability has a more significant impact on growth rates rather than temperature in this semi-arid study region.

2. Materials and Methods

2.1. Study Region

Our sampling site was located in the eastern Qilian Mountains, northeastern Tibetan Plateau (36°43′21″ N, 102°37′59″ E, 2070 m asl), marking the transitional zone between the Tibetan Plateau and the Loess Plateau (Figure 1). This area is characterized by pristine mountainous forests with abundant vegetation, displaying a distinct vertical distribution pattern. Within this region, P. tabulaeformis is primarily distributed on dry and infertile south-facing slopes below 2600 m above sea level, representing the westernmost boundary of its distribution in China. The sampling site, centrally located along the forest’s edge on a south-facing slope, lies in a remote area with minimal human impact. The vegetation comprises open stands of mono-specific P. tabulaeformis, with sparse interspecific competition. All sampled trees are healthy and free from pests and diseases. Given the open canopy of the forest (less than 20% coverage), we expect limited tree-to-tree interactions, with abiotic factors predominantly influencing tree growth.

2.2. Meteorological Data

Daily temperature and precipitation data from 1958 to 2011 were sourced from Minhe meteorological station (36°19′48″ N, 102°51′06″ E, 1814 m asl), located approximately 40 km from our sampling site (Figure 1). These data were obtained from the China Meteorological Data Sharing Service System (https://data.cma.cn/, accessed on 18 February 2024) with no missing values. The records indicate a typical semi-arid climate for the study area, with an average annual precipitation of 350.0 mm and an average annual temperature of 8.2 °C for the period of 1958–2011 (Figure 2a). The mean annual temperature showed a significant decreasing trend of 0.30 °C per decade from 1958 to 1978, followed by a significant increasing trend of 0.45 °C per decade from 1978 to 2011 (Figure 2b). Over the study period of 1958–2011, total precipitation exhibited a decrease of 9.35 mm per decade (Figure 2c).

2.3. Tree-Ring Data

Increment core samples were obtained at breast height (1.3 m) from each tree using increment borers, with two or three samples per tree. These samples were processed in the laboratory, mounted, air-dried, and sanded following standard procedures [31]. Tree-rings were cross-dated and measured to a precision of 0.001 mm using a Velmex measuring system and checked by the COFECHA program for quality control of cross-dating [32]. To remove age-related trends and preserve low-frequency signals, conservative detrending curves such as linear or negative exponential curves were applied to the cross-dated raw indices. The indices were standardized by dividing the raw width of the fitted curve value, then averaged to create standard (STD) chronologies using the ARSTAN program [33]. The representativeness of the common signal in the chronology was assessed by calculating the running mean inter-series correlation coefficient (Rbar) and expressed population signal (EPS) over a 50-year window with 25-year overlaps. The subsample signal strength (SSS), with a threshold value of 0.85, was used to determine the reliable period of the chronology (Figure 3a; Table 1). To examine the characteristics of the STD chronology, parameters such as mean sensitivity (MS), standard deviation (SD), mean inter-series correlation (MR), EPS, and signal-to-noise ratio (SNR) were also calculated (Table 1). The STD chronology was used to analyze the climate–growth relationships and validate the simulation results.

2.4. VS Model

The Vaganov–Shashkin (VS) model is a process-based model that operates on several assumptions. Firstly, it posits that cambial activity, and thus the anatomical features of newly formed tracheids, such as radial diameter and cell wall thickness, are primarily influenced by external factors of temperature, photoperiod, and soil moisture. Secondly, the model applies the principle of limiting factors to growth rate calculations, meaning that the growth rate at any given time is constrained by the most limiting factor. Lastly, the VS model simulates variations in growth rate and tree-ring structure in response to current season climatic changes [12].
The model calculates intra-annual radial tree growth as a function of three main factors [12]:
Gr(t) = GrE (t) × min [GrT (t), GrW (t)]
where GrE (t), GrT (t), and GrW (t) represent the partial growth rates on day t of the year, determined independently by solar irradiation (photoperiod), temperature, and soil water content, respectively.
Solar Irradiation: GrE is modeled as a harmonic function that considers latitude, declination angle, and hour angles.
Temperature: Growth is influenced by temperature through a piecewise linear function. Growth initiates when accumulated temperature surpasses the threshold value, Tbeg. Growth is restricted when the temperature falls below a minimum threshold (Tmin). The growth rate increases linearly up to the first optimum temperature (Topt1), remains constant between Topt1 and the second optimum temperature (Topt2), and then decreases linearly until the maximum limiting temperature (Tmax), where growth ceases.
Soil Moisture: Similarly to temperature, the impact of soil moisture on tree growth is quantified using a piecewise function with parameters Wmin, Wopt1, Wopt2, and Wmax. The daily soil water content (dW) is determined by a water balance equation accounting for precipitation, f (P), transpiration (Er), and runoff (Q). The actual precipitation (P) is multiplied by a coefficient (k1) to determine soil infiltration:
dW = f (P) − Er − Q
f (P) = min [k1 × P, Pmax]
The VS model is a multi-parametric model that presents significant challenges in determining the appropriate parameter values to create a simulated chronology that closely matches the observed tree-ring chronology [34]. In this study, certain parameters were established based on observed data, including the coefficient for water infiltration from soil (kr), the depth of the root system (Droot), the maximum daily precipitation entering the soil (Pmax), and the cumulative temperature for growth initiation (Tbeg) [29]. Other critical parameters, including soil moisture (Wmin, Wopt1, Wopt2, Wmax), temperature (Tmin, Topt1, Topt2, Tmax), and coefficients (K1, K2, and K3), were calibrated through an iterative process of refinement and comparison with the observed chronology. Initial values for these parameters were informed by previous studies conducted in similar regions [16]. The chronology was divided into two periods: calibration (1958–1984) and verification (1985–2011). The Pearson correlation coefficient between simulated and observed chronologies served as the metric for evaluating the accuracy of the VS model simulations.

2.5. Data Analysis and Statistics

To examine the relationship between tree growth and climate, correlation analysis was conducted by using STD chronology (observed tree-ring index). The correlation coefficients were calculated between the STD chronology and monthly regional climate data from 1958 to 2011. To consider the impact of climatic conditions on tree growth for both the current and previous years [35], the correlation analysis spanned a 14-month period from the previous August to the current October.
To assess the growth rate during different periods of the growing season, we computed the daily mean growth rates (Gr, GrE, GrT, and GrW) between years. Furthermore, to understand the impact of climatic variability on ring width formation, we compared intra-annual growth variations between years with wide and narrow rings, defined as one standard deviation above or below the mean ring width, respectively.
Additionally, both simulated and observed tree-ring width chronologies were correlated with regional gridded annual mean precipitation data from CRU TS 4.05 (http://climexp.knmi.nl, accessed on 18 February 2024; range of 34°–38° N, 101°–106° E) for the period from the previous August to the current July. This analysis aimed to validate the accuracy of the simulation and to explore the influence of climate on tree growth within the study region.

3. Results

3.1. Tree-Ring and Climate Responses

We found significant positive correlations between the ring-width index and precipitation in the previous August (p < 0.05) and the current July (p < 0.01). In contrast, there were significant negative correlations between the ring-width index and temperature in the previous September and the current May, June, and August (p < 0.05). The strength of these correlations was more significant when considering the annual mean precipitation, with the correlation coefficient between tree growth and precipitation from the previous August to the current July exceeding 0.6 (p < 0.01). Similarly, tree growth exhibited a significant negative correlation with the annual mean temperature for the same period. These results suggest that water availability, as indicated by precipitation, likely plays a crucial role in tree growth of the study region. Spatial correlation analysis further confirmed these findings, showing strong positive correlations between the tree-ring width index and precipitation from the previous August to the current July across extensive areas of the study region (Figure 4a), highlighting the significance of water availability in influencing tree growth across a wide geographic extent.

3.2. Process Model Analysis

Through iterative refinement and correlation comparison between simulated and observed chronologies, we identified the most appropriate physiological parameters for this study (Table 2). Using these optimized parameters, we observed a significant correlation (r = 0.68, p < 0.001, n = 27) between the observed and simulated chronologies during the calibration period 1958–1984. A similar level of agreement (r = 0.64, p < 0.001, n = 27) was found during the verification period 1985–2011 (Figure 5a). When combining these two periods, the simulated chronology explained 45.9% of the variance in the observed chronology (Figure 5b). Additionally, both the observed and simulated chronologies exhibited strong positive correlations with precipitation from the previous August to the current July across extensive regions of the study region (Figure 4). These results underscore the effectiveness of the process-based VS model in simulating tree-ring formation and in capturing large-scale patterns of tree growth responses to climate change in semi-arid conditions, such as those in the eastern Qilian Mountains.

3.3. Model Outputs of Phenology and Growth Rate of Trees

According to the VS model, the overall growth rate (Gr) of P. tabuliformis was primarily impacted by the lower growth rates associated with temperature (GrT) or soil moisture (GrW). The model outputs indicated that both the onset and end of radial growth in P. tabuliformis were determined by temperature-dependent growth rates (Figure 6). The mean onset and end of growth in P. tabuliformis occurred on 13 April (DOY 103) and 11 October (DOY 284), respectively, resulting in an average growing season of 182 days during the study period (1958–2011).
Although not significant (p = 0.21), the onset of the growing season exhibited a delayed trend from 1958 to 1978, with a delay rate of 6.5 days per decade during these two decades (Figure 7a). However, after 1978, the onset of the growing season occurred significantly earlier (p < 0.05), with an advancement rate of 5.1 days per decade (Figure 7a). No distinct trend was observed at the end of the growing season throughout the study period, although a non-significant advancement trend (p = 0.68) from 1958 to 1978 and a delayed trend (p = 0.26) from 1978 to 2011 was found, respectively (Figure 7b). These simulated results suggest an extension of the growing season in P. tabuliformis after 1978.
However, the extension of the growing season did not translate into an increased growth rate. The model outputs revealed that once the growth commenced, the growth rate was predominantly determined by water conditions (GrW) and remained consistent throughout the growing season (Figure 6). When analyzing intra-annual tree growth differences between years with wide and narrow rings using the model, no significant difference was found in the average temperature-dependent growth rate between wide and narrow rings (Figure 8a). This result indicates that temperature was not the primary limiting factor for the formation of wide or narrow rings. However, the mean relative growth rate due to soil moisture was significantly lower in narrow rings than in wide rings during the growing season, highlighting the crucial role of water availability in tree growth during this period (Figure 8b).

4. Discussion

4.1. Impacts of Precipitation on Tree Growth

The tree-ring index exhibited significant positive correlations with monthly precipitation, indicating that water availability is a key determinant of P. tabuliformis growth in the eastern Qilian Mountains. This semi-arid region, characterized by steep slopes, receives low annual precipitation (averaging 350 mm over the past six decades), which is insufficient to meet the trees’ water demands [36]. The higher temperatures in the region likely contributed to increased potential evapotranspiration, further decreasing soil moisture. This is supported by the significant negative correlations between the tree-ring index and summer temperatures during the current growing season, which underscore the importance of water availability for tree growth in this region. Therefore, moisture availability is identified as the primary growth-limiting factor in this area.
The significant positive correlation with precipitation and negative correlation with temperature in the previous August suggests that carbohydrates produced at the end of the previous growing season may contribute to the current year’s growth. This is supported by the monitoring study of photosynthates in the same species within our study region, which reported high photosynthetic rates in August and September at the end of the growing season [37]. Dendroclimatic studies in the eastern Qilian Mountains have also indicated the sensitivity of ring widths to moisture availability in various tree species, including Picea crassifolia [38,39], Picea wilsonii [40], and Juniperus przewalskii [41]. These findings are further supported by intra-annual monitoring studies in the region, which found that the intra-annual radial growth of trees is mainly correlated with the growing season precipitation [29,30,42].

4.2. Effects of Temperature on Phenology of Tree Growth

The VS model outputs provide valuable insights into the growth phenology over the study period. The simulated mean growth onset (13 April, DOY 103) during 1958–2011 aligns with the same period (mid-April) as we monitored in the same region in P. tabuliformis during the 2012–2021 growing seasons [29]. Thus, the VS model provides significant insights into the onset of tree growth during the past decades. However, the simulated mean growth end (11 October, DOY 284) is half a month later than the observed data (late September) [29]. This suggests that the VS model more accurately simulates growth onset than cessation. Previous studies have also highlighted better performance of the VS model in predicting growth onset than growth end [17,18,20], possibly due to a lack of modules accounting for carbon storage processes at the end of the growing season [18]. Additionally, the model may not fully capture the impact of water availability on the end of the growing season in arid or semi-arid regions, resulting in an overestimation of the timing of growth end [29]. Despite these limitations, the VS model demonstrates trends similar to those observed in phenology and is widely used for assessing long-term growth phenologies [14].
The model indicates that temperature is the primary factor determining the growth onset in P. tabuliformis over the eastern Qilian Mountains, with a minimum threshold temperature of 6 °C. This is consistent with our monitoring results that identify temperature as the main factor determining the onset of cell production in coniferous species over the eastern Qilian Mountains, including P. tabuliformis [29,30]. This finding is also in line with studies in temperate and boreal forests, where temperature significantly impacts the onset of tree growth [8,30,43]. For example, the mean threshold temperature for the onset of radial growth in J. przewalskii in the eastern Qilian Mountains is 6.1 °C, similar to findings at treeline in the eastern Italian Alps [44,45].
The model also suggests that temperature controls the end of radial tree growth, with the growth terminating as temperatures decrease to the threshold. This aligns with the monitoring studies in cold regions where temperature plays an important role in the end tree growth [8,46]. However, monitoring studies in the eastern Qilian Mountains indicate that summer drought conditions, driven by high temperature and soil moisture deficits, may determine the end of wood formation before the temperature decreases to the threshold [29,30]. Cambial cell division in a dry environment is closely related to water potential [47,48,49], and water deficit may lead to a decrease in cell turgor, resulting in a termination in cell division and enlargement toward the end of the growing season [50]. Therefore, further research is necessary to better understand the limitations of the VS model and to develop the model to account for the complex interactions between water availability and tree growth in semi-arid regions.
The VS model provides comprehensive growth phenologies for the entire simulated period, surpassing the limitations of radial growth monitoring studies that are typically based on one or a few years [14]. This enables us to analyze the changing trends in phenologies over recent decades and estimate the potential advancement of tree growth under future climate change. Our study indicates a delayed growth onset in the initial decades (1958–1978), followed by advancement in later years (1978–2011), corresponding with observed temperature trends. These findings reinforce the critical role of temperature in influencing the onset of radial growth in the eastern Qilian Mountains.
Predictions based on this study suggest that if temperatures continue to rise at the current rate in the study region (0.45 °C decade−1), the onset of radial tree growth in P. tabuliformis will occur approximately 5.1 days earlier per decade, while the end of tree growth will be delayed by approximately 2.2 days per decade. These projections are similar to those for J. przewalskii in the same region, where the onset advanced by 3.1 days and the end by 1.8 days per decade [30]. However, the uncertainty of the simulation regarding the growth end must be acknowledged, as the model’s consideration of the influence of water on the end of tree growth requires further refinement.

4.3. Influence of Water Availability on Intra-Annual Tree Growth

Our study revealed that an extended growing season, driven by rising temperatures, did not necessarily enhance radial growth. Instead, the growth rate of P. tabuliformis remained dependent on water conditions in the eastern Qilian Mountains. The significant difference in model outputs between years with wide and narrow rings further emphasizes the importance of water availability over temperature in determining tree growth rates in our study region. This is supported by the correlations between the observed STD chronology and climatic responses, indicating that tree growth is influenced by precipitation in July and the previous August. Similar findings have been reported in the region for different species based on microcoring studies, which emphasize the critical role of precipitation during the growing season in influencing growth rate [30,42]. Consistently, both intra-annual [51,52] and inter-annual [41,53,54] tree-ring growth studies have demonstrated a positive correlation between early summer precipitation and drought conditions on the northeastern Tibetan Plateau, suggesting that water availability is the primary factor influencing growth rate and wood production in alpine semi-arid environments.
Spatial correlation analysis revealed strong positive correlations between both the observed and simulated tree-ring width chronologies and precipitation from the previous August to the current July across the study region, highlighting the significance of water availability during this period for regional tree growth. However, when comparing the two correlation fields, the observed tree-ring index exhibited a higher correlation coefficient and a broader correlation field than the simulated index, suggesting that the observed index contains more climate information than the simulated index. This result was not found in previous simulations in J. przewalskii near the study region [16,55]. This discrepancy may arise because the VS model only considers climate information for the current year, omitting previous year influences, which could explain why the simulated index does not fully reflect climate information present in the observed index.

5. Conclusions

The process-based VS model was employed to simulate the phenology and growth dynamics of P. tabuliformis in the eastern Qilian Mountains using temperature and precipitation data from meteorological stations. Our study demonstrates that the simulated tree-ring index closely aligns with the observed tree-ring chronology, indicating the accuracy of the simulations in reflecting the climatic influences on tree growth in P. tabuliformis in this semi-arid region. The outputs from this model agree with the statistical relationships between tree-ring growth and climate factors. The findings suggest that growth phenologies are predominantly determined by threshold temperatures, while the growth rate is primarily influenced by soil moisture availability during the growing season. Notably, temperature was not found to have a significant impact on the rate of radial growth. These insights are in line with previously observed growth patterns for P. tabuliformis and other species within the same region. However, the VS model does not account for the effects of climate conditions from the previous growing season on the subsequent year’s tree growth, which could lead to discrepancies in the simulation results. Future enhancements to the model should incorporate this factor to improve the precision and reliability of growth simulations.

Author Contributions

Conceptualization, J.Z.; methodology, Z.M. and J.Z.; software, Z.M. and J.Z.; validation, Z.M., J.Z., J.L. and J.Y.; formal analysis, Z.M., L.L. and J.Z.; investigation, Z.M., Z.C. and J.Z.; writing—original draft preparation, Z.M. and J.Z.; writing—review and editing, J.Z.; visualization, Z.M. and J.Z.; supervision, J.Z.; project administration, J.Z.; funding acquisition, J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No. U21A2006 and 42001043).

Data Availability Statement

Data are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Deslauriers, A.; Fonti, P.; Rossi, S.; Rathgeber, C.; Gričar, J. Ecophysiology and Plasticity of Wood and Phloem Formation. In Dendroecology, Tree-Ring Analyses Applied to Ecological Studies, AmorosoLori, M.M., Daniels, L.D., Baker, P.J., Camarero, J.J., Eds.; Springer International Publishing: Berlin/Heidelberg, Germany, 2017; pp. 13–33. [Google Scholar]
  2. Pan, Y.; Birdsey, R.A.; Fang, J.; Houghton, R.; Kauppi, P.E.; Kurz, W.A.; Phillips, O.L.; Shvidenko, A.; Lewis, S.L.; Canadell, J.G.; et al. A large and persistent carbon sink in the world’s forests. Science 2011, 333, 988–993. [Google Scholar] [CrossRef] [PubMed]
  3. Pugh, T.A.M.; Rademacher, T.; Shafer, S.L.; Steinkamp, J.; Barichivich, J.; Beckage, B.; Haverd, V.; Harper, A.; Heinke, J.; Nishina, K.; et al. Understanding the uncertainty in global forest carbon turnover. Biogeosciences 2020, 17, 3961–3989. [Google Scholar] [CrossRef]
  4. IPCC. Summary for Policymakers. In Climate Change 2021: The Physical Science Basis; Working Group I Contribution to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change; IPCC: Paris, France, 2021. [Google Scholar]
  5. Adams, H.D.; Macalady, A.K.; Breshears, D.D.; Allen, C.D.; Stephenson, N.L.; Saleska, S.R.; Huxman, T.E.; McDowell, N.G. Climate-induced tree mortality: Earth system consequences. Eos Trans. Am. Geophys. Union 2010, 91, 153–154. [Google Scholar] [CrossRef]
  6. Zhang, J.; Gou, X.; Manzanedo, R.D.; Zhang, F.; Pederson, N. Cambial phenology and xylogenesis of Juniperus przewalskii over a climatic gradient is influenced by both temperature and drought. Agric. For. Meteorol. 2018, 260–261, 165–175. [Google Scholar] [CrossRef]
  7. Fernandez-de-Una, L.; Aranda, I.; Rossi, S.; Fonti, P.; Canellas, I.; Gea-Izquierdo, G. Divergent phenological and leaf gas exchange strategies of two competing tree species drive contrasting responses to drought at their altitudinal boundary. Tree Physiol. 2018, 38, 1152–1165. [Google Scholar] [CrossRef] [PubMed]
  8. Rossi, S.; Anfodillo, T.; Čufar, K.; Cuny, H.E.; Deslauriers, A.; Fonti, P.; Frank, D.; Gricar, J.; Gruber, A.; Huang, J.G.; et al. Pattern of xylem phenology in conifers of cold ecosystems at the Northern Hemisphere. Glob. Chang. Biol. 2016, 22, 3804–3813. [Google Scholar] [CrossRef] [PubMed]
  9. Zhang, J.; Gou, X.; Pederson, N.; Zhang, F.; Niu, H.; Zhao, S.; Wang, F. Cambial phenology in Juniperus przewalskii along different altitudinal gradients in a cold and arid region. Tree Physiol. 2018, 38, 840–852. [Google Scholar] [CrossRef] [PubMed]
  10. Valentine, H.T.; Makela, A. Bridging process-based and empirical approaches to modeling tree growth. Tree Physiol. 2005, 25, 769–779. [Google Scholar] [CrossRef]
  11. Shishov, V.V.; Tychkov, I.I.; Popkova, M.I.; Ilyin, V.A.; Bryukhanova, M.V.; Kirdyanov, A.V. VS-oscilloscope: A new tool to parameterize tree radial growth based on climate conditions. Dendrochronologia 2016, 39, 42–50. [Google Scholar] [CrossRef]
  12. Evans, M.N.; Reichert, B.K.; Kaplan, A.; Anchukaitis, K.J.; Vaganov, E.A.; Hughes, M.K.; Cane, M.A. A forward modeling approach to paleoclimatic interpretation of tree-ring data. J. Geophys. Res.: Biogeosciences 2006, 111, G03008. [Google Scholar] [CrossRef]
  13. He, M.; Shishov, V.; Kaparova, N.; Yang, B.; Bräuning, A.; Grießinger, J. Process-based modeling of tree-ring formation and its relationships with climate on the Tibetan Plateau. Dendrochronologia 2017, 42, 31–41. [Google Scholar] [CrossRef]
  14. Yang, B.; He, M.; Shishov, V.; Tychkov, I.; Vaganov, E.; Rossi, S.; Ljungqvist, F.C.; Brauning, A.; Griessinger, J. New perspective on spring vegetation phenology and global climate change based on Tibetan Plateau tree-ring data. Proc. Natl. Acad. Sci. USA 2017, 114, 6966–6971. [Google Scholar] [CrossRef] [PubMed]
  15. Gao, J.; Yang, B.; He, M.; Shishov, V. Intra-annual stem radial increment patterns of Chinese pine, Helan Mountains, Northern Central China. Trees 2019, 33, 751–763. [Google Scholar] [CrossRef]
  16. Zhang, J.; Gou, X.; Zhang, Y.; Lu, M.; Xu, X.; Zhang, F.; Liu, W.; Gao, L. Forward modeling analyses of Qilian Juniper (Sabina przewalskii) growth in response to climate factors in different regions of the Qilian Mountains, northwestern China. Trees-Struct Funct 2016, 30, 175–188. [Google Scholar] [CrossRef]
  17. Butto, V.; Shishov, V.; Tychkov, I.; Popkova, M.; He, M.; Rossi, S.; Deslauriers, A.; Morin, H. Comparing the Cell Dynamics of Tree-Ring Formation Observed in Microcores and as Predicted by the Vaganov-Shashkin Model. Front. Plant Sci. 2020, 11, 555055. [Google Scholar] [CrossRef] [PubMed]
  18. Jevšenak, J.; Tychkov, I.; Gričar, J.; Levanič, T.; Tumajer, J.; Prislan, P.; Arnič, D.; Popkova, M.; Shishov, V.V. Growth-limiting factors and climate response variability in Norway spruce (Picea abies L.) along an elevation and precipitation gradients in Slovenia. Int. J. Biometeorol. 2020, 65, 311–324. [Google Scholar] [CrossRef] [PubMed]
  19. Tumajer, J.; Begović, K.; Čada, V.; Jenicek, M.; Lange, J.; Mašek, J.; Kaczka, R.J.; Rydval, M.; Svoboda, M.; Vlček, L.; et al. Ecological and methodological drivers of non-stationarity in tree growth response to climate. Glob. Chang. Biol. 2022, 29, 462–476. [Google Scholar] [CrossRef] [PubMed]
  20. Tumajer, J.; Kašpar, J.; Kuželová, H.; Shishov, V.V.; Tychkov, I.I.; Popkova, M.I.; Vaganov, E.A.; Treml, V. Forward Modeling Reveals Multidecadal Trends in Cambial Kinetics and Phenology at Treeline. Front. Plant Sci. 2021, 12, 613643. [Google Scholar] [CrossRef] [PubMed]
  21. Breitenmoser, P.; Bronnimann, S.; Frank, D. Forward modelling of tree-ring width and comparison with a global network of tree-ring chronologies. Clim. Past 2014, 10, 437–449. [Google Scholar] [CrossRef]
  22. Touchan, R.; Shishov, V.V.; Meko, D.M.; Nouiri, I.; Grachev, A. Process based model sheds light on climate sensitivity of Mediterranean tree-ring width. Biogeosciences 2012, 9, 965–972. [Google Scholar] [CrossRef]
  23. Cao, Z.; Zhang, J.; Gou, X.; Wang, Y.; Sun, Q.; Yang, J.; Manzanedo, R.D.; Pederson, N. Increasing forest carbon sinks in cold and arid northeastern Tibetan Plateau. Sci. Total Environ. 2023, 905, 167168. [Google Scholar] [CrossRef] [PubMed]
  24. Gao, L.; Gou, X.; Deng, Y.; Wang, Z.; Gu, F.; Wang, F. Increased growth of Qinghai spruce in northwestern China during the recent warming hiatus. Agric. For. Meteorol. 2018, 260-261, 9–16. [Google Scholar] [CrossRef]
  25. Liu, L.-Y.; Wang, X.-J.; Gou, X.-H.; Yang, M.-X.; Zhang, Z.-H. Projections of surface air temperature and precipitation in the 21st century in the Qilian Mountains, Northwest China, using REMO in the CORDEX. Adv. Clim. Chang. Res. 2022, 13, 344–358. [Google Scholar] [CrossRef]
  26. Liu, Y.; Wang, Y.; Li, Q.; Song, H.; Linderholm, H.W.; Leavitt, S.W.; Wang, R.; An, Z. Tree-ring stable carbon isotope-based May–July temperature reconstruction over Nanwutai, China, for the past century and its record of 20th century warming. Quat. Sci. Rev. 2014, 93, 67–76. [Google Scholar] [CrossRef]
  27. Chen, F.; Yuan, Y.; Zhang, R.; Qin, L. A tree-ring based drought reconstruction (AD 1760–2010) for the Loess Plateau and its possible driving mechanisms. Glob. Planet. Chang. 2014, 122, 82–88. [Google Scholar] [CrossRef]
  28. Zeng, Q.; Rossi, S.; Yang, B. Effects of Age and Size on Xylem Phenology in Two Conifers of Northwestern China. Front. Plant Sci. 2017, 8, 2264. [Google Scholar] [CrossRef] [PubMed]
  29. Wang, L. Intra-Annual Radial Growth of Three Coniferous Species and Their Responses to Climate Factors in the Eastern Qilian Mountains; Lanzhou University: Lanzhou, Chian, 2023. [Google Scholar]
  30. Zhang, J.; Gou, X.; Alexander, M.R.; Xia, J.; Wang, F.; Zhang, F.; Man, Z.; Pederson, N. Drought limits wood production of Juniperus przewalskii even as growing seasons lengthens in a cold and arid environment. Catena 2021, 196, 104936. [Google Scholar] [CrossRef]
  31. Stokes, M.; Smiley, T. An Introduction to Tree-Ring Dating; UA Press: Tucson, AZ, USA, 1968. [Google Scholar]
  32. Holmes, R.L. Computer-assisted quality control in tree-ring dating and measurement. Tree-Ring Bull. 1983, 43, 69–78. [Google Scholar]
  33. Cook, E.R. A Time-Series Analysis Approach to Tree-Ring Standardization; The University of Arizona: Tucson, AZ, USA, 1985. [Google Scholar]
  34. Vaganov, E.A.; Hughes, M.K.; Shashkin, A.V. Introduction and Factors Influencing the Seasonal Growth of Trees; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  35. Fritts, H. Tree Rings and Climate; Academic Press: Cambridge, MA, USA, 1976. [Google Scholar]
  36. Liu, J.; Gou, X.; Gunina, A.; Long, X.-E.; Zhang, F.; Zhang, J. Soil nitrogen pool drives plant tissue traits in alpine treeline ecotones. For. Ecol. Manag. 2020, 477, 118490. [Google Scholar] [CrossRef]
  37. Wang, F.; Zhang, F.; Gou, X.; Fonti, P.; Xia, J.; Cao, Z.; Liu, J.; Wang, Y.; Zhang, J. Seasonal variations in leaf-level photosynthesis and water use efficiency of three isohydric to anisohydric conifers on the Tibetan Plateau. Agric. For. Meteorol. 2021, 308-309, 108581. [Google Scholar] [CrossRef]
  38. Gao, L.L.; Gou, X.H.; Deng, Y.; Yang, M.X.; Zhao, Z.Q.; Cao, Z.Y. Dendroclimatic Response of Picea crassifolia along an Altitudinal Gradient in the Eastern Qilian Mountains, Northwest China. Arct. Antarct. Alp. Res. 2013, 45, 491–499. [Google Scholar] [CrossRef]
  39. Zhang, W.; Gou, X.; Liu, W.; Li, J.; Su, J.; Dilawar, N.; Zhu, F.; Xia, J.; Du, M.; Wang, L.; et al. Divergent tree radial growth patterns of Qinghai spruce (Picea crassifolia) at the alpine timberline along a moisture gradient in the Qilian mountains, Northwest China. Agric. For. Meteorol. 2023, 328, 109240. [Google Scholar] [CrossRef]
  40. Fang, K.Y.; Gou, X.H.; Chen, F.H.; Yang, M.X.; Li, J.B.; He, M.S.; Zhang, Y.; Tian, Q.H.; Peng, J.F. Drought variations in the eastern part of northwest China over the past two centuries: Evidence from tree rings. Clim. Res. 2009, 38, 129–135. [Google Scholar] [CrossRef]
  41. Gou, X.; Deng, Y.; Gao, L.; Chen, F.; Cook, E.; Yang, M.; Zhang, F. Millennium tree-ring reconstruction of drought variability in the eastern Qilian Mountains, northwest China. Clim. Dyn. 2015, 45, 1761–1770. [Google Scholar] [CrossRef]
  42. Zhang, J.; Alexander, M.R.; Gou, X.; Deslauriers, A.; Fonti, P.; Zhang, F.; Pederson, N. Extended xylogenesis and stem biomass production in Juniperus przewalskii Kom. during extreme late-season climatic events. Ann. For. Sci. 2020, 77, 99. [Google Scholar] [CrossRef]
  43. Huang, J.G.; Ma, Q.; Rossi, S.; Biondi, F.; Deslauriers, A.; Fonti, P.; Liang, E.; Makinen, H.; Oberhuber, W.; Rathgeber, C.B.K.; et al. Photoperiod and temperature as dominant environmental drivers triggering secondary growth resumption in Northern Hemisphere conifers. Proc. Natl. Acad. Sci. USA 2020, 117, 20645–20652. [Google Scholar] [CrossRef] [PubMed]
  44. Deslauriers, A.; Rossi, S.; Anfodillo, T.; Saracino, A. Cambial phenology, wood formation and temperature thresholds in two contrasting years at high altitude in southern Italy. Tree Physiol. 2008, 28, 863–871. [Google Scholar] [CrossRef] [PubMed]
  45. Rossi, S.; Deslauriers, A.; Anfodillo, T.; Carraro, V. Evidence of threshold temperatures for xylogenesis in conifers at high altitudes. Oecologia 2007, 152, 1–12. [Google Scholar] [CrossRef]
  46. Rossi, S.; Girard, M.J.; Morin, H. Lengthening of the duration of xylogenesis engenders disproportionate increases in xylem production. Glob. Chang. Biol. 2014, 20, 2261–2271. [Google Scholar] [CrossRef] [PubMed]
  47. Ziaco, E.; Truettner, C.; Biondi, F.; Bullock, S. Moisture-driven xylogenesis in Pinus ponderosa from a Mojave Desert mountain reveals high phenological plasticity. Plant Cell Environ. 2018, 41, 823–836. [Google Scholar] [CrossRef] [PubMed]
  48. Balducci, L.; Cuny, H.E.; Rathgeber, C.B.; Deslauriers, A.; Giovannelli, A.; Rossi, S. Compensatory mechanisms mitigate the effect of warming and drought on wood formation. Plant Cell Environ. 2016, 39, 1338–1352. [Google Scholar] [CrossRef] [PubMed]
  49. Vieira, J.; Rossi, S.; Campelo, F.; Freitas, H.; Nabais, C. Xylogenesis of Pinus pinaster under a Mediterranean climate. Ann. For. Sci. 2014, 71, 71–80. [Google Scholar] [CrossRef]
  50. Steppe, K.; Sterck, F.; Deslauriers, A. Diel growth dynamics in tree stems: Linking anatomy and ecophysiology. Trends Plant Sci. 2015, 20, 335–343. [Google Scholar] [CrossRef] [PubMed]
  51. He, M.; Yang, B.; Wang, Z.; Bräuning, A.; Pourtahmasi, K.; Oladi, R. Climatic forcing of xylem formation in Qilian juniper on the northeastern Tibetan Plateau. Trees 2016, 30, 923–933. [Google Scholar] [CrossRef]
  52. Ren, P.; Rossi, S.; Gricar, J.; Liang, E.; Cufar, K. Is precipitation a trigger for the onset of xylogenesis in Juniperus przewalskii on the north-eastern Tibetan Plateau? Ann. Bot. 2015, 115, 629–639. [Google Scholar] [CrossRef] [PubMed]
  53. Liang, E.Y.; Shao, X.M.; Liu, X.H. Annual Precipitation Variation Inferred from Tree Rings since Ad 1770 for the Western Qilian Mts., Northern Tibetan Plateau. Tree-Ring Res. 2009, 65, 95–103. [Google Scholar] [CrossRef]
  54. Shao, X.; Xu, Y.; Yin, Z.Y.; Liang, E.; Zhu, H.; Wang, S. Climatic implications of a 3585-year tree-ring width chronology from the northeastern Qinghai-Tibetan Plateau. Quat. Sci. Rev. 2010, 29, 2111–2122. [Google Scholar] [CrossRef]
  55. Gou, X.; Zhou, F.; Zhang, Y.; Chen, Q.; Zhang, J. Forward modeling analysis of regional scale tree-ring patterns around the northeastern Tibetan Plateau, Northwest China. Biogeosciences Discuss. 2013, 10, 9969–9988. [Google Scholar] [CrossRef]
Plants 13 00980 g001
Figure 1. Location of sampling site (red triangle) and nearest national meteorological station (black circle).
Figure 1. Location of sampling site (red triangle) and nearest national meteorological station (black circle).
Plants 13 00980 g001
Plants 13 00980 g002
Figure 2. (a) The monthly mean temperature (symbol line) and precipitation (bars) during the period 1958–2011 recorded by the Minhe meteorological station. Annual mean temperature (b) and total precipitation (c) (grey solid line) and their trends (dash line) during the period 1958–2011.
Figure 2. (a) The monthly mean temperature (symbol line) and precipitation (bars) during the period 1958–2011 recorded by the Minhe meteorological station. Annual mean temperature (b) and total precipitation (c) (grey solid line) and their trends (dash line) during the period 1958–2011.
Plants 13 00980 g002
Plants 13 00980 g003
Figure 3. (a) Tree-ring width chronology (grey line) with 21-year running averages (red line) and sample depth (shaded area) of the chronology. The dashed line shows the year when SSS exceeds 0.85. (b) Correlations between the observed STD chronology and monthly temperature and precipitation; p denotes the previous year (e.g., p8 is the previous August), and C denotes the current year (e.g., C1 is the current January), p8C7 represents the mean value from previous August to current July. Bars surpassing the grey area signify correlations above the 95 % confidence level.
Figure 3. (a) Tree-ring width chronology (grey line) with 21-year running averages (red line) and sample depth (shaded area) of the chronology. The dashed line shows the year when SSS exceeds 0.85. (b) Correlations between the observed STD chronology and monthly temperature and precipitation; p denotes the previous year (e.g., p8 is the previous August), and C denotes the current year (e.g., C1 is the current January), p8C7 represents the mean value from previous August to current July. Bars surpassing the grey area signify correlations above the 95 % confidence level.
Plants 13 00980 g003
Plants 13 00980 g004
Figure 4. Spatial correlations between gridded precipitation from the previous August to the current July of CRU TS 4.05 and the (a) observed and (b) simulated tree-ring width chronology. Insignificant correlations (p > 0.05) are masked. The triangle and circle indicate the sampling site and Minhe meteorological station, respectively.
Figure 4. Spatial correlations between gridded precipitation from the previous August to the current July of CRU TS 4.05 and the (a) observed and (b) simulated tree-ring width chronology. Insignificant correlations (p > 0.05) are masked. The triangle and circle indicate the sampling site and Minhe meteorological station, respectively.
Plants 13 00980 g004
Plants 13 00980 g005
Figure 5. (a) Simulated (blue lines) and observed (red lines) chronologies for the calibration (1958–1984) and verification (1985–2011) periods. (b) The relationship between the observed and simulated chronologies during 1958–2011.
Figure 5. (a) Simulated (blue lines) and observed (red lines) chronologies for the calibration (1958–1984) and verification (1985–2011) periods. (b) The relationship between the observed and simulated chronologies during 1958–2011.
Plants 13 00980 g005
Plants 13 00980 g006
Figure 6. Simulated radial growth rates dependent on solar radiation (GrE), temperature (GrT), soil moisture (GrW), and the combined effect of all factors (Gr) during 1958–2011.
Figure 6. Simulated radial growth rates dependent on solar radiation (GrE), temperature (GrT), soil moisture (GrW), and the combined effect of all factors (Gr) during 1958–2011.
Plants 13 00980 g006
Plants 13 00980 g007
Figure 7. Simulated (a) onset and (b) end of growth (solid line) with their respective trends (dashed line) during 1958–2011.
Figure 7. Simulated (a) onset and (b) end of growth (solid line) with their respective trends (dashed line) during 1958–2011.
Plants 13 00980 g007
Plants 13 00980 g008
Figure 8. Simulated mean growth rate dependent on (a) temperature and (b) soil moisture for years with wide and narrow rings, respectively. The shaded areas represent plus/minus one standard deviation.
Figure 8. Simulated mean growth rate dependent on (a) temperature and (b) soil moisture for years with wide and narrow rings, respectively. The shaded areas represent plus/minus one standard deviation.
Plants 13 00980 g008
Table 1. Statistical characteristics and common interval analyses results of the STD chronology.
Table 1. Statistical characteristics and common interval analyses results of the STD chronology.
Statistical ParametersValue
Mean sensitivity (MS)0.507
Inter-series correlation0.827
Mean segment length (yr)84.3
Absent rings (%)0.742
All series Rbar0.613
Standard deviation0.160
Expressed population signal (EPS)0.978
Signal-to-noise ratio (SNR)44.441
Year of SSS > 0.85/yr (cores)1879 (4)
Table 2. Model parameters used in this study.
Table 2. Model parameters used in this study.
ParameterDescription (Units)Value
TminMinimum temperature for tree growth (°C)6
Topt1Lower limit of optimal temperatures (°C) 14
Topt2Upper limit of optimal temperatures (°C)18
TmaxMaximum temperature for tree growth (°C) 30
WminMinimum soil moisture for tree growth (v/v)0.04
Wopt1Lower limit of optimal soil moisture (v/v)0.29
Wopt2Upper limit of optimal soil moisture (v/v)0.8
WmaxMaximum soil moisture for tree growth (v/v)0.9
TbegTemperature sum for initiation of growth (°C), period = 10 days55
DrootDepth of root system (mm)1000
PmaxMaximum daily precipitation for saturated soil (mm)20
k1Fraction of precipitation penetrating soil 0.86
k2First coefficient for calculation of transpiration (mm/day)0.12
k3Second coefficient for calculation of transpiration (1/degree)0.176
krCoefficient for water infiltration from soil 0.002
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

Man, Z.; Zhang, J.; Liu, J.; Liu, L.; Yang, J.; Cao, Z. Process-Based Modeling of Phenology and Radial Growth in Pinus tabuliformis in Response to Climate Factors over a Cold and Semi-Arid Region. Plants 2024, 13, 980. https://doi.org/10.3390/plants13070980

AMA Style

Man Z, Zhang J, Liu J, Liu L, Yang J, Cao Z. Process-Based Modeling of Phenology and Radial Growth in Pinus tabuliformis in Response to Climate Factors over a Cold and Semi-Arid Region. Plants. 2024; 13(7):980. https://doi.org/10.3390/plants13070980

Chicago/Turabian Style

Man, Zihong, Junzhou Zhang, Junjun Liu, Li Liu, Jiqin Yang, and Zongying Cao. 2024. "Process-Based Modeling of Phenology and Radial Growth in Pinus tabuliformis in Response to Climate Factors over a Cold and Semi-Arid Region" Plants 13, no. 7: 980. https://doi.org/10.3390/plants13070980

APA Style

Man, Z., Zhang, J., Liu, J., Liu, L., Yang, J., & Cao, Z. (2024). Process-Based Modeling of Phenology and Radial Growth in Pinus tabuliformis in Response to Climate Factors over a Cold and Semi-Arid Region. Plants, 13(7), 980. https://doi.org/10.3390/plants13070980

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