Next Article in Journal
Hydraulic Conductivity Characteristics of a Clayey Soil Incorporating Recycled Rubber and Glass Granules
Previous Article in Journal
Basin-Scale Geochemical Assessment of Water Quality in the Ganges River during the Dry Season
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessing and Comparing Reference Evapotranspiration across Different Climatic Regions of China Using Reanalysis Products

1
Key Laboratory of Agricultural Soil and Water Engineering in Arid and Semiarid Areas, Ministry of Education, Northwest A&F University, Yangling District, Xianyang 712100, China
2
College of Water Resources and Architectural Engineering, Northwest A&F University, Yangling District, Xianyang 712100, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Water 2023, 15(11), 2027; https://doi.org/10.3390/w15112027
Submission received: 13 April 2023 / Revised: 20 May 2023 / Accepted: 22 May 2023 / Published: 26 May 2023

Abstract

:
This study aims to assess the accuracy of the reference evapotranspiration (ET0) estimated by CLDAS, ERA5 reanalysis products, and the quality of reanalysis weather variables required to calculate PM-ET0. For this purpose, the applicability of surface meteorological elements from the ERA5 reanalysis datasets provided by the European Centre for Medium-Range Weather Forecasts (ECMWF), and the second-generation China Meteorological Administration Land Data Assimilation System (CLDASV2.0) datasets are evaluated in China by comparison with local observations from 689 stations reported by the Chinese Meteorological Administration (CMA). Statistics including percent bias (PBias), coefficient of determination (R2), root mean square error (RMSE) and mean absolute error (MAE) are used to check the accuracy. The results show the highest correlation between reanalysis temperature and observations, with a mean R2 of 0.96, 0.90 for the CLDAS maximum and minimum air temperatures, and 0.87, 0.84 for ERA5. For the reanalysis of solar radiation (Rs) and relative humidity (RH), an overestimation trend is shown for Rs and an underestimation trend is shown for RH. For reanalysis of wind speed, a relatively low accuracy is shown. The accuracy of ET0 estimated by the two reanalysis products is acceptable in China, but the spatial and temporal consistency between the CLDAS estimates and site observations is higher, with a mean RMSE R2 of 0.91, 0.82 for CLDAS and 1.42, 0.70 for ERA5, respectively. Moreover, CLDAS reanalysis products are more effective in describing the boundary details of the study area.

1. Introduction

Evapotranspiration (ET) is an important meteorological element in the hydrological cycle, which is directly related to the energy and water balance of the Earth’s surface, determining the formation and evolution of the geographical environment [1], and also an important basis and a key link in the evaluation of agricultural water use efficiency [2,3]. As the global competition for water resources intensifies and water resources become depleted, coping with water scarcity requires a more accurate knowledge of ET [4].
Irrigation management in agricultural practices requires the accurate estimation of crop water consumption, which in turn requires an accurate estimation of crop evapotranspiration (ETc), and its forecasting is significant in developing crop irrigation systems and real-time irrigation scheduling [5,6]. A commonly used ETc estimation method at the field scale consists of using the Kc-ET0 approach, as proposed in FAO56 [7], where a crop coefficient (Kc) for the considered vegetation is multiplied by the crop reference evapotranspiration (ET0) to estimate the ETc, although the calculation of Kc is based on experience, and applications worldwide are successful [8,9,10].
At present, ET0 can also be estimated from pan evaporation, but such methods are expensive, technically complex, and not practical. Therefore, the Penman–Monteith model is usually used for the estimation ET0 [11,12,13,14]. However, the model requires the input of several climate variables, including the maximum and minimum air temperature (Tmax and Tmin), relative humidity (RH), solar radiation (Rs), and wind speed at 2 m height (U2). These variables are not observed, or the price of obtaining these meteorological data from the relevant meteorological services is prohibitive in many regions and locations in China. In addition, the common method of obtaining these meteorological data is station monitoring, which is time-consuming and makes it difficult to obtain long time series and high-quality meteorological variables.
Since the 1960s, the application of meteorological satellite data made up for the shortage of station observations, and the spatial and temporal continuity of meteorological data made great progress. After the 1990s, satellite data entered the data assimilation system, which further improved the accuracy of meteorological data [15,16]. Reanalysis information assimilates data obtained from remote sensing, ground-based observations and numerical simulations [17,18], which has the characteristics of high spatial and temporal resolutions and is an ideal driving data source for distributed models. The main advantage of reanalysis products is the free access to continuous meteorological data and its great advantage and potential to become an alternative to observed data.
Several global reanalysis datasets are available for ET0 estimation, such as WorldClim, which provides the global historical and future climate data and elevation data [19]; The World Meteorological Organization’s Climate Explorer weather data retrieval platform, with a wealth of global or regional climate data [20]; the European Centre for Medium-range Weather Forecasts (ECMWF), which provides ERA5 reanalysis products [21]; the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR), known as NCEP/NCAR Reanalysis I [22]; the National Aeronautics and Space Administration (NASA) Global Simulation; and the Assimilation Office for a Second Review of Modern Research and Applications (MERRA-2) [23,24], etc.
The main advantages of the reanalysis product are the high spatial and temporal resolutions, the ability to scale the product appropriately, and the stable improved model resolution and bias [25]. Currently, many authors have used reanalysis products to represent the spatiotemporal variability of surface climate variables. ERA5_land reanalysis data were applied to the observations of Peruvian glaciers verified by Martin et al., and were found to be appropriate for characterizing 2 m air temperature and relative humidity, particularly in the wet outer tropics [26]. Song et al. (2020) [27] compared the applicability of various soil moisture data in Inner Mongolia, and found that the ERA5 simulation capability is optimal. Srivastava et al. (2016) [28] compared the ET0 estimation from ERA-interim and NCEP for all the seasons, and revealed a similar performance to the seasonal assessment, with a higher agreement for the ERA-interim. The only drawback of using a downscaling ERA-interim reanalysis dataset with the WRF is that it is time-consuming. Martins et al. (2016) [29] used the ERA-interim and NOAA NCEP reanalysis datasets for the calculation of PM-ET0 in the Iberian Peninsula and compared it with the ET0 Observed at 130 meteorological stations in the Iberian Peninsula. The results supported the quality of the ERA-interim reanalysis data to calculate the ET0 computations. The ECMWF ERA5 reanalysis products were selected for this study because they had previously performed well and showed great a potential for estimating ET0. However, due to the bias reported for ERA products, the use of these data requires assessing the impact of bias on the performance of the estimations of reanalysis-based ET0 when compared with ET0 computed with observational data [30,31].
CLDAS-V2.0 is the latest version of the land surface assimilation system developed by the China National Meteorological Information Center, providing a variety of atmospheric driving field products (URL: https://data.cma.cn/, accessed on 1 May 2021). At the same time, some scholars have evaluated the applicability of different land surface model data. Liu et al. (2021) [32] compared CLDAS air temperature data with hourly air temperature data from 48,708 ground automatic meteorological stations in China for 2017–2018, and found that CLDAS air temperature data better reflects the interannual variability of air temperature in Chinese regions, with an average correlation coefficient of more than 0.991. Shi et al. (2018) [33] analyzed the temporal variation characteristics of average soil moisture in China by comparing CLDAS soil moisture with hourly observations of automatic soil moisture observation stations, concluding that the simulated values were very close to the observed values, and the time-series variation of soil moisture was captured in each study area. Using GLDAS and CLDAS forcing data, Yang et al. (2017) [34] evaluated precipitation and shortwave radiation from the ITPCAS, GLDAS, CLDAS and the gridded analysis based on CMA gauge observations (CN05.1) over mainland China during 2008–2014, and showed that CLDAS provides more realistic precipitation and radiation approximations than the other datasets. The CLDAS forcing dataset contains more in situ and remote sensing observations of China, so it is more accurate over mainland China. Liu et al. (2019) [35] compared and evaluated the simulation results of monthly soil moisture and monthly evapotranspiration, and reported that CLDAS_Noah-MP significantly improved the simulation results for mainland China and eight river basins; the R-value increased from 0.451 to 0.534 and the RMSE reduced from 0.078 to 0.068, which confirmed the superiority of CLDAS, as previously reported by Cui et al. (2018) [36] and Shi et al. (2018) [37].
The most recent studies have focused on assessing one or several meteorological variables for the reanalysis of atmospheric field-driven products, but relatively few have conducted comprehensive assessments of the meteorological variables required to calculate ET0 and reanalysis for ET0 estimation, especially for small temporal scales. In addition, the numerical–physical approach, horizontal and vertical resolution, and temporal variability of the reanalysis model will bring uncertainties in the reanalysis simulation process [38,39]. The resolution and accuracy of these reanalysis products appear to vary with the location and timescale of the study [40]. Although research indicates that the spatial resolution and reliability of the new generation of reanalysis data products have improved, there is still a need to evaluate the region-specific applicability of CLDASV2.0, ERA5 reanalysis weather variables and the estimation of ET0 in different climate zones of China, and find the best method for researchers and users to carry out irrigation practices.
The main objective of this research is to apply the ET0 calculated from reanalysis data to irrigation scheduling and management, for which it can be used in those places where the weather data quality is low or (and) weather variables are missing [41,42]. The correlations and bias of the five reanalyzed meteorological variables from CLDASV2.0 and ERA5 are compared based on the daily weather variable data obtained from 689 ground automatic meteorological stations in China from 2017–2019. Any biases in these products can be identified by comparison so that the bias correction method can be easily applied. Moreover, a comprehensive set of the statistical indicators is used to assess how calculated reanalysis PM-ET0 compared with the corresponding station observations, which also offers a reference for the selection and application of CLDASV2.0 and ERA5 reanalysis meteorological data across different climatic regions of China.

2. Materials and Methods

2.1. Meteorological Data

CLDAS uses fusion and assimilation technology to fuse data from various sources, such as ground observation, satellite observation, numerical model products and numerical model products, then output the land surface-driven products with high spatial and temporal resolution, including the maximum and minimum air temperatures, 2 m specific humidity, 10 m wind speed, etc. It covers the Asian region (60°–160° E, 0°–65° N), with a spatial resolution of 0.0625° × 0.0625° and a temporal resolution of 1 h. The dataset is formed by using the ECMWF numerical analysis/forecasting product as the background field and the topographic adjustment, multigrid variational assimilation, optimal interpolation and other techniques in the Chinese region, to fuse the ground-based automatic station observations and interpolate them to the reanalysis grid points [33,43].
ERA5 is the fifth generation of ECMWF atmospheric reanalysis of the global climate, created by the EU-funded Copernicus Climate Change Service (C3S). It assimilates remote sensing information, upper atmosphere, and near-surface conventional observation data, including different regions and sources, on a global scale, spanning from 1979 to the present, achieving real-time updates with a spatial resolution of 0.25° × 0.25° latitude-longitude [44]. ERA5 reanalysis data provide many kinds of meteorological elements, including 2 m air temperature, 2 m dew point temperature, and 10 m wind speed, etc. The data sources can be found at http://cdds.climate.copernicus.eu/.

2.2. FAO-56 Penman–Monteith Equation

The reference crop evapotranspiration was calculated using the FAO-56 Penman–Monteith equation recommended by the Food and Agriculture Organization of the United Nations (FAO) in 1998 and the PM-ET0 Equation (1). It is not an empirical equation, but a mechanistic combination equation derived from the definition of the Penman–Monteith equation, parameterized for the reference crops [45]. ET0 (mm day−1) is described as follows.
E T 0 = 0.408 Δ R n G + γ 900 T a + 273 u 2 e s e a Δ + γ 1 + 0.34 u 2
where E T 0 is the reference crop evapotranspiration (mm d−1), R n is the net radiation at the crop surface (MJ m−2 day−1) and G is the soil heat flux density (MJ m−2 day−1). T a is the mean daily air temperature at 2 m height (°C) calculated from the maximum and minimum air temperature (Tmax and Tmin, °C), u 2 is the wind speed at 2 m height (m s−1), e s is the saturated vapor pressure (kPa), e a is actual water vapor pressure (kPa), e s e a is the vapor pressure deficit (kPa) calculated from the saturation vapor pressure and actual vapor pressure, Δ is the slope of the vapor pressure curve (kPa °C−1) and γ is the psychrometric constant (kPa °C−1). In the daily calculation scale, G is small relative to R n , and is therefore neglected [46,47].

2.3. Data Sources

Station observation data were collected from 689 meteorological stations of the China Meteorological Administration (https://data.cma.cn/). The list of variables retrieved is defined by the Penman–Monteith model FAO56 framework, including the maximum and minimum air temperatures at 2 m, surface radiation, sunshine duration, relative humidity at 2 m and wind speed at 10 m. The data are subjected to quality control such as climate threshold value check, temporal consistency, and spatial consistency check.
ET0 is significantly affected by seasonal variation and regional differences. The two types of data from the different time series and different study regions are compared in the study, and the detailed analysis and comparison are presented in Section 3. Tmax and Tmin, RH, Rs, U2 and the estimated PM-ET0 from CLDAS2.0, ERA5, and station observations in 2017–2019 in the Chinese region are selected as the research objects. China was divided into 7 climatic zones according to climatic characteristics and heat indexes [48,49]. Regions 1–7 are the Northwest desert zone, Inner Mongolia grassland zone, Northeast humid and semi-temperate zone, humid and semi-humid warm temperate zone, humid subtropical zone, humid tropical zone and the Qinghai Tibet Plateau zone (Figure 1). Bias, correlation, and error of the estimated PM-ET0 across different climate zones and the CLDAS and ERA5 reanalysis weather variables used to estimate ET0 are analyzed.

2.4. Spatial Interpolation Method

The meteorological variables related to the calculation of PM-ET0 have been extracted from the CLDAS2.0, ERA5 reanalysis grid points, and then these grid point meteorological data were interpolated to the same latitude and longitude positions as the 689 ground observation sites, and the meteorological data obtained from the interpolation are evaluated as true values for the test. Grid data from four grid points around it were selected and interpolated to the station by the inverse distance weight (IDW) method. The formula is as follows:
Z = i = 1 n Z i D i i = 1 n 1 D i
where Z is the true value, Z i is the value of the control point, and D i is the weight coefficient.

2.5. Accurate Statistical Indicators

The assessment focused on the pairwise comparison between the weather stations’ observed values ( Q i ) and the corresponding values ( P i ) for the two reanalysis datasets for each meteorological variable used in ET0 computations and for the ET0 values [50]. The following several statistical indicators were used to assess the performance of the reanalysis datasets in representing the spatial and temporal variation of the meteorological variables.
The percent bias (PBias) measures the degree of deviation between the reanalysis estimated data and the corresponding site observations [51]. It is simply a normalized difference between the means of both sets of values: positive values indicate overestimation, and negative values indicate underestimation.
PBias = i = 1 n P i Q i i = 1 n Q i
The coefficient of determination (R2) is used to assess the degree of dispersion between the reanalysis estimated value and the site observations. The closer the R2 is to 1, the closer the actual site measurements are to the estimated and the higher the correlation, calculated as:
R 2 = i = 1 n Q i Q i ¯ P i P i ¯ 2 i = 1 n Q i Q i ¯ 2 i = 1 n P i P i ¯ 2
The root means square error (RMSE) measures the overall difference between the reanalysis estimated values and the site observations. Because of its high sensitivity to outliers, RMSE can better reflect the accuracy of the data.
RMSE = 1 n i = 1 n Q i P i 2
The mean absolute error (MAE) is a measure of the error between the paired reanalysis estimated values and site observations, computed as:
MAE = 1 n i = 1 n Q i P i 2

3. Results

3.1. Assessing the Accuracy of the Reanalysis Weather Variables

The assessment indicators PBias, RMSE, MAE and R2 for the reanalysis estimated air temperatures (Tmin CLDAS, Tmax CLDAS, Tmin ERA5, and Tmax ERA5) and the corresponding observed values are summarized in Table 1 and Table 2. Reanalysis Tmin shows a trend of overestimation, with the average PBias of Tmin CLDAS and Tmin ERA5 being 0.07 and 0.06, respectively. On the contrary, Tmax reveals a trend of underestimation, and the mean PBias of Tmax CLDAS and Tmax ERA5 are −0.06 and −0.34, respectively. When comparing the reanalysis Tmax, the correlation between reanalysis Tmin and site observations is higher, the average R2 of the two reanalysis Tmin being 0.96, 0.87 and the corresponding R2 values for reanalysis Tmax 0.90, 0.84, but with a larger PBias in the case of Tmin.
Tmax and Tmin are the key meteorological factors to calculate ET0. Therefore, RMSE, MAE and R2 are further used to assess the error and correlation of the reanalysis air temperature. It can be seen from Table 1 and Table 2 that the air temperature estimates of CLDAS and ERA5 products are consistent with the site observations, and the correlation values are between 0.76 and 0.98. Among the two products, CLDAS estimations are much closer to the observed values, with lower root means square error values (average range from 2.45 to 3.52) and higher correlation values (average range from 0.96 and 0.90) in most climate zones compared to those obtained from ERA5 (average range from 9.32 to 9.66 and 0.87 to 0.84, respectively).
Four indicators of statistical accuracy of reanalysis solar radiation (Rs) across seven different climate zones and mainland China are shown in Table 3. The overall trends of the PBias, RMSE, MAE and R2 for reanalysis Rs CLDAS and Rs ERA5 are consistent, with differences in the numerical values. The mean R2 of Rs CLDAS is between 0.50 and 0.65, and the mean R2 of Rs ERA5 is between 0.40 and 0.55, indicating that the reanalysis Rs agrees well with Rs Obs. However, reanalysis Rs CLDAS better reflects the variability of the site-observed Rs in the overwhelming of locations compared to Rs ERA5. The bias of Rs CLDAS and Rs ERA5 ranges from 0.02 to 0.17 and 0.06 to 0.24, demonstrating that both the reanalysis products slightly overestimate the observed Rs.
The spatial distribution of PBias, RMSE, MAE and R2 of Rs across the different climatic regions of China are presented in Figure 2. The regional distribution of the statistical indicators of Rs is characterized clearly. The estimation bias and error are generally lower in the northwest and northeast regions than in the southeast; in other words, the mean PBias, RMSE and MAE are small in arid and semi-arid regions (climate zones 1, 2) and cold regions (climate zone 3), with distributions ranging from 0.02 to 0.13, 4.04 to 4.43 and 2.87 to 3.51, respectively. The average biases and estimation errors are large for humid and semi-humid areas (climate zones 4,5 and climate zone 6), with distributions of the PBias, RMSE and MAE ranging from 0.10 to 0.24, 4.33 to 5.23 and 3.16 to 4.01, respectively. The R2 of Rs is the opposite of the previous three indicators. The R2 of Rs CLDAS and Rs ERA5 is between 0.65 and 0.58, and 0.55 and 0.50 in climatic zones 1, 2 and 3, while the R2 is between 0.51 and 0.54 and 0.38 and 0.46 in climatic zones 4, 5 and 6. As a whole, the reanalysis of Rs CLDAS and Rs ERA5 reflected the temporal and spatial variation characteristics of Rs Observed well at most stations.
The statistical indicators of reanalysis relative humidity are shown in Table 4. The PBias of RHCLDAS is negative (ranging from −0.07 to −0.30) across the 1–6 climate zones. The results indicate a slight underestimation, with the PBias of climate zone 7 being positive, with a value of 0.35. The PBias of RHERA5 are 0.04 and 0.20 in climate zones 4 and 7, respectively, indicating a slight overestimation trend, and the remaining climate zones have a negative PBias, with values from −0.17 to −0.34. As for the error, compared with other climate zones, the estimation errors for climate zone 7 are larger, the RMSE and MAE for RHCLDAS are 25.94% and 22.51%, respectively, and those of ERA5 are 22.70% and 18.84%. Reanalysis RH and the observed RHObs are in better agreement in climate zone 2, and the correlation values are 0.68 and 0.53 for RHCLDAS and RHERA5, but the estimation error is relatively high, with the corresponding RMSE values of 11.63% and 18.51%. Furthermore, the estimation errors for climate zones 1–3 (RMSE between 10.82 and 21.75%) are generally larger than those for climate zones 4–6 (RMSE between 7.91 and 17.47%), partly because the latter are located in the humid zones with little variation in daily relative humidity, while the former are in the arid desert and grassland zones with a low vegetation cover and more frequent changes in daily relative humidity.
In contrast to the performance of the previous weather variables, both U2 CLDAS and U2 ERA5 show a relatively low accuracy, with a high dispersion of pairs of 2 m wind speed derived from the reanalysis datasets and the station observations along the regression line, with R2 averages of 0.25 for U2 CLDAS and 0.18 for U2 ERA5. For PBias, U2 CLDAS is underestimated relative to RHObs., and PBias ranged from −0.31 to −0.63 in seven different climatic regions. U2 ERA5 is underestimated in climate zones 1–3, with PBias ranging from −0.33 to −0.78, and overestimated in the rest of the climate zones with PBias ranging from 0.21 to 0.41. In terms of RMSE and MAE, the mean values of the RMSE and MAE are 1.09, 0.84 for U2 CLDAS and 1.23, 0.99 for U2 ERA5 at the national scale, respectively, the details can be found in Table 5.

3.2. Comparison of the Annual Mean of the Observation-Calculated ET0 and Reanalysis-Estimated ET0

A spatial distribution of the annual mean values of ET0 calculated from the stations observations (ET0 Obs), ET0 estimated by CLDAS (ET0 CLDAS), and ET0 estimated by ERA5(ET0 ERA5) in China, is shown in Figure 3. The frequency distribution of the average annual evapotranspiration between different value ranges is shown in the histogram of Figure 3. The spatial distribution of the annual mean values of ET0 estimated from the two reanalysis datasets are approximately the same; however, from the perspective of frequency distribution, ET0 CLDAS is more consistent with ET0 Obs. For the observed calculated ET0, only 15.14% of stations have an ET0 Obs < 500 mm, 23.41% of stations an ET0 Obs between 650 and 800 mm, 16.15% of stations an ET0 Obs between 950 and 1000 mm, and 3.34% of stations an ET0 Obs > 1000 mm. For the CLDAS estimated ET0, 13.10% of the stations showed ET0 CLDAS <500 mm, ET0 CLDAS ranged from 650 mm to 800 mm for 25.21% of the stations, from 800 mm to 900 mm for 15.97% of stations, and 5.82% of stations exhibited an ET0 CLDAS >1000 mm. The frequency distribution of different intervals of the annual mean evapotranspiration estimated by CLDAS is similar to that of the site-observed ET0, but in most cases, ET0 CLDAS is slightly overestimated.
For ERA5 estimated ET0, 11.79% of stations have an ET0 ERA5 between 500 and 650 mm, which is 3.64% higher than the frequency of ET0 Obs observations in this range. Of the stations, 21.69% have ET0 ERA5 between 950 and 1000 mm, an increase of 5.53% over the percentage of ET0 Obs in this interval, while the proportion of ET0 ERA5 with an annual average evapotranspiration interval in the range of 800–950 mm is lower than that of ET0 Obs. The frequency distribution of ET0 ERA5 in different evapotranspiration intervals differs from that of ET0 Obs observations, and ET0 ERA5 is overall overestimated relative to the ET0 Obs site-observed values.
The statistical indicators of ET0 from station observations and two reanalysis estimates show significant regional variations in different climatic regions of China (Figure 4). The annual mean ET0 in climate zones 1 and 2 is mainly concentrated between 800 and 950 mm, with a higher ET0 in Xinjiang; for example, ET0 is greater than 1000 mm in Urumqi and other stations. The annual mean ET0 in climate zone 3 is low, with the lower values appearing in the eastern and northern parts of climate zone 3, and the lowest values of ET0 in the whole of China also occur in this region, which may be caused by the lower average temperature. The majority of the stations in climate zone 4 have a small trend of annual mean ET0 fluctuations and are concentrated in the range of 800–950 mm; however, ET0 changes in the west of climate zone 4 and the southwest of climate zone 5, and its accuracy statistical indicators also change abruptly at this boundary location, and the change in ET0 CLDAS is more significant, mainly because these two locations are susceptible to topography. The high fluctuation of ET0 in climatic zones 6 and 7, where the higher PBias and RMSE values also occur (Figure 4), is due to the influence of the topography of the Tibetan Plateau region. In addition, the short sunshine time and high relative humidity cause a lower ET0 in the southwestern part of climate zone 6 and the eastern part of climate zone 7. Compared with the ERA5 reanalysis data, CLDAS can better simulate the detailed characteristics of ET0 in different climate zones.
To explore the accuracy of estimating the daily ET0 from CLDAS, ERA5 reanalysis data for different climatic zones is shown. The paired comparisons between ET0 Obs and ET0 CLDAS (1st line, 3rd line), and ET0 Obs and ET0 ERA5 are shown in scatter plots in Figure 5, in which all stations are from seven different climate zones. The analysis in Figure 5 complements the one in Figure 4 in detail. The results show that except for climate zone 7, the simulated ET0 of the reanalysis datasets is closer to the ET0 Obs of the stations in the arid and semi-arid areas; that is, the scatter points are concentrated on the y = x line, and the lowest value of the ET0 CLDAS correlation value is 0.91, and the lowest value of the ET0 ERA5 correlation value is 0.85. The R2 values of the two reanalysis products are smaller in humid, and sub-humid areas, ranging from 0.71 to 0.86, and 0.63 to 0.72 for ET0 CLDAS and ET0 ERA5. The results indicate that the variability of meteorological variables is higher in humid and semi-humid areas than in arid and semi-arid areas. Moreover, the accuracy of the reanalysis products for extreme weather estimation is low, and the extreme values of ET0 Obs tend to be reconciled by the reanalysis datasets. For example, when the ET0 Obs value at site 50,557 is greater than 6 mm d−1, the reanalysis estimates ET0 to be less than the corresponding site-observed ET0 value. When the ET0 Obs value at the sites is 59, 133 is less than 2 mm d−1, and the ET0 simulated by ERA5 is greater than the ET0 Observed at the corresponding site. These phenomena can be explained in the reanalysis temperature variables in Section 2.2.

3.3. Assessing the Accuracy of ET0 CLDAS and ET0 ERA5 Estimates Using Reanalysis Products

The statistical accuracy indicators for the reanalysis estimation of ET0 and observed ET0 in the seven climate zones are summarized in Table 6. The results show that the performance of CLDAS in estimating ET0 is better than that of ERA5 reanalysis products: the average RMSE = 0.91 mm d−1 is smaller than that for ERA5 (RMSE = 1.42 mm d−1), while the average determination coefficient is larger for the CLDAS reanalysis products (R2 = 0.82) than for ERA5 (R2 = 0.70). Similar results have been shown by Liu et al. (2019) [35], who reported that the simulation results of CLDAS in most climate zones maintained a consistent and reasonable performance. The PBias distribution of ET0 estimated for two reanalysis products shows that the number of positive bias stations is more than the number of negative bias stations. The ET0 CLDAS shows a negative bias in climate zones 1,6 and climate zone 7, with a PBias of −0.05, −0.03 and −0.14, respectively; ET0 Obs is underestimated in these climate zones. ET0 CLDAS shows a positive bias in all other climate zones, with a PBias distribution range of 0.01–0.03; ET0 Obs is overestimated, but to a lesser extent. ET0 ERA5 shows negative bias in climate zones 1,2 and climate zone 7, with PBias ranging from −0.08 to −0.18, and the ET0 Obs observations are underestimated, the remaining climate zones have a positive PBias and the ET0 Obs observations are overestimated.
The results of the ET0 estimation generated from the two reanalysis datasets revealed that the reanalysis product presented a superior performance in the arid and semi-arid regions, with a higher correlation at most sites, with R2 > 0.85. Like the Northwest desert zone, the Inner Mongolia grassland zone and R2 of ET0 CLDAS are 0.91, 0.90, and the R2 of ET0 ERA5 are 0.87, 0.85, respectively, while the correlations were slightly worse in humid and semi-humid areas such as climate zones 5 and 6. The R2 distributions of ET0 CLDAS range from 0.69 to 0.77 and 0.54 to 0.62 for ET0 ERA5 (Table 6).To further evaluate the accuracy of the CLDAS, the ERA5 reanalysis datasets in estimating ET0, RMSE and MAE for different climate zones are compared. ET0 CLDAS and ET0 ERA5 have the largest estimation errors in climate zone 7, with an RMSE and MAE of 1.22, 1.06 and 1.67, 1.52, respectively. Except for climate zone 7, the RMSE and MAE of ET0 CLDAS range from 0.78 to 0.91 and 0.52 to 0.68, respectively, the distribution of the RMSE for ET0 ERA5 ranged from 1.12 to 1.56 and MAE from 1.04 to 1.39. Moreover, the higher values of the RMSE and MAE for ET0 are mainly concentrated in coastal humid areas such as climate zones 4, 5 and 6 (Figure 6), with an RMSE > 0.88 and MAE > 0.65 for ET0 CLDAS, and ET0 ERA5 with an RMSE > 1.42 and MAE > 1.28. Overall, the RMSE and MAE are shorter and closer to the optimal target when estimating ET0 using the CLDAS reanalysis products. The effects of topography, slope, monsoon and other factors on meteorological variables in different climatic zones are adopted and integrated into the CLDAS reanalysis datasets, which makes the estimation of ET0 CLDAS more accurate [35,52].

3.4. Comparison of Time Series Changes in the Monthly Mean Values of the Observed ET0 and Reanalysis of Estimated ET0

The monthly time series variation of ET0 for the seven climate zones across the country in 2017–2019 is shown in Figure 7. The two reanalysis estimates of ET0 agree well with the observed ET0 values in most cases, indicating that the reanalysis products capture the seasonal cycle and temporal evolution based on the observed ET0 to some extent. However, ET0 ERA5 overestimates the fluctuation amplitudes and the estimated ET0 values are all much higher than the observed values in climate zones 4 and 6. ET0 ERA5 underestimates the fluctuation amplitudes and the estimated ET0 values are all much lower than the observed values in climate zone 7. In contrast, the ET0 CLDAS estimates are closer to the observed values in most climate zones.
The calculation of ET0 is affected by meteorological variables such as air temperature, relative humidity and solar radiation, and these meteorological factors vary with time and region. ET0 also has obvious seasonal variation and spatial distribution characteristics. The monthly mean ET0 is large in summer and autumn in different climatic zones, and statistics found that the cumulative ET0 from May to October each year accounts for about 75–80% of the whole year. ET0 grows with increased precipitation and temperature during the summer season. Due to the influence of geographical location and meteorological factors, the peak of monthly mean ET0 occurs at different times, with climate zones 1 and 2 showing a single peak pattern of monthly mean ET0, with the peak occurring in July each year and the trough occurring in November, lasting until April the following year. There are differences between the annual intervals of monthly mean ET0 in climate zones 3 and 6, where 2017 showed a single peak type, with a peak in June, and 2018–2019 a double peak type, with a peak in May and a low peak in July. The monthly mean ET0 in climate zones 4, 5 and 7 show a double peak pattern, with peaks in April and July, and a valley from October to February of the next year.

4. Discussion

4.1. Analysis of Meteorological Variables Related to ET0 Estimation

In this paper, the key meteorological factors for calculating ET0 are evaluated. It was found that Tmin shows an overestimation trend and Tmax shows an underestimation trend, and the consistency of the reanalysis minimum temperature of the two datasets (average R2 0.96 and 0.87) is higher than that of the maximum temperature (0.90 and 0.84), but the PBias of the former is larger. Similar conclusions were found by Wu et al. (2022) [53] and Qian et al. (2022) [54] when they used reanalysis data to assess air temperature in the Chinese region. In addition, Simmons et al. (2010) [55] and Paredes et al. [56] also reported a tendency for the overestimation of Tmin and underestimation of Tmax by the ERA-Interim in Continental Portugal [57]. The authors hypothesized that CLDAS, ERA5 reanalysis datasets are capturing warming over the land more than over the sea, which makes Tmin CLDAS and Tmin ERA5 larger than Tmin Obs data; Tmax CLDAS, Tmax ERA5 showed the opposite performance. For the reanalysis Rs, both CLDAS and ERA5 slightly overestimate the site-observed Rs, the PBias of Rs CLDAS and Rs ERA5 range from 0.02 to 0.17 and 0.06 to 0.24, which is similar to the reports of Bojanowski et al. (2014) [57] and Urraca et al. (2017) [58]. Sheffield et al. (2006) [25] also reported a similar conclusion and combined this behavior of overestimation of global radiation for NCEP1 with the fact that the reanalysis dataset does not consider the blocking of shortwave radiation by the clouds. As for reanalysis of the relative humidity, the PBias of RHCLDAS is negative (ranging from −0.07 to −0.30) in the 1–6 climate zones, with a slight underestimation. The PBias of climate zone 7 is positive, and the value is 0.35. This is consistent with the conclusion of Simmons et al. (2010) [55] that the relative humidity of ERA5 reanalysis is underestimated in the mid–low latitudes. Fu et al. (2015) [59] also concluded that the RH of NCEP1 reanalysis was slightly underestimated in Australia. Differently, Martins et al. (2016) [29] reported no tendency for over- or underestimation of the RH when using the reanalysis monthly products for Iberia, which can be attributed to the different time scales of the study. Different from the previous four meteorological factors, the analysis of the wind speed showed a lower accuracy across the seven climatic regions. The average R2 of U2 CLDAS was 0.25, and the average R2 of U2 ERA5 0.18. Many scholars have reanalyzed the wind speed and obtained similar results, the low accuracy of the wind speed derived from reanalysis was reported by Carvalho et al. (2014) [60] studies applied to Continental Portugal. The same results were reported by Martins et al. (2016) [29] in the Iberian Peninsula for the problem of low accuracy in analyzing the wind speed. The two wind speed estimation alternatives proposed by Allen et al. (1998) [7] are to replace the daily wind speed with the annual average wind speed over the local area or use the world average wind speed value U2 = 2 m s−1 as the default wind speed, obtained from over 2000 weather stations worldwide. Jabloun and Sahli (2008) [61] compared the above methods and found that both methods obtained wind speeds with high applicability, but the ET0 was better estimated using the local area’s annual average wind speed.

4.2. Analysis of the Spatial Distribution of the Observation-Calculated ET0 and Reanalysis-Estimated ET0

The spatial and frequency distribution of the annual mean values of ET0 calculated from the station observations, ET0 estimated by CLDAS and ET0 estimated by ERA5 in China were analyzed. It was found that the frequency distribution and spatial distribution of different intervals of annual mean evapotranspiration estimated by CLDAS and ERA5 are similar to the observed ET0, but in most cases, the reanalysis ET0 is slightly overestimated. The overestimation of the reference evapotranspiration derived from reanalysis has been discovered in many relevant studies. Betts et al. (2009) [62] compared ET0 at the watershed scale using ERA-Interim reanalysis data and observations, and found that the reference crop evapotranspiration from reanalysis data was significantly larger than the observed value. Moreover, other studies using ERA-Interim reanalysis products also reported an overall overestimation bias for ET0 [28]. Dee and Uppala (2009, 2011) [44,63] considered that the overestimation of reanalysis data is mainly related to the utilization of small timescale (daily or smaller timescale) reanalysis products, water vapor movement processes and variational bias correction.

4.3. Analysis of the Accuracy of ET0 Estimated by CLDAS and ERA5 across Different Climatic Regions

Comparing the accuracy of ET0 estimation across different climatic regions, it was found that the variability of ET0 in humid and semi-humid areas is higher than that in arid and semi-arid areas. The RMSE of ET0 CLDAS ranged from 0.78 to 0.82 in climate zones 1–3, the RMSE ranged from 0.89 to 0.91 in climate zones 4–6, the RMSE of ET0ERA5 ranged from 1.12 to 1.33 in climate zones 1–3 and the RMSE ranged from 1.42 to 1.59 in climate zones 4–6. Su et al. (2015) [64] also found that the reanalysis estimates of ET0 in arid and extremely cold regions are more consistent with the site observations of ET0. The same trend was found by Zhang et al. (2018) [65] when they used ERA-Interim reanalysis data to study the characteristics of ET0 in global arid and semi-arid regions. Furthermore, Paredes et al. (2018) [57], using the ERA-Interim reanalysis product for the daily estimation of grass reference evapotranspiration, also found that the higher values of RMSE were mostly concentrated in coastal areas of Portugal. These findings are consistent with the conclusions of this paper. This can be attributed to the meteorological factors related to the estimation of ET0 in result Section 2.1. We also reported that the accuracy of reanalysis Rs in arid and semi-arid regions is higher than that in humid regions, which is partly because the calculation of solar radiation is a complex task in itself, and the high uncertainty of atmospheric turbidity in humid climatic conditions, which further exacerbates errors between the reanalysis of solar radiation data and station observations, Rs is the most significant meteorological factor affecting ET0 [66]. The high error of the reanalysis Rs leads to a high variability of ET0 in humid and semi-humid climate zones.

5. Summary and Conclusions

It is difficult to obtain long timeseries meteorological variables required for PM-ET0 calculation from ground stations. The reanalysis datasets provide all the weather variables required for calculating ET0, but needs to be evaluated for accuracy. In this paper, meteorological factors such as Tmin and Tmax, Rs, RH and U2 derived from CLDAS and ERA5 are statistically analyzed. Moreover, based on the daily meteorological data of 689 surface meteorological stations covering seven different climatic regions in China, ET0 CLDAS and ET0 ERA5 at different timescales are compared and evaluated.
The quality of meteorological factors for calculating PM-ET0 obtained from CLDAS and ERA5 reanalysis products is acceptable, and the correlation between reanalysis temperature and station observation temperature is high, with R2 > 0.90 for Tmin CLDAS, Tmax CLDAS, and R2 > 0.84 for Tmin ERA5, Tmax ERA5; however, there is a tendency for an underestimation of Tmax, and conversely, a tendency for an overestimation of Tmin. Rs and RH can be accurately estimated from reanalysis datasets, but reanalysis Rs has a slight overestimation trend for site observation Rs and an underestimation trend for site observation RH. Compared with the ERA5 reanalysis products, Tmin, Tmax, Rs and RH derived from the CLDAS products are closer to the site observations, and the overall RMSE and MAE values are smaller. In contrast, the accuracy of both U2 CLDAS and U2 ERA5 is coarser, but the reanalysis of wind speed does not cause a large variability in ET0.
Two reanalysis products capture the seasonal cycle and monthly time evolution based on the observed ET0. During the year, 75–80% of ET0 is concentrated in April–October, but the monthly mean ET0 in different climate zones shows single-peak and double-peak type variations, with different peak occurrence times. In addition, the agreement between the reanalysis estimated ET0 and site-observed ET0 is high in arid and semi-arid areas, with the lowest value of ET0 CLDAS correlation coefficient being 0.90 and ET0 ERA5 0.85. The consistency is slightly worse in wet and semi-humid areas, with the R2 variation range of 0.69–0.88 for ET0 CLDAS and 0.54–0.79 for ET0 ERA5. CLDAS reanalyzed meteorological variables which are used to compute the PM-ET0, and estimated ET0 closer to the site observation ET0 Obs, comparing the corresponding values of ERA5.
These results suggest that reanalysis weather products (Rs, Tmax, Tmin, RH and U2) can be used to estimate ET0 when the observed weather data are unavailable. However, the causes of the impact of errors in each weather variable on the accuracy of the estimated ET0 are not clear, and further studies are required to clarify the quantitative sensitivity of the errors in the reanalysis weather variables on the estimated ET0 as a way to decrease the reanalysis product estimation errors obtained when using daily time steps. Moreover, studies are needed to combine long past series data with forecasts, to achieve the real-time operation of irrigation scheduling models.

Author Contributions

Conceptualization, X.Y. and L.Q.; methodology, X.Y., W.W. and L.Q.; data curation, X.Y. and X.H. (Xuefei Huo); writing—original draft preparation, X.Y., L.Q., X.H. (Xiaotao Hu) and Y.W.; writing—review and editing, X.Y., L.Q., W.W. and X.H. (Xuefei Huo); project administration, X.Y. and L.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (52079113, U2243235), and the National key research and development program (2022YFD1900402-01).

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank the European Centre for Medium-Range Weather Forecasts (ECMWF) and Chinese Meteorological Administration/National Meteorological Information Center (CMA/NMIC) for providing ERA5 and CLDAS reanalysis meteorological products, and the CMA/NMIC for providing surface meteorological station data.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analysis, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Holdridge, L.R. Life Zone Ecology; Tropical Science: San Jose, Costa Rica, 1967. [Google Scholar]
  2. Pereira, L.S.; Perrier, A.; Allen, R.G.; Alves, I. Evapotranspiration: Review of concepts and future trends. J. Irrig. Drain Eng. 1999, 125, 45–51. [Google Scholar] [CrossRef]
  3. Pereira, L.S.; Allen, R.G.; Smith, M.; Raes, D. Crop evapotranspiration estimation with FAO56: Past and future. Agric. Water Manag. 2015, 147, 4–20. [Google Scholar] [CrossRef]
  4. Scanlon, B.R.; Faunt, C.C.; Longuevergne, L.; Reedy, R.C.; Alley, W.M.; McGuire, V.L.; McMahon, P.B. Groundwater depletion and sustainability of irrigation in the US High Plains and Central Valley. Proc. Natl. Acad. Sci. USA 2012, 109, 9320–9325. [Google Scholar] [CrossRef] [PubMed]
  5. Feng, Y.; Peng, Y.; Cui, N.; Gong, D.; Zhang, K. Modeling reference evapotranspiration using extreme learning machine and generalized regression neural network only with temperature data. Comput. Electron. Agric. 2017, 136, 71–78. [Google Scholar] [CrossRef]
  6. Dong, Y.; Zhao, Y.; Zhai, J.; Zhao, J.S.; Han, J.Y.; Wang, Q.M.; He, G.H. Changes in reference evapotranspiration over the non-monsoon region of China during 1961–2017: Relationships with atmospheric circulation and attributions. Int. J. Clim. 2020, 41, E734–E751. [Google Scholar] [CrossRef]
  7. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration: Guidelines for Computing Crop Water Requirements; Food and Agriculture Organization of the United Nations: Rome, Italy, 1998. [Google Scholar]
  8. Valiantzas, J.D. Simplified forms for the standardized FAO-56 Penman–Monteith reference evapotranspiration using limited weather data. J. Hydrol. 2013, 505, 13–23. [Google Scholar] [CrossRef]
  9. Pereira, L.S.; Paredes, P.; Rodrigues, G.C.; Neves, M. Modeling malt barley water use and evapotranspiration partitioning in two contrasting rainfall years. Assessing AquaCrop and SIMDualKc models. Agric. Water Manag. 2015, 159, 239–254. [Google Scholar] [CrossRef]
  10. Jensen, M.E.; Allen, R.G. Evaporation, Evapotranspiration, and Irrigation Water Requirements; American Society of Civil Engineers: Reston, VA, USA, 2016. [Google Scholar] [CrossRef]
  11. Monteith, J.L. Evaporation and the Environment in the State and Movement of Water in Living Organisms. Symp. Soc. Exp. Biol. 1965, 19, 205–234. [Google Scholar]
  12. Popova, Z.; Kercheva, M.; Pereira, L.S. Validation of the FAO methodology for computing ETo with limited data. Application to south Bulgaria. Irrig. Drain. 2006, 55, 201–215. [Google Scholar] [CrossRef]
  13. McMahon, T.A.; Finlayson, B.L.; Peel, M.C. Historical developments of models for estimating evaporation using standard meteorological data. Wiley Interdiscip. Rev. Water 2016, 3, 788–818. [Google Scholar] [CrossRef]
  14. Martínez-Cob, A.; Tejero-Juste, M. A wind-based qualitative calibration of the Hargreaves ET0 estimation equation in semiarid regions. Agric. Water Manag. 2004, 64, 251–264. [Google Scholar] [CrossRef]
  15. Eyer, J.R. Progress Achieved on Assimilation of Satellite Data in Numerical Weather Prediction over the Last 30 Years/Proceeding of ECMWF Seminar on Recent Developments in of Satellite Observations in Numerical Weather Prediction; ECMWF Publication: Reading, UK, 2007; pp. 1–27. [Google Scholar]
  16. Eyre, J.R.; Kelly, G.A.; McNally, A.P.; Andersson, E.; Persson, A. Assimilation of TOVS radiance information through one-dimensional variational analysis. Q. J. R. Meteorol. Soc. 1993, 119, 1427–1463. [Google Scholar] [CrossRef]
  17. Harris, B.A.; Kelly, G. A satellite radiance-bias correction scheme for data assimilation. Q. J. R. Meteorol. Soc. 2001, 127, 1453–1468. [Google Scholar] [CrossRef]
  18. Wang, Y.; Han, W.; Xue, J.; Li, J.; Lei, B.; Yang, J. Adaptive tuning of background error and satellite radiance observation error for operational variational assimilation. In MIPPR 2007: Remote Sensing and GIS Data Processing and Applications; and Innovative Multispectral Technology and Applications; SPIE: Washington, DC, USA, 2007; Volume 6790, pp. 1063–1071. [Google Scholar] [CrossRef]
  19. Fick, S.E.; Hijmans, R.J. WorldClim 2: New 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 2017, 37, 4302–4315. [Google Scholar] [CrossRef]
  20. Trouet, V.; Van Oldenborgh, G.J. KNMI Climate Explorer: A Web-Based Research Tool for High-Resolution Paleoclimatology. Tree-Ring Res. 2013, 69, 3–13. [Google Scholar] [CrossRef]
  21. Balsamo, G.; Albergel, C.; Beljaars, A.; Boussetta, S.; Brun, E.; Cloke, H.; Dee, D.; Dutra, E.; Muñoz-Sabater, J.; Pappenberger, F.; et al. ERA-Interim/Land: A global land surface reanalysis data set. Hydrol. Earth Syst. Sci. 2015, 19, 389–407. [Google Scholar] [CrossRef]
  22. Kanamitsu, M.; Ebisuzaki, W.; Woollen, J.; Yang, S.-K.; Potter, G.L. NCEP–DOE AMIP-II Reanalysis (R-2). Bull. Am. Meteorol. Soc. 2002, 83, 1631–1644. [Google Scholar] [CrossRef]
  23. Kazutoshi, O.; Junichi, T.; Hiroshi, K.; Masami, S.; Shinya, K.; Hiroaki, H.; Takanori, M.; Nobuo, Y.; Hirotaka, K. The JRA-25 Reanalysis. J. Meteorol. Soc. Jpn. 2007, 85, 369–432. [Google Scholar] [CrossRef]
  24. Kobayashi, S.; Ota, Y.; Harada, Y.; Ebita, A.; Moriya, M.; Onoda, H.; Onogi, K.; Kamahori, H.; Kobayashi, C.; Endo, H.; et al. The JRA-55 Reanalysis: General Specifications and Basic Characteristics. J. Meteorol. Soc. Jpn. Ser. II 2015, 93, 5–48. [Google Scholar] [CrossRef]
  25. Sheffield, J.; Goteti, G.; Wood, E.F. Development of a 50-Year High-Resolution Global Dataset of Meteorological Forcings for Land Surface Modeling. J. Clim. 2006, 19, 3088–3111. [Google Scholar] [CrossRef]
  26. Martí, B.; Jose, U.; Giovanni, L.; Philipp, K.; Álvaro, N.; Rolando, C. Validation of ERA5-Land temperature and relative humidity on four Peruvian glaciers using on-glacier observations. J. Mt. Sci. 2022, 19, 1849–1873. [Google Scholar] [CrossRef]
  27. Song, H.Q.; Sun, X.L.; Li, Y.P. Evaluation of ERA5 reanalysis soil moisture over inner mongolia. Sci. Technol. Eng. 2020, 20, 2161–2168. [Google Scholar]
  28. Srivastava, P.K.; Han, D.; Islam, T.; Petropoulos, G.P.; Gupta, M.; Dai, Q. Seasonal evaluation of evapotranspiration fluxes from MODIS satellite and mesoscale model downscaled global reanalysis datasets. Theor. Appl. Clim. 2016, 124, 461–473. [Google Scholar] [CrossRef]
  29. Martins, D.S.; Paredes, P.; Raziei, T.; Pires, C.; Cadima, J.; Pereira, L.S. Assessing reference evapotranspiration estimation from reanalysis weather products. An application to the Iberian Peninsula. Int. J. Clim. 2016, 37, 2378–2397. [Google Scholar] [CrossRef]
  30. Zhao, T.; Guo, W.; Fu, C. Calibrating and Evaluating Reanalysis Surface Temperature Error by Topographic Correction. J. Clim. 2008, 21, 1440–1446. [Google Scholar] [CrossRef]
  31. Hwang, S.; Graham, W.D.; Geurink, J.S.; Adams, A. Hydrologic implications of errors in bias-corrected regional reanalysis data for west central Florida. J. Hydrol. 2014, 510, 513–529. [Google Scholar] [CrossRef]
  32. Liu, Y.; Shi, C.X.; Wang, H.J.; Han, S. Applicability assessment of CLDAS temperature data in China. Atmos. Sci. 2021, 44, 540–548. [Google Scholar] [CrossRef]
  33. Shi, C.; Xie, Z.; Qian, H.; Liang, M.; Yang, X. China land soil moisture EnKF data assimilation based on satellite remote sensing data. Sci. China Earth Sci. 2011, 54, 1430–1440. [Google Scholar] [CrossRef]
  34. Yang, F.; Lu, H.; Yang, K.; He, J.; Wang, W.; Wright, J.S.; Li, C.; Han, M.; Li, Y. Evaluation of multiple forcing data sets for precipitation and shortwave radiation over major land areas of China. Hydrol. Earth Syst. Sci. 2017, 21, 5805–5821. [Google Scholar] [CrossRef]
  35. Liu, J.; Shi, C.; Sun, S.; Liang, J.; Yang, Z.-L. Improving Land Surface Hydrological Simulations in China Using CLDAS Meteorological Forcing Data. J. Meteorol. Res. 2019, 33, 1194–1206. [Google Scholar] [CrossRef]
  36. Cui, Y.Y.; Qin, J.; Jing, W.; Tan, G. Applicability evaluation of merged soil moisture in GLDAS and CLDAS products over Qinghai-Tibetan Plateau. Plateau Meteorol. 2018, 37, 123–136. [Google Scholar] [CrossRef]
  37. Shi, C.; Zhang, S.; Sun, S. Effect of Improved Precipitation CLDAS on Snow Simulation in China. Meteor. Mon. 2018, 44, 985–997. [Google Scholar] [CrossRef]
  38. Vanella, D.; Longo-Minnolo, G.; Belfiore, O.R.; Ramírez-Cuesta, J.M.; Pappalardo, S.E.; Consoli, S.; D’urso, G.; Chirico, G.B.; Coppola, A.; Comegna, A.; et al. Comparing the use of ERA5 reanalysis dataset and ground-based agrometeorological data under different climates and topography in Italy. J. Hydrol. 2022, 42, 101182. [Google Scholar] [CrossRef]
  39. Shea, J.G.; Worley, J.; Stern, R.N. An Introduction to Atmospheric and Oceanographic Datasets. (No.NCAR/TN-404+IA); University Corporation for Atmospheric Research: Boulder, CO, USA, 1994. [Google Scholar] [CrossRef]
  40. Sheffield, J.; Wood, E.F.; Roderick, M.L. Little change in global drought over the past 60 years. Nature 2012, 491, 435–438. [Google Scholar] [CrossRef]
  41. Boulard, D.; Castel, T.; Camberlin, P. Capability of a regional climate model to simulate climate variables requested for water balance computation: A case study over northeastern France. Clim. Dyn. 2015, 46, 2689–2716. [Google Scholar] [CrossRef]
  42. Paredes, P.; Rodrigues, G.C.; Alves, I.; Pereira, L.S. Partitioning evapotranspiration, yield prediction and economic returns of maize under various irrigation management strategies. Agric. Water Manag. 2014, 135, 27–39. [Google Scholar] [CrossRef]
  43. Shi, C.X.; Jiang, L.P.; Zhang, T.; Xu, B.; Han, S. Status and plans of CMA land data assimilation system (CLDAS) project. Geophys. Res. Abstr. 2014, 16, EGU2014-5671. [Google Scholar]
  44. Dee, D.P.; Uppala, S.M.; Simmons, A.J. The ERA-Interim reanalysis: Configuration and performance of the data assimilation system. Q. J. R. Meteorol. Soc. 2011, 137, 553–597. [Google Scholar] [CrossRef]
  45. Allen, R.G.; Smith, M.; Perrier, A.; Pereira, L.S. An update for the definition of reference evapotranspiration. ICID Bull. 1994, 43, 1–35. [Google Scholar]
  46. Takakura, T.; Kubota, C.; Sase, S.; Hayashi, M.; Ishii, M.; Takayama, K.; Nishina, H.; Kurata, K.; Giacomelli, G.A. Measurement of evapotranspiration rate in a single-span greenhouse using the energy-balance equation. Biosyst. Eng. 2009, 102, 298–304. [Google Scholar] [CrossRef]
  47. Mobilia, M.; Longobardi, A. Prediction of Potential and Actual Evapotranspiration Fluxes Using Six Meteorological Data-Based Approaches for a Range of Climate and Land Cover Types. ISPRS Int. J. Geo Inf. 2021, 10, 192. [Google Scholar] [CrossRef]
  48. Yin, Y.H.; Li, B.Y. A New Scheme for Climate Regionalization in China. Acta Geogr. Sin. 2010, 65, 3–12. [Google Scholar] [CrossRef]
  49. Zheng, J.Y.; Bian, J.J.; Ge, Q.S.; Hao, Z.X.; Yin, Y.H.; Liao, Y.M. The climate regionalization in China for 1981–2010. Geogr. Res. 2013, 32, 987–997. [Google Scholar]
  50. Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptual models part I—A discussion of principles. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
  51. Gupta, H.V.; Sorooshian, S.; Yapo, P.O. Status of Automatic Calibration for Hydrologic Models: Comparison with Multilevel Expert Calibration. J. Hydrol. Eng. 1999, 4, 135–143. [Google Scholar] [CrossRef]
  52. Tie, R.; Shi, C.; Wan, G.; Hu, X.J.; Kang, L.H.; Ge, L.L. CLDASSD: Reconstructing Fine Textures of the Temperature Field Using Super-Resolution Technology. Adv. Atmos. Sci. 2022, 39, 117–130. [Google Scholar] [CrossRef]
  53. Wu, L.; Qian, L.; Huang, G.; Liu, X.-G.; Wang, Y.-C.; Bai, H.; Wu, S.-F. Assessment of Daily of Reference Evapotranspiration Using CLDAS Product in Different Climate Regions of China. Water 2022, 14, 1744. [Google Scholar] [CrossRef]
  54. Qian, L.; Wu, L.F.; Liu, X.G.; Cui, Y.K.; Wang, Y.W. Comparison of CLDAS and Machine Learning Models for Reference Evapotranspiration Estimation under Limited Meteorological Data. Sustainability 2022, 14, 14577. [Google Scholar] [CrossRef]
  55. Simmons, A.J.; Willett, K.M.; Jones, P.D.; Thorne, P.W.; Dee, D.P. Low-frequency variations in surface atmospheric humidity, temperature, and precipitation: Inferences from reanalyses and monthly gridded observational data sets. J. Geophys. Res. 2010, 115, D01110. [Google Scholar] [CrossRef]
  56. Paredes, P.; Martins, D.S.; Pereira, L.S.; JorgePires, C. Accuracy of daily estimation of grass reference evapotranspiration using ERA-Interim reanalysis products with assessment of alternative bias correction schemes. Agric. Water Manag. 2018, 210, 340–353. [Google Scholar] [CrossRef]
  57. Bojanowski, J.S.; Vrieling, A.; Skidmore, A.K. A comparison of data sources for creating a long-term time series of daily gridded solar radiation for Europe. Sol. Energy 2014, 99, 152–171. [Google Scholar] [CrossRef]
  58. Urraca, R.; Martinez-De-Pison, E.; Sanz-Garcia, A.; Antonanzas, J.; Antonanzas-Torres, F. Estimation methods for global solar radiation: Case study evaluation of five different approaches in central Spain. Renew. Sustain. Energy Rev. 2017, 77, 1098–1113. [Google Scholar] [CrossRef]
  59. Fu, G.; Charles, S.P.; Timbal, B.; BranislavaOuyang, F. Comparison of NCEP-NCAR and ERA-Interim over Australia. Int. J. Clim. 2015, 36, 2345–2367. [Google Scholar] [CrossRef]
  60. Carvalho, D.; Rocha, A.; Gómez-Gesteira, M.; Santos, C.S. WRF wind simulation and wind energy production estimates forced by different reanalyses: Comparison with observed data for Portugal. Appl. Energy 2013, 117, 116–126. [Google Scholar] [CrossRef]
  61. Jabloun, M.; Sahli, A. Evaluation of FAO-56 methodology for estimating reference evapotranspiration using limited climatic data. Agric. Water Manag. 2008, 95, 707–715. [Google Scholar] [CrossRef]
  62. Betts, A.K.; Köhler, M.; Zhang, Y. Comparison of river basin hydrometeorology in ERA-Interim and ERA-40 reanalyses with observations. J. Geophys. Res. 2009, 114, D02101. [Google Scholar] [CrossRef]
  63. Dee, D.P.; Uppala, S. Variational bias correction of satellite radiance data in the ERA-Interim reanalysis. Q. J. R. Meteorol. Soc. 2009, 135, 1830–1841. [Google Scholar] [CrossRef]
  64. Su, T.; Feng, T.; Feng, G. Evaporation variability under climate warming in five reanalyses and its association with pan evaporation over China. J. Geophys. Res. Atmos. 2015, 120, 8080–8098. [Google Scholar] [CrossRef]
  65. Zhang, X.; Li, M.; Ma, Z. Evapotranspiration variability over global arid and semi-arid regions from 1982 to 2011. Chin. J Atmos. Sci. 2018, 42, 251–267. [Google Scholar] [CrossRef]
  66. Luo, H.; Cui, Y.; Duan, Z. Analysis of the Sensitivity of ET0 and the Main Meteorological Factors in Major Agricultural Regions in Tibet. Mod. Water Sav. Effic. Agric. Ecol. Irrig. Area Constr China 2010, 8, 130–136. [Google Scholar]
Figure 1. Spatial distribution of seven climatic zones in China and 689 utilized meteorological stations (The dots with different colors in the figure represent meteorological stations located in different climatic regions).
Figure 1. Spatial distribution of seven climatic zones in China and 689 utilized meteorological stations (The dots with different colors in the figure represent meteorological stations located in different climatic regions).
Water 15 02027 g001
Figure 2. Spatial distribution of solar radiation percent bias (PBias), root mean square error (RMSE), mean absolute error (MAE) and the correlation coefficient (R2) estimation with adopting CLDAS, ERA5, and station observation statistics.
Figure 2. Spatial distribution of solar radiation percent bias (PBias), root mean square error (RMSE), mean absolute error (MAE) and the correlation coefficient (R2) estimation with adopting CLDAS, ERA5, and station observation statistics.
Water 15 02027 g002
Figure 3. Spatial distribution of the annual mean PM-ET0 calculated from station observations and CLDAS, ERA5 reanalysis products in China.
Figure 3. Spatial distribution of the annual mean PM-ET0 calculated from station observations and CLDAS, ERA5 reanalysis products in China.
Water 15 02027 g003
Figure 4. Spatial distribution of PBias, RMSE and R2 comparing observations and reanalysis estimated PM-ET0 in various climates with different environments.
Figure 4. Spatial distribution of PBias, RMSE and R2 comparing observations and reanalysis estimated PM-ET0 in various climates with different environments.
Water 15 02027 g004
Figure 5. Comparing ET0 Obs with ET0 CLDAS and ET0 ERA5 on a daily scale, and using CLDAS, ERA5 reanalysis products of seven different climate zones (The blue circle in the figure represents the observed and estimated values, the dashed line is a 1:1 line, and the solid line is a scatter plot fitting curve).
Figure 5. Comparing ET0 Obs with ET0 CLDAS and ET0 ERA5 on a daily scale, and using CLDAS, ERA5 reanalysis products of seven different climate zones (The blue circle in the figure represents the observed and estimated values, the dashed line is a 1:1 line, and the solid line is a scatter plot fitting curve).
Water 15 02027 g005
Figure 6. Radar chart of statistical indicators of the PM-ET0 estimated by CLDAS and ERA5 reanalysis across seven climatic regions of China.
Figure 6. Radar chart of statistical indicators of the PM-ET0 estimated by CLDAS and ERA5 reanalysis across seven climatic regions of China.
Water 15 02027 g006
Figure 7. Time series of the monthly PM-ET0 from station observations, CLDAS and ERA5 reanalysis products during the period of 2017–2019 within the seven climates of China.
Figure 7. Time series of the monthly PM-ET0 from station observations, CLDAS and ERA5 reanalysis products during the period of 2017–2019 within the seven climates of China.
Water 15 02027 g007
Table 1. Statistical indicators for the minimum air temperature when comparing CLDAS and ERA5 reanalysis products with the observed values across seven climatic zones and mainland China.
Table 1. Statistical indicators for the minimum air temperature when comparing CLDAS and ERA5 reanalysis products with the observed values across seven climatic zones and mainland China.
Tmin CLDASTmin ERA5
Climate ZonesPBiasRMSEMAER2PBiasRMSEMAER2
(°C)(°C)
10.243.222.760.960.367.296.340.90
20.312.001.480.980.425.915.050.92
30.102.421.830.971.385.174.360.91
40.042.121.680.97−0.659.478.510.88
5−0.021.791.430.96−0.5713.5412.330.86
6−0.011.511.200.94−0.5316.8515.400.80
7−0.164.083.620.930.046.985.660.83
Average0.072.452.000.960.069.328.240.87
Table 2. Statistical indicators for the maximum air temperature when comparing CLDAS and ERA5 reanalysis products with the observed values across seven climatic zones and mainland China.
Table 2. Statistical indicators for the maximum air temperature when comparing CLDAS and ERA5 reanalysis products with the observed values across seven climatic zones and mainland China.
Tmax CLDASTmax ERA5
Climate ZonesPBiasRMSEMAER2PBiasRMSEMAER2
(°C)(°C)
1−0.144.243.670.95−0.328.407.350.87
2−0.033.252.500.93−0.296.315.410.90
3−0.013.192.450.94−0.285.564.640.89
4−0.033.262.590.92−0.3710.148.860.83
5−0.022.892.230.88−0.4112.7511.010.79
6−0.022.301.790.85−0.4315.0913.070.76
7−0.185.484.890.83−0.289.387.880.82
Average−0.063.522.870.90−0.349.668.320.84
Table 3. Statistical indicators of solar radiation originated from CLDAS and ERA5 reanalysis products in different climate zones and mainland China.
Table 3. Statistical indicators of solar radiation originated from CLDAS and ERA5 reanalysis products in different climate zones and mainland China.
Rs CLDASRs ERA5
Climate ZonesPBiasRMSEMAER2PBiasRMSEMAER2
(MJ m2 d−1)(MJ m2 d−1)
10.024.042.870.650.064.313.470.53
20.064.212.950.600.114.433.500.51
30.044.162.910.580.134.223.310.55
40.144.843.380.530.155.233.920.46
50.104.783.400.540.245.084.010.41
60.174.333.160.510.184.633.580.38
70.024.463.250.500.094.543.470.40
Average0.084.403.140.560.144.633.610.45
Table 4. Statistical indicators of relative humidity originated from CLDAS and ERA5 reanalysis products across different climate zones and mainland China.
Table 4. Statistical indicators of relative humidity originated from CLDAS and ERA5 reanalysis products across different climate zones and mainland China.
RHCLDASRHERA5
Climate ZonesPBiasRMSEMAER2PBiasRMSEMAER2
(%)(%)
1−0.3012.8410.290.58−0.3421.7518.800.45
2−0.2911.639.790.68−0.2618.5113.990.53
3−0.2010.828.230.59−0.2711.099.470.48
4−0.148.717.130.630.0415.4611.600.46
5−0.1410.248.270.51−0.1717.4713.080.41
6−0.077.916.320.62−0.2114.8310.170.48
70.3525.9422.510.370.2022.7018.840.32
Average−0.1212.5810.360.57−0.1417.4013.710.45
Table 5. Statistical indicators of wind speed originated from CLDAS and ERA5 reanalysis products across different climate zones of China.
Table 5. Statistical indicators of wind speed originated from CLDAS and ERA5 reanalysis products across different climate zones of China.
U2 CLDASU2 ERA5
Climate ZonesPBiasRMSEMAER2PBiasRMSEMAER2
(%)(%)
1−0.501.310.980.26−0.781.631.450.22
2−0.471.521.080.30−0.331.210.940.23
3−0.630.900.680.31−0.461.080.840.18
4−0.310.900.660.230.411.150.900.15
5−0.461.201.050.250.331.110.870.16
6−0.421.130.940.210.211.391.100.18
7−0.580.680.520.190.231.040.810.14
Average−0.481.090.840.25−0.061.230.990.18
Table 6. Comparing the statistical indicators relative to the estimation of ET0 by CLDAS. ERA5 reanalysis for 689 weather stations across the seven climatic regions and mainland China.
Table 6. Comparing the statistical indicators relative to the estimation of ET0 by CLDAS. ERA5 reanalysis for 689 weather stations across the seven climatic regions and mainland China.
ET0 CLDASET0 ERA5
Climate ZonesPBiasRMSEMAER2PBiasRMSEMAER2
mm d−1mm d−1
1−0.050.860.620.91−0.081.331.160.87
20.020.820.560.900.071.291.120.85
30.010.780.520.88−0.121.121.040.79
40.030.890.670.850.161.561.390.70
50.010.880.650.770.111.491.300.62
6−0.030.910.680.690.131.421.280.54
7−0.141.221.060.75−0.181.671.520.59
Average−0.010.910.680.820.011.421.280.70
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yu, X.; Qian, L.; Wang, W.; Huo, X.; Hu, X.; Wang, Y. Assessing and Comparing Reference Evapotranspiration across Different Climatic Regions of China Using Reanalysis Products. Water 2023, 15, 2027. https://doi.org/10.3390/w15112027

AMA Style

Yu X, Qian L, Wang W, Huo X, Hu X, Wang Y. Assessing and Comparing Reference Evapotranspiration across Different Climatic Regions of China Using Reanalysis Products. Water. 2023; 15(11):2027. https://doi.org/10.3390/w15112027

Chicago/Turabian Style

Yu, Xingjiao, Long Qian, Wen’e Wang, Xuefei Huo, Xiaotao Hu, and Yafei Wang. 2023. "Assessing and Comparing Reference Evapotranspiration across Different Climatic Regions of China Using Reanalysis Products" Water 15, no. 11: 2027. https://doi.org/10.3390/w15112027

APA Style

Yu, X., Qian, L., Wang, W., Huo, X., Hu, X., & Wang, Y. (2023). Assessing and Comparing Reference Evapotranspiration across Different Climatic Regions of China Using Reanalysis Products. Water, 15(11), 2027. https://doi.org/10.3390/w15112027

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