The angular distribution of the light reflected by the Earth's surface influences top-of-atmosphere (TOA) reflectance values. This surface reflectance anisotropy has implications for UV/Vis satellite retrievals of albedo, clouds, and trace gases such as nitrogen dioxide (NO 2 ). These retrievals routinely assume the surface to reflect light isotropically. Here we show that cloud fractions retrieved from GOME-2A and OMI with the FRESCO and OMCLDO2 algorithms have an east-west bias of 10 % to 50 %, which are highest over vegetation and forested areas, and that this bias originates from the assumption of isotropic surface reflection. To interpret the across-track bias with the DAK radiative transfer model, we implement the bidirectional reflectance distribution function (BRDF) from the Ross-Li semi-empirical model. Testing our implementation against state-of-the-art RTMs LIDORT and SCIATRAN, we find that simulated TOA reflectance generally agrees to within 1 %. We replace the assumption of isotropic surface reflection in the equations used to retrieve cloud fractions over forested scenes with scattering kernels and corresponding BRDF parameters from a daily high-resolution database derived from 16 years' worth of MODIS measurements. By doing this, the east-west bias in the simulated cloud fractions largely vanishes. We conclude that across-track biases in cloud fractions can be explained by cloud algorithms that do not adequately account for the effects of surface reflectance anisotropy. The implications for NO 2 air mass factor (AMF) calculations are substantial. Under moderately polluted NO 2 and backwardscattering conditions, clear-sky AMFs are up to 20 % higher and cloud radiance fractions up to 40 % lower if surface anisotropic reflection is accounted for. The combined effect of these changes is that NO 2 total AMFs increase by up to 30 % for backward-scattering geometries (and decrease by up to 35 % for forward-scattering geometries), which is stronger than the effect of either contribution alone. In an unpolluted troposphere, surface BRDF effects on cloud fraction counteract (and largely cancel) the effect on the clearsky AMF. Our results emphasise that surface reflectance anisotropy needs to be taken into account in a coherent manner for more realistic and accurate retrievals of clouds and NO 2 from UV/Vis satellite sensors. These improvements will be beneficial for current sensors, in particular for the recently launched TROPOMI instrument with a high spatial resolution.

We present the first comparison between new fAPAR and LAI products derived from the GlobAlbedo dataset and the widely-used MODIS fAPAR and LAI products. The GlobAlbedo-derived products are produced using a 1D two-stream radiative transfer (RT) scheme designed explicitly for global parameter retrieval from albedo, with consistency between RT model assumptions and observations, as well as with typical large-scale land surface model RT schemes. The approach does not require biome-specific structural assumptions (e.g., cover, clumping, understory), unlike more detailed 3D RT model approaches. GlobAlbedo-derived values of fAPAR and LAI are compared with MODIS values over 2002-2011 at multiple flux tower sites within selected biomes, over 1200ˆ1200 km regions and globally. GlobAlbedo-derived fAPAR and LAI values are temporally more stable than the MODIS values due to the smoothness of the underlying albedo, derived via optimal estimation (assimilation) using an a priori estimate of albedo derived from an albedo "climatology" (composited multi-year albedo observations). Parameters agree closely in timing but with GlobAlbedo values consistently lower than MODIS, particularly for LAI. Larger differences occur in winter (when values are lower) and in the Southern hemisphere. Globally, we find that: GlobAlbedo-derived fAPAR is~0.9-1.01ˆMODIS fAPAR with an intercept of~0.03; GlobAlbedo-derived LAI is~0.6ˆMODIS LAI with an intercept of~0.2. Differences arise due to the RT model assumptions underlying the products, meaning care is required in interpreting either set of values, particularly when comparing to fine-scale ground-based estimates. We present global transformations between GlobAlbedo-derived and MODIS products.

Surface albedo is of crucial interest in land–climate interaction studies, since it is a key parameter that affects the Earth’s radiation budget. The temporal and spatial variation of surface albedo can be retrieved from conventional satellite observations after a series of processes, including atmospheric correction to surface spectral bi-directional reflectance factor (BRF), bi-directional reflectance distribution function (BRDF) modelling using these BRFs, and, where required, narrow-to-broadband albedo conversions. This processing chain introduces errors that can be accumulated and then affect the accuracy of the retrieved albedo products. In this study, the albedo products derived from the multi-angle imaging spectroradiometer (MISR), moderate resolution imaging spectroradiometer (MODIS) and the Copernicus Global Land Service (CGLS), based on the VEGETATION and now the PROBA-V sensors, are compared with albedometer and upscaled in situ measurements from 19 tower sites from the FLUXNET network, surface radiation budget network (SURFRAD) and Baseline Surface Radiation Network (BSRN) networks. The MISR sensor onboard the Terra satellite has 9 cameras at different view angles, which allows a near-simultaneous retrieval of surface albedo. Using a 16-day retrieval algorithm, the MODIS generates the daily albedo products (MCD43A) at a 500-m resolution. The CGLS albedo products are derived from the VEGETATION and PROBA-V, and updated every 10 days using a weighted 30-day window. We describe a newly developed method to derive the two types of albedo, which are directional hemispherical reflectance (DHR) and bi-hemispherical reflectance (BHR), directly from three tower-measured variables of shortwave radiation: downwelling, upwelling and diffuse shortwave radiation. In the validation process, the MISR, MODIS and CGLS-derived albedos (DHR and BHR) are first compared with tower measured albedos, using pixel-to-point analysis, between 2012 to 2016. The tower measured point albedos are then upscaled to coarse-resolution albedos, based on atmospherically corrected BRFs from high-resolution Earth observation (HR-EO) data, alongside MODIS BRDF climatology from a larger area. Then a pixel-to-pixel comparison is performed between DHR and BHR retrieved from coarse-resolution satellite observations and DHR and BHR upscaled from accurate tower measurements. The experimental results are presented on exploring the parameter space associated with land cover type, heterogeneous vs. homogeneous and instantaneous vs. time composite retrievals of surface albedo.

Abstract. The angular distribution of the light reflected by the Earth's surface influences top-of-atmosphere (TOA) reflectance values. This surface reflectance anisotropy has implications for UV/Vis satellite retrievals of albedo, clouds, and trace gases such as nitrogen dioxide (NO2). These retrievals routinely assume the surface to reflect light isotropically. Here we show that cloud fractions retrieved from GOME-2A and OMI with the FRESCO and OMCLDO2 algorithms have an East-West bias of 10 % to 50 % over rugged terrain, and that this bias originates from the assumption of isotropic surface reflection. To interpret the across-track bias with the DAK radiative transfer model, we implement the Bidirectional Reflectance Distribution Function (BRDF) from the Ross-Li semi-empirical model. Testing our implementation against state-of-art RTMs LIDORT and SCIATRAN, we find that simulated TOA reflectance generally agrees to within 1 %. By replacing the assumption of isotropic surface reflection in the cloud retrievals over vegetated scenes with scattering kernels and corresponding BRDF parameters from a daily, high-resolution database derived from 16 years' worth of MODIS measurements, the East-West bias in the retrieved cloud fractions largely vanishes. We conclude that across-track biases in cloud fractions can be explained by cloud algorithms not adequately accounting for the effects of surface reflectance anisotropy. The implications for NO2 air mass factor (AMF) calculations are substantial. Under moderately polluted NO2 and backscatter conditions, clear-sky AMFs are up to 20 % higher and cloud radiance fractions up to 40 % lower if surface anisotropic reflection is accounted for. The combined effect of these changes is that NO2 total AMFs increase by up to 30 % for backscatter geometries (and decrease by up to 35 % for forward scattering geometries), stronger than the effect of either contribution alone. In an unpolluted troposphere, surface BRDF effects on cloud fraction counteract (and largely cancel) the effect on the clear-sky AMF. Our results emphasize that surface reflectance anisotropy needs to be taken into account in a coherent manner for more realistic and accurate retrievals of clouds and NO2 from UV/Vis satellite sensors. These improvements will be beneficial for current sensors, in particular for the recently launched TROPOMI instrument with a high spatial resolution.

CAR (Cloud Absorption Radiometer) is a multi-angular and multi-spectral airborne radiometer instrument, whose radiometric and geometric characteristics are well calibrated and adjusted before and after each flight campaign. CAR was built by NASA (National Aeronautics and Space Administration) in 1984. On 16 May 2008, a CAR flight campaign took place over the well-known calibration and validation site of Railroad Valley in Nevada, USA (38.504°N, 115.692°W). The campaign coincided with the overpasses of several key EO (Earth Observation) satellites such as Landsat-7, Envisat and Terra. Thus, there are nearly simultaneous measurements from these satellites and the CAR airborne sensor over the same calibration site. The CAR spectral bands are close to those of most EO satellites. CAR has the ability to cover the whole range of azimuth view angles and a variety of zenith angles depending on altitude and, as a consequence, the biases seen between satellite and CAR measurements due to both unmatched spectral bands and unmatched angles can be significantly reduced. A comparison is presented here between CAR's land surface reflectance (BRF or Bidirectional Reflectance Factor) with those derived from Terra/MODIS (MOD09 and MAIAC), Terra/MISR, Envisat/MERIS and Landsat-7. In this study, we utilized CAR data from low altitude flights (approx. 180 m above the surface) in order to minimize the effects of the atmosphere on these measurements and then obtain a valuable ground-truth data set of surface reflectance. Furthermore, this study shows that differences between measurements caused by surface heterogeneity can be tolerated, thanks to the high homogeneity of the study site on the one hand, and on the other hand, to the spatial sampling and the large number of CAR samples. These results demonstrate that satellite BRF measurements over this site are in good agreement with CAR with variable biases across different spectral bands. This is most likely due to residual aerosol effects in the EO derived reflectances.