Introduction
Haze is an atmospheric phenomenon which is associated with horizontal visibilities of less than l0 km and atmospheric relative
humidity (RH) less than 90 % (Liu et al., 2013). It is well known that haze occurs as a result of pollution. For example, the release of sulfur dioxide
(SO2), nitrogen oxides (NOx), and volatile organic compounds/hydrocarbons coupled with sunlight leads to aerosol formation, particles or the photochemical production of atmospheric particles (Sezer
et al., 2005; Pudasainee et al., 2006). These particles are called aerosol. Aerosol has a variety of effects on climate and
environment both directly and indirectly. The direct effect is through scattering, which cools the atmosphere and surface system or by
absorption of incoming solar radiation, which also cools the surface but warms the atmosphere. Indirectly, aerosol impacts on cloud
formation and the microphysical properties of clouds, which in turn influence cloud albedo and precipitation (Li et al., 2011) adding
to their negative health impacts. Aerosols are also the carriers of toxic substances such as heavy metals and polycyclic aromatic
hydrocarbons (Wilkomirski et al., 2011). In Beijing, under high pollution conditions, the concentrations of sulfate and nitrate have
been shown to account for one-third of the particle matter (PM10) mass and two-thirds of the PM2.5 mass, a part of which is
attributed to the additional secondary conversion of SO4-2 from SO2 and NO3-
from NOx (Ji et al., 2012). Haze has a significant effect on regional climatic phenomena, such as monsoon (Chung
et al., 2002; Evan et al., 2011), and on the environment, e.g. air quality (Lin et al., 2012) and visibility (Zhao et al., 2011). Aerosol
can adversely affects human health (Evan et al., 2011), especially for the elderly, children (American Academy of Pediatrics Committee
on Environmental Health, 1993), and even newborn children (Dadvand et al., 2013).
A thick smoke haze enveloped the eastern and northern part of China in December 2016. Pictures taken by cameras onboard the satellite
TERRA/AQUA show that the area of China affected by haze exceeded about 1.5 million square kilometres. The poor visibility resulted
in several highways and regional airports being closed for extended periods. The situation deteriorated significantly during the haze
event and became a matter of public concern.
Satellite observations of the reflectance of solar radiation at the top of the atmosphere are used to determine aerosol optical
thickness (AOT), which is used as an indicator of air quality (Kaufman et al., 2002). There are numerous attempts for the retrieval of
aerosol properties from satellite observations. AOT retrieval algorithms have been developed for use with the measurements of Moderate
Resolution Imaging Spectroradiometer (MODIS) (e.g. Dark-Target (Levy et al., 2013), DeepBlue (Hsu et al., 2013), the Multiangle
implementation of atmospheric correction (MAIAC) (Lyapustin et al., 2011)), Advanced Along-Track Scanning Radiometer (AATSR)
(e.g. AATSR Dual-Viewing (ADV) (Kolmonen et al., 2016; Sogacheva et al., 2017), Oxford-RAL Aerosol and Cloud (ORAC) (Thomas et al.,
2009), and Swansea University (SU) (North et al., 1999) algorithms). AOT is also derived from observations of the Multi-angle Imaging
SpectroRadiometer (MISR) (Diner et al., 2005), PARASOL's Polarization and Directionality of the Earth's Reflectances (POLDER) (Dubovik
et al., 2014), Sea-Viewing Wide Field-of-View Sensor (SeaWiFS) (Sayer et al., 2012) etc.
One challenge for the derivation of AOT long-term datasets from satellite observation is to generate comparable AOT data products from
the different instruments, which have limited lifetimes. Consequently, mature aerosol algorithms, which can be applied to data from
instruments on different platforms, are required. For example, the three MODIS aerosol algorithms have been applied to the Visible
Infrared Imaging Radiometer Suite (VIIRS) instrument and the three AATSR algorithms have been proposed to be applied to the
observations of the Sea and Land Surface Temperature Radiometer (SLSTR) instrument (Popp et al., 2016).
The MERIS instrument onboard Environmental Satellite (Envisat) provided valuable information for different applications (Verstraete
et al., 1999). There are several previous attempts to
develop AOT retrieval algorithms for MERIS, e.g. the Bremen AErosol Retrieval (BAER; von Hoyningen-Huene et al., 2003, 2011), and the
European Space Agency (ESA) standard aerosol retrieval (Santer et al., 2007). These had mixed success (Mei et al., 2017a). BAER has
limited accuracy away from dark-vegetated surfaces and primarily for non-absorbing aerosols (de Leeuw et al., 2015; Holzer-Popp et al.,
2013), while the ESA standard AOT retrieval tends to overestimate AOT (de Leeuw et al., 2015). The recently developed eXtensible Bremen
AErosol (XBAER) algorithm (Mei et al., 2017a, b) has been internally validated in the Aerosol-Climate Change Initiative (Aerosol-CCI)
project (Popp et al., 2016), and shows very promising results.
The newly launched (on 16 February, 2016) instrument Ocean Land Colour Instrument (OLCI) continues the work of MERIS as it contains all
MERIS channels. Theoretically it is possible to transfer the mature MERIS retrieval algorithms to the OLCI instrument. In this paper,
the XBAER algorithm has been applied to OLCI instrument for the first time. To our best knowledge, this is the first publication of AOT
retrieved from OLCI. Although Sentinel-3 has only recently been launched, applying XBAER to OLCI data we have identified a haze event
over Beijing, China, during December 2016. We use observations by OLCI during this episode to test our retrieval of AOT. This study is a necessary first step to observing the aerosol in the Arctic, which is an overarching long-term objective.
In this paper, the characteristics of OLCI and MERIS instruments are presented and compared in Sect. 2. The XBAER algorithm is
briefly explained in Sect. 3. Section 4 shows the comparison between OLCI and MERIS instruments – first XBAER OLCI-derived AOT results
and a comparison with AOT from MODIS/MISR and AERONET observations is shown and discussed from a global point of view. The AOT
retrieved during the regional haze event is also presented and discussed in Sect. 4. Conclusions are given in Sect. 5.
XBAER algorithm
The XBAER algorithm was designed for the retrieval of AOT from MERIS and similar observations. It has its own cloud-screening approach,
aerosol type selection and surface parameterization (Mei et al., 2017a, b). The cloud-screening algorithm minimizes cloud
contamination for aerosol retrieval in XBAER. The XBAER cloud-masking algorithm determines the presence of cloud by using (i) the
brightness of the scene, (ii) the homogeneity or variability of the top of the atmosphere reflectance, and (iii) cloud height
information (Mei et al., 2017b). The threshold values in the XBAER cloud-masking algorithm are selected by a two-step process. The
ranges for the thresholds were determined by using accurate radiative transfer modelling with different surface and atmospheric
scenarios. A histogram analysis has been used for different cloud, aerosol, and surface scenarios to estimate the optimal threshold
values for each criterion.
The XBAER algorithm uses a generic one-parametric surface parameterization for both land and ocean. XBAER uses a set of space–time-dependent spectral coefficients to describe surface properties. The spatial and temporal resolutions are 10 km and monthly,
respectively. The surface spectral reflectance can be determined simultaneously with AOT in an iterative procedure (Mei et al.,
2017a). This approach assumes that the wavelength-dependent properties of surface spectral reflectance are constrained by space- and
time-dependent spectral coefficients. The wavelength-independent single parameters (soil-adjusted vegetation index (SAVI) for land
retrieval and normalized differential pigment index (NDPI) for ocean retrieval) have been used as the “tuning” parameters. The
definitions of SAVI and NDPI are
SAVI=R(λ14)-R(λ7)R(λ14)+R(λ7)+L(1+L),L=1-2R(λ14)+1-(2R(λ14)+1)2-8(R(λ14)-R(λ7))2,
where R is the SSR and the subscript for the wavelength denotes the MERIS channel numbers defined in Table 1, and
NDPI=R(λ2)-R(λ5)R(λ3).
In this manner, XBAER is not limited to dark surfaces (ocean, vegetation) and also retrieves AOT over bright surfaces (e.g. desert,
semiarid, and urban areas).
XBAER uses MODIS Dark-Target aerosol type assumptions and the expected aerosol type for a given region and season is taken from an
analysis of Aerosol Robotic Network (AERONET) and Maritime Aerosol Network (MAN) observations for both land and ocean. AOT and surface
reflectance are retrieved by minimizing the difference between simulated and measured top-of-the-atmosphere (TOA) reflectance using a look-up table (LUT),
created by the radiative transfer software package SCIATRAN (Rozanov et al., 2014). Details of the XBAER algorithm can be found in Mei
et al. (2017a, b). A post-processing technique used in Aerosol-CCI project and the MODIS monthly snow fraction dataset have been
additionally applied to avoid unresolved clouds/snow (Popp et al., 2016).
Spectral response function of MERIS (dashed lines) and OLCI (solid lines) for overlap channels.
![]()
(a) Surface reflectances of the three selected surface types; (b) comparisons of the simulated TOA reflectance for different
combinations of MERIS and OLCI SRF values. Green, orange, and blue colours in (a) and (b) represent vegetation, soil, and water simulations. Filled circles in (b) are differences for simulations using MERIS and OLCI SRFs. Circles in (b) are differences for simulations with and without convolution.
![]()
Global comparison of OLCI XBAER AOT with AERONET observations for 2016 December. R and “match_ups” refer to the Pearson
correlation coefficient and the number of locations used in the validation respectively. The dashed lines are ±15%τ±0.10.
![]()
Comparison of the retrieved global monthly mean AOT at 0.55 µm for December 2016. (a) MODIS fire
product, (b) MISR, (c) MODIS (Dark-Target and DeepBlue combined), and (d) OLCI (XBAER).
![]()
Time series of meteorological parameters and pollutants during December 2016. (a) Wind direction and wind speed (kmh-1), (b) temperature (∘C) and relative humidity (%), (c) atmospheric pressure (hPa) and visibility (m), (d) SO2 and NO2 concentration (µgm-3), (e) AOT and Ångström coefficient (440–870 nm) (Alpha), and (f) PM2.5 hourly and daily concentration (µgm-3). The atmospheric components and meteorological data are from https://www.aqistudy.cn/historydata/index.php and https://www.wunderground.com/.
![]()
Daily MODIS RGB and AOT for East China (100–125∘ E, 25–45∘ N) during the period 16–21 December 2016 (from top left to bottom right).
![]()
Results
Verification
One important characteristic investigated is the instrument spectral response function (SRF) because it is the major difference between
MERIS and OLCI for overlap channels. Figure 1 shows the SRF for the MERIS and OLCI overlap channels. The OLCI SRF mean dataset
(https://sentinel.esa.int/web/sentinel/technical-guides/sentinel-3-olci/olci-instrument/spectral-response-function-data) has been
used. Differences between MERIS and OLCI SRF are identified but have negligible impact on the retrieved AOT.
In order to quantitatively investigate the impact of different SRFs, the TOA reflectances have been simulated with and without taking
SRF into account. The simulations have been determined by undertaking radiative transfer simulations using SCIATRAN for atmospheric
and surface conditions (Rozanov et al., 2014). The MERIS observation geometry for 2 July 2009 over Paris was used to perform
a forward simulation. In particular, the solar zenith angle, viewing angle, and relative azimuth were set to (32.32∘,
28.7∘, 30.65∘) as suggested in Mei et al. (2017a).
In order to design representative simulated scenarios, we define a comprehensive set of aerosol optical parameters, surface spectral
reflectances, and other atmospheric properties comprising temperature and pressure profiles, the profiles of the concentration of
gaseous absorbers and scattering. Suitable ranges of values for all relevant inputs for the radiative transfer model are obtained by statistical analysis of
corresponding global products (Mei et al., 2016a). For this purpose, we use the following parameters.
Surface reflectance: three typical surface types representing vegetation, soil and water, i.e. relatively dark land (vegetation-covered
city), bright land (desert), and water surface (ocean surface), were used. The typical vegetation and soil spectra are adapted from von
Hoyningen-Huene et al. (2011), the liquid water spectrum comes from the SCIATRAN database (see references in Rozanov et al.,
2014). Figure 2 shows the corresponding surface reflectance spectra for selected surface types.
Aerosol scenarios: within the ESA Aerosol-CCI project, a representative value for global mean AOT of 0.25 has been selected
(Holzer-Popp et al., 2013; de Leeuw et al., 2015). Thus an AOT of 0.25 was selected for the simulation of “vegetation” and “water”
cases. An AOT value of 0.5 was used for the “soil” scenario to represent a “real” case for the Sahara region. Moderately absorbing
(fine-mode radius rv, f=0.150 µm, coarse-mode radius rv, c=3.19 µm, fine-mode
variance σf=0.408, coarse-mode variance σc=0.754, fine/coarse-mode volumes
(µm3µm-3) are 0.055 and 0.038), pure maritime type (rv, f=0.150 µm, rv, c=3.19 µm, σf=0.408, σc=0.754, fine/coarse-mode volumes
(µm3µm3-) are 0.04 and 0.296) and dust aerosol model (rv, f=0.140 µm, rv, c=1.74 µm, σf=0.454, σc=0.687, fine/coarse-mode volumes
(µm3µm-3) are 0.02 and 0.157) were used for aerosol types.
Other atmospheric parameters: the profiles of temperature, pressure, and concentration of the gases ozone, O3, nitrogen
dioxide, NO2, and molecular oxygen, O2, and water vapour, H2O, which all absorb in the 400–900 nm
spectral region were provided by the Bremen 2-D chemical transport model (Sinnhuber et al., 2009).
In Table 1 the spectral channels of OLCI and MERIS are given. Figure 2a presents the surface spectral reflectance for the three surface
types selected. Figure 2b presents the simulated TOA differences for the above scenarios. The differences for all surface/atmospheric
conditions are less than 1.5 %. These are similar to the simulation with and without convolution for MERIS with the exception of
the O2A and water vapour channels. However, the potential impacts of different SRFs may also introduce some uncertainties to the XBAER
cloud mask due to the relatively strong impact of SRF to the O2A channels (about 20 % difference).
First XBAER AOT retrieval for OLCI and its validation
AERONET observations are considered to be the “ground truth” for satellite validation (Holben et al., 1998). Here, we collocate the
XBAER OLCI aerosol retrievals with the AERONET Version 3.0
(https://aeronet.gsfc.nasa.gov/new_web/Documents/AERONET-V3_News_Final.pdf, last access: 15 May 2017), Level 1.5 (level 2.0 for
both AERONET Version 2.0 and 3.0 are not available until 15 May 2017) (Holben et al., 1998; Smirnov et al., 2000). As AERONET does not
provide AOT at 0.55 µm, data are interpolated to 0.55 µm using quadratic fits on a log–log scale (Eck et al.,
1999). Since AERONET provides a point measurement with high temporal resolution while satellite observations represent a “regional”
measurement depending on the satellite spatial resolution for a particular overpass time, spatial statistics for the OLCI data are
calculated and compared to the temporal statistics of the AERONET observations taken within ±30 min of OLCI overpass
following the spatial–temporal technique of Ichoku et al. (2002).
Figure 3 is a plot which compares XBAER-derived and AERONET-observed AOT at 0.55 µm. The collocations of Fig. 3 contain
various surface and aerosol types, which ensure a wide representativeness of the validation. A total of 733 collocations were found for
December 2016. The colour of each ordered pair (0.025×0.025 increment) represents the number of such matchups. Negative AOT
(>-0.1) values are possible and reasonable as a result of the noise in satellite observations (Levy et al., 2007) and uncertainties of surface parameterization. The comparison here excluded
negative values and only AOT values between 0.0 and 2.5 are used following the validation method of other aerosol products (Sayer
et al., 2012; Levy et al., 2013). The validation contains various surface and aerosol types, which ensures a wide representativeness of
the validation. The regression equation is y=(0.81x±0.04)+(0.11±0.01) with slightly higher correlation compared to the
first MERIS validation (R=0.82 vs. R=0.78) (Mei et al., 2016a). The AOT is reasonably correlated between the two datasets
(R=0.82), with increased scatter for high aerosol loadings. The majority of the data (87.5 %) are for low aerosol loadings
(AOT <0.3). The comparison between XBAER AOTs and AERONET observations shows the acceptable quality of the first OLCI XBAER
results.
Figure 4 shows the global monthly AOT of December 2016 for MODIS collection 6 (Levy et al., 2013), MISR (Diner et al., 2005) and OLCI
(XBAER) algorithm. In order to identify biomass burning events, the active fire points of MODIS
(https://lance.modaps.eosdis.nasa.gov/cgi-bin/imagery/firemaps.cgi) are added to the figures. Please note that the MISR
“FIRSTLOOK” product is used because the monthly Land Surface and Aerosol products are not yet processed for December 2016 (NASA Langley ASDC, personal
communication, 2017). MODIS/MISR on board of TERRA and OLCI on board of Sentinel-3 have very
similar overpass time (within 30 min difference). Therefore, all four results should show similar patterns for large AOT from desert
dust events over the Sahara, biomass burning over West Africa and the Amazon region, and anthropogenic pollution over India and East Asia. In
Fig. 4, XBAER AOT from OLCI shows similar patterns as the AOT from MODIS and MISR for both land and ocean. However, there are
differences in the magnitude of the AOTs. Biomass burning over Africa, as observed in the MODIS active fire product, produces a
“plume belt” of high AOT near the Equator. This is observed in all three AOT products. The AOT distribution pattern over India, which
depends on the unique meteorological conditions and emissions, is captured by the three AOT data products as well. MODIS and OLCI show
similar pattern and magnitude of large AOT over eastern China while the values from MISR are slightly lower, which may be due to the
relatively small sampling compared to MODIS and OLCI. However, the retrievals from XBAER over Australia are higher than those of MODIS
and MISR. In addition to potential contamination by thin clouds observed in the RGB composite figures, the calibration uncertainties
associated with a new instrument may also contribute to the bias of XBAER-derived AOT. The large AOT differences over the Sahara may stem, in
part, from different assumptions in the different algorithms for bright surfaces (Lyapustin et al., 2011b; Mei et al.,
2016a). Different patterns over the Amazon can most likely be attributed to the use of different cloud-screening methods. The global
patterns obtained indicate that the generic XBAER algorithm works over both dark and bright surfaces using its flexible surface
parameterization approach. For relative dark surfaces, the one-parametric surface parameterization is dominated by the first term (SAVI
or NDPI tuned term) making XBAER behave like the Dark-Target-like retrieval algorithm. For bright surfaces such as desert, XBAER
becomes similar to the DeepBlue AOT retrieval algorithm.
Same as Fig. 6 but for MERRA AOT.
![]()
Same as Fig. 6 but for OLCI.
![]()
Same as Fig. 6 but for NO2 from the GOME2a/GOME2b combined IUP–UB product.
![]()
Beijing haze event observed by OLCI
In the following we show the ability of the retrievals of XBAER used with OLCI data to resolve spatial aerosol patterns on a synoptic
scale. A prolonged haze event was observed over Beijing during the period 16–21 December 2016. The intention of applying XBAER to
this event is to show the potential of the retrieval to resolve aerosol patterns at a local level and thus being able to support future
studies analysing such events. This event is investigated by both ground-based measurements and satellite observations. Figure 5a shows
that winds at the surface were weak, with a daily averaged wind speed lower than 3.5 ms-1 during the period, causing the
accumulation of pollutants on a regional scale. The low temperature and high surface pressure near the surface indicate relatively
stable atmospheric conditions in the vertical direction (Fig. 5b and c). The dispersion of pollutants out of the boundary layer is
therefore slow. The relative humidity remained high (Fig. 5b), causing the aerosol particles size to increase by the uptake of water
(Winkler, 1988), thus making the haze event stronger. Under these meteorological conditions, pollutants can accumulate over the North
China Plain (NCP) (Li et al., 2011).
Figure 5d shows the time series of concentration of SO2 and NO2 in the boundary layer provided by ground-based
measurements. The concentrations of SO2 and NO2 for haze periods are three to five times larger than those on
relatively clear days. We thus assume anthropogenic activities to be the major source of AOT. Figure 5e and f show the AOT from AERONET
sites and the time series of the daily mean concentration of PM2.5 in Beijing with clearly increased values in the same
timespan (16–21 December). The lack of larger AOT values observed by AERONET is likely going back to too strict cloud-screening
procedures. The daily mean concentrations of PM2.5 during 16–21 December 2016 ranged from 107.1 to
394.5 µgm-3, which is far above the daily PM2.5 limit of the new threshold value set as the Chinese Ambient
Air Quality standard (75 µgm-3)
(http://transportpolicy.net/index.php?title=China:_Air_Quality_Standards). A large Ångström coefficient (Fig. 5e) shows that fine
particles dominate during this period. In summary, we find that the cause of the haze event in Beijing and northeastern China can be traced back
to (1) the stable meteorological conditions (low wind speeds and temperature inversion), (2) local emissions, and (3) high relative
humidity.
Figure 6 shows the MODIS/Terra-derived AOT for the haze period. According to Fig. 6, this intense part of the haze episode has been
partly observed by MODIS. However, a large part of it (under cloud-free conditions) during the first 3 days is missing, mainly due
to cloud masking applied in the MODIS aerosol retrieval. Figure 7 shows the AOT from Modern-Era Retrospective analysis for Research and
Applications Version 2 (MERRA-2) simulation (Rienecker et al., 2011) in order to exclude the impact from cloud screening. According to
Fig. 7, the shape of the area covered by high AOT in MERRA remains stable except for 20 December, indicating the relative stable
meteorological condition during the haze period. Due to the narrower swath width of OLCI compared to that of MODIS (1270 km
vs. 2300 km), OLCI has a longer “revisit” time for a repetitive observation of the ground scene in Beijing. According to
Figs. 6 and 8, XBAER discards fewer clear-sky ground scenes than the MODIS retrieval, in particular on 18 and 19 December over
eastern China. For this period, the aerosol over eastern China has been characterized as “moderately absorbing aerosol”. The single-scattering albedo at
0.675 µm from AERONET has values between 0.88 and 0.91, indicating relatively strong absorption from anthropogenic
activities. The magnitude of AOT for the overlap regions between OLCI and MERRA are comparable according to Figs. 7 and 8. The regions
of high aerosol agree well with the areas having high NO2 columns in Global Ozone Monitoring Experiment 2 (GOME2) (Richter
et al., 2011) for the corresponding time period as presented in Fig. 9. Figure 8 illustrates that cloud masking, surface treatment, and
aerosol type selection in XBAER all work well for the detection of extreme haze events. Studies like the one by Zheng et al. (2015)
usually focus on the origin of such plumes and the speciation of aerosol particles on a city level. XBAER results utilizing
multi-spectral imagery such as provided by OLCI can support this kind of study to identify plume transport and extension. This
implies that the OLCI instrument can provide important data on AOT for atmospheric research.
Discussions
In this study, we have applied XBAER to data from the OLCI instrument onboard Sentinel-3 for the first time on both synoptic and global
scale. The potential differences caused by different spectral response functions for OLCI and MERIS have been investigated by using
SCIATRAN to generate representative simulated scenarios for dust aerosol type over desert, moderately absorbing aerosol over vegetation
regions, and maritime aerosol over water. The overall differences for all selected channels for XBAER are smaller than 1.5 %. This
implies that XBAER can be used to retrieve AOT from OLCI. Although relatively large differences caused by SRFs (approximately 20 %)
have been found for the O2A channels, the global retrieval of OLCI shows that the original MERIS cloud masking, which includes
the use of O2A channels, works well for OLCI and can potentially even be improved as only MERIS-heritage channels have been
used so far with OLCI.
The global monthly mean XBAER AOT maps for December 2016 show good agreement with those by MODIS and MISR. The comparison with AERONET
measurements reveals that XBAER can provide promising results over both dark and bright surface. The first comparison with AERONET
shows acceptable agreement between the two data sets, with a regression yielding y=(0.81x±0.04)+(0.11±0.01) and correlation
of R=0.82. The global retrievals confirm that XBAER is valid for both dark and bright surfaces because of its use of an optimized
monthly global SSR spectral coefficients dataset.
A significant haze event during December 2016 over Beijing has been analysed in this paper based on ground-based and satellite
observations to show the potential of the retrieval to resolve aerosol patterns at a local level and thus being able to support future
studies analysing such events. This large haze event has been attributed to the large local emissions under unfavourable meteorological
conditions (temperature inversion in vertical direction and no advection). The MODIS/Terra- and OLCI-derived AOT both detect the haze
event. However, due to cloud screening, the MODIS AOT partly misses it while the OLCI AOT is able to detect the main pattern of haze for
clear conditions. The overlap retrieval for both MODIS and OLCI has similar values, indicating that OLCI provides another useful data
source for air pollution monitoring.
Although the study shows that XBAER can be applied to OLCI observations for synoptic to global applications, several important issues
need to be addressed in the future work. Potential cloud contamination due to both the relative large calibration uncertainty of OLCI
compared to MERIS as well as the impact of SRF on O2A channel need to be investigated with the new version of level 1 TOA
reflectance dataset. Modification or improvement for OLCI cloud screening will be included, besides the criteria of brightness,
texturing/variability, and cloud altitude of the scenes (Mei et al., 2017b). The underestimation of AOT over regions like the Sahara could be explained by the spheroid dust model adapted from MODIS-DT
algorithm due to the impact of non-sphericity of dust particles on the aerosol phase function (Mei et al., 2016a) – a new spheroid model
accounting for aerosol particle non-sphericity will be included in the new version (Dubovik et al., 2006). The cloud-screening
evaluation shows that approximately 5–10 % clouds may be misclassified as retrievable clear cases for MERIS (Mei et al., 2016b),
which introduces both bias and potential patchiness of XBAER-derived AOT for OLCI. Thus a new cloud post-processing, following the
AATSR dual-view (ADV) algorithm (Sogacheva et al., 2017), will be applied to discard the pixels that might potentially be affected by
cloud (cloud edge, very thin cloud, and so on).