|
1.IntroductionEvapotranspiration (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.12–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. 2.Study Area and Datasets2.1.Study AreaOur 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.29–31 2.2.DatasetsThe 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. 2.2.1.Remote sensing dataTwo 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. Table 1Datasets used in this study.
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). 2.2.2.Data from the HiWATER experimentThe 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. Table 2Details of the LAS systems in HiWATER-MUSOEXE.
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 3.MethodologyWe 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. 3.1.Pixel Evapotranspiration AlgorithmThe SEB describes the transfer of energy between the land surface and the atmosphere. The energy budget is commonly expressed as follows: where is the net radiation (), is the soil heat flux (), is the sensible heat flux (), and is the latent heat flux (). Sensible heat flux is the flux caused by the temperature gradient between the land and the air, whereas LE is the change in energy due to ET or condensation. Atmospheric turbulence at the boundary layer causes the heat transfer at the land surface and atmospheric boundary. LE can be derived from the residual of the energy balance equation when , , and are known. 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 where is the downward shortwave radiation (), is the shortwave radiation reflected by the surface (), is the downward atmospheric longwave radiation (), and is the upward longwave radiation emitted by the surface (). is the surface albedo, is the emissivity of the land surface, , is the effective emissivity of the atmosphere, is the Stefan–Boltzmann constant, is the air temperature, and is the surface radiation temperature.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 where and represent the ratio of and with full canopy coverage and bare soil.is the turbulent heat transfer between the surface and the atmosphere that is derived using the temperature gradient between them. According to gradient diffusion theory where is the air density (), is the specific heat of air at a constant pressure, and is the radiative surface temperature (K). is the air temperature at a reference height (K), and is the aerodynamic resistance (), which is determined using the zero-plane displacement height and the roughness length. was calculated based on MOST using a stability correction function.38–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:3.2.Accuracy Evaluation IndexTo 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 () where represents the observed value, represents the estimated value, and represents the number of observations.The ground-based land surface temperature () was calculated using the Stefan–Boltzmann law from the AMS measurements of the longwave radiation fluxes25 where and are the in situ surface upwelling and atmospheric downwelling longwave radiations, respectively, and is the surface broadband emissivity.4.Results and Discussion4.1.Evaluation of Surface VariablesThe 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 Table 3The station validation of VIIRS LST.
Table 4The station validation of MODIS LST.
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. 4.2.Validation of Sensible Heat FluxesFigure 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 ). 4.3.Error and Sensitivity AnalysisMany 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.
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. 5.ConclusionWe 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:
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. AcknowledgmentsWe 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. ReferencesW. P. Kustas et al.,
“Estimating subpixel surface temperatures and energy fluxes from the vegetation index-radiometric temperature relationship,”
Remote Sens. Environ., 85
(4), 429
–440
(2003). http://dx.doi.org/10.1016/S0034-4257(03)00036-1 Google Scholar
V. M. Bindhu, B. Narasimhan and K. P. Sudheer,
“Development and verification of a non-linear disaggregation method (NL-DisTrad) to downscale MODIS land surface temperature to the spatial scale of Landsat thermal data to estimate evapotranspiration,”
Remote Sens. Environ., 135
(8), 118
–129
(2013). http://dx.doi.org/10.1016/j.rse.2013.03.023 Google Scholar
X. J. Li et al.,
“Estimating subpixel surface heat fluxes through applying temperature-sharpening methods to MODIS data,”
Remote Sens., 9
(8), 836
(2017). http://dx.doi.org/10.3390/rs9080836 Google Scholar
Y. Q. Wang et al.,
“Is scale really a challenge in evapotranspiration estimation? A multi-scale study in the Heihe oasis using thermal remote sensing and the three-temperature model,”
Agric. For. Meteorol., 230 128
–141
(2016). http://dx.doi.org/10.1016/j.agrformet.2016.03.012 0168-1923 Google Scholar
M. Bisquert et al.,
“Estimating high resolution evapotranspiration from disaggregated thermal images,”
Remote Sens. Environ., 187 423
–433
(2016). http://dx.doi.org/10.1016/j.rse.2016.10.049 Google Scholar
A. N. French, D. J. Hunsaker and K. R. Thorp,
“Remote sensing of evapotranspiration over cotton using the TSEB and METRIC energy balance models,”
Remote Sens. Environ., 158 281
–294
(2015). http://dx.doi.org/10.1016/j.rse.2014.11.003 Google Scholar
S. M. Liu et al.,
“Measurements of evapotranspiration from eddy-covariance systems and large aperture scintillometers in the Hai River Basin, China,”
J. Hydrol., 487
(9), 24
–38
(2013). http://dx.doi.org/10.1016/j.jhydrol.2013.02.025 JHYDA7 0022-1694 Google Scholar
Z. Z. Jia et al.,
“Validation of remotely sensed evapotranspiration over the Hai River Basin, China,”
J. Geophys. Res., 117
(D13), 13113
(2012). http://dx.doi.org/10.1029/2011JD017037 Google Scholar
W. M. L. Meijninger et al.,
“Determination of area-averaged water vapour fluxes with large aperture and radio wave scintillometers over a heterogeneous surface-Flevoland field experiment,”
Boundary Layer Meteorol., 105
(1), 63
–83
(2002). http://dx.doi.org/10.1023/A:1019683616097 Google Scholar
Q. Z. Mu et al.,
“Development of a global evapotranspiration algorithm based on MODIS and global meteorology data,”
Remote Sens. Environ., 111
(4), 519
–536
(2007). http://dx.doi.org/10.1016/j.rse.2007.04.015 Google Scholar
Q. Z. Mu, M. S. Zhao and S. W. Running,
“Improvements to a MODIS global terrestrial evapotranspiration algorithm,”
Remote Sens. Environ., 115
(8), 1781
–1800
(2011). http://dx.doi.org/10.1016/j.rse.2011.02.019 Google Scholar
Y. F. Jin, J. T. Randerson and M. L. Goulden,
“Continental-scale net radiation and evapotranspiration estimated using MODIS satellite observations,”
Remote Sens. Environ., 115
(9), 2302
–2319
(2011). http://dx.doi.org/10.1016/j.rse.2011.04.031 Google Scholar
J. B. Fisher, K. P. Tu and D. D. Baldocchi,
“Global estimates of the land-atmosphere water flux based on monthly AVHRR and ISLSCP-II data, validated at 16 FLUXNET sites,”
Remote Sens. Environ., 112
(3), 901
–919
(2008). http://dx.doi.org/10.1016/j.rse.2007.06.025 Google Scholar
Z. Q. Peng et al.,
“Remote sensing algorithm for surface evapotranspiration considering landscape and statistical effects on mixed pixels,”
Hydrol. Earth Syst. Sci., 20
(11), 4409
–4438
(2016). http://dx.doi.org/10.5194/hess-20-4409-2016 Google Scholar
Z. B. Su,
“The surface energy balance system (SEBS) for estimation of turbulent heat fluxes,”
Hydrol. Earth Syst. Sci., 6
(1), 85
–100
(2002). http://dx.doi.org/10.5194/hess-6-85-2002 Google Scholar
W. G. M. Bastiaanssen et al.,
“A remote sensing surface energy balance algorithm for land (SEBAL).: Part 2: validation,”
J. Hydrol., 212 213
–229
(1998). http://dx.doi.org/10.1016/S0022-1694(98)00254-6 JHYDA7 0022-1694 Google Scholar
W. J. Shuttleworth and J. S. Wallace,
“Evaporation from sparse crops-an energy combination theory,”
Q. J. R. Meteorolog. Soc., 111
(469), 839
–855
(1985). http://dx.doi.org/10.1002/qj.49711146910 QJRMAM 0035-9009 Google Scholar
J. M. Norman, W. P. Kustas and K. S. Humes,
“Source approach for estimating soil and vegetation energy fluxes in observations of directional radiometric surface temperature,”
Agric. For. Meteorol., 77
(3), 263
–293
(1995). http://dx.doi.org/10.1016/0168-1923(95)02265-Y 0168-1923 Google Scholar
T. Carlson,
“An overview of the ‘triangle method’ for estimating eurface evapotranspiration and soil moisture from satellite imagery,”
Sensors, 7
(8), 1612
–1629
(2007). http://dx.doi.org/10.3390/s7081612 SNSRES 0746-9462 Google Scholar
L. Jiang and S. Islam,
“Estimation of surface evaporation map over southern Great Plains using remote sensing data,”
Water Resour. Res., 37
(2), 329
–340
(2001). http://dx.doi.org/10.1029/2000WR900255 WRERAQ 0043-1397 Google Scholar
M. Yebra et al.,
“Evaluation of optical remote sensing to estimate actual evapotranspiration and canopy conductance,”
Remote Sens. Environ., 129
(2), 250
–261
(2013). http://dx.doi.org/10.1016/j.rse.2012.11.004 Google Scholar
K. Wang and S. L. Liang,
“An improved method for estimating global evapotranspiration based on satellite determination of surface net radiation, vegetation index, temperature, and soil moisture,”
J. Hydrometeorol., 9
(4), 712
–727
(2008). http://dx.doi.org/10.1175/2007JHM911.1 Google Scholar
S. M. Bateni and S. L. Liang,
“Estimating surface energy fluxes using a dual-source data assimilation approach adjoined to the heat diffusion equation,”
J. Geophys. Res., 117
(117), 17118
(2012). http://dx.doi.org/10.1029/2012JD017618 JGREA2 0148-0227 Google Scholar
Y. Liu et al.,
“Using data from Landsat, MODIS, VIIRS and PhenoCams to monitor the phenology of California oak/grass savanna and open grassland across spatial scales,”
Agric. For. Meteorol., 237 311
–325
(2017). http://dx.doi.org/10.1016/j.agrformet.2017.02.026 0168-1923 Google Scholar
H. Li et al.,
“Evaluation of the VIIRS and MODIS LST products in an arid area of Northwest China,”
Remote Sens. Environ., 142
(1), 111
–121
(2014). http://dx.doi.org/10.1016/j.rse.2013.11.014 Google Scholar
D. D. Wang et al.,
“Direct estimation of land surface albedo from VIIRS data: algorithm improvement and preliminary validation,”
J. Geophys. Res.: Atmos., 118
(22), 12,577
–12,586
(2013). http://dx.doi.org/10.1002/2013JD020417 Google Scholar
S. B. Duan and Z. L. Li,
“Spatial downscaling of MODIS land surface temperatures using geographically weighted regression: case study in Northern China,”
IEEE Trans. Geosci. Remote Sens., 54
(11), 6458
–6469
(2016). http://dx.doi.org/10.1109/TGRS.2016.2585198 IGRSD2 0196-2892 Google Scholar
X. Li et al.,
“Heihe watershed allied telemetry experimental research (HiWATER): scientific objectives and experimental design,”
Bull. Am. Meteorol. Soc., 94 1145
–1160
(2013). http://dx.doi.org/10.1175/BAMS-D-12-00154.1 BAMIAT 0003-0007 Google Scholar
X. Li et al.,
“Watershed allied telemetry experimental research,”
J. Geophys. Res., 114
(D22), 22103
(2009). http://dx.doi.org/10.1029/2008JD011590 Google Scholar
S. M. Liu et al.,
“Upscaling evapotranspiration measurements from multi-site to the satellite pixel scale over heterogeneous land surfaces,”
Agric. For. Meteorol., 230 97
–113
(2016). http://dx.doi.org/10.1016/j.agrformet.2016.04.008 0168-1923 Google Scholar
Z. W. Xu et al.,
“Intercomparison of surface energy flux measurement systems used during the HiWATER-MUSOEXE,”
J. Geophys. Res. Atmos., 118
(23), 140
–157
(2013). http://dx.doi.org/10.1002/2013JD020260 JGRDE3 0148-0227 Google Scholar
C. Conrad et al.,
“Mapping and assessing water use in a Central Asian irrigation system by utilizing MODIS remote sensing products,”
Irrig Drainage Syst, 21
(3), 197
–218
(2007). http://dx.doi.org/10.1007/s10795-007-9029-z Google Scholar
T. Islam et al.,
“A physics-based algorithm for the simultaneous retrieval of land surface temperature and emissivity from VIIRS thermal infrared data,”
IEEE Trans. Geosci. Remote Sens., 55
(1), 563
–576
(2017). http://dx.doi.org/10.1109/TGRS.2016.2611566 IGRSD2 0196-2892 Google Scholar
L. Belsma and E. Gottshall,
“Joint Polar Satellite System (JPSS) VIIRS surface reflectance Algorithm Theoretical Basis Document (ATBD),”
https://jointmission.gsfc.nasa.gov/sciencedocs/2015-06/474-00034_ATBD-VIIRS-Surface-Reflectance_A.pdf. Google Scholar
E. Valor and V. Caselles,
“Mapping land surface emissivity from NDVI: application to European, African, and South American areas,”
Remote Sens. Environ., 57
(3), 167
–184
(1996). http://dx.doi.org/10.1016/0034-4257(96)00039-9 Google Scholar
Y. H. Ran et al.,
“Spatial representativeness and uncertainty of eddy covariance carbon flux measurements for upscaling net ecosystem productivity to the grid scale,”
Agric. For. Meteorol., 230 114
–127
(2016). http://dx.doi.org/10.1016/j.agrformet.2016.05.008 0168-1923 Google Scholar
X. Li et al.,
“A multiscale dataset for understanding complex eco-hydrological processes in a heterogeneous oasis system,”
Sci. Data, 4 170083
(2017). Google Scholar
B. Choudhury and J. Monteith,
“A four-layer model for the heat budget of homogeneous land surfaces,”
Q. J. R. Meteorol. Soc., 114
(480), 373
–398
(1988). http://dx.doi.org/10.1002/(ISSN)1477-870X QJRMAM 0035-9009 Google Scholar
C. A. Paulson,
“The mathematical representation of wind speed and temperature profiles in the unstable atmospheric surface layer,”
J. Appl. Meteorol., 9
(6), 857
–861
(1970). http://dx.doi.org/10.1175/1520-0450(1970)009<0857:TMROWS>2.0.CO;2 Google Scholar
S. K. Ambast, A. K. Keshari and A. K. Gosain,
“An operational model for estimating regional evapotranspiration through surface energy partitioning (RESEP),”
Int. J. Remote Sens., 23
(22), 4917
–4930
(2002). http://dx.doi.org/10.1080/01431160110114501 IJSEDK 0143-1161 Google Scholar
B. J. Choudhury, R. J. Reginato and S. B. Idso,
“An analysis of infrared temperature observations over wheat and calculation of latent heat flux,”
Agric. For. Meteorol., 37
(1), 75
–88
(1986). http://dx.doi.org/10.1016/0168-1923(86)90029-8 0168-1923 Google Scholar
W. P. Kustas,
“Estimates of evapotranspiration with a one- and two-layer model of heat transfer over partial canopy cover,”
J. Appl. Meteorol., 29
(8), 704
–715
(1990). http://dx.doi.org/10.1175/1520-0450(1990)029<0704:EOEWAO>2.0.CO;2 Google Scholar
J. M. Chen et al.,
“Comparison of terrestrial evapotranspiration estimates using the mass transfer and Penman–Monteith equations in land surface models,”
J. Geophys. Res. Biogeosci., 118
(4), 1715
–1731
(2013). http://dx.doi.org/10.1002/2013JG002446 Google Scholar
Z. M. Wan,
“New refinements and validation of the collection-6 MODIS land-surface temperature/emissivity product,”
Remote Sens. Environ., 140
(1), 36
–45
(2014). http://dx.doi.org/10.1016/j.rse.2013.08.027 Google Scholar
A. Neftel, C. Spirig and C. Ammann,
“Application and test of a simple tool for operational footprint evaluations,”
Environ. Pollut., 152
(3), 644
–652
(2008). http://dx.doi.org/10.1016/j.envpol.2007.06.062 Google Scholar
W. M. L. Meijninger et al.,
“Determination of area-averaged sensible heat fluxes with a large aperture scintillometer over a heterogeneous surface-Flevoland field experiment,”
Boundary Layer Meteorol., 105
(1), 37
–62
(2002). http://dx.doi.org/10.1023/A:1019647732027 Google Scholar
BiographyXiaojun 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. |