Next Article in Journal
Assimilating FY-3D MWHS2 Radiance Data to Predict Typhoon Muifa Based on Different Initial Background Conditions and Fast Radiative Transfer Models
Next Article in Special Issue
The SOLAR-HRS New High-Resolution Solar Spectra for Disk-Integrated, Disk-Center, and Intermediate Cases
Previous Article in Journal
Building a Practical Multi-Sensor Platform for Monitoring Vessel Activity near Marine Protected Areas: Case Studies from Urban and Remote Locations
Previous Article in Special Issue
Editorial on New Challenges in Solar Radiation, Modeling and Remote Sensing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diurnal Cycle in Surface Incident Solar Radiation Characterized by CERES Satellite Retrieval

State Key Laboratory of Earth Surface Processes and Resource Ecology, College of Global Change and Earth System Science, Faculty of Geographical Science, Beijing Normal University, Beijing 100875, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(13), 3217; https://0-doi-org.brum.beds.ac.uk/10.3390/rs15133217
Submission received: 11 May 2023 / Revised: 14 June 2023 / Accepted: 15 June 2023 / Published: 21 June 2023
(This article belongs to the Special Issue New Challenges in Solar Radiation, Modeling and Remote Sensing II)

Abstract

:
Surface incident solar radiation (Rs) plays an important role in climate change on Earth. Recently, the use of satellite-retrieved datasets to obtain global-scale Rs with high spatial and temporal resolutions has become an indispensable tool for research in related fields. Many studies were carried out for Rs evaluation based on the monthly satellite retrievals; however, few evaluations have been performed on their diurnal variation in Rs. This study used independently widely distributed ground-based data from the Baseline Surface Radiation Network (BSRN) to evaluate hourly Rs from the Clouds and the Earth’s Radiant Energy System Synoptic (CERES) SYN1deg–1Hour product through a detrended standardization process. Furthermore, we explored the influence of cloud cover and aerosols on the diurnal variation in Rs. We found that CERES-retrieved Rs performs better at midday than at 7:00–9:00 and 15:00–17:00. For spatial distribution, CERES-retrieved Rs performs better over the continent than over the island/coast and polar regions. The Bias, MAB and RMSE in CERES-retrieved Rs under clear-sky conditions are rather small, although the correlation coefficients are slightly lower than those under overcast-sky conditions from 9:00 to 15:00. In addition, the range in Rs bias caused by cloud cover is 1.97–5.38%, which is significantly larger than 0.31–2.52% by AOD.

1. Introduction

Surface incident solar radiation (Rs) plays an important role in climate change as it is the major energy source in the Earth system [1,2]. The performance of Rs is closely related to water cycle, as it significantly affects the evaporation of surface water [3]. Robock et al. [4] have shown that the changes in soil temperature are consistent with those in Rs during the global brightening and dimming period. Changes in Rs affect the melting and growth of glaciers as well [5]; for example, the snow cover in the northern hemisphere did not change significantly before the 1980s but exhibited a sharp downward trend in the brightening period after the 1980s [6]. Several studies suggest that Rs also has an impact on ecosystems, driving crop growth as an energy factor [7,8]. In addition, the spatial and temporal distribution of Rs is the basis for energy policy decision making on solar power [9]. Therefore, decreases and increases in Rs (also known as global dimming and brightening, respectively) have received widespread attention [10,11,12]. The variations in cloud cover, aerosols, and other factors were also explored, as they contributed more to the change in Rs [13,14]. Cloud cover regulates the radiative energy balance of the atmosphere [15]. Compared with a clear-sky atmosphere, cloud cover can absorb and reflect a large amount of incident shortwave radiation, which has a cooling effect on the subsurface. Aerosols can absorb and scatter solar radiation, which can block incident solar radiation, especially by reducing the passage of ultraviolet rays, weakening the solar radiation reaching the ground [16].
The diurnal cycle in Rs due to the Earth’s rotation can cause diurnal changes in the surface and atmospheric state in all regions except the polar regions (where solar variability is more seasonal) through a variety of physical processes [17]. Ye et al. [18] noted that the diurnal temperature range (DTR) decreased rapidly in China from the 1960s to the 1980s when it became significantly dimmer, while in the early 1990s, the decline in DTR stopped during the period when Rs changed from dimming to brightening, which was confirmed by Du et al. [19]. Both sea surface temperature variability in the western Pacific warm pool [20] and the diurnal behavior of the spiral rainbands [21] are modulated by the diurnal cycle in Rs.
Currently, Rs can be obtained from ground-based observations, reanalyses datasets, and satellite-retrieved datasets [22,23,24]. The ground-based observations of Rs have the highest accuracy but are sparsely distributed [9]. Reanalyses datasets have better spatial continuity, but the bias in cloud and aerosol simulations can lead to low accuracy in Rs [25,26]. The accuracy of the satellite retrieved data is relatively higher than that of reanalyses, as clouds and aerosols observed by satellite were used in the radiative transfer model [27].
Most studies have evaluated satellite-retrieved Rs on monthly and annual time scales. These studies have shown that satellite-retrieved datasets are less biased and have better accuracy than reanalyses datasets at long time scales [14,28,29]. However, few evaluations have been performed on satellite-retrieved Rs on the diurnal scale. Therefore, a comprehensive and detailed evaluation of satellite-retrieved Rs on an hourly scale is needed. In this study, we compared Clouds and the Earth’s Radiant Energy System Synoptic (CERES)-retrieved hourly Rs with Baseline Solar Radiation Network (BSRN) observations to quantify their differences in the diurnal cycle of Rs. Furthermore, the influence of clouds and aerosols on the diurnal variation in Rs was explored.

2. Data and Methods

2.1. Ground-Based Observation Data

Ground-based observation data with high-quality instrumentation and long-term maintenance provide the most reliable and accurate Rs [9]. In 1988, the WMO proposed the establishment of a new international ground-based radiation baseline network (i.e., BSRN) under the World Climate Research Program (WCRP) [30]. BSRN provides high-temporal resolution ground-based radiation observations (1 min) that can be used to validate satellite-retrieved datasets, improve radiative transfer calculations in climate models, and support the detection and monitoring of long-term changes in ground-based radiation fluxes [31].
In this study, we use ground-based observations from the BSRN to evaluate the satellite-retrieved dataset. As CERES data began in March 2000, the research period in this study is March 2000–July 2021. Seventy–three stations observed Rs during this period, as shown in Figure 1. Figure 2 shows the length of the record at each site. Twenty sites marked in red were excluded because their recorded periods were less than 4 years. Finally, 53 sites marked in blue color were used for evaluation in this work. Of these, 38 are located on the continent, 8 are located on the continent along the island/coast, and 8 are located on the continent in the polar region.
Thermophile radiometers suffer a negative bias because of the infrared loss to the sky at nighttime [32]. Consequently, measurement data at nighttime were abandoned in this study. We focus on the period from 7:00 to 17:00 each day, as the sample size of the data in this period is larger than 200,000 (see Figure 3), accounting for more than 50% of the data for each hour.

2.2. Satellite-Retrieved Dataset

Clouds and the Earth’s Radiant Energy System (CERES) is an investigation to provide earth radiation budget data through satellite sensing, which is produced, archived, and made available to the scientific community by the Atmospheric Sciences Data Center (ASDC), the Langley Research Center (LaRC), and the National Aeronautics and Space Administration (NASA) [33].
The CERES SYN Ed4.1 product contains hourly Rs, cloud cover, and aerosol optical depths with a 1° × 1° spatial resolution. Cloud properties from high-resolution imagers on several satellites were precisely matched with broadband radiance data. CERES has been providing clouds since 2000 using the algorithm developed for the second edition of CERES (Ed2) until 2002. To improve the accuracy of the clouds, CERES Edition 4 (Ed4) applies the revised algorithm [34]. The aerosols are from the Multi-scale Atmospheric Transport and Chemistry (MATCH) model constituents and the NASA-GSFC Moderate-resolution Imaging Spectroradiometer (MODIS) MOD04_L2/MYD04_L2 products [35]. The clouds and aerosols are used as inputs of the Fu–Liou radiative transfer model. Additional inputs are pressure, temperature and water vapor profiles from the Global Modeling and Assimilation Office (GMAO) Goddard Earth Observing System Model (GEOS). Gaseous absorption in the shortwave region is treated by the method described in Kato et al. [36], and absorption by water vapor, carbon dioxide, ozone, methane and oxygen were considered. Then, the adjusted fluxes were derived by constraining the calculated at the top of the atmosphere (TOA) fluxes to the observed CERES TOA fluxes. Rs in the CERES-SYN Ed4.1 products were computed hourly in approximate equal-area grid boxes [37].

2.3. Methods

Taking into account that the Rs measurement stations are distributed throughout the Earth, we converted UTC to local time for subsequent evaluation. To reduce the diurnal cycle and seasonal cycle in the Rs on the evaluation results, the observational data and CERES retrieved data were standardized as follows:
S W i , s t d = S W i S W i , T O A
where SWi,std represents the standardized radiation value at hour i, SWi represents the original radiation value at hour i and SWi,TOA represents the radiation value at the top of the atmosphere at hour i. Solar radiation at the top of atmosphere (SWTOA) is affected by the solar zenith angle and Sun–Earth distance, so SWTOA contains diurnal and seasonal variation, as shown in Figure 4. Both the observational Rs data and CERES retrieved Rs data were standardized by dividing SWTOA from CERES. The standardized Rs data are greater than or equal to 0 and less than 1 and unitless. It is worth noting that the Rs appearing in the following is the standardized one.
In this study, the Bias, mean absolute bias (MAB), root-mean-square error (RMSE), and correlation coefficient (R) are used to evaluate satellite radiation data with ground-based observations, as shown in Equations (2)–(5).
Bias = 1 n i = 1 n ( S i O i )
MAB = 1 n i = 1 n | S i O i |
RMSE = i = 1 n ( S i O i ) 2 n
R = i = 1 n ( S i  ​ S ¯ ) ( O i  ​ O ¯ ) i = 1 n ( S i  ​ S ¯ ) 2 i = 1 n ( O i  ​ O ¯ ) 2
where Si is the satellite-retrieved Rs at hour i and Oi is the observed Rs at hour i. The mean absolute bias was used here to avoid the offsetting of positive and negative deviations.

3. Results

3.1. Difference between CERES-Retrieved and BSRN Hourly Rs

Most Rs data are in the range from 0.50 to 0.75, as shown in Figure 5l. Hourly Rs were underestimated by −1.10% in CERES. The MAB between satellite-retrieved and observed Rs is less than 10%, while the RMSE is 12.88% for the whole period. As the seasonal and diurnal cycles were removed, their correlation coefficient was 0.83, which is much less than 0.92 for calculation with absolute values. This is consistent with the value of 0.95 noted by David et al. [35].
Evaluations were also conducted each hour, as shown in Figure 5a–k. Hourly Rs retrieved by CERES show negative Bias from 7:00 to17:00, with the largest Bias of −2.29% at 7:00 and the smallest Bias of −0.43% at 13:00. The MAB ranges from 8.77% to 8.92% for 10:00–13:00 and is relatively large, 10.03% and 9.89%, for 7:00 and 17:00, respectively. For other hours, MAB ranges from 9.09% to 9.52%. It decreases until midday and increases afterward. The RMSE is smallest at 10:00 (12.51%) and largest at 7:00 (13.66%). The diurnal variation of RMSE is similar to that of MAB. R is higher than 0.8 most of the time, except at 7:00 and 17:00. R is higher at midday and lower at 7:00 and 17:00. These results suggest that the CERES-retrieved Rs performs better at midday than at 7:00–8:00 and 16:00–17:00. The smaller difference between satellite-retrieved and observed Rs from 10:00 to 13:00 may be attributed to the orbit overpass times related to the two sensors, Terra and Aqua. Terra is in a descending sun-synchronous orbit with an equator-crossing time at 10:30, while Aqua is in ascending sun-synchronous orbits with an equator-crossing time at 13:30 [38].
Most of the stations show negative Bias, especially at 7:00 and 17:00, and it is particularly evident for stations located along the island/coast (Figure 6 and Table A1). MAB and RMSE are especially larger at 7:00 and 17:00 with relatively smaller R for most stations, as shown in Figure 7, Figure 8 and Figure 9 and Table A1. They perform better at midday. These are consistent with those shown in Figure 6. Notably, the largest Bias (−21.95% for Bias, 22.93% for MAB, 25.79% for RMSE, and 0.56 for R) was found at the IZA station located in Tenerife (Canary Islands, Spain) [39], which is consistent with the results of Hao et al. [40] and Tang et al. [22]. This may be attributed to its high altitude above the subtropical inversion layer, as cloud cover only affects the lower part of the area (below 2000 m), while the upper part of the island is cloudless [41].
Figure 10 and Table 1 summarize the statistical parameters for different types of stations. The negative Bias for the continental stations (−0.01–−1.38%) is smaller than that for the island/coastal stations (−1.84–−8.16%), and they are both underestimated more at 7:00 and 17:00 and less at midday. The stations located in the polar regions show the opposite variation, with negative Bias increasing until 9:00, then decreasing throughout the day and finally showing a positive Bias after 16:00. The MAB for continental stations ranges from 8.14% to 9.66%, which is smaller than that for island/coastal (10.67–13.54%) and polar (9.84–12.11%) stations. The RMSE for continental stations ranges from 11.44% to 12.51%, which again is smaller than that for island/coastal (13.94–16.93%) and polar (13.13–15.74%) stations. The R value for continental stations ranges from 0.74 to 0.80, which is higher than that for island/coastal (0.59–0.70) and polar (0.66–0.76) stations. The CERES-retrieved Rs performs better at most continental stations, although the spatial variabilities of these statistical parameters are relatively larger as they cover different land cover types, climate zones, and surface topography. Rapid weather changes and the presence of both land and water within the grid (edge effects) may lead to the worst performance along the island/coastal stations [40]. CERES-retrieved Rs data also perform relatively poorly at polar stations, which may be caused by the failure of cloud detection as more ice and snow exist in this region, and the temperature of clouds is usually not lower than that of surface snow and ice [31]. Additionally, Urraca et al. [41] pointed out the failure of most radiation products over polar regions, where strong intra-annual variations are present due to low solar elevation angles in winter, seasonal snowfall and the low viewing angle of satellites. The Rs difference over the continent is observed to be smaller than that along the island/coast and polar regions, which is confirmed by the findings in other studies [42,43].

3.2. Effect of Clouds and AOD on the Bias in CERES-Retrieved Hourly Rs

Cloud cover and aerosol optical depth (AOD) are two important factors that regulate Rs [13,44]. Some studies have shown that CERES-retrieved cloud cover is relatively accurate [34,45,46,47] and AOD has high accuracy with AERONET observations [48,49]. This section focuses on the CERES-retrieved Rs bias under different cloud cover and AOD (at 550 nm) categories.
Here, we define cloud cover of less than 20% as clear-sky, greater than 80% as overcast-sky and everything else as cloudy-sky. The Bias is smallest under clear-sky conditions as shown in Figure 11a, with the smallest value of approximately 0% at 9:00, slightly underestimated at midday (−0.16%) and overestimated at 14:00–17:00 (0.60%). Large underestimations in Rs were found under other two conditions especially for cloudy-sky conditions, with smaller negative Bias at midday (−1.36% for cloudy-sky, −0.09% for overcast-sky) and larger negative Bias at 7:00 and 17:00 (−3.84% for cloudy-sky, −1.92% for overcast-sky).
MAB and RMSE for all times under clear-sky conditions are significantly (3.33–7.01% for MAB and 5.84–10.31% for RMSE) smaller than those under cloudy-sky (9.96–11.91% for MAB and 13.42–15.20% for RMSE) and overcast-sky conditions (10.32–11.09% for MAB and 13.66–14.67% for RMSE), especially around midday (Figure 11b,c). Compared with cloudy-sky conditions, MAB and RMSE under overcast-sky conditions are 0.26–1.10% and 0.10–0.33% larger for 9:00–15:00. The diurnal cycle of MAB and RMSE is strongest (1.30% variation and 1.59% variation) under clear-sky conditions with low values (3.33% for MAB and 5.84% for RMSE) at midday and high values (7.01% for MAB and 10.31% for RMSE) at 7:00–9:00 and 15:00–17:00. They are rather stable for cloudy-sky (0.63% variation for MAB and 0.59% variation for RMSE) and overcast-sky (0.22% variation for MAB and 0.29% variation for RMSE) conditions.
Figure 11d shows that R under overcast-sky conditions (0.75–0.77) is higher than that under clear-sky conditions (0.69–0.74) for 9:00–15:00. CERES-retrieved Rs under cloudy-sky conditions has the lowest R for each hour, ranging from 0.49 at midday to ~0.60 at 7:00–9:00 and 16:00–17:00. MAB and RMSE under cloudy-sky conditions are similar to overcast-sky and all-sky conditions, but R is much lower. This may be attributed to the Bias under cloudy-sky conditions being 1.44% and 1.17% smaller than those under overcast-sky and all-sky conditions.
R is lower at midday and higher at 7:00 and 17:00 under clear-sky and cloudy-sky conditions, while it is the opposite under overcast-sky conditions. Under clear-sky and cloudy-sky conditions, MAB and RMSE decrease/increase while R also decreases/increases correspondingly, which is the opposite of the situation presented in Figure 10 (MAB and RMSE increase while R decreases and vice versa). This may be because more clouds are likely to appear at midday [50], and there are many broken clouds in low clouds. Broken clouds are smaller in size and more disperse in the sky, and their area may be much smaller than the single pixel of CERES satellite data. Thus, many broken clouds at these hours can easily be misjudged as clear-sky or cloudy-sky conditions. It is easier for the satellite to detect the circumstances under overcast-sky conditions. R is higher for all-sky conditions with more fluctuation than that for other conditions because Rs is more stable in the same cloud category.
Furthermore, we classified cloud cover more specifically into “0–20”, “20–40”, “40–60”, “60–80”, and “80–100” categories. The AOD values were divided into “0–0.05”, “0.05–0.1”, “0.1–0.15”, “0.15–0.3” and “0.3–8” categories (see Figure 12 and Figure 13). A negative Bias of Rs were found at almost all hours. Bias was the smallest, nearly zero, in the “0–20” and “80–100” cloud cover categories. This means that the Rs Bias is small under clear-sky conditions and cloudy-sky conditions. Bias, calculated by the mean method, is largest in the “60–80” cloud cover category for all hours, with the largest value of −4.71% at 7:00. Median values show that a larger negative Bias is shown in the “40–60” and “60–80” cloud cover categories, with the highest value of −5.38% at 7:00 for the “40–60” interval. At 17:00, the difference between the mean and the median values is the smallest. From Table 2, we found changes in cloud cover contribute to 1.97–5.38% changes in Rs Bias, and the variation in Rs Biases caused by cloud cover is in the range of 0.86–2.43%. The change/variation in the median values (4.39–5.38%/1.98–2.43%) is larger than that in the mean values (1.97–4.03%/0.86–1.75%) as shown in Table 2, which may be due to the positive values being offset by the negative values for mean calculations.
For the mean, at 7:00, the negative Bias decreases significantly with increasing AOD; from 8:00 to 14:00, the negative Bias decreases with increasing AOD for the “0–0.3” range and then increases for the “0.3–8” range; from 16:00 to 17:00, the negative Bias increases with increasing AOD. The median is generally similar to the mean except for the “0–0.05” interval. AOD changes lead to a change of 0.31–2.52% and a variation of 0.12–1.01% in the Rs Biases (Table 2). The use of median and mean values derived similar results.
In general, the change in Rs Bias caused by AOD at each hour is significantly smaller than that caused by cloud cover, which in turn supports that the diurnal variation in Biases in CERES-retrieved Rs is more sensitive to cloud cover than AOD.

4. Discussion

Based on the observed Rs data at different sites in the BSRN observational network, this study evaluates CERES-retrieved hourly Rs data and explores the impact of cloud cover and AOD on the bias of CERES-retrieved Rs data. The impact of diurnal and seasonal variation on the Rs evaluation was removed by dividing the solar radiation at TOA.
The Bias, MAB and RMSE values over the continent are 2.37%, 2.36% and 2.55% smaller than those over the island/coast and polar regions, respectively, while R is 0.1 higher. This may be attributed to rapid weather changes and the presence of both land and water within the grid (edge effects) for island/coast regions, and the failure of cloud detection as more ice and snow exist for polar regions. The spatial distribution of MAB and RMSE is consistent with the results of Yang et al. [27] and Tang et al. [22]. The Bias, MAB, RMSE at 7:00 are relatively large, then decrease at midday and increase again at 17:00. R is higher at midday than at 7:00 and 17:00. The diurnal variation in R is identical to that in Hao et al. [40], but other statistical parameters are different because we removed the diurnal cycle and seasonal cycle in the Rs. MAB and RMSE for all times under clear-sky conditions are significantly 3.31–7.71% and 3.35–8.59% smaller than those under other conditions, especially around midday. Bias under clear-sky conditions is about 1.47% smaller for overcast-sky conditions (except for nearly the same at 11:00–13:00) and 2.13% smaller for cloudy-sky conditions. However, R under clear-sky conditions is lower than that under overcast-sky conditions for 9:00–15:00. R under cloudy-sky conditions is significantly lower than those under other conditions. This may be because more clouds tend to occur at midday, and during these hours, many broken clouds can easily be misjudged as clear-sky or cloudy-sky conditions. Satellites are more likely to detect circumstances under overcast-sky conditions. The change in Rs bias caused by AOD is 0.31–2.52% at all hours, which is significantly smaller than 1.97–5.38% caused by cloud cover. This is consistent with our previous studies as cloud cover drive the short-term variation in Rs and aerosols modulate its long-term variation [51,52].

5. Conclusions

Satellite-retrieved datasets have relatively high accuracy and continuity in spatial distribution compared to traditional ground-based observations and reanalyses datasets. However, the evaluations of satellite-retrieved Rs data on the diurnal scale are still lacking. In this study, CERES-retrieved Rs on the diurnal scale was evaluated by comparison with the BSRN observations during a 21-year period from 2000 to 2021, and the influence of clouds and aerosols on the diurnal variation in Rs was investigated. The findings of this study included the following:
  • CERES-retrieved Rs performs better at 11:00–13:00 (−0.55% for Bias, 8.87% for MAB, 12.58% for RMSE, and 0.84 for R) than at other hours (1.26% for Bias, 10.00% for MAB, 13.50% for RMSE, and 0.81 for R).
  • For spatial distribution, CERES-retrieved Rs performs better over the continent (−0.42% for Bias, 8.89% for MAB, 12.12% for RMSE, and 0.83 for R) than over the island/coast (−1.01% for Bias, 9.38% for MAB, 13.00% for RMSE, and 0.74 for R) and polar (−1.70% for Bias, 10.85% for MAB, 14.30% for RMSE and 0.72 for R) regions.
  • The Bias, MAB, and RMSE in CERES-retrieved Rs under clear-sky conditions are rather small, although the correlation coefficients are slightly lower than those under overcast-sky conditions from 9:00 to 15:00. R in CERES-retrieved Rs under cloudy-sky conditions are the lowest.
  • The change in Rs bias caused by cloud cover is 1.97–5.38%, significantly larger than 0.31–2.52% by AOD.
In this study, the BSRN network stations used to evaluate the hourly Rs data retrieved by CERES have a high accuracy. However, the observation quality of the BSRN station, the urbanization level of the city where the station is located, and whether there are natural or artificial factors that affect the observation near the station have not been further evaluated and discussed in detail. Additionally, in the BSRN observation network, there are many missing and null values at some stations, which makes the ground-based observation sequence discontinuous. Although some studies have shown that the accuracy of cloud covers and AOD in CERES is relatively high [45,46,47,48,49], further evaluation of their diurnal variations is needed. It is beyond the scope of this study but is essential for our future work.

Author Contributions

Conceptualization, Q.M.; Methodology, Q.M.; Software, L.L.; Validation, L.L. and Q.M.; Formal Analysis, Q.M.; Investigation, Q.M.; Resources, L.L.; Data Curation, L.L.; Writing—Original Draft Preparation, L.L.; Writing—Review and Editing, Q.M.; Visualization, L.L.; Supervision, Q.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Science Foundation of China, grant number 41930970.

Data Availability Statement

Baseline Surface Radiation Network (BSRN) data are available at https://dataportals.pangaea.de/bsrn/?q=LR0100, accessed on 1 July 2021; Clouds and the Earth’s Radiant Energy System for CERES SYN data are available at https://ceres.larc.nasa.gov/order_data.php, accessed on 1 July 2021.

Acknowledgments

We thank the following institutions for sharing their data freely: the NASA Langley Research Center Atmospheric Science Data Center for CERES SYN data; the World Data Center PANGAEA for BSRN observation data.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. The information (station name, latitude (°), longitude (°) and altitude (m)) and evaluation of CERES hourly Rs against the ground-based measurements for 53 stations. N is sample size of participating in the calculation.
Table A1. The information (station name, latitude (°), longitude (°) and altitude (m)) and evaluation of CERES hourly Rs against the ground-based measurements for 53 stations. N is sample size of participating in the calculation.
SiteLatLongElStatistical
Parameters
7:008:009:0010:0011:0012:0013:0014:0015:0016:0017:00
ALE
(2004.08–2014.03)
−23.8133.9547Bias−2.04−2.18−1.69−0.650.11.231.962.853.124.154.82
MAB11.1211.0110.9310.439.859.799.829.279.409.499.40
RMSE14.2114.3914.2913.5813.2813.1813.1712.3112.6112.5612.35
R0.720.720.720.740.760.770.770.790.780.790.80
N16381715178718481886191218991839177917091633
ASP
(2000.03–2020.07)
71.32−156.68Bias−1.69−0.460.060.420.400.701.211.181.852.45−2.20
MAB6.755.114.193.884.284.674.815.35.896.379.32
RMSE9.688.147.287.027.668.068.158.749.259.4712.69
R0.780.790.820.830.820.800.840.830.820.810.70
N67096757676867786775676467586756675267546753
BAR
(2000.02–2017.08)
32.27−64.678Bias−3.05−3.10−3.36−3.37−2.91−1.91−0.730.792.374.295.75
MAB11.1810.8710.2710.129.699.369.168.968.999.459.86
RMSE14.8514.3813.7113.5712.8612.4812.0611.6711.5011.9012.22
R0.720.750.780.790.810.820.820.830.830.820.82
N34263804412643814555466746154448420639123570
BER
(2000.03–2017.08)
36.61−97.52317Bias−1.480.330.640.480.500.240.380.33−0.57−0.43−1.63
MAB8.188.258.598.578.958.949.078.788.958.699.13
RMSE10.9011.0011.5911.6012.0912.1212.2811.7911.9211.5412.07
R0.810.810.780.770.760.760.760.770.770.780.73
N47064704469747044697469546974700470247043994
BIL
(2000.03–2019.07)
40.07−88.37213Bias−0.261.881.731.551.291.351.251.391.772.690.45
MAB7.756.305.755.585.465.635.465.665.926.847.76
RMSE11.038.858.418.358.098.418.008.228.449.4510.83
R0.830.900.890.870.870.870.870.860.870.870.80
N60136775677767796781677967776780677967795774
BON
(2009.01–2020.04)
40.13−105.21689Bias0.490.50−0.94−1.23−1.46−1.09−0.81−0.81−0.620.812.09
MAB8.248.818.548.368.418.618.518.698.478.488.80
RMSE11.1312.1212.3212.1812.1312.2212.1312.2212.0211.6811.49
R0.830.850.860.870.860.860.860.850.850.860.86
N51837231724072557269727072707273727372707164
BOS
(2009.01–2020.04)
40.05−1051577Bias0.46−0.47−1.18−1.48−1.41−1.11−0.93−0.47−0.450.982.63
MAB10.9810.6010.4810.5711.5311.6812.3112.4612.2711.9111.44
RMSE15.3315.0314.8914.6715.6915.7116.3216.3715.9515.2514.55
R0.730.720.690.670.640.680.690.690.690.690.68
N72807315731673297332733373287326732773475107
BOU
(2000.02–2016.06)
−15.6−47.711023Bias−1.77−0.72−1.28−2.09−2.75−2.76−2.68−1.280.180.953.30
MAB12.3110.8810.2910.2811.2311.5912.2912.3012.2812.1411.50
RMSE16.8715.2614.5614.2515.2415.6316.3116.0515.7315.4714.56
R0.640.690.670.640.610.630.640.660.670.650.66
N57675805580958305826584258475846584658494083
BRB
(2006.02–2019.04)
51.974.930Bias−1.58−1.00−0.070.942.303.082.891.691.311.470.41
MAB7.406.456.156.938.339.4510.3910.4510.7010.6710.87
RMSE9.828.928.9910.0611.8813.1814.1814.1714.2914.2413.90
R0.800.810.800.750.710.690.690.700.720.720.63
N30853086309130893086308330863086308230843082
CAB
(2005.02–2021.04)
50.22−5.3288Bias−2.56−0.751.300.890.921.430.931.040.640.14−1.47
MAB9.629.188.768.438.408.328.458.518.608.929.15
RMSE12.3511.8811.4411.2911.2111.1511.3111.3111.3111.5811.97
R0.800.840.860.870.870.870.860.860.860.850.82
N39395045586458695866587158705869587158714323
CAM
(2001.01–2017.07)
44.085.06100Bias−0.16−0.031.511.571.711.961.521.150.29−0.88−1.91
MAB9.369.459.329.629.829.919.619.489.159.469.61
RMSE12.0612.1912.1612.5712.7012.9412.6612.4611.9212.1912.45
R0.800.810.830.830.830.830.830.840.850.830.81
N42415696569556965698569756975697569750633933
CAR
(2000.03–2018.12)
36.91−75.7137Bias−0.69−0.60−1.02−1.53−2.66−2.76−4.05−4.70−5.25−6.22−7.12
MAB7.728.558.848.899.059.119.8710.3010.9411.2911.75
RMSE10.7311.9012.6112.8112.8212.8913.6214.1915.0715.2215.71
R0.830.830.790.760.740.730.720.720.720.760.72
N49696664666466636663666266596658665866595296
CLN
(2000.05–2016.11)
42.82−1.6471Bias−5.33−3.04−2.14−1.48−1.84−1.93−2.14−3.13−4.91−7.60−9.58
MAB9.507.947.106.766.837.277.318.5310.1712.5614.20
RMSE12.1110.409.679.519.5810.1510.1911.4413.1215.7217.79
R0.860.880.880.880.870.870.870.850.820.780.66
N56475635561956165609559756055602562856474069
CNR
(2009.07–2021.05)
−12.1996.846Bias2.574.644.724.484.214.292.883.132.51−0.26−1.59
MAB9.289.629.169.218.848.808.588.758.969.169.25
RMSE12.4713.3812.8513.0412.6112.5512.0212.2312.1812.4012.29
R0.830.860.870.860.860.860.860.860.860.830.80
N40834328432043174317432243244325433243312978
COC
(2004.10–2020.05)
−30.6723.991287Bias−3.20−0.91−0.210.210.240.891.051.181.190.86−0.21
MAB10.139.278.708.458.648.458.528.648.668.748.39
RMSE13.1512.0411.6311.5511.8611.7011.7111.7211.5611.5810.84
R0.610.670.670.670.650.680.680.680.700.710.72
N42094214421742094197419041864179417241624146
DAA
(2000.07–2020.01)
−12.43130.930Bias−3.470.370.690.941.432.372.701.881.871.76−1.11
MAB9.605.954.904.584.614.965.375.746.537.739.24
RMSE13.539.938.648.428.769.509.919.7210.2311.2612.60
R0.750.810.800.790.780.780.800.790.760.750.77
N34083478350335093507350935133516351935163506
DAR
(2002.07–2015.01)
−75.1123.43233Bias−4.08−2.22−1.001.321.751.251.211.47−0.44−0.37−0.99
MAB8.357.737.808.229.319.208.097.977.697.517.90
RMSE10.9310.4710.7111.9113.1012.9211.9911.8210.9910.7510.56
R0.640.670.690.670.610.640.670.660.690.690.58
N45324533453345344533453345134516451745174518
DOM
(2006.01–2021.02)
36.63−1161007Bias−8.62−11.24−12.71−10.59−6.39−5.17−4.39−4.74−7.75−6.87−3.86
MAB12.5012.5413.6011.687.796.736.226.7610.1210.439.27
RMSE16.6816.9818.2115.3910.689.669.259.9514.4515.4813.95
R0.470.480.480.530.630.690.720.700.650.680.73
N25372752296331013271339434393378326331162919
DRA
(2009.01–2020.04)
−12.42130.932Bias0.891.290.58−0.56−0.77−0.48−0.51−0.310.632.275.03
MAB6.966.025.244.874.965.415.425.735.826.708.89
RMSE10.619.618.928.318.219.048.909.299.309.7612.42
R0.820.870.850.830.800.780.790.790.810.820.86
N56927215721072137214721672217223722372407016
DWN
(2008.04–2020.07)
36.61−97.49318Bias−5.32−2.63−0.381.932.632.052.202.28−0.52−0.86−2.53
MAB9.157.757.698.139.028.897.777.857.537.568.31
RMSE11.8010.1710.6611.9313.1612.7711.7811.8110.6410.6910.92
R0.640.670.690.670.590.610.660.660.680.700.60
N39994024402540304037404240414039403540274026
E13
(2000.02–2019.07)
−27.6−48.5211Bias−0.591.581.421.120.760.840.720.861.201.99−0.52
MAB7.776.265.705.545.395.585.375.565.796.617.74
RMSE10.998.768.348.308.018.367.898.148.309.2610.89
R0.830.900.890.870.870.870.870.870.870.870.80
N61256884689468906892689368946893689268925987
FLO
(2000.04–2021.03)
48.32−105.1634Bias0.130.000.30−0.140.510.19−0.45−1.38−0.59−0.570.25
MAB8.918.728.568.428.278.108.458.869.029.218.88
RMSE11.6711.6211.6911.4511.3810.9311.3211.7512.2412.3712.08
R0.810.850.860.860.860.870.860.860.850.830.72
N44854501450945104509450845104509449544764470
FPE
(2009.01–2020.05)
33.58130.43Bias1.261.210.060.560.690.480.00−0.790.482.313.33
MAB9.008.948.738.028.028.278.609.089.419.779.90
RMSE12.5012.4912.4611.4711.3511.6611.9912.6112.7613.0713.05
R0.830.840.820.830.820.810.800.800.800.790.76
N56316701673567376747675867786775674656944363
FUA
(2010.04–2021.04)
34.25−89.8798Bias−0.760.710.620.550.610.750.800.40−0.010.21−2.24
MAB7.407.227.207.387.667.887.717.557.487.598.10
RMSE9.619.589.7410.0310.3610.5510.4410.219.9610.0610.67
R0.880.910.910.900.900.900.900.900.900.890.84
N40144013401340124010400740104010401240143232
GCR
(2009.01–2020.04)
−23.5615.04407Bias3.193.021.921.090.730.630.510.080.020.681.21
MAB7.357.446.545.946.136.456.636.957.177.548.40
RMSE9.589.928.928.378.719.089.459.8510.0810.4211.35
R0.860.890.890.900.900.890.890.880.850.850.83
N52857093715872037203720472067205720571937150
GOB
(2012.05–2021.04)
−70.65−8.2542Bias2.004.001.261.090.750.350.62−0.29−1.54−4.38−7.75
MAB13.9110.987.234.312.942.592.492.763.546.069.20
RMSE17.3915.5011.527.945.765.014.925.365.918.6611.82
R0.330.530.700.780.780.730.760.780.700.620.66
N29033211321332143214321332133213321532153214
GVN
(2000.03–2021.01)
24.34124.25.7Bias−13.25−13.34−13.45−12.15−10.83−8.67−8.06−6.71−4.80−2.66−1.34
MAB16.6816.3016.1915.0114.0612.6012.2411.5410.9410.4510.06
RMSE20.9720.4020.0218.6917.8716.0515.8414.8514.1713.6413.21
R0.540.590.620.650.670.700.700.720.710.700.70
N38314381496353655736592559425640529848604364
ISH
(2010.04–2021.04)
28.31−16.52372.9Bias−1.25−0.85−0.55−0.270.440.861.160.94−0.23−1.18−2.66
MAB7.127.457.578.118.098.568.728.678.658.448.66
RMSE9.5110.0310.2711.1011.1911.9212.1512.0511.7411.3511.56
R0.850.870.870.850.840.830.830.830.840.850.84
N40094007400640034006400640044005400640084008
IZA
(2009.03–2021.06)
8.72167.710Bias−25.76−22.14−19.26−16.82−15.28−14.21−14.72−14.97−17.67−22.17−27.90
MAB26.6023.3620.9119.2618.6018.2618.5718.8121.0424.5528.93
RMSE28.9825.8323.4321.6820.8020.5120.8521.0123.2426.9331.66
R0.480.240.190.190.220.230.260.340.390.530.52
N44414435443944424445444444434441444444473662
KWA
(2000.03–2017.08)
−45.05169.7350Bias−3.74−0.82−0.49−0.31−0.34−0.95−0.24−0.020.910.771.11
MAB10.129.118.628.138.068.568.138.127.977.767.75
RMSE13.0912.0811.6011.1911.1411.7811.1711.2010.9110.4410.22
R0.550.670.710.720.710.710.720.740.740.750.75
N59455934587256695585548055965645590459585965
LAU
(2000.03–2018.12)
60.14−1.1880Bias−7.58−7.80−5.14−4.02−3.64−3.58−3.14−3.10−3.35−1.59−1.00
MAB15.9316.0514.4413.5913.6613.2813.0513.1013.2212.1411.38
RMSE20.5921.0018.3017.2517.5317.2417.0617.1817.4316.0715.02
R0.570.580.680.690.670.680.680.680.670.690.70
N43196132631963646391641364116409641164205318
LER
(2001.01–2017.07)
52.2114.12125Bias0.12−0.63−1.74−1.01−1.18−0.48−0.27−0.540.581.422.20
MAB11.4211.3011.4610.8510.7310.7110.8011.5011.6911.6711.60
RMSE14.1314.2414.8613.9213.8713.7313.8014.6914.8314.6614.34
R0.710.740.730.780.780.790.780.730.720.690.66
N37484378537553785380537953815381448438723235
LIN
(2000.03–2018.12)
−2.06147.46Bias−1.67−1.50−1.11−1.96−1.77−1.67−2.20−2.38−3.57−5.05−6.92
MAB8.819.199.229.099.239.359.569.5810.6911.2212.20
RMSE11.2411.9312.0912.1012.0812.4312.6712.7314.0414.6715.58
R0.860.850.860.860.860.850.840.830.790.750.70
N30353893389238943893389338933894389430312418
LRC
(2014.12–2021.05)
24.291547.1Bias−3.87−2.66−1.69−1.04−0.270.691.751.991.58−0.50−2.30
MAB10.3710.2310.2011.0311.6812.3912.5612.2411.6710.839.78
RMSE13.2313.2013.4114.5715.3616.1716.6016.1815.2113.9512.56
R0.700.730.730.710.710.700.680.690.700.690.60
N48874891489648894890488948894890488848864886
MAN
(2000.03–2013.10)
−0.52166.97Bias−8.90−3.62−1.82−0.690.040.360.370.630.23−0.71−0.93
MAB12.749.558.397.477.167.327.027.207.287.367.07
RMSE16.4612.6111.5010.6710.6310.9110.4610.5410.3210.109.37
R0.580.700.710.700.710.690.710.720.740.760.82
N39753976397739753969396439583954395939663972
MNM
(2010.04–2021.04)
78.9311.9311Bias−4.64−1.01−0.460.721.762.443.483.964.384.232.99
MAB9.368.137.757.908.398.679.249.509.539.368.97
RMSE12.2311.0110.8211.5212.4512.6913.3313.5013.1212.5511.77
R0.550.640.610.590.570.560.550.550.580.600.63
N48574860486148654863486648694865486448654866
NAU
(2000.03–2013.09)
48.712.21156Bias−9.97−8.79−7.43−6.45−4.67−3.18−1.93−2.12−3.84−3.68−3.65
MAB15.2214.2513.6213.0312.2311.9612.0611.2110.9910.7910.75
RMSE20.0318.5417.4816.8815.9415.7115.7914.5814.4514.1914.09
R0.600.640.670.680.690.680.640.700.730.730.72
N42994601486450575183526752005040483945794278
NYA
(2000.03–2021.07)
46.826.94491Bias−0.650.772.351.871.211.080.660.430.47−0.24−2.34
MAB8.238.528.748.188.648.679.048.879.169.6910.09
RMSE10.6611.3111.7411.1111.6211.4512.0011.7412.1012.5813.12
R0.850.860.870.880.860.870.850.860.840.810.77
N33634435499750005003500150014998499649953950
PAL
(2005.10–2019.12)
40.72−77.93376Bias1.393.313.282.721.892.010.740.660.720.41−1.12
MAB9.2610.3810.8910.6710.7410.3710.2310.2510.2010.059.51
RMSE12.2513.8114.7814.7114.6614.2313.8613.8413.6713.1312.56
R0.830.810.800.800.790.800.800.800.810.820.81
N55487500749774977496750375027505750875075674
PAY
(2000.03–2020.12)
−9.07−40.32387Bias4.792.612.161.591.451.260.740.530.711.230.51
MAB9.378.178.208.138.408.628.638.578.769.109.36
RMSE12.0811.2511.3711.2811.5111.8011.8211.7311.7612.0512.38
R0.850.880.890.890.880.880.870.870.860.840.78
N72447228722372287226722772327230723272255216
PSU
(2009.01–2020.04)
50.21−104.7578Bias1.412.984.365.075.094.623.502.361.750.882.49
MAB9.269.569.539.178.708.367.877.667.827.968.27
RMSE11.9612.4912.3711.9911.4010.8910.2410.0110.3110.5311.32
R0.640.690.710.730.730.730.690.660.640.630.61
N32463275328632873289328732853284328932883287
PTR
(2006.12–2018.07)
43.06141.317.2Bias−6.64−5.18−4.08−3.52−2.97−2.30−2.61−2.87−3.85−2.54−1.69
MAB12.5411.5610.909.949.238.989.279.7810.269.449.09
RMSE16.3015.7014.7813.8713.0112.6713.3013.8814.6213.2812.58
R0.720.740.770.800.810.810.790.780.780.790.78
N27053448427742894298430443074307430843063471
REG
(2000.02–2011.12)
30.8634.78500Bias−5.96−7.76−7.38−6.64−5.81−4.25−3.80−3.64−3.97−4.46−4.29
MAB11.8113.3313.6713.5413.5713.1413.0412.9212.5111.7010.83
RMSE15.1917.5118.1618.1318.2817.7417.5917.5616.9915.7514.29
R0.720.680.680.680.660.660.670.670.680.690.72
N30243872387338753872387538753873387338783089
SAP
(2010.04–2020.11)
−29.44−53.82489Bias1.762.842.512.401.872.442.202.473.464.877.86
MAB7.356.085.464.854.454.314.334.685.256.398.89
RMSE10.078.998.618.007.457.637.477.858.288.9511.80
R0.800.830.810.810.820.820.840.850.880.890.92
N31693472346934743469346834713472347434713469
SBO
(2003.01–2012.12)
47.0512.963108.9Bias−3.09−2.04−2.26−1.66−0.680.280.77−0.43−0.76−1.67−4.34
MAB8.197.246.826.276.386.717.077.227.417.719.12
RMSE10.649.769.418.939.239.9210.3410.3510.6610.8712.17
R0.890.910.910.910.900.890.890.880.870.880.86
N37573759375937583759375837593761376137623762
SMS
(2006.04–2017.06)
−89.98−24.82800Bias−15.77−15.15−15.26−13.77−12.10−9.58−8.98−7.89−7.96−8.05−7.42
MAB19.6019.0719.7319.5419.2719.1619.4418.7618.1117.8016.64
RMSE23.9123.2223.5623.5423.2523.1323.4422.8922.1122.0120.80
R0.680.670.640.600.580.530.500.490.500.480.49
N23302885288228552854286828662888291726051959
SPO
(2000.02–2017.03)
43.73−96.62473Bias−4.56−5.83−6.66−7.01−6.80−6.11−5.08−3.81−2.83−1.51−0.22
MAB8.329.079.399.649.318.908.407.587.477.467.71
RMSE11.1812.0912.4412.7912.5111.9511.3910.1210.099.9010.12
R0.700.690.690.680.700.710.730.760.750.760.75
N28002794277127622719273327522756277328072824
SXF
(2009.01–2019.11)
−69.0139.5918Bias3.392.540.751.281.361.351.451.251.352.732.54
MAB8.538.748.247.637.487.387.467.828.298.698.69
RMSE11.5612.3712.0811.1411.0210.7710.9511.3811.9011.9311.80
R0.860.850.850.860.860.860.860.850.850.850.83
N47865916593259335942594659435936594059384538
SYO
(2000.03–2021.04)
22.795.531385Bias6.695.494.674.284.735.295.645.846.888.529.45
MAB13.6412.9012.5812.6112.5712.2212.2111.9912.4612.5813.16
RMSE18.3117.7117.3417.2617.1316.7616.7816.6817.2117.3017.59
R0.600.620.640.650.660.690.690.690.680.700.69
N40844616506554565756591057835455506845864063
TAM
(2000.03–2021.06)
36.06140.125Bias−4.23−3.37−3.02−2.54−2.50−2.31−0.730.671.230.460.28
MAB7.025.765.134.985.286.046.486.957.337.967.96
RMSE9.027.616.856.867.248.399.7310.7011.1911.2110.73
R0.760.740.720.670.620.650.720.770.800.790.76
N72987314731873227324732173267324732673237323
TAT
(2000.03–2021.05)
58.2526.4670Bias−1.540.130.590.891.020.48−0.14−1.15−1.99−2.88−6.33
MAB8.857.967.647.427.648.448.408.738.949.3611.96
RMSE11.6910.6210.5210.2510.4211.3311.1511.5011.7512.1916.18
R0.790.860.880.890.880.860.850.840.840.830.70
N64437725772777267726772477257724772377247246
TOR
(2000.03–2020.11)
39.7511732Bias−1.07−1.11−1.06−0.41−0.75−0.49−0.200.01−0.86−1.32−2.18
MAB10.5810.079.408.578.558.828.989.7410.7111.1811.41
RMSE13.4513.0512.3711.3611.3011.6711.8412.6413.9314.3414.33
R0.800.840.850.880.880.870.870.840.790.750.72
N52636247756775677564756975697571675056574735
XIA
(2005.01–2015.10)
82.49−62.42127Bias−1.602.072.342.271.60−2.43−2.01−2.23−1.61−1.21−1.65
MAB8.937.537.356.966.626.726.596.847.096.836.79
RMSE11.6710.0710.069.539.159.048.878.989.178.838.87
R0.780.870.880.880.880.900.900.890.870.870.84
N36353650365036503651365136523653365336532784

References

  1. Zhao, R.Z.; Wang, K.C.; Wu, G.C.; Zhou, C.L. Temperature annual cycle variations and responses to surface solar radiation in China between 1960 and 2016. Int. J. Climatol. 2021, 41, E2959–E2978. [Google Scholar] [CrossRef]
  2. Kodera, K.; Kuroda, Y. Dynamical response to the solar cycle. J. Geophys. Res. Atmos. 2002, 107, 4749. [Google Scholar] [CrossRef] [Green Version]
  3. Papaioannou, G.; Kitsara, G.; Athanasatos, S. Impact of global dimming and brightening on reference evapotranspiration in Greece. J. Geophys. Res. Atmos. 2011, 116, D9. [Google Scholar] [CrossRef]
  4. Robock, A. Forty-five years of observed soil moisture in the Ukraine: No summer desiccation (yet). Geophys. Res. Lett. 2005, 32, L03401. [Google Scholar] [CrossRef] [Green Version]
  5. Pellicciotti, F.; Brock, B.; Strasser, U.; Burlando, P.; Funk, M.; Corripio, J. An enhanced temperature-index glacier melt model including the shortwave radiation balance: Development and testing for Haut Glacier d’Arolla, Switzerland. J. Glaciol. 2005, 51, 573–587. [Google Scholar] [CrossRef]
  6. Wild, M.; Trüssel, B.; Ohmura, A.; Long, C.N.; König-Langlo, G.; Dutton, E.G.; Tsvetkov, A. Global dimming and brightening: An update beyond 2000. J. Geophys. Res. Atmos. 2009, 114, D10. [Google Scholar] [CrossRef] [Green Version]
  7. Ballare, C.L.; Caldwell, M.M.; Flint, S.D.; Robinson, A.; Bornman, J.F. Effects of solar ultraviolet radiation on terrestrial ecosystems. Patterns, mechanisms, and interactions with climate change. Photochem. Photobiol. Sci. 2011, 10, 226–241. [Google Scholar] [CrossRef]
  8. Caldwell, M.M.; Bornman, J.F.; Ballare, C.L.; Flint, S.D.; Kulandaivelu, G. Terrestrial ecosystems, increased solar ultraviolet radiation, and interactions with other climate change factors. Photochem. Photobiol. Sci. 2007, 6, 252–266. [Google Scholar] [CrossRef]
  9. Korany, M.; Boraiy, M.; Eissa, Y.; Aoun, Y.; Abdel Wahab, M.M.; Alfaro, S.C.; Blanc, P.; El-Metwally, M.; Ghedira, H.; Hungershoefer, K.; et al. A database of multi-year (2004–2010) quality-assured surface solar hourly irradiation measurements for the Egyptian territory. Earth Syst. Sci. Data 2016, 8, 105–113. [Google Scholar] [CrossRef] [Green Version]
  10. Wang, K.; Dickinson, R.E.; Ma, Q.; Augustine, J.A.; Wild, M. Measurement Methods Affect the Observed Global Dimming and Brightening. J. Clim. 2013, 26, 4112–4120. [Google Scholar] [CrossRef]
  11. Allen, R.J.; Norris, J.R.; Wild, M. Evaluation of multidecadal variability in CMIP5 surface solar radiation and inferred underestimation of aerosol direct effects over Europe, China, Japan, and India. J. Geophys. Res. Atmos. 2013, 118, 6311–6336. [Google Scholar] [CrossRef]
  12. Xia, X. A closer looking at dimming and brightening in China during 1961–2005. Ann. Geophys. 2010, 28, 1121–1132. [Google Scholar] [CrossRef]
  13. Obregón, M.Á.; Serrano, A.; Costa, M.J.; Silva, A.M. Global Spatial and Temporal Variation of the Combined Effect of Aerosol and Water Vapour on Solar Radiation. Remote Sens. 2021, 13, 708. [Google Scholar] [CrossRef]
  14. Wang, K.; Ma, Q.; Li, Z.; Wang, J. Decadal variability of surface incident solar radiation over China: Observations, satellite retrievals, and reanalyses. J. Geophys. Res. Atmos. 2015, 120, 6500–6514. [Google Scholar] [CrossRef]
  15. Li, J.; You, Q.; He, B. Distinctive spring shortwave cloud radiative effect and its inter-annual variation over southeastern China. Atmos. Sci. Lett. 2020, 21, e970. [Google Scholar] [CrossRef] [Green Version]
  16. Wang, Y.; Lyu, R.; Xie, X.; Meng, Z.; Huang, M.; Wu, J.; Mu, H.; Yu, Q.R.; He, Q.; Cheng, T. Retrieval of gridded aerosol direct radiative forcing based on multiplatform datasets. Atmos. Meas. Tech. 2020, 13, 575–592. [Google Scholar] [CrossRef]
  17. Zhang, Y.; Long, C.N.; Rossow, W.B.; Dutton, E.G. Exploiting diurnal variations to evaluate the ISCCP-FD flux calculations and radiative-flux-analysis-processed surface observations from BSRN, ARM, and SURFRAD. J. Geophys. Res. 2010, 115, D15. [Google Scholar] [CrossRef] [Green Version]
  18. Ye, J.; Li, F.; Sun, G.; Guo, A. Solar dimming and its impact on estimating solar radiation from diurnal temperature range in China, 1961–2007. Theor. Appl. Climatol. 2010, 101, 137–142. [Google Scholar] [CrossRef]
  19. Du, J.; Wang, K.; Wang, J.; Jiang, S.; Zhou, C. Diurnal Cycle of Surface Air Temperature within China in Current Reanalyses: Evaluation and Diagnostics. J. Clim. 2018, 31, 4585–4603. [Google Scholar] [CrossRef]
  20. Shinoda, T. Impact of the Diurnal Cycle of Solar Radiation on Intraseasonal SST Variability in the Western Equatorial Pacific. J. Clim. 2005, 18, 2628–2636. [Google Scholar] [CrossRef]
  21. Zhou, S.; Ma, Y.; Ge, X. Impacts of the Diurnal Cycle of Solar Radiation on Spiral Rainbands. Adv. Atmos. Sci. 2016, 33, 1085–1095. [Google Scholar] [CrossRef]
  22. Tang, W.; Yang, K.; Qin, J.; Li, X.; Niu, X. A 16-year dataset (2000–2015) of high-resolution (3 h, 10 km) global surface solar radiation. Earth Syst. Sci. Data 2019, 11, 1905–1915. [Google Scholar] [CrossRef] [Green Version]
  23. Jiang, H.; Lu, N.; Qin, J.; Yao, L. Hourly 5-km surface total and diffuse solar radiation in China, 2007–2018. Sci. Data 2020, 7, 311. [Google Scholar] [CrossRef] [PubMed]
  24. Letu, H.; Nakajima, T.Y.; Wang, T.; Shang, H.; Ma, R.; Yang, K.; Baran, A.J.; Riedi, J.; Ishimoto, H.; Yoshida, M.; et al. A New Benchmark for Surface Radiation Products over the East Asia–Pacific Region Retrieved from the Himawari-8/AHI Next-Generation Geostationary Satellite. Bull. Am. Meteorol. Soc. 2022, 103, E873–E888. [Google Scholar] [CrossRef]
  25. He, Y.; Wang, K.; Feng, F. Improvement of ERA5 over ERA-Interim in Simulating Surface Incident Solar Radiation throughout China. J. Clim. 2021, 34, 3853–3867. [Google Scholar] [CrossRef]
  26. Träger-Chatterjee, C.; Müller, R.W.; Trentmann, J.; Bendix, J. Evaluation of ERA-40 and ERA-interim re-analysis incoming surface shortwave radiation datasets with mesoscale remote sensing data. Meteorol. Z. 2010, 19, 631–640. [Google Scholar] [CrossRef]
  27. Yang, D.; Bright, J.M. Worldwide validation of 8 satellite-derived and reanalysis solar radiation products: A preliminary evaluation and overall metrics for hourly data over 27 years. Sol. Energy 2020, 210, 3–19. [Google Scholar] [CrossRef]
  28. Ma, Q.; Wang, K.; Wild, M. Impact of geolocations of validation data on the evaluation of surface incident shortwave radiation from Earth System Models. J. Geophys. Res. Atmos. 2015, 120, 6825–6844. [Google Scholar] [CrossRef] [Green Version]
  29. White, J.W.; Hoogenboom, G.; Wilkens, P.W.; Stackhouse, P.W., Jr.; Hoel, J.M. Evaluation of Satellite-Based, Modeled-Derived Daily Solar Radiation Data for the Continental United States. Agron. J. 2011, 103, 1242–1251. [Google Scholar] [CrossRef] [Green Version]
  30. Driemel, A.; Augustine, J.; Behrens, K.; Colle, S.; Cox, C.; Cuevas-Agullo, E.; Denn, F.M.; Duprat, T.; Fukuda, M.; Grobe, H.; et al. Baseline Surface Radiation Network (BSRN): Structure and data description (1992–2017). Earth Syst. Sci. Data 2018, 10, 1491–1501. [Google Scholar] [CrossRef] [Green Version]
  31. Tang, W.; Qin, J.; Yang, K.; Zhu, F.; Zhou, X. Does ERA5 outperform satellite products in estimating atmospheric downward longwave radiation at the surface? Atmos. Res. 2021, 252, 105453. [Google Scholar] [CrossRef]
  32. Dutton, E.G.; Michalsky, J.J.; Stoffel, T.; Forgan, B.W.; Hickey, J.; Nelson, D.W.; Alberta, T.L.; Reda, I. Measurement of Broadband Diffuse Solar Irradiance Using Current Commercial Instrumentation with a Correction for Thermal Offset Errors. J. Atmos. Ocean. Technol. 2001, 18, 297–314. [Google Scholar] [CrossRef]
  33. Wielicki, B.A.; Barkstrom, B.R.; Harrison, E.F.; Lee, R.B.; Smith, G.L.; Cooper, J.E. Clouds and the Earth’s Radiant Energy System (CERES): An Earth Observing System Experiment. Bull. Am. Meteorol. Soc. 1996, 77, 853–868. [Google Scholar] [CrossRef]
  34. Trepte, Q.Z.; Minnis, P.; Sun-Mack, S.; Yost, C.R.; Chen, Y.; Jin, Z.H.; Hong, G.; Chang, F.L.; Smith, W.L.; Bedka, K.M.; et al. Global Cloud Detection for CERES Edition 4 Using Terra and Aqua MODIS Data. IEEE Trans. Geosci. Remote Sens. 2019, 57, 9410–9449. [Google Scholar] [CrossRef]
  35. Rutan, D.A.; Kato, S.; Doelling, D.R.; Rose, F.G.; Nguyen, L.T.; Caldwell, T.E.; Loeb, N.G. CERES Synoptic Product: Methodology and Validation of Surface Radiant Flux. J. Atmos. Ocean. Technol. 2015, 32, 1121–1143. [Google Scholar] [CrossRef]
  36. Kratz, D.P.; Rose, F.G. Accounting for molecular absorption within the spectral range of the CERES window channel. J. Quant. Spectrosc. Radiat. Transf. 1999, 61, 83–95. [Google Scholar] [CrossRef]
  37. Collins, W.D.; Rasch, P.J.; Eaton, B.E.; Khattatov, B.V.; Lamarque, J.-F.; Zender, C.S. Simulating aerosols using a chemical transport model with assimilation of satellite aerosol retrievals: Methodology for INDOEX. J. Geophys. Res. Atmos. 2001, 106, 7313–7336. [Google Scholar] [CrossRef] [Green Version]
  38. Loeb, N.G.; Doelling, D.R.; Wang, H.; Su, W.; Nguyen, C.; Corbett, J.G.; Liang, L.; Mitrescu, C.; Rose, F.G.; Kato, S. Clouds and the Earth’s Radiant Energy System (CERES) Energy Balanced and Filled (EBAF) Top-of-Atmosphere (TOA) Edition-4.0 Data Product. J. Clim. 2018, 31, 895–918. [Google Scholar] [CrossRef]
  39. García, R.D.; Cuevas, E.; Ramos, R.; Cachorro, V.E.; Redondas, A.; Moreno-Ruiz, J.A. Description of the Baseline Surface Radiation Network (BSRN) station at the Izaña Observatory (2009–2017): Measurements and quality control/assurance procedures. Geosci. Instrum. Method Data Syst. 2019, 8, 77–96. [Google Scholar] [CrossRef] [Green Version]
  40. Hao, D.; Asrar, G.R.; Zeng, Y.; Zhu, Q.; Wen, J.; Xiao, Q.; Chen, M. DSCOVR/EPIC-derived global hourly and daily downward shortwave and photosynthetically active radiation data at 0.1°  ×  0.1° resolution. Earth Syst. Sci. Data 2020, 12, 2209–2221. [Google Scholar] [CrossRef]
  41. Urraca, R.; Huld, T.; Gracia-Amillo, A.; Martinez-de-Pison, F.J.; Kaspar, F.; Sanz-Garcia, A. Evaluation of global horizontal irradiance estimates from ERA5 and COSMO-REA6 reanalyses using ground and satellite-based data. Sol. Energy 2018, 164, 339–354. [Google Scholar] [CrossRef]
  42. Boland, J.; David, M.; Lauret, P. Short term solar radiation forecasting: Island versus continental sites. Energy 2016, 113, 186–192. [Google Scholar] [CrossRef]
  43. Wang, H.; Pinker, R.T. Shortwave radiative fluxes from MODIS: Model development and implementation. J. Geophys. Res. Atmos. 2009, 114, D19. [Google Scholar] [CrossRef]
  44. Redemann, J.W.; Wood, R.; Zuidema, P.; Doherty, S.J.; Luna, B.; LeBlanc, S.E.; Diamond, M.S.; Shinozuka, Y.; Chang, I.Y.; Ueyama, R.; et al. An overview of the ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) project: Aerosol–cloud–radiation interactions in the southeast Atlantic basin. Atmos. Chem. Phys. 2021, 21, 1507–1563. [Google Scholar] [CrossRef]
  45. Levy, R.C.; Mattoo, S.; Sawyer, V.; Shi, Y.; Colarco, P.R.; Lyapustin, A.I.; Wang, Y.; Remer, L.A. Exploring systematic offsets between aerosol products from the two MODIS sensors. Atmos. Meas. Tech. 2018, 11, 4073–4092. [Google Scholar] [CrossRef] [Green Version]
  46. Yost, C.R.; Minnis, P.; Sun-Mack, S.; Chen, Y.; Smith, W.L. CERES MODIS Cloud Product Retrievals for Edition 4—Part II: Comparisons to CloudSat and CALIPSO. IEEE Trans. Geosci. Remote Sens. 2021, 59, 3695–3724. [Google Scholar] [CrossRef]
  47. Minnis, P.; Sun-Mack, S.; Chen, Y.; Chang, F.L.; Yost, C.R.; Smith, W.L.; Heck, P.W.; Arduini, R.F.; Bedka, S.T.; Yi, Y.; et al. CERES MODIS Cloud Product Retrievals for Edition 4—Part I: Algorithm Changes. IEEE Trans. Geosci. Remote Sens. 2021, 59, 2744–2780. [Google Scholar] [CrossRef]
  48. Kaufman, Y.J.; Tanré, D.; Boucher, O. A satellite view of aerosols in the climate system. Nature 2002, 419, 215–223. [Google Scholar] [CrossRef]
  49. Quaas, J.; Boucher, O.; Bellouin, N.; Kinne, S. Satellite-based estimate of the direct and indirect aerosol climate forcing. J. Geophys. Res. Atmos. 2008, 113, 947–955. [Google Scholar] [CrossRef]
  50. An, N.; Wang, K.; Zhou, C.; Pinker, R.T. Observed Variability of Cloud Frequency and Cloud-Base Height within 3600 m above the Surface over the Contiguous United States. J. Clim. 2017, 30, 3725–3742. [Google Scholar] [CrossRef]
  51. Ma, Q.; Wang, K.C.; He, Y.Y.; Su, L.Y.; Wu, Q.Z.; Liu, H.; Zhang, Y.R. Homogenized century-long surface incident solar radiation over Japan. Earth Syst. Sci. Data 2022, 14, 463–477. [Google Scholar] [CrossRef]
  52. Wang, K.C.; Dickinson, R.E.; Su, L.; Trenberth, K.E. Contrasting trends of mass and optical properties of aerosols over the Northern Hemisphere from 1992 to 2011. Atmos. Chem. Phys. 2012, 12, 9387–9398. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Geographical distribution of observation sites used for the evaluation of satellite-retrieved Rs data. The blue, red, and magenta circles indicate the site location in the polar (Arctic or Antarctic), continent, and island/coast, respectively. The black crosses indicate the sites excluded in this study as fewer data were available.
Figure 1. Geographical distribution of observation sites used for the evaluation of satellite-retrieved Rs data. The blue, red, and magenta circles indicate the site location in the polar (Arctic or Antarctic), continent, and island/coast, respectively. The black crosses indicate the sites excluded in this study as fewer data were available.
Remotesensing 15 03217 g001
Figure 2. Time duration of each site in BSRN for 2000–2021. Red indicates that sites with short measurement periods are excluded from this work; blue indicates that sites are used in this work.
Figure 2. Time duration of each site in BSRN for 2000–2021. Red indicates that sites with short measurement periods are excluded from this work; blue indicates that sites are used in this work.
Remotesensing 15 03217 g002
Figure 3. The sample size of the observation data at different times. The red line indicates 50% of the total sample size.
Figure 3. The sample size of the observation data at different times. The red line indicates 50% of the total sample size.
Remotesensing 15 03217 g003
Figure 4. (a) Monthly solar radiation at TOA at XIA station in 2005; (b) hourly solar radiation at TOA at XIA station on 25 April 2005.
Figure 4. (a) Monthly solar radiation at TOA at XIA station in 2005; (b) hourly solar radiation at TOA at XIA station on 25 April 2005.
Remotesensing 15 03217 g004
Figure 5. Scatter plots of annual average of hourly CERES retrieved and observed Rs from 2000 to 2021 for each hour from 7:00 to 17:00 and all hours. The statistical parameters in the upper left corner of (al) were calculated using Rs at 53 stations in Figure 1, and N is sample hours participating in the calculation.
Figure 5. Scatter plots of annual average of hourly CERES retrieved and observed Rs from 2000 to 2021 for each hour from 7:00 to 17:00 and all hours. The statistical parameters in the upper left corner of (al) were calculated using Rs at 53 stations in Figure 1, and N is sample hours participating in the calculation.
Remotesensing 15 03217 g005
Figure 6. Bias between hourly CERES-retrieved Rs and observed Rs in each station for individual hours from 7:00 to 17:00. Unit: %.
Figure 6. Bias between hourly CERES-retrieved Rs and observed Rs in each station for individual hours from 7:00 to 17:00. Unit: %.
Remotesensing 15 03217 g006
Figure 7. Mean absolute bias (MAB) between hourly CERES-retrieved Rs and observed Rs in each station for individual hours from 7:00 to 17:00. Unit: %.
Figure 7. Mean absolute bias (MAB) between hourly CERES-retrieved Rs and observed Rs in each station for individual hours from 7:00 to 17:00. Unit: %.
Remotesensing 15 03217 g007
Figure 8. Root-mean-square error (RMSE) between hourly CERES-retrieved Rs and observed Rs in each station for individual hours from 7:00 to 17:00. Unit: %.
Figure 8. Root-mean-square error (RMSE) between hourly CERES-retrieved Rs and observed Rs in each station for individual hours from 7:00 to 17:00. Unit: %.
Remotesensing 15 03217 g008
Figure 9. Correlation coefficients (R) between hourly CERES-retrieved Rs and observed Rs in each station for individual hours from 7:00 to 17:00.
Figure 9. Correlation coefficients (R) between hourly CERES-retrieved Rs and observed Rs in each station for individual hours from 7:00 to 17:00.
Remotesensing 15 03217 g009
Figure 10. Diurnal variations of statistical parameters between hourly CERES-retrieved Rs and observed Rs for different types of sites. (a) Bias %; (b) Mean absolute bias (MAB) %; (c) Root-mean-square error (RMSE) %; and (d) Correlation coefficient (R).
Figure 10. Diurnal variations of statistical parameters between hourly CERES-retrieved Rs and observed Rs for different types of sites. (a) Bias %; (b) Mean absolute bias (MAB) %; (c) Root-mean-square error (RMSE) %; and (d) Correlation coefficient (R).
Remotesensing 15 03217 g010
Figure 11. Diurnal variations of statistical parameters for CERES retrieved hourly Rs at all ground-based sites for different cloud cover conditions. (a) Bias %; (b) Mean absolute bias (MAB) %; (c) Root-mean-square error (RMSE) %; and (d) Correlation coefficient (R).
Figure 11. Diurnal variations of statistical parameters for CERES retrieved hourly Rs at all ground-based sites for different cloud cover conditions. (a) Bias %; (b) Mean absolute bias (MAB) %; (c) Root-mean-square error (RMSE) %; and (d) Correlation coefficient (R).
Remotesensing 15 03217 g011
Figure 12. Median and mean values of differences between CERES and BSRN hourly Rs under different cloud cover categories from 7:00 to 17:00. The differences refer to the values that CERES data minus BSRN data.
Figure 12. Median and mean values of differences between CERES and BSRN hourly Rs under different cloud cover categories from 7:00 to 17:00. The differences refer to the values that CERES data minus BSRN data.
Remotesensing 15 03217 g012
Figure 13. Median and mean values of differences between CERES and BSRN hourly Rs under different AOD categories from 7:00 to 17:00. The differences refer to the values that CERES data minus BSRN data.
Figure 13. Median and mean values of differences between CERES and BSRN hourly Rs under different AOD categories from 7:00 to 17:00. The differences refer to the values that CERES data minus BSRN data.
Remotesensing 15 03217 g013
Table 1. Evaluation of CERES hourly Rs against the ground-based measurements for different regions. Unit: Bias %, MAB %, RMSE %.
Table 1. Evaluation of CERES hourly Rs against the ground-based measurements for different regions. Unit: Bias %, MAB %, RMSE %.
RegionStatistical Parameters7:008:009:0010:0011:0012:0013:0014:0015:0016:0017:00
ContinentalBias−1.38−0.43−0.33−0.15−0.060.07−0.01−0.17−0.42−0.56−1.28
MAB9.308.728.368.148.248.388.438.598.809.119.66
RMSE12.3612.0511.8211.5211.4411.5611.6011.7411.8611.9712.51
R0.770.800.800.800.780.780.790.790.790.780.74
Island/coastalBias−8.16−5.59−4.27−3.36−2.67−2.25−1.84−1.76−2.34−2.84−3.90
MAB13.5412.2511.3810.8010.7110.6810.7110.6710.9611.1711.55
RMSE16.9315.5114.5414.0414.0614.0814.0813.9414.1214.1814.38
R0.590.620.630.630.620.620.630.650.670.700.70
PolarBias−3.17−3.76−4.07−3.58−2.76−1.84−1.28−0.74−0.250.891.81
MAB12.1111.9211.9711.2710.4810.019.939.8410.3910.5410.49
RMSE15.6915.6315.7414.7913.9013.3313.2713.1313.8414.0213.83
R0.660.680.690.720.740.760.760.760.740.740.74
TotalBias−2.68−1.71−1.49−1.15−0.86−0.57−0.48−0.50−0.68−0.69−1.21
MAB10.319.679.298.968.918.958.979.079.349.6110.05
RMSE13.4412.8912.5712.2612.2112.2712.3112.3712.5812.7413.13
R0.730.760.760.760.760.760.760.760.760.760.74
Table 2. Ranges (largest bias minus smallest bias) in median Biases or average Biases (shown in Figure 12 and Figure 13) from CERES-retrieved hourly Rs at all ground-based sites for each hour under different cloud cover and AOD conditions. Variation (the standard deviation of the biases) in the median Biases or average Biases shown in parentheses. Unit: %.
Table 2. Ranges (largest bias minus smallest bias) in median Biases or average Biases (shown in Figure 12 and Figure 13) from CERES-retrieved hourly Rs at all ground-based sites for each hour under different cloud cover and AOD conditions. Variation (the standard deviation of the biases) in the median Biases or average Biases shown in parentheses. Unit: %.
Impact
Factor
Range
(Variation)
7:008:009:0010:0011:0012:0013:0014:0015:0016:0017:00
cloud
cover
Median
Biases
5.38
(2.43)
4.83
(2.27)
4.50
(2.13)
4.51
(2.01)
4.49
(1.99)
4.39
(1.98)
4.46
(2.04)
4.39
(2.02)
4.43
(2.14)
4.68
(2.31)
5.09
(2.42)
Average
Biases
3.68
(1.48)
3.87
(1.46)
3.12
(1.18)
2.25
(0.94)
2.12
(0.94)
1.97
(0.86)
2.11
(0.92)
2.31
(0.94)
3.26
(1.27)
4.01
(1.57)
4.03
(1.75)
AODMedian
Biases
1.90
(0.74)
0.36
(0.14)
0.50
(0.19)
0.38
(0.16)
0.31
(0.12)
0.77
(0.30)
1.04
(0.39)
1.06
(0.41)
1.50
(0.60)
2.17
(0.93)
2.09
(0.86)
Average
Biases
2.52
(1.01)
1.79
(0.72)
2.27
(0.89)
2.19
(0.86)
1.80
(0.70)
1.23
(0.45)
0.79
(0.29)
0.57
(0.22)
0.46
(0.17)
2.08
(0.80)
2.10
(0.83)
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

Lu, L.; Ma, Q. Diurnal Cycle in Surface Incident Solar Radiation Characterized by CERES Satellite Retrieval. Remote Sens. 2023, 15, 3217. https://0-doi-org.brum.beds.ac.uk/10.3390/rs15133217

AMA Style

Lu L, Ma Q. Diurnal Cycle in Surface Incident Solar Radiation Characterized by CERES Satellite Retrieval. Remote Sensing. 2023; 15(13):3217. https://0-doi-org.brum.beds.ac.uk/10.3390/rs15133217

Chicago/Turabian Style

Lu, Lu, and Qian Ma. 2023. "Diurnal Cycle in Surface Incident Solar Radiation Characterized by CERES Satellite Retrieval" Remote Sensing 15, no. 13: 3217. https://0-doi-org.brum.beds.ac.uk/10.3390/rs15133217

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