Articles | Volume 13, issue 4
Data description paper
06 Apr 2021
Data description paper |  | 06 Apr 2021

Standardized flux seasonality metrics: a companion dataset for FLUXNET annual product

Linqing Yang and Asko Noormets

Phenological events are integrative and sensitive indicators of ecosystem processes that respond to climate, water and nutrient availability, disturbance, and environmental change. The seasonality of ecosystem processes, including biogeochemical fluxes, can similarly be decomposed to identify key transition points and phase durations, which can be determined with high accuracy, and are specific to the processes of interest. As the seasonality of different processes differ, it can be argued that the interannual trends and responses to environmental forcings can be better described through the fluxes' own temporal characteristics than through correlation to traditional phenological events like bud break or leaf coloration. Here we present a global dataset of seasonality or phenological metrics calculated for gross primary productivity (GPP), ecosystem respiration (RE), latent heat (LE), and sensible heat (H) calculated for the FLUXNET2015 Dataset of about 200 sites and 1500 site years of data. The database includes metrics (i) on an absolute flux scale for comparisons with flux magnitudes and (ii) on a normalized scale for comparisons of change rates across different fluxes. Flux seasonality was characterized by fitting a single-pass double-logistic model to daily flux integrals, and the derivatives of the fitted time series were used to extract the phenological metrics marking key turning points, season lengths, and rates of change. Seasonal transition points could be determined with a 90 % confidence interval of 6–11 d for GPP, 8–14 d for RE, 10–15 d for LE, and 15–23 d for H. The phenology metrics derived from different partitioning methods diverged, at times significantly.

This Flux Seasonality Metrics Database (FSMD) can be accessed at the US Department of Energy's (DOE) Environmental Systems Science Data Infrastructure for a Virtual Ecosystem (ESS-DIVE,; Yang and Noormets, 2020). We hope that it will facilitate new lines of research, including (1) validating and benchmarking ecosystem process models, (2) parameterizing satellite remote sensing phenology and PhenoCam products, (3) optimizing phenological models, and (4) generally expanding the toolset for interpreting ecosystems responses to changing climate.

Please read the corrigendum first before continuing.

1 Introduction

Phenology, the timing of life-cycle events and phases of plants and animals, and their relationship with the environment, especially climate (Lieth, 1974; Piao et al., 2019), is an important indicator of ecosystem dynamics. It is an integrating record of the effects of global warming and other environmental changes on biological processes (Noormets et al., 2009; Post and Stenseth, 1999; Weltzin et al., 2020). Current phenology studies focus primarily on structural changes such as bud break, flowering, leaf coloring, and leaf fall. However, the functional aspects of plant activities, although invisible, also provide quantitative measures of plant responses to changes in environmental conditions and underlie the structural changes (Fitzjarrald et al., 2001; Schwartz, 2003; Schwartz and Crawford, 2013). Ecosystem processes such as biogeochemical fluxes also show seasonal changes and can be decomposed to key transition dates and phase durations, which characterize the exchanges of mass and energy between plants and the environment and may exert mutual feedback (Baldocchi et al., 2001; Freedman et al., 2001).

Currently, the phenology datasets mainly have three sources: (i) ground-based observations of plant structural changes, (ii) camera-based observations of canopy reflectance (or near-surface remote sensing observations), and (iii) satellite-based observations of land surface reflectance. The ground-based phenology is a traditional but very useful method in phenology studies. For example, the USA National Phenology Network (USA-NPN) has collated observations of first bloom and first leaf of lilac and honeysuckle from the 1960s across the contiguous United States (Schwartz et al., 2012; CONUS, United States territory, not including Hawaii or Alaska; Betancourt et al., 2007; Glynn and Owen, 2015). The USA-NPN was established in part to assemble long-term phenology datasets for a broad array of species across the United States, which can be used to determine the extent to which species, populations, and communities are vulnerable to ongoing and projected future changes in climate (Glynn and Owen, 2015; Schwartz et al., 2012). The camera-based phenology observations such as the PhenoCam network (, last access: 21 March 2021; Richardson et al., 2018; Richardson, 2019) use high-resolution digital cameras to characterize canopy phenology through the color information from the images (Brown et al., 2016; Richardson et al., 2018). The remote sensing has been used to detect vegetation green-up and canopy development (Ganguly et al., 2010; Julien and Sobrino, 2009; Zhang et al., 2003, 2018). While the remote-sensing-based phenology product estimates transition dates from a continuous reflectance time series, it is compared against ground-based event data. The seasonality of ecosystem processes, including biogeochemical fluxes, can similarly be decomposed to identify key transition points and phase durations, which can be determined with high accuracy, and are specific to the processes of interest. As the seasonalities of different processes differ, it can be argued that the interannual trends and responses to environmental forcings can be better described through the fluxes' own temporal characteristics than through correlation to traditional phenological events like bud break or leaf coloration.

Therefore, the objective of this paper is to generate an objective and standardized flux seasonality metrics dataset, which can act as the companion dataset for the FLUXNET product. This study aims to develop a comprehensive framework for studying the seasonality of ecosystem processes systematically with eddy covariance flux measurements including gross primary productivity (GPP), ecosystem respiration (RE), latent heat (LE), and sensible heat (H). As ecosystem and Earth system models are increasingly tested and developed using the very high temporal and increasing spatial coverage of eddy covariance sites, the added information of explicit and standardized flux-specific transition times offers unprecedented opportunity to refine the process representations in models even further (Baldocchi et al., 2001; Falge et al., 2002; Noormets, 2009; Wofsy et al., 1993). The remainder of this paper is organized to the dataset description, the summary of estimating the seasonality metrics from idealized seasonal curves of different fluxes, description of the model performance, uncertainties of the flux seasonality metrics, and conclusions.

2 Data

FLUXNET is a global network of regional networks of eddy covariance sites that measure the exchange of CO2, water vapor, and energy between vegetation and the atmosphere (Baldocchi et al., 2001; Baldocchi, 2008). Recently, harmonized data processing protocols have been developed (Pastorello et al., 2020), and the growing global coverage of these observations has become the de-facto ground-truthing tool for both mechanistic ecosystem models as well as global planetary circulation models (Baldocchi, 2003, 2020). The data include continuous (i.e., gap-filled) measurements of net ecosystem exchange of CO2 (NEE), latent and sensible heat fluxes (LE and H), and microclimate data (air temperature, humidity, wind speed and direction, solar radiation, soil temperature, and soil water content), all at a 30 min time step. Estimates of canopy photosynthesis and ecosystem respiration, derived from the data using an empirical model, are also typically available. The data undergo quality assurance, and missing half-hourly averages are filled using standardized methods to provide continuous data records (Foken and Wichura, 1996; Mauder et al., 2008; Pastorello et al., 2019). The current study uses the FLUXNET2015 Dataset (, last access: 21 March 2021), which includes over 200 sites and around 1500 site years of data (Fig. 1). The gap-filled 30 min data series of fluxes and micrometeorological conditions were aggregated to daily totals. The example sites for each biome were selected based on the following boundary conditions: (1) distinct seasonality of all fluxes; (2) data coverage of observed and high-quality gap-filled data >75 % (defined by variable-specific data quality flags in the FLUXNET database; Reichstein et al., 2005). Therefore, the coverage of different fluxes is different, in which GPP has the highest coverage and H has the lowest coverage. The final dataset included 170 sites and 1035 site years for GPP, 173 sites and 1013 site years for RE, 130 sites and 773 site years for LE, and 110 sites and 446 site years for H. Even though the FLUXNET sites are mainly distributed in the Northern Hemisphere and temperate ecosystems, they still have high spatial and temporal representativeness (Yu et al., 2019).

Figure 1Global distribution of eddy covariance flux sites included in this study, color coded by their International Geosphere-Biosphere Programme (IGBP) biome classes.

Although there is a broad agreement between different flux partitioning approaches (Moffat et al., 2007), and many approaches have converged recently, the current FLUXNET data product still includes a couple of alternative estimates of derived fluxes (RE, GPP). The latest interpretation of respiration fluxes, in particular, is that the night- and daytime estimates may represent different facets of the same phenomenon (Keenan et al., 2019). Different researchers may choose different partitioning approaches for different purposes. Hence, the phenological metrics dataset includes metrics for all biogeochemical fluxes reported in the FLUXNET dataset. Debating the strengths of different partitioning approaches, and indeed the quality of the underlying dataset, is beyond the scope of this study. We assume that the data reported has passed certain thresholds, even though some additional screening has been necessary, and not all site years are of sufficient coverage and quality to estimate the flux seasonality metrics. The metrics are reported both in absolute flux scale (to allow comparisons against commonly reported values) and in relative, normalized scale (to allow comparisons of development rates among different fluxes).

3 Method

3.1 Phenology metrics and uncertainties from flux observation

The seasonal dynamics of the ecosystem fluxes generally have five distinctive phases, which results from the interaction between the inherent biological and ecological processes and the changes in environmental conditions, and reflects the unique functioning of plant community at different stages of the growing season (Gu et al., 2009). The five phases are (1) pre-phase, baseline dormant season flux, before leaf development; (2) flux development period, a rapid increase in flux rate, concurrent or immediately following leaf emergence and expansion; (3) peak flux period, a relatively steady stage in the middle of the growing season; (4) flux recession period, a rapidly declining stage to the baseline; and (5) termination phase, the onset of a new dormant season, following leaf senescence and abscission. Sites with non-standard flux seasonality, where the R2 of model fit was below 0.75, were filtered out during general model fit assessment. Most filtering related to poor quality gap filling with obviously distorted seasonality, but seven sites with multiple peaks during the same year were also excluded from the current study. Southern Hemisphere sites were analyzed by shifting the calendar year cutoffs by 180 d.

Here, we fit the single-pass double-logistic model as first described by Gu et al. (2009) to fit the flux time series. This function exhibits broad structural flexibility and robust convergence, both of which are important for automated processing. The temporal variation in eddy-flux data for an entire growing cycle can be modeled using the function

(1) F m t = f 0 + a 1 1 + e - b 1 t - t 1 c 1 - a 2 1 + e - b 2 t - t 2 c 2 ,

where Fm is the flux value in a given day of year (DOY) t, f0 is the dormant season base flux, and a1 and a2 are parameters about the flux magnitude. Parameters b1, b2, c1, and c2 are related with the transitions or curvature parameters. The model was fit to daily integrated fluxes, following iterative procedures:

  • a.

    Fit Eq. (1) to the flux time series and calculate a predicted value for each DOY.

  • b.

    For each point in the time series, compute the ratio of the observed to predicted flux.

  • c.

    Conduct the Grubbs test (Grubbs, 1969) to identify outliers in the obtained ratios.

  • d.

    If an outlier is detected, remove this outlier and go to step c.

  • e.

    If no outliers are found, remove the data points whose ratios are more than 1 standard deviation (1 s) below the mean ratio.

  • f.

    Fit Eq. (1) to the time series of the daily flux measurements.

The DOYs at which the fitted logistic curve showed characteristic curvature changes were identified with the formula shown in Table 1 derived analytically from the seven parameters of Eq. (1) corresponding to the minimum and maximum values of the first and second derivatives. The first derivative of Eq. (1) is given by

(2) F m t = a 1 c 1 e - ( b 1 t - t 1 ) b 1 ( 1 + e - b 1 t - t 1 ) c 1 + 1 - a 2 c 2 e - ( b 2 t - t 2 ) b 2 ( 1 + e - b 2 t - t 2 ) c 2 + 1 .

The first derivative can be considered as the rate of change of the flux. The maximum of the first derivative occurs early and the minimum late in the growing season (Fig. 2). The day on which the maximal growth rate of each flux occurs is termed the “midpoint of flux development period” (DOYMFD; point B), the day on which the minimal growth rate occurs is termed the “midpoint of flux recession period” (DOYMFR; point E), and the interval between these two transition dates is termed the “length of flux midpoint” (LFM=E-B).

Figure 2An example of the seasonal dynamics of gross primary productivity (GPP), and metrics of transition points of different phases derived from the extremes of the first (Fmt) and second (Fm′′t) derivatives of the fitted logistic function (Eq. 1). For visual clarity, the scales of the first and second derivatives are enhanced 20-fold and 200-fold, respectively (orange and purple lines). The blue line indicates the double-logistic model (Eq. 1) fitted to the observed flux time series (black dots). The slope of the green dash lines indicates the rate of change during the flux development, peak flux, and flux recession period. The phenological transition points are marked with the vertical dashed lines, and the bootstrap estimates of 90 % confidence intervals of these metrics are indicated with the horizontal red error bars about the seven key transition points.


The second derivative of Eq. (1) is given by


The second derivative can be considered as the rate of the growth rate of the flux. The spring and fall maxima of the second derivative mark the start of flux development period (DOYSFD; point A) and end of the flux recession period (DOYEFR; point F), whereas the minima mark the transition from the flux development to peak flux period (point C=DOYEFD=DOYSPF) and the transitions from the peak flux to flux recession period (point D=DOYEPF=DOYSFR). Periods between AC, CD, and DF mark the length of flux development, peak flux, and flux recession periods (LFD, LPF, and LFR, respectively). We also calculated the peak flux (Fmax), date of peak flux (DOYFmax), the rate of the flux development period (RFD;l1), the rate of peak flux period (RPF;l2), and the rate of the flux recession period (RFR;l3). Period AF is the length of active season (LAS).

The fitted daily fluxes observations were normalized to 0–1, and then the rates of change were calculated (namely NRFD, NRPF, and NRFR), which can be used for comparison of development rates among different fluxes.

All the seasonality metrics were summarized in Table 1.

Table 1Seasonality metrics estimated for biogeochemical fluxes including gross primary production, ecosystem respiration, latent heat, and sensible heat.

Note that “–” means units vary for phenology metrics of different fluxes.

Download Print Version | Download XLSX

3.2 Evaluation of parameter quality

3.2.1 Model fit statistics

The fit of Eq. (1) to flux time series was characterized through the coefficient of determination (R2), root mean square error (RMSE), empirical BIAS, and agreement index (AI).

The R2 value of a regression is a measure of the portion of the variance of the dependent variable accounted for by the explanatory variables and characterizes the goodness of fit of the fitted model,

(4) R 2 = t = 1 n ( F o ( t ) - F o ¯ ) 2 - t = 1 n ( F o ( t ) - F m ( t ) ) 2 t = 1 n ( F o ( t ) - F o ¯ ) 2 ,

where Fm and Fo are the predicted and observed values, respectively, Fo¯ is the mean value of the observations, and n is the sample size or the number of days in the year.

The RMSE was estimated as the square root of the mean value of the squared residuals:

(5) RMSE = t = 1 n ( F o ( t ) - F m ( t ) ) 2 n .

The BIAS was calculated as the mean value of the model's residuals:

(6) BIAS = t = 1 n ( F o ( t ) - F m ( t ) ) n .

The agreement index (AI; Willmott, 2013) provides a measure of the relative error in model estimates, combining the information contained in the correlation coefficient (R) and RMSE, and is popular in model assessments and calibration (Gu et al., 2002; Zhou et al., 2016). It is calculated as

(7) AI = 100 - 100 t = 1 n ( F o ( t ) - F m ( t ) ) 2 t = 1 n ( F m ( t ) - F o ¯ + F o ( t ) - F o ¯ 2 .

AI is dimensionless and ranges from 0 (complete disagreement) to 100 (perfect fit). The AI is also sensitive to differences between observed and modeled means (Willmott, 2013). Thus, the AI is well suited for comparing model fits across different biomes and climates.

3.2.2 Uncertainties estimation

The uncertainties in the flux seasonality metrics estimates arise from two sources: (i) the day-to-day variability of fluxes, particularly during the transition periods, that affect the overall goodness of fit of Eq. (1) (see Sect. 3.2.1) and (ii) the consistency of change in climatological drivers during the transition periods that can manifest as early or late cold or warm spells providing conflicting signals to plant development and can affect specific metrics without affecting others. The overall model fit statistics can be used to identify the suitability of different data sources for flux seasonality assessment, but they are not good indicators of the confidence in specific seasonality metrics. Assessing the quality of the underlying flux data is beyond the scope of the current study, and all reported flux values are assumed “true” and the best possible estimates. The uncertainties in the seasonality metrics were estimated similar to Elmore et al. (2012), using Monte Carlo bootstrapping (Efron, 1979). Bootstrapping is a statistical procedure that resamples a single dataset to create many simulated samples. The advantage of the bootstrapping is that parameters can be estimated without assumptions about the normal distribution and using also a small sample size. The distribution of parameter estimates for these bootstrap models provides valuable information about parameter uncertainty and correlation that is free of assumptions about the underlying data distributions. In this study, random uniform sampling with replacement was conducted 500 times for each site year, and the seasonality metrics were estimated for each iteration of the bootstrapped dataset (Elmore et al., 2012; Klosterman et al., 2014). The 5th and 95th percentiles of the 500 bootstrapped phenology metrics estimates were taken as the confidence interval of the mean estimated from the original dataset. The uncertainties of the seven key transition dates are shown in Fig. 2, with red horizontal error bars to indicate the uncertainty intervals.

4 Results

4.1 Model fit statistics

The double-logistic model (Eq. 1) captured the temporal dynamics of widely divergent flux time series (Fig. 3). Although the model fit statistics do not directly translate to the quality of the seasonality metrics estimates (see Sect. 3.2.2), the general fit statistics deserve a brief review. Table 2 shows the overall performance of the fitted model for the different fluxes. The primary explanatory variable behind the fit statistics, as well as the differences between different fluxes, was the range of day-to-day variability in the flux time series. For example, H was generally more variable than LE, both of which are much more variable than RE and GPP, which resulted in lower fit statistics for H and LE than for GPP and RE (Table 2, Fig. 3). The model fit statistics reported in Table 2 were largely consistent in ranking the goodness of fit among biomes. Biomes with well expressed seasonal flux magnitude differences (mixed forest, deciduous broadleaf forest, evergreen needleleaf forest, croplands) exhibited consistently higher fit statistics than biomes with weak seasonality of fluxes (CSH and EBF). The latter also exhibited relatively greater day-to-day variability of fluxes, resulting in lower fit statistics. For GPP, the fit statistics were practically indistinguishable for time series partitioned based on daytime and nighttime partitioning methods (GPP-DT: AI=98.81, R2=0.96; GPP-NT: AI=98.73; R2=0.95), whereas for RE the nighttime partitioning method showed a marginally better fit than the daytime method (AI=98.66 and R2=0.95 versus AI=97.82 and R2=0.92).

Figure 3Examples of the seasonal dynamics of different fluxes for 10 sites representative of different biomes. One biome, open shrubland, was left off because of space limitations on a single page. The blue line indicates the double-logistic model (Eq. 1) fitted to the observed flux time series (black dots). The phenological transition points are marked with the vertical dashed lines, and the bootstrap estimates of 90 % confidence intervals of these metrics are indicated with the horizontal red error bar for corresponding transition points. Note that the confidence intervals are not always symmetrical to the best estimate.


Table 2Model fit statistics to different flux time series in different biomes.

R2: coefficient of determination; RMSE: root mean square error; AI: agreement index; GPP: gross primary production; RE: ecosystem respiration; H: sensible heat flux; LE: latent heat flux; DT: daytime partitioning method; NT: nighttime partitioning method; CORR: corrected H_F_MDS or LE_F_MDS; MDS: marginal distribution sampling; DBF: deciduous broadleaf forests; EBF: evergreen broadleaf forests; ENF: evergreen needle leaf forests; CRO: croplands; OSH: open shrublands; CSH: closed shrublands; GRA: grasslands; MF: mixed forests; SAV: savannas; WSA: woody savannas; WET: permanent wetlands.

Download Print Version | Download XLSX

4.2 Uncertainties of seasonality metrics

The uncertainties of individual flux seasonality metrics (Table 1) were estimated as the 5th and 95th percentiles of 500 Monte Carlo bootstrapping samples ranging from about a week to several weeks, and the uncertainties of phase durations tended to generally be greater than those of individual transition dates. Generally, the uncertainties were the lowest for the phenology metrics of GPP and highest for H (shown as horizontal red lines on Fig. 3). The average uncertainties of transition dates ranged from 6–11 d for GPP, 8 to 14 d for RE, 10 to 15 d for LE, and 15 to 23 d for H. The average uncertainties of duration length ranged from 12–20 d for GPP, 14 to 23 d for RE, 16 to 25 d for LE, and 23 to 32 d for H.

For all fluxes, deciduous broadleaf forest always showed the lowest uncertainties among all biomes. Uncertainties of flux development midpoints were always lower than those of start and end points. Meanwhile, the uncertainties of duration dates are larger than those of transition dates, indicating the compounding effect of the uncertainties of the start and end dates of the active season. The length of the dormant season also affects the uncertainties: the longer the dormant season, the lower the uncertainties.

4.3 Alternative data sources

4.3.1 Daily peak versus daily total fluxes

Even with standardized data sources like FLUXNET2015 Dataset, some user discretion of data aggregation remains, which may affect the reproducibility of different analyses. We mentioned earlier that the FLUXNET2015 Dataset contains GPP and RE estimates from different gap-filling algorithms. In addition, the daily time series can be assembled as either daily integrated fluxes (more common) or as daily peak fluxes, as proposed by Gu et al. (2009), reasoning that daily peaks may be less sensitive to cloudiness and thus better capture the seasonal dynamics of flux capacity. To test this line of argumentation, and identify the best representation of flux seasonality, we started the current study by comparing the model fit statistics for either of these data types. Except for H, all other fluxes had higher model fit statistics with daily integrated than daily peak fluxes (Fig. 4). Although only R2 values are shown, all other fit statistics confirmed the same pattern (Table 2). Therefore, we adopted the daily integrated fluxes for the following phenology metrics dataset generation. Although the day- and nighttime partitioning methods yielded different flux estimates, and different model fits, choosing the best between them is not appropriate at this point. Both partitioning methods have their uses, and their respective merits have not been conclusively proven. The differences between the seasonality metrics of each dataset will be discussed in Sect. 4.2.3.

Figure 4Comparison of the coefficients of determination (R2) of fitted seasonality curves between daily integrated fluxes (x axis) and daily peak fluxes (y axis) for GPP (a, b), RE (c, d), LE (e), and H (f). Both daytime (DT) and nighttime (NT) flux partitioning models are shown for GPP and RE.


Although daily flux totals were identified as the preferred scalar for deriving seasonality metrics from, we will report here the differences between the metrics estimated from the daily peak ecosystem respiration and daily integrated ecosystem respiration data (Fig. 5). This can be viewed as the minimum methodological uncertainty in a best-case scenario in the sense that the correlation between the two sets of metrics was much higher (average R2=0.81; Fig. 5) than between other sources of variability (e.g., the daytime and nighttime partitioning models resulted in seasonality metrics with an average R2=0.27; Fig. 7). The differences between the metric derived from daily integrated and peak flux values are generally smaller than the confidence intervals of individual estimates, except in WSA and SAV. Annual peak flux from daily peak ecosystem respiration exhibited consistently greater values than daily integrals (Fig. 5o), as would be expected.

Figure 5The scatter plots of different phenology metrics from RE-NT daily integrated data and RE-NT daily peak data.


4.3.2 Comparison of different flux partitioning methods

The significance of the assumptions made by partitioning methods to fill flux time series has been emphasized from the perspective of flux integrals (Kruijt et al., 2004). Here we show that the choice of the partitioning model can also affect the seasonality of daily integrated fluxes and thus the seasonality metrics (Fig. 6). We report here the differences between the daytime and nighttime models of flux partitioning as exemplified by the RE and GPP time series, but the lessons apply for all partitioning approaches. Most importantly, mixing of time series filled with different models should not be done.

Figure 6Seasonal dynamics of RE from daytime and nighttime partitioning methods for 11 sites representative of different biomes.


Figure 7The scatter plots of seasonality metrics from RE data using daytime and nighttime partitioning methods.


Ecosystem respiration (RE)

The nighttime and daytime flux partitioning methods can yield similar or dissimilar daily RE, and sometimes even the seasonality can diverge significantly between them (e.g., US-PFa 2014 in Fig. 6).

Figure 8Seasonal dynamics of GPP from daytime and nighttime partitioning method for 11 sites representative of different biomes.


Figure 9The scatter plots of different phenology metrics from GPP data using daytime and nighttime partitioning methods.


The differences in flux seasonality metrics based on the partitioning method (Fig. 7) are sizable, non-systematic, and generally greater for season length metrics than transition date metrics (because the season length is determined by two transition dates, each of which is subject to deviation among methods). Interestingly, the seasonality metrics of flux development period (Fig. 7a, d, e, f) are more consistent than those of the peak flux period (Fig. 7g, h, i) and flux recession period (Fig. 7b, j, k, l). A more detailed analysis of the performance of the consistency of the different partitioning methods is the subject of future studies, and we may find answers more from information pertaining more to the accuracy of the flux estimates than the seasonality estimates.

Gross primary production (GPP)

The seasonality of GPP time series differed less between the two flux partitioning methods than did RE (Fig. 8). It is also obvious that the seasonalities of RE and GPP for the same sites differed significantly in terms of the seasonal timing, symmetry, peak duration, and other aspects. A more detailed assessment of these differences is the subject of a forthcoming analysis (Yang and Noormets, 2021).

As suggested by the extent of overlap in the flux estimates on Fig. 8, the seasonality metrics of GPP from the two partitioning methods were also more consistent compared to RE (Fig. 9). The scatter was distributed around the 1:1 line and with little bias. However, similar to RE, the metrics of the flux development period of GPP were also more consistent among the partitioning methods than those of peak flux period and flux recession period. Yet, the rate of flux development was more variable between the methods than the rate of flux recession (Fig. 9p, r, s, u). And like for ER, the phase duration metrics of GPP were more variable than the transition dates.

5 Data availability

This database can be obtained at the US Department of Energy's (DOE) Environmental Systems Science Data Infrastructure for a Virtual Ecosystem (ESS-DIVE,; Yang and Noormets, 2020).

6 Significance

The Flux Seasonality Metrics Database (FSMD) presented here summarizes important latent features contained in the land surface biogeochemical fluxes and is likely to allow novel insights into the functioning of the biosphere and assist in the development and validation of novel functionality in Earth system models. Improving the predictive capabilities of ecosystem biogeochemistry models on interannual and decadal scales remains a challenge, and variability in the seasonality of different fluxes has been recognized as a key uncertainty. Importantly, the seasonality of GPP in models is often forced to match observations with arbitrary coefficients (Straube et al., 2018), as the divergence of LE and GPP seasonality is not captured in common LAI-driven models (Wu et al., 2017; Restrepo-Coupe et al., 2017). By developing process-specific seasonality references, explicit validation of these fluxes becomes possible. In addition, the ability to discern shifts in seasonality from those in flux capacity or vegetation structure can also be important in correcting the attributions of observed changes.

A standardized flux seasonality metrics dataset can also support other seasonality assessment tools, like near-surface (e.g., PhenoCam network; Richardson et al., 2012) and remote optical sensors (Broich et al., 2015; White and Nemani, 2006; Ganguly et al., 2010; Gamon et al., 2016). However, as these methods purport to infer GPP from the greenness information, they are vulnerable to the same lags between leaf development, LE, and GPP that have undermined the model development mentioned above. In addition, pixel heterogeneity and clouds can reduce the potential of remote sensing approaches.

The FSMD was designed to capture and depict the seasonality of different ecosystem processes. FSMD will be updated within 6 months of each major release of the FLUXNET database, making it possible to quantify the differences and similarities between different ecosystem processes in their responses to changes in climatic conditions. Some potential applications of this dataset have been mentioned before, and there are likely many others. We anticipate that the FSMD will stimulate new research in global change and Earth science disciplines where land–atmosphere exchange dynamics play a central role.

Author contributions

LY and AN designed the research; LY implemented the research and wrote the manuscript; AN edited and revised the manuscript.

Competing interests

The authors declare that they have no conflict of interest.


This work was supported, in part, by the US Department of Energy's Office of Science through the AmeriFlux Management Project (award no. M1800063) and the USDA McIntire-Stennis Cooperative Forestry Program (project no. 121209). LY was supported by a Texas A&M AgriLife Research Strategic Initiative Assistantship. We thank everybody who has contributed data and effort to producing the FLUXNET database.

Review statement

This paper was edited by David Carlson and reviewed by Geoffrey Henebry and one anonymous referee.


Baldocchi, D.: “Breathing” of the terrestrial biosphere: lessons learned from a global network of carbon dioxide flux measurement systems, Aust. J. Bot., 56, 1–26,, 2008. 

Baldocchi, D. D.: Assessing the eddy covariance technique for evaluating carbon dioxide exchange rates of ecosystems: past, present and future, Glob. Change Biol., 9, 479–492,, 2003. 

Baldocchi, D. D.: How eddy covariance flux measurements have contributed to our understanding of Global Change Biology, Glob. Change Biol., 26, 242–260,, 2020. 

Baldocchi, D., Falge, E., Gu, L., Olson, R., Hollinger, D., Running, S., Anthoni, P., Bernhofer, C., Davis, K., Evans, R., Fuentes, J., Goldstein, A., Katul, G., Law, B., Lee, X., Malhi, Y., Meyers, T., Munger, W., Oechel, W., Paw, K. T., Pilegaard, K., Schmid, H. P., Valentini, R., Verma, S., Vesala, T., Wilson, K., and Wofsy, S.: FLUXNET: A New Tool to Study the Temporal and Spatial Variability of Ecosystem–Scale Carbon Dioxide, Water Vapor, and Energy Flux Densities, B. Am. Meteorol. Soc., 82, 2415–2434,<2415:FANTTS>2.3.CO;2, 2001. 

Betancourt, J. L., Schwartz, M. D., Breshears, D. D., Brewer, C. A., Frazer, G., Gross, J. E., Mazer, S. J., Reed, B. C., and Wilson, B. E.: Evolving plans for the USA National Phenology Network, Eos, 88, 211–211,, 2007. 

Broich, M., Huete, A., Paget, M., Ma, X., Tulbure, M., Coupe, N. R., Evans, B., Beringer, J., Devadas, R., Davies, K., and Held, A.: A spatially explicit land surface phenology data product for science, monitoring and natural resources management applications, Environ. Modell. Softw., 64, 191–204,, 2015. 

Brown, T. B., Hultine, K. R., Steltzer, H., Denny, E. G., Denslow, M. W., Granados, J., Henderson, S., Moore, D., Nagai, S., SanClements, M., Sánchez-Azofeifa, A., Sonnentag, O., Tazik, D., and Richardson, A. D.: Using phenocams to monitor our changing Earth: toward a global phenocam network, Front. Ecol. Environ., 14, 84–93,, 2016. 

Efron, B.: Bootstrap Methods: Another Look at the Jackknife, Ann. Stat., 7, 1–26,, 1979. 

Elmore, A. J., Guinn, S. M., Minsley, B. J., and Richardson, A. D.: Landscape controls on the timing of spring, autumn, and growing season length in mid-Atlantic forests, Glob. Change Biol., 18, 656–674,, 2012. 

Falge, E., Baldocchi, D., Tenhunen, J., Aubinet, M., Bakwin, P., Berbigier, P., Bernhofer, C., Burba, G., Clement, R., Davis, K. J., Elbers, J. A., Goldstein, A. H., Grelle, A., Granier, A., Guðmundsson, J., Hollinger, D., Kowalski, A. S., Katul, G., Law, B. E., Malhi, Y., Meyers, T., Monson, R. K., Munger, J. W., Oechel, W., Paw U, K. T., Pilegaard, K., Rannik, Ü., Rebmann, C., Suyker, A., Valentini, R., Wilson, K., and Wofsy, S.: Seasonality of ecosystem respiration and gross primary production as derived from FLUXNET measurements, Agr. Forest Meteorol., 113, 53–74,, 2002. 

Fitzjarrald, D. R., Acevedo, O. C., and Moore, K. E.: Climatic Consequences of Leaf Presence in the Eastern United States, J. Climate, 14, 598–614,<0598:CCOLPI>2.0.CO;2, 2001. 

Foken, T. and Wichura, B.: Tools for quality assessment of surface-based flux measurements, Agr. Forest Meteorol., 78, 83–105,, 1996. 

Freedman, J. M., Fitzjarrald, D. R., Moore, K. E., and Sakai, R. K.: Boundary Layer Clouds and Vegetation–Atmosphere Feedbacks, J. Climate, 14, 180–197,<0180:BLCAVA>2.0.CO;2, 2001. 

Gamon, J. A., Huemmrich, K. F., Wong, C. Y., Ensminger, I., Garrity, S., Hollinger, D. Y., Noormets, A., and Penuelas, J.: A remotely sensed pigment index reveals photosynthetic phenology in evergreen conifers, P. Natl. Acad. Sci. USA, 113, 13087–13092,, 2016. 

Ganguly, S., Friedl, M. A., Tan, B., Zhang, X., and Verma, M.: Land surface phenology from MODIS: Characterization of the Collection 5 global land cover dynamics product, Remote Sens. Environ., 114, 1805–1816,, 2010. 

Glynn, P. D. and Owen, T. W.: Review of the USA National Phenology Network, USGS Numbered Series, U.S. Geological Survey, Reston, VA, 26 pp.,, 2015. 

Grubbs, F. E.: Procedures for Detecting Outlying Observations in Samples, Technometrics, 11, 1–21,, 1969. 

Gu, L., Baldocchi, D., Verma, S. B., Black, T. A., Vesala, T., Falge, E. M., and Dowty, P. R.: Advantages of diffuse radiation for terrestrial ecosystem productivity, J. Geophys. Res.-Atmos., 107, ACL 2-1-ACL 2–23,, 2002. 

Gu, L., Post, W. M., Baldocchi, D. D., Black, T. A., Suyker, A. E., Verma, S. B., Vesala, T., and Wofsy, S. C.: Characterizing the Seasonal Dynamics of Plant Community Photosynthesis Across a Range of Vegetation Types, in: Phenology of Ecosystem Processes, 1st Edn., Springer, New York, NY, 35–58,, 2009. 

Julien, Y. and Sobrino, J. A.: Global land surface phenology trends from GIMMS database, Int. J. Remote Sens., 30, 3495–3513,, 2009. 

Keenan, T. F., Migliavacca, M., Papale, D., Baldocchi, D., Reichstein, M., Torn, M., and Wutzler, T.: Widespread inhibition of daytime ecosystem respiration, Nat. Ecol. Evol., 3, 407–415,, 2019. 

Klosterman, S. T., Hufkens, K., Gray, J. M., Melaas, E., Sonnentag, O., Lavine, I., Mitchell, L., Norman, R., Friedl, M. A., and Richardson, A. D.: Evaluating remote sensing of deciduous forest phenology at multiple spatial scales using PhenoCam imagery, Biogeosciences, 11, 4305–4320,, 2014. 

Kruijt, B., Elbers, J. A., von Randow, C., Araújo, A. C., Oliveira, P. J., Culf, A., Manzi, A. O., Nobre, A. D., Kabat, P., and Moors, E. J.: The Robustness of Eddy Correlation Fluxes for Amazon Rain Forest Conditions, Ecol. Appl., 14, 101–113,, 2004. 

Lieth, H.: Phenology and Seasonality Modeling, Springer-Verlag Berlin, Heidelberg, New York,, 1974. 

Moffat, A. M., Papale, D., Reichstein, M., Hollinger, D. Y., Richardson, A. D., Barr, A. G., Beckstein, C., Braswell, B. H., Churkina, G., Desai, A. R., Falge, E., Gove, J. H., Heimann, M., Hui, D., Jarvis, A. J., Kattge, J., Noormets, A., and Stauch, V. J.: Comprehensive comparison of gap-filling techniques for eddy covariance net carbon fluxes, Agr. Forest Meteorol., 147, 209–232,, 2007. 

Mauder, M., Foken, T., Clement, R., Elbers, J. A., Eugster, W., Grünwald, T., Heusinkveld, B., and Kolle, O.: Quality control of CarboEurope flux data – Part 2: Inter-comparison of eddy-covariance software, Biogeosciences, 5, 451–462,, 2008. 

Noormets, A.: Phenology of Ecosystem Processes, 1st Edn., Springer, New York, NY,, 2009. 

Noormets, A., Chen, J., Gu, L., and Desai, A.: The Phenology of Gross Ecosystem Productivity and Ecosystem Respiration in Temperate Hardwood and Conifer Chronosequences, in: Phenology of Ecosystem Processes: Applications in Global Change Research, edited by: Noormets, A., Springer, New York, 59–85,, 2009. 

Pastorello, G., Trotta, C., Ribeca, A., Elbashandy, A., Barr, A., and Papale, D.: ONEFlux: Open Network-Enabled Flux processing pipeline, AmeriFlux Management Project, European Ecosystem Fluxes Database, available at: (last access: 21 March 2021), 2019. 

Pastorello, G., Trotta, C., Canfora, E., Chu, H., Christianson, D., Cheah, Y. W., Poindexter, C., Chen, J., Elbashandy, A., Humphrey, M., Isaac, P., Polidori, D., Ribeca, A., van Ingen, C., Zhang, L., Amiro, B., Ammann, C., Arain, M. A., Ardo, J., Arkebauer, T., Arndt, S. K., Arriga, N., Aubinet, M., Aurela, M., Baldocchi, D., Barr, A., Beamesderfer, E., Marchesini, L. B., Bergeron, O., Beringer, J., Bernhofer, C., Berveiller, D., Billesbach, D., Black, T. A., Blanken, P. D., Bohrer, G., Boike, J., Bolstad, P. V., Bonal, D., Bonnefond, J. M., Bowling, D. R., Bracho, R., Brodeur, J., Brummer, C., Buchmann, N., Burban, B., Burns, S. P., Buysse, P., Cale, P., Cavagna, M., Cellier, P., Chen, S., Chini, I., Christensen, T. R., Cleverly, J., Collalti, A., Consalvo, C., Cook, B. D., Cook, D., Coursolle, C., Cremonese, E., Curtis, P. S., D'Andrea, E., da Rocha, H., Dai, X., Davis, K. J., De Cinti, B., de Grandcourt, A., De Ligne, A., De Oliveira, R. C., Delpierre, N., Desai, A. R., Di Bella, C. M., di Tommasi, P., Dolman, H., Domingo, F., Dong, G., Dore, S., Duce, P., Dufrene, E., Dunn, A., Dusek, J., Eamus, D., Eichelmann, U., ElKhidir, H. A. M., Eugster, W., Ewenz, C. M., Ewers, B., Famulari, D., Fares, S., Feigenwinter, I., Feitz, A., Fensholt, R., Filippa, G., Fischer, M., Frank, J., Galvagno, M., Gharun, M., Gianelle, D., Gielen, B., Gioli, B., Gitelson, A., Goded, I., Goeckede, M., Goldstein, A. H., Gough, C. M., Goulden, M. L., Graf, A., Griebel, A., Gruening, C., Grunwald, T., Hammerle, A., Han, S., Han, X., Hansen, B. U., Hanson, C., Hatakka, J., He, Y., Hehn, M., Heinesch, B., Hinko-Najera, N., Hortnagl, L., Hutley, L., Ibrom, A., Ikawa, H., Jackowicz-Korczynski, M., Janous, D., Jans, W., Jassal, R., Jiang, S., Kato, T., Khomik, M., Klatt, J., Knohl, A., Knox, S., Kobayashi, H., Koerber, G., Kolle, O., Kosugi, Y., Kotani, A., Kowalski, A., Kruijt, B., Kurbatova, J., Kutsch, W. L., Kwon, H., Launiainen, S., Laurila, T., Law, B., Leuning, R., Li, Y., Liddell, M., Limousin, J. M., Lion, M., Liska, A. J., Lohila, A., Lopez-Ballesteros, A., Lopez-Blanco, E., Loubet, B., Loustau, D., Lucas-Moffat, A., Luers, J., Ma, S., Macfarlane, C., Magliulo, V., Maier, R., Mammarella, I., Manca, G., Marcolla, B., Margolis, H. A., Marras, S., Massman, W., Mastepanov, M., Matamala, R., Matthes, J. H., Mazzenga, F., McCaughey, H., McHugh, I., McMillan, A. M. S., Merbold, L., Meyer, W., Meyers, T., Miller, S. D., Minerbi, S., Moderow, U., Monson, R. K., Montagnani, L., Moore, C. E., Moors, E., Moreaux, V., Moureaux, C., Munger, J. W., Nakai, T., Neirynck, J., Nesic, Z., Nicolini, G., Noormets, A., Northwood, M., Nosetto, M., Nouvellon, Y., Novick, K., Oechel, W., Olesen, J. E., Ourcival, J. M., Papuga, S. A., Parmentier, F. J., Paul-Limoges, E., Pavelka, M., Peichl, M., Pendall, E., Phillips, R. P., Pilegaard, K., Pirk, N., Posse, G., Powell, T., Prasse, H., Prober, S. M., Rambal, S., Rannik, U., Raz-Yaseef, N., Reed, D., de Dios, V. R., Restrepo-Coupe, N., Reverter, B. R., Roland, M., Sabbatini, S., Sachs, T., Saleska, S. R., Sanchez-Canete, E. P., Sanchez-Mejia, Z. M., Schmid, H. P., Schmidt, M., Schneider, K., Schrader, F., Schroder, I., Scott, R. L., Sedlak, P., Serrano-Ortiz, P., Shao, C., Shi, P., Shironya, I., Siebicke, L., Sigut, L., Silberstein, R., Sirca, C., Spano, D., Steinbrecher, R., Stevens, R. M., Sturtevant, C., Suyker, A., Tagesson, T., Takanashi, S., Tang, Y., Tapper, N., Thom, J., Tiedemann, F., Tomassucci, M., Tuovinen, J. P., Urbanski, S., Valentini, R., van der Molen, M., van Gorsel, E., van Huissteden, K., Varlagin, A., Verfaillie, J., Vesala, T., Vincke, C., Vitale, D., Vygodskaya, N., Walker, J. P., Walter-Shea, E., Wang, H., Weber, R., Westermann, S., Wille, C., Wofsy, S., Wohlfahrt, G., Wolf, S., Woodgate, W., Li, Y., Zampedri, R., Zhang, J., Zhou, G., Zona, D., Agarwal, D., Biraud, S., Torn, M., and Papale, D.: The FLUXNET2015 dataset and the ONEFlux processing pipeline for eddy covariance data, Sci. Data, 7, 225,, 2020. 

Piao, S., Liu, Q., Chen, A., Janssens, I. A., Fu, Y., Dai, J., Liu, L., Lian, X., Shen, M., and Zhu, X.: Plant phenology and global climate change: Current progresses and challenges, Glob. Change Biol., 25, 1922–1940,, 2019. 

Post, E. and Stenseth, N. C.: Climatic Variability, Plant Phenology, and Northern Ungulates, Ecology, 80, 1322–1339,[1322:Cvppan]2.0.Co;2, 1999. 

Reichstein, M., Falge, E., Baldocchi, D., Papale, D., Aubinet, M., Berbigier, P., Bernhofer, C., Buchmann, N., Gilmanov, T., Granier, A., Grunwald, T., Havrankova, K., Ilvesniemi, H., Janous, D., Knohl, A., Laurila, T., Lohila, A., Loustau, D., Matteucci, G., Meyers, T., Miglietta, F., Ourcival, J.-M., Pumpanen, J., Rambal, S., Rotenberg, E., Sanz, M., Tenhunen, J., Seufert, G., Vaccari, F., Vesala, T., Yakir, D., and Valentini, R.: On the separation of net ecosystem exchange into assimilation and ecosystem respiration: review and improved algorithm, Glob. Change Biol., 11, 1424–1439,, 2005. 

Restrepo-Coupe, N., Levine, N. M., Christoffersen, B. O., Albert, L. P., Wu, J., Costa, M. H., Galbraith, D., Imbuzeiro, H., Martins, G., da Araujo, A. C., Malhi, Y. S., Zeng, X., Moorcroft, P., and Saleska, S. R.: Do dynamic global vegetation models capture the seasonality of carbon fluxes in the Amazon basin? A data-model intercomparison, Glob. Change Biol., 23, 191–208,, 2017. 

Richardson, A. D., Anderson, R. S., Arain, M. A., Barr, A. G., Bohrer, G., Chen, G., Chen, J. M., Ciais, P., Davis, K. J., Desai, A. R., Dietze, M. C., Dragoni, D., Garrity, S. R., Gough, C. M., Grant, R., Hollinger, D. Y., Margolis, H. A., McCaughey, H., Migliavacca, M., Monson, R. K., Munger, J. W., Poulter, B., Raczka, B. M., Ricciuto, D. M., Sahoo, A. K., Schaefer, K., Tian, H., Vargas, R., Verbeeck, H., Xiao, J., and Xue, Y.: Terrestrial biosphere models need better representation of vegetation phenology: results from the North American Carbon Program Site Synthesis, Glob. Change Biol., 18, 566–584,, 2012. 

Richardson, A. D., Hufkens, K., Milliman, T., Aubrecht, D. M., Chen, M., Gray, J. M., Johnston, M. R., Keenan, T. F., Klosterman, S. T., Kosmala, M., Melaas, E. K., Friedl, M. A., and Frolking, S.: Tracking vegetation phenology across diverse North American biomes using PhenoCam imagery, Sci. Data, 5, 180028,, 2018. 

Richardson, A. D.: Tracking seasonal rhythms of plants in diverse ecosystems with digital camera imagery, New Phytol., 222, 1742–1750,, 2019. 

Schwartz, M. D.: Phenology: An Integrative Environmental Science, 1st Edn., Tasks for Vegetation Science, 39, Springer, the Netherlands,, 2003. 

Schwartz, M. D. and Crawford, T. M.: Detecting Energy-Balance Modifications at the Onset of Spring, Phys. Geogr., 22, 394–409,, 2013. 

Schwartz, M. D., Betancourt, J. L., and Weltzin, J. F.: From Caprio's lilacs to the USA National Phenology Network, Front. Ecol. Environ., 10, 324–327,, 2012. 

Straube, J. R., Chen, M., Parton, W. J., Asso, S., Liu, Y.-A., Ojima, D. S., and Gao, W.: Development of the DayCent-Photo model and integration of variable photosynthetic capacity, Front. Earth Sci., 12, 765–778,, 2018. 

Weltzin, J. F., Betancourt, J. L., Cook, B. I., Crimmins, T. M., Enquist, C. A. F., Gerst, M. D., Gross, J. E., Henebry, G. M., Hufft, R. A., Kenney, M. A., Kimball, J. S., Reed, B. C., and Running, S. W.: Seasonality of biological and physical systems as indicators of climatic variation and change, Climatic Change, 163, 1755–1771,, 2020. 

White, M. A. and Nemani, R. R.: Real-time monitoring and short-term forecasting of land surface phenology, Remote Sens. Environ., 104, 43–49,, 2006. 

Willmott, C. J.: On the Validation of Models, Phys. Geogr., 2, 184–194,, 2013. 

Wofsy, S. C., Goulden, M. L., Munger, J. W., Fan, S. M., Bakwin, P. S., Daube, B. C., Bassow, S. L., and Bazzaz, F. A.: Net Exchange of CO2 in a Mid-Latitude Forest, Science, 260, 1314–1317,, 1993. 

Wu, J., Serbin, S. P., Xu, X., Albert, L. P., Chen, M., Meng, R., Saleska, S. R., and Rogers, A.: The phenology of leaf quality and its within-canopy variation is essential for accurate modeling of photosynthesis in tropical evergreen forests, Glob. Change Biol., 23, 4814–4827,, 2017. 

Yang, L. and Noormets, A.: Flux Seasonality Metrics Database: A companion dataset for FLUXNET annual product, ESS-DIVE,, 2020. 

Yang, L. and Noormets, A.: Asynchrony of the seasonal dynamics of gross primary productivity (GPP) and ecosystem respiration (RE): Asynchrony of GPP and RE, Glob. Change Biol., submitted, 2021.  

Yu, R., Ruddell, B. L., Kang, M., Kim, J., and Childers, D.: Anticipating global terrestrial ecosystem state change using FLUXNET, Glob. Change Biol., 25, 2352–2367,, 2019. 

Zhang, X., Liu, L., Liu, Y., Jayavelu, S., Wang, J., Moon, M., Henebry, G. M., Friedl, M. A., and Schaaf, C. B.: Generation and evaluation of the VIIRS land surface phenology product, Remote Sens. Environ., 216, 212–229,, 2018. 

Zhang, X. Y., Friedl, M. A., Schaaf, C. B., Strahler, A. H., Hodges, J. C. F., Gao, F., Reed, B. C., and Huete, A.: Monitoring vegetation phenology using MODIS, Remote Sens. Environ., 84, 471–475,, 2003. 

Zhou, Y., Wu, X., Ju, W., Chen, J. M., Wang, S., Wang, H., Yuan, W., Andrew Black, T., Jassal, R., Ibrom, A., Han, S., Yan, J., Margolis, H., Roupsard, O., Li, Y., Zhao, F., Kiely, G., Starr, G., Pavelka, M., Montagnani, L., Wohlfahrt, G., D'Odorico, P., Cook, D., Arain, M. A., Bonal, D., Beringer, J., Blanken, P. D., Loubet, B., Leclerc, M. Y., Matteucci, G., Nagy, Z., Olejnik, J., Paw U, K. T., and Varlagin, A.: Global parameterization and validation of a two-leaf light use efficiency model for predicting gross primary production across FLUXNET sites, J. Geophys. Res.-Biogeo., 121, 1045–1072,, 2016. 


The requested paper has a corresponding corrigendum published. Please read the corrigendum first before downloading the article.

Short summary
We present a flux seasonality metrics database (FSMD) depicting a set of standardized metrics of ecosystem biogeochemical fluxes of CO2, water, and energy, including transition dates, phase lengths, and rates of change with uncertainty estimates. FSMD allows assessment of spatial and temporal patterns in developmental dynamics, validation of novel aspects of phenology product, and process models. It is calculated from FLUXNET2015 data product and will be updated with new FLUXNET data releases.