Validation of Two MODIS Aerosols Algorithms with SKYNET and Prospects for Future Climate Satellites Such as the GCOM-C / SGLI

1 Earth Observation Research Center/JAXA, 2-1-1 Sengen, Tsukuba, Ibaraki 305-8505, Japan 2 Center for Environmental Remote Sensing (CEReS), Chiba University, 1-33 Yayoi-Cho, Inage-Ku, Chiba 263-8522, Japan 3 Center for Environmental Measurement and Analysis, National Institute for Environmental Studies, 16-2 Onogawa, Tsukuba, Ibaraki 305-8506, Japan 4Department of Network and Computer Engineering, Tokai University, 2-28-4 Tomigaya, Shibuya-ku, Tokyo 151-0063, Japan


Introduction
The determination of the optical properties of aerosols and their size distribution around the globe has been a significant contemporary research effort of late [1].Some of the major factors that have enabled this progress are the better spectral and spatial capacities of satellites and ground-based radiometers, the improvement of the aerosol signal filtering methods, the better knowledge of the aerosols particles' shapes, and so forth.The qualitative and quantitative importance of the data collected through the continuous monitoring of aerosols and their daily global coverage, by various satellites, has permitted a better characterization of the role of the aerosols in the climate dynamics and the understanding of their temporal and local variation.The most common aerosols optical and physical properties used for this characterization are the AOT, the AE, the particle size distribution, the single scattering albedo (SSA), the aerosol phase function, the asymmetry factor (AF), the refractive index (RI), and so forth.To obtain these aerosols characteristics from satellites' observations, a detailed model of the aerosols properties is required [2].The increasing number and spatial expansion of field measurement campaigns and ground sites' coverage have helped to gather statistically meaningful data for the development of aerosols models.Ground collected data not only serve as preliminary input for aerosols models but also are reliable evaluation and validation tools for aerosols products derived from satellites (GMS-5/SVISSR, NOAA-AVHRR, Terra-and Aqua-MODIS and MISR, OrbView-2/SeaWiFS, ENVISAT-MERIS, etc.).The objective of these validation studies is to improve the quality of the aerosols properties retrieved from satellites' observations as well as the confidence level in these products.Since the launch of Terra-and Aqua-MODIS, ground data from AERONET have been used to validate and improve the MODIS AOT over ocean [3] and land [4,5].NOAA-AVHRR AOT data have been also evaluated using sun photometers of AERONET and others [6,7].The retrieval quality over the ocean of POLDER/ADEOS aerosols measurements has been also conducted through comparisons with AERONET sun photometer data [8].
The increasing accuracy of aerosols products, as a result of all these evaluation/validation and correction studies, could be used to reduce the uncertainties associated with the aerosol radiative forcing of the global climate [9].Due to the complex and diverse composition, size, structure, and origin of aerosols and though a lot of progress has been made in the satellites' observations, there are still various challenges for a globally acceptable accuracy of the aerosol optical and microphysical properties of different types of aerosols.The high spectral and spatial resolutions of climate-oriented satellites such as the GCOM-C/SGLI, scheduled to be launched around 2015 by the Japan Aerospace Exploration Agency (JAXA), aim at accomplishing better distinctions between aerosols particles and other atmospheric components.However, the programmed long lifespan of such satellites (3 series) and the connectivity between similar/different satellites' series will pose other accuracy issues proper to long-term climate observations [10][11][12] that may only be alleviated with better accuracy retrieval algorithms.
In the present study, the potential performance of future satellites' products is discussed through a validation study of aerosol retrievals by present satellites' (Terra-and Aqua-MODIS satellites) algorithms using spectrally compatible channels with those of the coming GCOM-C/SGLI.This validation is conducted against a three-year (2008-2010) daily ground-truth dataset from three SKYNET observation sites (Fukue, Hedo, and Miyako).These ground sites are located in the Pacific East Asian region, an area at the confluence of actively changing and mixed aerosols (natural and anthropogenic particles resulting from rapid industrialization).The retrievals examined are from two algorithms using similar observations, Terra-and Aqua-MODIS satellites' calibrated radiances (MOD/MYD021, version 5.1).The first algorithm is named here as MODIS-GLI, as it was previously applied to the Global Imager (GLI) onboard the Advanced Earth Observation Satellite-II (ADEOS-II).The second algorithm is the MODIS-NASA algorithm, used for MODIS standard products' retrievals by the National Aeronautics and Service Administration (NASA).The aerosols properties investigated are the aerosol optical thickness (AOT) and the Ångström exponent (AE).The differences between the satellites' retrievals and ground observations are analyzed, issues plaguing the accuracy of each of these satellites' algorithms are discussed, and lessons for future satellites' algorithms improvements, such as the algorithm of the future GCOM-SGLI, are drawn.To conduct this work, the following steps are adopted: (1) the introduction, highlighting the basic settings and motives of the study; (2) the description of the ground data collection approach and the satellites' algorithms main features; (3) the presentation of the main characteristics of the aerosols properties of the study area; (4) the validation scheme of the satellites' retrievals, and, discussion on the uncertainties plaguing these retrievals; (5) the lessons for future satellites' aerosols products; (6) the summary of the results obtained.

Ground Data Collection and Satellites' Algorithms
The analyses conducted in this study aim at validating tropospheric aerosols' retrievals from two algorithms, using similar observations from both Terra-and Aqua-MODIS satellites.These satellites have spectrally compatible channels with the future GCOM-C/SGLI satellite.This validation study is conducted through a comparison of the satellites' retrievals with ground-truth measurements from three SKYNET sites: Fukue (32.752 ∘ N latitude, 128.682 ∘ E longitude, and 50 m altitude), Hedo (26.867 ∘ N latitude, 128.249 ∘ E longitude, and 65 m altitude), and Miyako (24.737 ∘ N latitude, 125.327 ∘ E longitude, and 50 m altitude).These sites are islands located in the Pacific East Asian region, south of Japan.The choice of these three locations for the satellites' retrieved aerosols properties' evaluation study derives from the fact that they are at the confluence of one of the most complex mixtures of aerosols types.It is known that the Asian atmospheric environment has been worsened rapidly by Asian dust, anthropogenic gases, and aerosols in recent years [13]; it has been also the subject of various aerosols projects.Some of the examples are the Asian-Pacific Regional Aerosol Characterization Experiment (ACE-Asia) [14][15][16] which was a multiplatform project where aerosols, radiative fluxes measurements were conducted over the East China Sea; and the Atmospheric Brown Cloud East Asian Regional Experiment 2005 [17,18], with the study of aerosols radiative characteristics and aerosol direct radiative forcing.The SKYNET, whose data are used to validate against the satellite retrievals in this study, is a well-developed and maintained network of climate radiation stations with ground sites spread all over Japan and the East Asian region [17,19].It serves for the monitoring of aerosol and cloud properties as well as other weather/climate variables.The basic aerosols properties at the SKYNET sites are obtained from radiation measurements performed with sky radiometers.These radiometers record direct sky radiance intensities at 11 predefined wavelengths (315, 340, 380, 400, 500, 675, 870, 940, 1020, 1600, and 2200 nm) and scattering angles.The measurement frequency is 10 or 15 minutes.The solar directbeam transmission and sky radiance obtained from sky radiometer measurements are used for aerosols properties such as the AOT, the aerosol column size distributions, the complex refractive indices, the single scattering albedo, and the asymmetry factor, at multiple wavelengths.A program package named Skyrad.pack[20] is used for the retrieval of these parameters.The SKYNET helps to extend the global aerosol monitoring from the ground to locations not always covered by the better known AErosol RObotic NETwork (AERONET).Previous studies have shown a good match between the SKYNET and the AERONET aerosols measurements in areas where both networks coexist, confirming the quality and consistency of the SKYNET measurements.One of these is a study conducted over the Loess Plateau of Northwest China.The study compares the aerosol optical depths at two common wavelengths (675 and 870 nm) observed by the CIMEL sun photometer of the AERONET and those observed by the POM-01 sky radiometer of the SKYNET, based on 703 clear-sky measurements during the year 2009 [21].The study shows a high correlation coefficient (0.99) and a low relative difference of about 3% between the measurements of the two instruments.Bi et al. [22] also showed a good match between the AOT retrieved by both networks at four major channels (440, 675, 870, and 1020).
Both the SKYNET and AERONET aerosols networks have been very often used to validate and refine aerosols models for the satellites' retrievals [3,4,9,[23][24][25][26].The multiplication of ground stations and the improvement of the quality of the data obtained with the introduction of new parameterization schemes in the radiative transfer derivation codes and a more rigorous cloud mask application have increased the confidence level in the use of these groundbased data for satellite products' validation purpose [25,27].
In this study, daily aerosol retrievals from both Terraand Aqua-MODIS at 0.5 m wavelength are analyzed.The matching time and space area between these satellite data and those from the ground stations are ±30 min and 25 km around each ground site (time and space average are applied), respectively.The parameters studied are the AOT at 0.5 m and the AE from visible channels.The MODIS-NASA data (MOD/MYD04, version 5.1) produce the AOT at 0.55 m.Based on the wavelength relationship between the AE and the AOT, the latter are interpolated to 0.5 m, in order to match the retrievals from the MODIS-GLI algorithm and the ground data AOT.
The basic retrieval scheme of the MODIS-GLI algorithm over the ocean, described in Nakajima and Higurashi [23], Higurashi and Nakajima [28], and Higurashi et al. [29], uses the visible (0.660 m) and near-infra-red (0.870 m) channels to retrieve the AOT.The satellite-received radiances are synthesized in four look-up tables (LUTs) derived from radiative transfer calculations under the Mie theory scheme.The algorithm assumes a complex refractive index (RI = 1.5-0.005)and a bimodal size distribution expressed as where  is the number of lognormal aerosol modes,   is the median or geometric mean radius at each mode,   is the variance or width of the mode,   represents the particle peak volume concentration of the mode; this is an undetermined parameter to be retrieved from satellite-received radiances.
The AOT and the peak ratio  of the size distribution are retrieved by comparing measured and theoretical (precalculated LUTs) reflectance values.The optimal values of AOT and , that minimize the root mean square deviation between observed and theoretical reflectance, are iteratively searched.The  parameter, representing the contribution of the large particles in the aerosol total load, is used to calculate the AE.
For the MODIS-NASA algorithm, the retrieval scheme over the ocean, described by Kaufmann and Tanré [30], uses seven wavelengths from 0.47 to 2.13 m to retrieve the AOT and the volume distribution (in the range of 0.08-5 m radius) from measured radiances.LUTs from radiative transfer calculations of the MODIS radiance spectral bands are used to retrieve aerosols parameters at each particle size distribution mode: the small particle mode or accumulation mode particles (5 cases) and the large or coarse mode particles (6 cases).The AOT and the volume distribution are the primarily derived parameters from the algorithm, while the AE is (as in the MODIS-GLI algorithm) a secondarily derived parameter.The AOT is obtained from the measured radiances, by linear interpolation between five optical thickness values of each aerosol particle mode combination.The two modes that give the best selected value (ratio between the modes) are the selected aerosol models for the retrieval.

Aerosols Characteristics
Prior to the evaluation of the satellite retrievals against ground data, a quick description of some of the major characteristics of the aerosols distribution observed in the region covering the ground sites selected for this study may be useful to understand the aerosols of the region.Figure 1 shows the prevailing pattern of the AOT and AE distributions in the region, as retrieved by the MODIS-GLI algorithm during winter (month of January 2010), and the geographical location of the SKYNET ground sites selected.In this figure, the aerosol distribution shows high AOTs near the coastal area of the southwest (adjacent to China) and in the northwest area representing the sea of Japan.But higher AE (finer particles) are visible in the former than the latter.Closer ranges of both AOT and AE can be seen in a north-south band running along the SKYNET sites, with however thicker layers and finer particles in the north than the south.For the three-year data set matching the satellites' overpass time, the northern site (Fukue) has the thickest aerosols (mean AOT of 0.30 ± 0.23) as well as the finest particles (mean AE of 0.97 ± 0.35).The other 2 sites (Hedo and Miyako) located in the south have mean AOT values of 0.26 ± 0.20 and 0.27 ± 0.15, respectively and AE of 0.95 ± 0.35 and 0.89 ± 0.36, respectively.As we move from the north to the south, that is, further from the major coastal area and towards the central parts of the Pacific Ocean (Figure 1), the aerosol particle size increases (AE values get smaller).aerosols measurements from the ground (SKYNET sites) timely matching the satellite retrievals (from MODIS-GLI and MODIS-NASA algorithms) for the AOT and AE, respectively, during the 3-year study period.Analyses of these plots show that the MODIS-GLI AOT values are systematically overestimated compared to ground data, while the MODIS-NASA AOT values show a better match with these ground data, as shown by the higher correlation coefficient compared to that of MODIS-GLI.In both MODIS-GLI and MODIS-NASA cases, the dispersion around the 1 : 1 correlation line tends to increase at large AOT values (above 1).

Evaluation of Satellites' Derived Aerosols Data
For the AE, the graphs of Figure 3 show that the satellites' retrievals from both algorithms are clearly underestimated.This underestimation is slightly larger with the MODIS-GLI data than with the MODIS-NASA, and its magnitude seems to vary with the proportion of small/large particles' size of the aerosol load present in the ground data.The distribution of ground AE data appears as two overlapping datasets (around a cutoff value of 0.8) and, the increase in the number of ground data with AE > 0.8 against those with AE < 0.8 leads to an increasing underestimation of the AE by the satellite retrievals.In the large particle group, the satellite particles are smaller in size than the ground data, while in the small particle group, the satellite data are larger in size.And as, in general, small particles tend to be in higher number than large particles, among the ground data, the AE would be underestimated.The determination coefficient in the AE comparison is lower than that in the AOT comparison for both MODIS-GLI and MODIS-NASA cases.Among the 3 ground sites, the best match between the ground and both satellites' algorithms AOT and AE is at Fukue (northernmost site).
Table 1 shows the discrepancies between the ground measurements and the satellites' retrievals, expressed as the simple difference and the RMSE (root mean square error) of the AOT and AE.At the 3 sites, there are higher AOT RMSEs and biases (against ground observations) with the MODIS-GLI algorithm (ranges 0.214-0.234and −0.194-−0.203,resp.)compared to the MODIS-NASA algorithm (ranges 0.089-0.147and −0.008-0.009,resp.).The AOT bias of the MODIS-GLI algorithm is consistently negative, that is, the average AOT from this algorithm is overestimated (by an average of around 0.2), while the bias of the MODIS-NASA shows only a slightly under/overestimated value of around ±0.010.These AOT validation biases confirm the conclusions of Remer et al. [9] indicating that the MODIS-NASA ocean algorithm has virtually no offset and little bias except for a slight underprediction at high AOT, and the linear regression line follows the 1 : 1 line closely.The uncertainty in the AE satellite estimation is expected to be larger than that of the AOT.An expression of this uncertainty is [31] is the Ångström exponent, and  1 and  2 are the aerosol optical thickness at wavelengths  1 and  2 ; Δ, Δ 1 , Δ 2 are the errors in the Ångström exponent and the aerosol optical thickness at the wavelengths  1 and  2 , respectively.As a rule, satellite retrievals of aerosol optical thickness are more robust than the corresponding retrievals of aerosol size; consequently, the fraction of the fine mode contribution to the total optical thickness magnifies relative errors in the retrievals of individual path radiances by up to 4 to 8 times the original percentages [9].The AE uncertainty values are much closer than the AOT errors; between the two satellites' algorithms, the RMSE and bias vary from 0.414-0.566and 0.145-0.248,respectively, for the MODIS-GLI, then 0.321-0.388and 0.065-0.150,respectively, for the MODIS-NASA (Table 1).Both algorithms show consistently a positive bias in the AE, that is, smaller AE (larger size of particles) than ground data. of the ground data with however a consistently low sensitivity to low AOT values (<0.1).The principal mode occurs at the same ranges as that of the ground data, that is, around 0.2-0.3.The difference in the data distribution between this algorithm and the ground AOT distribution at most of the ranges is below 5%.Contrary to the MODIS-NASA algorithm, the MODIS-GLI algorithm shows its highest AOT peak range at >0.7, contrary also to the ground data peak, which is around 0.3.This MODIS-GLI higher range peak depletes the number of data in the lower AOT ranges and contributes to the AOT overestimation by this algorithm.Concerning the AE, it shows a unimodal distribution (mode at 1.2) at the Fukue ground site while the Hedo and Miyako ground data distributions appear to be bimodal.This distribution contrasts with that of the two satellites' algorithms, as they show at all sites a unimodal distribution with the modes consistently at 0.8 (for the MODIS-GLI) and 1.0 (for the MODIS-NASA), and fewer data are available beyond these modes.This contributes to the underestimation of the AE in the high AE ranges by the 2 satellites' algorithms, overwhelming the overestimation in the low AE ranges.
The comparison of the AOT, then AE seasonal distribution between the ground and satellite data, is shown at Figure 5, as the monthly average of each parameter for the 3year study period.The MODIS-GLI algorithm overestimates the AOT nearly equally each month.A slight inflection of this overestimation is seen in summer, around July.There is a better match each month between the MODIS-NASA and the ground data.
For the AE, winter months (January, February, March) and summer months (July, August, September), respectively, for the northern site (Fukue) and the southern sites (Hedo and Miyako), show closer values between the ground and both the MODIS-GLI and MODIS-NASA data.The AE is mostly underestimated by the satellites' algorithms.But this underestimation is generally within the ground data monthly variation range (standard deviation represented by vertical blue bars on the map).

Impact of External Factors and Aerosol Model Assumptions on Satellites' Retrievals.
As mentioned before, a prior step in the retrieval of aerosols properties is the choice of an aerosol model.This model is based on various assumptions.The uncertainties in these assumptions will likely affect the accuracy of the aerosols properties retrieved from satellite observations.Some of these uncertainties are the particle size, refractive index, the single scattering albedo, the particle sphericity [30], the satellite sensors' calibration, the cloud signal (thin clouds or adjacent cloud to the aerosols) interference, the asymmetry factor, and so forth.In this section we will review the contribution of these uncertainties on the accuracy of the satellite aerosols retrievals.Initial evaluations of Terra-and Aqua-MODIS ocean aerosols products, made by Kaufman and Tanré [30], suggest that the expected uncertainty for MODIS small AOTs ranges from 0.05 to 0.1 and increases to 20-30% for high AOTs.In Remer et al. [9], the uncertainty is given as ±0.03 ± 0.05 AOT.

Cloud Contamination Effect.
Atmosphere areas where aerosols properties are retrieved are supposed to be cloudfree.Examination of Terra-and Aqua-MODIS cloud masks suggest that, even after rigorous cloud screening, the resulting clear-sky dataset could still be contaminated (sometimes by 15% or more) by small clouds [32,33].Also, in cloud-free areas, adjacent cloud pixels to clear-sky areas can increase (brightening effect) or decrease (shadowing effect) the intensity of the aerosol radiative signal.The brightening effect is expected to play a more prominent role in these aerosols areas.To discuss the possible effects of cloud contamination on the satellite retrievals, we will examine the variation of an independent variable, the average cloud fraction (CF) in each spatial subarea (the CF is from the MOD/MYD06 version 5.1 MODIS cloud product) where the MODIS-NASA algorithm retrieval was performed, against the change in AE and AOT (both at the ground and satellite levels) then the bias of the satellite retrievals, that is, the difference (Ground-MODIS-GLI or Dif GLI and Ground-MODIS-NASA or Dif NASA). Figure 6 shows the results of these relationships.On this figure, the CF data are sorted in 0.1 ranges (i.e., 10% CF).The maximum range is >0.7 (i.e., >70% CF).Let us examine first the possibility of cloud contamination on the AE.The AE ground variation shows that it is nearly steady (no specific trend) along the CF variation ranges, while both the MODIS-NASA and MODIS-GLI AE show a decreasing trend (larger sizes) with the CF increase.Consequently, the AE difference ground-satellite (Dif GLI and Dif NASA) tends to increase with the CF increase.A stronger AE decreasing trend at high CF is observable with the MODIS-GLI than the MODIS-NASA at all sites.These trends show that both satellites' algorithm retrieved AE appear to be contaminated by clouds.
The AOT as the AE shows no specific trend in the ground data variation with the CF.In contrast, the MODIS-GLI AOT gradually increases with the CF increase while the increase of the latter with the MODIS-NASA AOT is slighter.Consequently, the difference in absolute value between the ground and satellite AOT (Dif GLI for the difference with MODIS-GLI or Dif NASA for the difference with MODIS-NASA) tends to increase, substantially (Dif GLI) or slightly (Dif NASA) with the CF increase.The AOT bias observed here may surpass 0.3 (e.g., at the CF peak of 0.7).The slight trend between the Dif NASA and the CF, found in this study confirms the remarks of Remer et al. [9] that the MODIS monthly mean AOT over the ocean is not significantly cloud contaminated (based on the good agreement between the MODIS-NASA ocean algorithm and 4 AERONET sites' AOT).
For the 3-year dataset, the average cloud amount is higher in the south (58.5% and 60.2% at Hedo and Miyako, resp.)than the north (44.2% at Fukue).And, coincidently, the cloud contamination and the negative bias (Table 1) are lower at the latter area compared to the former areas.
The AE and AOT larger biases in the MODIS-GLI than the MODIS-NASA, in relation to the cloud contamination, can be explained by the pixel size selection scale differences between the two algorithms.The MODIS-GLI algorithm retrieves aerosols properties on a pixel-by-pixel basis while the MODIS-NASA algorithm retrievals are made on 10 × 10-pixel boxes (only pixels with the highest noncloudy probability are selected and averaged).This selection process unfortunately reduces the number of data to be analyzed but produces less cloud-contaminated products compared to the pixel-by-pixel retrieval approach adopted by the MODIS-GLI algorithm.Kaufman and Tanré [30] suggest that when the aerosol product is not given on a pixel by pixel basis (0.5 × 0.5 km 2 ) but rather over a grid of 10 × 10 km 2 or 50 × 50 km 2 , clear sky areas used for the determination of aerosols can lead to a cloud contamination of aerosols and therefore contribute to measurement errors of these aerosols properties.Decreasing the spatial resolution for the aerosol product could help to reduce cloud contamination but also enhance the signal-to-noise-ratio (SNR) for the AE retrieval [30,34].

Particle Size and Shape Distribution Assumption.
In order to convert the aerosol mass measurements to size distributions, a model of the aerosol size distribution as a function of aerosol mass is needed [35].If "actual" aerosol properties are too different from those assumed by the model, then the satellite retrieved products will be less accurate [36].For the 3-year study period and at the 3 ground sites examined the distribution is consistently bimodal: small particle mode (accumulation) and large particle mode (coarse) radii.It was noticed that the small mode radius tends to move to higher values with the increase of the particle volume.This creates a larger dispersion (as will be shown later with the standard deviation) of the mode.The average radii of the fine mode ± standard deviation (STD) are 0.15 ± 1.75 m, 0.16 ± 1.63 m, and 0.16±1.72 m and the coarse modes 4.18 ± 2.35 m, 4.92± 2.27 m, and 4.72 ± 2.21 m for Fukue, Hedo, and Miyako, respectively.The accumulation mode radius is nearly similar at the 3 sites (0.15-0.16), and the coarse mode is larger in the south than the north.The STDs within each mode range are also very comparable.The accumulation and coarse modes radii used in the MODIS-GLI algorithm are, respectively, 0.17 ± 1.3, 3.44 ± 2.37.These values are, respectively, higher and lower than those of the ground sites, while their STDs Figure 6: Cloud fraction (CF) variation against aerosols properties retrievals from ground (AE and AOT ground), MODIS-GLI (AE and AOT GLI), and MODIS-NASA (AE and AOT NASA).Then CF variation against satellite bias, that is, difference ground-satellite (ground AE or AOT-MODIS-GLI AE or AOT), (ground AE or AOT-MODIS-NASA AE or AOT), during the 3 years of observations.The AE or AOT scales are on the right while the scale of the Dif is on the left.

AE and Dif
Accumulation Figure 7: AE and AE difference (ground-satellite) against the accumulation mode volume ratio (  /  ) variation.
are lower and close to those of the ground observations.The MODIS-NASA algorithm uses a combination of models whose median radius ± STD varies from 0.035 ± 0.40 to 0.10 ± 0.60 for the accumulation mode and 0.40 ± 0.60 to 1.00 ± 0.80 for the coarse mode, in the wavelengths ranges 0.47-0.86m.These particle size mode values and STDs are lower than those of the corresponding ground data.Single process originating particles have smaller STDs, while multiple processes originating particles have wider log-normal distributions [37].The smaller accumulation mode STDs of both satellites' algorithms compared to the ground data may imply that the aerosols models retrievals do not take into account most of the processes at the origin of these particles.As a consequence, the particle size of the satellite retrieved aerosols may differ from that of the ground.In some areas like semiarid areas, the effective radius of the accumulation mode particle may increase with increasing AOT [38].This tendency was also noticed in the ground data of the present study.However no relationship could be established between the increase of ground-satellite discrepancies with either the accumulation or the coarse mode radii.Figure 7 shows the histograms of the AE and AE difference (ground-satellite) against the ratio accumulation (small) mode volume against large mode volume (  /  ).The AE difference tends to increase with the   /  increase.As noticed before both satellites' algorithms use smaller STDs at the accumulation mode than the ground observed STDs, though the mode radius itself is not very different from that of the ground.The STD of the accumulation mode may increase the uncertainty in the satellites retrieved AE by both algorithms.
The particle shape (nonsphericity) can also be a potential contributor to the inaccuracy of satellite aerosols retrievals.According to Levy et al. [39] the MODIS-NASA ocean algorithm does not perform well in dust-laden atmosphere, and this problem is attributed to poor assumptions of the nonsphericity of the particles in the dust-aerosol phase functions.The peak in dust (where particles are less spherical) amount occurring in the aerosol load of the study area during spring does not particularly increase the difference between the ground AOT or AE with the satellite derived data (as shown by Figure 5).Also the 1 : 1 good alignment of the AOT of the MODIS-NASA algorithm with ground data does not show a specific seasonal dependency.All of these may lead to the conclusion that there is a less significant impact of the sphericity in the satellites' algorithms examined compared to other algorithms.Furthermore, tests conducted with a nonspherical model, by Nakajima and Higurashi [23], on the MODIS-GLI algorithm suggest that the nonspherical scattering effect is limited in this algorithm.

Impact of the Refractive Index.
In a given location, atmospheric aerosols are characterized by their concentration, their size distribution, their shape, their vertical profile, and by their chemical composition [2].The latter is indicative of the magnitude of the real and imaginary parts of the refractive index (RI) of aerosols.The real part of the RI (RIr) is spectrally related to the scattering of the radiation while the imaginary part (RIi) is a consequence of the absorption of radiation, with a null value meaning the particle does not absorb radiation.In the satellite aerosols retrieval scheme, the aerosol model assumes specific realistic values for the real and imaginary parts of the RI.In the wavelength range 0.47-0.86m, the MODIS-GLI algorithm uses aerosol models with a single RI (1.5-i0.005),based on field data analyses.The real part RI of this algorithm is close to general tropospheric aerosols indices [40].Meanwhile, the MODIS-NASA algorithm uses a combination of models where the RIs are set to vary according to the area and the volume size distribution mode ratio.The values used by this algorithm are derived from the sun/sky ground photometers of AERONET [9], 1.45-i0.0035or 1.40-i0.0020for the fine mode and 1.35-i0.001for the coarse mode, of all aerosols types except dustlike types where the RIs are 1.53-i0.003and 1.53-i0.000for 0.47 m and 0.86 m wavelengths, respectively [41].The average RIs at 0.5 m wavelength obtained from ground data for the whole study period are 1.44-i0.0043,1.43-i0.0016,1.44-i0.0029for Fukue, Hedo, and Miyako sites, respectively.The MODIS-NASA real and imaginary parts are closer or within the ranges of the ground RI values than the MODIS-GLI are to the latter.The following distribution of RI at the ground sites was observed (in brackets are the percentage of data at Fukue, Hedo, and Miyako, respectively, for the RI range considered): the real part RI: <1.40 (33%, 41%, 36%), 1.4-1.5 (37%, 29%, 32%), >1.5 (30%, 30%, 32%); and the imaginary part RI: <0.002 (54%, 84%, 75%), 0.002-0.005(22%, 7%, 10%), >0.005 (24%, 9%, 15%).Only one-third of the data are beyond RIr 1.5, and less than 25% have RIi above 0.005 (RI values considered for the MODIS-GLI).
The comparison with the AOT shows that (Figure 8) the RIr tends to increase with the AOT increase, for the ground and both satellites' algorithms data.The AOT RMSE is high at high RIr and vice versa for low RIr values such as 1.4 (closer to the values used in the MODIS-NASA algorithm).For the MODIS-GLI, the AOT RMSE is high in the 1.4-1.45range and tends to decrease beyond 1.5 (closer to the values used by the MODIS-GLI algorithm).All these variations show that the RIr is an important contributor to the enhancement of the satellite AOT uncertainty.No trend was seen between the RIi and the AOT RMSE nor between the RIr or RIi and the AE RMSE variation.4.2.4.Impact of Scattering Parameters.The main property representing the scattering by aerosols is generally expressed by the single scattering albedo (SSA).The scattering preferential direction may be either forward (+1 for pure forward scattering) or backward (-1 for pure backscattering).The proportion of light scattered in the forward direction is defined as the asymmetry factor (AF).For satellites' aerosol retrievals, the SSA and AF are precomputed using a Miebased theory (spherical particles) scattering model.Irregularly shaped tropospheric particles may have larger SSA and smaller AF than their equal volume spherical counterparts [40].The average SSA ± STD and AF ± STD observed by the SKYNET ground sites are, respectively, 0.95 ± 0.06 and 0.69 ± 0.04 for Fukue, 0.98 ± 0.05 and 0.70 ± 0.04 for Hedo, and 0.96 ± 0.07 and 0.71 ± 0.05 for Miyako.Figure 9 shows the variation of the AF against the AE and the AE Dif.The ground AE consistently decreases with the increase of the AF at the three ground sites; that is, the increase in the forward scattering of light results from the increase of the particle size.No such clear trend is observed between the ground AF and both, the satellites' AE of the GLI and NASA algorithms.As a result, the ground-satellite AE discrepancies decrease as the AF increases that is, as the scattering of light in the forward direction increases.The uncertainty in the determination of the AF appears then to impact both satellites' algorithms AE retrievals.The ground AF slope change with the AF is not as steep in the northern Fukue as in the southern Hedo and Miyako.This is consistent with the increase in coarse particles from the north to the south (increased forward radiation scattering) and consequently the AE RMSE.No consistent correlation was found between the SSA and either, the ground AE, the AOT or the error in the determination of the satellite retrieved AOT and AE.This could be related to the fact that the influence of error in the retrieved SSA is important only for thick layers [36].And the data used here are mostly thin aerosols layers.

Impact of Solar and Satellite Zenith
Angles.An ideal aerosol algorithm would retrieve AOTs of equal quality, independent of solar and observing geometry [42].No consistent trend showing possible or significant satellite bias was observed between both the AOT and AE (AOT and AE RMSE) and the solar zenith angle or the satellite zenith angle.Using collection 5 MODIS dark target data over land, Levy et al. [42] suggest that there is only a slight over-and underestimate of the AOT by ∼0.01, respectively, on the sunglint side and the sun-shadow side of the MODIS swath.

Implications for Future Satellites' Algorithms
The RMSE and Dif for the AOT retrieved from the MODIS-GLI algorithm are higher compared to those of the  MODIS-NASA algorithm.The accuracies of both algorithms for the AE are much closer to each other as well as lower than the AOT accuracies.Some of the factors explaining these accuracy differences range from the cloud contamination, the particle size and shape distribution, the refractive index, and the scattering parameters to the sun-satellite geometry.The lessons drawn from the possible impact of these parameters on satellite retrievals will help to improve the accuracy of future satellites' aerosols products such as those from the coming GCOM-C/SGLI satellite.The cloud contamination was found to be stronger in the AOT of MODIS-GLI than that of MODIS-NASA.This was explained by the differences in the spatial resolution of the retrievals: pixel-by-pixel retrieval in the MODIS-GLI and pixels' selection retrieval from boxes of 10 × 10 pixels in the MODIS-NASA.To avoid this contamination in the future and in case the pixel-by-pixel approach is chosen, a cloud fraction analysis should accompany the aerosols retrievals and be used to correct the aerosols properties' inaccuracies.The MODIS-NASA approach looks attractive as retrievals are made only on the high probability clear pixels of a wider pool of pixels, but the inconvenience is the lower amount of output data.Changes in aerosols mixtures of the region would need a better representation.For example, this study showed that the aerosol particle radii deviations from the mode should be wider in order to integrate the dynamic processes at the origin of the aerosols of this region.Also, a wider range (than that already used by the MODIS-NASA algorithm) for the refractive indices would be necessary to reduce the uncertainty in the AE.The high spatial and spectral resolution of the GCOM-C/SGLI satellite (e.g., 250 m at all visible and thermal infrared channels) should help to better discriminate clouds from cloud-free areas than actual satellites do and therefore reduce possible cloud contamination of the aerosols products.The improvement in the AOT accuracy will certainly increase that of the AE.The increase in the number of ground sites measurements of aerosols such as in the SKYNET, AERONET, and other networks, providing better quality data and quantitative characterization of aerosols, would be of great asset for the improvement of the assumptions used in aerosols models: predetermined particle size distribution (modes' radii and volumes), refractive index, asymmetry factor, single scattering albedo, ambient humidity and particle growth, and so forth.

Conclusion
The quality of aerosols retrievals from future satellites was discussed in this study, through a validation study of aerosols products from actual satellites (Terra-and Aqua-MODIS) having compatible channels with those of the coming GCOM-C/SGLI.Aerosols products derived from two algorithms (MODIS-GLI and MODIS-NASA) using these satellites were evaluated against ground-truth data from 3 SKYNET sites located in the Pacific East Asian region.The results obtained showed a systematic overestimation of the AOT by the MODIS-GLI and a better agreement between the ground and MODIS-NASA AOT and AE underestimation by both algorithms.The large particle AE overestimation was overwhelmed by the small particle underestimation.The seasonal variation of the aerosols parameters (AOT and AE) based on monthly averages showed similar trends/variations though the direct values differed, sometimes with MODIS-NASA or more often with MODIS-GLI.Also nearly constant monthly AOT differences between the ground aerosols data and the satellite retrievals were observed.The AE showed closer values between the ground and satellite data during winter and spring for the Fukue site and summer and autumn for the southern sites (Hedo and Miyako).The effect of cloud contamination on the satellites retrievals, examined through the use of an independent variable, the cloud fraction, showed that cloud contamination was more pronounced on the AOT MODIS-GLI than that of MODIS-NASA.The AE analysis showed that this parameter decreased (particles get larger) with the CF for both algorithms.The difference in the AOT contamination between the two satellites' algorithms was found to be mostly due to the retrieval pixel resolution difference: pixel-by-pixel for the MODIS-GLI, therefore higher probability of contamination compared to the pixels' selection from larger boxes for the MODIS-NASA.Also, this difference could be explained by the fact that the MODIS-NASA uses wider spectral information to estimate the aerosol model and, may perform more accurate retrieval of spectral AOT [36].The AE distribution showed that when the ground data have a bimodal distribution, the uncertainty in the AE retrievals by the two satellites' algorithms is higher.The shape of the particle size distribution examined showed larger STDs than those considered by the satellites' algorithms, leading to limited mixtures in the aerosols sizes retrieved by the satellites' algorithms.Another contributor to the ground-satellite differences was the uncertainty on the satellite retrievals, due to the asymmetry factor (AF).A negative correlation was found with both the ground AE and the AE difference (ground-satellite) against the AF.The increase of the light scattering in the forward direction is a consequence of the increasing size of the aerosols particles.The real part of the refractive index was found to be correlated with the AOT RMSE of both satellites' algorithms and not with the AE RMSE.To end this study, some measures aimed at improving the accuracy of future satellites' algorithms were discussed.These involve remedies to cloud contamination and corrections of the particle size distribution and scattering models.

4. 1 .Figure 1 :
Figure 1: Geographical location of the three SKYNET radiation sites (Fukue, Hedo, Miyako) used in this study to evaluate the satellite retrievals, over an image of satellite derived aerosols properties from the MODIS-GLI algorithm, for the month of January.

Figure 2 :
Figure2: Scatter diagrams between ground (Fukue, Hedo, Miyako) and satellites' derived aerosols AOT.Ground data are on the -axis while satellite derived data (from both algorithms, MODIS-GLI and MODIS-NASA) are on the -axis.AOT is the aerosol optical thickness.

Figure 4 Figure 3 :
Figure 3: Scatter diagrams between ground (Fukue, Hedo, Miyako) and satellites' derived aerosols AE.Ground data are on the -axis while satellite derived data (from both algorithms, MODIS-GLI and MODIS-NASA) are on the -axis.The AE is the Ångström exponent.

Figure 4 :
Figure 4: Comparison between the frequency distribution of ground measurements and satellite derived (GLI for the MODIS-GLI algorithm, NASA for the MODIS-NASA algorithm) aerosols properties (AE (a), AOT (b)).

Figure 5 :
Figure 5: Average monthly variation of the ground (Ground) AOT at the 3 SKYNET sites (Fukue, Hedo, Miyako) and the corresponding satellite retrievals: MODIS-GLI (GLI) and MODIS-NASA (NASA).The errors bars (blue) represent the standard deviation of the monthly ground data.

Figure 8 :
Figure 8: Real part of the refractive index (RIr) variation at the ground against ground AOT, and the satellite retrieved AOT and AOT RMSE.

Figure 9 :
Figure 9: Aerosol asymmetry factor (AF) histogram distribution against AE and the satellites' algorithms AE, and the difference groundsatellite.

Table 1 :
Error analyses showing the root mean square error (RMSE) between the ground data and the satellite derived data for the MODIS-GLI algorithm (GLI RMSE) and the MODIS-NASA algorithm (NASA RMSE), and the simple differences between the ground data and the MODIS-GLI (GLI Dif), then the MODIS-NASA (NASA Dif) satellites' algorithms.