Ahmad Iqbal1, Naveed Ahmad1, Hassan Mohy ud Din1, Michel Van Roozendael2, Muhammad Shehzaib Anjum1, Muhammad Zeeshan Ali Khan1, Muhammad Fahim Khokhar This email address is being protected from spambots. You need JavaScript enabled to view it.1 1 Institute of Environmental Science and Engineering (IESE), National University of Science and Technology (NUST), Islamabad-44000, Pakistan
2 Royal Belgian Institute for Space Aeronomy (B), Brussels, Belgium
Received:
December 28, 2021
Copyright The Author(s). This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are cited.
Revised:
April 22, 2022
Accepted:
April 30, 2022
Download Citation:
||https://doi.org/10.4209/aaqr.210398
Iqbal, A., Ahmad, N., Mohy ud Din, H., Van Roozendael, M., Anjum, M.S., Zeeshan Ali Khan, M., Khokhar, M.F. (2022). Retrieval of NO2 Columns by Exploiting MAX-DOAS Observations and Comparison with OMI and TROPOMI Data during the Time Period of 2015–2019. Aerosol Air Qual. Res. 22, 210398. https://doi.org/10.4209/aaqr.210398
Cite this article:
Nitrogen dioxide (NO2)—a criteria major air pollutant, is of paramount importance due to its role in atmospheric chemistry and tropospheric ozone formation. Exposure to high concentrations of NO2 has been reported to cause various health issues in humans. This study presents an intercomparison of NO2 retrieval settings using the Differential Optical Absorption Spectroscopy (DOAS) and based on the literature published over the last 20 years. Comparison between results of various settings, as reported in the literature, and settings used for this study show a good correlation with R2 > 0.97. This paper also presents validation of satellite observation through ground-based MAX-DOAS measurements from September 2015 to September 2019. Daily MAX-DOAS measurements have shown a strong positive correlation of 70.75% and 77.74% with OMI and TROPOMI, respectively. The average monthly correlation for OMI and TROPOMI with MAX-DOAS is 88.39% and 91.91% respectively. The comparison of the slopes of regression plots for daily and monthly datasets of OMI and TROPOMI vs. MAX-DOAS reveals that TROPOMI data is more synonymous to MAX-DOAS than OMI data. The error analysis indicates that for TROPOMI measurements calculated biases are significantly improved in case of TROPOMI as compared to OMI measurements. It is pertinent to mention that TROPOMI measurements can capture the local NO2 pollution much better than OMI and its predecessor instruments like GOME-2, SCIAMACHY and GOME. Seasonal trends of NO2 column densities have shown a peak in the winter season (November–January) while lowest NO2 column density is recorded in monsoon season.HIGHLIGHTS
ABSTRACT
Keywords:
NO2 retrieval settings, DOAS Algorithm, Satellite observations, Validations
Nitrogen dioxide (NO2) is of prime importance due to its active role in the photochemistry of the atmosphere (Ma et al., 2013; Chan et al., 2015). NO2 is also one of the major criteria air pollutants, which poses serious threats to human health (WHO, 2006). Exposure to high concentrations of nitrogen dioxide has been reported to cause damage to the respiratory system and aggravates bronchitis in some cases (Valavanidis et al., 2008; Ma et al., 2013; Khokhar et al., 2017). Furthermore, exposure to NO2 has been considered the root cause of reduced lungs efficiency in infants (Trasande and Thurston, 2005) and exacerbating asthma among primary school children (Youssef Agha et al., 2012). However, higher incidents of lung diseases have been associated with long-term exposure to NO2 (Soto-Martinez and Sly, 2010). Especially in developing countries like Pakistan, lack of political will, inadequate resources, and air quality monitoring infrastructure combined with weak enforcement of air pollution standards have primarily increased the vulnerability of humans to air pollution-related diseases (Ali et al., 2021; Anjum et al., 2021). The importance of nitrogen dioxide in the atmosphere is not only due to its health effects but also because (a) it absorbs solar radiation and reduces atmospheric visibility (Majewski et al., 2014); (b) it affects the oxidizing capacity of the troposphere by reacting with radical species, like hydroxyl radical, present in the atmosphere, (c) it plays a crucial role in maintaining ozone concentration in the troposphere as it initiates the catalytic formation of ozone (U.S. EPA, 2006). Sources of NO2 include both natural and anthropogenic activities. But most of these sources are related to the burning of biomass, fossil fuel consumption, and different industrial emissions (Tack et al., 2015; Gratsea et al., 2016; Chan, 2017; Malik et al., 2020). On the other hand, it is mainly removed from the atmosphere through its reaction with free radicals like hydroxyl radicals (OH) and ozone (O3), which ultimately results in Nitric acid and NO3 (Crutzen, 1970; Athanasopoulou et al., 2008; Wang et al., 2012; Yang et al., 2019). In recent years, optical methods were used to perform NO2 measurements in the atmosphere. The advantages of this method are high sensitivity, fast measurements and simultaneous detection of several trace gases (Lohberger et al., 2004). In particular, multi-axis DOAS (MAX-DOAS) instruments (Hönninger et al., 2004; Platt and Stutz, 2008; Vlemmix et al., 2015) have been used for the retrieval of vertical column densities of various tropospheric trace gases (NO2, O3, CHOCHO, SO2, BRO, IO, HCHO and HONO) (Leser et al., 2003; Coburn et al., 2011; Hendrick et al., 2014; Wang et al., 2017; Khan et al., 2018; Chan et al., 2019; Kreher et al., 2019; Javed et al., 2019a, b). MAX-DOAS measurements can be performed from different platforms like ground-based, vehicle-based, ship-based and satellite-based (Halla et al., 2011; Peters et al., 2012; Lin et al., 2014; Shaiganfar et al., 2015; Dix et al., 2016; Shabbir et al., 2016). Cheng et al. (2019) explored the temporal variation of tropospheric NO2 vertical column densities (VCDs) using MAX-DOAS measurements in Shangdianzi region of China and compared them with OMI data. In another study, MAX-DOAS has been used for the observation of tropospheric SO2, NO2, and HCHO in Hefei, Nanjing, and Shanghai (Yangtze River Delta area) (Tian et al., 2018). Longer-term studies covering more than 5 years were also conducted using Max-DOAS (Chan et al., 2018). In Pakistan, air pollution is a predictable key environmental issue as it is ranked in the top polluted countries in the region. The recurrence of smog episodes in Pakistan and neighboring regions is a clear manifestation of worsened air quality and consequent disruption in socio-economic circles (Khokhar et al., 2021; Anwar et al., 2021). The menace of degraded air quality is aided by economic development, unprecedented increase in the vehicular fleet, poor quality of fuel, and lack of any satisfactory measures due to the absence of a proper air quality monitoring network. Large-scale mitigation measures are direly needed to be taken to maintain ambient air quality within permissible limits. Very few studies have been conducted in Pakistan to measure the atmospheric levels of NO2 in Pakistan (Ali and Athar, 2008; Tariq et al., 2014; Khokhar et al., 2016b; Zeb et al., 2019; Bilal et al., 2021). For instance, Zeb et al. (2019) exploited multiple satellite instruments to identify the spatial and decadal variation of trace gases in Pakistan. In a similar study, Shabbir et al. (2016) have used Car-mounted MAX-DOAS to measure the NO2 concentration over the NH5 highway from Islamabad to Lahore, in Pakistan. However, so far, no validation has been performed in the moderately polluted city of Islamabad, where temporal gradients of NO2 are small. Thus, the objective of this study is to present (a) the assessment of a DOAS fit NO2 retrieval suitable for measurements at Islamabad (b) an intercomparison between a large number of NO2 DOAS retrieval settings consistently applied to the same spectral data, and (c) the use of ground-based tropospheric NO2 measurements performed at Islamabad to validate satellite observations from the OMI and TROPOMI instruments. This paper is organized as follows In Section 2, the design of the intercomparison study is discussed including the instrument, study site, and measurement protocol. In Section 3, results of a regression analysis of various NO2 retrieval settings are presented. In Section 4, the validation of OMI and TROPOMI satellite observations using ground-based MAX-DOAS data are discussed. Finally, Section 5 provides a conclusion of the study. Various studies have been conducted previously for the retrieval of NO2 measurements by the MAX-DOAS technique. In this study, settings extracted from the literature were used to conduct an intercomparison of DOAS fit settings performed on data recorded on 1 November 2017. We used results from this exercise, to define optimal settings for NO2 measurements at Islamabad, Pakistan. Table S1 contains a list of the settings that were used for DOAS analysis on 1 November 2017. The table lists the studies included, Wavelength ranges for each study, and absorber fitted used in a specific study. It also includes various parameters retrieved e.g., Maximum dSCD (molecules cm–2), Average dSCD (molecules cm–2), RMS, spectrum shift and dSCD error (molecules cm–2). The DOAS technique used in this study is based on the Lambert-Beer law which states that the transmittance of light through an absorbing medium is directly proportional to the concentration of the substance present in the optical path. Here, in Eq. (1) “Io” means incident intensity whereas “I” is referred to the measured intensity. Various trace gases can be obtained simultaneously within a selected fitting interval. The retrieval method considers laboratory absorption cross-sections and by using spectra of the desired gases, slant column densities (SCDs) can be retrieved using a linear-least squares minimization. The DOAS technique is especially used for the retrieval of trace gases that have a short lifetime in the atmosphere like NO2. QDOAS is retrieval software designed by the Royal Belgium Institute for Space and Aeronomy for the analysis of UV-visible spectra. Different cross-sections: NO2 at 298 K (Vandaele et al., 1996), NO2 a at 220 K (Vandaele et al., 1996), O3 at 223 K (Serdyuchenko et al., 2014), O4 at 293 K (Thalman and Volkamer, 2013), H2O (Rothman, 2010) and Ring were used in the fitting interval of 409–445 nm for the development of the optimized NO2 retrieval settings considered in this study (see Table S1). Measurements were carried out at elevation angles 2°, 4°, 5°, 10°, 15°, 30°, 45°, and 90° and polynomial order of 4 was selected. To minimize the stratospheric influence, the zenith spectrum of each cycle was used. The residual spectrum that remains after subtracting all absorption cross-sections through least-squares minimization is characterized by its root mean square (RMS). It is generally of the order of 10–4 or higher for good instruments. Different trace gases cross-sections and aerosols retrieval play important role in the uncertainty of the vertical column densities (VCDs). Sites with low trace gases can have higher uncertainty in retrieved VCDs. The uncertainties introduced by the geometric approximation depend on the viewing geometry and the aerosol load. However, instrumental conditions can also affect the uncertainty and accuracy of the results. Many studies conducted previously have revealed that typical error in VCDs retrieved by MAX-DOAS is around 20% (Hendrick et al., 2014; Wang et al., 2014; Franco et al., 2015). Fig. 1 shows an example of a DOAS fit for NO2. Red lines represent the reference spectrum and the black line shows the measured spectrum after subtracting all other absorbers. It was measured on 23 May 2018 at 12:36 UTC. SZA was 73° and the elevation angle was 5°. The mean dSCD value for this fit was 3.55 × 1016 molecules cm–2. Differential slant column densities depend upon the light pathways which is characterized by an air mass factor (AMF) usually calculated using a radiative transfer model. In this study, however the AMF was estimated using a simple geometric approximation as described in (Hönninger et al., 2004; Wittrock et al., 2004; Wagner et al., 2010). This method allows the retrieval of column densities of trace gases even in cloudy conditions (Wagner et al., 2010). A differential air mass factor (dAMF) was used because the slant column densities obtained were differential SCDs corresponding to the difference between SCDs measured at low elevations and dSCDs recorded at the zenith for the same scan (Liu et al., 2016). dAMF is difference in AMF obtained at a certain angle and AMF obtained at 90°: So, Eq. (3) will be: Using geometric approximation, AMF can be found as: Then Eq. (5) will be: This method has been used widely to derive the tropospheric vertical column (VCD) of trace gases. Despite its simplicity, it works with good accuracy for NO2 when applied to measurements at 30° of elevation. To quantify the tropospheric trace gases, the MAX-DOAS uses different elevation angles, from 0 to 90°, to capture the scattered sunlight (Bobrowski et al., 2003). MAX-DOAS can measure several trace gases, at a time, in the visible and Ultra-Violet Spectral ranges and the resultant less residual allows its use in even less polluted environment. The MAX-DOAS instrument used in this study is based on a light weight spectrometer with dimensions of 13 cm × 19 cm × 14 cm. It contains a 40 mm lens mounted at the front, coupled with the spectrograph by fiber optics and an electronic circuit. Czerny-Turner spectrometer (USB 2000+, Ocean Optics In.), having a spectral resolution of 0.7 nm, records the solar scattered light between 320–460 nm. A stepper motor, installed in it, is used to move the entire instrument to achieve the desired elevation angle. Scattered light enters the instrument via a lens that is tightly sealed to avoid any condensation of water vapors and prevent dust particles to damage the interior optics. In this study, DOAS fit retrieval settings were used to analyze the MAX-DOAS measurements from September 2015 to September 2019. The OMI observations (https://disc.gsfc.nasa.gov/datasets/OMNO2d_003/summary) for this period were obtained from GIOVANNI and the data product used was “NO2 Tropospheric column 30% cloud screened (OMNO2d v003)” for better agreement of tropospheric NO2. However, later on, the tropospheric monitoring instrument (TROPOMI) data (https://developers.google.com/earth-engine/datasets/catalog/COPERNICUS_S5P_NRTI_L3_NO2) was also available starting from August 2018. TROPOMI product of tropospheric NO2 column number density (having a spatial resolution of 7 × 3.5 km was taken through the Google Earth Engine platform (ee.ImageCollection("COPERNICUS/S5P/NRTI/L3_NO2")) from September 2018 to September 2019. TROPOMI is a satellite instrument onboard Sentinel 5P satellite. Slant Column Densities obtained from OMI and their uncertainties have been affected by the instrument’s degradation (radiometric), row irregularity and stripes. According to a study, temporal analysis of NO2 Slant column Density uncertainty over the period of 2005–2016 has been 0.013 × 1015 ± 0.001 × 1015 molec. cm–2 year–1 (Zara et al., 2018). Uptil now, less than 5% optical degradation in visible channel of OMI has been observed (Boersma et al., 2011; Müller et al., 2016; Schenkeveld et al., 2017). Similarly, noise is also increasing gradually over the period of time owing to the long-term use of CCD detector (Schenkeveld et al., 2017). Furthermore, the details on the slant column retrieval method, stability, uncertainties for TROPOMI NO2 measurements and its comparison with OMI are available in a study conducted by Van Geffen et al. (2020). Islamabad, the capital of Pakistan was selected as a site for this study. Pakistan’s climate is described by four seasons based on the occurrence of monsoon: Winter (December–February), Pre-Monsoon (March–May), Monsoon (June–September) and Post Monsoon (October–November) (Naheed et al., 2013). Islamabad is located at the latitude of 33.349°N; and longitude of 72.324°E, at an elevation above sea level of 450 meters (1827 ft). Islamabad is a mountainous region located in the Pothohar Plateau. The instrument was mounted on the rooftop of the Institute of Environmental Sciences and Engineering (IESE), National University of Sciences and Technology (NUST) Islamabad, Pakistan. The lowest elevation angle used in this study was 2°. A sensitive study was constructed to find optimal settings for retrieval of NO2 using MAX-DOAS for local climatic conditions. The efforts for this study are inspired by the QA4ECV project (http://www.qa4ecv.eu/) in which algorithms were improved for satellite NO2 retrievals. Similarly, this study attempts to find optimum settings for NO2 retrieval through ground-based MAX-DOAS for local conditions in Islamabad. Though sensitivity to NO2 spectral lines is generally known to be highest in the 400–450 nm region, analysis were carried out to find the best fitting sub-window in this range that could further improve the results and reduce the RMS and Slant column error of the data obtained at the site. Moreover, different O4 cross-sections were also analyzed to find out the most suitable one for the NO2 retrieval in the study area. These were Thalman and Volkamer (2013) at 293 K, Hermans (2011) at 298 K, Thalman and Volkamer (2013) at 203 K and Burkholder and Talukdar (1994). Among these, Burkholder cross-section improved the O4 best fit. Moreover, RMS and dSCDs were also improved by using Burkholder cross-section. In addition to this, a wavelength range of 320 nm to 445 nm was used keeping in view the spectral range of an available MAX-DOAS instrument. Initially, the upper boundary of the fitting interval was kept constant at 445nm whereas the lower boundary was increased up to 420 nm with an interval of 10 nm each. Then lower boundary was kept constant at 320 nm and the upper boundary was reduced from 445nm up to 355 nm with an interval of 10 nm each. Consequently, RMS, Slant error and dSCDs of each fitting interval were compared and the least RMS and Slant error was observed at 410 nm to 445nm interval as the lower and upper boundary of fitting interval respectively. Further analysis with 409 nm–445 nm as fitting window and Burkholder O4 cross-section for NO2 retrieval shows reduction of RMS and NO2 slant column error. Thus, this setting was considered as optimum for the study area i.e., Islamabad. Moreover, to compare the settings of this study, literature was consulted with studies having fitting windows within the range of fitting interval of the instrument used in this study i.e., 320 nm–465 nm. Thus, the reported settings of 19 studies (for more details, see Table S1) were used to analyze and compare the results with the settings found in this study. Fig. 2(a) shows dSCDs of NO2 obtained at an elevation angle 10° on 1 November 2017. It can be observed that dSCDs obtained using optimal settings of this study are among the highest as compared to dSCDs obtained by using the settings of previous studies reported in the literature. dSCDs obtained using Shaiganfar et al. (2011) settings show the highest values but corresponding RMS are larger than those obtained using the settings of this study. The comparison of the NO2 RMS obtained using the various settings considered for this study is shown in Fig. 2(b). It clearly shows that the RMS for dSCDs obtained by using the settings of this study is low as compared to that of most of the others. Similarly, Fig. 2(c) shows the slant column error for the NO2 dSCDs measured in this study is among the lowest. The DOAS fit retrieval settings for NO2 used in this study have not only reduced Slant error and RMS but also improved dSCDs which depicts that these settings are suitable for the study site. To intercompare the different data sets, linear regressions were performed using data sets obtained with standard (optimized) settings of this study as a reference. Accordingly, dSCDs retrieved from different settings found in the literature were used to draw correlation plots as shown in Fig. 3. As can be seen, a tight correlation is observed for most of the cases. However, the highest correlation (R2 = 0.999) is found with the settings of Drosoglou et al. (2018) which includes a strong NO2 absorption band of 411–445 nm. Whereas, Lee et al. (2009) have shown the lowest value (R2 = 0.972) which includes an absorption band of 364–383 nm (for more details see Table S1). It can be attributed to the wavelength dependency of tropospheric NO2 AMFs. For instance, Lee et al. (2009) has used a spectral fitting window of (364–383 nm) which falls in the UV-A spectral range and identified a longer light path and higher extinction coefficient. Furthermore, it might be largely impacted by the Rayleigh scattering dominant towards a shorter wavelength. However, wavelength window 409–445nm, was used in this study for further retrieval of NO2 dSCDs. The values for the slope, intercept and RMS are calculated as a part of the regression analysis. Whereas, Spectrum shift is obtained from the DOAS retrieval algorithm used for MAX-DOAS data analysis. Fig. 4 presents a different view of the data set already presented in Fig. 3. It shows RMS, slope, intercept, and spectrum shift calculated from regression analysis for all the measurements obtained by different settings with settings found in this study. It was done for different elevation viewing angles (2°, 4°, 5°, 10°, 15°, 30°, and 45°). Values are colored according to the elevation viewing angles. The dash lines indicate the standard deviation for slope, intercept, and spectrum shift whereas in RMS it represents the average value. A close relation of the parameters is observed among all the settings used for different elevation angles. Three parameters of regression used in Fig. 4 (RMS, slope, and intercept) do not display any significant angular dependency. Some data display larger deviations irrespective of elevation angles. For instance, RMS for Tian et al. (2018), Chan et al. (2015), and Lee et al. (2009) lie above the mean value. Whereas, for slope, Shaiganfar et al. (2011), Wang et al. (2017), Leser et al. (2003), Chan et al. (2020), Dimitropoulou et al. (2020) and Kreher et al. (2020) showed largest deviations. However, a larger Spectrum Shift is observed at smaller angles as shown in Fig. 4. The spatial extent and magnitude of tropospheric NO2 columns over Islamabad are highlighted in Fig. 5. The lower two panels show the OMI obtained area-averaged maps for Islamabad. Where the bottom map represents the time periods between October 2004 to September. 2019 and the second panel represents the period between September 2018 and September 2019. A very minute difference can be seen between the two maps showing that there has been very little significant change in the spatial extent of NO2. This could also be a result of the coarse resolution of the OMI. The first panel of Fig. 5 highlights the TROPOMI Tropospheric NO2 columns, also for the time period September 2018 to 2019. It is apparent that the better spatial resolution (7 × 3.5 km2 during the study duration for which data is used) of the TRPOMI instrument helps differentiate between hotspots more easily. In this case, the higher levels of NO2 are observed near congested roads and densely packed areas of Islamabad city while peri-urban and rural areas in the vicinity of Islamabad are showing relatively lower NO2 columns. However, it is important to note that there is not a stark difference between the distributions of pollution levels for the long-term averaged map and the yearly map from OMI instrument. This urges the need for comprehensive ground-based observations for the validation of satellite products that can be used for pollution management and to devise abatement strategies. Thus, this study having NO2 columns observed by the MAX-DOAS instrument at the IESE-NUST site in combination with satellite data can serve such purposes to a certain extent. Especially, areas without having a regular air pollution monitoring network. The comparison of available data for all three instruments i.e., ground-based, and two satellites is presented in Figs. 6(a) and 6(b). The data for ground-based MAX-DOAS and OMI was available from September 2015 to September 2019 except for January and February 2016. The data for TROPOMI was available from September 2018 to September 2019. Initially, the on-site ground-based MAX-DOAS VCDs for the study period were averaged over a local time of 1300 hrs to 1400 hrs on daily basis as the TROPOMI and OMI over passing time are 1330 hrs and 1340 hrs respectively. So, averaged MAX-DOAS VCDs (1300 hrs–1400 hrs local time) were used for satellite validation. Further, the data was filtered on daily basis keeping in view the cloud cover of less than 20% and to make sure the comparison of three instruments for same-day i.e., the comparison is shown for all those days for which data is available for all three instruments (see Fig. 6(a)). The comparison for daily measurements of three instruments i.e., ground-based MAX-DOAS, OMI, and TROPOMI is presented in Fig. 6(a) while the monthly comparison for these instruments is shown in Fig. 6(b). Fig. 6(a) shows that daily satellite observations (OMI and TROPOMI) are corresponding to MAX-DOAS observations and showing a similar trend for the entire period but values are lower as compared to the measurements obtained by ground-based MAX-DOAS. However, the values obtained for TROPOMI are higher and showed good agreement with MAX-DOAS as compared to OMI. It is mainly because of the higher spatial resolution of TROPOMI (7 × 3.5 km2) as compared to OMI (13 × 24 km2) in addition to better spectral resolution i.e., 0.54 nm and 0.63 nm for TROPOMI and OMI respectively, and improved signal to noise ratio of around 1500 for TROPOMI and around 500 for OMI. Therefore, TROPOMI measurements are able to capture the local NO2 pollution much better than OMI. However, satellite instruments quantitively underestimated NO2 columns (see Fig. 6(a)) for the entire period as compared to ground-based MAX-DOAS amidst to the fact that satellites average over larger ground pixel sizes and are relatively less sensitive to boundary layer pollution. Moreover, satellite observations through OMI and TROPOMI are showing the same trend as that of ground observations through MAX-DOAS (see Fig. 6(b)). Also, owing to its high spatial resolution, it is pertinent to mention that TROPOMI measurements are able to capture the local NO2 pollution much better than OMI and its predecessor instruments like GOME-2, SCIAMACHY, GOME. Considering the data from three instruments, linear regressions were performed for daily as well as monthly measurements of these instruments. Linear regression was performed for ground-based MAX DOAS vs. OMI and ground-based MAX-DOAS vs. TROPOMI. The correlation for daily measurements between MAX-DOAS and OMI has shown a very good positive correlation with R2 = 0.50 (r = 0.70) as shown in Fig. 7(a). Fig. 7(b) shows a positive correlation between daily NO2 measurements of ground-based MAX-DOAS and TROPOMI. This shows an improved correlation with R2 = 0.60 (r = 0.77) which could be due to improved spatial resolution of TROPOMI as compared to OMI. The monthly correlation of ground-based MAX-DOAS with OMI and TROPOMI is shown in Figs. 7(c) and 7(d) respectively. It has shown a high positive correlation with R2 = 0.78 (r = 0.88) for OMI (see Fig. 7(c)) While it also depicts an improved and high positive correlation with R2 = 0.84 (r = 0.91) for TROPOMI (see Fig. 7(d)). Through linear regression and correlation plots, it can be seen that the instruments i.e., ground-based MAX-DOAS vs. OMI and ground-based MAX-DOAS vs. TROPOMI have shown a strong positive correlation for daily as well as for monthly NO2 measurements. The slope of the regression analysis between MAX-DOAS and TROPOMI depicts a linear positive relationship between the two (see Fig. 7(b)). The value of slope for this model is 0.1893 for TROPOMI. The slope value shows underestimation by TROPOMI as the slope is nearer to 0 than 1 (i.e., 0.1). On the other hand, values obtained from the OMI satellite shows the same trend. But its slope value is indicating an increase of 0.1492 as compared to 0.1893 of TROPOMI (see Fig. 7(a)). However, the monthly correlation of a MAX-DOAS with OMI and TROPOMI data not only improves R2 significantly but also improves the slope. The slope with OMI data improved from 0.14 to 0.18 whereas, the slope of TROPOMI improved significantly from 0.18 to 0.40 as shown in Figs. 7(c) and 7(d) for OMI and TROPOMI respectively. The comparison of the slopes of OMI and TROPOMI vs. MAX-DOAS plots reveals that TROPOMI data is more synonymous to that of data obtained by MAX-DOAS. It can be attributed mainly owing to a higher and improved spatial and spectral resolution of TROPOMI than that of OMI. The slope of the correlation plots also exhibits underestimation of the satellite data and it is mainly due to restrictions of satellite sensors and retrieval algorithms. It can be further attributed to uncertainties caused by the coarse resolution of satellite sensors as compared to ground-based MAX-DOAS. Although satellite data underestimates the NO2 column densities, it shows similar variation in NO2 columns and can be speculated that there is a good agreement of MAX-DOAS measurements with TROPOMI data, instead of OMI. Furthermore, the seasonal variations of NO2 column densities for all three instruments are identical. As shown in Figs. 6(a) and 6(b) NO2 column density peaks in the winter season while the lower NO2 columns are observed in the monsoon season. Depending upon various factors i.e., OH concentration, photolysis rate and varying meteorological conditions, NO2 has a lifetime of a few hours to about one day (Beirle et al., 2003; Seinfeld and Pandis, 1997; Sheel et al., 2010). The largest known sink of NOx is a reaction of NO2 with hydroxyl (OH) radical (Stavrakou et al., 2013). As there is low temperature and limited solar flux is available during the winter season, there is less NO2 loss due to reaction with OH radical that leads to higher NO2 lifetime and column densities. The other contributing factor is fewer precipitation rates that hinder frequent OH production due to less availability of water vapors (Sheel et al., 2010). This leads to relatively limited removal of NO2 molecules. Furthermore, during the winter season increased usage of wood, cow dung and coal burning for cooking and space heating results in increased NO2 emissions. Due to the longer NO2 lifetime and more emissions of NO2 during the winter season, higher NO2 values are observed in winter season (Martin et al., 2010). A similar trend can be observed in Fig. 6(b). Grey rectangles in Fig. 6(b) show the monsoon season in Pakistan and longer and heavier spells of monsoon rains. Therefore, the monsoon exhibits lower NO2 columns in the seasonal cycle mainly due to the larger availability of OH radical causing a decreased NO2 lifetime. This could also be attributed to the advection of moist clean air masses, heavy precipitation (wet scavenging) and strong actinic fluxes during the monsoon season (Wang and Jacob, 1998; Jacob, 2003; Khokhar et al., 2015). The observations from satellite as well as ground-based instruments recorded the highest NO2 column densities when days are longer than night and photolysis of NO2 occurs for a much longer time. Error analysis between the MAX-DOAS measurements and satellite observations i.e., OMI and TROPOMI were performed by finding the root mean square error (RMSE) and the mean absolute error (MAE) given by Eqs. (8) and (9), respectively. where n is the number of observations. Besides computing differences in averages, overestimation or underestimation of a model can be quantified by calculating the mean bias (MB), which is given by Eq. (10): Error analysis computed over a different time duration under the study period on daily as well as a monthly basis is presented in Table 1. The daily RMSE, MAE and MB between MAX-DOAS and OMI for the study period (2015–2019) were 1.08E+16, 8.82E+15 and 8.59E+15 respectively while on monthly basis these were 9.46E+15, 8.74E+15 and 8.74+15 respectively. The daily and monthly analysis of RMSE, MAE and MB for the period of 2018–2019 were also performed to check the synchronization of satellite observations (TROPOMI and OMI) with ground-based MAX-DOAS measurements. The daily RMSE, MAE and MB between MAX-DOAS and OMI for the period of 2018–2019 were 8.87E+15, 7.34E+15 and 7.24E+15 respectively while on monthly basis these were 9.42E+15, 8.09E+15 and 8.09+15 respectively. Similarly, the daily RMSE, MAE and MB between MAX-DOAS and TROPOMI for the period of 2018-2019 were 4.49E+15, 3.40E+15 and 2.15E+15 respectively while on monthly basis these were 4.34E+15, 3.04E+15 and 2.76E+15 respectively. These error analysis of MAX-DOAS vs. OMI and TROPOMI on a daily as well as monthly basis indicates that TROPOMI measurements are more synchronized with ground-based MAX-DOAS measurements. It is pertinent to note that TROPOMI’s finer spatial resolution has significantly improved the biases. We acknowledge the free use of Tropospheric NO2 column data from the TROPOMI and OMI sensors. We acknowledge NUST Islamabad, Pakistan for providing partial financial support as a Master thesis research fund to conduct this study. Very special gratitude goes to the satellite group of Max-Planck Institute for Chemistry Mainz, Germany for providing the mini MAX-DOAS instrument.1 INTRODUCTION
2 DESIGN OF INTERCOMPARISON STUDY AND MEASUREMENT PROTOCOL
2.1 Introduction
2.2 Differential Optical Absorption Spectrometry (DOAS) Technique
2.3 Measurement and Data Analysis Protocol for the Development of New Retrieval Settings for NO2Fig. 1. The DOAS fit observed for NO2 on 23 May 2018 at 5° elevation angle and SZA of 73°.
2.4 Instrument
2.5 Measurement Site
3 RESULTS
3.1 Comparison of dSCD, RMS, and Slant Error Retrieved for Intercomparison of NO2 DOAS Fit Settings Listed in Table S1Fig. 2. Comparison of (a) NO2 dSCDs, (b) RMS and (c) Slant Column error retrieved by MAX-DOAS through various DOAS fit settings used in literature and settings used for this study.
3.2 Regression AnalysisFig. 3. Regression analysis of NO2 dSCDs for various DOAS fit retrieval settings found in literature plotted against the reference DOAS fit retrieval settings used for this study.
Fig. 4. Slope, intercept, spectrum shift and RMS of regression analysis for dSCDs of various DOAS fit retrieval settings found in literature plotted against the reference DOAS fit retrieval settings of this study.
4 VALIDATION OF SATELLITE TROPOSPHERIC NO2 OBSERVATIONS
Fig. 5. Spatial and temporal extent of NO2 columns over Islamabad using OMI and TROPOMI data.
Fig. 6. Comparison of (a) daily (b) monthly vertical column densities (VCDs) retrieved through ground-based MAX-DOAS, tropospheric NO2 column densities by OMI, and tropospheric NO2 column densities by TROPOMI from September 2015 to September 2019. Grey rectangles in lower panel (b) indicate the monsoon period in Pakistan.
Fig. 7. Correlation plots between (a) daily MAX-DOAS NO2 and OMI NO2 (b) daily MAX-DOAS NO2 and TROPOMI NO2 (c) average monthly MAX-DOAS NO2 (average of 1300–1400 h, UTC+5) and OMI NO2 in molecules cm–2 from September 2015 to September 2019, and (d) average monthly MAX-DOAS NO2 and TROPOMI NO2 from September 2018 to September 2019.
4.1 Error Analysis
5 CONCLUSIONS
ACKNOWLEDGMENTS
REFERENCES