Abstract:Leaf phenology, represented at ecosystem scale by the seasonal dynamics of leaf area index (LAI), is a key control on the exchanges of CO2, energy and water between the land and atmosphere. Robust simulation of leaf phenology is thus important for both dynamic global vegetation models (DGVMs) and land-surface representations in climate and Earth System models. There is no general agreement on how leaf phenology should be modelled. However, a recent theoretical advance posits a universal relationship between the time course of 'steady-state' gross primary production (GPP) and LAI - that is, the mutually consistent LAI and GPP that would pertain if weather conditions were held constant. This theory embodies the concept that leaves should be displayed when their presence is most beneficial to plants, combined with the reciprocal relationship of LAI and GPP via (a) the Beer's law dependence of GPP on LAI, and (b) the requirement for GPP to support the allocation of carbon to leaves. Here we develop a global prognostic leaf phenology model, combining this theoretical approach with a parameter-sparse terrestrial GPP model (the P model) that achieves good fit to GPP derived from flux towers in all biomes, and a scheme based on the P model that predicts seasonal maximum LAI as the lesser of an energy-limited rate (maximizing GPP) and a water-limited rate (maximizing the use of available precipitation). The exponential moving average method is used to represent the time lag between leaf allocation and modelled steady-state LAI. The model captures satellite-derived LAI dynamics across biomes at both site and global levels. Since this model outperforms the 15 DGVMs used in the TRENDY project, it could provide a basis for improved representation of leaf-area dynamics in vegetation and climate models.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to establish a global Leaf Area Index (LAI) dynamic model. This model can be independent of satellite data, simulate LAI changes on time scales from seasons to multiple years, and accurately reflect the LAI dynamics both at flux tower sites and globally. Specifically, the paper aims to: 1. **Develop a global leaf area dynamic model**: This model does not rely on satellite data and can simulate LAI changes on time scales from seasons to multiple decades. 2. **Evaluate the performance of the model**: Use in - situ measurements and satellite data to evaluate the model to verify its applicability and accuracy in different biomes around the world. ### Background and Motivation Leaf phenology (i.e., the annual life cycle of leaves, including budding, unfolding, growth, senescence, and abscission) has an important impact on the carbon, energy, and water exchanges between terrestrial ecosystems and the climate. Therefore, accurately simulating leaf phenology is crucial for the terrestrial surface representation in Dynamic Global Vegetation Models (DGVMs) and climate and earth system models. Currently, although there are many models that attempt to simulate leaf phenology, most of these models are based on empirical parameters and lack a clear theoretical basis. In addition, existing models perform poorly in simulating LAI dynamics on a global scale, especially in high - latitude and tropical regions. ### Research Methods The paper adopts a new theoretical approach, combining the inter - relationship between GPP (Gross Primary Productivity) and LAI. The specific steps are as follows: 1. **Derivation of steady - state LAI dynamics**: Through Beer's law (which describes the dependence of GPP on LAI) and the requirement that GPP supports LAI development, the dynamics of steady - state LAI are derived. 2. **Simulation of actual LAI dynamics**: Considering that leaf allocation lags behind steady - state LAI, an exponentially weighted moving - average method is used to simulate the dynamics of actual LAI. 3. **Prediction of the annual maximum LAI**: Based on the energy - limitation and water - limitation equations, predict the annual maximum LAI. ### Key Equations of the Model - **Energy - limitation equation**: \[ f_{\text{APAR}}^{\text{opt}}=\min \left\{1-\frac{z}{k A_0}, \frac{c_a(1 - \chi)}{1.6 D}\left[f_0 P / A_0\right]\right\} \] where: - \( f_{\text{APAR}}^{\text{opt}} \) is the optimal fraction of absorbed photosynthetically active radiation. - \( z \) is a constant, representing the cost of building and maintaining leaves. - \( k \) is the light attenuation coefficient, set to 0.5. - \( A_0 \) is the potential GPP. - \( c_a \) is the partial pressure of atmospheric CO2. - \( \chi \) is the ratio of the partial pressure of CO2 inside the leaf to the partial pressure of atmospheric CO2. - \( D \) is the vapor pressure deficit during the growing season. - \( P \) is the total annual precipitation. - \( f_0 \) is the part of the annual precipitation that can be transpired by plants and is related to the climate aridity index (AI). - **Optimal LAI**: \[ \text{LAI}_{\text{opt}} = -\frac{1}{k} \ln (1 - f_{\text{APAR}}^{\text{opt}}) \] - **Simulation of actual LAI**: \[ \text{LAI}_{\text{sim}}=\alpha \times L_s[t]+(1 - \alpha) \times L_s[t - 1] \] where: - \( \text{LAI}_{\text{sim}} \) is the simulated actual LAI. - \( L_s \) is the steady - state LAI. - \( \alpha \) is the weighting factor, set to 0.067, corresponding to a memory effect of approximately 15 days. ### Conclusion Through the above methods, the paper has developed a model that can be independent of satellite data.

A general model for the seasonal to decadal dynamics of leaf area

A Semi-Prognostic Phenology Model for Simulating Multi-Decadal Dynamics of Global Vegetation Leaf Area Index

A Steady-State Approximation Approach to Simulate Seasonal Leaf Dynamics of Deciduous Broadleaf Forests Via Climate Variables

A Prognostic Vegetation Phenology Model to Predict Seasonal Maximum and Time Series of Global Leaf Area Index Using Climate Variables

Evaluation of Leaf Phenology of Different Vegetation Types From Local to Hemispheric Scale in CLM

Leaf Area Index in Earth System Models: How the Key Variable of Vegetation Seasonality Works in Climate Projections

A leaf age-dependent light use efficiency model for remote sensing the gross primary productivity seasonality over pantropical evergreen broadleaved forests

Seasonality of Leaf Area Index and Photosynthetic Capacity for Better Estimation of Carbon and Water Fluxes in Evergreen Conifer Forests

A New Temperature-Photoperiod Coupled Phenology Module in LPJ-GUESS Model V4.1: Optimizing Estimation of Terrestrial Carbon and Water Processes

Predicting vegetation type through physiological and environmental interactions with leaf traits: evergreen and deciduous forests in an earth system modeling framework

Using Satellite Data to Improve the Leaf Phenology of A Global Terrestrial Biosphere Model

The phenology of leaf quality and its within-canopy variation is essential for accurate modeling of photosynthesis in tropical evergreen forests

Evaluation of key parameters controlling phenology-induced variability of surface fluxes in land surface models

A Seasonally Varying Phenology for High Resolution Simulations with the COSMO-CLM Model

Improved Modeling of Canopy Transpiration for Temperate Forests by Incorporating a LAI-based Dynamic Parametrization Scheme of Stomatal Slope

Developing a parsimonious canopy model (PCM v1.0) to predict forest gross primary productivity and leaf area index of deciduous broad-leaved forest

Influence of vegetation phenology on modelling carbon fluxes in temperate deciduous forest by exclusive use of MODIS time-series data

Phenology-induced energy and carbon fluxes in land surface models

Improved process representation of leaf phenology significantly shifts climate sensitivity of ecosystem carbon balance

Terrestrial biosphere models need better representation of vegetation phenology: results from the North American Carbon Program Site Synthesis

Modeling land surface phenology in a mixed temperate forest using MODIS measurements of leaf area index and land surface temperature