Several studies have shown that changes in incoming solar radiation and variations of the diffuse fraction can significantly modify the vegetation carbon uptake. Hence, monitoring the incoming solar radiation at large scale and with high temporal frequency is crucial for this reason along with many others. The European Organization for the Exploitation of Meteorological Satellites (EUMETSAT) Satellite Application Facility for Land Surface Analysis (LSA SAF) has operationally disseminated in near real time estimates of the downwelling shortwave radiation at the surface since 2005. This product is derived from observations provided by the SEVIRI instrument onboard the Meteosat Second Generation series of geostationary satellites, which covers Europe, Africa, the Middle East, and part of South America. However, near real time generation of the diffuse fraction at the surface level has only recently been initiated. The main difficulty towards achieving this goal was the general lack of accurate information on the aerosol particles in the atmosphere. This limitation is less important nowadays thanks to the improvements in atmospheric numerical models. This study presents an upgrade of the LSA SAF operational retrieval method, which provides the simultaneous estimation of the incoming solar radiation and its diffuse fraction from the satellite every 15 min. The upgrade includes a comprehensive representation of the influence of aerosols based on physical approximations of the radiative transfer within an atmosphere-surface associated medium. This article explains the retrieval method, discusses its limitations and differences with the previous method, and details the characteristics of the output products. A companion article will focus on the evaluation of the products against independent measurements of solar radiation. Finally, the access to the source code is provided through an open access platform in order to share the expertise on the satellite retrieval of this variable with the community.

The present paper is devoted to the quality assessment of two global land surface albedo products developed by Meteo France in the frame of the Copernicus Climate Change Service (C3S) and the LSA-SAF (Satellite Application Facility on Land Surface Analysis), herein called, respectively, VGT (VeGeTation) (the C3Sv1 dataset, derived from VGT sensors onboard Satellites for the Observation of the Earth, also called SPOT) and ETAL (European polar system Ten-day surface ALbedo, derived from Advanced Very High Resolution Radiometers (AVHRR) onboard METeorological OPerational (METOP) satellites). The evaluation study inter-compared these products with measurements at 33 ground stations and two independent operational products, MTAL-R/NRT (Meteosat second generation Ten-day ALbedo Reprocessed/Near Real-Time) and MODIS (MODerate resolution Imaging Spectroradiometer), over two distinct four-year periods. In accordance with the prescription from the Land Product Validation group of the joint Committee on Earth Observation Satellites (LPV/CEOS), the evaluation was addressed per land cover; furthermore, two albedo regimes were considered throughout the evaluation to distinguish between high (over 0.15) and low (below 0.15) surface albedo behaviors. First, we show that both VGT and ETAL products agree well with the measurements and the other satellite products at the ground stations. Second, when inter-compared with MODIS, the results are noteworthy for ETAL as opposed to VGT, with 11 out of 13 land cover types passing the Global Climate Observing System (GCOS) requirements for more than 80% of the sites for albedo values less than 0.15 (compared with none for VGT) and 10 out of 14 land cover types passing the GCOS requirements for more than 50% of the sites for albedo values greater than 0.15 (compared with 5 for VGT). Finally, a pixel-by-pixel analysis reveals that VGT overestimates the surface albedo as compared with MODIS by about 0.02 in absolute value for albedo values less than 0.15 and by about 22% in relative value for albedo values greater than 0.15. The root-mean-square-deviation (RMSD) in absolute value is about 0.015 for albedo values less than 0.15 and 51.5% in relative value for albedo values greater than 0.15. In contrast, the bias for ETAL when compared with MODIS remains very small. Over the four-year period, ETAL overestimates the surface albedo as compared with MODIS by 0.001 in absolute value for the regime of surface albedo less than 0.15 and by about 5.8% in relative value for albedo values greater than 0.15. The RMSD in absolute value is about 0.014 for albedo values less than 0.15 and 19.4% in relative value for albedo values greater than 0.15. Assuming that the MODIS product is a good reference, a relative bias of around 6% can be judged satisfactory for ETAL surface albedo. The lower performance of the VGT (C3Sv1) product is currently the subject of investigation. Work is ongoing to upgrade it further towards the final C3S product.

High frequency knowledge of the spatio-temporal distribution of the downwelling surface shortwave flux (DSSF) and its diffuse fraction (fd) at the surface is nowadays essential for understanding climate processes at the surface–atmosphere interface, plant photosynthesis and carbon cycle, and for the solar energy sector. The European Organization for the Exploitation of Meteorological Satellites (EUMETSAT) Satellite Application Facility for Land Surface Analysis operationally delivers estimation of the MDSSFTD (MSG Downwelling Surface Short-wave radiation Fluxes—Total and Diffuse fraction) product with an operational status since the year 2019. The method for retrieval was presented in a companion paper. Part 2 now focuses on the evaluation of the MDSSFTD algorithm and presents a comparison of the corresponding outputs, i.e., total DSSF and diffuse fraction (fd) components, against in situ measurements acquired at four Baseline Surface Radiation Network (BSRN) stations over a seven-month period. The validation is performed on an instantaneous basis. We show that the satellite estimates of DSSF and fd meet the target requirements defined by the user community for all-sky (clear and cloudy) conditions. For DSSF, the requirements are 20 Wm⁻² for DSSF < 200 Wm⁻², and 10% for DSSF ≥ 200 Wm⁻². The mean bias error (MBE) and relative mean bias error (rMBE) compared to the ground measurements are 3.618 Wm⁻² and 0.252%, respectively. For fd, the requirements are 0.1 for fd < 0.5, and 20% for fd ≥ 0.5. The MBE and rMBE compared to the ground measurements are −0.044% and −17.699%, respectively. The study also provides a separate analysis of the product performances for clear sky and cloudy sky conditions. The importance of representing the cloud–aerosol radiative coupling in the MDSSFTD method is discussed. Finally, it is concluded that the quality of the aerosol optical depth (AOD) forecasts currently available is accurate enough to obtain reliable diffuse solar flux estimates. This quality of AOD forecasts was still a limitation a few years ago.