Next Article in Journal
Rural Roads Are Paving the Way for Land-Use Intensification in the Uplands of Laos
Previous Article in Journal
Effect of Organic Fertilizer on Soil Bacteria in Maize Fields
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulation of Gross Primary Productivity Using Multiple Light Use Efficiency Models

1
School of Environment and Municipal Engineering, Lanzhou Jiaotong University, Lanzhou 730000, China
2
Gansu Academy of Eco-Environmental Sciences, Lanzhou 730000, China
3
Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou 730000, China
*
Author to whom correspondence should be addressed.
Land 2021, 10(3), 329; https://doi.org/10.3390/land10030329
Submission received: 6 February 2021 / Revised: 13 March 2021 / Accepted: 19 March 2021 / Published: 23 March 2021
(This article belongs to the Section Land Systems and Global Change)

Abstract

:
Gross primary productivity (GPP) is the most basic variable in a carbon cycle study that determines the carbon that enters the ecosystem. The remote sensing-based light use efficiency (LUE) model is one of the primary tools that is currently used to estimate the GPP at the regional scale. Many remote sensing-based GPP models have been developed in the last several decades, and these models have been well evaluated at some sites. However, an accurate estimation of the GPP remains challenging work using LUE models because of uncertainties in the model caused by model parameters, model forcing, and vegetation spatial heterogeneity. In this study, five widely used LUE models, Glo-PEM, VPM, EC-LUE, the MODIS GPP algorithm, and C-fix, were selected to simulate the GPP of the Heihe River Basin forced using in situ measurements. A multiple-model averaging method, Bayesian model averaging (BMA), was used to combine the five models to obtain a more reliable GPP estimation. The BMA was trained using carbon flux data from five eddy covariance towers located at dominant vegetation types in the study area. Generally, the BMA method performed better than any single LUE model. From the case study in the study area, it is indicated that the trained BMA is an efficient method to combine multiple LUE models and can improve the GPP simulation accuracy.

1. Introduction

Gross primary productivity (GPP) is the rate of carbon fixation through vegetation photosynthesis. GPP is a key measure of the carbon mass flux in carbon cycle studies [1]. Accurately estimating the GPP at the regional scale is of importance to understand the carbon cycle mechanism in a specific region. The light use efficiency (LUE) model, which uses remote sensing observations as inputs, is one the most effective tools to estimate regional GPP [2]. In the LUE model, GPP is a function of the amount of photosynthetically active radiation (PAR), the fraction of PAR absorbed by the plant canopy (fAPAR), and environmental factors (e.g., temperature, soil moisture, and vapor pressure deficit) [3]. Additionally, the GPP can be obtained at the ecosystem scale by using the eddy covariance (EC) technique. A lot of networks have been established using the EC technique to monitor the carbon flux between the vegetation and the atmosphere [4]. The EC measurements provide in-situ data to validate and calibrate the LUE models, and this enhances the development of LUE models.
In recent years, many LUE models have been developed, demonstrating different performances in different regions [5,6,7,8]. To reduce the uncertainties in regional GPP estimation, the simple average of the LUE model ensemble is widely used in many studies [9]. Thus, how to integrate the outputs of multiple models to obtain a more reliable GPP is valuable for the exploration of ecosystem function variations with climate change. The Bayesian model averaging (BMA) method [10] has been witnessed to be an effective method to integrate the output of multiple models, and it has been used to simulate water flux using the model ensemble [11]. However, how the BMA performs for the LUE model’s ensemble integration has not been well studied. Thus, the primary objectives of this study are (1) to evaluate the BMA method at dominant vegetation types in the study area and (2) to obtain a more accurate GPP estimation in the study area using the ensemble of LUE models.

2. Materials and Methods

2.1. Study Area

The Heihe River Basin is the second-largest inland river basin in China and is located in the northwest arid region of China. The elevation ranges from 5300 m upstream to 900 m downstream. It contains diverse landscapes, such as glaciers, permafrost, alpine steppe/meadow, and forests in the upstream and irrigated croplands, riparian vegetation, wetlands, and gobi desert in the middle stream and downstream [12]. An EC measurement network was established in 2013 in the Heihe River Basin [13,14]. The flux between the ecosystem and the atmosphere and ancillary meteorological variables were measured in the dominant vegetation types of this region. The measured data was collected from five sites from 2014 to 2015. The detailed information of these sites is listed in Table 1, while the spatial distribution of these sites is shown in Figure 1. The Arou site and the Dashalong site are located in the upstream, and the vegetation is the alpine meadow. The Daman site and the Shidi site are located in the middle stream, and the vegetation is maize (cropland) and reed (wetland), respectively. The Sidaoqiao site is located in the riparian zone of the downstream, and the vegetation at this site is tamarisk.
The EC measured 10 Hz data was processed using despike, coordinate rotation, frequency response correction, and WPL correction to obtain the half-hourly data [15]. Then the half-hourly EC data was gap-filled and partitioned into the GPP and ecosystem respiration using the REddyProc package (https://cran.r-project.org/web/packages/REddyProc (accessed on 10 January 2019)) [16] in the R environment. The measured air temperature, downward shortwave radiation, relative humidity, soil moisture, and GPP daily value were used to force or validate the LUE model. The MODIS reflectance at these sites was extracted from the MOD09A1 collection 6. The normalized difference vegetation index (NDVI), the enhanced vegetation index (EVI), and the land surface water index (LSWI) were calculated from the MODIS reflectance data at these sites.

2.2. LUE Models

The LUE model is an effective tool to estimate the GPP using remote sensing data. LUE models can be generalized as the following formula:
GPP   =   FAPAR   ×   PAR   ×   L U E m a x ×   f ( T ,   S ,   --- )
where FAPAR is the fraction of absorbed photosynthetically active radiation (PAR). L U E m a x is the maximum light using efficiency, and f(T, S, ---) is an environmental scalar on the photosynthesis rate. According to different schemes to quantify the environmental stress on photosynthesis, many LUE models have been developed. In this study, five widely used LUE models were selected to simulate the GPP seasonal dynamic at dominant vegetation types of the study area. They were Global Production Efficiency Model (GLO-PEM) [17], Vegetation Photosynthesis Model (VPM) [18], Carbon fixation model (C-Fix) [5], Eddy Covariance-Light Use Efficiency (EC-LUE) [6], and MODIS GPP product algorithm (MODIS-PSN) [7]. All the models were forced using in-situ measured meteorological data and default parameters for each site. Information regarding these models is listed in Table 2.

2.3. Bayesian Model Averaging

Bayesian model averaging is a scheme for the integration of outputs from multiple models. More detail regarding the theory can be found in the literature [10]. In the BMA method, the posterior distribution of the simulation variable, y, can be calculated using the following expression:
p ( y | D ) = i = 1 k p ( f i | D ) p i ( y | f i , D )
where y is the variable to be simulated, D = [ y 1 o b s , y 2 o b s , , y T o b s ] is the training data with length T, f = [ f 1 , f 2 , , f k ] is the ensemble of outputs from k different models, p ( f i | D ) is the posterior probability of a simulation given the observation D, and p i ( y | f i , D ) is the posterior distribution of y given the model simulation f i and the observed dataset D.
Finally, the posterior mean is estimated using the following formula:
E [ y | D ] =   i = 1 k p ( f i | D ) · E [ p i ( y | f i , D ) ] = i = 1 k w i ( a i + b i f i )
where w i is the weight for model i and i = 1 k w i = 1 ; a i and b i are bias correction terms that are derived using a simple linear regression of y on f i .
Measurements obtained from the five sites in 2014 were used to train the BMA method, and measurements obtained in 2015 were used to validate the performance of the BMA method.

2.4. Model Accuracy Evaluation

The determinant coefficient, R2, and the root mean square error (RMSE) were used to evaluate the model performance. The formulas for R2 and the RMSE are as follows:
R 2 = 1   i = 1 n ( P i O i ) 2 i = 1 n ( O i O ¯ ) 2
R M S E = 1 N i = 1 n ( P i O i ) 2
where P i is the simulated GPP, O i is the EC-measured GPP in this study,   O ¯ is the average of all the measurements, and N is the number of measurements.

3. Results

3.1. Carbon Flux Dynamics in the Heihe River Basin

The seasonal carbon flux dynamics of the five dominant vegetation types in the Heihe River Basin are shown in Figure 2. The annual carbon flux is summarized in Table 3. Carbon flux at these sites had an obvious seasonal cycle. The carbon flux began increasing in early May, reached the maximum in July and returned to the minimum in later October (Figure 2). Alpine grassland is the dominant vegetation type in the upper stream of the study area. Both Arou and Dashalong are covered with alpine grass. From the carbon flux measurement, the annual GPP was 862.76 gC/m2/year at Arou and 477.91 gC/m2/year at Dashalong. The annual Net Ecosystem Exchange (NEE) was −144.08 gC/m2/year at Arou and −309.89 gC/m2/year at Dashalong. The annual Ecosystem Respiration (ER) was 718.69 gC/m2/year at Arou and 168.02 gC/m2/year at Dshalong. Artificial and natural oases are the dominant landscapes in the middle stream of the Heihe River Basin. Carbon flux at two sites in the oases, a cropland (Daman) and a wetland (Shidi), was measured. In the middle stream, the annual GPP for the cropland was 1364 gC/m2/year and 1087 gC/m2/year for the wetland, the annual NEE for the cropland was −688 gC/m2/year and −585 gC/m2/year for the wetland, and the annual ER for the cropland was 676 gC/m2/year and 502 gC/m2/year for the wetland. Natural oasis dominates the vegetated area in the downstream of the study area. The Sidaoqiao site is located in the natural oasis region of the downstream. The vegetation type is tamarisk at the Sidaoqiao site. This site is close to the main river channel. At the Sidaoqiao site, the annual GPP was 709 gC/m2/year, the annual NEE was −201 gC/m2/year, and the annual ER was 508 gC/m2/year. All of these sites were carbon sinks and absorbed more than 100 gC/m2/year.

3.2. The BMA Based GPP Estimation

The determinant coefficient and RMSE were calculated between the simulated GPP and the EC-measured GPP (Table 4). As expected, the results from the trained BMA method had the highest accuracy. Averaging the R2 and RMSE among all of the sites indicated that BMA had the highest R2 and the lowest RMSE among all the models during the training and validation stages (Table 4). However, if one looks at one of the five sites, the BMA was not always the highest-accuracy method, and some LUE models performed better than the BMA. The LUE models performed differently at different sites. At a specific site, the performance of the LUE models varied greatly. For example, the R2 of the C-fix model in 2014 at the Dashalong site was 0.73, but the R2 was 0.97 in 2015 at the Daman site. At the Arou site, the MODIS-PSN model performed better than other LUE models, but the C-fix model performed better than other LUE models at the Sidaoqiao site. Overall, at each site, the BMA method was close to the best one of the LUE models. To further explore the performance of these methods, the residual is plotted for each method (Figure 3). The residuals of the BMA method were uniformly distributed around 0 and were less than the standard deviation (SD) of the measurements. The residuals of the BMA method were smaller than all the LUE models. All the LUE models had values outside the 1x standard deviation range. Even GLO-PEM and EC-LUE had values outside the 2x standard deviation range. The convergence trajectory of the coefficients for LUE models at the Arou site is shown in Figure 4. The initial coefficients for each LUE model were assigned a fixed value, and the value changed during the iteration and quickly reached a steady-state (Figure 4).

4. Discussion

The high spatial heterogeneity of GPP is an obstacle to accurately obtain regional GPP. We found that the carbon flux obtained using the eddy covariance method varied greatly among the different vegetation types. The artificial oasis dominated by croplands showed the highest GPP and carbon sequestration capacity. This indicates that both natural factors and human activity have a great impact on carbon flux and cannot be neglected. Eddy covariance is one of the commonly used techniques to measure ecosystem-scale carbon flux, but still it involves many uncertainties that should be noted when analyzing the data. The uncertainties primarily result from three aspects. First, the eddy covariance data quality is affected by sensor configurations and meteorological conditions [19]. Hence, strict quality control is required for the raw data before its use. Second, due to instrument malfunction and spikes, data gaps are inevitable in the raw data, and gap-filling is a necessary step to obtain continuous carbon flux [20]. This process also leads to uncertainties in the results. Third, the eddy covariance method only directly measures the NEE, and the GPP and ER are partitioned from the NEE with some hypothesis, a fact that can also result in uncertainties in the GPP and ER [21,22].
The LUE model performance changed greatly with vegetation types. The LUE model is a widely used method to estimate the GPP, especially at the regional scale. However, GPP simulated using the LUE model contains uncertainties that result from the model parameters, meteorological input data, and remote sensing input data [23,24]. Model parameters vary with species and vegetation community compositions, and the LUE models usually only provide default parameters for the major vegetation types. When applying the LUE models in a small region, the default parameters may result in great uncertainties. The uncertainties that result from meteorological input data are primarily caused by gap filling and the interpolation of measured data. Remote sensing input data (such as FPAR, NDVI, EVI, LSWI and others) in the LUE model also contain uncertainties that result from bad data quality [25]. Some studies reported that model parameters have a larger impact on the simulation accuracy than the meteorological input data and remote sensing input data [26].
BMA is an effective method to improve modeling accuracy. Every model has some limitations caused by its hypothesis and parameterization scheme [27]. The average of multiple model outputs can overcome this problem and be thought of as more reliable than a single model output. In many studies, the classical average of a model ensemble has been used to conduct the research [28,29]. BMA provides a way to obtain the model ensemble weighted average [30]. By training the BMA method with measurement data, different weights are assigned to different models, and the weights also vary with spatial and vegetation type changes. This is more reasonable than the classical average because the performances of the LUE models are different at different locations and with different vegetation types [31,32]. Many studies have reported that BMA can improve the simulation accuracy than the classical average in stream-flow forecasting models [33] and evapotranspiration models [11]. This is consistent with the finding of the current study, where the RMSE of the BMA method was 5.8 gC/m2/8 day during the training period and 8.1 gC/m2/8 day during the validation period, while the classical averages were 12.6 gC/m2/8 day and 13.98 gC/m2/8 day, respectively. This demonstrates that BMA has the potential to improve GPP estimation by integrating multiple LUE models. Meanwhile, training process is vital to BMA application. The short data period for BMA training may result in some uncertainties in the simulation in this study. The training data should be as comprehensive as possible. Additionally, BMA is also able to calculate the uncertainties in its estimated GPP, which is useful when conducting the analysis with remote sensed GPP. In this study, we take the Heihe River Basin as an example to demonstrate the performance of the BMA method. As a universal method, if BMA is well-calibrated at the global scale, it can integrate the global GPP products and output high quality global GPP estimation.

5. Conclusions

In this study, five LUE models and the multi-model ensemble method were evaluated using measured daily carbon flux and meteorological data in the Heihe River Basin. The measured carbon flux data indicate that the GPP varied greatly among the different ecosystems in the Heihe River Basin from upstream to downstream. The cropland and wetland in the middle stream had the highest carbon uptake capacities. Performance of the LUE models indicated that the LUE models contain uncertainties resulting from model parameters and model struct. Single LUE model performance is not robust in different vegetation types. The LUE models without parameter calibration exhibited great inconsistencies in GPP simulations. The BMA method is an effective tool to integrate multi-model ensemble output. The RMSE between the BMA-simulated GPP and EC-measured GPP was the smallest. The use of the BMA method allows for the combination of multiple LUE models to yield a high-accuracy GPP.

Author Contributions

Conceptualization, X.W. and J.Z.; formal analysis, X.W. and J.Z.; writing—original draft preparation, X.W.; writing—review and editing, X.W., J.Z. and J.R.; All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No. 41771466), and the Special Fund for the Key Program of Science and Technology of Qinghai Province (Grant no. 2017-SF-A6).

Data Availability Statement

All data and codes used in this study are available by requesting the corresponding author ([email protected]).

Acknowledgments

This work was supported by the Special Fund for Key Program of Science and Technology of Qinghai Province (Grant no. 2017-SF-A6). We would like to thank the National Tibetan Plateau Data Center (http://data.tpdc.ac.cn (accessed on 8 October 2018)) for providing the carbon flux and meteorological data. We would also like to thank the Land Processes Distributed Active Center (LP DAAC) for providing MODIS reflectance product MOD09A1 collection6. We thank LetPub (www.letpub.com (accessed on 8 March 2021)) for its linguistic assistance during the preparation of this manuscript. We thank the three anonymous reviewers for their constructive suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Beer, C.; Reichstein, M.; Tomelleri, E.; Ciais, P.; Jung, M.; Carvalhais, N.; Rödenbeck, C.; Arain, M.A.; Baldocchi, D.; Bonan, G.B.; et al. Terrestrial Gross Carbon Dioxide Uptake: Global Distribution and Covariation with Climate. Science 2010, 329, 834–838. [Google Scholar] [CrossRef] [Green Version]
  2. Turner, D.P.; Ritts, W.D.; Cohen, W.B.; Gower, S.T.; Zhao, M.; Running, S.W.; Wofsy, S.C.; Urbanski, S.; Dunn, A.L.; Munger, J.W. Scaling Gross Primary Production (GPP) over boreal and deciduous forest landscapes in support of MODIS GPP product validation. Remote Sens. Environ. 2003, 88, 256–270. [Google Scholar] [CrossRef] [Green Version]
  3. Yuan, W.; Cai, W.; Xia, J.; Chen, J.; Liu, S.; Dong, W.; Merbold, L.; Law, B.; Arain, A.; Beringer, J.; et al. Global comparison of light use efficiency models for simulating terrestrial vegetation gross primary production based on the LaThuile database. Agric. For. Meteorol. 2014, 192–193, 108–120. [Google Scholar] [CrossRef]
  4. Baldocchi, D.; Falge, E.; Gu, L.; Olson, R.; Hollinger, D.; Running, S.; Anthoni, P.; Bernhofer, C.; Davis, K.; Evans, R.; et al. FLUXNET: A New Tool to Study the Temporal and Spatial Variability of Ecosystem-Scale Carbon Dioxide, Water Vapor, and Energy Flux Densities. Bull. Am. Meteorol. Soc. 2001, 82, 2415–2434. [Google Scholar] [CrossRef]
  5. Veroustraete, F.; Sabbe, H.; Eerens, H. Estimation of carbon mass fluxes over Europe using the C-Fix model and Euroflux data. Remote Sens. Environ. 2002, 83, 376–399. [Google Scholar] [CrossRef]
  6. Yuan, W.; Liu, S.; Zhou, G.; Zhou, G.; Tieszen, L.L.; Baldocchi, D.; Bernhofer, C.; Gholz, H.; Goldstein, A.H.; Goulden, M.L.; et al. Deriving a light use efficiency model from eddy covariance flux data for predicting daily gross primary production across biomes. Agric. For. Meteorol. 2007, 143, 189–207. [Google Scholar] [CrossRef] [Green Version]
  7. Running, S.W.; Nemani, R.R.; Heinsch, F.A.; Zhao, M.; Reeves, M.; Hashimoto, H. A Continuous Satellite-Derived Measure of Global Terrestrial Primary Production. Bioscience 2004, 54, 547–560. [Google Scholar] [CrossRef]
  8. Xiao, X.; Hollinger, D.; Aber, J.; Goltz, M.; Davidson, E.A.; Zhang, Q.; Moore, B. Satellite-based modeling of gross primary production in an evergreen needleleaf forest. Remote Sens. Environ. 2004, 89, 519–534. [Google Scholar] [CrossRef]
  9. Ballantyne, A.; Smith, W.; Anderegg, W.; Kauppi, P.; Sarmiento, J.; Tans, P.; Shevliakova, E.; Pan, Y.; Poulter, B.; Anav, A.; et al. Accelerating net terrestrial carbon uptake during the warming hiatus due to reduced respiration. Nat. Clim. Chang. 2017, 7, 148–152. [Google Scholar] [CrossRef]
  10. Duan, Q. Multi-model ensemble hydrologic prediction using Bayesian model averaging. Adv. Water Resour. 2007, 30, 1371–1386. [Google Scholar] [CrossRef] [Green Version]
  11. Zhu, G.; Li, X.; Zhang, K.; Ding, Z.; Han, T.; Ma, J.; Huang, C.; He, J.; Ma, T. Multi-model ensemble prediction of terrestrial evapotranspiration across north China using Bayesian model averaging. Hydrol. Process. 2016, 30, 2861–2879. [Google Scholar] [CrossRef]
  12. Li, X.; Li, X.; Li, Z.; Ma, M.; Wang, J.; Xiao, Q.; Liu, Q.; Che, T.; Chen, E.; Yan, G.; et al. Watershed Allied Telemetry Experimental Research. J. Geophys. Res. 2009, 114. [Google Scholar] [CrossRef] [Green Version]
  13. Li, X.; Cheng, G.D.; Liu, S.M.; Xiao, Q.; Ma, M.G.; Jin, R.; Che, T.; Liu, Q.H.; Wang, W.Z.; Qi, Y.; et al. Heihe Watershed Allied Telemetry Experimental Research (HiWATER): Scientific Objectives and Experimental Design. Bull. Am. Meteorol. Soc. 2013, 94, 1145–1160. [Google Scholar] [CrossRef]
  14. Liu, S.; Li, X.; Xu, Z.; Che, T.; Xiao, Q.; Ma, M.; Liu, Q.; Jin, R.; Guo, J.; Wang, L.; et al. The Heihe Integrated Observatory Network: A Basin-Scale Land Surface Processes Observatory in China. Vadose Zone J. 2018, 17. [Google Scholar] [CrossRef]
  15. Zhang, Z.; Wang, W.; Ma, M.; Xu, Z.; Wu, Y.; Huang, G.; Tan, J. Data Processing and Product Analysis of Eddy Covariance Flux Data for WATER. Remote Sens. Technol. Appl. 2010, 25, 788–796. [Google Scholar]
  16. Wutzler, T.; Lucas-Moffat, A.; Migliavacca, M.; Knauer, J.; Sickel, K.; Šigut, L.; Menzer, O.; Reichstein, M. Basic and extensible post-processing of eddy covariance flux data with REddyProc. Biogeosciences 2018, 15, 5015–5030. [Google Scholar] [CrossRef] [Green Version]
  17. Prince, S.D.; Goward, S.N. Global Primary Production: A Remote Sensing Approach. J. Biogeogr. 1995, 22, 815–835. [Google Scholar] [CrossRef]
  18. Xiao, X.; Zhang, Q.; Braswell, B.; Urbanski, S.; Boles, S.; Wofsy, S.; Moore, B.; Ojima, D. Modeling gross primary production of temperate deciduous broadleaf forest using satellite images and climate data. Remote Sens. Environ. 2004, 91, 256–270. [Google Scholar] [CrossRef]
  19. Foken, T.; Wichura, B. Tools for quality assessment of surface-based flux measurements. Agric. For. Meteorol. 1996, 78, 83–105. [Google Scholar] [CrossRef]
  20. Papale, D.; Reichstein, M.; Aubinet, M.; Canfora, E.; Bernhofer, C.; Kutsch, W.; Longdoz, B.; Rambal, S.; Valentini, R.; Vesala, T.; et al. Towards a standardized processing of Net Ecosystem Exchange measured with eddy covariance technique: Algorithms and uncertainty estimation. Biogeosciences 2006, 3, 571–583. [Google Scholar] [CrossRef] [Green Version]
  21. LASSLOP, G.; REICHSTEIN, M.; PAPALE, D.; RICHARDSON, A.D.; ARNETH, A.; BARR, A.; STOY, P.; WOHLFAHRT, G. Separation of net ecosystem exchange into assimilation and respiration using a light response curve approach: Critical issues and global evaluation. Glob. Chang. Biol. 2010, 16, 187–208. [Google Scholar] [CrossRef] [Green Version]
  22. Reichstein, M.; Falge, E.; Baldocchi, D.; Papale, D.; Aubinet, M.; Berbigier, P.; Bernhofer, C.; Buchmann, N.; Gilmanov, T.; Granier, A.; et al. On the separation of net ecosystem exchange into assimilation and ecosystem respiration: Review and improved algorithm. Glob. Chang. Biol 2005, 11, 1424–1439. [Google Scholar] [CrossRef]
  23. Wang, X.; Ma, M.; Huang, G.; Veroustraete, F.; Zhang, Z.; Song, Y.; Tan, J. Vegetation primary production estimation at maize and alpine meadow over the Heihe River Basin, China. Int. J. Appl. Earth Obs. Geoinf. 2012, 17, 94–101. [Google Scholar] [CrossRef]
  24. Wang, X.; Ma, M.; Li, X.; Song, Y.; Tan, J.; Huang, G.; Zhang, Z.; Zhao, T.; Feng, J.; Ma, Z.; et al. Validation of MODIS-GPP product at 10 flux sites in northern China. Int. J. Remote Sens. 2013, 34, 587–599. [Google Scholar] [CrossRef]
  25. Wang, X.; Xiao, J.; Li, X.; Cheng, G.; Ma, M.; Che, T.; Dai, L.; Wang, S.; Wu, J. No Consistent Evidence for Advancing or Delaying Trends in Spring Phenology on the Tibetan Plateau. J. Geophys. Res. 2017, 122, 3288–3305. [Google Scholar] [CrossRef]
  26. Wang, H.; Li, X.; Ma, M.; Geng, L. Improving Estimation of Gross Primary Production in Dryland Ecosystems by a Model-Data Fusion Approach. Remote Sens. 2019, 11, 225. [Google Scholar] [CrossRef] [Green Version]
  27. Parrish, M.A.; Moradkhani, H.; DeChant, C.M. Toward reduction of model uncertainty: Integration of Bayesian model averaging and data assimilation. Water Resour. Res. 2012, 48. [Google Scholar] [CrossRef]
  28. Zhang, Y.; Xiao, X.; Guanter, L.; Zhou, S.; Ciais, P.; Joiner, J.; Sitch, S.; Wu, X.; Nabel, J.; Dong, J.; et al. Precipitation and carbon-water coupling jointly control the interannual variability of global land gross primary production. Sci. Rep. 2016, 6, 39748. [Google Scholar] [CrossRef] [Green Version]
  29. Dormann, C.F.; Calabrese, J.M.; Guillera-Arroita, G.; Matechou, E.; Bahn, V.; Bartoń, K.; Beale, C.M.; Ciuti, S.; Elith, J.; Gerstner, K.; et al. Model averaging in ecology: A review of Bayesian, information-theoretic, and tactical approaches for predictive inference. Ecol. Monogr. 2018, 88, 485–504. [Google Scholar] [CrossRef] [Green Version]
  30. Raftery, A.E.; Gneiting, T.; Balabdaoui, F.; Polakowski, M. Using Bayesian Model Averaging to Calibrate Forecast Ensembles. Mon. Weather Rev. 2005, 133, 1155–1174. [Google Scholar] [CrossRef] [Green Version]
  31. Wu, C.; Munger, J.W.; Niu, Z.; Kuang, D. Comparison of multiple models for estimating gross primary production using MODIS and eddy covariance data in Harvard Forest. Remote Sens. Environ. 2010, 114, 2925–2939. [Google Scholar] [CrossRef]
  32. Sun, Z.; Wang, X.; Zhang, X.; Tani, H.; Guo, E.; Yin, S.; Zhang, T. Evaluating and comparing remote sensing terrestrial GPP models for their response to climate variability and CO2 trends. Sci. Total Environ. 2019, 668, 696–713. [Google Scholar] [CrossRef] [PubMed]
  33. Vrugt, J.A.; Diks, C.G.H.; Clark, M.P. Ensemble Bayesian model averaging using Markov Chain Monte Carlo sampling. Environ. Fluid Mech. 2008, 8, 579–595. [Google Scholar] [CrossRef]
Figure 1. Distribution of the study sites in the Heihe River Basin.
Figure 1. Distribution of the study sites in the Heihe River Basin.
Land 10 00329 g001
Figure 2. Seasonal carbon flux dynamics at the five study sites.
Figure 2. Seasonal carbon flux dynamics at the five study sites.
Land 10 00329 g002
Figure 3. Residual plots of different models at the five sites (a: VPM; b: EC-LUE; c: C-fix; d: MODIS-PSN; e: GLO-PEM; f: BMA). The dashed line is the 1*SD and 2*SD (SD = 25 gC/m2/8 day for the GPP measurement at Arou site).
Figure 3. Residual plots of different models at the five sites (a: VPM; b: EC-LUE; c: C-fix; d: MODIS-PSN; e: GLO-PEM; f: BMA). The dashed line is the 1*SD and 2*SD (SD = 25 gC/m2/8 day for the GPP measurement at Arou site).
Land 10 00329 g003
Figure 4. The weights convergence trajectory of the LUE models during the BMA method training.
Figure 4. The weights convergence trajectory of the LUE models during the BMA method training.
Land 10 00329 g004
Table 1. Information of the Eddy Covariance (EC) sites in the Heihe River Basin.
Table 1. Information of the Eddy Covariance (EC) sites in the Heihe River Basin.
Site NameLand CoverLatitudeLongitudeElevation (m)
ArouGrassland38.0473100.46433033
DashalongGrassland38.839998.94063739
DamanCropland (Maize)38.85551100.37221556
ShidiWetland38.97514100.44641460
SidaoqiaoShrubland42.0012101.1374873
Table 2. The Light Use Efficiency (LUE) models used in this study.
Table 2. The Light Use Efficiency (LUE) models used in this study.
ModelMain FormulaInput Variables
GLO-PEMGPP = APAR × LUEmax× f(T) × f(SM) × f(VPD)PAR, NDVI, T, SM, VPD
VPMGPP = APAR × LUEmax× f(T) × f(LSWI) × f(P)PAR, EVI, T, LSWI
C-FixGPP = APAR × LUEmax× f(T) × f(CO2)PAR, NDVI, T, CO2
EC-LUEGPP = APAR × LUEmax× f(T) × f(EF)PAR, NDVI, T, EF
MODIS-PSNGPP = APAR × LUEmax× f(TMIN) × f(VPD)PAR, NDVI, TMIN, VPD
T: air temperature; SM: soil moisture; VPD: vapor pressure deficit; EF: evaporative fraction.
Table 3. Annual carbon flux at the five study sites in the Heihe River Basin.
Table 3. Annual carbon flux at the five study sites in the Heihe River Basin.
Site NameGPP (gC/m2/yr)NEE (gC/m2/yr)ER (gC/m2/yr)
201420152014201520142015
Arou845.90879.63−175.61−112.54670.29767.08
Dashalong506.58449.24−314.61−305.16191.97144.08
Daman1329.961397.60−726.94−648.23603.02749.36
Shidi1008.571165.77−559.52−609.76449.05556.02
Sidaoqiao664.82753.58−190.00−211.96474.82541.62
Table 4. Determinant coefficients and RMSEs between the simulated GPP and the EC-measured GPP. (Data measured during 2014/2015 was used to train/validate the BMA method. The highest R2 and the smallest RMSE were bolded.).
Table 4. Determinant coefficients and RMSEs between the simulated GPP and the EC-measured GPP. (Data measured during 2014/2015 was used to train/validate the BMA method. The highest R2 and the smallest RMSE were bolded.).
Site R2RMSE (gC/m2/8 day)
VPMEC-LUEC-FixGLO-PEMMODIS_PSNBMAVPMEC-LUEC-FixGLO-PEMMODIS_PSNBMA
ArouTraining0.950.950.950.950.970.966.78.18.426.06.04.8
Validation0.950.940.940.940.950.9510.17.010.923.810.36.9
DashalongTraining0.830.790.730.790.790.88.87.110.710.99.66.0
Validation0.880.910.840.820.760.877.13.98.213.28.44.7
DamanTraining0.950.940.950.940.960.9511.513.69.623.910.98.4
Validation0.960.960.970.950.980.9715.78.48.241.56.78.6
ShidiTraining0.830.850.870.810.820.849.915.514.739.513.68.6
Validation0.750.820.790.830.820.8118.416.117.040.913.515.1
SidaoqiaoTraining0.930.910.940.930.910.955.77.24.824.38.53.4
Validation0.850.840.870.850.780.879.39.77.425.711.67.6
All sitesTraining0.900.900.900.890.900.919.511.010.221.610.85.8
Validation0.880.900.890.880.850.9012.79.910.925.111.38.1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, J.; Wang, X.; Ren, J. Simulation of Gross Primary Productivity Using Multiple Light Use Efficiency Models. Land 2021, 10, 329. https://doi.org/10.3390/land10030329

AMA Style

Zhang J, Wang X, Ren J. Simulation of Gross Primary Productivity Using Multiple Light Use Efficiency Models. Land. 2021; 10(3):329. https://doi.org/10.3390/land10030329

Chicago/Turabian Style

Zhang, Jun, Xufeng Wang, and Jun Ren. 2021. "Simulation of Gross Primary Productivity Using Multiple Light Use Efficiency Models" Land 10, no. 3: 329. https://doi.org/10.3390/land10030329

APA Style

Zhang, J., Wang, X., & Ren, J. (2021). Simulation of Gross Primary Productivity Using Multiple Light Use Efficiency Models. Land, 10(3), 329. https://doi.org/10.3390/land10030329

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