Next Article in Journal
Irrigation and Precipitation Hydrological Consistency with SMOS, SMAP, ESA-CCI, Copernicus SSM1km, and AMSR-2 Remotely Sensed Soil Moisture Products
Next Article in Special Issue
Mapping Forest Aboveground Biomass Using Multisource Remotely Sensed Data
Previous Article in Journal
Assessing Scale Dependence on Local Sea Level Retrievals from Laser Altimetry Data over Sea Ice
Previous Article in Special Issue
A New Remote Sensing Method to Estimate River to Ocean DOC Flux in Peatland Dominated Sarawak Coastal Regions, Borneo
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Combining Phenological Camera Photos and MODIS Reflectance Data to Predict GPP Daily Dynamics for Alpine Meadows on the Tibetan Plateau

1
College of Geography and Environment Sciences, Northwest Normal University, Lanzhou 730070, China
2
Key Laboratory of Remote Sensing of Gansu Province, Heihe Remote Sensing Experimental Research Station, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou 730000, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(22), 3735; https://doi.org/10.3390/rs12223735
Submission received: 19 October 2020 / Revised: 11 November 2020 / Accepted: 11 November 2020 / Published: 13 November 2020
(This article belongs to the Special Issue Remote Sensing of Carbon Cycle Science)

Abstract

:
Gross primary production (GPP) is the overall photosynthetic fixation of carbon per unit space and time. Due to uncertainties resulting from clouds, snow, aerosol, and topography, it is a challenging task to accurately estimate daily GPP. Daily digital photos from a phenological camera record vegetation daily greenness dynamics with little cloud or aerosol disturbance. It can be fused with satellite remote sensing data to improve daily GPP prediction accuracy. In this study, we combine the two types of datasets to improve the estimation accuracy of GPP for alpine meadow on the Tibetan Plateau. To examine the performance of different methods and vegetation indices (VIs), three experiments were designed. First, GPP was estimated with the light use efficiency (LUE) model with the green chromatic coordinate (GCC) from the phenological camera and vegetation index from MODIS, respectively. Second, GPP was estimated with the Backpropagation neural network machine learning algorithm (BNNA) method with GCC from the phenological camera and vegetation index from MODIS, respectively. Finally, GPP was estimated with the BNNA method using GCC and vegetation index as inputs at the same time. Compared with eddy covariance GPP, GPP predicted by the BNNA method with GCC and vegetation indices as inputs at the same time showed the highest accuracy of all the experiments. The results indicated that GCC had a higher accuracy than NDVI and EVI when only one vegetation index data was used in the LUE model or the BNNA method. The R2 of GPP estimated by BNNA and GPP from eddy covariance increased by 0.12 on average, RMSE decreased by 1.13 g C·m−2·day−1 on average, and MAD decreased by 0.87 g C·m−2·day−1 on average compared with GPP estimated by the traditional LUE model and GPP from eddy covariance. This study puts forth a new way to improve the estimation accuracy of GPP on the Tibetan Plateau. With the emergence of a large number of phenological cameras, this method has great potential for use on the Tibetan Plateau, which is heavily affected by clouds and snow.

Graphical Abstract

1. Introduction

Gross primary production (GPP) refers to the overall photosynthetic fixation of carbon per unit space and time [1]. GPP is extremely important for studying the terrestrial ecosystem carbon cycle and global climate change, because GPP reflects not only the vegetation activity, but also the response of the ecosystem to climate change [2,3,4,5]. Therefore, accurately estimating the temporal and spatial dynamics of GPP on a regional and global scale is of great significance for understanding the terrestrial ecosystem carbon cycle, assessing grain production, predicting future climate change scenarios and making climate-policy decisions [6,7,8]. Although there are many GPP models have been developed in the past few decades, it is still a challenging work to accurately estimate daily GPP.
Monteith first proposed the method of estimating vegetation productivity by using photosynthetic effective radiation absorbed by vegetation and established the basis of light use efficiency model [1]. This method is widely used in regional GPP estimation with remote sensing data [9,10]. However, owing to the influence of atmospheric conditions, cloud and sensor degradation, there were large uncertainties in GPP estimated by remote sensing vegetation index data and the light use efficiency (LUE) model [11], such as the CASA model [12], the global production efficiency model (GLO-PEM) [13], the vegetation photosynthesis model (VPM) [14] or the eddy covariance-light use efficiency (EC-LUE) model [15]. However, it is reported that the satellite remote sensing-based GPP contains many problems caused by snow and cloud contaminants and low temporal resolution. For example, GPP was underestimated at the majority of sites in African ecosystems [16] and has unstable performance at various sites [17,18].
Phenological cameras can provide almost time-continuous images at a small spatial scale; meanwhile, they are rarely disturbed by clouds [19,20,21]. Phenological cameras have the advantages of low cost, small size, convenient installation, and easy handling, and they are widely used to monitor plant dynamics [22]. The phenological camera-based color index was also used to estimate canopy GPP in many studies [23,24,25]. The performance of GPP estimated by color indexes from phenological camera photos was of varying vegetation types, and the estimation effect of GPP on grassland was better than that of evergreen coniferous forests and deciduous broad-leaved forests [11]. Thus, the phenology camera can provide more reliable daily greenness dynamics and can be fused with satellite remote sensing data to improve daily GPP estimate for alpine meadow ecosystems.
In recent years, data-based machine learning has been widely used in remote sensing data fusion and model integration. Many studies have indicated that machine learning works well in vegetation productivity estimation [26,27,28]. Machine learning models find rules by training sample data and use these rules to predict unknown results. GPP estimation based on machine learning algorithms can make full use of existing multisource remote sensing data, and they can provide a long time series and high-precision estimation of GPP. Therefore, machine learning has advantage of combining multisource data by machine learning to estimate GPP. The machine learning method predicted GPP with an RMSE of 1.87 g C·m−2·day−1 and an R2 of 0.71 using meteorological and flux data from 33 AmeriFlux sites between 2000 and 2003 in America [29]. The machine learning method was also used to estimate eight days GPP (at the site level) with superior performance on the European spatial scale [30]. Lee successfully estimated the daily GPP of the forest ecosystem using the machine learning method with MODIS data at site scale, and the RMSE reached 1.03 g C·m−2·day−1 [31]. These studies mentioned above indicated that the machine learning methods is an effective and widely used method to correctly predict GPP.
The rapid use of phenology cameras in the Tibetan Plateau provides high-quality data that reflect alpine seasonal changes. However, in what extent it can be used to improve GPP estimation for alpine meadows has not been studied. In view of the advantages and disadvantages of phenology camera digital photos and MODIS data, combing digital camera and satellite remote sensing data is an effective way to improve GPP prediction accuracy on the Tibetan Plateau because digital cameras and satellite remote sensing provide complementary information. Thus, the objective of this study are to: (1) compare the performance of color index from phenological photos and vegetation indices from satellite remote sensing data in GPP estimation for alpine meadows on the Tibetan Plateau; (2) compare the performance of the LUE model and a machine learning method in GPP estimation, and (3) explore to what extent the GPP estimation accuracy can be improved by combining multisource remote sensing data for alpine meadows on the Tibetan Plateau.

2. Materials and Methods

2.1. Study Area

The Tibetan Plateau is the largest and highest plateau on Earth [32]. The ecosystems of the Tibetan Plateau are sensitive and fragile, and they are sensitive to global climate change. Alpine meadows are a dominant vegetation type on the Tibetan Plateau. The A’rou (AR) eddy covariance (EC) observation station (Figure 1) is located in the eastern part of the Tibetan Plateau. The land surface of the AR station is alpine meadow. The longitude and latitude of the AR station are 100°27′52.9″ E and 38°02′39.8″ N and the altitude is 3033 m. The mean annual temperature and annual precipitation are 0.7 °C and 400 mm, respectively. We divided the seasons on the Tibetan Plateau, with spring as May, summer as June to July, and autumn as August to September [33,34].

2.2. Meteorological Data

The meteorological data from 2017 to 2018 observed at the AR station were used as input to simulate GPP, including global solar radiation (RG, w/m2), saturated vapor pressure deficit (VPD, Pa), air temperature (T_2m, air, °C) at a height of 2 m above ground level, and photosynthetically active radiation (PAR, MJ/d/m2). These variables are observed every 10 min, and daily RG, VPD, PAR, and minimum temperature (Tmin) were calculated from raw 10 min data in this study (Figure 2).

2.3. Remote Sensing Data

MOD09GA is the daily MODIS surface reflectance product, which was atmospherically corrected [35]. The band reflectance data of MOD09GA can be obtained from Google Earth Engine (https://earthengine.google.com/) and its spatial resolution and temporal resolution are 250 m and daily, respectively [36]. The 3 × 3 pixels around the AR station were extracted from MOD09GA data (Figure 3). Daily NDVI and EVI were calculated from the reflectance data using the following formulas:
NDVI = ρ NIR ρ RED ρ NIR + ρ RED
EVI = 2.5 × ρ NIR ρ RED ρ NIR + 6 × ρ RED 7.5 × ρ BLUE + C
where ρRED, ρNIR, and ρBLUE are the band reflectance of the red, near-infrared, and blue band, respectively. C is the adjustment parameter, and the value of C is set to 1.

2.4. GPP from Eddy Covariance

GPP data from AR eddy covariance flux site (Figure 1) on the Tibetan Plateau were used to evaluate the GPPs derived from remote sensing data and phenological camera photos. Eddy covariance raw flux data were processed rigidly, including despike, time lag correction, coordinate rotation, frequency response correction, and WPL correction [37]. After gap-filling and partition, daily GPP was calculated based on ecosystem respiration (Re) and net ecosystem exchange (NEE) data by using Equations (3) and (4) [38]:
NEE = ε max × PAR × GEE max ε max × PAR + GEE max Re
GPP = NEE + Re
where εmax is the maximum LUE in g C/MJ APAR, GEEmax is the maximum gross ecosystem exchange in g C/m2/s [39].

2.5. Phenological Camera Photos and Color Index Extraction

Phenological cameras (NV201E digital cameras, Haoxin Netscape Network Technology Co., Ltd., Shenzhen, China) were installed at the AR station on 15th April, 2017. They continually record the surface status by taking photos almost vertically and toward the south to avoid the tower shadow. These phenological cameras took photos with a resolution of 1600 × 1200 and can specify the shooting time and frequency (Figure 4). During the observation period, photos were saved in JPEG format. These phenological cameras were operated in automatic exposure mode with the white balance adjustment in “automatic mode”. In 2017 year, the times the photos were taken were 7:30, 7:45, 8:00, 11:30, 11:45, 12:00, 17:30, 17:45, and 18:00. In 2018 year, the time was adjusted and the times were 8:30, 8:45, 9:00, 11:30, 11:45, 12:00, 17:30, 17:45, and 18:00 respectively. We deleted the photos when the light was too strong or too weak from the photo set, so as to get the photos of real recorded features under normal weather conditions. The green chromatic coordinate (GCC) was calculated to reflect the canopy greenness, which was widely used to monitor canopy development and was strongly related to GPP [40,41]. GCC was calculated using following formulas:
GCC 1 = DN G DN R + DN G + DN B
GCC 2 = DN G DN R + DN B
where DN is the digital number and R, G, and B were the pixel values of red, green, and blue channels, respectively. GCC1 and GCC2 are both green chromatic coordinate (GCC), but the denominators in the calculation Equations (5) and (6) were different.
To ensure the reliability of GCC data, we set some conditions on the quality of the photos as following: (1) Make sure that there were no underexposed photos. (2) No shadows of the the flux tower appeared in the photos. (3) No obvious animals and people appeared in the photos. We extracted digital RGB data from the ROI (400 × 600) of the images and calculated GCC1 and GCC2 for all photos. We found that the R2 of ROI_GCC1 and GPP from eddy covariance were higher than the R2 of ROI_GCC2 and GPP from eddy covariance (Figure 5), so ROI_GCC1 was used for all GPP estimations in this study.

2.6. GPP Model Based on VIs and GCC

The LUE model was used in study to estimate GPP, as follows:
GPP model = ε × FPAR × PAR
ε = ε max × f ( Tmin ) × f ( VPD )
FPAR = α × NDVI
FPAR = α × EVI
FPAR = α ×   CI norm
where ε is light use efficiency in photosynthesis (g C/MJ) and PAR (MJ/d/m2) is photosynthetically active radiation. PAR = n × RG. n is the proportion of PAR in the global total solar radiation and was not a fixed value. For calculations, n = 0.45 was used. FPAR is the fraction of absorbed photosynthetically active radiation by vegetation. In this study, FPAR was estimated from the normalized color index (CInorm) of GCC, NDVI and EVI by Equations (9)–(11). α is the coefficient constant and is assumed to be 1.0 [14]. The value of εmax was set to 1.6 g C·m−2·day−1·MJ−1 according to literature for alpine meadow [39].
f(Tmin) and f(VPD) are simple linear piecewise slope functions of daily minimum temperature (Tmin) and daily average VPD [42]. f(Tmin) and f(VPD) reflect the environmental stress on vegetation photosynthesis and are expressed as:
f ( Tmin ) = { 0 , Tmin < Tmin min Tmin Tmin min Tmin max Tmin min , Tmin min < Tmin < Tmin max 1 , Tmin > Tmin max
f ( VPD ) = { 0 , VPD > VPD max VPD max VPD VPD max VPD min , VPD min < VPD < VPD max 1 , VPD < VPD min
where VPDmax, VPDmin, Tminmax, and Tminmin are model parameters and the values were from the default parameters for grass in the BPLUT table [43].
The GCC values of RGB color photos of the AR station from 2017 to 2018 were calculated in Python 3.5, and the CInorm index was calculated with the GCC values (Equation (14)). The CImax and CImin were the maximum and minimum value of GCC in each year, respectively. The GPPs based on CInorm and the remote sensing VIs was conducted by Equations (15)–(17), and the equations are given by:
CI norm = CI CI min CI max CI min
GPP model = ε × CI norm × PAR
GPP model = ε × NDVI × PAR
GPP model = ε × EVI × PAR

2.7. BP Neural Network Algorithm

The BP neural network is a multilayer feedforward network trained by error back propagation. The BP neural network algorithm (BNNA) contains an input layer, hidden layer, and output layer, and each layer includes a specific number of neurons. We chose BNNA for this study because it had the following four capabilities: 1. BNNA has a strong nonlinear mapping ability; 2. BNNA has a highly self-learning and self-adaptive ability; 3. BNNA has the ability to apply learning results to new knowledge; 4. BNNA has a certain fault tolerance. We used Keras based on TensorFlow high-level Application Programming Interface (Google Artificial Intelligence Team Google Brain, Mountain view, Santa Clara County, CA, USA) and BNNA to achieve machine learning fusion of meteorological data, MOD09GA VIs, and phenological cameras GCC to predict GPP. First, the input datasets were normalized and mapped to the range [0,1], and then BNNA was trained with data in 2017. Training was stopped when BNNA reached the iteration stop condition (Terror = 1 × 10−7) and number of iterations (epochs = 100). The data in 2018 were used to validate the BNNA to test the performance of established model. R2, RMSE, Mean Absolute Deviation (MAD) indices between the GPP from eddy covariance and simulated GPP were calculated to evaluate the accuracy during training and validation periods. The structure and key parameters of BNNA are reported in Figure 6 and Table 1. The process of machine learning GPP modeling consisted of four groups of experiments (Table 2).
Too many or too few iterations of the training model in the BNNA will cause uncertainty in training errors. Only repeated experiments can sum up a relatively minimum number of trainings and obtain better training results. Finally, we set n = 100. If the activation function was not set properly, the BP neural network will not converge. We compared the convergence effects of different activation functions and found that the effects of the “tanh” activation was better than others for the training. The learning rate was generally set to 0.1 for BNNA.

2.8. Model Output Accuracy Evaluation

We calculated the R2, RMSE, and MAD in this study to evaluate the performances of the GPPs (GPP_NDVI, GPP_EVI, GPP_CInorm) estimated by the remote sensing GPP model and the GPPs (GPP_train1, GPP_train2, GPP_train3, GPP_train4, GPP_valid1, GPP_valid2, GPP_valid3, GPP_valid4) estimated by the BNNA (Figure 7).
P and O, represent the GPP simulation values and the GPP from eddy covariance, respectively, and i and n represent the days of the year and the total number of days in the year, respectively. The R2, RMSE, and MAD are expressed as follows:
  R 2 = 1 i = 1 n ( P i O i ) 2 i = 1 n O i 2
  RMSE = 1 n i = 1 n ( P i O i ) 2
  MAD = i = 1 n | P i O i | n

3. Results

3.1. Simulating GPP with GCC

The seasonal GPP dynamics estimated from GCC are shown in Figure 8 and Figure A1. The total growing season GPP was 790.48 g C·m−2·day−1 in 2017 and 809.24 g C·m−2·day−1 in 2018, respectively. Generally, GPP_CInorm showed the same seasonal dynamics with the GPP from eddy covariance in 2017 and 2018. The R2, MAD, and RMSE of GPP_CInorm and the GPP from eddy covariance in the growing seasons of 2017 were 0.68, 2.11 g C·m−2·day−1, and 2.47 g C·m−2·day−1, respectively. In the year 2018, R2, MAD, and RMSE between GPP_CInorm and GPP from eddy covariance were 0.78, 1.63 g C·m−2·day−1, and 1.94 g C·m−2·day−1, respectively. There was a systematic deviation between GPP from GCC and GPP from EC in spring and autumn, and these phenomena appeared both in 2017 and 2018. In 2017, the GPP_CInorm had the highest accuracy in autumn, followed by summer, and spring was the worst in the LUE model (Table 3).
In 2018, the GPP_CInorm had the highest accuracy in summer, followed by autumn, and spring was the worst in the LUE model (Table 3). More information about the changes of GPP from eddy covariance and GPP_CInorm in 2017 and 2018 can be found in the appendix (Appendix A, Figure A1).

3.2. Simulating GPP with MOD09GA VIs

The comparison between GPP from MODIS NDVI and GPP from eddy covariance in 2017 and 2018 at the AR station are shown in Figure 9a,b. In the growing season of 2017, the R2, RMSE, and MAD between GPP_NDVI and GPP from eddy covariance were 0.36, 3.30 g C·m−2·day−1, and 2.46 g C·m−2·day−1 respectively. In 2017 growing season, the GPP_NDVI had the highest accuracy in autumn, followed by spring and summer (Table 3). The R2, RMSE, and MAD between GPP_NDVI and GPP from eddy covariance were 0.46, 3.29 g C·m−2·day−1, and 2.45 g C·m−2·day−1 in the growing season of 2018, respectively. In 2018, the GPP_NDVI had the highest accuracy in autumn, followed by summer and spring (Table 3).
The comparison between GPP from MODIS EVI and GPP from eddy covariance in 2017 and 2018 at the AR station are shown in Figure 9c,d. In the 2017 growing season, the R2, RMSE, and MAD of GPP_EVI and GPP from eddy covariance were 0.53, 2.84 g C·m−2·day−1, and 2.19 g C·m−2·day−1, respectively. The LUE model with EVI as input (GPP_EVI) had the highest accuracy in spring, followed by autumn and summer (Table 3). During growing season in 2018, the R2, RMSE, and MAD of GPP_EVI and GPP from eddy covariance were 0.59, 2.85 g C·m−2·day−1, and 2.22 g C·m−2·day−1, respectively. The LUE model with EVI as input (GPP_EVI) performed best in autumn, followed by summer and spring (Table 3).
GPPs estimated from MODIS NDVI (GPP_NDVI) and EVI (GPP_EVI) were able to successfully capture the season dynamics of GPP from eddy covariance (Figure 9a–d). However, throughout the growing season, daily GPP_NDVI/GPP_EVI was dramatically underestimated at some dates comparing with GPP from eddy covariance, this was mainly caused by cloud or snow contamination of the satellite remote sensing data. The number of cloud contaminated data during the growing season in 2017 was 56, and 49 in 2018. The cumulative MAD caused by GPP_NDVI and GPP_EVI were 3.794 g C·m−2·day−1–3.849 g C·m−2·day−1 and 3.151 g C·m−2·day−1–3.193 g C·m−2·day−1, respectively. More information about the changes of GPP from eddy covariance and GPP_NDVI and GPP_EVI in two years (2017 and 2018) can be found in the appendix (Appendix A, Figure A2).

3.3. Simulating GPP with Multisource Data Using Machine Learning

In machine learning GPP estimation experiment, data in 2017 were used for training, and data in 2018 were used for validation. Four experiments were designed to test the contribution of different input datasets. The time period of the four GPP training experiments was from 26 April to 30 September, 2017. The time period of the four GPP validation experiments was from 26 April to 30 September, 2018. Figure 10a–d and Figure 11a–d are the results of training and validation in machine learning Exp. 1, Exp. 2, Exp. 3 and Exp. 4, respectively.
During the growing season, autumn had the highest training accuracy, followed by summer and spring corresponding to GPP_train1. For the GPP_train2, spring had the highest training accuracy, followed by autumn and summer. In the GPP_train3 experiment, spring had the highest training accuracy, followed by summer and autumn. In the GPP_train4 experiment, autumn had the highest training accuracy, followed by summer and spring (Table 4). In the four validation experiments (Figure 11), it is found that autumn had the highest validation accuracy for all of these experiments (GPP_valid1, GPP_valid2, GPP_valid3 and GPP_valid4), followed by summer and spring.

3.4. Comparing the Performance of Different GPP Estimation Methods

The statistical indicators of GPPs estimated by MOD09GA and GCC according to the LUE model and BNNA are shown in Table 5. The seasonal patterns between the GPP from eddy covariance and GPP from machine learning method were generally consistent, except for some large fluctuations between May and June. We compared the R2, RMSE, and MAD of the training results of GPP and found that the training accuracy of the GPP_train4 was the highest followed by GPP_train1, GPP_train3, and GPP_train2. The training accuracy of the GPP_train4 indicated that the training results can be improved by combining multisource data compared with the GPP_train1, GPP_train2, and GPP_train3 in machine learning. In addition, we found that the validation accuracy of GPP_valid4 was the highest followed by GPP_valid1, GPP_valid3, and GPP_valid2. Therefore, the effect of GPP_train4 and GPP_valid4 was the best of machine learning. Among the LUE models forced by different vegetation indices, the model accuracy from high to low were GPP_CInorm, GPP_EVI, and GPP_NDVI. When compared with the LUE models, the R2, MAD, and RMSE of machine learning were clearly better than those of the LUE models. Therefore, the machine learning method was an effective approach and can significantly improve the simulation accuracy of daily GPP as compared to the LUE model.

4. Discussion

4.1. GPP Estimation Based on GCC

The GCC can be used to quantitatively describe the vegetation canopy greenness in alpine meadows on the Tibetan Plateau, and the CInorm derived from GCC can be used to estimate the vegetation GPP. This was consistent with the finding that the GCC can quantify the total GPP of the grassland and deciduous broad-leaved forest [11]. In this study, we found that GPP_CInorm has the highest accuracy, followed by GPP_EVI, and GPP_NDVI was the worst among them. This consequence illustrates that long-term phenological photos can replace remote sensing data thereby obtaining higher-precision regional GPP. In addition, we found that there was no deviation in the GCC peak value obtained from the AR station, but GCC had systematic deviation in spring and autumn. However, previous studies reported that the GCC peak date in growing season is much different from observed GPP in forests and grasslands [21]. The GCC peak value of the phenology photo of the TKY site (deciduous broad-leaved forest) and the TKC site (evergreen coniferous forest) in Japan appeared about 20–30 days before the GPP peaks [26]. GCC and GPP of the Tongyu (meadow steppe) site in the Songnen Plain of China and the Maodeng Ranch (typical steppe) site in Inner Mongolia showed consistent and noticeable seasonal changes. GCC and GPP had large values during the growing season, and the values during the dormant period were very close to zero. The peak of GCC of meadow steppe lagged behind the peak of GPP by about ten days, and the peak of GCC of typical steppe lagged behind by GPP by about 20 days [44].
The advantages of GPP estimated by GCC include high continuity of phenology photos and that phenology photos can be obtained multiple times a day with less affects from clouds. In addition, GCC will be affected by some factors when estimating GPP, such as the exposure rate of the photo, the spatial scale of the photo, the height angle of the sun, the presence of shadows or animals in the photo, temperature, and photosynthesis. These factors will affect the values of GCC and thus GPP estimation. To our knowledge, there have been few studies on the use of phenological photos (GCC) to estimate GPP on the Tibetan Plateau. Therefore, the performance of GCC in GPP estimation for alpine meadow on the Tibetan Plateau was tested in this study.

4.2. GPP Estimation Based on Remote Sensing Vegetation Indices

The remote sensing data had great advantages of wide coverage and high spatiotemporal resolution, which can be used to estimate GPP for regional and even global over a long period of time. However, the complexity of cloud cover, snow, moisture, temperature, and vegetation growth will all affect the accuracy of GPP estimation by remote sensing [45]. We used the remote sensing vegetation index MOD09GA NDVI and EVI to estimate the GPPs of the AR station on an annual scale. For NDVI, we found that R2 was 0.50 and RMSE was 2.78 g C·m−2·day−1 in 2017; R2 was 0.57 and RMSE was 2.76 g C·m−2·day−1 in 2018. For EVI, R2 was 0.85 and RMSE was 1.65 g C·m−2·day−1 in 2017; R2 was 0.87 and RMSE was 1.56 g C m−2day−1 in 2018 (Appendix A, Figure A2). Li estimated the GPPs of AR station in 2009 and reported that R2 of LUE models were up to 0.8 at AR station at annual scale [46]. This is consistent with the performance of the LUE model in this study. This study found that at the AR station, using EVI to estimate GPP was better than using NDVI. This was also reported by a previous study [47]. Compared with NDVI, the blue band was used to correct the background effects when calculating EVI. This difference may be one of the reasons that EVI had a better correlation with GPP than NDVI [48]. These findings were consistent with conclusions that the remote sensing GPP model can accurately estimate GPP in alpine meadows.
The GPP_NDVI and GPP_EVI estimated by the LUE model based on the MOD09GA vegetation indices fluctuated greatly during the growing seasons, which were lower than the GPP from eddy covariance (Figure 9). GPP_NDVI and GPP_EVI had many cloud contaminated data points close to zero during the growing season, which may be caused by some noise points in the MOD09GA product. For GPP_NDVI, the number of cloud-contaminated data in the growing season of 2017 year was 56, R2 of GPP_NDVI and GPP from eddy covariance was 0.22, the number of without cloud contaminated data was 92, and the R2 of GPP_NDVI and GPP from eddy covariance was 0.61. The number of cloud contaminated data in the growing season of 2018 year was 49, the R2 of GPP_NDVI and GPP from eddy covariance was 0.26, the number of without cloud contaminated data was 99, and the R2 of GPP_NDVI and GPP from eddy covariance was 0.60. For GPP_EVI, the number of cloud contaminated data in the growing season of 2017 year was 48, R2 of GPP_EVI and GPP from eddy covariance was 0.40, the number of without cloud contaminated data was 100, and the R2 of GPP_EVI and GPP from eddy covariance was 0.64. The number of cloud contaminated data in the growing season of 2018 year was 46, the R2 of GPP_EVI and GPP from eddy covariance was 0.43, the number of without cloud contaminated data was 102, and the R2 of GPP_EVI and GPP from eddy covariance was 0.71. For the non-growth seasons, because all the alpine grasses were withered and photosynthesis was weak, GPP_NDVI and GPP_EVI were very small and close to zero, which was consistent with the GPP from eddy covariance (Appendix A, Figure A2).

4.3. The Impact of εmax for GPP Estimation

The maximum light energy utilization rate (εmax) of the remote sensing GPP model is an important parameter and has an important impact on the accuracy of GPP modeling. Turner revealed that the MODIS GPP algorithm may cause εmax to be underestimated [49]. MODIS GPP products in the midwestern U.S. croplands were underestimated due to the inaccuracy of εmax estimation. Therefore, the LUE values measured in the field should be used consistently in large-scale GPP modeling [50]. Lin conducted a monthly parameterization of the one-leaf LUE model and the two-leaf LUE model of eighteen typical plant function types. The results indicated that the εmax of most plant function types change with the season and had a significant positive correlation with the leaf area index and environmental temperature changes [17]. Wei developed a data-driven RFR-GPP model using random forest regression method, which showed LUE varied largely with latitude [51]. Zheng examined different sources of uncertainty in GPP simulated by LUE models and found that model structures (e.g., FPAR, Ws, and Ts) influenced both model parameters and model performance [52]. The εmax in alpine meadows at different altitudes on the Tibetan Plateau estimated to be 1.83 g C/MJ and 1.33 g C/MJ [47,53]. The εmax of this study was equal to 1.6 g C/MJ in alpine meadows [39]. εmax also vary with spatial and even in different year. Using a constant εmax for a certain biome may cause errors in LUE-oriented GPP modeling [54]. Therefore, obtaining reasonable εmax is vital to LUE models.

4.4. BNNA Performance Evaluation for GPP Estimation

In recent years, machine learning methods have been widely used in GPP estimation for its flexibility in application. In this work, it was found that machine learning algorithm (BNNA) performed well in alpine meadow GPP prediction. The machine learning algorithm (BNNA) was more accurate than the LUE models. The advantages of using machine learning method to fuse multisource data to estimate GPP were that machine learning methods can easily fuse data from different sources together and can combine the advantages of different data for GPP estimation. However, the main shortages were the lack of mechanism process, the need for a large number of training samples, and the results of training in different regions cannot be applied to other regions. The implication of this study was to provide a new idea to improve the GPP estimation accuracy on the Tibetan Plateau. With the large number of phenological cameras, the method will have good application potential on the Tibetan Plateau, which is greatly affected by clouds and snow.

5. Conclusions

The GPP_CInorm has the highest accuracy, followed by GPP_EVI, and GPP_NDVI was the worst in the LUE model. However, it should be noted that there is a problem with the seasonal dynamics of CInorm results, manifested as a significant deviation existing in spring and autumn. In several sets of machine learning experiments, the accuracy of GPP was found to be the highest when both the CInorm based on phenological camera photo and MODIS vegetation indices (NDVI and EVI) are input at the same time, indicating that the two types of data can complement each other’s defects and improve estimation accuracy of GPP. Compared with the GPP accuracy estimated by EVI and NDVI in the LUE model, the GPP accuracy estimated by combining multisource data is increased by 27% and 59%, respectively. These results illustrate that the overall accuracy of machine learning is significantly higher than the traditional LUE model. However, compared with the GPP accuracy estimated by CInorm in the LUE model, it is reduced by 3%, which is caused by the lower quality of GPP data from eddy covariance and the shorter time series of GCC data. Due to the late installation of phenological camera resulting in the quantity of photos available now is very small. Thus, we can get longer time series GCC data later, the purpose is to reduce the impact of the amount of data on the accuracy of machine learning estimation GPP.

Author Contributions

X.W. designed the research; X.Z., Y.Z., X.B., and X.W. collected and processed data; X.Z. and X.W. analyzed the data and wrote the draft; S.Z. provided import input to the manuscript; X.Z. and X.W. revised the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key R&D Program of China (Grant No. 2017YFA0604801), the National Natural Science Foundation of China (Grant No. 41771466), the Special Fund for Key Program of Science and Technology of Qinghai Province (Grant No. 2017-SF-A6), and the Youth Innovation Promotion Association CAS to X.W. (No. 2020422).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. The annual changes of the GPP_CInorm and the GPP from eddy covariance, time from (a) 26 April to 16 December 2017, and (b) from 5 January to 31 December 2018.
Figure A1. The annual changes of the GPP_CInorm and the GPP from eddy covariance, time from (a) 26 April to 16 December 2017, and (b) from 5 January to 31 December 2018.
Remotesensing 12 03735 g0a1
Figure A2. The year-round changes of GPP_NDVI in 2017 (a), GPP_NDVI in 2018 (b), GPP_EVI in 2017 (c), GPP_EVI in 2018 (d) along with the GPP from eddy covariance in the same year.
Figure A2. The year-round changes of GPP_NDVI in 2017 (a), GPP_NDVI in 2018 (b), GPP_EVI in 2017 (c), GPP_EVI in 2018 (d) along with the GPP from eddy covariance in the same year.
Remotesensing 12 03735 g0a2

References

  1. Monteith, J.L. Solar Radiation and Productivity in Tropical Ecosystems. J. Appl. Ecol. 1972, 9, 747–766. [Google Scholar] [CrossRef] [Green Version]
  2. Wang, S.; Huang, K.; Yan, H.; Yan, H.; Zhou, L.; Wang, H.; Zhang, J.; Yan, J.; Zhao, L.; Wang, Y.; et al. Improving the light use efficiency model for simulating terrestrial vegetation gross primary production by the inclusion of diffuse radiation across ecosystems in China. Ecol. Complex. 2015, 23, 1–13. [Google Scholar] [CrossRef]
  3. Jia, W.; Liu, M.; Wang, D.; He, H.; Shi, P.; Li, Y.; Wang, Y. Uncertainty in simulating regional gross primary productivity from satellite-based models over northern China grassland. Ecol. Indic. 2018, 88, 134–143. [Google Scholar] [CrossRef]
  4. Lobell, D.B.; Hicke, J.A.; Asner, G.P.; Field, C.B.; Tucker, C.J.; Los, S.O. Satellite estimates of productivity and light use efficiency in United States agriculture, 1982–1998. Glob. Chang. Biol 2002, 8, 722–735. [Google Scholar] [CrossRef]
  5. Yu, G.-R.; Zhu, X.-J.; Fu, Y.-L.; He, H.-L.; Wang, Q.-F.; Wen, X.-F.; Li, X.-R.; Zhang, L.-M.; Zhang, L.; Su, W.; et al. Spatial patterns and climate drivers of carbon fluxes in terrestrial ecosystems of China. Glob. Chang. Biol. 2013, 19, 798–810. [Google Scholar] [CrossRef]
  6. Zhang, W.J.; Wang, H.M.; Yang, F.T.; Yi, Y.H.; Wen, X.F.; Sun, X.M.; Yu, G.R.; Wang, Y.D.; Ning, J.C. Underestimated effects of low temperature during early growing season on carbon sequestration of a subtropical coniferous plantation. Biogeosciences 2011, 8, 1667–1678. [Google Scholar] [CrossRef] [Green Version]
  7. Mekonnen, Z.A.; Grant, R.F.; Schwalm, C. Contrasting changes in gross primary productivity of different regions of North America as affected by warming in recent decades. Agric. For. Meteorol. 2016, 218–219, 50–64. [Google Scholar] [CrossRef] [Green Version]
  8. Jahan, N.; Gan, T.Y. Developing a gross primary production model for coniferous forests of northeastern USA from MODIS data. Int. J. Appl. Earth Obs. Geoinf. 2013, 25, 11–20. [Google Scholar] [CrossRef]
  9. Xiao, J.; Zhuang, Q.; Baldocchi, D.D.; Law, B.E.; Richardson, A.D.; Chen, J.; Oren, R.; Starr, G.; Noormets, A.; Ma, S.; et al. Estimation of net ecosystem carbon exchange for the conterminous United States by combining MODIS and AmeriFlux data. Agric. For. Meteorol. 2008, 148, 1827–1847. [Google Scholar] [CrossRef] [Green Version]
  10. Wu, C.; Niu, Z.; Gao, S. Gross primary production estimation from MODIS data with vegetation index and photosynthetically active radiation in maize. J. Geophys. Res. 2010, 115. [Google Scholar] [CrossRef] [Green Version]
  11. Toomey, M.; Friedl, M.A.; Frolking, S.; Hufkens, K.; Klosterman, S.; Sonnentag, O.; Baldocchi, D.D.; Bernacchi, C.J.; Biraud, S.C.; Bohrer, G.; et al. Greenness indices from digital cameras predict the timing and seasonal dynamics of canopy-scale photosynthesis. Ecol. Appl. 2015, 25, 99–115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Potter, C.S.; Randerson, J.T.; Field, C.B.; Matson, P.A.; Vitousek, P.M.; Mooney, H.A.; Klooster, S.A. Terrestrial ecosystem production: A process model based on global satellite and surface data. Glob. Biogeochem. Cycle 1993, 7, 811–841. [Google Scholar] [CrossRef]
  13. Prince, S.D.; Goward, S.N. Global Primary Production: A Remote Sensing Approach. J. Biogeogr. 1995, 22, 815–835. [Google Scholar] [CrossRef]
  14. 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]
  15. 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]
  16. Sjöström, M.; Ardö, J.; Arneth, A.; Boulain, N.; Cappelaere, B.; Eklundh, L.; de Grandcourt, A.; Kutsch, W.L.; Merbold, L.; Nouvellon, Y.; et al. Exploring the potential of MODIS EVI for modeling gross primary production across African ecosystems. Remote Sens. Environ. 2011, 115, 1081–1089. [Google Scholar] [CrossRef]
  17. Lin, X.; Chen, B.; Chen, J.; Zhang, H.; Sun, S.; Xu, G.; Guo, L.; Ge, M.; Qu, J.; Li, L.; et al. Seasonal fluctuations of photosynthetic parameters for light use efficiency models and the impacts on gross primary production estimation. Agric. For. Meteorol. 2017, 236, 22–35. [Google Scholar] [CrossRef] [Green Version]
  18. Montero, P.; Daneri, G.; González, H.E.; Iriarte, J.L.; Tapia, F.J.; Lizárraga, L.; Sanchez, N.; Pizarro, O. Seasonal variability of primary production in a fjord ecosystem of the Chilean Patagonia: Implications for the transfer of carbon within pelagic food webs. Cont. Shelf Res. 2011, 31, 202–215. [Google Scholar] [CrossRef]
  19. Richardson, A.D.; Braswell, B.H.; Hollinger, D.Y.; Jenkins, J.P.; Ollinger, S.V.J.E.A. Near-surface remote sensing of spatial and temporal variation in canopy phenology. Ecol. Appl. 2009, 19, 1417–1428. [Google Scholar] [CrossRef]
  20. Sonnentag, O.; Hufkens, K.; Teshera-Sterne, C.; Young, A.M.; Friedl, M.; Braswell, B.H.; Milliman, T.; O’Keefe, J.; Richardson, A.D.J.A.; Meteorology, F. Digital repeat photography for phenological research in forest ecosystems. Agric. For. Meteorol. 2012, 152, 159–177. [Google Scholar] [CrossRef]
  21. Delpierre, N.; Soudani, K.; Berveiller, D.; Dufrêne, E.; Hmimina, G.; Vincent, G. “Green pointillism”: Detecting the within-population variability of budburst in temperate deciduous trees with phenological cameras. Int. J. Biometeorol. 2020, 64, 663–670. [Google Scholar] [CrossRef] [PubMed]
  22. Alberton, B.; Torres, R.d.S.; Cancian, L.F.; Borges, B.D.; Almeida, J.; Mariano, G.C.; Santos, J.D.; Morellato, L.P.C. Introducing digital cameras to monitor plant phenology in the tropics: Applications for conservation. Perspect. Ecol. Conser 2017, 15, 82–90. [Google Scholar] [CrossRef]
  23. Westergaard-Nielsen, A.; Lund, M.; Hansen, B.U.; Tamstorf, M.P. Camera derived vegetation greenness index as proxy for gross primary production in a low Arctic wetland area. ISPRS J. Photogramm. Remote Sens 2013, 86, 89–99. [Google Scholar] [CrossRef]
  24. Wang, H.; Jia, G.; Epstein, H.E.; Zhao, H.; Zhang, A. Integrating a PhenoCam-derived vegetation index into a light use efficiency model to estimate daily gross primary production in a semi-arid grassland. Agric. For. Meteorol. 2020, 288–289, 107983. [Google Scholar] [CrossRef]
  25. Zhao, J.; Zhang, Y.; Tan, Z.; Song, Q.; Liang, N.; Yu, L.; Zhao, J. Using digital cameras for comparative phenological monitoring in an evergreen broad-leaved forest and a seasonal rain forest. Ecol. Infor. 2012, 10, 65–72. [Google Scholar] [CrossRef]
  26. Saitoh, T.M.; Nagai, S.; Saigusa, N.; Kobayashi, H.; Suzuki, R.; Nasahara, K.N.; Muraoka, H. Assessing the use of camera-based indices for characterizing canopy phenology in relation to gross primary production in a deciduous broad-leaved and an evergreen coniferous forest in Japan. Ecol. Infor. 2012, 11, 45–54. [Google Scholar] [CrossRef]
  27. Wang, M.; Sun, R.; Zhu, A.; Xiao, Z. Evaluation and Comparison of Light Use Efficiency and Gross Primary Productivity Using Three Different Approaches. Remote Sens. 2020, 12, 1003. [Google Scholar] [CrossRef] [Green Version]
  28. Wolanin, A.; Camps-Valls, G.; Gómez-Chova, L.; Mateo-García, G.; van der Tol, C.; Zhang, Y.; Guanter, L. Estimating crop primary productivity with Sentinel-2 and Landsat 8 using machine learning methods trained with radiative transfer simulations. Remote Sens. Environ. 2019, 225, 441–457. [Google Scholar] [CrossRef]
  29. Yang, F.; Ichii, K.; White, M.A.; Hashimoto, H.; Michaelis, A.R.; Votava, P.; Zhu, A.X.; Huete, A.; Running, S.W.; Nemani, R.R. Developing a continental-scale measure of gross primary production by combining MODIS and AmeriFlux data through Support Vector Machine approach. Remote Sens. Environ. 2007, 110, 109–122. [Google Scholar] [CrossRef]
  30. Tramontana, G.; Ichii, K.; Camps-Valls, G.; Tomelleri, E.; Papale, D. Uncertainty analysis of gross primary production upscaling using Random Forests, remote sensing and eddy covariance data. Remote Sens. Environ. 2015, 168, 360–373. [Google Scholar] [CrossRef]
  31. Lee, B.; Kim, E.; Lim, J.-H.; Kang, M.; Kim, J. Application of Machine Learning Algorithm and Remote-sensed Data to Estimate Forest Gross Primary Production at Multi-sites Level. Korean J. Remote Sens 2019, 35, 1117–1132. [Google Scholar] [CrossRef]
  32. Du, M.; Kawashima, S.; Yonemura, S.; Zhang, X.; Chen, S. Mutual influence between human activities and climate change in the Tibetan Plateau during recent years. Glob. Planet. Chang. 2004, 41, 241–249. [Google Scholar] [CrossRef]
  33. Zheng, Z.; Zhu, W.; Chen, G.; Jiang, N.; Fan, D.; Zhang, D. Continuous but diverse advancement of spring-summer phenology in response to climate warming across the Qinghai-Tibetan Plateau. Agric. For. Meteorol. 2016, 223, 194–202. [Google Scholar] [CrossRef]
  34. An, S.; Chen, X.; Zhang, X.; Lang, W.; Ren, S.; Xu, L. Precipitation and Minimum Temperature are Primary Climatic Controls of Alpine Grassland Autumn Phenology on the Qinghai-Tibet Plateau. Remote Sens. 2020, 12, 431. [Google Scholar] [CrossRef] [Green Version]
  35. Huete, A.R.; Liu, H.Q.; Batchily, K.; van Leeuwen, W. A comparison of vegetation indices over a global set of TM images for EOS-MODIS. Remote Sens. Environ. 1997, 59, 440–451. [Google Scholar] [CrossRef]
  36. Shokr, M. Potential directions for applications of satellite earth observations data in Egypt. Egyptian. J. Remote Sens. Space Sci. 2011, 14, 1–13. [Google Scholar] [CrossRef] [Green Version]
  37. Liu, S.M.; Xu, Z.W.; Wang, W.Z.; Jia, Z.Z.; Zhu, M.J.; Bai, J.; Wang, J.M. A comparison of eddy-covariance and large aperture scintillometer measurements with respect to the energy balance closure problem. Hydrol. Earth Syst. Sci. 2011, 15, 1291–1306. [Google Scholar] [CrossRef] [Green Version]
  38. Gilmanov, T.G.; Soussana, J.F.; Aires, L.; Allard, V.; Ammann, C.; Balzarolo, M.; Barcza, Z.; Bernhofer, C.; Campbell, C.L.; Cernusca, A.; et al. Partitioning European grassland net ecosystem CO2 exchange into gross primary productivity and ecosystem respiration using light response function analysis. Agric. Ecosyst. Environ. 2007, 121, 93–120. [Google Scholar] [CrossRef]
  39. 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]
  40. Ahrends, H.E.; Etzold, S.; Kutsch, W.L.; Stoeckli, R.; Bruegger, R.; Jeanneret, F.; Wanner, H.; Buchmann, N.; Eugster, W. Tree phenology and carbon dioxide fluxes: Use of digital photography for process-based interpretation at the ecosystem scale. Clim. Res. 2009, 39, 261–274. [Google Scholar] [CrossRef] [Green Version]
  41. Knox, S.H.; Dronova, I.; Sturtevant, C.; Oikawa, P.Y.; Matthes, J.H.; Verfaillie, J.; Baldocchi, D. Using digital camera and Landsat imagery with eddy covariance data to model gross primary production in restored wetlands. Agric. For. Meteorol. 2017, 237, 233–245. [Google Scholar] [CrossRef] [Green Version]
  42. Zhou, R.; Zhang, Y.; Song, Q.; Lin, Y.; Sha, L.; Jin, Y.; Liu, Y.; Fei, X.; Gao, J.; He, Y.; et al. Relationship between gross primary production and canopy colour indices from digital camera images in a rubber (Hevea brasiliensis) plantation, Southwest China. Forest. Ecol. Manag. 2019, 437, 222–231. [Google Scholar] [CrossRef]
  43. Heinsch, F.A.; Reeves, M.; Votava, P.; Kang, S.; Milesi, C.; Zhao, M.S. User’s Guide GPP and NPP (MOD17A2/A3) Products NASA MODIS Land Algorithm, Version 2.0. 2002. Available online: http://www.ntsg.umt.edu/modis/MOD17UsersGuide.pdf (accessed on 24 December 2008).
  44. Zhao, H.; Jia, G.; Wang, H.; Zhang, A.; Xu, X. Seasonal and interannual variations in carbon fluxes in East Asia semi-arid grasslands. Sci. Total Environ. 2019, 668, 1128–1138. [Google Scholar] [CrossRef] [PubMed]
  45. Chai, X.; Shi, P.; Zong, N.; He, Y.; Zhang, X.; Xu, M.; Zhang, J. A growing season climatic index to simulate gross primary productivity and carbon budget in a Tibetan alpine meadow. Ecol. Indic. 2017, 81, 285–294. [Google Scholar] [CrossRef]
  46. Li, F.; Wang, X.F.; Zhao, J.; Zhang, X.Q.; Zhao, Q.J. A method for estimating the gross primary production of alpine meadows using MODIS and climate data in China. Int. J. Remote Sens. 2013, 34, 8280–8300. [Google Scholar] [CrossRef]
  47. Li, Z.; Yu, G.; Xiao, X.; Li, Y.; Zhao, X.; Ren, C.; Zhang, L.; Fu, Y. Modeling gross primary production of alpine ecosystems in the Tibetan Plateau using MODIS images and climate data. Remote Sens. Environ. 2007, 107, 510–519. [Google Scholar] [CrossRef]
  48. 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]
  49. Turner, D.P.; Ritts, W.D.; Zhao, M.; Kurc, S.A.; Running, S.W.J.I.T.o.G.; Sensing, R. Assessing Interannual Variation in MODIS-Based Estimates of Gross Primary Production. IEEE Trans. Geosci. Remote Sens. 2006, 44, 1899–1907. [Google Scholar] [CrossRef]
  50. Xin, Q.; Broich, M.; Suyker, A.E.; Yu, L.; Gong, P. Multi-scale evaluation of light use efficiency in MODIS gross primary productivity for croplands in the Midwestern United States. Agric. For. Meteorol. 2015, 201, 111–119. [Google Scholar] [CrossRef] [Green Version]
  51. Wei, S.; Yi, C.; Wei, F.; Hendrey, G. A global study of GPP focusing on light-use efficiency in a random forest regression model. Ecosphere 2017, 8. [Google Scholar] [CrossRef]
  52. Zheng, Y.; Zhang, L.; Xiao, J.; Yuan, W.; Yan, M.; Li, T.; Zhang, Z. Sources of uncertainty in gross primary productivity simulated by light use efficiency models: Model structure, parameters, input data, and spatial resolution. Agric. For. Meteorol. 2018, 263, 242–257. [Google Scholar] [CrossRef]
  53. Xu, L.; Zhang, X.; Shi, P.; Yu, G. Response of canopy quantum yield of alpine meadow to temperature under low atmospheric pressure on Tibetan Plateau. Sci. China Ser. D Earth Sci. 2006, 49, 219–225. [Google Scholar] [CrossRef]
  54. Wang, H.; Jia, G.; Fu, C.; Feng, J.; Zhao, T.; Ma, Z. Deriving maximal light use efficiency from coordinated flux measurements and satellite data for regional gross primary production modeling. Remote Sens. Environ. 2010, 114, 2248–2258. [Google Scholar] [CrossRef]
Figure 1. Location of the A’rou (AR) station and vegetation distribution on the Tibetan Plateau.
Figure 1. Location of the A’rou (AR) station and vegetation distribution on the Tibetan Plateau.
Remotesensing 12 03735 g001
Figure 2. (a) RG, (b) VPD, (c) Tmin and (d) PAR of the A’rou (AR) station on the Tibetan plateau from 2017 to 2018.
Figure 2. (a) RG, (b) VPD, (c) Tmin and (d) PAR of the A’rou (AR) station on the Tibetan plateau from 2017 to 2018.
Remotesensing 12 03735 g002
Figure 3. (a) NDVI and (b) EVI of the A’rou (AR) station on the Tibetan plateau from 2017 to 2018.
Figure 3. (a) NDVI and (b) EVI of the A’rou (AR) station on the Tibetan plateau from 2017 to 2018.
Remotesensing 12 03735 g003
Figure 4. Example photos at different moments of the day by the phenological camera at A’rou (AR) station on 22 July 2018.
Figure 4. Example photos at different moments of the day by the phenological camera at A’rou (AR) station on 22 July 2018.
Remotesensing 12 03735 g004
Figure 5. Changes of ROI_GCC1 (a,b) and ROI_GCC2 (c,d) with GPP from eddy covariance, time from (a) 15 April to 16 December, 2017, (b) from 5 January to 31 December 2018, (c) from 15 April to 16 December 2017, and (d) from 5 Januaryto 31 December 2018.
Figure 5. Changes of ROI_GCC1 (a,b) and ROI_GCC2 (c,d) with GPP from eddy covariance, time from (a) 15 April to 16 December, 2017, (b) from 5 January to 31 December 2018, (c) from 15 April to 16 December 2017, and (d) from 5 Januaryto 31 December 2018.
Remotesensing 12 03735 g005
Figure 6. Structure of the BP neural network algorithm.
Figure 6. Structure of the BP neural network algorithm.
Remotesensing 12 03735 g006
Figure 7. Flow diagram of data processing and analyses.
Figure 7. Flow diagram of data processing and analyses.
Remotesensing 12 03735 g007
Figure 8. The daily dynamics of the GPP_CInorm and the GPP from eddy covariance during the growing seasons (1 May to 1 October), (a) 2017, (b) 2018.
Figure 8. The daily dynamics of the GPP_CInorm and the GPP from eddy covariance during the growing seasons (1 May to 1 October), (a) 2017, (b) 2018.
Remotesensing 12 03735 g008
Figure 9. The daily dynamics of the GPP_NDVI, GPP_EVI and the GPP from eddy covariance during the growing seasons (1 May to 1 October), (a) 2017, (b) 2018, (c) 2017, and (d) 2018.
Figure 9. The daily dynamics of the GPP_NDVI, GPP_EVI and the GPP from eddy covariance during the growing seasons (1 May to 1 October), (a) 2017, (b) 2018, (c) 2017, and (d) 2018.
Remotesensing 12 03735 g009
Figure 10. Training results of the four experiments for the BNNA method. (a), (b), (c) and (d) represent the predicted GPP and the GPP from eddy covariance in GPP_train1, GPP_train2, GPP_train3, and GPP_train4, respectively.
Figure 10. Training results of the four experiments for the BNNA method. (a), (b), (c) and (d) represent the predicted GPP and the GPP from eddy covariance in GPP_train1, GPP_train2, GPP_train3, and GPP_train4, respectively.
Remotesensing 12 03735 g010
Figure 11. Validation of the four experiments for BNNA method. (a), (b), (c), and (d) represent the year-round changes of GPP_valid and the GPP from eddy covariance in GPP_valid1, GPP_valid2, GPP_valid3, and GPP_valid4, respectively.
Figure 11. Validation of the four experiments for BNNA method. (a), (b), (c), and (d) represent the year-round changes of GPP_valid and the GPP from eddy covariance in GPP_valid1, GPP_valid2, GPP_valid3, and GPP_valid4, respectively.
Remotesensing 12 03735 g011
Table 1. The parameter settings in the BP neural network algorithm.
Table 1. The parameter settings in the BP neural network algorithm.
Parameters and ValuesExplanatory Notes
epochs = 100Number of iterations
batch_size = 6Number of feeding data for each training model
learnrate = 0.1Learning rate
Activation = “tanh”Activation function
Terror = 1 × 10−7Iteration stop condition
Loss = ‘mean_squared_error’Loss function
Dropout = 0.2Rejection rate of each layer
model Hierarchical model
Table 2. Machine learning-based GPP estimation experiments.
Table 2. Machine learning-based GPP estimation experiments.
ExperimentsInput LayerOutput Layer
Exp 1f(VPD), PAR, f(Tmin), CInormGPP
Exp 2f(VPD), PAR, f(Tmin), NDVIGPP
Exp 3f(VPD), PAR, f(Tmin), EVIGPP
Exp 4f(VPD), PAR, f(Tmin), CInorm, NDVI, EVIGPP
GPP refers to the GPP from eddy covariance.
Table 3. Statistical indicators of the GPP from eddy covariance with the GPP simulation values of remote sensing GPP model in spring, summer, and autumn.
Table 3. Statistical indicators of the GPP from eddy covariance with the GPP simulation values of remote sensing GPP model in spring, summer, and autumn.
GPPR2RMSEMAD
GPP_CInorm 201720182017201820172018
Spring0.530.381.561.161.240.95
Summer0.550.752.772.072.441.80
Autumn0.700.712.552.162.201.88
GPP_NDVISpring0.320.241.121.340.851.08
Summer0.290.343.804.072.943.12
Autumn0.360.493.523.102.842.46
GPP_EVISpring0.510.380.881.250.610.99
Summer0.420.493.183.462.572.82
Autumn0.490.573.142.762.612.23
GPP refers to the GPP from eddy covariance. The seasons on the Tibetan Plateau, with spring as May, summer as June to July, and autumn as August to September.
Table 4. Statistical indicators of the GPP from eddy covariance with the GPP simulation of the remote sensing GPP model in spring, summer, and autumn.
Table 4. Statistical indicators of the GPP from eddy covariance with the GPP simulation of the remote sensing GPP model in spring, summer, and autumn.
ExperimentSeasonR2RMSEMAD
TrainingValidationTrainingValidationTrainingValidation
1Spring0.170.151.902.221.641.90
Summer0.640.581.501.611.191.27
Autumn0.660.781.481.311.121.03
2Spring0.530.292.622.182.381.83
Summer0.420.351.801.901.451.59
Autumn0.430.691.661.481.351.16
3Spring0.600.372.101.671.871.35
Summer0.480.411.691.861.391.56
Autumn0.470.651.631.631.261.26
4Spring0.360.301.371.581.201.33
Summer0.660.581.511.621.231.30
Autumn0.660.701.431.411.001.09
1, 2, 3, and 4 refer to Exp. 1, Exp. 2, Exp. 3, and Exp. 4, respectively. The seasons on the Tibetan Plateau, with spring as May, summer as June to July, and autumn as August to September. Training refers to the training period (from 26 April to 30 September 2017), Validation refers to the validation period (from 26 April to 30 September 2018).
Table 5. Statistical indicators of the GPP from eddy covariance with the GPP simulation of the remote sensing GPP model and machine learning from 26 April to 30 September 2017 and from 26 April to 30 September 2018.
Table 5. Statistical indicators of the GPP from eddy covariance with the GPP simulation of the remote sensing GPP model and machine learning from 26 April to 30 September 2017 and from 26 April to 30 September 2018.
Date/ExperimentGPPR2 MADRMSE
LUE model2017GPP_CInorm0.692.062.46
2018GPP_CInorm0.781.641.95
2017GPP_NDVI0.372.453.29
2018GPP_NDVI0.472.413.28
2017GPP_EVI0.542.182.83
2018GPP_EVI0.592.192.84
Machine Learning NormalizedTrainingGPP_train10.671.281.61
TrainingGPP_train20.541.611.96
TrainingGPP_train30.621.451.77
TrainingGPP_train40.731.151.47
ValidationGPP_valid10.731.341.69
ValidationGPP_valid20.661.491.83
ValidationGPP_valid30.691.411.74
ValidationGPP_valid40.751.251.56
2017 (Training) refers to a period (from 26 April to 30 September 2017), 2018 (Validation) refers to a period (from 26 April to 30 September 2018). GPP refers to the GPP from eddy covariance. Bold fields and values are the best performers in LUE model and Machine Learning, respectively.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhou, X.; Wang, X.; Zhang, S.; Zhang, Y.; Bai, X. Combining Phenological Camera Photos and MODIS Reflectance Data to Predict GPP Daily Dynamics for Alpine Meadows on the Tibetan Plateau. Remote Sens. 2020, 12, 3735. https://doi.org/10.3390/rs12223735

AMA Style

Zhou X, Wang X, Zhang S, Zhang Y, Bai X. Combining Phenological Camera Photos and MODIS Reflectance Data to Predict GPP Daily Dynamics for Alpine Meadows on the Tibetan Plateau. Remote Sensing. 2020; 12(22):3735. https://doi.org/10.3390/rs12223735

Chicago/Turabian Style

Zhou, Xuqiang, Xufeng Wang, Songlin Zhang, Yang Zhang, and Xuejie Bai. 2020. "Combining Phenological Camera Photos and MODIS Reflectance Data to Predict GPP Daily Dynamics for Alpine Meadows on the Tibetan Plateau" Remote Sensing 12, no. 22: 3735. https://doi.org/10.3390/rs12223735

APA Style

Zhou, X., Wang, X., Zhang, S., Zhang, Y., & Bai, X. (2020). Combining Phenological Camera Photos and MODIS Reflectance Data to Predict GPP Daily Dynamics for Alpine Meadows on the Tibetan Plateau. Remote Sensing, 12(22), 3735. https://doi.org/10.3390/rs12223735

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