Terrestrial Planet Optical Phase Curves. I. Direct Measurements of the Earth

NASA’s EPOXI mission used the Deep Impact spacecraft to observe the disk-integrated Earth as an analog to terrestial exoplanets’ appearance. The mission took five 24 hr observations in 2008–2009 at various phase angles (57.°7–86.°4) and ranges (0.11–0.34 au), of which three equatorial (E1, E4, E5) and two polar (P1, North and P2, South). The visible data taken by the HRIV instrument ranges from 0.3 to 1.0 μm, taken trough seven spectral filters that have spectral widths of about 100 nm, and which are centered about 100 nm apart, from 350 to 950 nm. The disk-integrated, 24 hr averaged signal is used in a phase angle analysis. A Lambertian-reflecting, spherical planet model is used to estimate geometric albedo for every observation and wavelength. The geometric albedos range from 0.143 (E1, 950 nm) to 0.353 (P2, 350 nm) and show wavelength dependence. The equatorial observations have similar values, while the polar observations have higher values due to the ice in view. Therefore, equatorial observations can be predicted for other phase angles, but (Earth-like) polar views (with ice) would be underestimated.


Introduction
After almost three decades of exoplanet discoveries, we now know that possibly up to 50% of solar-type stars have a small, terrestrial-type planet in its habitable zone (Bryson et al. 2020). However, even though we know such small, potentially habitable planets are out there, little is known about the properties of their atmospheres, surface, climate, and with that, their habitability. To obtain such information on small planets in the habitable zones of solar-type stars, direct observations of the visible to near-infrared starlight that these planets reflect and/or the thermal radiation that these planets emit should be measured. Except for planets around the nearest stars, which might be imaged by large ground-based telescopes such as the European Extremely Large Telescope, such measurements would require dedicated space telescopes beyond current capabilities.
Space telescopes that could perform these measurements are being proposed and designed with specific interest in exoplanet observations. The Habitable Exoplanet Observatory (Habex) proposed mission focuses on habitable zone Earth-like planets around Sun-like stars. The various instruments would cover a wavelength range from ultraviolet 0.115 μm to near-infrared 1.8 μm (Gaudi et al. 2020). The Large UV/Optical/IR Surveyor (LUVOIR) is another proposed mission that will have the capabilities to observe and characterize a more diverse range of exoplanets. It will have a wavelength range from 0.1 to 2.5 μm (The LUVOIR Team 2019). Both these missions show the interest in exoplanets and the desire to have space telescopes committed to this topic. They are also pushing the technological limits, Habex, by using a free-flying starshade as an external occulter and LUVOIR by using very a large telescope. These missions will only fly by the end of the next decade (∼2040). By then, the Earth's appearance as a model exoplanet (and its dominant factors) must be well understood to compare with results and interpret conditions on spatially unresolved exoplanets.
Even at this time, only a disk-integrated signal can be anticipated. Most Earth observation data comes from near-Earth/Earth orbiting spacecraft (NASA EODIS, ESA GMES; Pfeifer et al. 2012;Aschbacher & Milagro-Pérez 2012;Kansakar & Hossain 2016). Those data are often highly detailed, large volume and recreating Earth exoplanet-like data can be very tedious as it entails data reduction/selection, numerical modeling to reconstruct a simulated view of the whole disk from a distance and/or require special treatments (data gaps, spacecraft maneuvers; Hearty et al. 2009;Robinson et al. 2011;Mettler et al. 2020). Certain physics important for local phenomena that maybe captured by high-resolution instruments is not significant in a full Earth observation. It is essential to simplify modeling as much as possible while retaining the greatest amount of environmental information that can be retrieved unambiguously, a task for which empirical whole-disk observations can provide essential constraints. Data from deep space is needed, capturing the full Earth. Since the distance is not large on astronomical scale, the signal-to-noise ratio is very high. This is advantageous for the analysis compared to real exoplanet observations. The EPOXI mission took such data to allow the reconstruction of Earth as an exoplanet. The EPOXI mission was an extended mission after Deep Impact had completed its mission at 9P/Tempel 1 (Blume 2005). The extended mission included a second comet flyby (Thomas et al. 2013), photometric observations of a number of the then relatively small number of known transiting exoplanets (Christiansen et al. 2010), and observations of Earth as an exoplanet analog . Those data can be used in the future to compare with exoplanets and allow an initial assessment of the habitability of the planet. This article will investigate the change in brightness of Earth with changing phase of illumination comparable to observing a planet move through its orbit, with the outline of this paper is as follows. In Section 2, we describe the EPOXI data that we use for our analysis. In Section 3, a phase angle analysis is performed. Finally, Section 4 contains our conclusions.

The EPOXI Data
Within the Extrasolar Planetary Observation and Characterization part of the EPOXI mission, the Deep Space Impact spacecraft was used to observe the Earth from afar. A detailed description of the targeting, acquisition, selection, processing, and calibration of the data can be found in the paper by Livengood et al. (2011), which serves as the main reference. Other investigations have been done with subsets of the data such as azimuthal mapping of Earth's surface units and polar regions from multicolor light curves and Cowan et al. (2009Cowan et al. ( , 2011; Using surface and atmospheric properties to constrain models for Earth's visible spectrum and rotational light curve (Fujii et al. 2011;Robinson et al. 2011); and using visible colors to empirically categorize Earth among the planets of our Solar System (Crow et al. 2011).
Since the publication of Livengood et al. (2011), the data has been recalibrated. The same techniques are used and the same steps are followed as described by Livengood et al. (2011), except translated from Interactive Data Language to Python code. Note that in this paper, we include polar observations that were not reported in the paper by Livengood et al. (2011) but were used in Cowan et al. (2011). The data set used here is a calibrated second version retrieved from the NASA Planetary Data System. 6 The measurements by the High-Resolution Instrument Visible camera (HRIV) are whole-disk images of the Earth on a 512 × 512 pixel subframe of the whole detector array, taken through seven spectral filters (located in a filter wheel) that have spectral widths of about 100 nm, and which are centered about 100 nm apart, from 350 to 950 nm. Livengood et al. (2011) addressed the convolution of the actual filter functions with the solar spectrum to estimate the upper and lower limits of each filter's bandwidth. Table 1 lists the characteristics of the data, and Table 2 the center wavelengths and the wavelengths that are covered by the seven filters.
Four of the filters (with nominal center wavelengths of 450, 550, 650, and 850 nm; see Table 2) were sampled every 15 minutes over a period of 24 hr, resulting in 97 measurements for each filter. The measurements at the start and the end of the 24 hour period cover the same orientation of the Earth. The other three filters (with nominal center wavelengths of 350, 750, and 950 nm; see Table 2) were sampled every 60 minutes for the same 24 hour period, resulting in 25 measurements for each filter. Here, the measurements at the start and the end also cover the same orientation of the Earth. Note that the data also includes near-infrared spectroscopy of the Earth at 1.0-4.8 μm, that was collected every two hr, but those data will be Note. In the columns (from left to right): the observation, the date (in UTC), the phase angle α, the illuminated fraction a f of the Earth, the subobserver latitude, the subsolar latitude, and the spacecraft-Earth range. a The illuminated fraction f is related to phase angle α through ( ) a = + f 1 cos 1 2 (Cowan et al. 2011). Note. In the columns (from left to right): the observation, and the reflected solar flux for the seven spectral bands. The reflected fluxes have been averaged over the Earth's disk, have been averaged over 24 hr, and have been scaled to a distance of 1 au from the spacecraft and the Sun. The reflected fluxes and their standard deviations are in units of 10 −7 W/m 2 /μm. The wavelength listed in the top row is the nominal center wavelength of each filter. The wavelength below is the average wavelength, and the range below (in nm) is the full-width-at-half-maximum response to the solar spectrum  addressed in a future paper. The sampling strategy was designed to fit the limited spacecraft memory. Due to the changing position of the spacecraft with respect to the Earth, there are different measurement sets. We will analyze three so-called "equatorial" observations, which were taken more or less over Earth's equator. They will be referred to as E1, E4, and E5. We will also analyze two "polar" observations, which were taken at high latitude (60°); one over the North pole (P1) and one over the South pole (P2). The Moon transits the Earth during E4, the Moon's signal has been removed from the affected images. Note that P1 was acquired shortly after the northern spring equinox (i.e., on 2009 March 27-28,), thus slightly more than half of the polar cap was illuminated, while P2 was acquired only a few weeks after the southern spring equinox (i.e., on 2009 October 4-5), thus slightly more than half of the southern polar cap was illuminated. The P1 and P2 measurements have been used by Cowan et al. (2011) in their principal component analysis of the visible spectroscopy of Earth's rotational light curve.
The first Earth image of four of the observational runs is shown in Figure 1, showing the different size of the disk depending on the geocentric distance. The spacecraft-Earth distance is 0.18 au for E1, 0.34 au for E4 and E5 (the latter image is not shown), and 0.11 au for P1 and P2. The Earth covers around half the image for E1 (there, the Earth's diameter is 233 pixels) on a 28,000 × 28,000 km image. For E4 and E5, the Earth is smaller (the Earth's diameter 129 and 126 pixels, respectively) on a 51,000 × 51,000 km and a 52,000 × 52,000 km image, respectively. For P1 and P2, the Earth takes up most of the image (the diameters are 373 and 372 pixels, respectively) on a 17,000 × 17,000 km and a 18,000 × 18,000 km image, respectively.
As can be seen, the images are not sharp as the HRIV instrument was out-of-focus, since the final design has no mechanism to focus the instrument (Klaasen et al. 2008). The lack of focus is no limitation for our application, as an exo-Earth would only be observed as a single pixel. In order to reduce the pixelated imaged to a single pixel value, we integrate the signals across the image. For this, first, the Earth's centroid is determined in every image by finding the highest correlation between the image and a half disk. Then, using aperture photometry, the Earth signal is retrieved, see Figure 2. The aperture radius starts at the Earth's edge and is increased until the collected signal no longer changes. A circular band (annulus) around the aperture is used to collect the local background signal, which is subtracted from the aperture signal. The outer radius is a factor 2 4 times the aperture radius. Table 2 shows the disk-integrated reflected fluxes as derived from the images with the global background level subtracted. The diurnally averaged reflected fluxes as shown in Table 2, can be compared with the values in Table 2 in Livengood et al. 2011. The values for the 950 nm band are slightly lower than the values originally published. The difference is due to calibration improvements (L. M. Feaga 2020, personal communication).
The values of the standard deviation in the data represent the signal variability due to the Earth's rotation over the 24 hour period of each set of observations, rather than the measurement uncertainty ( Figure 3). Precision uncertainty estimated from pixel-to-pixel variations in the background is negligible compared to the variability in the reflected flux that is due to Earth's heterogeneous surface and cloud coverage over the 24 hour period over which the data was acquired. The diskintegrated reflected fluxes with their variability at the various phase angles and filters have been plotted in Figure 3, the details of which will be discussed in Section 3.

Phase Angle Analysis
The flux that is reflected by the Earth and that is measured by the spacecraft, depends on the coverage of the disk at the time of the observation (clouds, ocean, continents, ...) and also, strongly, on the illumination and viewing geometries, and the fraction f of the disk that is illuminated and visible, thus on the phase angle. Phase angle α is defined as the angle between the direction to the Sun and the direction to the spacecraft as measured from the center of the planet (the Sun-Earthspacecraft angle). Note that we ignore the local variations in the angle between the directions to the Sun and the spacecraft that are due to parallaxes. At a phase angle α of 0°, the Earth's disk that is visible for the spacecraft would be fully illuminated by the Sun, while at α = 180°, the night-side of the Earth would be in view. Table 1 shows that for our observations, phase angle α ranges from 57°.7 (E1) to 86°.4 (P2).
The dependence of the signal on the planetary phase angle due to the illuminated fraction f and the illumination and viewing directions of the planet needs to be taken into account when analyzing the observations and when trying to predict observations at other phase angles, since the goal of the analysis of disk-integrated data like ours is to retrieve planet characteristics like surface coverage and weather patterns. If not taken into account correctly, efforts to recover surface properties, rotational variability and weather patterns would obtain false results.
A simple model to predict the phase angle variation of the signal due to the variation in f and the illumination and viewing geometries is a model of a sphere with a Lambertian (thus isotropically), horizontally homogeneous reflecting surface. This model is an approximation because the Earth's atmosphere has a wavelength-dependent influence on the signal such as through scattering and absorption, and including phenomena such as reddened Sun-sets, the surface reflection is wavelength dependent, and the scattering and reflection is usually nonisotropic, thus non-Lambertian. However, as a first order estimation of the Earth's phase angle behavior, this model should allow us to gain insight into the phase angle variation without invoking detailed parameters that might not be representative for exoplanets.
The disk-integrated reflected flux of a Lambertian-reflecting, spherical planet at a phase angle α is given by ( where A s is the surface albedo, r is the planet radius, d is the distance between the planet and the observer, and πF 0 is the incident solar flux measured perpendicular to the direction of incidence. The geometric albedo A G is defined as the flux that is reflected by the planetary body at a phase angle α of 0°divided by the flux that would have been reflected by a white, Lambertian (thus isotropically) reflecting disk of the same cross-sectional area (and with its normal pointing toward the observer) that is located at the same distance from the observer (Shepard 2017). The geometric albedo A G of a Lambertian-reflecting planet (Equation (1) where j(α) are the observational data at 1 au geocentric distance. The incoming solar flux, πF 0 , is wavelength dependent and with the r d 2 2 factor it is scaled from the Earth's surface to 1 au geocentric distance. We use the ASTM 2000 Standard Extraterrestrial Spectrum Reference E-490-00 (Livengood et al. 2011) for πF 0 . The geometric albedos are averaged over the 24 hr period.
The average geometric albedos of the Earth computed using Equation (2), thus assuming a Lambertian-reflecting planet (without an atmosphere), and using the data for j(α), in the seven spectral bands have been listed in Table 3. These albedos range from 0.143 (E1, 950 nm) to 0.353 (P2, 350 nm). Figure 4 shows these derived geometric albedos as functions of the wavelength for the five observations. The five spectra have similar shapes with relatively high values at the shortest wavelengths, consistent with Rayleigh scattering, a smooth decrease to relatively low values around 600 nm, and increase to a local maximum at 850 nm, followed by a sharp decrease to the lowest values at 950 nm, consistent with absorption by water vapor at ∼1 μm (Crow et al. 2011). From the figure, it is clear that the polar observations (P1 and P2) have very similar fluxes that are significantly larger than those of the equatorial observations, which are also very similar to each other. The higher fluxes for the poles could be due to the polar ice caps that are in view of the spacecraft (see Figure 1; visible by direct inspection) or to a larger cloud coverage fraction. This will be further discussed in a later paper. The E4 and E5 values are consistently increasingly greater than E1 at larger wavelength. This is consistent with the smaller illumination fraction having more solid surface in view and less ocean surface. Note that the E1 was around equinox and the E5 observation was around solstice, so seasonal effects could also cause a difference. Nonetheless, the equatorial observations are quite similar too each other.
The geometric albedos listed in Table 3 (Robock 1980). They are closer to the observations but they are either too low or too high due to the observations wavelength dependence or due to the higher polar observations. For the equatorial observations, the values are for most wavelengths just too high. Mallama et al. (2017) reports geometric albedo values for the Johnson-Cousins system or Sloan system. In this article values range from 0.392 (R C , 0.641 μm) to 0.688 (U, 0.360 μm) for Johnson-Cousins and from 0.388 (r′, 0.616 μm) to 0.722 (u′, 0.355 μm) for Sloan. These values are calculated using the same EPOXI data presented here but a different phase function is used to extrapolate to zero phase angle. They use the (steeper) phase function for the "realistic clouds" case plotted in Figure 7 of Tinetti et al. (2006b) and therefore they have higher values than the ones reported here. The geometric albedo values presented here are not necessarily an accurate representation of Earth's geometric albedo. In fact, the Lambertian scattering model cannot fully represent Earth's phase function, since it cannot account for the difference in retrieved geometric albedo from the polar data compared to the equatorial data, even though the geometry at the subsolar point must be the same for both models and thus should represent the same albedo. The values retrieved here serve as a means to extrapolate the phase curve to a limited extent beyond the modest range of phase angle directly observed and to compare the degree of consistency between the differing phase angle observations. These observations also provide a cautionary case study for future retrievals from actual exoplanets.
It is important to note that the phase angle dependence of the flux that is reflected by a planet with a Lambertian-reflecting surface and no atmosphere is generally less steep than that of a planet with an atmosphere. In particular, for a given geometric albedo A G , and at α = 90°, the flux of a Lambertian planet with a gaseous atmosphere is only 75% of that of a Lambertian planet without atmosphere, and at α = 60°, it is 78% (Stam et al. 2006). This would mean that a geometric albedo derived assuming a Lambertian-reflecting planet would be an underestimation compared to the real value. A model comparison to show the effect of the Lambertian model assumption is shown in Appendix.
As a next step, we used the wavelength dependent geometric albedos that we derived from the E1 data (the first row of Table 3) and computed the phase curves for a Lambertianreflecting planet (see Equation (1)), to allow a comparison with the data taken at the different phase angles. The comparison is shown in Figure 3. The phase curves based on the data from the E1 observations fit the other two equatorial observations, E4 and E5, very well, confirming the similarity between the three equatorial observations that we noted before when inspecting the geometric albedos (thus at a phase angle of 0°, Figure 4) over the whole phase angle range. The phase curves also confirm that the polar data (P1 and P2) are very different from the equatorial data. In particular, the reflected fluxes are much higher. The absolute differences between the hence computed phase curves and the observations data are listed in Table 4. Table 4 contains the absolute mean differences in flux (in the table "mean deviation") between all data points that were taken during the 24 hour period over which the data for each observation was acquired (25 for the 350, 750, and 950 nm   Table 2. The variability in the data is mainly due to the rotation of the Earth and to variations in the weather patterns during the 24 hour period over which the data is acquired. The solid curves represent the reflected fluxes as computed assuming a Lambertian-reflecting planet (without atmosphere) with (for each filter band) a geometric albedo A G derived from the E1 data. Note that while these curves fit not only the E1 data, they also fit the E4 and E5 data, while the P1 and P2 data are clearly higher than the curves. 7 Taken from https://nssdc.gsfc.nasa.gov/planetary/factsheet/earthfact. html, acquired on 01/12/2020. filters; and 97 for the 450, 550, 650, and 850 nm filters) and the computed value of the phase curve for the Lambertianreflecting planet at the phase angle pertaining to the observation (cf . Table 1). Because the Earth rotates within the 24 hour period and because local weather patterns, in particular clouds, change over time, the mean deviation also reflects these actual variations in time. The deviation is larger for the polar observations as these are brighter than the equatorial observations and are thus not fitted by the phase curve computed for E1. Table 4 also contains the standard deviation of the mean deviation, which also reflects the actual variations in time, such as the Earth's rotation. The standard error of the mean is given to assess the confidence in the mean. The values are given for Equatorial observation 1 (E1) as well, even though those data were used to calculate the phase curves. It can be seen that the mean errors for E4 and E5 are close to those for E1. Thus we can be confident that the differences between the data points of E4 and E5 and the phase curves are also due to time variations, such as the Earth's rotation. As the values for the standard deviation are mainly driven by time variations such as due to Earth's rotation, they are very similar for the equatorial observations as they cover the same region of the Earth (apart from the clouds).
For the polar observations, the mean errors are larger than the standard deviations indicating that there there is indeed an offset in the data that is not due to the Earth's rotation, i.e., the polar observations are significantly brighter. The standard deviations for the polar data are similar to those for the equatorial data, indicating that the Earth's rotation and cloud variations have a similar effect, except for the longest wavelengths in the north polar observations (P1). Those standard deviations are smaller than those for P2, which can also be seen in Figure 3, especially for the 950 nm data points of P1. This would indicate that at the time of observations ('09 March 27-28 for P1, see Table 1), the north polar region appeared to be more rotationally symmetric than the south polar region at the time of observation ('09 October 4-5 for P2, see Table 1). This is consistent with direct inspection of North American and Asian land masses in the images.
Finally, Table 4 also contains the standard errors of the mean, in order to proof that the limited number of data points (25 images for the 350, 750, and 950 nm filters; and 97 images for the other filters) do not undermine the above conclusions. They are all significantly smaller than the mean error, showing that the latter is relatively certain. 8 Therefore, the mean errors are expected to remain smaller (for E1, E4, and E5) and larger (for P1 and P2) than the standard deviations even when more Note. The last row shows the mean P/E, which is the ratio of the mean polar albedo (P1, P2) and the mean equatorial albedo (E1, E4, and E5). The uncertainty of this ratio is found by using propagation of error.  (350,450,550,650,750,850, and 950 nm) as derived from the measurements E1, E4, E5, P1, and P2 (see Table 2) assuming a Lambertianreflecting planet (without atmosphere). data would be added (note that the standard errors of the mean are somewhat larger for the filters with less observations, i.e., 350, 750, and 950 nm, than for the other filters. The observations can also be compared with different phase curves. Tinetti et al. (2006b) presents a phase light curve for the 0.5-0.9 μm wavelength range in Figure 7(a) (the general Earth model is presented in Tinetti et al. 2006a) that can be used. The albedos of the equatorial observations (of the relevant wavelength filters) match the "realistic clouds" model closely. They are slightly higher than the model but this is hard to assess with the low sampling rate of the model lines. The polar observations clearly do not agree with the "realistic clouds" model. They are closer to the "strato-cumulus" and "altostratus" model lines (low and middle clouds; Tinetti et al. 2006b), although slightly lower than those lines. They are lower because the model lines are for a fully cloud-covered disk-averaged spectra. From visual inspection of Figure 1, it can be seen that the EPOXI polar observations have more clouds in view than the equatorial observations. The midlatitudes are more cloudy than the tropics (Cowan et al. 2011), which explains why the polar observations have higher albedo values than the Tinetti et al. 2006b "realistic clouds" model due to the increased solid angle presented in the field of view. Modeling has been done, which confirms that the polar observations have slightly higher cloud coverage, this work will be published later.

Conclusion
The EPOXI mission imaged the Earth's full disk from afar and at five occasions: three equatorial and two polar views (the latter included one north and one south pole observation) at visible wavelengths (from 0.3 to 1.0 μm). The images were taken with the High-Resolution Instrument Visible onboard NASA's Deep Impact spacecraft and vary in distance to the Earth and in phase angle at the time of observation. The complete and calibrated version of the visible data is presented in this article along with a phase angle analysis.
Additional near-infrared data that was also acquired is not used in this article but has been presented before by Livengood et al. (2011) and that description serves as our main reference. A detailed description of the targeting, acquisition, selection, and processing (including calibration) along with a presentation of the old data (newly calibrated data is used in this article) can be found in that paper.
The Earth images have been taken through seven spectral filters centered 100 nm apart, covering wavelengths from 350 to 950 nm, and they were sampled every 15 minutes for four filters (450, 550, 650, and 850 nm) and every 60 minutes for the other three filters (350, 750, and 950 nm) each over a period of 24 hr. Using aperture photometry, disk-integrated Earth fluxes are retrieved for every filter, averaged over 24 hr and scaled to a 1 au range from the spacecraft. These fluxes are compared with a Lambertian-reflecting planet model (see Equation (1)), which is often used as a first order assumption or reference as comparison (Fujii et al. 2010(Fujii et al. , 2011Cowan et al. 2009Cowan et al. , 2011Robinson et al. 2011).
The fluxes of one of the equatorial observations are used to determine the geometric albedo of the Earth at every filter, assuming a Lambertian-reflecting planet (thus without scattering atmosphere), and to analyze the phase angle behavior for the other four observations. This approach allows us to match the other two equatorial observations over the phase angles, confirming the similarity between the three observations (recall that each set of observations was integrated over 24 hr, thus each covering a whole rotation of the Earth), but it underestimates the geometric albedos derived from the two polar observations. That could be due to the ice that is visible in those two observations.
Lambertian-reflecting planet (thus without scattering atmosphere; see Table 3). The hence derived geometric albedos range from 0.143 to 0.353 and are, as expected, wavelength dependent. The albedos derived from the polar observations are distinctly higher than those from the equatorial observations, with the difference being larger for longer wavelengths, where the influence of the atmosphere is smaller. The geometric albedo values that we retrieve are smaller than the classic (visible v-band, 0.549 μm) geometric albedo value of 0.434. 9 In this article values for geometric albedo range from 0.392 (R C , 0.641 μm) to 0.688 (U, 0.360 μm) for Johnson-Cousins and from 0.388 (r′, 0.616 μm) to 0.722 (u′, 0.355 μm) for Sloan. These values are calculated using the same EPOXI data presented here but a different phase function is used to extrapolate to zero phase angle. They use the (steeper) phase function for the "realistic clouds" case plotted in Figure 7 of Tinetti et al. (2006b).  The zero phase angle model is the correct geometric albedo and the difference with the other model lines shows the effect of using a Lambertian phase function to retrieve the geometric albedo. B1 has ocean surface albedo without cloud and B2 with clouds. A vegatation surface without cloud is shown in B3 and with clouds in B4. B5 is an ice surface without clouds. 9 Taken from https://nssdc.gsfc.nasa.gov/planetary/factsheet/earthfact. html, acquired on 01/12/2020.