Next Article in Journal
Analysis of the Influence of Deforestation on the Microphysical Parameters of Clouds in the Amazon
Next Article in Special Issue
Joint Characterization of Sentinel-2 Reflectance: Insights from Manifold Learning
Previous Article in Journal
Assessing Through-Water Structure-from-Motion Photogrammetry in Gravel-Bed Rivers under Controlled Conditions
Previous Article in Special Issue
The Impact of NPV on the Spectral Parameters in the Yellow-Edge, Red-Edge and NIR Shoulder Wavelength Regions in Grasslands
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Accuracy Assessment of Photochemical Reflectance Index (PRI) and Chlorophyll Carotenoid Index (CCI) Derived from GCOM-C/SGLI with In Situ Data

1
Graduate School of Science and Technology, University of Tsukuba, 1-1-1 Tennoudai, Tsukuba 305-8572, Ibaraki, Japan
2
Geological Survey of Japan, National Institute of Advanced Industrial Science and Technology, Central 7, 1-1-1 Higashi, Tsukuba 305-8567, Ibaraki, Japan
3
Earth Observation Research Center, Japan Aerospace Exploration Agency, 2-1-1 Sengen, Tsukuba 305-8505, Ibaraki, Japan
4
Earth System Division, National Institute for Environmental Studies, 16-2 Onogawa, Tsukuba 305-8506, Ibaraki, Japan
5
Teshio Experimental Forest, Field Science Center for Northern Biosphere, Hokkaido University, 131 Toikanbetsu, Horonobe 098-2943, Hokkaido, Japan
6
Kansai Research Center, Forestry and Forest Products Research Institute, 68 Nagai-kyutaro, Momoyama-cho, Fushimi-ku, Kyoto 612-0855, Japan
7
Uryu Experimental Forest, Field Science Center for Northern Biosphere, Hokkaido University, Moshiri, Horokanai, Uryu 074-0741, Hokkaido, Japan
8
Faculty of Life and Environmental Sciences, University of Tsukuba, 1-1-1 Tennoudai, Tsukuba 305-8572, Ibaraki, Japan
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(21), 5352; https://doi.org/10.3390/rs14215352
Submission received: 24 August 2022 / Revised: 8 October 2022 / Accepted: 17 October 2022 / Published: 26 October 2022
(This article belongs to the Special Issue Feature Papers for Section Biogeosciences Remote Sensing)

Abstract

:
The photochemical reflectance index (PRI) and the chlorophyll carotenoid index (CCI) are carotenoid-sensitive vegetation indices, which can monitor vegetation’s photosynthetic activities. One unique satellite named “Global Change Observation Mission-Climate (GCOM-C)” is equipped with a sensor, “Second Generation Global Imager (SGLI)”, which has the potential to frequently and simultaneously observe PRI and CCI over a wide swath. However, the observation accuracy of PRI and CCI derived from GCOM-C/SGLI remains unclear in forests. Thus, we demonstrated their accuracy assessment by comparing them with in situ data. We collected in situ spectral irradiance data at four forest sites in Japan for three years. We statistically compared satellite PRI with in situ PRI, and satellite CCI with in situ CCI. From the obtained results, the satellite PRI showed poor agreement (the best: r = 0.294 ( p < 0.05 )) and the satellite CCI showed good agreement (the best: r = 0.911 ( p < 0.001 )). The greater agreement of satellite CCI is possibly because satellite CCI contained fewer outliers and satellite CCI was more resistant to small noise, compared to satellite PRI. Our results suggest that the satellite CCI is more suitable for practical use than the satellite PRI with the latest version (version 3) of GCOM-C/SGLI’s products.

1. Introduction

The photochemical reflectance index (PRI) is a narrow-band vegetation index proposed by Gamon et al. [1,2] as a proxy of vegetation’s photosynthetic activities. PRI is affected by environmental factors such as soil moisture, precipitation, and air temperature [3,4,5,6,7,8,9,10,11,12,13,14]. It is defined as
P R I = ρ ( 531 ) ρ ( 570 ) ρ ( 531 ) + ρ ( 570 ) ,
where ρ ( λ ) represents the reflectance at λ   nm wavelength. The PRI is sensitive to the short-term (such as diurnal) changes in the composition of xanthophyll pigments, namely, zeaxanthin, antheraxanthin, and violaxanthin [15,16]. These pigments convert from one to another in the xanthophyll cycle, and the changes are closely related to the photosynthetic activity. On the other hand, the PRI is also sensitive to the long-term (such as seasonal) changes in the ratio between chlorophyll and carotenoid pigment pools [15,16,17,18]. Various studies have indicated relationships between PRI and light use efficiency (LUE) [19,20], which is the ratio between the gross primary production (GPP) and the absorbed photosynthetic active radiation (APAR) at both the leaf scale [21,22,23,24] and the canopy scale [25,26,27,28]. Garbulsky et al. [29] and Zhang et al. [30] have conducted meta-analyses and found significant relationships between the PRI and the LUE over the various species and the different spatiotemporal scales. At the ecosystem scale, the PRI has been observed by satellites, aircrafts, or unmanned aerial vehicles equipped with hyperspectral imagers [5,31,32,33,34]. However, in these observations, the PRI can be irregularly obtained with only a narrow swath. Thus, some alternative indices, which are regularly available over a wide swath with satellites, have been developed for the ecosystem-scale observation of PRI.
One of the widely used alternative indices is the chlorophyll carotenoid index (CCI), which is also proposed by Gamon et al. [35]. The CCI is defined from the reflectance at 531 nm and the reflectance at red color regions as
C C I = ρ ( 531 ) ρ ( red ) ρ ( 531 ) + ρ ( red ) .
For example, at the ecosystem scale, a sensor named Moderate Resolution Imaging Spectroradiometer (MODIS) carried on Terra and Aqua satellites has been used to obtain the CCI [35,36]. Middleton et al. [37] obtained the CCI from MODIS band 11 (531 nm) and band 1 (645 nm). Drolet et al. [36] used the combination of MODIS band 11 and band 14 (678 nm) to obtain the CCI. CCI has been used for monitoring the changes in the LUE [8,38] and phenology [12,39]. Although CCI was developed as an alternative index of PRI, several studies have reported the CCI has the advantage over the PRI. Springer et al. [8] and Wong et al. [40] found that CCI is more sensitive to GPP than LUE at the leaf and canopy scale. Therefore, more detailed analysis of PRI and CCI based on the regular and simultaneous observation of them at the ecosystem scale was needed to understand photosynthetic activities and phenology of forests [41].
Currently, both the PRI and the CCI can be regularly obtained over a wide swath with a new satellite: Global Change Observation Mission-Climate (GCOM-C), launched by Japan Aerospace Exploration Agency (JAXA) on 23 December 2017 [42,43]. GCOM-C is equipped with a sensor named Second-Generation Global Imager (SGLI). GCOM-C/SGLI observes with 1150 k m swath range for its visible and near-infrared (NIR) bands. The observation interval is once every two days around Japan (latitude = 35°N) with 14 days’ revisit time. GCOM-C/SGLI’s band5 (VN05: 529.7 nm), band6 (VN06: 566.1 nm), and band8 (VN08: 672.4 nm) can detect the reflectance at 531 nm, 570 nm, and red color regions, respectively. Hence, GCOM-C/SGLI has the potential to observe the PRI and the CCI simultaneously at the ecosystem scale. In addition, the spatial resolution of the PRI and the CCI obtained from GCOM-C/SGLI is 250 m, which is four times finer than that of the MODIS CCI [35,36,37].
Nevertheless, to the author’s best knowledge, no study has conducted even the accuracy assessments of the PRI and the CCI yet in forests. Therefore, the purpose of this study is to conduct the accuracy assessments of the PRI and the CCI derived from GCOM-C/SGLI by comparing them with in situ observation data of forests.

2. Materials and Methods

2.1. Study Sites

The current study was conducted at four forest sites (Teshio: TSE, Takayama: TKY, Fuji Yoshida: FJY, Fuji Hokuroku: FHK) in Japan between 2018 and 2020 (Figure 1 and Figure 2, Table 1). These four sites belong to AsiaFlux [44], Japanese Long Term Ecological Research Network (JaLTER) [45], and Phenological Eyes Network (PEN) [46].
TSE is located in the northern part of Hokkaido Prefecture. The dominant species of the canopy is the young Japanese larch (a hybrid of Larix kaempferi and L. gmelinii), which is a deciduous needleleaf tree. The dominant species of the understory is the dwarf bamboo (Sasa senanensis or S. kurilensis). TSE is generally covered by snow from November to April. TKY is located on the northwestern slope of Mt. Norikura, which is a part of the Hida Mountain Range. The dominant species of the canopy are a variety of Mongolian oak (Quercus crispula), the Japanese white birch (Betula platyphylla Sukatchev var. japonica Hara), and the Erman’s birch (B. ermanii). All these species are deciduous broadleaf trees. The understory is dominated by the dwarf bamboo (S. senanensis). TKY is usually covered by snow from December to March. FJY is located about 10 k m north of Mt. Fuji. The dominant species of the canopy is the red pine (Pinus densiflora), which is an evergreen needle leaf tree. A variety of Mongolian oak (Q. crispula) and the jolcham oak (Q. serrata) dominate the understory. FJY is occasionally covered by snow in winter. FHK is also located about 10 k m north of Mt. Fuji, and FHK is approximately 1.3   k m away from FJY. The canopy is dominated by the grown Japanese larch (L. kaempferi), and co-dominated by red pine (P. densiflora), Cornus controversa, and a variety of Mongolian oak (Q. crispula). The ferns (Dryopteris crassirhizoma or D. expansa) and the dwarf bamboo (Sasamorpha borealis) dominate the understory. FHK is also occasionally covered by snow in winter.

2.2. In Situ Data

2.2.1. In Situ Data Collection

At each study site, we collected spectral irradiance data and fisheye images. The spectral irradiance data were used for the accuracy assessment of PRI and CCI, and the fisheye images were used for observing the conditions of vegetation, such as snow covers, leafing, autumn colors, and leaf falling. We installed the observation instruments for spectral irradiance data and fisheye images around the top of the each observation tower (Figure 3 and Table 2).
In situ spectral irradiance data were measured with hemispherical spectroradiometers (HSSR): MS-700 (EKO Instrument Co., Ltd., Tokyo, Japan). The step of data acquisition was 3.3   n m from 350 n m to 1050 n m , and the full width at half maximum (FWHM) was 10 n m . MS-700 has been used in many studies [49,50,51,52,53,54]. At each of TSE, FJY, and FHK, two MS-700 were installed around the top of the each observation tower (Figure 3a). One was fixed upward to measure the incident light from the sky, and the other was fixed downward to measure the reflected light from the vegetation (Figure 3a). At TKY, one MS-700 was installed, and the whole MS-700 was rotated by an external motor upward and downward by turns to measure the incident light and the reflected light (Figure 3b). The duration of one observation maneuver at TKY to measure the incident light and the reflected light was ten minutes. The MS-700 used in this study was occasionally (approximately every two years) calibrated by the manufacturer, using a National Institute of Standards and Technology (NIST) traceable halogen lamp.
At TSE, TKY, and FHK, we attached a masking device [55] to MS-700 for excluding the near-horizontal light from the sky and the reflected light from the body part of the observation towers (Figure 3). At FJY, the masking device was not installed.
We also collected fisheye images of the vegetation with time-lapse cameras named Automatic-capturing Digital Fisheye Camera (ADFC) (Figure 3). ADFC is a combined system of a digital camera (COOLPIX4300 or COOLPIX4500, Nikon Corp., Tokyo, Japan), a fisheye lens (FC-E8, Nikon Corp., Tokyo, Japan), and a waterproof housing case. Fisheye images taken by ADFC have been used in many studies [50,52,56,57]. Figure 4 shows examples of the fisheye images.
Acquisition intervals of the spectral irradiance data and fisheye images differed in the four study sites. The intervals are shown in Table 3. The spectral irradiance data and fisheye images were sometimes not collected. At TSE and FHK, spectral irradiance data were not collected in winter because we stopped and removed MS-700, being wary of damages caused by low temperature and snow. At TKY, the reflectance was not calculated in 2020 because of a malfunction of the masking device. At FJY, the fisheye images were not fully collected because of the problem of ADFC. However, FJY is closely located to FHK (see Figure 1), so the snow condition at FJY was referred to the snow condition at FHK in this study.

2.2.2. In Situ Data Processing

From the in situ spectral irradiance data (MS-700 data), we calculated the time series of each vegetation index in two ways, namely, “original” and “simulated”, for each day and each site. The “original” indices followed the original definition of PRI and CCI (Equations (1) and (2), respectively). The “simulated” indices used weighted-average reflectance within each band using GCOM-C/SGLI’s relative spectral response (RSR) (see Figure 5 and Table 4). We describe their details in the followings.
To derive the original PRI ( PRI original ) and the original CCI ( CCI original ), the monochromatic reflectance (the reflectance in each band of MS-700) ρ ( λ ) was firstly calculated as follows:
ρ ( λ ) = g ( λ ) f ( λ )
where f ( λ ) and g ( λ ) indicate the spectral irradiance data of incident light and reflected light at λ   nm , respectively.
The spectral irradiance data precisely at 531 n m , 570 n m , and  645 n m were estimated by liner interpolation, because the spectral resolution of MS-700 was 3.3 nm and the reflectance precisely at these three wavelengths were not available. Thus, the spectral irradiance at these three wavelengths were estimated by linear interpolation. The center wavelengths of the MS-700 bands, which were the first and second closest to three wavelengths and used for linear interpolation, are listed in Table A1. PRI original and CCI original were calculated as
PRI original = ρ ( 531 ) ρ ( 570 ) ρ ( 531 ) + ρ ( 570 )
and
CCI original = ρ ( 531 ) ρ ( 645 ) ρ ( 531 ) + ρ ( 645 ) .
To derive the simulated PRI ( PRI simulated ) and the simulated CCI ( CCI simulated ), we simulated bands’ values of GCOM-C/SGLI as follows. First, we resampled the spectral irradiance measured by MS-700 from 3.3   n m spectral resolution to 0.1   n m spectral resolution with linear interpolation. Then, we calculated the following values:   
VN x simulated = Λ g ( λ ) R S R x ( λ ) d λ Λ f ( λ ) R S R x ( λ ) d λ
where x represents the band number, VN x simulated represents simulated VN x value, Λ indicates the integration section (300 n m – 1100 n m ), R S R x ( λ ) is RSR of the band x, and  d λ = 0.1   n m . The RSR was provided at 0.1   n m spectral resolution on the JAXA’s official webpage [58]. PRI simulated and CCI simulated were calculated as:
PRI simulated = VN 05 simulated VN 06 simulated VN 05 simulated + VN 06 simulated
and
CCI simulated = VN 05 simulated VN 08 simulated VN 05 simulated + VN 08 simulated .

2.3. Satellite Data

2.3.1. Satellite Data Collection

GCOM-C/SGLI level 2 atmospheric corrected land surface reflectance (RSRF) daily products, whose version was 3, were used as the satellite data. The spatial resolution of the bands in visible and NIR regions of the products was 250   m . The details of each band are described in Table 4. Note that some parts of the products are not released from JAXA at present because the major update for the products (from version 2 to version 3) started in November 2021, and it has not been finished yet. Hence, the authors created the RSRF products by applying version 3 algorithm for version 2 input data: the top of atmosphere radiance (LTOA) products. Certainly, version 2 LTOA products were not the latest; however, the difference in LTOA products between version 2 and version 3 was only tiny debugs. Therefore, the RSRF products used in our study can be regarded as the same as the latest RSRF daily products, which will be freely available on JAXA’s FTP server [60] (HDF5 format, WGS84 datum, sinusoidal projection). They included a quality assessment (QA) flag whose details are shown in Table 5.

2.3.2. Satellite Data Processing

We used the highest quality RSRF data: the RSRF data whose QA flag equaled 2 (only the bit 1 equaled 1 and the others equaled 0). The satellite PRI ( PRI satellite ) and the satellite CCI ( CCI satellite ) were calculated using the following equations:
PRI satellite = VN 05 satellite VN 06 satellite VN 05 satellite + VN 06 satellite
and
CCI satellite = VN 05 satellite VN 08 satellite VN 05 satellite + VN 08 satellite ,
where VN x satellite represents the reflectance value of RSRF products of the band x. We calculated PRI satellite and CCI satellite from the RSRF products with a Python package named “h5py” (version 2.10.0). The data processing tools developed by the authors are freely available on GitHub [62]. For the accuracy assessment of PRI satellite and CCI satellite , the pixel whose center was the nearest to each site location was extracted from the RSRF products. To select the nearest pixel, we reprojected the RSRF products from sinusoidal projection to the equirectangular projection based on the GCOM-C/SGLI manual [59]. The location of each study site and the nearest pixel were displayed in Figure 6.

2.4. Comparison and Statistical Analysis

For the sake of comparison and accuracy assessment, we extracted the original and simulated indices measured with MS-700 at the closest time to the observation time of GCOM-C/SGLI for each day and each site from their daily time series data calculated in Section 2.2.2. The observation time of GCOM-C/SGLI was obtained from the RSRF products.
Then, we displayed and compared the seasonal changes in the original, simulated, and satellite indices. The vegetation’s conditions (snow covers, leafing, autumn colors, and leaf falling) were interpreted from the downward fisheye images taken by ADFC, and the conditions were simultaneously displayed. In addition, we examined the agreement between the simulated indices and the satellite indices with the scatter plots and statistics. We calculated the coefficient of correlation (r), root mean square error ( RMSE ), and mean absolute error ( MAE ) as follows:
r = i = 1 n ( x i x ¯ ) ( y i y ¯ ) i = 1 n ( x i x ¯ ) 2 i = 1 n ( y i y ¯ ) 2 ,
RMSE = 1 n i = 1 n ( x i y i ) 2 ,
and
MAE = 1 n i = 1 n | x i y i |
where n represents the sample size; x i and y i indicate the ith sample data of simulated index and satellite index; x ¯ and y ¯ were their sample means, respectively.

3. Results

3.1. Accuracy Assessment of PRI

Focusing on PRI original and PRI simulated , Figure 7 shows similar seasonal trends between them. During the growing seasons (from leafing to leaf falling), however, PRI simulated showed a smaller range of seasonal variation than PRI original . This was because they were approximately the same around the leafing and autumn colors, but PRI simulated was smaller than PRI original at their peak (around DOY = 200 ).
As for PRI satellite , Figure 7 and Figure 8 showed that PRI satellite had some significant errors (hereinafter called “outliers”), even if data screening with the QA flag, which can extract satellite data in the best condition, was applied. The outliers occurred in all four study sites, providing poor agreement.
In addition to the outliers, the relatively tiny errors of PRI satellite , which fluctuated around the PRI simulated (hereinafter called “small noise”), also affected the results of the accuracy assessments. For example, PRI satellite in 2018 at TSE fairly matched with PRI simulated (Figure 7a); however, the scatter plot, including 2018, 2019, and 2020, shows the range of small noise in PRI satellite was approximately the same as the range of seasonal variations in PRI simulated (Figure 8a). Due to the small noise in PRI satellite , the features of seasonal variations in PRI became unclear.
Therefore, owing to the outliers and small noise, the results of the accuracy assessments were not as good as expected. The best agreement between PRI simulated and PRI satellite was 0.294 ( p < 0.05 ) at FHK (Table 6).

3.2. Accuracy Assessment of CCI

Figure 9 shows that CCI simulated had roughly similar seasonal trends to CCI original . However, in detail, there were slight differences between them. The range of seasonal variation in CCI simulated was larger than that of CCI original . At TSE and TKY, they were approximately the same around the leafing and the leaf falling seasons, but around their peak, CCI simulated was larger than CCI original . At FJY, CCI simulated was slightly higher than CCI original except for the beginning of winter. At FHK, they were almost the same from the leafing to the middle of summer, and then CCI simulated was larger than CCI original . Around the autumn colors season, CCI simulated was smaller than CCI original . Taking together, we found that CCI simulated showed slightly larger seasonal variation than CCI original because the wavelength at the red region band of CCI simulated was longer than that of CCI original (see Figure 5).
About the CCI satellite , Figure 9 shows that CCI satellite contained many outliers without the QA flag, but most of them were eliminated with the QA flag. Small noise was also contaminated with CCI satellite , but the seasonal trends of CCI satellite were clear even with small noise. The reason was that the range of seasonal variation of CCI was relatively large and it was not strongly affected by the small noise (Figure 10).
Especially at TKY, there were few outliers and small noise, resulting in the best agreement between CCI simulated and CCI satellite : r = 0.911 ( p < 0.001 ), R M S E = 0.079 , and M A E = 0.058 (Table 7). At both FHK and TSE, which is the deciduous needleleaf forest, CCI satellite reasonably matched with CCI simulated , but CCI satellite showed positive bias at TSE (Figure 10a). At FJY, the coefficient of correlation was relatively low, and the relationships between CCI simulated and CCI satellite were seasonally changed; CCI satellite underestimated CCI simulated around winter but overestimated CCI simulated around summer. Overall, the results of the accuracy assessment of CCI satellite were better than those of PRI satellite .

4. Discussion

The PRI and the CCI are carotenoid-sensitive vegetation indices derived from the narrow green band reflectance [1,2,35,36,37]. The two indices are considered practical tools for monitoring the photosynthetic activities of vegetation [29,30,39,40,41]. According to the specifications of GCOM-C/SGLI, the satellite has the potential to derive the PRI and the CCI. Therefore, this study validated the accuracy of the satellite PRI and CCI derived from GCOM-C/SGLI by comparing them with in situ spectral data at four forest sites. As a result, we found that PRI satellite was poorly matched with PRI simulated (the best: r = 0.294 ( p < 0.05 ) at FHK), and by contrast, CCI satellite matched well with CCI simulated (the best: r = 0.911 ( p < 0.001 ) at TKY).
Compared to CCI satellite , PRI satellite was more strongly affected by outliers and small noise, which contributed to the remarkable differences in the agreements between PRI satellite and CCI satellite . Therefore, in the following section, we will discuss the reasons for the outliers and small noise. Then, we will discuss the differences in the agreements of CCI satellite and CCI simulated for each site, focusing on the footprint of the sensors and forest structures of the study sites because the agreements between CCI satellite and CCI simulated differed for each study site.

4.1. Common Outliers for PRI satellite and CCI satellite

Some outliers were found in both PRI satellite and CCI satellite on the same day. Figure 7 and Figure 9 showed that the common outliers sometimes occurred in snow seasons. These outliers must not be correctly removed by data screening with the QA flag (especially bit 5 for detecting snow or ice (see Table 5)).
To investigate the effects of snow contamination on PRI satellite and CCI satellite , we extracted snow-free PRI satellite and CCI satellite , and created the scatter plots with statistical analysis. We manually removed the PRI satellite and CCI satellite observed during the snow seasons (the gray bands in Figure 7 and Figure 9), and one day before and after the snow seasons by checking the downward fisheye images taken by ADFC. Figure 11 shows the scatter plots of the snow-free PRI and CCI. At FJY, where the in situ data were continuously collected even in winter, the positive outliers of PRI satellite and CCI satellite around the winter seasons were successfully eliminated (Figure 11c,g). In terms of CCI, the seasonal dependency of the relationships between CCI satellite and CCI simulated at FJY became clear. As for the statistical analysis, the correlation coefficient at FJY was improved for both PRI and CCI (Table 8). Nevertheless, at three other sites, the agreements between simulated and satellite indices were not substantially improved. The results suggest that we still need to remove the outliers which occurred in non-snow seasons.
We also found common outliers in non-snow seasons. We illustrated one example of common outliers at FHK with true color images and spatial distributions of PRI satellite and CCI satellite in Figure 12. The figure shows that the observation site was located on the boundary region between the screened and non-screened areas. In order to reveal the actual sky condition at the observation time, we displayed the sky images taken by upward ADFC around the observation time of GCOM-C/SGLI in Figure 13. The figure implies that it must be cloudy at the observation time. This finding suggests that the cloud screening with the QA flag (bit 6 and bit 7 for detecting cloud) might be insufficient, and the overlooking of the cloud might occur around the boundary region.
As Motohka et al. [50] reported, such cloud contamination is a severe problem for the analysis of vegetation indices. Therefore, the QA flag of GCOM-C/SGLI should be improved in future updates of the products.

4.2. Unique Outliers for PRI satellite

In addition to the common outliers for PRI satellite and CCI satellite , we found two kinds of outliers unique to PRI satellite . One is the “striping” outliers and the other is the “cluster” outliers.

4.2.1. Striping Outliers

One example of the striping outliers is displayed in Figure 14. Figure 14 shows that the striping outliers occurred parallel from northeast to southwest only in PRI satellite . We also illustrated the spatial distribution of the reflectance of VN05, VN06, and VN08, which were used for calculating PRI satellite and CCI satellite in Figure 15. Figure 15 shows that no striping outliers occurred in VN05 and VN06. In short, the striping outliers must be due to the position of detectors on SGLI for each band and the nearest neighbor resampling for geometric correction of RSRF products. According to the official document of GCOM-C/SGLI [59], SGLI is equipped with the detectors of each band in a line; each detector has different relationships between the spatial coordinates on the satellite data and spatial coordinates on the Earth’s surface. The detectors of VN05 and VN06 are installed relatively far from each other, whereas the detectors of VN05 and VN08 are installed next to each other. Hence, the difference in the relationships of spatial coordinates between VN05 and VN06 was larger than that between VN05 and VN08. Then, because of the geometric correction with nearest neighbor resampling, the correspondence of pixels between VN05 and VN06 must sometimes be displaced by more than one pixel. As a result, after the geometric correction, some pixels in VN05 and VN06 must refer to the different pixels in the input data of geometric correction. This difference possibly caused the striping outliers in PRI satellite when PRI satellite was calculated. By contrast, the displacement between VN05 and VN08 must be less than one pixel. Hence, pixels in VN05 and VN08 may refer the same pixels in the input data of geometric correction, and no striping outliers was caused in CCI satellite .

4.2.2. Cluster Outliers

In addition to the striping outliers, one example of the cluster outliers is illustrated in Figure 16. We found that the outliers were distributed like clusters only in PRI satellite (Figure 16c), but not in CCI satellite (Figure 16d). We also displayed the spatial distribution of the reflectance of VN05, VN06, and VN08 in Figure 17. According to Figure 17, VN06 had the cluster outliers around where PRI satellite contained the cluster outliers. Thus, the cluster outliers in VN 06 satellite may be one factor of the cluster outliers in PRI satellite . As for the cluster outliers in VN06, however, we found no significant reasons. Thus, we will continue to pursue the reasons and solutions for the cluster outliers in future work.

4.3. Demonstrations of Removing the Outliers

As we discussed in Section 4.1 and Section 4.2, both PRI satellite and CCI satellite contained the common outliers, and PRI satellite uniquely included the striping and cluster outliers. Hence, we need to remove these outliers for the application of PRI satellite and CCI satellite . Here, we propose to calculate the spatial mean as one method to remove the outliers.
According to the JAXA’s official report, the uncertainty of the geometric correction was less than 0.5 pixels [63]. Thus, we calculated the mean and standard deviation of four neighbor pixels, including the nearest neighbor pixel shown in Figure 6. For the calculation of the mean and standard deviation, if at least one pixel in four neighbor pixels was screened with the QA flag, we did not calculate the mean and standard deviation for accuracy assessments.
First, we show the results of PRI in Figure 18a–d and Table 9. As a result of calculating the spatial mean, many outliers were removed, and r, RMSE , and MAE were improved in many cases. At TKY, the results were drastically improved, and at FJY, the outliers, which we could not remove by using snow-free data set in Section 4.4, were eliminated.
Next, the results of CCI are shown in Figure 18e–h and Table 9. The outliers of CCI were also removed, and especially at FJY and FHK, we successfully removed the significant outliers. By contrast, the statistical result at TKY was not improved. The reason may be the spatial resolution of CCI satellite became too coarse, by calculating the mean value with four neighbor pixels.

4.4. Small Noise and Index Design

In addition to the outliers, PRI satellite and CCI satellite included small noise. One of the factors of small noise can be the uncertainty of the geometric correction. JAXA reported that the uncertainty of the geometric correction in version 3 RSRF products was of less than 0.5 pixels [63]. However, uncertainty of less than 0.5 pixels can be insufficient, particularly at TSE, where the homogeneous vegetation area was limited (see Figure 6a–d). The uncertainty possibly caused the contamination of other vegetation types, which must induce the small noise.
Moreover, the uncertainty of atmospheric correction can also be responsible for the small noise. We demonstrated the accuracy assessment of VN05, VN06, and VN08 in Figure A1 and Table A2. The results showed insufficient accuracy, especially for VN05 and VN06.
The small noise affected the results of the accuracy assessment of PRI satellite more strongly than CCI satellite , as seen in Figure 8 and Figure 10. The difference in the effects of small noise between PRI satellite and CCI satellite may be attributed to the design of each index. As can be seen in Figure 7 and Figure 9, the range of seasonal variation in PRI satellite was smaller than that in CCI satellite . Hence, PRI satellite was easily and strongly affected by the small noise in comparison to CCI satellite . At the leaf and the canopy scale, such a small noise might not occur frequently and PRI worked well to monitor vegetation’s conditions. However, small noise easily happened for satellite remote sensing, which requires atmospheric correction, geometric correction, and bidirectional reflectance distribution function (BRDF) correction. Therefore, compared to PRI, CCI may be more suitable for monitoring the vegetation with satellite remote sensing because CCI is “resistant” and PRI is “susceptible” to the small noise.

4.5. Footprint Effects for Accuracy Assessment of CCI satellite

As listed in Table 6 and Table 7, the accuracy of CCI satellite was higher than that of PRI satellite at all study sites. However, the agreements between CCI satellite and CCI simulated differed for each site. Especially at TSE and FJY, we found the distinctive relationships between CCI satellite and CCI simulated (Figure 10):
At TSE, indeed, the agreement was reasonably good ( r = 0.868 ( p < 0.001 )), but CCI satellite showed a positive bias. The positive bias was probably due to the differences in the footprint between in situ sensor (MS-700) and the satellite sensor (SGLI). As mentioned in Section 2.1, the dominant species at TSE was “young” Japanese larch, and its population density at the canopy was not high, as shown in Figure 19a. Hence, around the edge of the observation area of MS-700, MS-700 laterally observe the sides of trees (Japanese larch), rather than from straight above, and the trees tend to hide the understory (dwarf bamboo) behind them. In contrast, GCOM-C/SGLI observes the trees from above. Therefore, the contribution of Japanese larch to the reflectance observed by MS-700 might be larger than that observed by the GCOM-C/SGLI, and it may cause the positive bias.
At FJY, the relationships between CCI satellite and CCI simulated were seasonally changed, as seen in Figure 10c (clearer in Figure 11g). The seasonal dependency might be affected by the difference in footprints between MS700 and SGLI. As we mentioned in Section 2.2.1, the masking device, which excludes the reflected light from the body part of the observation tower, was not installed at FJY. Thus, as seen in Figure 19c, around ten percent of the footprint of MS-700 was occupied by the observation tower, whereas the area of the observation tower in one pixel of SGLI was tiny. Then, the observation tower contributed the reflectance measured by MS-700 more significantly than that observed by SGLI. The reflected light from the observation tower can be affected by some variables, such as Sun elevation, which changes seasonally. Therefore, the relationships between CCI simulated and CCI satellite may be seasonally changed.
As discussed above, the difference in the footprint might be responsible for the distinctive relationships between CCI simulated and CCI satellite at TSE and FJY. We are aware that the difference in the footprint was the limitation of this study. In contrast, at TKY and FHK, the forests were relatively close to climax (Figure 19b,d), and the vegetation’s condition in the footprint of MS-700 may be similar to that of GCOM-C/SGLI. Indeed, the results of accuracy assessment at TKY and FHK were better than those at TSE and FJY.

5. Conclusions

The accuracy assessments of PRI and CCI derived from GCOM-C/SGLI were conducted by comparing them with in situ data. As a result, GCOM-C/SGLI provided PRI satellite with poor accuracy (the best: r = 0.294 ( p < 0.05 ) at FHK) and CCI satellite with good accuracy (the best: r = 0.911 ( p < 0.001 ) at TKY). Thus, CCI satellite must be suitable for monitoring photosynthetic activities of vegetation with RSRF products (version 3) of GCOM-C/SGLI rather than PRI satellite .
We found that some outliers, possibly caused by insufficient QA flags, affected the results of the accuracy assessments of both PRI satellite and CCI satellite . Moreover, there were two kinds of outliers unique to PRI satellite . The first was the striping outliers, which might be caused by the combination of the position of detectors on SGLI and the nearest neighbor sampling method for geometric correction. The second was the cluster outliers, which might be caused by the cluster outliers in VN06. The two kinds of outliers may account for the lower accuracy of PRI satellite in comparison to CCI satellite .
In addition to the outliers, small noise also affected the accuracy of PRI satellite more significantly than CCI satellite , because of the smaller range of seasonal variation in PRI satellite . For vegetation monitoring with satellite remote sensing techniques, small noise may be unavoidable, so we should consider the “susceptibility” of the vegetation index to small noise as well as the ability to represent the target phenomena.
In the future, we should develop a method to remove or correct the outliers and continue to validate the accuracy of PRI satellite and CCI satellite in other vegetation types: paddy field, cropland, and grassland, for instance, as partially investigated by Bayarsaikhan et al. [64]. Furthermore, we need to assess the accuracy over enough footprints; at least one pixel of RSRF products of GCOM-C/SGLI. For the application of PRI satellite and CCI satellite , we plan to validate the their ability to track photosynthetic activities by comparing them with in situ eco-physiological data, such as LUE and GPP. We also pay attention to the availability of the green-red vegetation index (GRVI) [56] derived from other satellites, such as MODIS, Sentinel-2, and Himawari-8. The GRVI does not require “narrow” green band reflectance, and Yin et al. [65] have reported that GRVI performed similarly to CCI. Such combinations with eco-physiological data and collaborations over the satellites must provide a better understanding of terrestrial ecosystems with PRI and CCI.

Author Contributions

Conceptualization, T.S., T.K.A. and K.N.N.; methodology, T.S. and T.K.A.; software, T.S.; formal analysis, T.S.; field work and data curation, T.S., T.K.A., R.I., K.T., S.T., T.N. and K.N.N.; writing—original draft preparation, T.S.; writing—review and editing, T.S., T.K.A., R.I., S.T. and K.N.N.; visualization, T.S.; supervision, T.S., T.K.A. and K.N.N.; project administration, K.N.N.; funding acquisition, T.S., T.K.A., T.N. and K.N.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by Japan Science and Technology Agency (JST) Support for Pioneering Research Initiated by the Next Generation (SPRING) Grant Numbers JPMJSP2124; Japan Aerospace Exploration Agency (JAXA) Global Change Observation Mission (GCOM) Grant Numbers ER2GCF103 and ER3GCF104; Japan Society for the Promotion of Science (JSPS) KAKENHI Grant Numbers JP19K20433 and JP22H05739.

Data Availability Statement

The data used in this study are publicly available from the following websites with a user account or a user request: GCOM-C/SGLI data: https://gportal.jaxa.jp/gpr/?lang=en (accessed on 23 August 2022) (see Section 2.3.1) and in situ observation data: http://www.pheno-eye.org/ (accessed on 23 August 2022). The satellite data processing tools are available via https://github.com/tigersasagawa/sgli-tools (accessed on 23 August 2022).

Acknowledgments

The authors thank the Japan Aerospace Exploration Agency (JAXA) for providing the satellite data, the Phenological Eyes Network (PEN) for providing in situ data, and Murakami, H. (JAXA); Kobayashi, T. (University of Tsukuba: UT); Phan, C.D. (UT); Iwao, K. (National Institute of Advanced Industrial Science and Technology: AIST); Tsuchida, S. (AIST); and Mizuochi, H. (AIST) for their advice on this study.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. The center wavelengths of the MS-700 bands which were the nearest or the second nearest to 531 nm, 570 nm, and 645 nm and used for the liner interpolation.
Table A1. The center wavelengths of the MS-700 bands which were the nearest or the second nearest to 531 nm, 570 nm, and 645 nm and used for the liner interpolation.
Site ID
and
Period
Direction
of
MS-700
The
Nearest
Neighbor
Band’s Peak
to 531 nm
The Second
Nearest
Neighbor
Band’s Peak
to 531 nm
The
Nearest
Neighbor
Band’s Peak
to 570 nm
The Second
Nearest
Neighbor
Band’s Peak
to 570 nm
The
Nearest
Neighbor
Band’s Peak
to 645 nm
The Second
Nearest
Neighbor
Band’s Peak
to 645 nm
TSE
2018-01-01

2020-12-31
 
upward
 
529.70 nm
 
533.05 nm
 
569.93 nm
 
573.28 nm
 
643.51 nm
 
646.84 nm
 
downward
 
531.35 nm
 
528.00 nm
 
571.52 nm
 
568.18 nm
 
645.12 nm
 
641.78 nm
TKY
2018-01-01

2019-05-07
 
upward
 
529.64 nm
 
532.94 nm
 
569.14 nm
 
572.42 nm
 
644.65 nm
 
647.93 nm
 
downward
 
529.64 nm
 
532.94 nm
 
569.14 nm
 
572.42 nm
 
644.65 nm
 
647.93 nm
TKY
2019-05-08

2020-12-31
 
upward
 
529.88 nm
 
533.17 nm
 
569.33 nm
 
572.62 nm
 
644.93 nm
 
648.22 nm
 
downward
 
529.88 nm
 
533.17 nm
 
569.33 nm
 
572.62 nm
 
644.93 nm
 
648.22 nm
FJY
2018-01-01

2020-12-31
 
upward
 
529.85 nm
 
533.15 nm
 
569.40 nm
 
572.69 nm
 
645.04 nm
 
641.76 nm
 
downward
 
529.76 nm
 
533.06 nm
 
569.34 nm
 
572.64 nm
 
644.99 nm
 
648.28 nm
FHK
2018-01-01

2018-12-31
 
upward
 
532.47 nm
 
529.18 nm
 
568.65 nm
 
571.94 nm
 
644.19 nm
 
647.47 nm
 
downward
 
530.03 nm
 
533.33 nm
 
569.54 nm
 
572.83 nm
 
645.08 nm
 
641.81 nm
FHK
2019-01-01

2020-12-31
 
upward
 
531.79 nm
 
528.45 nm
 
568.58 nm
 
571.92 nm
 
645.43 nm
 
642.09 nm
 
downward
 
531.86 nm
 
528.51 nm
 
568.69 nm
 
572.04 nm
 
645.69 nm
 
642.34 nm
Figure A1. The scatter plots of simulated band value and satellite band value used for calculation of PRI and CCI. The shape of each point represents the year. The color of each point corresponds to the DOY. Note that the range of the x-axis and y-axis is limited from 0 to 0.1, so a few points out of the range are not displayed.
Figure A1. The scatter plots of simulated band value and satellite band value used for calculation of PRI and CCI. The shape of each point represents the year. The color of each point corresponds to the DOY. Note that the range of the x-axis and y-axis is limited from 0 to 0.1, so a few points out of the range are not displayed.
Remotesensing 14 05352 g0a1
Table A2. The statistics of the accuracy assessment of VN05, VN06, and VN08.
Table A2. The statistics of the accuracy assessment of VN05, VN06, and VN08.
Site IDnrRMSEMAE
VN05TSE84 0.856 ( p < 0.001 )0.0160.010
TKY40 0.839 ( p < 0.001 )0.0320.026
FJY146 0.136 ( p = 0.101 )4.9484.908
FHK65 0.343 ( p < 0.01 )0.0120.010
VN05TSE84 0.831 ( p < 0.001 )0.0170.012
TKY40 0.794 ( p < 0.001 )0.0370.030
FJY146 0.199 ( p < 0.05 )5.5945.550
FHK65 0.250 ( p < 0.05 )0.0150.013
VN05TSE84 0.882 ( p < 0.001 )0.0190.014
TKY40 0.899 ( p < 0.001 )0.0310.024
FJY146 0.481 ( p < 0.001 )4.6974.621
FHK65 0.738 ( p < 0.001 )0.0170.013

References

  1. Gamon, J.A.; Field, C.B.; Bilger, W.; Björkman, O.; Fredeen, A.L.; Peñuelas, J. Remote sensing of the xanthophyll cycle and chlorophyll fluorescence in sunflower leaves and canopies. Oecologia 1990, 85, 1–7. [Google Scholar] [CrossRef] [PubMed]
  2. Gamon, J.A.; Peñuelas, J.; Field, C.B. A narrow-waveband spectral index that tracks diurnal changes in photosynthetic efficiency. Remote Sens. Environ. 1992, 41, 35–44. [Google Scholar] [CrossRef]
  3. Gamon, J.A.; Serrano, L.; Surfus, J.S. The photochemical reflectance index: An optical indicator of photosynthetic radiation use efficiency across species, functional types, and nutrient levels. Oecologia 1997, 112, 492–501. [Google Scholar] [CrossRef] [PubMed]
  4. Thenot, F.; Méthy, M.; Winkel, T. The Photochemical Reflectance Index (PRI) as a water-stress index. Int. J. Remote Sens. 2002, 23, 5135–5139. [Google Scholar] [CrossRef]
  5. Suárez, L.; Zarco-Tejada, P.J.; Sepulcre-Cantó, G.; Pérez-Priego, O.; Miller, J.R.; Jiménez-Muñoz, J.C.; Sobrino, J. Assessing canopy PRI for water stress detection with diurnal airborne imagery. Remote Sens. Environ. 2008, 112, 560–575. [Google Scholar] [CrossRef]
  6. Porcar-Castell, A.; Garcia-Plazaola, J.I.; Nichol, C.J.; Kolari, P.; Olascoaga, B.; Kuusinen, N.; Fernández-Marín, B.; Pulkkinen, M.; Juurola, E.; Nikinmaa, E. Physiology of the seasonal relationship between the photochemical reflectance index and photosynthetic light use efficiency. Oecologia 2012, 170, 313–323. [Google Scholar] [CrossRef] [PubMed]
  7. Hmimina, G.; Dufrêne, E.; Soudani, K. Relationship between photochemical reflectance index and leaf ecophysiological and biochemical parameters under two different water statuses: Towards a rapid and efficient correction method using real-time measurements. Plant Cell Environ. 2014, 37, 473–487. [Google Scholar] [CrossRef]
  8. Springer, K.R.; Wang, R.; Gamon, J.A. Parallel Seasonal Patterns of Photosynthesis, Fluorescence, and Reflectance Indices in Boreal Trees. Remote Sens. 2017, 9, 691. [Google Scholar] [CrossRef] [Green Version]
  9. Xu, S.; Liu, Z.; Zhao, L.; Zhao, H.; Ren, S. Diurnal Response of Sun-Induced Fluorescence and PRI to Water Stress in Maize Using a Near-Surface Remote Sensing Platform. Remote Sens. 2018, 10, 1510. [Google Scholar] [CrossRef] [Green Version]
  10. Eitel, J.U.H.; Maguire, A.J.; Boelman, N.; Vierling, L.A.; Griffin, K.L.; Jensen, J.; Magney, T.S.; Mahoney, P.J.; Meddens, A.J.H.; Silva, C.; et al. Proximal remote sensing of tree physiology at northern treeline: Do late-season changes in the photochemical reflectance index (PRI) respond to climate or photoperiod? Remote Sens. Environ. 2019, 221, 340–350. [Google Scholar] [CrossRef]
  11. Yang, J.C.; Magney, T.S.; Yan, D.; Knowles, J.F.; Smith, W.K.; Scott, R.L.; Barron-Gafford, G.A. The photochemical reflectance index (PRI) captures the ecohydrologic sensitivity of a semiarid mixed conifer forest. J. Geophys. Res. Biogeosci. 2020, 125, e2019JG005624. [Google Scholar] [CrossRef]
  12. Pierrat, Z.; Nehemy, M.F.; Roy, A.; Magney, T.; Parazoo, N.C.; Laroque, C.; Pappas, C.; Sonnentag, O.; Grossmann, K.; Bowling, D.R.; et al. Tower-based remote sensing reveals mechanisms behind a two-phased spring transition in a mixed-species boreal forest. J. Geophys. Res. Biogeosci. 2021, 126, e2020JG006191. [Google Scholar] [CrossRef]
  13. Kohzuma, K.; Tamaki, M.; Hikosaka, K. Corrected photochemical reflectance index (PRI) is an effective tool for detecting environmental stresses in agricultural crops under light conditions. J. Plant Res. 2021, 134, 683–694. [Google Scholar] [CrossRef]
  14. Tsujimoto, K.; Hikosaka, K. Estimating leaf photosynthesis of C3 plants grown under different environments from pigment index, photochemical reflectance index, and chlorophyll fluorescence. Photosynth. Res. 2021, 148, 33–46. [Google Scholar] [CrossRef] [PubMed]
  15. Filella, I.; Porcar-Castell, A.; Munné-Bosch, S.; Bäck, J.; Garbulsky, M.F.; Peñuelas, J. PRI assessment of long-term changes in carotenoids/chlorophyll ratio and short-term changes in de-epoxidation state of the xanthophyll cycle. Int. J. Remote Sens. 2009, 30, 4443–4455. [Google Scholar] [CrossRef]
  16. Gamon, J.A.; Kovalchuck, O.; Wong, C.Y.S.; Harris, A.; Garrity, S.R. Monitoring seasonal and diurnal changes in photosynthetic pigments with automated PRI and NDVI sensors. Biogeosciences 2015, 12, 4149–4159. [Google Scholar] [CrossRef] [Green Version]
  17. Gitelson, A.A.; Gamon, J.A.; Solovchenko, A. Multiple drivers of seasonal change in PRI: Implications for photosynthesis 1. Leaf level. Remote Sens. Environ. 2017, 191, 110–116. [Google Scholar] [CrossRef] [Green Version]
  18. Gitelson, A.A.; Gamon, J.A.; Solovchenko, A. Multiple drivers of seasonal change in PRI: Implications for photosynthesis 2. Stand level. Remote Sens. Environ. 2017, 190, 198–206. [Google Scholar] [CrossRef] [Green Version]
  19. Monteith, J.L. Solar Radiation and Productivity in Tropical Ecosystems. J. Appl. Ecol. 1972, 9, 747–766. [Google Scholar] [CrossRef] [Green Version]
  20. Monteith, J.L.; Moss, C.J.; Cooke, G.W.; Pirie, N.W.; Bell, G.D.H. Climate and the efficiency of crop production in Britain. Philos. Trans. R. Soc. Lond. B Biol. Sci. 1977, 281, 277–294. [Google Scholar] [CrossRef]
  21. Penuelas, J.; Filella, I.; Gamon, J.A. Assessment of photosynthetic radiation-use efficiency with spectral reflectance. New Phytol. 1995, 131, 291–296. [Google Scholar] [CrossRef]
  22. Peñuelas, J.; Llusia, J.; Pinol, J.; Filella, I. Photochemical reflectance index and leaf photosynthetic radiation-use-efficiency assessment in Mediterranean trees. Int. J. Remote Sens. 1997, 18, 2863–2868. [Google Scholar] [CrossRef]
  23. Winkel, T.; Méthy, M.; Thénot, F. Radiation use efficiency, chlorophyll fluorescence, and reflectance indices associated with ontogenic changes in water-limited Chenopodium quinoa leaves. Photosynthetica 2002, 40, 227–232. [Google Scholar] [CrossRef]
  24. Nakaji, T.; Oguma, H.; Fujinuma, Y. Seasonal changes in the relationship between photochemical reflectance index and photosynthetic light use efficiency of Japanese larch needles. Int. J. Remote Sens. 2006, 27, 493–509. [Google Scholar] [CrossRef]
  25. Nichol, C.J.; Huemmrich, K.F.; Andrew Black, T.; Jarvis, P.G.; Walthall, C.L.; Grace, J.; Hall, F.G. Remote sensing of photosynthetic-light-use efficiency of boreal forest. Agric. For. Meteorol. 2000, 101, 131–142. [Google Scholar] [CrossRef] [Green Version]
  26. Nichol, C.J.; Lloyd, J.; Shibistova, O.; Arneth, A.; Röser, C.; Knohl, A.; Matsubara, S.; Grace, J. Remote sensing of photosynthetic-light-use efficiency of a Siberian boreal forest. Tellus B Chem. Phys. Meteorol. 2002, 54, 677–687. [Google Scholar] [CrossRef]
  27. Nakaji, T.; Ide, R.; Oguma, H.; Saigusa, N.; Fujinuma, Y. Utility of spectral vegetation index for estimation of gross CO2 flux under varied sky conditions. Remote Sens. Environ. 2007, 109, 274–284. [Google Scholar] [CrossRef]
  28. Hall, F.G.; Hilker, T.; Coops, N.C.; Lyapustin, A.; Huemmrich, K.F.; Middleton, E.; Margolis, H.; Drolet, G.; Black, T.A. Multi-angle remote sensing of forest light use efficiency by observing PRI variation with canopy shadow fraction. Remote Sens. Environ. 2008, 112, 3201–3211. [Google Scholar] [CrossRef] [Green Version]
  29. Garbulsky, M.F.; Peñuelas, J.; Gamon, J.; Inoue, Y.; Filella, I. The photochemical reflectance index (PRI) and the remote sensing of leaf, canopy and ecosystem radiation use efficiencies: A review and meta-analysis. Remote Sens. Environ. 2011, 115, 281–297. [Google Scholar] [CrossRef]
  30. Zhang, C.; Filella, I.; Garbulsky, M.; Peñuelas, J. Affecting Factors and Recent Improvements of the Photochemical Reflectance Index (PRI) for Remotely Sensing Foliar, Canopy and Ecosystemic Radiation-Use Efficiencies. Remote Sens. 2016, 8, 677. [Google Scholar] [CrossRef]
  31. Huang, W.; Lamb, D.W.; Niu, Z.; Zhang, Y.; Liu, L.; Wang, J. Identification of yellow rust in wheat using in-situ spectral reflectance measurements and airborne hyperspectral imaging. Precis. Agric. 2007, 8, 187–197. [Google Scholar] [CrossRef]
  32. Stagakis, S.; Markos, N.; Sykioti, O.; Kyparissis, A. Monitoring canopy biophysical and biochemical parameters in ecosystem scale using satellite hyperspectral imagery: An application on a Phlomis fruticosa Mediterranean ecosystem using multiangular CHRIS/PROBA observations. Remote Sens. Environ. 2010, 114, 977–994. [Google Scholar] [CrossRef]
  33. Rossini, M.; Fava, F.; Cogliati, S.; Meroni, M.; Marchesi, A.; Panigada, C.; Giardino, C.; Busetto, L.; Migliavacca, M.; Amaducci, S.; et al. Assessing canopy PRI from airborne imagery to map water stress in maize. ISPRS J. Photogramm. Remote Sens. 2013, 86, 168–177. [Google Scholar] [CrossRef]
  34. Stagakis, S.; Markos, N.; Sykioti, O.; Kyparissis, A. Tracking seasonal changes of leaf and canopy light use efficiency in a Phlomis fruticosa Mediterranean ecosystem using field measurements and multi-angular satellite hyperspectral imagery. ISPRS J. Photogramm. Remote Sens. 2014, 97, 138–151. [Google Scholar] [CrossRef]
  35. Gamon, J.A.; Fred Huemmrich, K.; Wong, C.Y.S.; Ensminger, I.; Garrity, S.; Hollinger, D.Y.; Noormets, A.; Peñuelas, J. A remotely sensed pigment index reveals photosynthetic phenology in evergreen conifers. Proc. Natl. Acad. Sci. USA 2016, 113, 13087–13092. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Drolet, G.G.; Huemmrich, K.F.; Hall, F.G.; Middleton, E.M.; Black, T.A.; Barr, A.G.; Margolis, H.A. A MODIS-derived photochemical reflectance index to detect inter-annual variations in the photosynthetic light-use efficiency of a boreal deciduous forest. Remote Sens. Environ. 2005, 98, 212–224. [Google Scholar] [CrossRef]
  37. Middleton, E.M.; Huemmrich, K.F.; Landis, D.R.; Black, T.A.; Barr, A.G.; McCaughey, J.H. Photosynthetic efficiency of northern forest ecosystems using a MODIS-derived Photochemical Reflectance Index (PRI). Remote Sens. Environ. 2016, 187, 345–366. [Google Scholar] [CrossRef]
  38. Kim, J.; Ryu, Y.; Dechant, B.; Lee, H.; Kim, H.S.; Kornfeld, A.; Berry, J.A. Solar-induced chlorophyll fluorescence is non-linearly related to canopy photosynthesis in a temperate evergreen needleleaf forest during the fall transition. Remote Sens. Environ. 2021, 258, 112362. [Google Scholar] [CrossRef]
  39. Wang, R.; Gamon, J.A.; Emmerton, C.A.; Springer, K.R.; Yu, R.; Hmimina, G. Detecting intra- and inter-annual variability in gross primary productivity of a North American grassland using MODIS MAIAC data. Agric. For. Meteorol. 2020, 281, 107859. [Google Scholar] [CrossRef]
  40. Wong, C.Y.S.; D’Odorico, P.; Bhathena, Y.; Arain, M.A.; Ensminger, I. Carotenoid based vegetation indices for accurate monitoring of the phenology of photosynthesis at the leaf-scale in deciduous and evergreen trees. Remote Sens. Environ. 2019, 233, 111407. [Google Scholar] [CrossRef]
  41. Wong, C.Y.S.; D’Odorico, P.; Arain, M.A.; Ensminger, I. Tracking the phenology of photosynthesis using carotenoid-sensitive and near-infrared reflectance vegetation indices in a temperate evergreen and mixed deciduous forest. New Phytol. 2020, 226, 1682–1695. [Google Scholar] [CrossRef] [PubMed]
  42. Imaoka, K.; Kachi, M.; Fujii, H.; Murakami, H.; Hori, M.; Ono, A.; Igarashi, T.; Nakagawa, K.; Oki, T.; Honda, Y.; et al. Global Change Observation Mission (GCOM) for monitoring carbon, water cycles, and climate change. Proc. IEEE 2010, 98, 717–734. [Google Scholar] [CrossRef]
  43. Hori, M.; Murakami, H.; Miyazaki, R.; Honda, Y.; Nasahara, K.; Kajiwara, K.; Nakajima, T.Y.; Irie, H.; Toratani, M.; Hirawake, T.; et al. GCOM-C Data Validation Plan for Land, Atmosphere, Ocean, and Cryosphere. Trans. Jpn. Soc. Aeronaut. Space Sci. Aerosp. Technol. Jpn. 2018, 16, 218–223. [Google Scholar] [CrossRef] [Green Version]
  44. AsiaFlux. Available online: http://www.asiaflux.net/ (accessed on 23 August 2022).
  45. Japan Long Term Ecological Research Network (JaLTER). Available online: http://www.jalter.org/en/ (accessed on 23 August 2022).
  46. Phenological Eyes Network (PEN). Available online: http://www.pheno-eye.org/ (accessed on 23 August 2022).
  47. Kottek, M.; Grieser, J.; Beck, C.; Rudolf, B.; Rubel, F. World Map of the Köppen-Geiger climate classification updated. Meteorol. Z. 2006, 15, 259–263. [Google Scholar] [CrossRef]
  48. Rubel, F.; Brugger, K.; Haslinger, K.; Auer, I. The climate of the European Alps: Shift of very high resolution Köppen-Geiger climate zones 1800–2100. Meteorol. Z. 2017, 26, 115–125. [Google Scholar] [CrossRef]
  49. Nagai, S.; Nasahara, K.N.; Tsuchida, S.; Motohka, T.; Muraoka, H. Phenological eyes network (PEN) and ground-truthing activity for satellite remote sensing. In Proceedings of the 34th International Symposium on Remote Sensing of Environment, Sydney, Australia, 10–15 April 2011. [Google Scholar]
  50. Motohka, T.; Nasahara, K.N.; Murakami, K.; Nagai, S. Evaluation of sub-pixel cloud noises on MODIS daily spectral indices based on in situ measurements. Remote Sens. 2011, 3, 1644–1662. [Google Scholar] [CrossRef] [Green Version]
  51. Choi, J.P.; Kang, S.K.; Choi, G.Y.; Nasahara, K.N.; Motohka, T.; Lim, J.H. Monitoring canopy phenology in a deciduous broadleaf forest using the Phenological Eyes Network (PEN). J. Ecol. Environ. 2011, 34, 149–156. [Google Scholar] [CrossRef] [Green Version]
  52. Nasahara, K.N.; Nagai, S. Review: Development of an in situ observation network for terrestrial ecological remote sensing: The Phenological Eyes Network (PEN). Ecol. Res. 2015, 30, 211–223. [Google Scholar] [CrossRef]
  53. Nagai, S.; Nasahara, K.N.; Inoue, T.; Saitoh, T.M.; Suzuki, R. Review: Advances in in situ and satellite phenological observations in Japan. Int. J. Biometeorol. 2016, 60, 615–627. [Google Scholar] [CrossRef] [Green Version]
  54. Yan, D.; Zhang, X.; Nagai, S.; Yu, Y.; Akitsu, T.; Nasahara, K.N.; Ide, R.; Maeda, T. Evaluating land surface phenology from the Advanced Himawari Imager using observations from MODIS and the Phenological Eyes Network. Int. J. Appl. Earth Obs. Geoinf. 2019, 79, 71–83. [Google Scholar] [CrossRef]
  55. Ide, R.; Hirose, Y.; Oguma, H.; Saigusa, N. Development of a masking device to exclude contaminated reflection during tower-based measurements of spectral reflectance from a vegetation canopy. Agric. For. Meteorol. 2016, 223, 141–150. [Google Scholar] [CrossRef]
  56. Motohka, T.; Nasahara, K.N.; Oguma, H.; Tsuchida, S. Applicability of green-red vegetation index for remote sensing of vegetation phenology. Remote Sens. 2010, 2, 2369–2387. [Google Scholar] [CrossRef] [Green Version]
  57. Nagai, S.; Akitsu, T.; Saitoh, T.M.; Busey, R.C.; Fukuzawa, K.; Honda, Y.; Ichie, T.; Ide, R.; Ikawa, H.; Iwasaki, A.; et al. 8 million phenological and sky images from 29 ecosystems from the Arctic to the tropics: The Phenological Eyes Network. Ecol. Res. 2018, 33, 1091–1092. [Google Scholar] [CrossRef]
  58. SGLI Sensor Characterization. Available online: https://suzaku.eorc.jaxa.jp/GCOM_C/data/prelaunch/index.html (accessed on 23 August 2022).
  59. JAXA. GCOM-C “SHIKISAI” Data Users Handbook. 2018. Available online: https://gportal.jaxa.jp/gpr/assets/mng_upload/GCOM-C/GCOM-C_SHIKISAI_Data_Users_Handbook_en.pdf (accessed on 23 August 2022).
  60. JAXA G-Portal. Available online: https://gportal.jaxa.jp/gpr/?lang=en (accessed on 23 August 2022).
  61. Murakami, H. GCOM-C/SGLI Land Atmospheric Correction Algorithm. 2021. Available online: https://suzaku.eorc.jaxa.jp/GCOM_C/data/ATBD/ver3/V3ATBD_T1A_Atmcorr_murakami.pdf (accessed on 23 August 2022).
  62. Sasagawa, T. GCOM-C/SGLI Data Processing Tools. Available online: https://github.com/tigersasagawa/sgli-tools (accessed on 24 September 2022).
  63. JAXA. Summary of the SGLI Products—Validation Results (Ver. 3.00). 2021. Available online: https://suzaku.eorc.jaxa.jp/GCOM_C/data/files/V3_summary_en.pdf (accessed on 22 August 2022).
  64. Bayarsaikhan, U.; Akitsu, T.K.; Tachiiri, K.; Sasagawa, T.; Nakano, T.; Uudus, B.; Nasahara, K.N. Early validation study of the photochemical reflectance index (PRI) and the normalized difference vegetation index (NDVI) derived from the GCOM-C satellite in Mongolian grasslands. Int. J. Remote Sens. 2022, 43, 5145–5172. [Google Scholar] [CrossRef]
  65. Yin, G.; Verger, A.; Descals, A.; Filella, I.; Peñuelas, J. A broadband green-red vegetation index for monitoring gross primary production phenology. J. Remote Sens. 2022, 2022, 9764982. [Google Scholar] [CrossRef]
Figure 1. Locations of the four study sites where in situ data were collected.
Figure 1. Locations of the four study sites where in situ data were collected.
Remotesensing 14 05352 g001
Figure 2. Overviews of the four study sites. The silvery artificial structure in each photo is the observation tower. The date next to the site name indicates when we took each photo.
Figure 2. Overviews of the four study sites. The silvery artificial structure in each photo is the observation tower. The date next to the site name indicates when we took each photo.
Remotesensing 14 05352 g002
Figure 3. Examples of the instruments: MS-700, masking device, and Automatic-capturing Digital Fisheye Camera (ADFC) at (a) FHK and (b) TKY. At TSE and FJY, the instruments were installed basically in the same manner as (a) FHK, but FJY was not equipped with the masking device for MS-700. At TKY, an external motor rotates MS-700 to observe the incident and reflected light (b).
Figure 3. Examples of the instruments: MS-700, masking device, and Automatic-capturing Digital Fisheye Camera (ADFC) at (a) FHK and (b) TKY. At TSE and FJY, the instruments were installed basically in the same manner as (a) FHK, but FJY was not equipped with the masking device for MS-700. At TKY, an external motor rotates MS-700 to observe the incident and reflected light (b).
Remotesensing 14 05352 g003
Figure 4. Examples of fisheye images taken by ADFC. These images were taken in 2019 at TKY. The bottom numbers represent day of year (DOY).
Figure 4. Examples of fisheye images taken by ADFC. These images were taken in 2019 at TKY. The bottom numbers represent day of year (DOY).
Remotesensing 14 05352 g004
Figure 5. The relative spectral response (RSR) in the visible and near-infrared (NIR) range of Global Change Observation Mission-Climate (GCOM-C)/Second Generation Global Imager (SGLI). The original data was obtained from [58]. The solid black lines represent the wavelength at 531 n m , 570 n m , and 645 n m , originally used to derive the photochemical reflectance index (PRI) and the chlorophyll carotenoid index (CCI). The blue, green, and orange dotted lines are reflectance measured by MS-700 at FHK on 9 April 2020 (DOY = 100 ), 18 July 2020 (DOY = 200 ), and 26 October 2020 (DOY = 300 ), respectively. Each reflectance was observed at 10:31:00 (Japan standard time (JST)) on each day.
Figure 5. The relative spectral response (RSR) in the visible and near-infrared (NIR) range of Global Change Observation Mission-Climate (GCOM-C)/Second Generation Global Imager (SGLI). The original data was obtained from [58]. The solid black lines represent the wavelength at 531 n m , 570 n m , and 645 n m , originally used to derive the photochemical reflectance index (PRI) and the chlorophyll carotenoid index (CCI). The blue, green, and orange dotted lines are reflectance measured by MS-700 at FHK on 9 April 2020 (DOY = 100 ), 18 July 2020 (DOY = 200 ), and 26 October 2020 (DOY = 300 ), respectively. Each reflectance was observed at 10:31:00 (Japan standard time (JST)) on each day.
Remotesensing 14 05352 g005
Figure 6. The location of each study site and the nearest pixel of GCOM-C/SGLI. The background true color image was created from Sentinel-2 level 2 products. The cyan circle indicates the location of the observation tower. The red square indicates a range of the nearest pixel of GCOM-C/SGLI.
Figure 6. The location of each study site and the nearest pixel of GCOM-C/SGLI. The background true color image was created from Sentinel-2 level 2 products. The cyan circle indicates the location of the observation tower. The red square indicates a range of the nearest pixel of GCOM-C/SGLI.
Remotesensing 14 05352 g006
Figure 7. Time series of PRI original , PRI simulated , and PRI satellite from 2018 to 2020 at the four study sites. The blue square is PRI original , the red plus is PRI simulated , the black circle is PRI satellite not screened with the quality assessment (QA) flag, and the black cross is PRI original screened with the QA flag. The gray bands are snow seasons, the orange bands are autumn colors seasons, the green dotted lines are leafing timings, and the red dotted lines are leaf falling timings.
Figure 7. Time series of PRI original , PRI simulated , and PRI satellite from 2018 to 2020 at the four study sites. The blue square is PRI original , the red plus is PRI simulated , the black circle is PRI satellite not screened with the quality assessment (QA) flag, and the black cross is PRI original screened with the QA flag. The gray bands are snow seasons, the orange bands are autumn colors seasons, the green dotted lines are leafing timings, and the red dotted lines are leaf falling timings.
Remotesensing 14 05352 g007
Figure 8. Scatter plots between PRI satellite and PRI simulated . PRI satellite was screened with the QA flag (QA = 2 ). The dotted line represents the 1:1 line. The black solid line is the linear regression line. The shape of each point represents the year: the circle is 2018, the triangle is 2019, and the square is 2020. The color of each point corresponds to the DOY.
Figure 8. Scatter plots between PRI satellite and PRI simulated . PRI satellite was screened with the QA flag (QA = 2 ). The dotted line represents the 1:1 line. The black solid line is the linear regression line. The shape of each point represents the year: the circle is 2018, the triangle is 2019, and the square is 2020. The color of each point corresponds to the DOY.
Remotesensing 14 05352 g008
Figure 9. Time series of CCI original , CCI simulated , and CCI satellite from 2018 to 2020 at the four study sites. The blue square is CCI original , the red plus is CCI simulated , the black circle is CCI satellite not screened with the QA flag, and the black cross is CCI original screened with the QA flag. The gray bands are snow seasons, the orange bands are autumn colors seasons, the green dotted lines are leafing timings, and the red dotted lines are leaf falling timings.
Figure 9. Time series of CCI original , CCI simulated , and CCI satellite from 2018 to 2020 at the four study sites. The blue square is CCI original , the red plus is CCI simulated , the black circle is CCI satellite not screened with the QA flag, and the black cross is CCI original screened with the QA flag. The gray bands are snow seasons, the orange bands are autumn colors seasons, the green dotted lines are leafing timings, and the red dotted lines are leaf falling timings.
Remotesensing 14 05352 g009
Figure 10. Scatter plots between CCI satellite and CCI simulated . CCI satellite was screened with the QA flag (QA = 2 ). The dotted line represents the 1 : 1 line. The black solid line is the linear regression line. The shape of each point represents the year: the circle is 2018, the triangle is 2019, and the square is 2020. The color of each point corresponds to the DOY.
Figure 10. Scatter plots between CCI satellite and CCI simulated . CCI satellite was screened with the QA flag (QA = 2 ). The dotted line represents the 1 : 1 line. The black solid line is the linear regression line. The shape of each point represents the year: the circle is 2018, the triangle is 2019, and the square is 2020. The color of each point corresponds to the DOY.
Remotesensing 14 05352 g010
Figure 11. Scatter plots of snow-free PRI (ad) and snow-free CCI (eh). Satellite data were screened with the QA flag (QA = 2 ). Additionally, satellite data observed in snow seasons were manually removed. The dotted line represents the 1 : 1 line. The black solid line is the linear regression line. The shape of each point represents the year. The color of each point corresponds to the DOY.
Figure 11. Scatter plots of snow-free PRI (ad) and snow-free CCI (eh). Satellite data were screened with the QA flag (QA = 2 ). Additionally, satellite data observed in snow seasons were manually removed. The dotted line represents the 1 : 1 line. The black solid line is the linear regression line. The shape of each point represents the year. The color of each point corresponds to the DOY.
Remotesensing 14 05352 g011
Figure 12. The spatial distribution of PRI satellite and CCI satellite on 9 October 2019 (DOY = 282 ) around FHK. The cyan circles in each figure indicate the location of FHK. The red rectangle in (a) represents the range of (bd). (a) shows the true color image of GCOM-C/SGLI. (b) shows the true color image, (c) shows PRI satellite , and (d) shows CCI satellite . (a,b) are not screened with the QA flag and (c,d) are screened with the QA flag. The black area represents where the RSRF product was unavailable, and the white area represents the screened area with the QA flag.
Figure 12. The spatial distribution of PRI satellite and CCI satellite on 9 October 2019 (DOY = 282 ) around FHK. The cyan circles in each figure indicate the location of FHK. The red rectangle in (a) represents the range of (bd). (a) shows the true color image of GCOM-C/SGLI. (b) shows the true color image, (c) shows PRI satellite , and (d) shows CCI satellite . (a,b) are not screened with the QA flag and (c,d) are screened with the QA flag. The black area represents where the RSRF product was unavailable, and the white area represents the screened area with the QA flag.
Remotesensing 14 05352 g012
Figure 13. The sky images taken by upward ADFC around the observation time of GCOM-C/SGLI (10:46:44 (JST)) at FHK on 9 October 2019 (DOY = 282).
Figure 13. The sky images taken by upward ADFC around the observation time of GCOM-C/SGLI (10:46:44 (JST)) at FHK on 9 October 2019 (DOY = 282).
Remotesensing 14 05352 g013
Figure 14. The spatial distribution of PRI satellite and CCI satellite on 2 August 2019 (DOY = 214 ) around TKY. The cyan circles in each figure indicate the location of TKY. The red rectangle in (a) represents the range of (bd). (a) shows the true color image of GCOM-C/SGLI. (b) shows the true color image, (c) shows PRI satellite , and (d) shows CCI satellite . (a,b) are not screened with the QA flag and (c,d) are screened with the QA flag.
Figure 14. The spatial distribution of PRI satellite and CCI satellite on 2 August 2019 (DOY = 214 ) around TKY. The cyan circles in each figure indicate the location of TKY. The red rectangle in (a) represents the range of (bd). (a) shows the true color image of GCOM-C/SGLI. (b) shows the true color image, (c) shows PRI satellite , and (d) shows CCI satellite . (a,b) are not screened with the QA flag and (c,d) are screened with the QA flag.
Remotesensing 14 05352 g014
Figure 15. The spatial distribution of (a) VN 05 satellite , (b) VN 06 satellite , and (c) VN 08 satellite on 2 August 2019 (DOY = 214 ) around TKY. The range of each map is the same as that of Figure 14b–d. All maps are screened with the QA flag. The cyan circles in each figure indicate the location of TKY.
Figure 15. The spatial distribution of (a) VN 05 satellite , (b) VN 06 satellite , and (c) VN 08 satellite on 2 August 2019 (DOY = 214 ) around TKY. The range of each map is the same as that of Figure 14b–d. All maps are screened with the QA flag. The cyan circles in each figure indicate the location of TKY.
Remotesensing 14 05352 g015
Figure 16. The spatial distribution of PRI satellite and CCI satellite on 8 August 2018 (DOY = 220 ) around TSE. The cyan circles in each figure indicate the location of TSE. The red rectangle in (a) represents the range of (bd). (a) shows the true color image of GCOM-C/SGLI. (b) shows the true color image, (c) shows PRI satellite , and (d) shows CCI satellite . (a,b) are not screened with the QA flag and (c,d) are screened with the QA flag.
Figure 16. The spatial distribution of PRI satellite and CCI satellite on 8 August 2018 (DOY = 220 ) around TSE. The cyan circles in each figure indicate the location of TSE. The red rectangle in (a) represents the range of (bd). (a) shows the true color image of GCOM-C/SGLI. (b) shows the true color image, (c) shows PRI satellite , and (d) shows CCI satellite . (a,b) are not screened with the QA flag and (c,d) are screened with the QA flag.
Remotesensing 14 05352 g016
Figure 17. The spatial distribution of (a) VN 05 satellite , (b) VN 06 satellite , and (c) VN 08 satellite on 8 August 2018 (DOY = 220 ) around TSE. The range of each map is the same as that of Figure 16b–d. All maps are screened with the QA flag. The cyan circles in each figure indicate the location of TSE.
Figure 17. The spatial distribution of (a) VN 05 satellite , (b) VN 06 satellite , and (c) VN 08 satellite on 8 August 2018 (DOY = 220 ) around TSE. The range of each map is the same as that of Figure 16b–d. All maps are screened with the QA flag. The cyan circles in each figure indicate the location of TSE.
Remotesensing 14 05352 g017
Figure 18. Scatter plots of spatial mean of PRI satellite (ad) and CCI satellite (eh) calculated from four neighbor pixels. Satellite data were screened with the QA flag (QA = 2 ). The error bar represents standard deviation for each point. The dotted line represents the 1 : 1 line. The black solid line is the linear regression line. The shape of each point represents the year. The color of each point corresponds to the DOY.
Figure 18. Scatter plots of spatial mean of PRI satellite (ad) and CCI satellite (eh) calculated from four neighbor pixels. Satellite data were screened with the QA flag (QA = 2 ). The error bar represents standard deviation for each point. The dotted line represents the 1 : 1 line. The black solid line is the linear regression line. The shape of each point represents the year. The color of each point corresponds to the DOY.
Remotesensing 14 05352 g018
Figure 19. The downward fisheye images taken by ADFC in 2018 at (a) TSE, (b) TKY, (c) FJY, and (d) FHK. The numbers under the images are DOY when the images were taken. The images approximately display the approximately same observation area as downward MS-700 for each site.
Figure 19. The downward fisheye images taken by ADFC in 2018 at (a) TSE, (b) TKY, (c) FJY, and (d) FHK. The numbers under the images are DOY when the images were taken. The images approximately display the approximately same observation area as downward MS-700 for each site.
Remotesensing 14 05352 g019
Table 1. Specifications of four study sites.
Table 1. Specifications of four study sites.
Site
ID
Site
Name
Vegetation
Type
Latitude,
Longitude,
and Elevation
(WGS84)
Köppen–Geiger
Climate
Classification
[47,48]
Canopy
Height
Dominant Species
TSETeshioDeciduous
Needleleaf
Forest
45°3′20.99″N,
142°6′25.72″E,
70 m
Dfb10 m Hybrid larch
(Larix kaempferi ×
L. gmelinii),
Sasa senanensis,
and S. kurilensis
TKYTakayamaDeciduous
Broadleaf
Forest
36°8′42.79″N,
137°25′24.54″E,
1420 m
Dfb15–18  m Quercus crispula,
Betula platyphylla
Sukatchev var. japonica
Hara, B. ermanii,
and S. senanensis
FJYFuji YoshidaEvergreen
Needleleaf
Forest
35°27′16.36″N,
138°45′44.10″E,
1030 m
Cfb20 m Pinus densiflora,
Q. crispula, and Q. serrata
FHKFuji HokurokuDeciduous
Needleleaf
Forest
35°26′36.88″N,
138°45′52.93″E,
1100 m
Cfb25 m L. kaempferi,
P. densiflora,
Cornus controversa,
and Q. crispula
Table 2. The vertical positional information of the canopy and downward MS-700 at each site.
Table 2. The vertical positional information of the canopy and downward MS-700 at each site.
Site IDThe Height at Where
MS-700 Was Installed
The Distance between
the Canopy and MS-700
TSE23 m13 m
TKY18 m0–5  m
FJY28 m8 m
FHK32 m7 m
Table 3. Data acquisition intervals of MS-700 and Automatic-capturing Digital Fisheye Camera (ADFC) in the study sites. The time is described according to Japan Standard Time (JST).
Table 3. Data acquisition intervals of MS-700 and Automatic-capturing Digital Fisheye Camera (ADFC) in the study sites. The time is described according to Japan Standard Time (JST).
Site IDMS-700 Upward
(To the Sky)
MS-700 Downward
(To the Vegetations)
ADFC Downward
(To the Vegetations)
TSEEvery 1 min
04:00–19:59
Every 1 min
04:00–19:59
12:00
TKYLiner interpolation between
1 min before and 2 min after
the downward observation
Every 10 min
09:10–15:00
Every 15 min
07:00–16:45
FJYEvery 10 min
04:00–20:00
Every 10 min
04:00–20:00
None
FHKEvery 2 min
06:01–18:59
Every 4 min
06:03–18:59
Every 1 h
06:00–18:00
Table 4. Specifications of Visible and Near-Infrared Radiometer (VNR) of Second Generation Global Imager (SGLI) onboard Global Change Observation Mission-Climate (GCOM-C). The following information is based on GCOM-C data users handbook officially published by Japan Aerospace Exploration Agency (JAXA) [59].
Table 4. Specifications of Visible and Near-Infrared Radiometer (VNR) of Second Generation Global Imager (SGLI) onboard Global Change Observation Mission-Climate (GCOM-C). The following information is based on GCOM-C data users handbook officially published by Japan Aerospace Exploration Agency (JAXA) [59].
BandCenter WavelengthBand WidthSaturation LevelInstantaneous Field
Number[nm][nm][ W m 2 sr 1 μ m 1 ]of View (IFOV) [m]
VN01379.910.6240–241250
VN02412.310.3305–318250
VN03443.310.1457–467250
VN04490.010.3147–150250
VN05529.719.1361–364250
VN06566.119.895–96250
VN07672.322.069–70250
VN08672.421.9213–217250
VN09763.111.4351–359250
VN10867.120.937–38250
VN11867.420.8305–306250
Table 5. The details of the Quality Assessment (QA) flag of the level 2 atmospheric corrected land surface reflectance (RSRF) daily products of GCOM-C/SGLI (version 3) [61]. Bit 0 is the least significant bit.
Table 5. The details of the Quality Assessment (QA) flag of the level 2 atmospheric corrected land surface reflectance (RSRF) daily products of GCOM-C/SGLI (version 3) [61]. Bit 0 is the least significant bit.
Bit NumberDescriptionValue = 0Value = 1
0No dataNoYes
1Ocean or landOceanLand
2CoastNoYes
3Sun glint  > 0.005 NoYes
4Sun glint  > 0.12 NoYes
5Detection of snow or iceNoYes
6Cloud by target day estimationNoYes
7Probably cloud by multi day estimationNoYes
8Optical thickness  > 0.8 NoYes
9SaturatedNoYes
10The number of bidirectional reflectance factor (BRF) samples  3 NoYes
11Stray lightNoYes
12ShadowNoYes
13Detection of cloud or thick aerosol for polarization channelsNoYes
14Recovery of the data with previous days observation (for non-polarization bands)NoYes
15Recovery of the data with previous days observation (for polarization bands)NoYes
Table 6. The statistics of the accuracy assessment between PRI satellite and PRI simulated . n is the sample size, r is the coefficient of correlation, RMSE is the root mean square error, and MAE is the mean absolute error.
Table 6. The statistics of the accuracy assessment between PRI satellite and PRI simulated . n is the sample size, r is the coefficient of correlation, RMSE is the root mean square error, and MAE is the mean absolute error.
Site IDnrRMSEMAE
TSE84 0.289 ( p = 0.01 )0.0480.031
TKY40 0.245 ( p = 0.129 )0.1240.084
FJY146 0.180 ( p < 0.05 )0.0930.066
FHK65 0.294 ( p < 0.05 )0.0850.049
Table 7. The statistics of the accuracy assessment between CCI satellite and CCI simulated .
Table 7. The statistics of the accuracy assessment between CCI satellite and CCI simulated .
Site IDnrRMSEMAE
TSE84 0.868 ( p < 0.001 )0.1060.086
TKY40 0.911 ( p < 0.001 )0.0790.058
FJY146 0.775 ( p < 0.001 )0.0840.065
FHK65 0.874 ( p < 0.001 )0.1120.083
Table 8. The statistics of the accuracy assessment for snow-free PRI satellite and snow-free CCI satellite .
Table 8. The statistics of the accuracy assessment for snow-free PRI satellite and snow-free CCI satellite .
Site IDnrRMSEMAE
PRITSE81 0.166 ( p = 0.138 )0.0490.031
TKY29 0.274 ( p = 0.150 )0.1260.078
FJY115 0.228 ( p < 0.05 )0.0900.063
FHK65 0.294 ( p < 0.05 )0.0850.049
CCITSE81 0.863 ( p < 0.001 )0.1080.088
TKY29 0.915 ( p < 0.001 )0.0890.065
FJY115 0.847 ( p < 0.001 )0.0840.064
FHK65 0.874 ( p < 0.001 )0.1120.083
Table 9. The statistics of the accuracy assessment of PRI satellite and CCI satellite with the mean value of four neighbor pixels.
Table 9. The statistics of the accuracy assessment of PRI satellite and CCI satellite with the mean value of four neighbor pixels.
Site IDnrRMSEMAE
PRITSE69 0.339 ( p < 0.01 )0.0380.028
TKY18 0.596 ( p < 0.01 )0.0390.034
FJY120 0.217 ( p < 0.05 )0.0650.050
FHK53 0.088 ( p = 0.531 )0.0450.034
CCITSE69 0.874 ( p < 0.001 )0.1040.088
TKY18 0.884 ( p < 0.001 )0.0870.064
FJY120 0.886 ( p < 0.001 )0.0710.058
FHK53 0.901 ( p < 0.001 )0.0850.070
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sasagawa, T.; Akitsu, T.K.; Ide, R.; Takagi, K.; Takanashi, S.; Nakaji, T.; Nasahara, K.N. Accuracy Assessment of Photochemical Reflectance Index (PRI) and Chlorophyll Carotenoid Index (CCI) Derived from GCOM-C/SGLI with In Situ Data. Remote Sens. 2022, 14, 5352. https://doi.org/10.3390/rs14215352

AMA Style

Sasagawa T, Akitsu TK, Ide R, Takagi K, Takanashi S, Nakaji T, Nasahara KN. Accuracy Assessment of Photochemical Reflectance Index (PRI) and Chlorophyll Carotenoid Index (CCI) Derived from GCOM-C/SGLI with In Situ Data. Remote Sensing. 2022; 14(21):5352. https://doi.org/10.3390/rs14215352

Chicago/Turabian Style

Sasagawa, Taiga, Tomoko Kawaguchi Akitsu, Reiko Ide, Kentaro Takagi, Satoru Takanashi, Tatsuro Nakaji, and Kenlo Nishida Nasahara. 2022. "Accuracy Assessment of Photochemical Reflectance Index (PRI) and Chlorophyll Carotenoid Index (CCI) Derived from GCOM-C/SGLI with In Situ Data" Remote Sensing 14, no. 21: 5352. https://doi.org/10.3390/rs14215352

APA Style

Sasagawa, T., Akitsu, T. K., Ide, R., Takagi, K., Takanashi, S., Nakaji, T., & Nasahara, K. N. (2022). Accuracy Assessment of Photochemical Reflectance Index (PRI) and Chlorophyll Carotenoid Index (CCI) Derived from GCOM-C/SGLI with In Situ Data. Remote Sensing, 14(21), 5352. https://doi.org/10.3390/rs14215352

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