Evapotranspiration (ET) is one of the components in the water cycle and the surface energy balance systems. It is fundamental information for agriculture, water resource management, and climate change research. This study presents a scheme for regional actual evapotranspiration estimation using multi-source satellite data to compute key land and meteorological variables characterizing land surface, soil, vegetation, and the atmospheric boundary layer. The algorithms are validated using ground observations from the Heihe River Basin of northwest China. Monthly data estimates at a resolution of 1 km from the proposed algorithms compared well with ground observation data, with a root mean square error (RMSE) of 0.80 mm and a mean relative error (MRE) of −7.11%. The overall deviation between the average yearly ET derived from the proposed algorithms and ground-based water balance measurements was 9.44% for a small watershed and 1% for the entire basin. This study demonstrates that both accuracy and spatial depiction of actual evapotranspiration estimation can be significantly improved by using multi-source satellite data to measure the required land surface and meteorological variables. This reduces dependence on spatial interpolation of ground-derived meteorological variables which can be problematic, especially in data-sparse regions, and allows the production of region-wide ET datasets.

Net radiation plays an essential role in determining the thermal conditions of the Earth’s surface and is an important parameter for the study of land-surface processes and global climate change. In this paper, an improved satellite-based approach to estimate the daily net radiation is presented, in which sunshine duration were derived from the geostationary meteorological satellite (FY-2D) cloud classification product, the monthly empirical a_s and b_s Angstrom coefficients for net shortwave radiation were calibrated by spatial fitting of the ground data from 1997 to 2006, and the daily net longwave radiation was calibrated with ground data from 2007 to 2010 over the Heihe River Basin in China. The estimated daily net radiation values were validated against ground data for 12 months in 2008 at four stations with different underlying surface types. The average coefficient of determination (R²) was 0.8489, and the averaged Nash-Sutcliffe equation (NSE) was 0.8356. The close agreement between the estimated daily net radiation and observations indicates that the proposed method is promising, especially given the comparison between the spatial distribution and the interpolation of sunshine duration. Potential applications include climate research, energy balance studies and the estimation of global evapotranspiration.

As a vital role in the processes of the energy balance and hydrological cycles, actual evapotranspiration (ET) is relevant to many agricultural, ecological and water resource management studies. The available global or regional ET products provide ET estimations with various temporal ranges, spatial resolutions and calculation methods (algorithms, inputs and parameterization, etc.), leading to varying degrees of introduced uncertainty. Northern China is the main agriculturally productive region supporting the whole country; thus, understanding the spatial and temporal changes in ET is essential to ensure water resource and food security. We developed a synthesis ET dataset for Northern China at a 1000 m spatial resolution, with a monthly temporal resolution covering a period ranging from 1982 to 2017, using an in-depth assessment of several ET products. Specifically, assessments were performed using in situ measured ET from eddy covariance (EC) observation towers at the site-pixel scale over interannual months under the conditions of different land cover types, climatic zones and elevation levels to select the most optimally performing ET products to be used in the synthesized ET dataset. Eight indicators under 21 conditions were involved in the assessment sheet, while the statistics of the different ET product occurrences and corresponding ratios were analyzed to select the best-performing ET products to build the synthesis ET dataset using the weighted mean method. The weights were determined by the Taylor skill score (TSS), calculated with ET products and EC ET observation data. Based on the assessment results, the Penman–Monteith–Leuning (PML_v2), ETWatch and Operational Simplified Surface Energy Balance (SSEBop) datasets were selected for implementation in the synthesis ET dataset from 2003 to 2017, while Global Land Evaporation Amsterdam Model (GLEAM) v3.3a, complementary relationship (CR) ET, and Numerical Terradynamic Simulation Group (NTSG) datasets were chosen for the synthesis ET dataset from 1982 to 2002. The weighted mean synthesized results from 2003 to 2017 performed well when compared to the in situ measured EC ET values produced under all of the above conditions, while the synthesized results from 1982 to 2002 performed well through the water balance method in Heihe River Basin. These results can provide more stable ET estimations for Northern China, which can contribute to relevant agricultural, ecological and hydrological studies.

Sunshine duration is an important indicator of the amount of solar radiation received in a region and an important input parameter for the study of atmospheric energy balance, climate change, ecosystem evolution, and social sustainability. Currently, extrapolation and interpolation of data from meteorological stations are the most common methods used to calculate sunshine duration on a regional scale. However, it is difficult to obtain high precision sunshine duration in areas lacking ground observation or where sunshine duration is highly heterogeneous on the ground. In this paper, a new method is proposed to estimate sunshine duration with hourly total cloud amount (CTA) data from sunrise to sunset derived from the Fengyun-2G geostationary meteorological satellite (FY-2G). This method constructs a new index known as daytime mean total cloud coverage amount and provides quadratic equations relating daytime mean total cloud coverage amount to relative sunshine duration in different seasons. The method was validated with ground observation data for 2016 from 18 meteorological stations in the Three-River Headwaters Region of Qinghai Province, China. For individual stations, the coefficient of determination (R²) between estimated and measured sunshine was at least 0.894, the RMSE (root mean square error) was 0.977 h/day or less, the MAE (mean absolute error) was 0.824 h/day or less, the RE (relative error) was 0.150 or lower, and the value of d was 0.963 or greater, which validated that the proposed method can effectively predict daily sunshine duration. These equations can also provide higher precision estimates of regional-scale sunshine duration. This was demonstrated by comparing, for the entire study region, the spatial distribution of sunshine duration estimated from season-based equations with results from three different interpolation methods based on ground observations. Overall, the study confirms that total cloud amount measures from a geostationary satellite can be used to successfully estimate sunshine duration.

Sunshine duration is an important variable that is widely used in atmospheric energy balance studies, analysis of the thermal loadings on buildings, climate research, and the evaluation of agricultural resources. In most cases, it is calculated using an interpolation method based on regional-scale meteorological data from field stations. Accurate values in the field are difficult to obtain without ground measurements. In this paper, a satellite-based method to estimate sunshine duration is introduced and applied over the Heihe River Basin. This method is based on hourly cloud classification product data from the FY-2D geostationary meteorological satellite (FY-2D). A new index—FY-2D cloud type sunshine factor—is proposed, and the Shuffled Complex Evolution Algorithm (SCE-UA) was used to calibrate sunshine factors from different coverage types based on ground measurement data from the Heihe River Basin in 2007. The estimated sunshine duration from the proposed new algorithm was validated with ground observation data for 12 months in 2008, and the spatial distribution was compared with the results of an interpolation method over the Heihe River Basin. The study demonstrates that geostationary satellite data can be used to successfully estimate sunshine duration. Potential applications include climate research, energy balance studies, and global estimations of evapotranspiration.

Aerodynamic roughness length is an important parameter for surface fluxes estimates. This paper developed an innovative method for estimation of aerodynamic roughness length (z_0m) over farmland with a new vegetation index, the Hot-darkspot Vegetation Index (HDVI). To obtain this new index, the normalized-difference hot-darkspot index (NDHD) is introduced using a semi-empirical, kernel-driven bidirectional reflectance model with multi-temporal Proba-V 300-m top-of-canopy (TOC) reflectance products. A linear relationship between HDVI and z_0m was found during the crop growth period. Wind profiles data from two field automatic weather station (AWS) were used to calibrate the model: one site is in Guantao County in Hai Basin, in which double-cropping systems and crop rotations with summer maize and winter wheat are implemented; the other is in the middle reach of the Heihe River Basin from the Heihe Watershed Allied Telemetry Experimental Research (HiWATER) project, with the main crop of spring maize. The iterative algorithm based on Monin–Obukhov similarity theory is employed to calculate the field z_0m from time series. Results show that the relationship between HDVI and z_0m is more pronounced than that between NDVI and z_0m for spring maize at Yingke site, with an R² value that improved from 0.636 to 0.772. At Guantao site, HDVI also exhibits better performance than NDVI, with R² increasing from 0.630 to 0.793 for summer maize and from 0.764 to 0.790 for winter wheat. HDVI can capture the impacts of crop residue on z_0m, whereas NDVI cannot.