Next Article in Journal
A Self-Calibrated Non-Parametric Time Series Analysis Approach for Assessing Insect Defoliation of Broad-Leaved Deciduous Nothofagus pumilio Forests
Previous Article in Journal
Quantifying the Evapotranspiration Rate and Its Cooling Effects of Urban Hedges Based on Three-Temperature Model and Infrared Remote Sensing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Validation of OMI HCHO Products Using MAX-DOAS observations from 2010 to 2016 in Xianghe, Beijing: Investigation of the Effects of Aerosols on Satellite Products

1
State Key Laboratory of Remote Sensing Science, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing 100101, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
China Centre for Resources Satellite Data and Application, Beijing 100830, China
*
Authors to whom correspondence should be addressed.
Remote Sens. 2019, 11(2), 203; https://doi.org/10.3390/rs11020203
Submission received: 6 December 2018 / Revised: 15 January 2019 / Accepted: 15 January 2019 / Published: 21 January 2019
(This article belongs to the Section Atmospheric Remote Sensing)

Abstract

:
Formaldehyde (HCHO) is one of the most abundant hydrocarbons in the atmosphere. Its absorption features in the 320–360 nm range allow its concentration in the atmosphere to be retrieved from space. There are two versions of HCHO datasets derived from the Ozone Monitoring Instrument (OMI)—one provided by the Royal Belgian Institute for Space Aeronomy (BIRA-IASB) and one provided by the National Aeronautics and Space Administration (NASA)—referred to as OMI-BIRA and OMI-NASA, respectively. We conducted daily comparisons of OMI-BIRA and multi-axis differential optical absorption spectrometry (MAX-DOAS), OMI-NASA and MAX-DOAS, and OMI-BIRA and OMI-NASA and monthly comparisons of OMI-BIRA and MAX-DOAS and OMI-NASA and MAX-DOAS. Daily comparisons showed a strong impact of effective cloud fraction (eCF), and correlations were better for eCF < 0.1 than for eCF < 0.3. By contrast, the monthly and multi-year monthly mean values yielded correlations of R2 = 0.60 and R2 = 0.95, respectively, for OMI-BIRA and MAX-DOAS, and R2 = 0.45 and R2 = 0.78 for OMI-NASA and MAX-DOAS, respectively. Therefore, use of the monthly mean HCHO datasets is strongly recommended. We conducted a sensitivity test for HCHO air mass factor (AMF) calculations with respect to the HCHO profile, the aerosol extinction coefficient (AEC), the HCHO profile–AEC combination, the aerosol optical depth (AOD), and the single scattering albedo (SSA) to explicitly account for the aerosol optical effects on the HCHO AMF. We found that the combination of AEC and HCHO profiles can account for 23–39% of the HCHO AMF variation. Furthermore, a high load of absorptive aerosols can exert a considerable effect (−53%) on the AMF. Finally, we used the HCHO monthly mean profiles from Goddard Earth Observing System coupled to Chemistry (GEOS-Chem), seasonal mean AECs from Cloud-Aerosol LIDAR with Orthogonal Polarization (CALIOP) and monthly climatologies of AOD and SSA from the OMAERUV (OMI level-2 near UV aerosol data product) dataset at Xianghe station to determine the aerosol correction. The results reveal that aerosols can account for +6.37% to +20.7% of the HCHO monthly change. However, the changes are greatest in winter and are weaker in summer and autumn, indicating that the aerosol correction is more applicable under high-AAOD conditions and that there may be other reasons for the significant underestimation between satellite and MAX-DOAS observations.

1. Introduction

In recent years, with air pollution control, PM2.5 has declined significantly in China (http://kjs.mep.gov.cn/). O3 is now the only pollutant whose mass burden has continued to increase in the 74 experimental cities of China, and the mass concentration is expected to continue to increase [1,2]. Furthermore, surface ozone has become the primary pollutant in many cities in China in the summer [3]. Volatile organic compounds (VOCs), as an important precursor of ozone, play an important role in air quality [4,5]. Studies have shown that VOCs are the key to ozone control [2,6]. In addition, VOCs can also generate secondary organic aerosols [7]. Formaldehyde (HCHO) is an intermediate of the oxidation reaction of various VOCs in the atmosphere. Therefore, HCHO can be used as a tracer for VOCs in the absence of other VOC observations [8]. Background HCHO in the atmosphere is mainly derived from the oxidation of CH4. In continental regions, biomass burning, direct emissions from the industrial sector, and the oxidation of non-methane VOCs (NMVOCs) emitted by biological sources can result in regional features of the spatial and temporal distributions of HCHO [9,10]. Due to photochemical reactions and the oxidation of OH radicals, the life cycle of HCHO is only 1.5 h [11].
HCHO has a unique and sufficiently strong absorption signature in the 320–360 nm range that can be used to detect HCHO concentrations from space. UV-Vis hyperspectral sensors, such as Global Ozone Monitoring Experiment (GOME), prove that it is possible to monitor the HCHO column [12]. HCHO retrieval has also been conducted based on GOME-2 [13], Scanning Imaging Absorption spectrometer for Atmospheric ChartographY (SCIAMACHY) [10], the Ozone Measurement Instrument (OMI) [14,15], Ozone Mapping and Profiler Suite (OMPS) [16,17], and TROPOspheric Monitoring Instrument (TROPOMI) [18]. The OMI has been in orbit for 15 years, and its long-term observation of HCHO provides valuable HCHO information. There are two HCHO algorithms based on OMI data: the operational OMHCHO product published by the National Aeronautics and Space Administration (NASA) and the HCHO dataset released by the Belgian Institute for Space Aeronomy (BIRA). For convenience, these products are referred to as OMI-NASA and OMI-BIRA, respectively. With improvements in sensors, many users have already benefitted from these products for several applications, such as the study of long-term trends [10,19], top-down VOC inventory studies [20,21], sensitivity analysis of ozone production [8], and biomass burning source emission assessment [22].
According to previous studies, eastern China and the Pearl River Delta region in southern China are hotspots of HCHO. However, previous validation results in Xianghe show that the regression line between satellite observation and multi-axis differential optical absorption spectrometry (MAX-DOAS) strongly depends on the parameter settings when calculating the air mass factor (AMF). Clear underestimation is observed when HCHO profiles from IMAGES mode are used a priori (−2.6 × 1015 molec/cm2); by contrast, the correction of the multi-year average data during the period 2008–2013 increases from 0.6 to 0.9 when the HCHO profiles from local MAX-DOAS are used, indicating an underestimation of −0.9 × 1015 molec/cm2 [15]. The validation results by Wang et al. [23] in Wuxi show that HCHO is biased by approximately −20%, 20%, and 12% based on OMI, GOME-2A and GOME-2B, respectively. Under clear-sky conditions using the vertical profile of MAX-DOAS as the shape factor, the AMF of HCHO changes by 11%, and the correlation is improved by 35%. However, the results show that the HCHO profile cannot fully explain the underestimation phenomenon; other error sources should be identified to further improve the satellite retrieval algorithms.
The AMF calculation error is one of the main sources of uncertainty in the HCHO satellite retrieval products [18]. The AMF uncertainty is mainly due to assumptions about atmospheric conditions, especially for the HCHO shape factor and cloud-aerosol correction [24]. In existing HCHO inversion algorithms, the aerosol information is implicitly considered in the cloud information, whereas there is no explicit consideration of the aerosol contribution. However, the contribution of aerosols to trace gases (NO2) has received increasing attention. Lin et al. [25,26] quantitatively assessed the effects of aerosols on cloud parameters and the NO2 AMF during NO2 retrieval in China and developed the POMINO algorithm, which can calculate the NO2 AMF online without look-up tables and considers aerosols explicitly. The results show that the correlation coefficient (R2) of the POMINO algorithm increased from 0.72 to 0.96, indicating that an explicit treatment of aerosols is important for space-based NO2 retrievals. However, because the aerosol parameters (extinction coefficient profile, aerosol optical depth (AOD), single scattering albedo (SSA), asymmetry factor, etc.) are complex, it is difficult to quantify them using the conventional look-up table method and, thus, the effects of aerosols on the HCHO AMF have not been explicitly considered. Both HCHO retrieval algorithms implicitly exploit aerosol information in cloud information. To some extent, the scattering effect of non-absorptive aerosols implicit in cloud information can offset the contribution of some aerosols to trace gas retrieval [18,27] because non-absorptive aerosols and clouds have similar radiation contributions. However, absorptive aerosols have the opposite effect of clouds, which can reduce the spectral reflectance. Compared with considering the absorptive aerosol contribution (AMFa), not considering the absorptive aerosol contribution (AMFoff) will cause an overestimation of the AMF, which in turn will underestimate the concentration of HCHO in the vertical column.
In previous studies, the satellite-retrieved HCHO datasets produced by different research institutions were validated using the MAX-DOAS or model simulation individually. No comparative evaluation of the applicability of different algorithms in specific regions has been performed. In this study, we validate daily and monthly averages of the HCHO data derived from OMI by BIRA and NASA using the MAX-DOAS observations in Xianghe. To account for the impact of clouds on both MAX-DOAS and satellite retrieval, we adopt the method used by Wang et al. [23], which separately evaluates the cloud effects on both satellite and MAX-DOAS observations. Additionally, the aerosol optical proprieties at Xianghe station derived from OMI are investigated, and the aerosol effect on the HCHO AMF is examined using SCIATRAN. Finally, the monthly mean HCHO vertical column density (VCD) is corrected based on the contribution of aerosols.
The paper is organized as follows: in Section 2, we describe the MAX-DOAS observations in Xianghe and the satellite products involved in this study. In Section 3, we compare the daily and monthly HCHO VCDs derived from BIRA and NASA with MAX-DOAS observations. In addition, we evaluate the effect of aerosols on the HCHO AMF and correct its effect on the monthly mean HCHOBIRA and HCHONASA. Section 4 discusses the validation results. Section 5 presents the summary and conclusions.

2. MAX-DOAS Measurements and Satellite Datasets

2.1. MAX-DOAS Instrument and Data Analysis

The MAX-DOAS measurement technique has been developed to retrieve tropospheric trace gas total columns and profiles. The most recent generation of MAX-DOAS instruments allows the measurement of aerosols and a number of tropospheric pollutants, such as NO2, HCHO, SO2, HONO, O4 and CHOCHO [28,29,30]. The MAX-DOAS instrument used in this study was constructed and assembled by the Royal Belgian Institute for Space Aeronomy (BIRA-IASB) (see Clémer et al. [31]). It consists essentially of a telescope mounted on a sun-tracker (which can point at any angle and in any azimuthal direction) combined with two spectrographs: one for the UV spectrum (300–390 nm) and one for the visible spectrum (400–720 nm). HCHO retrieval is based on sequential observations at elevation angles of 2, 4, 6, 8, 10, 12, 15, 30 and 90°. During the period analysed in this study (2010–2016), the instrument was stationed approximately 55 km to the east-southeast of Beijing at the meteorological observatory in Xianghe (39.75°N, 116.96°E). The vertical profile retrievals can also be treated as optimal estimates because of the scientific-grade instruments installed in China [15]. Ground observations indicate a clear seasonal cycle in HCHO concentrations, with a winter minimum of approximately 5 × 1015 molec cm−2 and a summer maximum roughly five times higher [29]. In general, we expect that the VCDs from MAX-DOAS have much lower uncertainties than those from satellite observations, though the bias of MAX-DOAS HCHO VCDs is −23 ± 28% with respect to near-surface volume mixing ratios [23,29].

2.2. HCHO Versions from OMI

2.2.1. OMI-BIRA HCHO Product

The BIRA-IASB algorithm for the retrieval of HCHO columns was updated by De Smedt et al. [15]. To better account for O2-O2 and BrO absorption, the new version (V14) of the HCHO retrieval scheme includes three fitting intervals: 339–364 nm for the pre-fit of O4, 328.5–359 nm for the pre-fit of BrO and 328.5–346 nm for the fit of the HCHO slant column density (SCD). QDOAS software is used for the DOAS retrieval, with the earthshine radiances averaged in the equatorial Pacific as reference spectra. A destriping correction and background normalization are applied in the across-swath position. The 340-nm HCHO AMF is calculated using the LIDORT v3.3 radiative transfer model, which uses a priori HCHO profiles from a global 2.0° (latitude) × 2.5° (longitude) IMAGESv2 chemistry transport model simulation. Cloud effects are corrected using the independent pixel approximation [15]. No explicit correction is applied for aerosols, but the cloud correction scheme accounts for a large part of their scattering. All rows affected by an OMI row anomaly are identified using the fitting residuals as follows: fitting residuals 3 times larger than the average fitting residual, pixels with effective cloud fractions of >0.4, solar zenith angles larger than 70°, and individual vertical column errors larger than 3 times the column are assigned a poor-quality (<0) flag.

2.2.2. OMI-NASA HCHO Product

The official NASA HCHO product is provided by the updated Smithsonian Astrophysical Observatory (SAO) retrieval, as described in González Abad et al. [14] and evaluated in Zhu et al. [32]. HCHO SCDs are retrieved through a direct non-linear least-squares fitting of spectral radiances within the interval 328.5–356.5 nm. The retrieval algorithm includes dynamic calibration of solar and radiance wavelengths, an under-sampling correction, and a common mode residual spectrum. Daily radiances are used as reference spectra. The retrieved SCDs are converted to VCDs using AMFs taken from the look-up table pre-computed using the VLIDORT radiative transfer model, which uses a priori HCHO profiles from a global 2.0° (latitude) × 2.5° (longitude) GEOS-Chem chemistry transport model simulation. A daily post-processing normalization correction (as a function of latitude and detector row) is applied to reduce retrieval biases, minimize noise, and reduce cross-track striping. Only good pixels with a main data quality flag of 0 are used in this work. Barkley et al. [33] used the OMI-NASA HCHO data as a factor in the assessment of air quality changes over the Middle East during 2005–2014.

2.3. Statistical Evaluation Methods

In this study, the OMI-BIRA and OMI-NASA HCHO products are compared with MAX-DOAS observations and with each other. The root mean square error (RMSE) between the satellite products and MAX-DOAS columns is used to quantify the systematic biases and is calculated following Equation (1).
R M S E = 1 n ( V S i V M i ) 2 n

2.4. AMF Calculation

The AMF is defined as the ratio of the SCD to the VCD. A geometric AMF ( A M F G ) can be defined as A M F G = sec θ s + sec θ v , which is a function of the solar zenith angle θ s and the satellite viewing angle θ v in the absence of atmospheric scattering. In the actual atmosphere, the AMF also depends on the trace gas profile, air pressure, surface albedo, and aerosol profiles. In this work, the AMFs are calculated at 340 nm by SCIATRAN [34,35]. The HCHO profiles from GEOS-Chem, which feature maximum values concentrated near the ground and an exponential decrease in values with height, are used to specify the vertical distribution. High-resolution GEOS-Chem model runs were performed for the period 2007 (0.667° long. × 0.5° lat.) and 2014 (0.3125° long. × 0.25° lat.), and each run involved 36 layers, with ~10 layers below 3 km. The model profiles from this run were also used for NO2 retrievals over China [25,26]. The HCHO AMF is sensitive to the vertical distribution of atmospheric species, and a more sophisticated calculation of the HCHO AMF includes the effects of clouds [15,18] and aerosols [25,26]. The AMF is calculated by applying the following equation:
A M F i = 1 V O D i ln I i I + i ,
where I + i and I i are the intensities with and without the absorber i , respectively. The vertical optical density V O D i ( λ ) from the ground to the top of the atmosphere h 0 is defined as V O D i ( λ ) = 0 h 0 d h ρ i ( h ) σ i ( λ , h ) , where h is the altitude, ρ is the number density of the absorber, and σ is the absorption cross-section at a given wavelength. This method is similar to [36] and is valid only for an optically thin atmosphere.

3. Validation of Satellite Datasets

Figure 1 shows the spatial distribution of multi-year mean HCHO data from OMI-BIRA and OMI-NASA during the period 2010–2016. The daily averaged OMI-BIRA (in square dots) and OMI-NASA (in circle dots) data within 50 km around Xianghe station and the daily averaged MAX-DOAS data (in triangle dots) within the period 12:00–15:00 local time (LT) (corresponding to the satellite overpass time of 13:30) are shown in Figure 2. The black, red and green lines represent the monthly averaged OMI-BIRA, MAX-DOAS and OMI-NASA data, respectively. The vertical bars indicate one standard deviation. The comparison in the following section is based on these data.

3.1. Daily Comparisons

Based on the daily averaged data in Figure 2, comparisons were performed between the daily products of OMI-BIRA and MAX-DOAS, OMI-NASA and MAX-DOAS, and OMI-NASA and OMI-BIRA. To characterize the cloud effect on the comparisons, the effective cloud fraction (eCF), which is obtained from the OMCLDO2 cloud product [37], was classified using value ranges of 0–0.3, 0–0.2, and 0–0.1 for the three comparisons.
The coincident observation number of correlative days for eCF < 0.3, eCF < 0.2 and eCF < 0.1 amounts to 1003 (40.17%), 748 (29.96%) and 447 (17.90%), respectively, in Xianghe during the period 2010.03–2016.12. Notably, a previous study found that cloud effects on MAX-DOAS results are negligible for satellite validation activities and that cloud effects on the HCHO products become substantial for eCF > 0.3 [23]. Figure 3a–c displays scatter plots of the daily averaged OMI-BIRA and OMI-NASA datasets versus the corresponding MAX-DOAS measurements for different eCF bins of 0–0.3, 0–0.2 and 0–0.1, respectively. The correlation coefficient (R) for eCF < 10% is better than the R for eCF < 30% for the three groups of comparisons. The OMI-BIRA and OMI-NASA datasets have similar correlations with MAX-DOAS, i.e., R = 0.26 (R = 0.34) and R = 0.27 (R = 0.32), respectively, for eCF < 0.3 (eCF < 0.1). The correlations between OMI-NASA and OMI-BIRA are 0.48, 0.53 and 0.56 for eCF < 0.3, eCF < 0.2, and eCF < 0.1, respectively, based on the daily comparison.
Figure A1 to Figure A3 show daily comparisons between OMI-BIRA and MAX-DOAS, OMI-NASA and MAX-DOAS, and OMI-NASA and OMI-BIRA by season. Figure A1 shows that the correlation between OMI-BIRA and MAX-DOAS is the best in spring, with R = 0.25 (R = 0.26) for 0 < eCF < 0.2 (0 < eCF < 0.1). The changes in eCF can change the correlation between the satellite data and the MAX-DOAS observations when the upper limit of eCF changes from 0.3 to 0.2 for all seasons except winter; however, less variation is observed for stricter cloud screening (eCF < 0.1) for summer. In winter, there is a weak correlation or even a negative correlation between OMI-BIRA and MAX-DOAS, indicating that the winter satellite data need to be used cautiously. Figure A2 shows that OMI-NASA has a better correlation with MAX-DOAS in spring. When stricter eCF control is employed, the satellite dataset’s correlation with the ground observations can be improved for spring and summer, but the two datasets have a poor correlation in summer and winter. OMI-NASA and OMI-BIRA have similar correlations except for winter, with R values ranging from 0.25 to 0.48 for eCF < 0.3 (Figure A3). The differences in dependence on eCF of OMI-BIRA and MAX-DOAS, OMI-NASA and MAX-DOAS, and OMI-BIRA and OMI-NASA will be discussed in Section 4.

3.2. Monthly Comparisons

This section describes the monthly comparisons of OMI-BIRA and MAX-DOAS, and OMI-NASA and MAX-DOAS. The R2 values of OMI-BIRA and MAX-DOAS, and OMI-NASA and MAX-DOAS are 0.60 and 0.45, respectively, and the corresponding RMSEs are 2.20 × 1015 molec/cm2 and 2.90 × 1015 molec/cm2. These results are significantly better than the daily correlations in Section 3.1. The multi-year monthly average is shown in Figure 4b. The R2 of OMI-BIRA (OMI-NASA) and MAX-DOAS is 0.95 (0.78), with an RMSE of 0.72 × 1015 molec/cm2 (1.64 × 1015 molec/cm2). The correlation is significant within the 95% confidence interval. The multi-year monthly averages of the three datasets are shown in Figure 5. The multi-year monthly average of OMI-BIRA is significantly lower than that of MAX-DOAS. The largest difference occurs in June, reaching −9.472 × 1015 molec/cm2 (−59.4723%), and the smallest difference occurs in January, with a value of −1.337 × 1015 molec/cm2 (−15.05%), as indicated in Table 1. By contrast, the multi-monthly average of OMI-NASA is higher than that of MAX-DOAS from October to February, with a range of +9.70% to 33.74%, and is slightly lower in March–October, with the smallest difference occurring in September. October has the largest negative difference, reaching −16.28%. However, the difference is less significant than that of OMI-BIRA.

3.3. Aerosol Effects on the Satellite Results

To further understand the impact of aerosols on the HCHO AMF calculation, we performed a sensitivity analysis to determine the dependence of the HCHO AMF calculation on the vertical distribution of the HCHO profile, vertical distribution of the aerosol extinction coefficient (AEC), AOD, and SSA.

3.3.1. HCHO AMF Dependence on the HCHO Profile

Cluster analysis of the HCHO vertical profile from GEOS-Chem in 2007 over China was performed, and the results of the 5–8 categories (NHCHO) are shown in Figure 6a. The AMF corresponding to the vertical HCHO distribution is shown in Figure 6b. The results show that the position and magnitude of the HCHO maximum between different clusters are well-defined and can represent background and hotspot regions (e.g., cluster5 (blue line) and cluster6 (brick-red line) for NHCHO = 8 represent low and high HCHO concentrations near the surface, respectively], in addition to elevation information (e.g., cluster3 (magenta line) for NHCHO = 6, cluster7 (orange line) for NHCHO = 7 and cluster4 (green line) for NHCHO = 8 indicate HCHO profiles of high elevation). NHCHO = 8 can fully reflect the influence of the vertical HCHO distribution on the AMF calculation. The maximum difference in the AMF between different clusters is approximately 20%. The test results for NHCHO > 8 show that the effect on the AMF reflected by the different classes is basically equivalent to NHCHO = 8.

3.3.2. HCHO AMF Dependence on the AEC

In this study, the AEC provided by Regional Atmospheric Modeling System and Models-3 Community Multiscale Air Quality (RAMS-CMAQ) and Cloud-Aerosol LIDAR with Orthogonal Polarization (CALIOP) are used to test the effects of AEC on the calculation of the HCHO AMF. The air quality modelling system RAMS-CMAQ has been used to investigate regional atmospheric pollution and environmental issues in many studies [2,38,39]. The regional model CMAQs can simulate various chemical and physical processes valuable for understanding atmospheric trace gas transformation and distribution [40]. CALIOP, which is aboard the Calipso satellite, is an elastic backscatter LIDAR that transmits polarized laser light at 532 and 1064 nm [41]. The LIDAR signal inversion begins at approximately 30 km above the ground and continues to the surface. Its orbit repeats every 16 days [42]. Wu et al. [43] noted that RAMS-CMAQ can reflect the vertical distribution of aerosols in China. The AEC peak is ~0.5 km above the ground, but it is underestimated (0.1 km−1) in north-western China and overestimated (0.3 km−1) in central and eastern China. Considering the low temporal resolution of Calipso, the observations of the RAMS-CMAQ and CALIOP datasets complement each other to some extent. The two datasets over China in 2013–2015 were used to perform cluster analysis, and Figure 7a,b shows the results of 5-8 categories (NAEC) for RAMS-CMAQ and CALIOP, respectively. Figure 7 shows that the magnitude and height of the AEC peak between different clusters are different and that the differences in the AEC among the classes are significantly larger for CALIOP than for RAMS-CMAQ. The RAMS-CMAQ classes show a similar vertical shape factor, which consistently decreases with increasing altitude. By contrast, CALIOP has a finer vertical resolution than RAMS-CMAQ and can reflect more complex features, e.g., the peak position changes from the surface to 2 km. The dependence of the AMF on the two AEC datasets is tested with an AOD set to 1 and an SSA set to 0.92. We find that the changes in the AMF caused by changes in the AEC types are relatively insignificant. As shown in Figure 8, for RAMS-CMAQ, the change caused by the AEC profile type is small (between 5% and 8%), and for Calipso, the change is larger (between 12% and 23%). When more categories are used, more subtle differences can be reflected.

3.3.3. HCHO AMF Dependence on the AEC and HCHO Profile Combination

When the HCHO vertical profiles of NHCHO = 5 in Section 3.3.1 and the CALIOP AEC profiles of NAEC = 8 in Section 3.3.2 are selected, the range of AMF variation in the group is 23–39% (Figure 9), substantially larger than the variability induced by changes in the vertical HCHO profile and AEC individually. Therefore, the relative positions of the HCHO and AEC profiles affect the AMF results. Consequently, in the AMF calculation, it is necessary to consider the combination of both.

3.3.4. HCHO AMF Dependence on AOD and SSA

As shown in Figure 10, when the aerosol SSA remains unchanged, the HCHO AMF decreases with increasing AOD. Furthermore, the smaller the SSA, the greater the effect of a change in AOD on the AMF. For SSA = 0.82, when AOD changes from 0.1 to 2, the AMF decreases by 53%, and the AMF differences lessen with increasing AOD. When the SSA is 0.97, the AOD changes from 0.1 to 2, and the AMF decreases by 14%. When AOD remains unchanged, the AMF shows an increasing trend with increasing SSA. The larger the AOD, the more obvious the increasing trend of the AMF. Figure 10b shows the percentage change in the AMF with SSA = 0.82 as a reference. The change in the AMF is largest for SSA = 0.97 and AOD = 2, up to 83%.

3.4. Aerosol Correction

According to the results of the HCHO AMF sensitivity analysis, we extract monthly average HCHO profiles based on GEOS-Chem for 2007 and 2014 in the Xianghe area and the seasonal AEC averages during the period 2013–2015 as the initial input of the AMF calculation (Figure 11). At Xianghe station, the vertical distribution of HCHO decreases with altitude. The HCHO is mainly distributed in the height range of 2 km near the ground. The maximum is 2.8 ppb in August, and the minimum is less than 1 ppb in April. However, in December, the HCHO concentration near the ground is larger than the average of the adjacent months and is greater than 2 ppb. The vertical profile distribution of aerosols is similar in summer, autumn and winter. The maximum value occurs at 500 m. Except for the relatively high value at 500 m, the vertical distribution characteristics in spring differ from those in the other seasons in the range of 2–4 km. These differences may be related to dust transport in spring.
The AOD and SSA parameters are derived from the OMI level-2 near UV aerosol data product ‘OMAERUV’. OMAERUV products are derived from the measured reflectance data from OMI observations and climatological surface albedo data from Total Ozone Mapping Spectrometer (TOMS) at 354 nm and 388 nm [44]. The validation results at Xianghe station in 2005–2008 indicate that the correlation of AOD is 0.87, with an RMSE of 0.17 [45]. Furthermore, the validation results in 2005–2014 show that the correlation of AOD is 0.79, with an RMSE of 0.21 [46]. Despite the 13 × 24 km2 coarse spatial resolution, using near-UV observations for aerosol characterization has the advantage of reducing retrieval errors associated with land surface reflectance characteristics due to the low surface albedo in this spectral region. In this paper, monthly climatologies of aerosol AOD and SSA at Xianghe station are calculated. The regional heterogeneity of aerosol optical properties in China is large and varies with the seasons [47]. The monthly SSA ranges from 0.86 to 0.93. Starting in May, there was a marked rise in SSA, indicating a decrease in aerosol absorption. The seasonal cycle of AOD at 354 nm is the same as that in other regions: AOD increased significantly to a maximum of 1.30 in August and reached a minimum of 0.15 in October.
Using the monthly HCHO vertical profile as the reference scene, the seasonal vertical extinction coefficient profile at Xianghe station and the monthly AOD and SSA parameters are included in the AMF calculation. Compared with the reference scene, which does not consider aerosol parameters, the relative changes in November, December, April and March were large, i.e., 20.70%, 20.60%, 18.10%, and 10.31%, respectively. These changes were all positive, and the HCHO concentration increased after correction. The OMI-BIRA and OMI-NASA HCHO multi-month monthly mean data after aerosol correction are shown Figure 12. After the correction, the OMI-NASA results are higher than the ground-based MAX-DOAS results from October to May, and the monthly averages from May to September are closer to the MAX-DOAS data. OMI-BIRA increased in different proportions each month, and the correction reduced the discrepancy between the OMI-BIRA and MAX-DOAS data.

4. Discussion

In this work, the weak correlation values between the daily satellite datasets and the MAX-DOAS observations compared with NO2 [25,26,48] indicate that the uncertainties in the satellite daily HCHO columns are large, resulting in errors that cause these differences, especially in winter. The daily comparison results are in qualitative agreement with the validation results for Wuxi, where the OMI-BIRA and MAX-DOAS relationship has R2 = 0.12 for eCF < 0.1 [23]. The daily comparison shows a larger random error, and the cloud effect can partly explain the bias in the whole dataset. However, at the season scale, the random error remains significant, and the effect of cloud correction depends on the season. This dependence may be due to cloud properties in different seasons or the cloud retrieval algorithm. At low eCF (eCF < 0.1, close to the detection limit of the cloud fraction [37,49]), cloud top pressure retrieval is highly unstable, which can introduce uncertainly in HCHO retrieval [18]. This can also be seen in Figure A4, in which the correlation of monthly comparison for eCF < 0.3 is better than that for eCF < 0.1. In addition, because OMI-BIRA and OMI-NASA adopt similar cloud correction methods, the dependence of the correlation between the two on eCF is weaker than those of the correlations of the satellite datasets with MAX-DOAS measurements. The poor correction in winter can be attributed to the low winter seasonal character of HCHO over Xianghe and the detection limitations for OMI-BIRA (7 × 1015 molec/cm2) and OMI-NASA (10 × 1015 molec/cm2). In summary, for both OMI-BIRA and OMI-NASA, the daily data should be used with caution. Use of the daily average data in winter is not recommended.
The monthly comparison shows that the monthly average can reduce the random error in the satellite products and ground-based data. OMI-BIRA has a larger systematic underestimation error than OMI-NASA, consistent with previous studies [32]. The GEOS-Chem and IMAGES chemical transport models supply the HCHO profiles (shape factors) for OMI-NASA and OMI-BIRA, respectively. The HCHO profiles assumed for the HCHO AMF calculation may be a key factor accounting for the systematic bias [15,23]. The other main factor related to the systematic difference between OMI-BIRA and OMI-NASA is background correction (see Table 2 for more detail). The inconsistent simulation of the HCHO concentration in the chemical transport models in the reference sector will dominate the monthly systematic bias. Furthermore, compared with OMI-BIRA, no pre-processing related to BrO interference is conducted for OMI-NASA. The overestimate provided by OMI-NASA in winter months may be unrealistic and may be due to strong interference of BrO at high latitudes. Theys et al. [50] showed that the zonal mean of tropospheric BrO vertical columns was larger from September to April of the next year at 30–60°N. The monthly comparisons prove that the monthly mean value of satellite observations can be used for air quality assessment. In addition, the RMSE of OMI-BIRA is smaller than that of OMI-NASA. This difference can be attributed to the consistency of the algorithm between the ground-based MAX-DOAS and the satellite-based OMI-BIRA.
With respect to the aerosol effect on HCHO AMF, our results indicate that the combination of the HCHO profile and the AEC profile can make a considerable difference under pollution conditions. The SSA is obviously smaller on the North China Plain, where there are more absorptive aerosols. Therefore, in the calculation of the AMF, aerosol parameters such as AOD and SSA need to be taken into account. The aerosol correction can partly eliminate the systematic underestimation despite other more effective factors (HCHO shape factor and cloud effect). The explicit consideration of the contribution of aerosol is necessary, especially under a high load of absorptive aerosols.

5. Conclusions

In this work, HCHO VCD data derived from OMI by different algorithms provided by BIRA and NASA were validated using MAX-DOAS measurements at Xianghe station during the period 2010.03–2016.12. We compared the daily average OMI-BIRA and OMI-NASA data with the corresponding MAX-DOAS results under different eCF conditions and compared the OMI-BIRA and OMI-NASA datasets to each other. We found that the eCF has a significant effect on the daily comparison under the conditions of eCF < 0.3 (eCF < 0.1), with R = 0.26 (R = 0.34), R = 0.27 (R = 0.32) and R = 0.48 (R = 0.56) for OMI-BIRA and MAX-DOAS, OMI-NASA and MAX-DOAS, and OMI-BIRA and OMI-NASA, respectively. Furthermore, the results of the daily comparison by season show that both versions of the HCHO dataset have the best correlation with MAX-DOAS in spring and the poorest correlation in winter. The comparison of OMI-BIRA and OMI-NASA by season shows similar results for all seasons except for winter. The monthly (multi-year monthly) comparison shows a better correlation compared with the daily result, with R2 = 0.60 (R2 = 0.95) and R2 = 0.45 (R2 = 0.78) for OMI-BIRA and MAX-DOAS, and OMI-NASA and MAX-DOAS, respectively. In addition, the multi-year monthly RMSEs of OMI-BIRA and MAX-DOAS, and OMI-NASA and MAX-DOAS are significantly lower than the monthly RMSEs. Notably, the multi-year monthly average of OMI-BIRA is significantly lower than that of MAX-DOAS, ranging from −15.05% to −59.47%. OMI-NASA also tends to be lower than MAX-DOAS, but the differences are smaller than those between OMI-BIRA and MAX-DOAS.
The HCHO AMF sensitivity test showed that the difference in the AMF between different HCHO profile types is less than 20%. The differences in the AMF caused by AEC types are 5–8% and 12–23% for RAMS-CMAQ and CALIOP, respectively. However, the AMF variation is as high as 23–39% when the AEC and HCHO profile types are combined. In addition, we found that when the SSA is smaller and the AOD is larger, the calculated AMF is substantially smaller than that without the aerosol parameter (up to −53%). According to the result of the sensitivity test, we applied the aerosol correction to the multi-year monthly HCHO VCD data at Xianghe station. The results showed that after aerosol correction, the monthly HCHO changes by +6.37 to +20.7%.
In summary, we have provided HCHO usage recommendations for different algorithms at Xianghe station based on validation with MAX-DOAS data. The aerosol correction can be applied to the daily HCHO retrieval if the aerosol parameters are available.

Author Contributions

Conceptualization, L.C. (Liangfu Chen); Data curation, S.Z.; Formal analysis, Y.W.; Methodology, Y.W. and C.Y.; Project administration, Y.Z.; Resources, Z.W.; Validation, C.Y. and L.C. (Liangxiao Cheng); Writing—original draft, Y.W.

Funding

This research was funded by the National Key Research and Development Program of China (Grant No. 2017YFB0503901), the Strategic Priority Research Program of Chinese Academy of Sciences (Grant No. XDA19040201) and the National Natural Science Foundation of China (Grant No. 41771391).

Acknowledgments

We thank all the mission scientists and principal investigators who prepared and provided the OMI datasets used in this study. The authors are grateful to BIRA for open access to OMI HCHO retrieval results. We are sincerely grateful to Hendrick et al. for providing the MAX-DOAS HCHO dataset. We also thank Lin, J.T., et al. and Wu, T., et al. for providing the GEOS-Chem HCHO profile and CALIOP dataset.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Daily comparisons of OMI-BIRA and MAX-DOAS; left panels: 0 ≤ eCF ≤ 30%, middle panels: 0 ≤ eCF ≤ 20%, right panels: 0 ≤ eCF ≤ 10%. (a) Spring (MAM (March-April-May)), (b) Summer (JJA (June-July-August)), (c) Autumn (SON (September-October-November)), (d) Winter (DJF (December-January-February)).
Figure A1. Daily comparisons of OMI-BIRA and MAX-DOAS; left panels: 0 ≤ eCF ≤ 30%, middle panels: 0 ≤ eCF ≤ 20%, right panels: 0 ≤ eCF ≤ 10%. (a) Spring (MAM (March-April-May)), (b) Summer (JJA (June-July-August)), (c) Autumn (SON (September-October-November)), (d) Winter (DJF (December-January-February)).
Remotesensing 11 00203 g0a1
Figure A2. Daily comparisons of OMI-NASA and MAX-DOAS; left panels: 0 ≤ eCF ≤ 30%, middle panels: 0 ≤ eCF ≤ 20%, right panels: 0 ≤ eCF ≤ 10%. (a) Spring (MAM), (b) Summer (JJA), (c) Autumn (SON), (d) Winter (DJF).
Figure A2. Daily comparisons of OMI-NASA and MAX-DOAS; left panels: 0 ≤ eCF ≤ 30%, middle panels: 0 ≤ eCF ≤ 20%, right panels: 0 ≤ eCF ≤ 10%. (a) Spring (MAM), (b) Summer (JJA), (c) Autumn (SON), (d) Winter (DJF).
Remotesensing 11 00203 g0a2
Figure A3. Daily comparisons of OMI-NASA and OMI-BIRA; left panels: 0 ≤ eCF ≤ 30%, middle panels: 0 ≤ eCF ≤ 20%, right panels: 0 ≤ eCF ≤ 10%. (a) Spring (MAM), (b) Summer (JJA), (c) Autumn (SON), (d) Winter (DJF).
Figure A3. Daily comparisons of OMI-NASA and OMI-BIRA; left panels: 0 ≤ eCF ≤ 30%, middle panels: 0 ≤ eCF ≤ 20%, right panels: 0 ≤ eCF ≤ 10%. (a) Spring (MAM), (b) Summer (JJA), (c) Autumn (SON), (d) Winter (DJF).
Remotesensing 11 00203 g0a3
Figure A4. Monthly comparisons of OMI-BIRA and MAX-DOAS, OMI-NASA and MAX-DOAS, and OMI-NASA and OMI-BIRA. (a) 0 ≤ eCF ≤ 30%, (b) 0 ≤ eCF ≤ 20%, (c) 0 ≤ eCF ≤ 10%.
Figure A4. Monthly comparisons of OMI-BIRA and MAX-DOAS, OMI-NASA and MAX-DOAS, and OMI-NASA and OMI-BIRA. (a) 0 ≤ eCF ≤ 30%, (b) 0 ≤ eCF ≤ 20%, (c) 0 ≤ eCF ≤ 10%.
Remotesensing 11 00203 g0a4

References

  1. Deng, X.J.; Zhou, X.J.; Wu, D.; Tie, X.X.; Tan, H.B.; Li, F.; Bi, X.Y.; Deng, T.; Jiang, D.H. Effect of atmospheric aerosol on surface ozone variation over the Pearl River Delta region. Sci. China Earth Sci. 2011, 54, 744–752. [Google Scholar] [CrossRef]
  2. Han, X.; Zhu, L.Y.; Wang, S.L.; Meng, X.Y.; Zhang, M.G.; Hu, J. Modeling study of impacts on surface ozone of regional transport and emissions reductions over North China Plain in summer 2015. Atmos. Chem. Phys. 2018, 18, 12207–12221. [Google Scholar] [CrossRef]
  3. Ministry of Ecology and Environment of the People’s Republic of China. China Environmental Status Bulletin. Available online: http://www.mee.gov.cn/hjzl/zghjzkgb/lnzghjzkgb/ (accessed on 18 January 2019).
  4. Houweling, S.; Dentener, F.; Lelieveld, J. The impact of nonmethane hydrocarbon compounds on tropospheric photochemistry. J. Geophys. Res. Atmos. 1998, 103, 10673–10696. [Google Scholar] [CrossRef] [Green Version]
  5. Poisson, N.; Kanakidou, M.; Crutzen, P.J. Impact of non-methane hydrocarbons on tropospheric chemistry and the oxidizing power of the global troposphere: 3-dimensional modelling results. J. Atmos. Chem. 2000, 36, 157–230. [Google Scholar] [CrossRef]
  6. Cheng, N.; Chen, Z.; Sun, F.; Sun, R.; Dong, X.; Xie, X.; Xu, C. Ground ozone concentrations over Beijing from 2004 to 2015: Variation patterns, indicative precursors and effects of emission-reduction. Environ. Pollut. 2018, 237, 262–274. [Google Scholar] [CrossRef] [PubMed]
  7. Kanakidou, M.; Seinfeld, J.; Pandis, S.; Barnes, I.; Dentener, F.; Facchini, M.; Dingenen, R.V.; Ervens, B.; Nenes, A.; Nielsen, C. Organic aerosol and global climate modelling: A review. Atmos. Chem. Phys. 2005, 5, 1053–1123. [Google Scholar] [CrossRef]
  8. Duncan, B.N.; Yoshida, Y.; Olson, J.R.; Sillman, S.; Martin, R.V.; Lamsal, L.; Hu, Y.; Pickering, K.E.; Retscher, C.; Allen, D.J. Application of OMI observations to a space-based indicator of NOx and VOC controls on surface ozone formation. Atmos. Environ. 2010, 44, 2213–2223. [Google Scholar] [CrossRef]
  9. Parrish, D.; Ryerson, T.; Mellqvist, J.; Johansson, J.; Fried, A.; Richter, D.; Walega, J.; Washenfelder, R.D.; De Gouw, J.; Peischl, J. Primary and secondary sources of formaldehyde in urban atmospheres: Houston Texas region. Atmos. Chem. Phys. 2012, 12, 3273–3288. [Google Scholar] [CrossRef] [Green Version]
  10. De Smedt, I.; Müller, J.-F.; Stavrakou, T.; Van Der A, R.; Eskes, H.; Van Roozendael, M. Twelve years of global observations of formaldehyde in the troposphere using GOME and SCIAMACHY sensors. Atmos. Chem. Phys. 2008, 8, 4947–4963. [Google Scholar] [CrossRef] [Green Version]
  11. Anderson, L.G.; Lanning, J.A.; Barrell, R.; Miyagishima, J.; Jones, R.H.; Wolfe, P. Sources and sinks of formaldehyde and acetaldehyde: An analysis of Denver’s ambient concentration data. Atmos. Environ. 1996, 30, 2113–2123. [Google Scholar] [CrossRef]
  12. Chance, K.; Palmer, P.I.; Spurr, R.J.; Martin, R.V.; Kurosu, T.P.; Jacob, D.J. Satellite observations of formaldehyde over North America from GOME. Geophys. Res. Lett. 2000, 27, 3461–3464. [Google Scholar] [CrossRef] [Green Version]
  13. De Smedt, I.; Van Roozendael, M.; Stavrakou, T.; Müller, J.; Lerot, C.; Theys, N.; Valks, P.; Hao, N.; Van Der A, R. Improved retrieval of global tropospheric formaldehyde columns from GOME-2/MetOp-A addressing noise reduction and instrumental degradation issues. Atmos. Meas. Tech. 2012, 5, 2933–2949. [Google Scholar] [CrossRef] [Green Version]
  14. González Abad, G.; Liu, X.; Chance, K.; Wang, H.; Kurosu, T.; Suleiman, R. Updated Smithsonian Astrophysical Observatory Ozone Monitoring Instrument (SAO OMI) formaldehyde retrieval. Atmos. Meas. Tech. 2015, 8, 19–32. [Google Scholar] [CrossRef] [Green Version]
  15. De Smedt, I.; Stavrakou, T.; Hendrick, F.; Danckaert, T.; Vlemmix, T.; Pinardi, G.; Theys, N.; Lerot, C.; Gielen, C.; Vigouroux, C. Diurnal, seasonal and long-term variations of global formaldehyde columns inferred from combined OMI and GOME-2 observations. Atmos. Chem. Phys. 2015, 15, 12241–12300. [Google Scholar] [CrossRef]
  16. González Abad, G.; Vasilkov, A.; Seftor, C.; Liu, X.; Chance, K. Smithsonian Astrophysical Observatory Ozone Mapping and Profiler Suite (SAO OMPS) formaldehyde retrieval. Atmos. Meas. Tech. 2016, 9, 2797–2812. [Google Scholar] [CrossRef] [Green Version]
  17. Li, C.; Joiner, J.; Krotkov, N.A.; Dunlap, L. A new method for global retrievals of HCHO total columns from the Suomi National Polar-orbiting Partnership Ozone Mapping and Profiler Suite. Geophys. Res. Lett. 2015, 42, 2515–2522. [Google Scholar] [CrossRef] [Green Version]
  18. De Smedt, I.; Theys, N.; Yu, H.; Danckaert, T.; Lerot, C.; Compernolle, S.; Van Roozendael, M.; Richter, A.; Hilboll, A.; Peters, E.; et al. Algorithm theoretical baseline for formaldehyde retrievals from S5P TROPOMI and from the QA4ECV project. Atmos. Meas. Tech. 2018, 11, 2395–2426. [Google Scholar] [CrossRef] [Green Version]
  19. De Smedt, I.; Stavrakou, T.; Muller, J.F.; van der A, R.J.; Van Roozendael, M. Trend detection in satellite observations of formaldehyde tropospheric columns. Geophys. Res. Lett. 2010, 37, L18808. [Google Scholar] [CrossRef]
  20. Fu, T.M.; Jacob, D.J.; Palmer, P.I.; Chance, K.; Wang, Y.X.; Barletta, B.; Blake, D.R.; Stanton, J.C.; Pilling, M.J. Space-based formaldehyde measurements as constraints on volatile organic compound emissions in east and south Asia and implications for ozone. J. Geophys. Res. Atmos. 2007, 112. [Google Scholar] [CrossRef] [Green Version]
  21. Stavrakou, T.; Müller, J.-F.; Smedt, I.D.; Roozendael, M.V.; Van Der Werf, G.; Giglio, L.; Guenther, A. Global emissions of non-methane hydrocarbons deduced from SCIAMACHY formaldehyde columns through 2003–2006. Atmos. Chem. Phys. 2009, 9, 3663–3679. [Google Scholar] [CrossRef] [Green Version]
  22. Stavrakou, T.; Müller, J.-F.; Smedt, I.D.; Roozendael, M.V.; Van der Werf, G.; Giglio, L.; Guenther, A. Evaluating the performance of pyrogenic and biogenic emission inventories against one decade of space-based formaldehyde columns. Atmos. Chem. Phys. 2009, 9, 1037–1060. [Google Scholar] [CrossRef] [Green Version]
  23. Wang, Y.; Beirle, S.; Lampel, J.; Koukouli, M.; De Smedt, I.; Theys, N.; Li, A.; Wu, D.X.; Xie, P.H.; Liu, C.; et al. Validation of OMI, GOME-2A and GOME-2B tropospheric NO2, SO2 and HCHO products using MAX-DOAS observations from 2011 to 2014 in Wuxi, China: Investigation of the effects of priori profiles and aerosols on the satellite products. Atmos. Chem. Phys. 2017, 17, 5007–5033. [Google Scholar] [CrossRef]
  24. Lorente, A.; Boersma, K.F.; Yu, H.; Dörner, S.; Hilboll, A.; Richter, A.; Liu, M.; Lamsal, L.N.; Barkley, M.; De Smedt, I. Structural uncertainty in air mass factor calculation for NO2 and HCHO satellite retrievals. Atmos. Meas. Tech. 2017, 10, 759–782. [Google Scholar] [CrossRef]
  25. Lin, J.T.; Martin, R.V.; Boersma, K.F.; Sneep, M.; Stammes, P.; Spurr, R.; Wang, P.; Van Roozendael, M.; Clémer, K.; Irie, H. Retrieving tropospheric nitrogen dioxide over China from the Ozone Monitoring Instrument: Effects of aerosols, surface reflectance anisotropy and vertical profile of nitrogen dioxide. Atmos. Chem. Phys. 2014, 14, 1441–1461. [Google Scholar] [CrossRef]
  26. Lin, J.T.; Liu, M.Y.; Xin, J.Y.; Boersma, K.F.; Spurr, R.; Martin, R.; Zhang, Q. Influence of aerosols and surface reflectance on satellite NO2 retrieval: Seasonal and spatial characteristics and implications for NOx emission constraints. Atmos. Chem. Phys. 2015, 15, 12653–12714. [Google Scholar] [CrossRef]
  27. Boersma, K.F.; Eskes, H.J.; Dirksen, R.J.; Veefkind, J.P.; Stammes, P.; Huijnen, V.; Kleipool, Q.L.; Sneep, M.; Claas, J.; Leitão, J.; et al. An improved tropospheric NO2 column retrieval algorithm for the Ozone Monitoring Instrument. Atmos. Meas. Tech. 2011, 4, 2329–2388. [Google Scholar] [CrossRef]
  28. Hendrick, F.; Muller, J.F.; Clemer, K.; Wang, P.; De Maziere, M.; Fayt, C.; Gielen, C.; Hermans, C.; Ma, J.Z.; Pinardi, G.; et al. Four years of ground-based MAX-DOAS observations of HONO and NO2 in the Beijing area. Atmos. Chem. Phys. 2014, 14, 765–781. [Google Scholar] [CrossRef]
  29. Vlemmix, T.; Hendrick, F.; Pinardi, G.; De Smedt, I.; Fayt, C.; Hermans, C.; Piters, A.; Wang, P.; Levelt, P.; Van Roozendael, M. MAX-DOAS observations of aerosols, formaldehyde and nitrogen dioxide in the Beijing area: Comparison of two profile retrieval approaches. Atmos. Meas. Tech. 2015, 8, 941–963. [Google Scholar] [CrossRef]
  30. Wang, T.; Hendrick, F.; Wang, P.; Tang, G.; Clemer, K.; Yu, H.; Fayt, C.; Hermans, C.; Gielen, C.; Muller, J.F.; et al. Evaluation of tropospheric SO2 retrieved from MAX-DOAS measurements in Xianghe, China. Atmos. Chem. Phys. 2014, 14, 11149–11164. [Google Scholar] [CrossRef]
  31. Clemer, K.; Van Roozendael, M.; Fayt, C.; Hendrick, F.; Hermans, C.; Pinardi, G.; Spurr, R.; Wang, P.; De Maziere, M. Multiple wavelength retrieval of tropospheric aerosol optical properties from MAXDOAS measurements in Beijing. Atmos. Meas. Tech. 2010, 3, 863–878. [Google Scholar] [CrossRef] [Green Version]
  32. Zhu, L.; Jacob, D.J.; Kim, P.S.; Fisher, J.A.; Yu, K.; Travis, K.R.; Mickley, L.J.; Yantosca, R.M.; Sulprizio, M.P.; De Smedt, I.; et al. Observing atmospheric formaldehyde (HCHO) from space: Validation and intercomparison of six retrievals from four satellites (OMI, GOME2A, GOME2B, OMPS) with SEAC(4)RS aircraft observations over the southeast US. Atmos. Chem. Phys. 2016, 16, 13477–13490. [Google Scholar] [CrossRef]
  33. Barkley, M.; Abad, G.G.; Kurosu, T.P.; Spurr, R.; Torbatian, S.; Lerot, C. OMI air-quality monitoring over the Middle East. Atmos. Chem. Phys. 2017, 17, 4687–4709. [Google Scholar] [CrossRef] [Green Version]
  34. Rozanov, A.; Rozanov, V.; Buchwitz, M.; Kokhanovsky, A.; Burrows, J.P. SCIATRAN 2.0—A new radiative transfer model for geophysical applications in the 175–2400 nm spectral region. Adv. Space Res. 2005, 36, 1015–1019. [Google Scholar] [CrossRef]
  35. Rozanov, V.V.; Rozanov, A.V.; Kokhanovsky, A.A.; Burrows, J.P. Radiative transfer through terrestrial atmosphere and ocean: Software package SCIATRAN. J. Quant. Spectrosc. Radiat. Transf. 2014, 133, 13–71. [Google Scholar] [CrossRef]
  36. Wittrock, F. The Retrieval of Oxygenated Volatile Organic Compounds by remote Sensing Techniques. Doctoral dissertation, Universität Bremen, Bremen, Germany, 2006. [Google Scholar]
  37. Acarreta, J.; De Haan, J.; Stammes, P. Cloud pressure retrieval using the O2-O2 absorption band at 477 nm. J. Geophys. Res. Atmos. 2004, 109, D05204. [Google Scholar] [CrossRef]
  38. Xiao, H.; Zhang, M.; Zhu, L.; Liren, X.U. Model Analysis of Influences of Aerosol Mixing State upon Its Optical Properties in East Asia. Adv. Atmos. Sci. 2013, 30, 1201–1212. [Google Scholar]
  39. Han, X.; Zhang, M.; Gao, J.; Wang, S.; Chai, F. Modeling analysis of the seasonal characteristics of haze formation in Beijing. Atmos. Chem. Phys. 2014, 14, 10231–10248. [Google Scholar] [CrossRef] [Green Version]
  40. Byun, D.; Schere, K.L. Review of the Governing Equations, Computational Algorithms, and Other Components of the Models-3 Community Multiscale Air Quality (CMAQ) Modeling System. Appl. Mech. Rev. 2006, 59, 51–77. [Google Scholar] [CrossRef]
  41. Stephens, G.L.; Vane, D.G.; Boain, R.; Mace, G.; Sassen, K.; Wang, Z.; Illingworth, A.; Oconnor, E.; Rossow, W.; Durden, S.L. The Cloudsat Mission and the EOS Constellation: A New Dimension of Space-Based Observation of Clouds and Precipitation; NASA: Washington, DC, USA, 2001. [Google Scholar]
  42. Li, J.; Han, Z. Aerosol vertical distribution over east China from RIEMS-Chem simulation in comparison with CALIPSO measurements. Atmos. Environ. 2016, 143, 177–189. [Google Scholar] [CrossRef]
  43. Wu, T.; Fan, M.; Tao, J.; Su, L.; Wang, P.; Liu, D.; Li, M.; Han, X.; Chen, L. Aerosol Optical Properties over China from RAMS-CMAQ Model Compared with CALIOP Observations. Atmosphere 2017, 8, 201. [Google Scholar] [CrossRef]
  44. Torres, O.; Ahn, C.; Chen, Z. Improvements to the OMI near-UV aerosol algorithm using A-train CALIOP and AIRS observations. Atmos. Meas. Tech. 2013, 6, 3257–3270. [Google Scholar] [CrossRef] [Green Version]
  45. Ahn, C.; Torres, O.; Jethva, H. Assessment of OMI near-UV aerosol optical depth over land. J. Geophys. Res. Atmos. 2014, 119, 2457–2473. [Google Scholar] [CrossRef] [Green Version]
  46. Zhang, W.H.; Gu, X.F.; Xu, H.; Yu, T.; Zheng, F.J. Assessment of OMI near-UV aerosol optical depth over Central and East Asia. J. Geophys. Res. Atmos. 2016, 121, 382–398. [Google Scholar] [CrossRef]
  47. Kang, L.T.; Chen, S.Y.; Huang, J.P.; Zhao, S.M.; Ma, X.J.; Yuan, T.G.; Zhang, X.R.; Xie, T.T. The Spatial and Temporal Distributions of Absorbing Aerosols over East Asia. Remote Sens. 2017, 9, 20. [Google Scholar] [CrossRef]
  48. Lamsal, L.N.; Krotkov, N.A.; Celarier, E.A.; Swartz, W.H.; Pickering, K.E.; Bucsela, E.J.; Gleason, J.F.; Martin, R.V.; Philip, S.; Irie, H.; et al. Evaluation of OMI operational standard NO2 column retrievals using in situ and surface-based NO2 observations. Atmos. Chem. Phys. 2014, 14, 11587–11609. [Google Scholar] [CrossRef]
  49. Veefkind, J.P.; de Haan, J.F.; Sneep, M.; Levelt, P.F. Improvements to the OMI O2-O2 operational cloud algorithm and comparisons with ground-based radar-lidar observations. Atmos. Meas. Tech. 2016, 9, 6035–6049. [Google Scholar] [CrossRef]
  50. Theys, N.; Van Roozendael, M.; Hendrick, F.; Yang, X.; De Smedt, I.; Richter, A.; Begoin, M.; Errera, Q.; Johnston, P.V.; Kreher, K.; et al. Global observations of tropospheric BrO columns using GOME-2 satellite data. Atmos. Chem. Phys. 2011, 11, 1791–1811. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Location of Xianghe and the multi-year averaged formaldehyde (HCHO) over the North China Plain (2010–2016). Left: OMI-BIRA; right: OMI-NASA.
Figure 1. Location of Xianghe and the multi-year averaged formaldehyde (HCHO) over the North China Plain (2010–2016). Left: OMI-BIRA; right: OMI-NASA.
Remotesensing 11 00203 g001
Figure 2. Daily averaged and monthly averaged OMI-BIRA, OMI-NASA and MAX-DOAS data.
Figure 2. Daily averaged and monthly averaged OMI-BIRA, OMI-NASA and MAX-DOAS data.
Remotesensing 11 00203 g002
Figure 3. Daily comparisons of OMI-BIRA and MAX-DOAS, OMI-NASA and MAX-DOAS, and OMI-NASA and OMI-BIRA: (a) 0 ≤ eCF ≤ 30%, (b) 0 ≤ eCF ≤ 20%, (c) 0 ≤ eCF ≤ 10%.
Figure 3. Daily comparisons of OMI-BIRA and MAX-DOAS, OMI-NASA and MAX-DOAS, and OMI-NASA and OMI-BIRA: (a) 0 ≤ eCF ≤ 30%, (b) 0 ≤ eCF ≤ 20%, (c) 0 ≤ eCF ≤ 10%.
Remotesensing 11 00203 g003
Figure 4. Monthly comparisons of OMI-BIRA and MAX-DOAS, and OMI-NASA and MAX-DOAS: (a) 2010–2016 monthly average; (b) 2010–2016 multi-year monthly average.
Figure 4. Monthly comparisons of OMI-BIRA and MAX-DOAS, and OMI-NASA and MAX-DOAS: (a) 2010–2016 monthly average; (b) 2010–2016 multi-year monthly average.
Remotesensing 11 00203 g004
Figure 5. Multi-year monthly averaged HCHO concentrations for OMI-BIRA, MAX-DOAS and OMI-NASA.
Figure 5. Multi-year monthly averaged HCHO concentrations for OMI-BIRA, MAX-DOAS and OMI-NASA.
Remotesensing 11 00203 g005
Figure 6. (a) HCHO profile clusters (NHCHO: 5–8); (b) air mass factor (AMF_ corresponding to the HCHO profiles in panel (a). Each colour bar corresponds to the HCHO profile of the same colour.
Figure 6. (a) HCHO profile clusters (NHCHO: 5–8); (b) air mass factor (AMF_ corresponding to the HCHO profiles in panel (a). Each colour bar corresponds to the HCHO profile of the same colour.
Remotesensing 11 00203 g006
Figure 7. The aerosol extinction coefficient (AEC) clusters: (a) Regional Atmospheric Modeling System and Models-3 Community Multi-scale Air Quality (RAMS-CMAQ); (b) CALIPSO (NAEC: 5–8).
Figure 7. The aerosol extinction coefficient (AEC) clusters: (a) Regional Atmospheric Modeling System and Models-3 Community Multi-scale Air Quality (RAMS-CMAQ); (b) CALIPSO (NAEC: 5–8).
Remotesensing 11 00203 g007
Figure 8. AMF values corresponding to the AEC clusters. The colours correspond to the AEC profiles in Figure 7 for each cluster.
Figure 8. AMF values corresponding to the AEC clusters. The colours correspond to the AEC profiles in Figure 7 for each cluster.
Remotesensing 11 00203 g008
Figure 9. AMF values corresponding to the combination of the HCHO profile and AEC clusters. The colours correspond to the AEC profiles in Figure 7 for each cluster. The change in the x-axis corresponds to the HCHO profiles in Figure 6a for NHCHO = 5.
Figure 9. AMF values corresponding to the combination of the HCHO profile and AEC clusters. The colours correspond to the AEC profiles in Figure 7 for each cluster. The change in the x-axis corresponds to the HCHO profiles in Figure 6a for NHCHO = 5.
Remotesensing 11 00203 g009
Figure 10. AMF values corresponding to single scattering albedo (SSA) and aerosol optical depth (AOD): (a) absolute AMF; (b) relative change.
Figure 10. AMF values corresponding to single scattering albedo (SSA) and aerosol optical depth (AOD): (a) absolute AMF; (b) relative change.
Remotesensing 11 00203 g010
Figure 11. Monthly mean HCHO profile (a) and seasonal mean AEC (b) in Xianghe.
Figure 11. Monthly mean HCHO profile (a) and seasonal mean AEC (b) in Xianghe.
Remotesensing 11 00203 g011
Figure 12. Before and after aerosol correction.
Figure 12. Before and after aerosol correction.
Remotesensing 11 00203 g012
Table 1. Absolute and relative differences between OMI-BIRA, OMI-NASA and MAX-DOAS based on multi-year averages.
Table 1. Absolute and relative differences between OMI-BIRA, OMI-NASA and MAX-DOAS based on multi-year averages.
OMI-BIRA
(1015 molec/cm2)
MAX-DOAS
(1015 molec/cm2)
OMI-NASA
(1015 molec/cm2)
Difference
(OMI-BIRA& MAX_DOAS)
Difference
(OMI-NASA& MAX_DOAS)
HCHOHCHOHCHOAbs_DFRel_DF (%)Abs_DFRel_DF (%)
MON018.894.9610.227.7315.4312.62−1.345.21−15.0533.74
MON28.484.6510.517.2013.779.64−2.033.26−23.9723.66
MON38.956.0613.017.9312.686.21−4.06−0.33−45.39−2.62
MON49.505.9013.717.3712.345.00−4.20−1.37−44.22−11.10
MON513.817.2518.718.3516.268.14−4.90−2.45−35.45−15.05
MON615.936.7225.4011.5022.298.41−9.47−3.11−59.47−13.96
MON716.868.2423.677.2621.829.05−6.82−1.85−40.43−8.50
MON815.027.7522.667.6720.318.19−7.64−2.35−50.90−11.57
MON912.077.5218.517.6118.369.25−6.44−0.15−53.35−0.83
MON1011.848.9217.7810.7315.298.70−5.94−2.49−50.16−16.28
MON119.838.2413.928.9015.4110.55−4.081.49−41.529.69
MON129.175.1411.688.4014.939.73−2.513.25−27.4021.79
Table 2. OMI-BIRA and OMI-NASA algorithm setting comparison.
Table 2. OMI-BIRA and OMI-NASA algorithm setting comparison.
OMI-BIRAOMI-NASA
Chemical transport modelIMAGESGEOS-Chem
Background correctionDestriping correctionRow-dependent median HCHO in the equatorial Pacific (15°S–15°N,180–240°E)GEOS-Chem simulated HCHO in reference section (90°S–90°N, 140–160°W), with 2.0° (latitude) × 2.5° (longitude)
Zonal reference sector correctionA polynomial in the entire reference sector (90° S-90°N, 180–240°E)
BrO interference considerationPre-fit BrO slant columns, then fixed BrO values during HCHO retrievalImplicit in the GEOS-Chem simulation in reference section
Reference[15][14]

Share and Cite

MDPI and ACS Style

Wang, Y.; Wang, Z.; Yu, C.; Zhu, S.; Cheng, L.; Zhang, Y.; Chen, L. Validation of OMI HCHO Products Using MAX-DOAS observations from 2010 to 2016 in Xianghe, Beijing: Investigation of the Effects of Aerosols on Satellite Products. Remote Sens. 2019, 11, 203. https://doi.org/10.3390/rs11020203

AMA Style

Wang Y, Wang Z, Yu C, Zhu S, Cheng L, Zhang Y, Chen L. Validation of OMI HCHO Products Using MAX-DOAS observations from 2010 to 2016 in Xianghe, Beijing: Investigation of the Effects of Aerosols on Satellite Products. Remote Sensing. 2019; 11(2):203. https://doi.org/10.3390/rs11020203

Chicago/Turabian Style

Wang, Yapeng, Zifeng Wang, Chao Yu, Songyan Zhu, Liangxiao Cheng, Ying Zhang, and Liangfu Chen. 2019. "Validation of OMI HCHO Products Using MAX-DOAS observations from 2010 to 2016 in Xianghe, Beijing: Investigation of the Effects of Aerosols on Satellite Products" Remote Sensing 11, no. 2: 203. https://doi.org/10.3390/rs11020203

APA Style

Wang, Y., Wang, Z., Yu, C., Zhu, S., Cheng, L., Zhang, Y., & Chen, L. (2019). Validation of OMI HCHO Products Using MAX-DOAS observations from 2010 to 2016 in Xianghe, Beijing: Investigation of the Effects of Aerosols on Satellite Products. Remote Sensing, 11(2), 203. https://doi.org/10.3390/rs11020203

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop