Evapotranspiration (ET) is an important parameter in the interactions among soil, vegetation, and the atmosphere, and it is central part of surface energy and water budgets.1,2 Approximately 46% to 53% of the incoming solar radiation absorbed by the surface is globally consumed by the ET process as latent heat.3,4 ET can provide useful information about the heat and vapor exchange between the land surface and the atmosphere. A thorough understanding of ET can improve water resources management at different scales; in particular, detailed ET estimations in agricultural areas can improve the detection of water stresses and aid in irrigation scheduling.5,6
Currently, the general equipment used to measure ET includes a weighing lysimeter, a Bowen-ratio energy balance system, eddy covariance systems (ECs), and a large aperture scintillometer (LAS). ECs and LASs are widely used around the world.7 ECs can directly measure heat, carbon dioxide flux, and water vapor, with a measurement scale of hundreds of meters.8 LASs can measure the sensible heat flux, and kilometer-scale ET can be acquired by combining this data with available energy data (net radiation minus surface soil heat flux).9
All the mentioned measurements are point-based data that are location-specific; thus, extrapolating these values to cover large regions is difficult. Moreover, the instruments needed to perform this generalization are expensive. With the modern advances in technology, remotely sensed data have revealed unprecedented large-scale atmospheric, oceanic, and land data. Remote sensing-based ET values offer a suitable data source for regional- to global-scale applications. Five types of methods have been developed to estimate ET from satellite images and meteorological variables. (1) The Penman–Monteith equation is used in the moderate resolution imaging spectrometer (MODIS) ET product MOD16.5,10,11 The inputs used in this model include the following: land cover, albedo, fraction of vegetation cover, and meteorological variables. (2) The Priestley–Taylor models expand the range of the Priestley–Taylor coefficient in the Priestley–Taylor equation or combine the physiological force factors with the energy component of ET.1213.–14 (3) The latent heat flux (LE) is estimated as a residual term based on the surface energy balance (SEB) equation. According to the partitioning of the sources and sinks of the soil–plant–atmosphere continuum, SEB-based models can be classified as single-15,16 or dual-source models.17,18 (4) The land surface temperature-vegetation index (LST-VI) method describes the dry and wet edges of the LST-VI region as having the minimum and maximum ET, respectively. These methods interpolate the available data and use the Priestley–Taylor or Penman–Monteith equations to calculate the LE.19,20 (5) Additional methods include statistical/empirical methods,21 the use of complementary models,22 and land-process models with data assimilation schemes.23
MODIS data have a daily return time (one overpass per day at varying view angles). Currently, such data are ideal for simulating regional-scale changes in surface energy and water fluxes. The high-resolution spectral configuration and improved radiometric and geographic accuracy of MODIS facilitate a high-quality global ET product given a frequency of 8 days.10,11 However, the MODIS sensors on the Terra and Aqua satellite platforms are aging; therefore, continuing to build this dataset in the future will require relying on data from the daily weather satellites represented by the Suomi National Polar-orbiting Partnership (NPP) and its successor.24 The visible infrared imaging radiometer suite (VIIRS) onboard the Suomi NPP satellite and the future Joint Polar Satellite System (JPSS) missions have similar designs to MODIS and thus provide a similar complement of visible and near-infrared and shortwave infrared channels to monitor the LST, albedo, and vegetation index products.25,26 Many websites provide these products, which lay the foundation for regional ET estimations. Moreover, there are few reports about the application of VIIRS data to calculate regional ET values. The purpose of this work is to investigate the reliability of VIIRS-derived surface variables for regional ET estimation.
Study Area and Datasets
Our study was conducted in the middle stream of the Heihe River Basin (HRB) at the Zhangye oasis, which is located near the city of Zhangye in the arid region of Gansu Province in northwestern China (100.10°E to 100.66°E, 38.68°N to 39.15°N). The Zhangye oasis has a semiarid climate, with an annual precipitation of and an annual potential ET of to 1800 mm; the annual average temperature is .27,28 Areas of the Gobi desert and alpine vegetation are located near the study area in the Qilian Mountains (see Fig. 1). Around the oasis, the land surface is mainly covered by maize, which has a growing season that lasts from May to September. The artificial oasis is highly heterogeneous, which impacts its thermodynamic and hydraulic features. Consequently, its water use efficiency and ET are variable. The HRB has long served as a test area for integrated watershed studies and land surface and hydrological experiments.28 One major objective of the Heihe Watershed Allied Telemetry Experimental Research (HiWATER), namely the Multi-Scale Observation Experiment on ET (HiWATER-MUSOEXE), is to capture strong land surface heterogeneities and the associated uncertainties within a watershed.2930.–31
The data in this study were mainly derived from VIIRS and MODIS. We combined these two sensors datasets with ancillary data and in situ HiWATER-MUSOEXE data to estimate and validate the VIIRS and MODIS heat fluxes.
Remote sensing data
Two kinds of remote sensing data were used in this study, specifically MODIS and VIIRS data. The MODIS sensor launched on December 18, 1999 aboard the Terra satellite has been in operation for nearly 18 years and provides global data that are used for monitoring long-term regional land surface heat and water vapor fluxes as well as other Earth-related changes.32 The MODIS data have high temporal (i.e., at least once every day) and spectral resolutions (i.e., 36 bands in total). The NPP/VIIRS was launched in late 2011 and has imaging capabilities in 22 bands with a spectral coverage from 0.412 to . Among the 22 bands, there are five high-resolution imaging channels (-bands), 16 moderate resolution channels (-bands), and 1 day/night band (DNB). The -bands have spatial resolutions of 375 m at nadir, and -bands (and DNB) have a resolution of 750 m at nadir. In the cross-track direction, the scan swath width is . The data for one granule covers of data and comprises 48 scans.33 The images produced by VIIRS are stable with balanced band performances.
In this paper, we mainly use two kinds of VIIRS operational products [i.e., VIIRS/NPP surface reflectance daily L2G global 1 km and 500 m SIN grid (VNP09GA) and VIIRS/NPP daily gridded land surface temperature 1-km SIN grid day (NPP_DLSTD_L3D)], which can be freely downloaded.34 The parameters from the MODIS sensor used for estimating heat fluxes, such as normalized difference vegetation index (NDVI), albedo, LSTs, and emissivity, are freely available from the NASA website. We used the MODIS and VIIRS image data that cover the HRB region from 2012. Many algorithms for estimating ET using surface parameters require clear-sky conditions; therefore, we combined the available information on data quality with visual interpretations to select satellite images without clouds. We obtained data for 9 days during the period of ground observations: June 19 and 30, July 8 and 27, August 3, 22, and 29, and September 2 and 14. These two kinds of satellite products used in this study are shown in Table 1.
Datasets used in this study.
|Product||Sensor||Spatial resolution (m)||Data interpretation||Location||Application|
|VNP09GA||NPP/VIIRS||500, 1000||Daily surface reflectance||Tile h25v05||Calculating NDVI, FVC, etc.|
|NPP_DLSTD_L3D||NPP/VIIRS||1000||Daily land surface temperature||Tile h25v05||Providing LST data|
|MOD11A1||MODIS||1000||Daily land surface temperature/emissivity||Tile h25v05||Providing LST data|
|MOD09GA||MODIS||500, 1000||Daily surface reflectance||Tile h25v05||Calculating NDVI, FVC, etc.|
All of the above data are required for the corresponding preprocessed. The MODIS products that cover the study area were downloaded and corrected for geometric and projection effects. VNP09GA is a level-2G surface reflectance product obtained by adjusting the top-of-atmosphere reflectance to compensate for atmospheric effects. The effects of molecular gases, including ozone and water vapor, and the effects of atmospheric aerosols are corrected. Previous studies showed that VIIRS LST products have reasonable accuracies.25 The VIIRS data covering the HRB were preprocessed using geometric correction, radiometric calibration, and atmospheric correction.14 To estimate the regional ET, various surface variables, including the NDVI, the surface albedo,26 the emissivity,35 the leaf area index (LAI), and the fractional vegetation coverage (FVC), were retrieved by applying the lumped method to the VIIRS data.14 These two kinds of satellite data were organized and paired with each other on the Universal Transverse Mercator (UTM) projection (i.e., WGS 84 UTM ZONE 47N).
Data from the HiWATER experiment
The in situ HRB observational data were provided by HiWATER. From June to September 2012, HiWATER designed nested observation matrices of and within the middle stream oasis.
For the larger observation matrix, four EC systems and one superstation (EC15) were installed in the oasis–desert ecosystem. Each station was supplemented with an automatic meteorological station (AMS) to record the meteorological and soil variables and to monitor the spatiotemporal variations of ET and its associated factors,28 including air pressure, precipitation, wind speed/direction, four-component radiation, and soil moisture. Each AMS has a different observation date, and its specific observation time and sensor characteristics can be found in Ref. 30. Within the oasis, 17 AMSs and EC towers were used, including one in a vegetable field (EC01), one in a residential area (EC04), one in an orchard (EC17), and 14 in maize fields. Because the HHZ station (shown in Fig. 1) and EC16 station lacked net radiations and soil heat flux observation data, they were excluded from this study. Additionally, four groups of optical LASs (eight sets, with two sets per group) were installed in the and MODIS pixels within the kernel experimental area (three groups in three MODIS pixels, named LAS1 to LAS3, from west to east, respectively, and one group in a MODIS pixels, LAS4). Considering the observational scale of the LAS systems and the resolution of satellite data used in this paper, we mainly use the EC data to validate the surface variables (i.e., land surface temperature, net radiation, and soil heat flux) and the LAS data to validate the sensible heat flux. The details of these LAS systems are shown in Table 2.
Details of the LAS systems in HiWATER-MUSOEXE.
|Station||Latitude (°N)||Longitude (°E)||Path length (m)||Height (m)||Type, manufactures|
|LAS1-N||38.8841||100.3510||3256||33.45||BLS900, Scintec Germany|
|LAS1-S||38.8547||100.3524||3256||33.45||zzlas, RR9340, Rainroot, China|
|LAS2-N||38.8827||100.3625||2841||33.45||BLS900, Scintec Germany|
|LAS2-S||38.8571||100.3618||2841||33.45||BLS450, Scintec Germany|
|LAS3-N||38.8834||100.3733||3111||33.45||BLS900, Scintec Germany|
|LAS3-S||38.8555||100.3722||3111||33.45||LAS, Kipp&zonen, Netherland|
|LAS4-N||38.8607||100.3785||1854||22.45||BLS450, Scintec Germany|
|LAS4-S||38.84682||100.3685||1854||22.45||zzlas, RR9340, Rainroot, China|
An LAS is a device that derives the turbulent intensity through measuring the refractive index of air, . The sensible heat flux (LAS-H) was obtained using Monin–Obukhov similarity theory (MOST) via an iterative method that combines various meteorological data (e.g., wind speed, air temperature, and air pressure).7,30 Finally, the LAS-H data were averaged across 30-min intervals. The data were carefully screened to ensure the LAS data quality: (i) the data were rejected when exceeded the saturated criterion, which was determined according to Ref. 36, (ii) the data were rejected when the demodulated signal was less than the threshold,30 (iii) the data were rejected when collected on nights with weak turbulence, and (iv) the data were rejected when precipitation occurred. Additionally, to avoid advection conditions, we selected data at 30-min intervals during periods when the values of the EC and LAS measurements were . The soil heat flux was measured using three soil heat plates at a depth of 6 cm at each site and the surface soil heat flux was calculated based on the soil temperature and moisture above the plates.14 The surface meteorological variables, such as wind direction, wind speed, air pressure, and relative humidity, were used to interpolate images using inverse distance weighting. The data obtained during the HiWATER experiment can be freely obtained from the Heihe Plan Data Management Center.37
We have developed a single-source ET estimation algorithm that combines sensor-derived surface variables with meteorological data over the northwestern China HRB area. For comparison, we use the MODIS-derived surface variables to calculate the regional ET values and the two ETs are compared with ground-based observations; the errors in the model are then analyzed and discussed.
Pixel Evapotranspiration Algorithm
The SEB describes the transfer of energy between the land surface and the atmosphere. The energy budget is commonly expressed as follows:Figure 2 shows the flowchart for ET retrieval using VIIRS data. The models that permit the estimations of , , and using remote sensing data are introduced.
, also known as the radiation balance or radiation budget, is the main energy source of land surface heat and mass transfers and exchanges, such as evaporation, and is calculated as the difference between the incoming and outgoing radiation
is the heat exchange between the surface and the deeper soils and is related to the exchanges of water and heat within the soil. Generally, heat is transferred from deeper depths to the surface at night. Studies have found that soil heat flux is correlated with and can be derived using land surface parameters.15 Because canopies exert significant influences on , the fractional canopy coverage is used to determine the ratio of to 14
is the turbulent heat transfer between the surface and the atmosphere that is derived using the temperature gradient between them. According to gradient diffusion theory3839.–40 Because remote sensing-based methods cannot obtain , the value of is usually replaced by the radiative surface temperature (). It should be noted that is not strictly equal to the in Eq. (4). The difference between these terms for homogeneous and full-coverage vegetation is to 2 K41 and it can reach 10 K in sparsely vegetative areas.42 The method used by Chen et al.43 corrects for this discrepancy by adding an “excess resistance” and thus, was used in this study. The modified can be specified as follows:
Accuracy Evaluation Index
To control model inputs and analyze error sources, the land surface temperature, net radiation, and soil heat flux were evaluated using in situ data. The analyses were performed using the root mean square error (RMSE), the mean bias error (MBE), and the coefficient of determination ()
The ground-based land surface temperature () was calculated using the Stefan–Boltzmann law from the AMS measurements of the longwave radiation fluxes25
Results and Discussion
Evaluation of Surface Variables
The station validation results of two LSTs are shown in Tables 3 and 4. The , MBE, and RMSE values of the VIIRS LST data are 0.77, 0.55 K, and 2.40 K, respectively. As shown in Table 3, the VIIRS MBE varies from to 1.95 K; the RMSE varies from 1.16 to 7.10 K. The VIIRS LSTs are warmer than the ground LSTs for most stations. From Table 4, the MBE values for the MODIS LST are negative for most stations, and the MBEs and RMSEs are much larger, with two stations (GB and SSW) showing and . From Fig. 1, the land cover of GB and SSW is uncultivated land. In our study area, the Gobi region and desert area comprise uncultivated land. Therefore, the MODIS LST products underestimate the LSTs for barren surfaces during the daytime, consistent with the results of Refs. 25 and 44. There are two possible reasons for the large Terra MODIS LST errors over barren surface sites. First, the LST product algorithm does not cover the wide range of LSTs, i.e., the daytime LSTs are greater than the air temperature at the surface level by , which is beyond the upper limit. Second, the large errors are present in the surface emissivity values in MODIS bands 31 and 32 estimated from land cover types.44
The station validation of VIIRS LST.
|Station||R2||MBE (K)||RMSE (K)||Station||R2||MBE (K)||RMSE (K)|
The station validation of MODIS LST.
|Station||R2||MBE (K)||RMSE (K)||Station||R2||MBE (K)||RMSE (K)|
As shown in Tables 3 and 4, the EC4 accuracy of both LST products is low. The main causes of the large errors are as follows: (1) buildings and soil/vegetation are distinct materials, and the emissivity algorithm may not be suitable for buildings and (2) the EC4 foundation is nonuniform and is not suitable for validation. After removing the EC4 data, the , MBE, and RMSE values of the VIIRS LST data are 0.86, 0.75 K, and 1.97 K, respectively, and the , MBE, and RMSE values of the MODIS LST data are 0.82, , and 2.77 K, respectively. The SD errors for both were large due to large errors on particular days. For example, although it was briefly cloudy above station SD on July 27, this area was not identified as cloudy in the cloud detection process.
The validation of net radiation and soil heat flux based on scatterplots is shown in Fig. 3. As shown in this figure, the , MBE, and RMSE values of the calculated using the VIIRS data were 0.68, , and , respectively, whereas the errors in the estimates obtained using the MODIS data were 0.64, , and , respectively. According to the sensitivity analysis of Eq. (2), and are highly sensitive variables when calculating , whereas the albedo, LST are not as sensitive. However, the and values used in the ET model are derived from ground observations; thus, the and are relatively accurate. Therefore, LST and albedo are the main sources of error and are discussed in Sec. 4.3. A 1 K bias in LST would result in an error of , but an albedo bias of 0.03 can lead to an error of when the incoming solar radiation is large. Studies have shown that land surface albedos retrieved using VIIRS data are better than those retrieved using MODIS;26 thus, the values calculated using the VIIRS data show better results.
The , MBE, and RMSE values of calculated using the VIIRS data were 0.58, , and , respectively, whereas the errors in the values calculated using the MODIS data were 0.52, , and , respectively. Since is calculated using , and show the same sensitivity results. As shown in the validation of in Fig. 3, the sensor-retrieved is overestimated at many sites, which may be related to the soil temperature and moisture above the soil heat plates at each site. For example, at EC5, the soil temperature and moisture were the same at different depths after July 19, which resulted in the surface being equal to the at a depth of 6 cm. However, the values below the surface are usually less than those at the soil surface; thus, the validation results of at EC5 indicate that was overestimated. However, is still underestimated as a whole () because is also underestimated.
Validation of Sensible Heat Fluxes
Figure 4 provides the turbulent heat flux results calculated from VIIRS and MODIS on June 19, 2012. The differences in the distributions of the heat fluxes calculated from the two sensors are not obvious because they follow certain underlying rules. Specifically, the spatial distributions of the turbulent heat fluxes are related to their underlying land covers. values are higher in areas covered by buildings and uncultivated land, including land patches in the desert areas, barren areas, and Gobi regions. However, values are lower in the water-covered and agricultural areas within the oasis. Comparing Figs. 4(e) and 4(f) and Figs. 4(g) and 4(h), we see that the values (LE values) of the barren areas in the southwest are lower (higher) than those of the desert areas in the southeast, although both areas were classified as uncultivated land. Because the barren areas in the southwest border of the Qilian Mountains, the groundwater levels are high due to snow melt and the downward movement of water, leading to a soil moisture content of according to in situ measurements at a depth of 2 cm.14 As shown in Figs. 4(e) and 4(f), some areas in the distribution map have values less than zero due to inversion from the oasis effect or irrigation. The HiWATER soil moisture data show that irrigation occurred on June 19, 2012. Irrigation is the main source of water within the oasis and cools the land surface to temperatures below that of the air.
To further analyze the heat fluxes retrieved from the two sensors, the values were evaluated using in situ data. In recent years, many researchers have employed a footprint model to validate remote sensing-based ET estimates (specifically only and ). When combined with the flux contribution source areas, this model solves the problems associated with mismatches between the observed values and the remote sensing-based estimates and leads to better results.45,46 In this study, we employ the path-weighting function of the LAS and the footprint model for point fluxes to calculate the LAS source area.30,46 The station-validated results for the values are shown in Fig. 5.
The footprint validation results show that the , MBE, and RMSE values of the values calculated using the VIIRS data were 0.68, , and , respectively, and those of the MODIS-based values were 0.56, , and , respectively. Compared with the observed values, the calculated using VIIRS data is overestimated (), whereas the calculated using MODIS data is underestimated. The reason is related to the LST. According to Eq. (4), is most sensitive to the LST values. A 1 K bias in LST would result in an error in of 20.5%. Consequently, LST is the most important factor affecting . However, this sensitivity to LST is unstable and depends on the strength of the turbulence present, as discussed in Sec. 4.3. The validation of the VIIRS and MODIS surface temperature products in Sec. 4.1 shows that the VIIRS LSTs are warmer than the ground LSTs (with an average bias of 0.55 K), whereas the Terra MODIS LSTs are significantly underestimated (with an average bias of ).
Error and Sensitivity Analysis
Many sources of error are unavoidable, including the calibration coefficients of the onboard satellite sensors; the drifts in the sensors caused by the aging of the instruments; mismatches in the timescales of the observational data, which are averaged over time, and the remotely sensed turbulent heat fluxes, which are instantaneous; and the geometric corrections cause a half-pixel bias less than or equal to the deviation of the artificially subjective interpretation. The estimation of is the most important and difficult step in the remote sensing-based assessment of the SEB; thus, the sensitivity of was analyzed first. Land surface variables, including LST, canopy height, FVC, and LAI, and meteorological variables, including wind speed, air temperature, and relative humidity, are the major factors in the sensitivity analysis of . Figure 6 presents a sensitivity analysis of . In this case, the average LST is 304 K, with a range from 297.5 to 311 K, and has a step size of 0.5 K; meanwhile, LAI ranges from 0.1 to 2.5 and has a step size of 0.1. The canopy height is 1 m, with a range from 0.1 to 1.9 m. The wind speed ranges from 0.50 to , the air temperature is , , and the relative humidity is . In addition, the land cover is maize and the reference is .
The air pressure is stable over short periods and has little effect on the ET results. Although, is sensitive to meteorological variables, such as wind speed and air temperature (see Fig. 6), the meteorological data used in the ET model are derived from ground observations; thus, the meteorological factors are relatively accurate. As shown in Fig. 6, LAI, LST, and the canopy height are all influential variables. is sensitive to LAI when LAI is . The momentum roughness length increases as LAI increases and the turbulent exchange is enhanced. However, when LAI is , the plant canopy is regarded as a continuum that is not a sensitive variable. Because our study area is dominated by agriculture and the study period extended from June to September, the crops in the middle reaches of the HRB grew quickly; thus, LAI was usually . Hence, LST and canopy height are the main sources of error.
1. Errors in LST. As shown in Fig. 6, a 1 K bias in LST would result in 20.5% error in when is , the corresponding errors in and are both only . Therefore, in this paper, LST is mainly affected by the impact of on the role of heat fluxes. However, as we mentioned in Sec. 3, depends on the temperature gradient between the surface and the atmosphere. Hence, the LST sensitivity is unstable and depends on the strength of the local turbulence. The strength of the turbulence determines the mass and energy transports and the resistance to heat transfer, which influences the sensitivity of to the LST. A weaker turbulence corresponds to a weaker LST sensitivity and vice versa.
2. Errors in the canopy height. The canopy height was obtained from a phenophase and classification map. Thus, the accuracy of the canopy height depends mainly on the plant growth state and classification accuracies. Even within the same region, the canopy height of a crop can differ due to differences in seeding times and soil attributes, such as soil moisture and fertilization.
As discussed above, the accuracy of the LST and canopy height parameters greatly influences the estimated turbulent heat flux. Since the VIIRS sensor is still running, the algorithm to calculate its LST product is still being explored.25 Canopy heights are known a priori and are derived from phenophase classifications; these classifications influence the accuracies of the surface roughness calculations. Multisource remote sensing data, such as active microwave and LiDAR data, can be used to obtain canopy heights for future studies,14 which would lead to more robust canopy height values. Since a footprint model was used in the validation, discrepancies between the in situ measurements and remote sensing pixels still exist in the footprints.3
Additionally, to correct the discrepancies between the remotely sensed radiative surface temperatures and the aerodynamic temperatures at the source of heat transport, a brief and well-constrained parameterization scheme (under a uniform and flat plant surface) that included excess resistance was used to calculate the aerodynamic resistance of the heat transfer.14 Since the surface cover of the HRB is very heterogeneous,4,30 multiple parameterization methods should be compared to select an optimal method.43 Although these preliminary results show that VIIRS data have certain advantages over MODIS data when calculating surface heat fluxes, the limited ground observation data in this paper use only 9 days in calculating and validating the results. Future work will involve testing this single-source ET estimation algorithm and the practicality of VIIRS data for calculating the heat fluxes in the range of actual environmental conditions.
We use a single-source model and VIIRS data to calculate the components of SEB. To determine the usefulness VIIRS sensors in estimating regional heat fluxes, we also used MODIS data to calculate the regional ET. After data preparation, SEB component inversion, verification (, , and ), and error analysis, the following conclusions can be drawn:
1. We have developed a single-source ET estimation algorithm that combines sensor-derived surface variables with meteorological data over the HRB area of northwestern China. This method can be used to estimate the turbulent heat flux based on VIIRS sensor data, which has certain advantages (i.e., the MBE of the sensible heat flux is ) compared with using MODIS data (i.e., the MBE of the sensible heat flux is ).
2. To control model inputs and analyze error sources, the model parameters (including the land surface temperature, net radiation, and soil heat flux) were evaluated using in situ data. The validation results show that the VIIRS LST product exhibits a more consistent agreement with in situ measurements on barren surfaces compared with the Terra MODIS LST product.
3. To estimate ET using a single-source model and remote sensing data, the sensible heat flux is the most important component, but the calculation of this component is also complex. In our model, the accuracy of the sensible heat flux mainly depends on the inversion accuracy of the surface temperature and the canopy height.
The VIIRS data have considerable advantages due to their fine temporal resolutions and free access. Moreover, the expansion and improvement of MODIS data and the next-generation of JPSS platforms will provide a broader application space. Because the Suomi-NPP/VIIRS has been in orbit for several years, the long-term data are promising for applications in monitoring energy budgets.
We are grateful to all the scientists and engineers who took part in the HiWATER experiment. This study was jointly supported by the Chinese Natural Science Foundation Project (grant no. 41371360) and the Special Fund of the Chinese Academy of Sciences (grant no. KZZD-EW-TZ-18). Generous help in revising the paper was provided by the editors and reviewers.
Xiaojun Li received his BE degree in remote sensing science and technology in 2015 and is now an undergraduate student at the State Key Laboratory of Remote Sensing Science, jointly sponsored by the Institute of Remote Sensing and Digital Earth of Chinese Academy of Sciences and Beijing Normal University, Chinese Academy of Sciences. His current research interests are the theory and method of remote sensing of land surface radiation and energy balance.
Xiaozhou Xin received his PhD in geology from the Institute of Remote Sensing Applications, Chinese Academy of Sciences in 2003. Currently, he is an associate professor and has been researching in the theory and method of remote sensing of land surface radiation and energy balance for more than 10 years.
Zhiqing Peng received his MSc degree from the State Key Laboratory of Remote Sensing Science, jointly sponsored by the Institute of Remote Sensing and Digital Earth of Chinese Academy of Sciences and Beijing Normal University, Chinese Academy of Sciences in 2017. His current research interests include deep learning, remote sensing of land surface radiation, and energy balance.
Hailong Zhang received his PhD in cartography and geography information system from the Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing, China, in 2010. Currently, he is an assistant researcher at the Institute of Remote Sensing Applications, CAS. His research interests include shortwave radiation estimation from geostationary and polar satellites.
Li Li received her MSc degree in cartography and geography information system from Beijing Normal University, China, in 2005 and her PhD in cartography and geography information system from the Institute of Remote Sensing Applications, CAS, Beijing, China, in 2010. Her current research interests include PAR estimation, FPAR inversion based on remote sensing data, and its applications in agriculture and environment.
Shanshan Shao received her BE degree in educational science in 2015 and now is an undergraduate student at the College of Educational Science, Anhui Normal University. She is mainly focused on psychological research.
Qinhuo Liu received his PhD in atmospheric physics from Peking University, Beijing, China, in 1997. Currently, he is a professor and the deputy director of the State Key Laboratory of Remote Sensing Science with the Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences. His research interests include radiation transfer modeling for optical remote sensing and terrestrial parameter inversion from multisource remote sensing data.