Primary Evaluation of the GCOM-C Aerosol Products at 380 nm Using Ground-Based Sky Radiometer Observations

: The Global Change Observation Mission-Climate (GCOM-C) is currently the only satellite sensor providing aerosol optical thickness (AOT) in the ultraviolet (UV) region during the morning overpass time. The observations in the UV region are important to detect the presence of absorbing aerosols in the atmosphere. The recently available GCOM-C dataset of AOT at 380 nm for January to September 2019 were evaluated using ground-based SKYNET sky radiometer measurements at Chiba, Japan (35.62 ◦ N, 140.10 ◦ E) and Phimai, central Thailand (15.18 ◦ N, 102.56 ◦ E), representing urban and rural sites, respectively. AOT retrieved from sky radiometer observations in Chiba and Phimai was compared with coincident AERONET and multi-axis di ﬀ erential optical absorption spectroscopy (MAX-DOAS) AOT values, respectively. Under clear sky conditions, the datasets showed good agreement. The sky radiometer and GCOM-C AOT values showed a positive correlation ( R ) of ~0.73 for both sites, and agreement between the datasets was mostly within ± 0.2 (the number of coincident points at both sites was less than 50 for the coincidence criterion of ≤ 30 km). At Chiba, greater di ﬀ erences in the AOT values were primarily related to cloud screening in the datasets. The mean bias error (MBE) (GCOM-C – sky radiometer) for the Chiba site was − 0.02 for a coincidence criterion of ≤ 10 km. For a similar coincidence criterion, the MBE values were higher for observations at the Phimai site. This di ﬀ erence was potentially related to the strong inﬂuence of biomass burning during the dry season (Jan–Apr). The diurnal variations in AOT, inferred from the combination of GCOM-C and ozone monitoring instrument (OMI) observations, showed good agreement with the sky radiometer data, despite the di ﬀ erences in the absolute AOT values. Over Phimai, the AOT diurnal variations from the satellite and sky radiometer observations were di ﬀ erent, likely due to the large di ﬀ erences in the AOT values during the dry season.


Introduction
Atmospheric aerosols play a crucial role in controlling Earth's radiation budget and thus impact regional climates and hydrological cycles. High near-surface concentrations of aerosol particles impact air quality and human health [1,2]. Global characterization of aerosols (and their effects on the climate) is difficult due to large spatial and temporal variability in their abundance and properties [3]. Satellite observations of aerosols are useful for understanding the large variability of aerosols in space and time [4][5][6]. Since 1990, a variety of sensors have been launched on Terra, Aqua, Aura, CALIPSO, Table 1. The datasets used in the study and the uncertainties in the retrieved product.

GCOM-C Data
GCOM-C carries a multi-spectral optical sensor called the second-generation global imager (SGLI), which is the successor to the Global Imager Onboard Advanced Earth Observing Satellite-II [23]. The SGLI consists of two radiometers, the visible and near-infrared radiometer (VNR) and the infrared scanner (IRS), covering the radiation observations over the wavelength range from near-ultraviolet (UV) to infrared. The SGLI-VNR observes polarized and non-polarized radiance. The VNR uses a wide-swath (1150 km) push-broom scan with a line CCD detector. On the other hand, the IRS uses the conventional cross-track mirror scan system with a 1400 km swath width. SGLI has 19 channels, including two polarization channels in the visible and near-infrared (NIR) regions [24,25]. The specifications of the SGLI channels are given in Table 2. The relative spectral response of the GCOM-C VNIR channels is shown in Figure 1.
Remote Sens. 2020, 12   The effectiveness of polarization and UV wavelengths for aerosol data retrieval was demonstrated by the polarization and directionality of Earth's reflectances (POLDER) and total ozone mapping spectrometer (TOMS) instruments, respectively [26,27]. GCOM-C utilizes both the  The effectiveness of polarization and UV wavelengths for aerosol data retrieval was demonstrated by the polarization and directionality of Earth's reflectances (POLDER) and total ozone mapping spectrometer (TOMS) instruments, respectively [26,27]. GCOM-C utilizes both the polarization channel and near-UV data to retrieve aerosol information from SGLI observations.
To maintain data observation quality and estimate the instrumental error, the on-orbit and calibration maneuvers are performed periodically. Table 3 shows the type of calibration performed for the SGLI instrument. A more detailed description of the SGLI instrument and calibration can be found in [23,28,29].

GCOM-C Aerosol Retrieval
GCOM-C AOT data at 380 nm was obtained from the GCOM-C standard data platform (https://gportal.jaxa.jp/gpr/?lang=en). GCOM-C has 19 channels, including near-UV (380 nm) and violet (412 nm) wavelengths and polarization channels in the red (670 nm) and NIR (865) regions. All channels from UV to NIR without significant gas absorption are used for retrieval of aerosols over land. Only the channels at wavelengths longer than 800 nm are used for aerosol retrieval over the ocean. Different algorithms are applied to the polarization and non-polarization channels. Only the retrieval of aerosols in the non-polarization channel is discussed here.

(a) Pre-retrieval calculated parameters
At first, simulated top of atmosphere (TOA) radiances are calculated using the radiative transfer model STAR [30][31][32]. The TOA reflectance at a particular channel i is approximated by the following equation: where ρ sim i is the TOA reflectance at the channel i, ρ a i is the atmospheric path reflectance, t s i is the total transmittance from the sun to the surface, t v i is the total transmittance from the surface to the sensor, ρ s i is the surface reflectance, s i is the spectral albedo for the illumination of the atmosphere from the ground, θ o is the solar zenith angle, θ is the satellite zenith angle, and ϕ is the sun/satellite relative azimuth angle. Secondly, the parameters in Equation (1) are precomputed for each candidate aerosol model. The aerosol models were assumed as an external mixture of the fine and coarse particles. The fine-mode aerosol and dust particle properties were set according to Omar et al. [33], based on global observations of aerosols from the AERONET [17] network. The coarse-mode model settings were based on the study of Sayer et al. [34]. The monomodal lognormal volume size distribution (r d ) used for the aerosol size of the fine and coarse mode models was where C v is the particle volume concentration, r v is the volume median radius, and σ is the standard deviation. The r v (σ) values were set to 0.143(1.537), 2.59(2.054), and 2.834(1.908) for fine, coarse marine, and coarse dust, respectively. The aerosol shape was assumed to be spherical in the fine and coarse marine model, whereas the non-spherical shape was assumed in the coarse dust model. The non-spherical parameters of the dust model were calculated based on the work of Nakajima et al. [35]. The real part of the refractive imaginary index for fine, coarse marine, and coarse dust were set to 1.439, 1.362, and 1.452, respectively. The imaginary part of the refractive index (m i ) for coarse marine and dust was set to 3 × 10 −9 and 0.0036, respectively, whereas the m i for the fine mode aerosols was perturbed to account for absorbing and non-absorbing aerosols. The aerosol models are summarized in Table 4.

(b) Retrieval of aerosol optical and physical properties
In the first step of the retrieval, the clear-sky pixel was selected initially using the cloud detection algorithm developed by Ishida et al. [36,37]. Secondly, the precalculated TOA reflectance was corrected for gas absorption, specifically for ozone and water vapor. Total ozone columns from OMI and water vapor columns from the Japan meteorological agency (JMA) global analysis dataset were used. Finally, the optimal estimation method [38] was used to retrieve the aerosol parameter. The outputs of the retrieval algorithm were AOT, Angstrom exponent (AE), single scattering albedo (SSA), and quality assurance (QA) flag. The data were screened based on quality and cloud flags included within the dataset. Only data from clear days (based on the GCOM-C cloud flag) were used for comparison with sky radiometer observations. The product version 1 was used in the current analysis. The details of the non-polarization aerosol product are explained in the work of Yoshida et al. [21]. Moreover, a detailed Remote Sens. 2020, 12, 2661 6 of 21 description of the retrieval algorithm is provided in the algorithm's theoretical basis document available in the following link https://suzaku.eorc.jaxa.jp/GCOM_C/data/product_std.html.

OMI Aerosol Product
OMI-observed near-UV radiances at 354 and 388 nm and the inversion algorithm named OMAERUV (OMI near UV aerosol retrieval algorithm) were used to derive the aerosol optical depth and single scattering albedo at 388 nm. The Lambertian equivalent reflectivity at 388 nm, the difference between reflectivity at the TOA and surface, and aerosol index were used in the OMAERUV algorithm [8]. The details of the OMI aerosol products have been documented in the work of Torres et al. [39,40], Jethva and Torres [41], and Ahn et al. [16]. The OMI aerosol product at 388 nm used in this study was obtained from the NASA Goddard Earth Sciences and Information services center server (https://avdc.gsfc.nasa.gov/pub/data/satellite/Aura/OMI/V03/L2OVP/OMAERUV/txt/ (last browsed on 10 January 2020). The OMI overpass time occurred at 13:30 LT.

Sky Radiometer Observations of Aerosol Optical Properties
SKYNET sky radiometers conducted observations of direct radiance, angular sky radiance, and zenith sky radiance at 11 wavelengths with a temporal resolution of 10 min. The direct and diffuse radiances were measured within 160 • of the center of the sun with a field of view of 1 • . The AOT, single scattering albedo (SSA), and refractive indexes at 340, 380, 400, 500, 675, 870, and 1020 nm were retrieved using the sky radiometer analysis package, Center for Environmental Remote Sensing (SR-CEReS, version 1) [42]. All pre-and post-processing was conducted using SKYRAD.pack version 5 [43], implemented in SR-CEReS for near-real-time data output. SKYNET sky radiometers used an on-site calibration method, called the improved Langley (IL) method [44], to determine the calibration constant (Fo). The solar disk sextant (SDS) method [44][45][46] was employed to calculate the solid viewing angle (SVA). The cloud screening algorithm of the sky radiometer data consists of three steps: (1) test with global irradiance data, (2) spectral variability test, and (3) statistical analyses test. The first test detects cloud affected data, and the spectral variability algorithm is applied to detect clear sky data. The final test detects outliers (if any) from clear sky data detected by the previous two steps. The detailed cloud screening algorithm is explained in the work of Khatri and Takamura [47]. The number of clear sky and cloudy data points retrieved from the sky radiometer observations at both sites are summarized in Table 5. Retrievals from SR-CEReS were validated during the NASA KORUS-AQ (Korea-US Air Quality) campaign in 2016 [42].

MAX-DOAS Observations of Aerosol and Trace Gases
The MAX-DOAS instrument is equipped with an ultraviolet-visible (UV-VIS) spectrometer (Maya2000Pro; Ocean Optics), which is located indoors. A telescope unit consisting of a single fixed telescope and a movable 45 • inclined mirror on a rotary actuator is located outdoors and conducts reference and off-axis observations. Scattered sunlight was measured at six ELs of 2 • , 3 • , 4 • , 6 • , 8 • , and 70 • every 15 min. The off-axis ELs were limited to below 10 • . This limitation is expected to minimize systematic error in the oxygen collision complex (O 4 ) fittings while maintaining high measurement sensitivity in the lowest layers of the aerosol and trace gas profiles [20]. High-resolution spectra were recorded from 310 to 515 nm (full width at half maximum of 0.4 at 386 and 476 nm). Wavelength calibration was performed using a high-resolution solar spectrum from [48]. The measured spectra were analyzed using the Japanese vertical profile retrieval algorithm version 2 (JM2) [20,49]. Aerosol information was retrieved from the O 4 absorption data. JM2 retrieves the vertical aerosol profile in two steps. First, O 4 differential slant column density (O 4 ∆SCD) at 357 and 476 nm are retrieved using the DOAS method [50], which is a nonlinear least-squares spectral fitting technique. Significant O 4 absorption values in the wavelength ranges of 338-370 and 460-490 nm were used to determine O 4 ∆SCD values. Second, the aerosol vertical profile is retrieved using O 4 ∆SCD values, and the differential air mass factor is computed (accounting for all ELs). The second retrieval step of JM2 is based on the optimal estimation method [38]. The retrieval procedures and error estimates are explained in detail by Irie et al. [20,49,51,52].
In addition to the aerosol products, glyoxal (CHOCHO) vertical column density (VCD) retrieved from the MAX-DOAS observations were also utilized. VCD and CHOCHO concentration data were retrieved using the JM2 algorithm. The fitting window of 436-490 nm was used for the CHOCHO retrievals. The CHOCHO retrieval procedure is explained in detail by Hoque et al. [22]. MAX-DOAS observations represented a horizontal resolution of~10 km.
To minimize the impact of clouds on the MAX-DOAS retrievals, data screening was performed. First, retrieved AOT greater than 3 (the highest value in the lookup table) was excluded. Such large AOT values potentially occur due to optically thick clouds. Further data screening criterion was based on the fitting residuals of O 4 ∆SCD and the trace ∆SCD (glyoxal slant column density (CHOCHO ∆SCD) for this study). The screening criterion was-(1) O 4 ∆SCD residuals less than 10%, (2) CHOCHO ∆SCD less than 50%, and (3) degree of freedom greater than 1.02. The details of the cloud screening procedure are explained in the work of Irie et al. [20,49].

AERONET Aerosol Data
AERONET is a network of automatic robotic sun and sky-scanning radiometers, which has more than 100 sites around the globe. AERONET data provides quality-assured aerosol optical properties to assess and validate satellite retrievals. The automated sun and sky scanning radiometers (CIMEL) measured direct sunlight with a 1.2 • viewing geometry, every 15 minutes at 340, 380, 440, 500, 675, 870, 940, and 1020 nm. It took~8 s to scan all eight wavelengths, with a motor-driven filter wheel positioning each filter in front of the detector. Then, the solar extinction measurements were utilized to retrieve the AOT values at every wavelength except for the 940 nm channel, which is dedicated to perceptible water content retrievals. Multiple algorithm parameters and statistical tests were utilized to minimize cloud impact on AERONET observations. The detailed cloud screening procedures of the AERONET observations are explained by Giles et al. [53]. The details of the AERONET retrievals are explained in the work of [17].

Sky Radiometer and AERONET Observations in Chiba
The consistency of the sky radiometer AOT values in Chiba was checked using an independent aerosol dataset from the coincident AERONET observations. AERONET observations were also used as additional information to explain aerosol characteristics in Chiba. Figure 2 shows the comparison between the SKYNET and AERONET (level 1.5, Version 3.0) AOT values at 380, 500, and 675 nm for the observation period from January to December 2019. Only clear sky days (based on the cloud flag in both datasets) were plotted. The statistics of the comparison is given in Table 6. The datasets showed excellent agreement at all three wavelengths, with a correlation coefficient (R) of~0.99. Thus, the absolute values and temporal variations of AOT values from both sources were quite similar. At all three wavelengths, the datasets generally agreed within ±0.02. Greater differences in a few cases could be attributed to differences in the retrieval procedure and cloud influences on both datasets. Overall, the excellent agreement allowed confidence in the quality of the SR-CEReS retrieval products for evaluating GCOM-C observations. Remote Sens. 2020, 12, x FOR PEER REVIEW 8 of 22 in both datasets) were plotted. The statistics of the comparison is given in Table 6. The datasets showed excellent agreement at all three wavelengths, with a correlation coefficient (R) of ~ 0.99. Thus, the absolute values and temporal variations of AOT values from both sources were quite similar. At all three wavelengths, the datasets generally agreed within ±0.02. Greater differences in a few cases could be attributed to differences in the retrieval procedure and cloud influences on both datasets. Overall, the excellent agreement allowed confidence in the quality of the SR-CEReS retrieval products for evaluating GCOM-C observations.   Figure 3 shows the AOT at 476 nm retrieved from the sky radiometer and MAX-DOAS observations in Phimai. The sky radiometer AOT values at 500 nm were converted to those at 476 nm using the sky radiometer Angstrom exponent (AE) data. Coincident days with at least four observations were selected. The statistics of the comparison are shown in Table 7. The sky radiometer and MAX-DOAS AOT values at 476 nm showed good agreement with an R-value of 0.73. In some cases, larger differences were observed between the sky radiometer and MAX-DOAS values at 476 nm. The integration times of the two instruments differed. It was unlikely that the difference in integration time was the dominant factor driving the differences in the AOT values. Another possible cause was the impact of clouds. An example case using LIDAR observations in Phimai is shown in Figure 3.    Figure 3 shows the AOT at 476 nm retrieved from the sky radiometer and MAX-DOAS observations in Phimai. The sky radiometer AOT values at 500 nm were converted to those at 476 nm using the sky radiometer Angstrom exponent (AE) data. Coincident days with at least four observations were selected. The statistics of the comparison are shown in Table 7. The sky radiometer and MAX-DOAS AOT values at 476 nm showed good agreement with an R-value of 0.73. In some cases, larger differences were observed between the sky radiometer and MAX-DOAS values at 476 nm. The integration times of the two instruments differed. It was unlikely that the difference in integration time was the dominant factor driving the differences in the AOT values. Another possible cause was the impact of clouds. An example case using LIDAR observations in Phimai is shown in Figure 3.

Evaluation of GCOM-C Data
Scatter plots of AOT at 380 nm retrieved from GCOM-C and sky radiometer observations in Chiba and Phimai are illustrated in Figure 4. The days of overlap between GCOM-C and sky radiometer data were analyzed. Sky radiometer data within ±20 min of the satellite overpass time were selected. Satellite data within a 30 km radius of the observation sites were compared. The observation time and satellite observation radius were optimized to obtain sufficient coincident points for comparison. As the Chiba site is situated near the sea, only the data over land were selected using the land-sea mask information. The number of coincident points available for the Chiba and Phimai site was 47 and 31, respectively. radiometer data were analyzed. Sky radiometer data within ±20 min of the satellite overpass time were selected. Satellite data within a 30 km radius of the observation sites were compared. The observation time and satellite observation radius were optimized to obtain sufficient coincident points for comparison. As the Chiba site is situated near the sea, only the data over land were selected using the land-sea mask information. The number of coincident points available for the Chiba and Phimai site was 47 and 31, respectively. A good positive correlation (R) of 0.77 and 0.75 at Chiba and Phimai, respectively, was observed. Figure 4c,d show the absolute differences between GCOM-C and sky radiometer observations over both sites. At both sites, the AOT values mostly agreed within ±0.2, and the differences larger than ±0.2 were observed mostly at higher AOT values. One potential reason for these larger differences could be the influence of clouds. At Chiba, the difference ~0.39 was observed on April 22 (Table 8). The coincident sky camera images showed clouds around the GCOM-C overpass time on April 22 (http://atmos3.cr.chiba-u.jp/skyview/svchib201904.html; last accessed on February 13, 2020). Similarly, a few clouds were observed during the GCOM-C overpass time on March 8 (http://atmos3.cr.chiba-u.jp/skyview/svchib201903.html; last accessed on February 13, 2020). However, the agreement between the two datasets was very good on March 8 (Table 8). Although both datasets were subjected to cloud screening schemes, some cloud pixels might not have been removed, which could lead to the differences observed on April 22. A good positive correlation (R) of 0.77 and 0.75 at Chiba and Phimai, respectively, was observed. Figure 4c,d show the absolute differences between GCOM-C and sky radiometer observations over both sites. At both sites, the AOT values mostly agreed within ±0.2, and the differences larger than ±0.2 were observed mostly at higher AOT values. One potential reason for these larger differences could be the influence of clouds. At Chiba, the difference~0.39 was observed on 22 April (Table 8). The coincident sky camera images showed clouds around the GCOM-C overpass time on 22 April (http://atmos3.cr.chiba-u.jp/skyview/svchib201904.html; last accessed on 13 February 2020). Similarly, a few clouds were observed during the GCOM-C overpass time on March 8 (http://atmos3.cr.chiba-u.jp/skyview/svchib201903.html; last accessed on 13 February 2020). However, the agreement between the two datasets was very good on 8 March (Table 8). Although both datasets were subjected to cloud screening schemes, some cloud pixels might not have been removed, which could lead to the differences observed on 22 April. The mean bias error (MBE) values in Chiba and Phimai were -0.06 and −0.19, respectively, for the coincidence criterion of 30 km around the sites. Table 9 shows the estimated MBE, root mean square error (RMSE), slope, and R-values, at both sites, for GCOM-C data of various spatial resolutions. Table 9. The MBE, RMSE, slope, and R calculated for the different spatial resolution of the GCOM-C observations. The parameters were calculated for both wavelengths in Chiba and Phimai. The number of data points (N) for each coincidence criterion for both sites is also shown. At both sites, MBE didn't change significantly with changes in the spatial resolution of the satellite coincidence point. At Chiba, the R-values were almost similar for the satellite coincidence point less than 50 km around the Chiba site. This indicated the good performance of GCOM-C retrievals at 380 nm over heterogeneous surfaces. The MBE values at Phimai were higher compared to those of Chiba, irrespective of the GCOM-C coincidence criterion. The potential reasons for the MBE and differences between the datasets are discussed in the later sections.

Aerosol Composition in Chiba
The changes in the aerosol composition could impact the observed differences in the dataset. To investigate the aerosol types in Chiba, AE, AAE, EAE, and aerosol volume size distribution data were analyzed. AAE values were calculated following the procedure of Irie et al. [18], using the following equations ln[AAOD(λ)] = a − AAE × ln(λ) where a is the intercept in Equations (3) and (5). The wavelength ranges from 340-870 nm were used for the calculation. Aerosol volume size distribution data were obtained from the coincident AERONET sunphotometer measurements. The aerosol volume size distribution retrieved using SR-CEReS is not yet well understood. The coincident AERONET volume size distribution data were used as supporting information along with sky radiometer observations to infer aerosol composition in Chiba. Figure 5 shows AE, AAE, and EAE values, as well as the mean aerosol volume size distribution for the observation period in Chiba. The aerosol mean volume size distribution was shown for coincident days between the SKYNET and AERONET observations. The AE values are generally greater than 1, indicating that small particles are dominant. The mean size distribution also showed that fine-mode aerosols were dominant in Chiba. AAE values range between 1 and 2, with values between 1 and 1.5 representing dominant urban-industrial aerosols, larger values (~>1.5) representing biomass burning aerosols, and the largest values (~>2.00) representing dust aerosols [54]. that fine-mode aerosols were dominant in Chiba. AAE values range between 1 and 2, with values between 1 and 1.5 representing dominant urban-industrial aerosols, larger values (~>1.5) representing biomass burning aerosols, and the largest values (~>2.00) representing dust aerosols [54]. AAE values between 1.2 and 1.5 occurred with the highest frequency during the observation period in Chiba. The observed AAE values (1.2-1.5) were in the range reported by Russel et al. [54] for dominant urban-industrial aerosols. However, the upper limit of the AAE range (i.e., 1.5) was similar to the AAE value (1.57) for dominant biomass burning aerosols reported by Irie et al. [18]; however, while they reported values over their entire observation period, the values reported in this work were daily mean AAE values. On some days, AAE values were as high as ~2, which fell within the range of AAE values for dominant dust aerosols reported by Russel et al. [54]. Figure 5 shows a case (February 20) when the AAE value was 1.95. The aerosol volume size distributions, EAE values, and 48-h backward trajectories of air masses at the height of 200 m arriving in Chiba are depicted in Figure 5d,e. The backward trajectories were calculated using the HYSPLIT model [55]. The EAE value on February 20 was < 1, indicating the presence of larger particles in the air. Aerosol volume size distribution data also supported the presence of larger particles on that day (Figure 5c). Backward trajectories showed that air masses from the Changchun area in China, where sand fields are located [56], arrived at the Chiba site. Moreover, the weather in Harbin city (located near Changchun) is often affected by sand and dust from Northwest China and Inner Mongolia [57]. This result indicated that dust particles were transported to the Chiba site.
However, a high AAE value did not necessarily indicate the presence of dust aerosols in Chiba. The AAE, EAE, aerosol volume size distribution, and 48-h backward trajectories on March 12 are shown in Figure 5; on that date, the AAE and EAE values were 1.85 and 1.32, respectively. Besides, On some days, AAE values were as high as~2, which fell within the range of AAE values for dominant dust aerosols reported by Russel et al. [54].  Figure 5d,e. The backward trajectories were calculated using the HYSPLIT model [55]. The EAE value on February 20 was < 1, indicating the presence of larger particles in the air. Aerosol volume size distribution data also supported the presence of larger particles on that day (Figure 5c). Backward trajectories showed that air masses from the Changchun area in China, where sand fields are located [56], arrived at the Chiba site. Moreover, the weather in Harbin city (located near Changchun) is often affected by sand and dust from Northwest China and Inner Mongolia [57]. This result indicated that dust particles were transported to the Chiba site.
However, a high AAE value did not necessarily indicate the presence of dust aerosols in Chiba. The AAE, EAE, aerosol volume size distribution, and 48-h backward trajectories on March 12 are shown in Figure 5; on that date, the AAE and EAE values were 1.85 and 1.32, respectively. Besides, the aerosol volume size distribution on that date showed that fine-mode aerosols were dominant in Chiba. The trajectories also showed that no dust was transported to the site on 12 March. However, the EAE value on March 12 was high, indicating the dominance of fine-mode aerosols. In most cases, the presence of dominant coarse-mode aerosols was noted when the EAE value was <~1. Moreover, the transport of dust was infrequent, despite air masses regularly arriving from the Changchun region. Thus, interpreting the aerosol composition in Chiba using one parameter (i.e., AAE or EAE) might lead to incorrect conclusions. Thus, multiple variables (EAE, AAE, and aerosol volume size distribution) were used to interpret information on aerosol composition. Figure 6 shows the sky radiometer and GCOM-C AOT, AE values, and aerosol size distribution on days with different aerosol composition. Three coincident days were selected based on AAE and EAE values, as well as the aerosol volume size distribution. On April 5, dust was transported to the Chiba site, as supported by the AAE and EAE values and the volume size distribution. The AE values in both datasets showed good agreement despite the large uncertainties in GCOM-C AE values. On April 5, both sky radiometer and GCOM-C observations showed high AOT values, compared to the other two days. This indicated the GCOM-C observations were able to capture the presence of dust on April 5, despite the differences in the absolute AOT values. On April 4, the AOT values in both datasets didn't agree within the error ranges, despite a good agreement in the AE values. Figure 6c shows the differences in the AOT values as a function of differences in the AE values. The differences in the AOT values didn't show any clear dependence on the differences in the AE values. Thus, it is unlikely that the impact of changes in the aerosol composition would have a significant impact on the differences in the AOT values.

Impact of Aerosol Composition Change on the Observed Differences in the Datasets in Chiba
the presence of dominant coarse-mode aerosols was noted when the EAE value was < ~1. Moreover, the transport of dust was infrequent, despite air masses regularly arriving from the Changchun region. Thus, interpreting the aerosol composition in Chiba using one parameter (i.e., AAE or EAE) might lead to incorrect conclusions. Thus, multiple variables (EAE, AAE, and aerosol volume size distribution) were used to interpret information on aerosol composition. Figure 6 shows the sky radiometer and GCOM-C AOT, AE values, and aerosol size distribution on days with different aerosol composition. Three coincident days were selected based on AAE and EAE values, as well as the aerosol volume size distribution. On April 5, dust was transported to the Chiba site, as supported by the AAE and EAE values and the volume size distribution. The AE values in both datasets showed good agreement despite the large uncertainties in GCOM-C AE values. On April 5, both sky radiometer and GCOM-C observations showed high AOT values, compared to the other two days. This indicated the GCOM-C observations were able to capture the presence of dust on April 5, despite the differences in the absolute AOT values. On April 4, the AOT values in both datasets didn't agree within the error ranges, despite a good agreement in the AE values. Figure 6c shows the differences in the AOT values as a function of differences in the AE values. The differences in the AOT values didn't show any clear dependence on the differences in the AE values. Thus, it is unlikely that the impact of changes in the aerosol composition would have a significant impact on the differences in the AOT values.

Biomass Burning Influence at the Phimai Site
According to Figure 7a, larger differences in AOT values in Phimai occurred primarily from March to May. In Phimai, the periods from January to April and July to September are considered

Biomass Burning Influence at the Phimai Site
According to Figure 7a, larger differences in AOT values in Phimai occurred primarily from March to May. In Phimai, the periods from January to April and July to September are considered the dry and wet seasons, respectively [22]. During the dry season, the Phimai site is influenced by biomass burning, which has been described in detail in several studies [22,58,59]. The mean AAE value during Remote Sens. 2020, 12, x FOR PEER REVIEW 15 of 22 the dry and wet seasons, respectively [22]. During the dry season, the Phimai site is influenced by biomass burning, which has been described in detail in several studies [22,58,59]. The mean AAE value during The dry season in Phimai was 1.57, which lied within the range of dominant biomass burning aerosols [18]. Figure 7b shows the seasonal variations in AOT values at 380 nm and glyoxal (CHOCHO) levels at Phimai. CHOCHO, which is predominantly an oxidization product of biogenic volatile organic carbons, can also be emitted from biomass burning and biofuel use [60]. Strong agreement between CHOCHO and fire radiative power (FRP) has been reported in the region, where CHOCHO sources are predominantly pyrogenic [61]. Figure 7b also shows good agreement between the seasonality of CHOCHO VCD and sky radiometer AOT values at 380 nm. This result indicated that higher AOT values mostly occurred during the dry season due to the influence of biomass burning. Figure 7c also shows differences in AOT values at 380 nm as a function of CHOCHO VCD. CHOCHO VCD values for the day prior to the coincident date between MAX-DOAS, sky radiometer, and GCOM-C observations were used to construct Figure 7c. The influence of biomass burning at Phimai depends on the fire intensity (i.e., FRP), location of burning, and movement of air masses to the site. Based on CHOCHO data obtained around the time of GCOM-C overpass, the influence of biomass burning might not be well captured, as the time of fire occurrence was unknown. CHOCHO data with a lag of 1 day was more likely to capture the biomass burning effect. According to Figure  7c, some large differences in AOT values at 380 nm were associated with high CHOCHO VCD values. This means the large differences in the AOT values mostly occurred during the intense biomass burning period. An example using FRP and backward trajectories is shown in Figure 8. The dry season in Phimai was 1.57, which lied within the range of dominant biomass burning aerosols [18]. Figure 7b shows the seasonal variations in AOT values at 380 nm and glyoxal (CHOCHO) levels at Phimai. CHOCHO, which is predominantly an oxidization product of biogenic volatile organic carbons, can also be emitted from biomass burning and biofuel use [60]. Strong agreement between CHOCHO and fire radiative power (FRP) has been reported in the region, where CHOCHO sources are predominantly pyrogenic [61]. Figure 7b also shows good agreement between the seasonality of CHOCHO VCD and sky radiometer AOT values at 380 nm. This result indicated that higher AOT values mostly occurred during the dry season due to the influence of biomass burning. Figure 7c also shows differences in AOT values at 380 nm as a function of CHOCHO VCD. CHOCHO VCD values for the day prior to the coincident date between MAX-DOAS, sky radiometer, and GCOM-C observations were used to construct Figure 7c. The influence of biomass burning at Phimai depends on the fire intensity (i.e., FRP), location of burning, and movement of air masses to the site. Based on CHOCHO data obtained around the time of GCOM-C overpass, the influence of biomass burning might not be well captured, as the time of fire occurrence was unknown. CHOCHO data with a lag of 1 day was more likely to capture the biomass burning effect. According to Figure 7c, some large differences in AOT values at 380 nm were associated with high CHOCHO VCD values. This means the large differences in the AOT values mostly occurred during the intense biomass burning period. An example using FRP and backward trajectories is shown in Figure 8  As the air masses present on March 8 had barely passed over any area of fire activity, the impact of fires at the site was low, and the difference between datasets was −0. 04. This indicated that, in addition to clouds, biomass burning had a significant impact on the observed differences in the observed AOT values in Phimai.

Diurnal Variation of AOT at 380 nm Inferred from GCOM-C and OMI Observations
The difference in the GCOM-C (10:30 LT) and OMI (13:30 LT) overpass time could be used to study the diurnal variations in AOT at 380 nm. The OMI UV aerosol products were provided at 388 nm, which was similar to the wavelength (i.e., 380 nm) of the sky radiometer and GCOM-C UV aerosol product. Uncertainty related to the slight difference in wavelength is expected to have a nonsignificant impact in the analysis; thus, the OMI AOT at 388 nm was not converted to 380 nm. However, in the text and figure, these data are referred to as AOT at 380 nm to avoid confusion and maintain consistency. Sky radiometer observations within ±20 min of the OMI overpass time (13:30 LT) were selected. For both OMI and GCOM-C, points within a 20-km radius of the sites were used. OMI data with a cloud flag and cloud fraction of 0 and < 0.2, respectively, were used. Figure 9 shows the diurnal variations in AOT at 380 nm over both sites on some selected coincident dates between the sky radiometer and satellite (GCOM-C + OMI) observations. In Chiba, the sky radiometer AOT values were higher compared to the morning values. Such an increase in the AOT values was also seen in the satellite (GCOM-C + OMI) observations, despite the differences in the absolute AOT  As the air masses present on March 8 had barely passed over any area of fire activity, the impact of fires at the site was low, and the difference between datasets was −0. 04. This indicated that, in addition to clouds, biomass burning had a significant impact on the observed differences in the observed AOT values in Phimai.

Diurnal Variation of AOT at 380 nm Inferred from GCOM-C and OMI Observations
The difference in the GCOM-C (10:30 LT) and OMI (13:30 LT) overpass time could be used to study the diurnal variations in AOT at 380 nm. The OMI UV aerosol products were provided at 388 nm, which was similar to the wavelength (i.e., 380 nm) of the sky radiometer and GCOM-C UV aerosol product. Uncertainty related to the slight difference in wavelength is expected to have a non-significant impact in the analysis; thus, the OMI AOT at 388 nm was not converted to 380 nm. However, in the text and figure, these data are referred to as AOT at 380 nm to avoid confusion and maintain consistency. Sky radiometer observations within ±20 min of the OMI overpass time (13:30 LT) were selected. For both OMI and GCOM-C, points within a 20-km radius of the sites were used. OMI data with a cloud flag and cloud fraction of 0 and <0.2, respectively, were used. Figure 9 shows the diurnal variations in AOT at 380 nm over both sites on some selected coincident dates between the sky radiometer and satellite (GCOM-C + OMI) observations. In Chiba, the sky radiometer AOT values were higher compared to the morning values. Such an increase in the AOT values was also seen in the satellite (GCOM-C + OMI) observations, despite the differences in the absolute AOT values.
This indicated that quantitative information on the temporal variation of aerosols over heterogeneous surfaces could be obtained from the combination of GCOM-C and OMI datasets. diurnal variations. Moreover, large differences between satellite (GCOM-C + OMI) and sky radiometer AOT values were observed. The example dates for the Phimai site were during the dry season when the site is strongly influenced by biomass burning (discussed in the previous sections). The differences in the diurnal variations could be related to the impact of biomass burning on the observations. It would be interesting to see the agreement in the diurnal variations during the wet season when such conditions (biomass burning) don't prevail. However, no coincident dates during the wet season were available within the present datasets. This issue will be addressed in our future studies. Overall, the combination of GCOM-C and OMI observations provided a good opportunity to infer aerosol characteristics in the UV region.

Discussion
The Phimai site experiences biomass burning impact during the dry season. The biomass burning influence had a significant impact on the large differences in the AOT values at Phimai, in addition to cloud contamination. This was also reflected in the high MBE values in Phimai compared to Chiba for the same coincident criterion of the satellite observations. It is expected that GCOM-C and sky radiometer observation in Phimai will show better agreement during the wet season due to the absence of strong emission sources. Such a hypothesis can be tested with long term datasets for Phimai and other rural sites of SKYNET. Some potential reasons for the large differences observed in the datasets under biomass burning influence are -(1) The distance of the ground-based stations from the fire locations is shorter compared to the satellite sensors. Thus, the strong influence of the fires on the ground-station can infer large differences in the comparison results. (2) The time of occurrence and the location of the fires are unknown and random. The optical and physical properties (size distribution, chemical composition, hygroscopicity, etc.) of biomass burning particles change Over Phimai, the sky radiometer and satellite (GCOM-C + OMI) observations showed different diurnal variations. Moreover, large differences between satellite (GCOM-C + OMI) and sky radiometer AOT values were observed. The example dates for the Phimai site were during the dry season when the site is strongly influenced by biomass burning (discussed in the previous sections). The differences in the diurnal variations could be related to the impact of biomass burning on the observations. It would be interesting to see the agreement in the diurnal variations during the wet season when such conditions (biomass burning) don't prevail. However, no coincident dates during the wet season were available within the present datasets. This issue will be addressed in our future studies. Overall, the combination of GCOM-C and OMI observations provided a good opportunity to infer aerosol characteristics in the UV region.

Discussion
The Phimai site experiences biomass burning impact during the dry season. The biomass burning influence had a significant impact on the large differences in the AOT values at Phimai, in addition to cloud contamination. This was also reflected in the high MBE values in Phimai compared to Chiba for the same coincident criterion of the satellite observations. It is expected that GCOM-C and sky radiometer observation in Phimai will show better agreement during the wet season due to the absence of strong emission sources. Such a hypothesis can be tested with long term datasets for Phimai and other rural sites of SKYNET. Some potential reasons for the large differences observed in the datasets under biomass burning influence are -(1) The distance of the ground-based stations from the fire locations is shorter compared to the satellite sensors. Thus, the strong influence of the fires on the ground-station can infer large differences in the comparison results. (2) The time of occurrence and the location of the fires are unknown and random. The optical and physical properties (size distribution, chemical composition, hygroscopicity, etc.) of biomass burning particles change considerably during the first 2-4 h of their atmospheric transport [62]. Such changes might not reflect in the satellite measurements due to the limited observation period (one measurement a day during the overpass time).
Comparing the OMI and collocated AERONET AOT at 380 nm at more than twenty-two sites, Ahn et al. [9] reported that R, RMSE, and slope values were in the range of 0.41-0.91, 0.08-0.25, and 0.33-0.92, respectively. The statistics of the comparison between SKYNET and GCOM-C are within the range of the reported values, except for the RMSE value (~0.35) at the Phimai site, where biomass burning influences can incur large differences in the AOT values. However, the study of Ahn et al. [9] reported RMSE value around 0.18 for sites located in the biomass burning regions. The limited number of collocated data points can be one of the potential reasons for the higher RMSE value at the Phimai site. The comparisons of Ahn et al. [9] included at least~100 collocated data points over four years. Moreover, despite the high RMSE value at the Phimai site, the R and slope values were almost similar to or better in some cases than the biomass burning sites included in the study of Ahn et al. [9]. Furthermore, in the diurnal variations in AOT at Phimai (Figure 9), the difference in the GCOM-C and sky radiometer was lower compared to that of between the sky radiometer and OMI. Thus, utilizing a long-term dataset, the statistics of the comparison between GCOM-C and sky radiometer is expected to improve, and the impact of the improved spatial resolution of GCOM-C will be more evident.
In addition to the potential reasons discussed here, uncertainties in the sky radiometer inversion product may also contribute to the observed biases. For example, the solar disk scan method can underestimate the calibration constant, resulting in an overestimation of sky radiances and SSA [63]. Moreover, the different assumptions in both retrieval algorithms can also impact the differences in the AOT values.
Apart from the discussions in this work, qualitative assessment of the retrieval algorithms (i.e., different assumptions, error in the model calculations, etc.) is required to identify further potential reasons for the disagreement in the AOT values. A qualitative assessment of the GCOM-C AOT data utilizing long term datasets and including more representative ground sites will be addressed in our future study.

Conclusions
The newly available GCOM-C AOT at 380 nm was evaluated using sky radiometer observations at the SKYNET Chiba and Phimai sites. GCOM-C is currently the only satellite providing AOT at 380 nm during morning overpass time. The sky radiometer observations in Chiba were compared with coincident AERONET observations, and the agreement between the datasets was mostly within ±0.02. The sky radiometer and MAX-DOAS observations in Phimai also showed good agreement mostly under clear sky conditions. At both the sites, the agreement between the sky radiometer and GCOM-C AOT at 380 nm was mostly within ±0.2, with a positive correlation of~0.73. The number of coincident points for the Chiba and Phimai site was 47 and 31, respectively, for a coincidence criterion of <=30 km. Larger differences in AOT values occurred at Chiba, mostly due to the effect of clouds on the observations. At times, long-range dust transport occurred at Chiba. Such a change in the aerosol composition did not show any clear impact on the differences in the AOT. The Phimai site is affected by biomass burning during the dry season. The influence of biomass burning had a significant impact on the differences observed in the datasets at the Phimai site. The difference in the overpass time of GCOM-C and OMI was utilized to study the diurnal variations in the AOT at 380 nm. The combination of GCOM-C and OMI AOT captured the diurnal variations in AOT at Chiba only. At Phimai, the coincident days among GCOM-C, OMI, and sky radiometer were only available for the dry season. Thus, the strong influence of biomass burning potentially led to the poor agreement of the diurnal variations in Phimai. Overall, the AOT values retrieved in both datasets (sky radiometer and GCOM-C) were consistent, indicating the strong potential of the GCOM-C UV AOT product.