Dual-View Clear-Sky Top-of-Atmosphere Albedo From Meteosat Second Generation Satellites

Geostationary observations offer the unique opportunity to resolve the diurnal cycle of the Earth's Radiation Budget at the top of the atmosphere (TOA), crucial for climate-change studies. However, a drawback of the continuous temporal coverage of the geostationary orbit is the fixed viewing geometry. As a consequence, errors in the angular distribution models (ADMs) used in the radiance-to-flux conversion process can result in systematic errors of the estimated radiative fluxes. In this work, focusing on clear-sky reflected TOA observations, we compare the overlapping views from Meteosat Second Generation satellites at 0{\deg} and 41.5{\deg}E longitude which enable a quantification of viewing-angle-dependent differences. Using data derived from the Spinning Enhanced Visible and InfraRed Imager (SEVIRI), we identify some of the main sources of discrepancies, and show that they can be significantly reduced at the level of one month. This is achieved, separately for each satellite, via a masking procedure followed by an empirical fit at the pixel-level that takes into account all the clear-sky data from that satellite, calculated separately per timeslot of the day, over the month of November 2016. The method is then applied to each month of 2017, and gives a quadratic mean of the albedo root-mean squared difference over the dual-view region which is comparable from month to month, with a 2017 average value of 0.01. Sources of discrepancies include the difficulty to estimate the flux over the sunglint ocean region close to the limbs, the absence of dedicated angular distribution models for the aerosol-over-ocean case in the data processing, and the existence of an observer-dependent diurnal-asymmetry artefact affecting the clear-sky-albedo dependence on the solar zenith angle particularly over land areas.


Introduction
The Earth's Radiation Budget (ERB) at the top of the atmosphere (TOA) is a crucial observable and one of the Essential Climate Variables [1] defined by the Global Climate Observing System (GCOS) for climate-change studies. For the Earth system, it describes the energy balance between what comes in from the Sun and what leaves the Earth, both as reflected solar (shortwave) radiation and as outgoing thermal (longwave) emission. The global warming currently ongoing with unprecedented speed since the mid-20th century is a direct consequence of an overall imbalance in the ERB, predominantly caused by human activities from the industrial revolution to the present [2].
Characterising as best as possible the state of the ERB can only be done from space, and for that we crucially need observer-independent quantities to be retrieved, knowing that satellite instruments can only directly provide observer-dependent radiances L in a given direction at a given time. 1 Going from the measured observer-dependent radiances to such observer-independent quantities involves two important steps. The first is called unfiltering and compensates for the spectral-dependence of the engineered instrument [3]. The second is the inversion-also known as angular conversion -which, from that unfiltered radiance in a given direction then associates an estimated observer-independent irradiance appropriate for the observed scene [4], typically via a collection of empirical angular distribution models (ADMs), such as the Clouds and the Earth's radiant Energy System (CERES) Tropical Rainfall Measurement Mission (TRMM) ADMs [5]. Formally, the irradiance or 'flux' F(θ ) leaving an imaginary surface element at the top of the atmosphere (in units of W m −2 ) is the integral, over all the outgoing solid angles in a hemisphere, of the individual radiances L(θ , θ vz , φ rel ) (in units of W m −2 sr −1 ) leaving that surface element [6]: where θ is the solar zenith angle, θ vz is the viewing zenith angle, and φ rel is the relative azimuth angle. As well-known, this relation simply reduces to F(θ ) = πL(θ ) when the flux is isotropic (Lambertian). ADMs are then typically introduced at this point [6], and provide anisotropic factors which compare an equivalent Lambertian flux πL(θ , θ vz , φ rel ) to the actual flux F(θ ) given in Eq. (1): With a proper set of empirical ADMs, Eq. (3) can then be inverted to derive the flux from the radiance. Finally, the albedo (a pure number between 0 and 1) is defined as the ratio of the reflected and incoming solar fluxes: where F incoming (θ ) = E cos θ ; E being the solar constant [7] corrected for the Sun-Earth distance [8].
For climate models, it is moreover important to precisely monitor the internal structures of the ERB over the globe and throughout the day, seeing that these govern important aspects of the climate on our planet; in particular, its diurnal cycle, and especially the formation of clouds during the day, is a key component of the tropical climate [9]. Most dedicated ERB observations thus far have been made from Low-Earth-Orbit (LEO) satellites. Chiefly among these are Sun-synchronous polar orbiters, such as NASA's Terra and Aqua satellites, which have been observing over the last two decades with the CERES instrument [10] each location on Earth-though not more than twice per day outside of the polar regions. The joint EUMETSAT/ESA geostationary (GEO) Meteosat Second Generation (MSG) satellites are actually in a unique position to complement these observations and observe the diurnal cycle, and this for two reasons. The first is that their geostationary orbits offer the advantage of an excellent temporal sampling of the observed locations throughout the day. The second is that they embark both the multispectral Spinning Enhanced Visible and InfraRed Imager (SEVIRI) with 12 spectral channels at a nadir resolution of 3 km (1 km for the high-resolution visible channel) [11], and the broadband Geostationary Earth Radiation Budget (GERB) instrument with nadir resolution of 50 km [12]. Every 15 minutes both these instruments provide observations of the full Earth disk seen from the satellite viewpoint. They do have the disadvantage of fixed viewing angles though. Compared to LEO satellites, this makes the measurements of the broadband TOA radiative fluxes more sensitive to angular-dependent errors. It is therefore particularly important to address any remaining observer-dependent systematics in flux or albedo data, knowing that such unphysical artefacts would introduce errors, e.g. in the Earth's Radiation Budget determination.
Interestingly, we now actually have the opportunity to cross-check retrieved GEO products. Indeed, since the end of the year 2016, there have been two MSG satellites at different longitudes and with overlapping scenes. One of them is at 0°longitude (MSG-3, replaced by MSG-4 on 2018/02/20), and the other is at 41.5°E longitude, providing the 'Indian Ocean Data Coverage' (MSG-1, which should be replaced by MSG-2 at 45.5°E longitude in 2022; see e.g. Ref. [13]). Taking advantage of this dual view, our aim with this work is then to compare, quantify the discrepancies and try to address them in order to make GEO-derived products as useful as possible. For simplicity the focus will only be on the clear-sky top-of-atmosphere broadband shortwave albedo, and we will often simply refer to it as the 'albedo'. In the following, we are only going to use and consider pure SEVIRI synthetic products: the so-called 'GERB-like' GL-SEV products [14]. 2 Involving a narrowband-to-broadband procedure [15] and relying on the CERES TRMM ADMs [5] for the angular-conversion process, these products consist of HDF5 files available every 15 minutes (96 per day), containing 3 × 3 SEVIRI-pixel image layers (9 km × 9 km resolution at nadir).
This paper is organised as follows. In Section 2, we present the basic methodology that we use to compare the results from the different satellites. We then show in Section 3 how the monthly consistency between the two views can be improved for the month of November 2016, via masks and an empirical-fit procedure. Finally, Section 4 shows that similar results can be obtained for the entire year 2017, before we conclude in Section 5.

Dual-view comparison method
The same treatment is applied independently to GL-SEV images from MSG-3 (SEV3) and MSG-1 (SEV1), and is done one month at a time. For that month, and for each individual HHMM timeslot among the 15-minute timeslots, our basic treatment is the following: • for each day at that timeslot, we first derive the instantaneous clear-sky TOA albedo image after applying the appropriate masks-e.g. keeping only those pixels for which the cloud-cover layer is zero (as discussed in the following, we are going to apply a number of extra conditions); • then, similarly to what is done for monthly hourly products 3 , we use all these images to calculate the monthly "representative albedo image" at that timeslot: in this work, we considered both the mean and the median, and unless otherwise stated we will show results obtained using the median albedo (robust statistics); • in a common 0.5°× 0.5°latitude-longitude grid, we can then calculate the grid-box difference for that specific timeslot; in other words, each pixel location x will store the albedo difference which would be 0 for all pixel if the calibration of SEV1 and SEV3 was the same, if there were no narrowband-to-broadband nor angular-dependent errors, and if the impact of having different clear-sky atmospheric paths from a given scene to each of the two satellites can be reduced to a strictly angular issue -an assumption already implicit in the radiance-to-flux conversion.
Once this treatment is done, we can proceed and apply statistics over the whole set of 96 'SEV3 − SEV1' images: for each pixel x, using all of the (t = 1, ..., n t x ) timeslots for which we have 2 We use what are, at the time of writing, the very latest available products covering the dual-view period: GL-SEV HR V003. These can be obtained via the https://gerb.oma.be website. 3 For monthly hourly products, see for instance the CERES SYN1deg [16] and the CM SAF MMDC [17] and [18] products. non-masked data at that location, we calculate the root-mean-squared difference or root-mean-squared deviation (i.e. the square root of the mean squared albedo difference): the bias (i.e. the mean albedo difference): and finally the bias-corrected standard deviation σ, from the corresponding variance:

Raw GL-SEV products
Since our aim is to assess and compare the retrieved products from the two different MSG satellites, here we use the fluxes exactly as they are given in the GL-SEV products. Note in particular that, as part of the data processing of these products, a shortwave flux over the sunglint region is provided [19], and that the case of 'aerosol over ocean' pixels is treated there as 'clear-sky ocean' with the use of identical CERES TRMM ADMs. Figure 1 (a) shows the root-mean-squared-difference result that we obtain from our comparison when we naively use the shortwave fluxes whenever they are defined in the GL-SEV products. 4 Among the most salient features that we can see is that there are issues related to unidentified clouds, which actually occur mainly at large viewing and solar zenith angles. What is then probably the most striking difference are the large extended regions (dark blue) in the bottom left and the bottom right of the image, where the discrepancies can be larger than 0.10. These can actually be linked to the sunglint regions (when the viewing direction coincides with the specular sunlight reflection off the ocean surface) getting close to either of the limbs of the overlapping area: the left one corresponds to MSG-1 and the right one, to MSG-3. Other areas over the ocean with discrepancies which can reach ∼ 0.05 are also visible close to the west and east coasts of Africa (blue), and are actually due to the presence of aerosols. These are serious discrepancies, given how low the albedo over clear-sky ocean typically is (∼ 0.06 under overhead solar conditions). Since the ocean covers about 70% of the Earth, it is particularly important to correct these. Finally, further discrepancies are found over land.

Masks
We can now try to improve the consistency between the albedo values retrieved from the two satellites by first applying a number of masks and cuts.
First of all, what is frequently done in the literature is to cut on the zenith angles as the quality is expected to suffer from increasing errors at large values; see e.g. Refs. [20,21]. The exact threshold may slightly differ from paper to paper; here we use θ , θ vz < 70°. 5 After doing that, we can see that things are already improved (see Figure 1 (b)), but notice that the discrepancies which we identified for the 4 The albedo is set to 1 where it exceeds this value. Such pixels are not masked, to help identify and address potential issues. 5 For consistency with an aerosol-optical-depth mask applied in the following, as there is no retrieval at larger zenith angles. , but now adding a sunglint mask, requiring the tilt angle to be > 25°over ocean pixels. (d, bottom right) Same as (c), now adding an aerosol mask over ocean, requiring that the aerosol optical depth at 0.6 µm is < 0.1, and extending the sunglint mask to retain only tilt angles > 40°over ocean pixels.
sunglint and the aerosol areas in particular actually remain of the same order, respectively > 0.10 and ∼ 0.05. Since very large differences coming from the shortwave flux over the sunglint regions remain, we then remove the modelled fluxes provided in those regions by masking all the ocean pixels for which the tilt angle is smaller than 25°. Doing so gives Figure 1 (c).
After applying these different masks, the discrepancies related to aerosols over the ocean are still clearly visible and still of the order of ∼ 0.05. Due to the lack of dedicated ADMs for that case 6 , we are going to mask these aerosol-loaded regions over ocean pixels. 7 There is an aerosol-optical-depth retrieval over the ocean in the GERB/GL-SEV processing [22]; we can therefore use the aerosol-optical-depth observations at 0.6 µm, and only keep those pixels for which this optical depth is smaller than 0.1. For consistency with the aerosol-retrieval algorithm and to avoid introducing artefacts in the images, the sunglint mask used in the previous step is then enlarged to only keep tilt angles greater than 40°over ocean pixels. This is our final result using masks; it is shown in Figure 1 (d). 6 Note that the CERES team proposes a method to account for aerosols in the clear-sky ocean ADM; see Ref [5]. However, this correction is not actually applied in the GL-SEV processing. 7 The presence of aerosols results in a quite diffuse reflection; very different from the strong specular reflection in the clear-sky ocean case. The use of clear-sky ocean ADMs that do not take aerosols into account is therefore particularly problematic.
Its decomposition in terms of standard deviation and bias is given in Figure 2, from which one can see that the standard deviation typically contributes much more to the RMSD than does the bias. Now, if we furthermore compare the statistics over the entire image obtained when naively using the fluxes as provided in the GL-SEV products (Figure 1 (a)) to the ones that we obtained after the masking procedures we just described (Figure 1 (d)), we obtain the results given in Table 1. These summary statistics provide a single value after looping over all the (x = 1, . . . , N x ) non-masked pixels in the corresponding image. More precisely, denoting the RMSD, σ, and bias of a given pixel via a · x subscript, the quantities given in that table are the square root of the average mean squared difference: the square root of the average variance: and finally, the mean absolute bias (or mean absolute difference): In particular, we see in Table 1 that the quadratic mean of the RMSD over all pixels, Eq. (9), is significantly reduced, going from 0.055 to 0.014. The consistency is much improved after this masking procedure, but we stress once again that we actually had to mask the aerosol over ocean case in order to obtain these results, and that having access to correct fluxes also in such cases is clearly a high priority; see e.g. Ref [23].

Diurnal-asymmetry artefact
We now move on to discuss a well-known issue with geostationary data: due to the fixed viewing geometry from GEO, there can be a systematic diurnal-asymmetry artefact in GEO shortwave fluxes [24]; see also Ref. [25].
In the majority of cases however, no asymmetry in the clear-sky top-of-atmosphere albedo is actually expected between the (local time) morning and afternoon. That is, at least if the surface properties do not change over time [26]. Indeed, the albedo over land essentially evolves as a function of the solar zenith angle, which is symmetrical with respect to the local noon [27,28]: with parameters a 60 = a(θ = 60°) and d, a dimensionless parameter tied to how exactly the albedo changes with cos θ . Equation (12) is extensively found in the literature, both for TOA and surface albedo studies as a function of the solar zenith angle; see Ref. [29] for a very detailed study which also reviews the literature. This functional form was first introduced in Ref. [27] for the case of infinite canopies. It was then used in Ref. [28] for all types of land surfaces. As shown in Figure 3, one can check against the CERES ADMs that the same functional form can be used over ocean pixels (rms of residuals: 0.005); this is using the average wind speed model. Knowing that the albedo, just like the flux, is expected to be only a function of the solar zenith angle at any single time and knowing the general functional form of this dependence, we can then take benefit of the GEO temporal sampling and consider all the observations at different times in order to obtain a more robust and consistent result at the level of one month. Following this method will provide a corrected monthly mean (level-3) GL-SEV product; it does not directly correct the instantaneous level (level-2). Note that this assumes that the surface properties and the atmospheric transmission affecting the clear-sky top-of-atmosphere radiation do not significantly change over that period of one month.
In Figure 4 (left), we show an example of the diurnal-asymmetry artefact in a case of a rarely cloudy sample pixel taken in the Sahara (18.62°latitude, 25.52°longitude), for which we can therefore get a very detailed view of the clear-sky-TOA-albedo temporal evolution throughout the day. This is seen from MSG-3 and we can clearly see that there are in fact two branches for similar values of the cosine of the solar zenith angle, the higher branch in this case corresponding to the morning, and the afternoon branch being the lower one. In this figure, each pair of colour and small symbols corresponds to the instantaneous clear-sky-TOA albedo for different days at the same timeslot. For each of these, we also overlay a summary representative (mean or median) value for that timeslot as a larger symbol  Figure 4. Diurnal-asymmetry artefact for the same location in the Sahara, seen from MSG-3 (left) and MSG-1 (right); see main text. For each timeslot (i.e. each colour + symbol pair), the large symbol inside a diamond shape shows the median albedo versus the mean cos θ for that timeslot, while the corresponding smaller symbols denote the instantaneous data used to compute those. inside a diamond shape-unless otherwise stated all such results shown here were obtained using the median albedo versus the mean cosine of the solar zenith angle at the corresponding timeslot.
In Figure 4 (right), we show the same case, but now seen from MSG-1. Notice how the two branches are then actually swapped: the higher branch corresponding to the afternoon in this case, and the morning branch being the lower one. This complete contradiction at the qualitative level highlights the fact that this asymmetry cannot be of physical origin. 8 The problem is rather related to imperfections in the ADMs used to convert radiances into flux, which, due to the fixed view of GEO observations, are then turned into systematic errors. This well-known situation is not specific to this given pixel, but truly an ubiquitous observation.
Since the existence of these branches is not physical 9 , one thing that we can do is to fit the albedo as a function of the solar zenith angle according to the functional form of Eq. (12), and then use that to calculate the representative mean or median albedo, before projecting. These fits are done by non-linear least-squares optimisations taking into account all the available information over the n t x timeslots: 10 min a 60 ,d n tx ∑ t=1 a t − a 60 1 + d 1 + 2dµ t 2 independently for each GEO pixel; (13) where for each timeslot t, µ t = cos θ ,t and a t are respectively the instantaneous values for the cosine of the solar zenith angle and for the albedo derived from the GL-SEV products. We do stress that the objective here is not to provide a predictive model of what the albedo might be, based on the scene; rather, this is merely an empirical fit, using only all the data for a given pixel over a complete month.
In the case of our example, the results are shown in Figure 5. Figure 6 shows another sample pixel, this time in the Horn of Africa and with remaining unidentified clouds from the MSG-3 viewpoint, especially in the morning (latitude 7.45°, longitude 48.70°). Note that rather similar results are obtained 8 While the branches themselves cannot be physical, one can notice the presence of overimposed instantaneous patterns. In contrast to the branches, those are actually consistent between the two satellites at corresponding times and likely of physical origin. They are in fact already present in SEVIRI narrowband radiances. See Appendix A. 9 Further note that any actual observer-independent physical effect with a significant impact on the albedo should remain clearly visible even with swapped branches from MSG-3 to MSG-1. If there were dew [26] in the early morning for instance, the shape of the two morning branches would then be similarly skewed upwards and the apparent symmetry seen here when swapping the morning and afternoon branches would be lost. 10 This is done using the NLopt library [30]. Having tested several of the available algorithms, we find that the derivative-free NEWUOA [31] and gradient-based SLSQP [32] algorithms are especially reliable, fast, and accurate for the problem at hand. when using the mean instead of the median albedo, but the fit is then necessarily more sensitive to the presence of unidentified clouds, systematically leading to a larger a 60 and therefore larger albedo values for as a function of cos(θ ) in those cases. Another interesting thing to notice is that, as this second example pixel is close to being at the same longitude as the MSG-1 satellite, seen from there, the branches are then not so spread but close to being indistinguishable from one another. In Figure 7, we show the root-mean square difference in clear-sky albedo retrieved from SEV3 and SEV1, where we both use all the masks that were introduced in Section 3.2 and the imposed angular consistency discussed here before projecting onto the latitude-longitude grid. We request at least three timeslots with defined albedo in order to make the fit and otherwise mask the corresponding pixel. 11 The corresponding summary statistics over all defined pixels are given in Table 2. In particular, we can see that the quadratic mean of the RMSD over the whole image is now as low as 0.008 in the overlapping region, and that no large discrepancy remains. Within the RMSD, the standard deviation shown in Figure 8 (left) is much decreased as a result of the fit (compare to Figure 2 (left)). Notice that the bias shown in Figure 8 (right), remains small, even though the average absolute bias is slightly larger overall after applying the fit: over the northern part of the African continent and the Arabian peninsula in particular we can see a slightly strengthened trend with a positive bias (light red colour) to the west of 20.75°E (halfway between 0°and 41.5°E) and a negative bias to the east (light blue colour).  Note that, alternatively, one could instead already proceed with a simple binning in cos θ . An approach of this type was actually adopted in a preliminary study with just two bins [33], and the results obtained for the same month ( RMSD 2 = 0.012, σ 2 = 0.008, and |bias| = 0.005) are comparable to those we report here in Table 2, giving confidence that the general method itself is quite robust.

Merged overhead albedo
Using false colour rendering, Figure 9 shows the merged overhead albedo a 0 in the different cases presented in Figure 1. The colours were chosen so that typical values over the ocean (a 0 ∼ 0.06) would be blue; vegetation (a 0 ∼ 0.12), green; dark and bright deserts (a 0 ∼ 0.24 and ∼ 0.36), brown and yellow. Aerosol-loaded regions over the ocean (with typical a 0 ∼ 0.08) appear in purple, and can be seen on panel (c) close to India and the Gulf of Guinea in particular. For each pixel, we used the fitted parameters obtained when imposing the functional form discussed in the previous section. For an overhead Sun, a 0 can be simply calculated from Whenever the albedo is defined for both SEV3 and SEV1 for a given latitude-longitude pixel, an average is taken. 12 The original sharp discontinuities between satellite views are significantly decreased as the masks described in Sec. 3.2 (zenith-angles (b), sunglint (c), and aerosols (d)) are being applied.

Visualising the diurnal-asymmetry artefact from the viewpoint of each satellite
Before applying the method described above to a full year of data, let us take a closer look at the spatial distribution of the observer-dependent diurnal-asymmetry artefact. It is indeed interesting to be able to visualise it, not only one pixel at a time, but for all pixels at once, from the viewpoint of each satellite. We are going to fit the morning and afternoon branches separately and compare them. 13 This will highlight where, according to each satellite, morning albedos appear larger than afternoon albedos and vice versa. Though the shapes of the individual branches are not strictly given by Eq. (12), this fitting procedure will provide instructive results. Figure 10 shows at pixel-level the result of such fits for both SEV3 and SEV1, for the rarely cloudy Sahara location discussed in Sec. 3.3. To keep as many pixels as possible while making sure that each branch is well-sampled and not overly sensitive to the curvature change close to the tip around the local noon, we find that a good trade-off is to only keep pixels with non-masked data spanning at least 12 different timeslots for each branch.  Figure 10. Same case as Figure 5, with each branch fitted to enable the visualisation shown in Figure 11.  Figure 11. Visualisation of the diurnal asymmetry artefact. This compares for each pixel the a 60 fit parameters (albedo corresponding to θ = 60°) obtained for each of the branches (a 60,AM − a 60,PM ). The asymmetry, seen from MSG-3 (left) or MSG-1 (right), is obviously strongly observer dependent. Figure 11 then compares for each pixel the a 60 parameter obtained when fitting the morning branch (AM) to the one obtained for the afternoon branch (PM). The gaps in ocean pixels are largely due to the sunglint and aerosol masks on the ocean, which prevent meeting our fit-quality requirements. As we focus on the common region seen from both satellites and centred on Africa, we see that the two satellites strongly disagree on which branch would be larger than the other virtually everywhere. From the point of view of MSG-3 (left panel), compared to the satellite position (centre of the image), the albedo tends to be systematically larger during the morning than in the afternoon for pixels on the right of the image (i.e greener on the right), while the situation tends to be reversed to the left of the image, where the afternoon albedo is then larger than the morning albedo (i.e redder on the left). In other words, the retrieved albedos for pixels to the east of the satellite tend to be brighter in the morning, and those to the west, brighter in the afternoon, as seen also in Figure 12. Seen from MSG-1 (right panel), there seems to be a similar left/right trend with respect to the satellite position at the centre of the image. This again highlights that this diurnal asymmetry is a viewing-angle issue. Note that the rare exceptions where both satellites agree appear to be mostly linked to overimposed cloud-identification issues, such as in Gabon or near the Horn of Africa in the morning (see also Figure 12). 15 Without correcting this artefact, SEV3 would tell that most of Africa has a larger albedo in the morning, while SEV1 would tell that it is larger in the afternoon instead. 16

One year
Finally, we now simply apply the treatment discussed in Secs. 3.2 and 3.3 to each of the months of the year 2017. Again, this is done independently for MSG-3 and MSG-1. As before, we compare the results and compute the summary statistics; these are given in Table 3. As expected, the results are comparable to those given in Tables 1 and 2 and the consistency is systematically much improved. For completeness, note that in the 2017 data which was available for us to compute those, there were six entire days with missing GL-SEV data from MSG-1: from January 17th to January 22th (MSG-1 SEVIRI decontamination), and on November 12th. Further averaging these monthly results for the quadratic mean of the RMSD, the quadratic mean of the standard deviation, and the mean absolute bias over all the months of 2017 in the masks-only case respectively gives 0.017, 0.015, and 0.005. When empirical fits are also used, one then respectively obtains 0.010, 0.006, and 0.005.

Conclusions
In this paper, we have had a closer look at the clear-sky top-of-atmosphere broadband shortwave albedo retrieved from MSG geostationary satellites positioned at 0°and 41.5°E longitude, and compared the results obtained over the very substantial overlapping region observed by both.
This was first done for the month of November 2016. We could identify discrepancies which exist at present in the GL-SEV products. Among the main sources of discrepancies, some in particular were identified to be due to the fluxes retrieved in the sunglint regions over ocean pixels close to the limbs, where the RMSD can locally well exceed values as large as 0.10; another important source of issues being related to aerosol-loaded regions, where the RMSD can be as large as ∼ 0.05. For the aerosol, the observer-dependent differences can be tied to the lack of a dedicated ADM processing. One of the key points stressed in this study is therefore that ADMs appropriate for the case of 'aerosol over ocean', are direly needed. We then aimed at improving the situation. This was achieved first by using masks to remove pixels for which it appears at present difficult to obtain a sufficiently reliable flux, and then by using an empirical physical fit of the expected general dependence of the albedo with the solar zenith angle to apply level-3 corrections at the level of one month-taking advantage of the high temporal resolution characteristic of geostationary satellites. The imposed angular consistency helped address what clearly turned out to be an unphysical observer-dependent diurnal asymmetry artefact, which we also quite extensively discussed. Note that at no point was the additional information coming from the dual view used to correct the results from either of the satellites. They both remained treated independently throughout. These led to sizeable improvements in terms of consistency.
The same approach was then applied over one full year, for all the months of 2017, and we similarly obtained that the agreement between the retrieved albedos from each satellite at the level of one month is much improved with this method. In average, over the twelve months of 2017, the quadratic mean of the albedo root-mean squared difference indeed goes down to 0.01, while the corresponding average quadratic mean of the standard deviation and average mean absolute bias are 0.006 and 0.005 respectively.
Funding: This research was funded by the Solar-Terrestrial Centre of Excellence (STCE). The STCE is a collaboration between the Royal Meteorological Institute of Belgium (RMIB), the Royal Observatory of Belgium (ROB), and the Belgian Institute for Space Aeronomy (BISA).
Acknowledgments: It is our pleasure to thank Christine Aebi for a careful reading of and various comments on the manuscript.

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

Abbreviations
The following abbreviations are used in this manuscript:

Appendix A
Here, we show that the patterns visible in Figure 4 in the instantaneous albedo data at each timeslot are already present in the broadband GL-SEV radiances. To remove any influence from the ADMs, Figure A1 shows this time the pseudo-albedo as a function of the day of the month for the 10:00 UTC timeslot. Since this is independent of the ADMs, there is no angular correction applied at all for the surface type with respect to an ideal Lambertian surface; it is therefore expected that the exact pseudo-albedo values from MSG-3 or MSG-1 cannot be directly meaningfully compared. What matters however is that a similar pattern/shape is clearly seen from each satellite over that timeslot.
Such a pattern is actually even visible in the SEVIRI narrowband radiance data and the two independent satellites again appear to agree; these fluctuations are therefore very likely physical. Figure A2 shows the VIS0.8 and VIS0.6 channels, and VIS0.8 in particular is quite clearly reminiscent of the pattern in broadband GL-SEV synthetic data.
Let us stress that the presence of such patterns is a general observation, not limited to this specific example case (though the pattern shape itself of course changes for different pixels and periods of time). They can sometimes surprisingly remain essentially unchanged at different times of the day, as is the case for the current Sahara example (seemingly same pattern repeated over different timeslots).  30  29  28  27  26  25  24  23  22  21  20  19  18  17  16  15  14  13  12  11  10  9  8  7  6  5  4  3  2  1 0.38 0.36 0.34 Figure A1. Consistent patterns in instantaneous data seen from MSG-3 (squares) and MSG-1 (dots). This shows the pseudo-albedo πL/F incoming , calculated directly from the broadband GL-SEV radiances at 10:00 UTC, for every day of the month of 2016/11 in the Sahara sample case shown in Figure 4. Note that in this case the pattern is inverted along the x-axis when shown as a function of the day of the month instead of cos θ (i.e. cos θ values are larger at the beginning of the month, smaller at the end).  Figure A1, but this time using directly SEVIRI narrowband radiance data for the corresponding pixels at 10:00 UTC, in the same 3 × 3 (S3) grid as the GL-SEV products.