the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Global datasets of hourly carbon and water fluxes simulated using a satellite-based process model with dynamic parameterizations
Jiye Leng
Jing M. Chen
Wenyu Li
Xiangzhong Luo
Mingzhu Xu
Rong Wang
Cheryl Rogers
Bolun Li
Yulin Yan
Abstract. Diagnostic terrestrial biosphere models (TBMs) forced by remote sensing observations have been a principal tool to provide benchmarks on global gross primary productivity (GPP) and evapotranspiration (ET). However, these models often estimate GPP and ET at coarse daily or monthly step, hindering analysis of ecosystems dynamics at the diurnal (hourly) scales, and prescribe some essential parameters (i.e., the Ball-Berry slope (m) and the maximum carboxylation rate at 25 °C (Vcmax25)) as constant, inducing uncertainties in the estimates of GPP and ET. In this study, we present hourly estimation of global GPP and ET datasets at a 0.25° resolution from 2001 to 2020 simulated with a widely used diagnostic TBM – Biosphere-atmosphere Exchange Process Simulator (BEPS). We employed eddy covariance observations and machine learning approaches to derive and upscale the seasonally varied m and Vcmax25 for carbon and water fluxes. The estimated hourly GPP and ET are validated against flux observation, remote sensing, and machine learning-based estimates across multiple spatial and temporal scales. The correlation coefficients (R2) and slopes between hourly tower-measured and modeled fluxes are: R2 = 0.83, regression slope = 0.92 for GPP and, R2 = 0.72, regression slope = 1.04 for ET. At the global scale, we estimated a global mean GPP of 137.78 ± 3.22 Pg C yr-1 (mean ± 1 SD) with a positive trend of 0.53 Pg C yr-2 (p < 0.001), and ET of 89.03 ± 0.82 × 103 km3 yr-1 with a slight positive trend of 0.10 × 103 km3 yr-2 (p < 0.001) from 2001 to 2020. The spatial pattern of our estimates agrees well with other products, with R2 = 0.77–0.85 and R2 = 0.74–0.90 for GPP and ET, respectively. Overall, this new global hourly dataset serves as a `handshake’ among process-based models, remote sensing, and the eddy covariance flux network, providing a reliable long-term estimate of global GPP and ET with diurnal patterns and facilitating studies related to ecosystem functional properties, global carbon, and water cycles. The hourly and accumulated daily GPP and ET estimates are available at https://doi.org/10.5281/zenodo.8240492 (Leng et al., 2023).
- Preprint
(3031 KB) - Metadata XML
-
Supplement
(366 KB) - BibTeX
- EndNote
Jiye Leng et al.
Status: final response (author comments only)
-
RC1: 'Comment on essd-2023-328', Anonymous Referee #1, 14 Sep 2023
This study reported a global product of hourly GPP and ET from 2001 to 2020. In general, this study is important and valuable because none of previous GPP and ET products focused on the hourly scale. However, the authors need introduce the methods more clear especially for the parameterization, which is quite important to judge this dataset is reliable. Besides, I also hope the authors can show the more details of parameters’ distribution and validations. Please see my detailed comments below. Vegetation transpiration is simulated by the supplementary Equ. 8. Did the authors set the same Rn-G to shaded and sunlit leaf? Obviously, there are large differences of latent heat for shaded and sunlit leaves. Line 115: why did you use GLOBMAP LAI data? The spatial resolution of GLOBMAP is quite coarse. Line 150: it is very important to know the details of optimization algorithm. It seems the parameterization method has not been published. Although the authors mentioned the supplementary materials, but I still did not get how the authors optimize the model parameters. Especially, you mentioned the parameters were optimized for each month at each site-year. From 2.3.1 and 2.3.2, I assumed the authors first inversed two model parameters m and Vcmax at eddy covariance sites, and then used machine-learning method to generate global gridded dataset of m and Vcmax, and finally, BEPS model was run based on the gridded m and Vcmax to estimate global GPP? if it is so, why the authors did not just upscale GPP from towers to global scale, just like Jung et al. 2009. Besides, the authors did not show the global patterns of m and Vcmax at all, and we cannot judge if their distributions are reliable. And it is also necessary to show the performance of machine learning method to simulate m and Vcmax at eddy covariance towers, which is quite important than GPP. In addition, I am curious that if the authors did not use gridded parameters, and just used site-based inversed parameters to simulate global GPP, how is the performance? Fig. 4 showed the better performance of hourly simulations than daily simulations. Is it possible? I am curious how the authors examine the performance of hourly simulations, and if the authors included all simulations of night and daytime together, which will result in a false high correlation. As this study aimed to produce an hourly GPP and ET, if there are large difference of parameters diurnal scale. Figure 9: the comparison does not make sense as the products used different LAI datasets, and the different trends basically depends on the trend of LAI datasets.
Citation: https://doi.org/10.5194/essd-2023-328-RC1 -
RC2: 'Comment on essd-2023-328', Anonymous Referee #2, 26 Sep 2023
This study presented hourly estimation of global GPP and ET dataset at a 0.25° resolution from 2001 to 2020. This product could be of great importance in understanding the diurnal ecosystem functionality, and promote monitoring of ecosystem responses to extreme climate events. However, there are several concerns regarding the methods explanation and result presentation that could dampen the science of this study.
Methods: 1. need particular explanation on the newly revised BEPS regarding how exactly the hourly GPP and ET were simulated but not from the previous version. If BEPS was designed to simulate the hourly products, were if just because the hourly inputs were not available before?
- Need a better presentation and validation on the newly optimized key photosynthesis and stomatal conductance model parameters (i.e., 𝑉cmax and and 𝑚). Was it only revised for the flux sites (as presented in Figure 2)? How were they interpolated into the global scales and what are their uncertainties? How about the spatial and temporal variations of these parameters globally? Also, what will be the differences in the simulated GPP and ET between the new dynamic parameters and original fixed parameters in terms of accuracy and spatial pattern?
- Figure 1. The bottom-up order is quite counter-intuitive to readers. Suggest using top-down order.
Results: 1. Figure 4, suggest adding the label of 1 in the slope subplots as a reference of good fitting.
- Need better presentations on the diurnal patterns of GPP and ET for different vegetation function types. For example, providing hourly curves for different vegetation function types and validated against flux site observations. It still unclear whether these products can capture the diurnal variations of GPP and ET.
Citation: https://doi.org/10.5194/essd-2023-328-RC2
Jiye Leng et al.
Data sets
Global datasets of hourly carbon and water fluxes simulated using a satellite-based process model with dynamic parameterizations Jiye Leng, Jing M. Chen, Wenyu Li, Xiangzhong Luo, Mingzhu Xu, Jane Liu, Rong Wang, Cheryl Rogers, Bolun Li, Yulin Yan https://doi.org/10.5281/zenodo.8240492
Model code and software
Biosphere-atmosphere Exchange Process Simulator (BEPS) Jiye Leng, Jing M. Chen, Xiangzhong Luo, Jane Liu https://github.com/JChen-UToronto/BEPS_hourly_site
Jiye Leng et al.
Viewed
HTML | XML | Total | Supplement | BibTeX | EndNote | |
---|---|---|---|---|---|---|
389 | 114 | 15 | 518 | 29 | 7 | 13 |
- HTML: 389
- PDF: 114
- XML: 15
- Total: 518
- Supplement: 29
- BibTeX: 7
- EndNote: 13
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1