El uso de datos satelitales en combinación con mediciones realizadas en superficie puede proporcionar información valiosa acerca de la química atmosférica y la calidad del aire. En este estudio se comparan mediciones en superficie de dióxido de nitrógeno (NO2) realizadas mediante la técnica de espectroscopia óptica de absorción diferencial (DOAS, por sus siglas en inglés) con mediciones del instrumento satelital para la medición de ozono (OMI, por sus siglas en inglés) realizadas de 2006 a 2011. Las mediciones realizadas desde la superficie presentaron grandes variaciones diarias y fueron en promedio tres veces más altas que las columnas medidas desde el espacio. La diferencia se atribuye a una fuerte heterogeneidad horizontal presente en la capa inferior de las columnas de NO2, las cuales fueron muestreadas por el instrumento satelital a partir de un área extensa; de igual manera, esta discrepancia se atribuye a la sensibilidad reducida del satélite cerca de la superficie, donde se encuentran las mayores concentraciones. A partir de los datos del OMI analizados se reconstruyeron mapas de distribución de NO2 sobre el centro de México, y se identificaron tres áreas principales de interés: la zona metropolitana de la Ciudad de México, que fue el área predominante; la zona altamente industrializada de Tula, al norte, y el valle de Cuernavaca, al sur. En este análisis se detectaron de igual forma variaciones estacionales de columnas de NO2 sobre el centro de México: se encontraron columnas más altas durante la estación fría y seca, seguidas por las de la estación caliente y seca; las columnas más bajas se encontraron durante la época de lluvias. Este conjunto de datos evidencia el transporte de contaminación de este gas desde Tula hasta la Ciudad de México, así como al Valle de Cuernavaca.
The use of satellite data in combination with ground-based measurements can provide valuable information about atmospheric chemistry and air quality. In this study, ground-based Differential Optical Absorption Spectroscopy (DOAS) measurements of nitrogen dioxide (NO2) conducted in central Mexico are compared with the space-borne Ozone Monitoring Instrument (OMI) dataset of 2006-2011. Ground-based measurements exhibited large day-to-day variations and were on average three times higher than the space-borne derived average over the observation site. This difference is attributed to strong horizontal inhomogeneity of the lower layer of the measured NO2 columns, sampled over a large footprint from the satellite instrument. Also, a reduced sensitivity of the satellite observation near the surface, where the largest concentrations are expected, could be responsible for this large discrepancy. From the analyzed OMI dataset, distribution maps of NO2 above central Mexico were reconstructed, allowing to identify three main areas with increased NO2 column densities: The dominating metropolitan area of Mexico City, the heavily industrialized region of Tula to the north and the Cuernavaca valley to the south. In this analysis, seasonal variability of NO2 columns over central Mexico was detected, finding higher NO2 columns during the dry and cold season, followed by the dry and warm period, and finally the lowest NO2 columns were found during the rainy season. Pollution transport of this gas from Tula into Mexico City, as well as towards the Cuernavaca valley, is evident from this dataset.
Mexico City is the third largest city of the world, with more than 20 million inhabitants (UN, 2012). It is located at 19.4° N latitude and approximately at an elevation of 2200 m above sea level (masl). For several years Mexico City has experienced severe air quality problems, which have been addressed by authorities trough the implementation of various programs. While the countermeasures and technological improvements have effectively decreased the amount of some pollutants, the city continues to experience pollution episodes resulting on the frequent violation of air quality standards.
Nitrogen dioxide (NO2) plays a major role in tropospheric and stratospheric chemistry (Crutzen, 1979). In polluted regions such as Mexico City, tropospheric NO2 concentrations can be highly variable in time and space and are influenced by both natural and anthropogenic emissions. In Mexico City, most of the NO2 present in the atmosphere comes from fossil fuel combustion from transportation and partly from industrial activities. NO2 levels in Mexico City have decreased by 30% between 1990 and 2008, and the average surface concentration registered in 2008 by four stations around the Universidad Nacional Autónoma de México (UNAM) campus in Mexico City—where the ground-based measurements presented in this study were conducted—was 56.8 ppb (SMA-GDF, 2010).
Satellites provide global observations, allowing to retrieve information from nearly every region in the globe (Rees, 2001). There are a number of instruments measuring from space, covering a wide range of wavelengths, with different spatial resolution and viewing geometries. In this work we have focused on the Ozone Monitoring Instrument (OMI) (Levelt et al. 2006a, b), a UV-VIS spectrometer on board NASA’s Earth Observing System-Aura satellite, able to provide a daily global coverage.
Data from the OMI instrument have been used in an extensive number of studies involving different atmospheric gases as NO2, O3, SO2, HCHO, among others (Bhartia, 2002; Chance, 2002). The NO2 data products provide information about total and tropospheric NO2 columns. The NO2 dataset in particular, has been validated in different studies. Celarier et al. (2006) compared NO2 columns from OMI with NO2 columns measured by the Systéme d’Analyse par Observation Zénithale (SAOZ) network, finding good agreement between the two of them although the measurements were not done at the same time. As a continuation, the different measurements and validation activities for OMI NO2 stratospheric, tropospheric and total columns using ground- and aircraft-based measurements were summarized by Celarier et al. (2008). Good agreement was reported between OMI and nearby instruments, finding correlations higher than 0.6 with ground-based instruments, whereas OMI NO2 stratospheric and total columns were found to be underestimated by 14 and 15-30%, respectively.
During the Dutch Aerosol and Nitrogen Dioxide Experiments for Validation of OMI and SCIA-MACHY (DANDELIONS) 2005 and 2006 field experiments in the Netherlands, Brinksma et al. (2008) found good agreement between tropospheric from OMI and MAX-DOAS (Multi-AXis), as well as between total NO2 from OMI and direct sun observations. Boersma et al. (2008) compared OMI NO2 tropospheric columns with in situ aircraft measurements during the INTEX-B campaign, the results showed good correlation with no significant bias between the two data sets. Irie et al. (2008) used MAX-DOAS measurements to validate OMI NO2 columns during the Mount Tai Experiment 2006 finding that OMI data may have a positive bias of 20% over the North China Plain. Gruzdev and Elokhov (2010) found that tropospheric OMI NO2 columns measured over Zvenigorod, Russia were approximately 40% lower than ground-based twilight measurements. Vlemmix et al. (2010) compared tropospheric ground-based measured NO2 columns with OMI-satellite tropospheric data finding no significant difference and a correlation of 0.88.Ionov et al. (2008) conducted a geophysical validation of OMI NO2 vertical column measurements with collocated ground-based observations at two stations in central Asia and Europe finding that mid-latitude OMI stratospheric NO2 column data underestimates ground-based measurements.
Direct-Sun mode measurements have been conducted as well in order to validate OMI NO2 columns. Wenig et al. (2008) used direct-Sun mode Brewer measurements to conduct a validation study with OMI tropospheric NO2 columns concluding that monthly averages of coinciding measurements yielded a correlation of 0.9, however OMI data was 25% lower than the ground-based Brewer measurements. Herman et al. (2009) conducted intercomparisons of NO2 columns measured with ground-based direct solar irradiance measurements to OMI NO2 columns finding good agreement (correlation of 0.73).
Regional and global models have also been used to compare and validate OMI NO2 tropospheric columns (Kim et al, 2009; Herron-Thorpe et al., 2010; Huijnen et al, 2010; Dirksen et al, 2011) and to estimate NOx emissions from different sources (Lin, 2012). Recently, Beirle et al. (2011) studied NOx emissions and lifetimes from megacities and power plants, while Valin et al. (2011) used the super-zoom mode of OMI to observe NO2slant columns from point sources and two cities, capturing significant spatial variability at fine scales.
In this work, the OMI NO2 data product (OMI/ Aura nitrogen dioxide total and tropospheric column 1-orbit L2 swath 13×24 km, data set version 003) has been analyzed and compared to ground-based measurements conducted at UNAM campus in Mexico City. NO2 column maps were reconstructed from the OMI data product using an algorithm based on the optimal estimation approach with a weighted assignment of the measurement onto a fine grid, in order to identify the spatial features of its distribution with higher resolution. Focus has been given to the study of NO2 mainly to its role as a catalyst in the production of tropospheric ozone (O3), which is responsible for one of the most severe air quality problems in the Mexico City metropolitan area (MCMA).
2Data set2.1OMIOMI is a nadir viewing spectrometer covering a spectral region of 264-504 nm, conducting atmospheric chemistry measurements with a resolution of 0.42-0.63 nm since 2004 on board NASA’s Earth Observation System (EOS)-Aura satellite. Aura follows a sun-synchronous polar orbit (705 km altitude), with an ascending local equator crossing time at 13:45. OMIs observations provide complete global coverage in one day with a nominal ground footprint of 13×24 km2 at nadir (OMI Team, 2009). It measures a number of air quality components such as NO2, SO2, BrO, HCHO, and aerosols (Levelt et al. 2006a, b).
The OMI operational algorithm is described in detail by Bucsela et al. (2006). For the evaluation of spectra recorded by OMI, a similar DOAS method (Platt and Stutz, 2008) is used as for the ground instruments. In this method a spectral fit is applied to all measured spectra. The evaluation of spectra results on slant column densities (SCDs) for each viewing pixel. The computed SCDs are then converted to vertical column densities (VCDs) using an air mass factor (AMF). The AMF is the ratio between the retrieved SCD and the VCD. The AMF concept was introduced by Solomon et al. (1987) and is dependent on clouds and aerosols, as well as viewing geometries and surface albedo.
In this work, the OMI/Aura nitrogen dioxide total and tropospheric column 1-orbit L2 swath 13×24 km, data set version 003, level 2, algorithm version 1.1.4.4 data product was downloaded from NASA’s Goddard Earth Sciences Data and Information Services Center (GESDISC) using the MIRADOR 1.42 tool (http://disc.sci.gsfc.nasa.gov/Aura/data-holdings/OMI). We have analyzed the dataset and only considered total NO2 columns from cloud fraction below 20% over the period 2006-2011. This dataset has been reprocessed and released on December 2011 (Claas, 2012).
2.1.1Reconstruction of the mean column NO2 distributionThe mean column NO2 distribution over a fine grid (of the order of about 2 km) can be reconstructed from measurements done on a larger footprint (here 13×24 km2 for OMI), solving the mathematically ill-posed problem Y = Kx + ε. In this equation, Y represents the vector containing all OMI NO2 measurements, x is the solution vector which contains the estimation of the mean column NO2 distribution on the chosen grid, ε represents the deviation between a column measured/retrieved from the space measurement and the estimated mean value. ε is mainly dependent on the day-to-day variance, but also on the seasonal variance and measurement precision. The matrix K describes how the measurements are related with the mean gas horizontal distribution. It represents a forward model that reproduces mainly the footprint of the OMI instrument (13×24 km2 at nadir). The reconstruction uses a Tihkonov-smoothing-constraint, similar as the one used for profile retrievals or in image restoration. Further details of this method can be found in Stremme et al. (2013).
2.2Zenith scattered sunlight DOASGround-based DOAS measurements were conducted at the UNAM campus in Mexico City (19.33° N, 99.18° W) between August 2010 and December 2011. Scattered light was collected by means of a telescope and transferred into an Ocean Optics TR2000 spectrometer with a spectral resolution of ~0.6 nm and wavelength coverage of 291-501 nm. Spectra were taken every minute using the DOASIS software (Kraus, 2006) and were evaluated using the QDOAS (version 1.00) software (Fayt et al., 2011). Differential cross-sections of NO2 at 298 K (Vandaele et al., 1998), O3 at 221 K and 241 K (Burrows et al, 1999), the oxygen dimer O4 (Hermans et al, 1999) and a Ring spectrum, generated at 273 K from a high resolution Kurucz file using the QDOAS software (Fayt et al., 2011) were included in the analysis. The oxygen dimer O4 was included in the evaluation procedure in order to detect enhancement on the light path of the measurements since it is expected to have a constant concentration in the atmosphere. Evaluations were conducted in the 405 and 465 nm wavelength ranges with a fixed reference spectrum measured on 21 December 2011 at 14:21 local time (LT) under clear sky and high ventilation conditions for the whole time series. Geometrical AMFs were calculated and used to convert SCDs to VCDs for further comparison with OMI results. For this it was considered that the light path through a trace gas layer varies with 1/cos θ, where θ is the solar zenith angle (Chen et al, 2009). It should be noted that this simplified approach without the use of a full radiative transfer calculation that effectively incorporates scattering processes, reflection on the earth’s surface, refraction, the curvature of the earth, the vertical distribution of trace gases and the influence of aerosols leads to additional errors discussed in detail by Platt and Stutz (2008).
2.3Surface measurementsSince 1986, a network of 34 stations run by the local government (Red Automática de Monitoreo Atmosférico, RAMA) reports continuous and permanent measurements of concentrations of O3, sulfur dioxide (SO2), nitrogen oxides (NOx), carbon monoxide (CO) and particulate matter PM10 and PM2.5 to evaluate air quality in the metropolitan area. All data generated by RAMA are freely available online at http://www.calidadaire.df.gob.mx/. NO2 is determined using a nitrogen oxides analyzer based on the chemilumi-nescence measurement technique. For this study NO2 concentration results generated during 2010 and 2011 by Pedregal (19.32° N, 99.20° W), Tlalpan (19.25 N, 99.18 W), Santa Úrsula (19.31°N, 99.14° W), Merced (19.42° N, 99.11° W) and FES Acatlán (19.48° N, 99.24° W) stations situated close to the study site were used to calculate NO2 columns, using mixing layer heights (section 2.4) as well as pressure and temperature data (section 2.5).
2.4LIDARA Vaisala CL31 ceilometer, an active instrument based on the LIDAR (light detection and ranging) measurement technique, was used to measure back-scattered aerosol density at the UNAM campus, allowing the estimation of mixing layer heights. The measurement principle involves the transmission of pulsed laser radiation into the atmosphere, which is further backscattered and detected with a certain time delay. The measurement range of the instrument used in this field study was from 0 up to 7500 m above the ground with a 10 m resolution. The LIDAR wavelength of 905 nm was provided by an InGaAs MOCVD laser diode, 110 ns, 1.2 μJ per pulse, and repetition rate of 8192 Hz. The mixing height was retrieved using the CL31 MLH software based on the gradient method, which is focused on selecting the maximum of the negative gradient of the backscatter coefficient (β) to be the top of the mixed layer (-dβ/dx). A more detailed description of this methodology can be found in Münkel et al. (2007).
2.5Meteorological parametersMeteorological parameters were recorded by a Davis Vantage Pro 2 station installed at the UNAM campus, forming part of the Programa de Estaciones Meteorológicas del Bachillerato Universitario (PEMBU) network. From the dataset recorded by the meteorological station, pressure and temperature were used for the calculation of ground-based NO2 columns based on NO2 concentrations measured by the RAMA network.
3Results and discussion3.1Ground-based NO2 columnsFigure 1 shows differential NO2columns measured at UNAM between August 2010 and December 2011. The gaps in the graph represent periods where no measurements were conducted, either due to experimental problems, or periods of cloudiness or rain identified through the analysis of the data itself using evaluation results of the oxygen dimer. Only measurements conducted between 14:00 and 15:00 LT are plotted in order to depict columns coincident with OMI overpasses (~14:21 LT) over Mexico City. Large variability on day-to-day NO2columns was measured, illustrating the dynamics of fresh emissions and rapid conversion of pollutants occurring in the study area. The ground-based average NO2 column values and standard deviation measured from August 2010 to December 2011 between 14:00 and 15:00 LT was (5.61 ± 3.32) ×1016 molecules/cm2.
Differential NO2 columns measured at UNAM between (a) August-December 2010, (b) February-June 2011, and (c) August-December 2011. The vertical gray lines represent the columns evaluation error. Only measurements conducted between 14:00 and 15:00 LT are shown. The data points showed in this figure are zenith measurements.
The average diurnal behavior of differential NO2 columns measured at UNAM on the same measurement period are presented in Figure 2. The maximum of the measured NO2 differential columns is found around noon, which is consistent with the reported hourly NOx emissions for the MCMA by the official inventory (SMA-GDF, 2010).
3.2NO2 spatial distribution
The average distribution of NO2 total columns derived from OMI measurements conducted between June 2006 and December 2011 over the Mexico basin is depicted in Figure 3. Higher NO2 total columns are observed in the downtown area of Mexico City, which gradually decrease with distance from the city center. The average and standard deviation of all OMI NO2 columns measured over UNAM during 2010-2011 at its overpass time, ~14:21, is (1.18 ± 1.1)×1016 molecules/cm2. While Figure 3 is simply the average of the columns falling on each grid point, Figure 4 shows a reconstruction of the mean column gas distribution according to the methodology described in section 2.1.1. The reconstructed value of the average over UNAM for the same time period is 1.7×1016 molecules/cm2. A clear difference can be observed between the average and reconstructed column gas distributions with the average column being smaller over UNAM. This bias is present during all seasons, it might originate from the sensitivity of the space-borne measurements to pollution near the surface and is discussed in more detail in section 3.3.
NO2 total columns distribution reconstructed from OMI measurements (see section 2.2.1) conducted between June 2006 and December 2011. Same scale as Figure 3.
From Figures 3 and 4, three main areas of interest in central Mexico can be identified: the Tula industrial area to the north, a dominating blur centered around Mexico City, and the Cuernavaca valley to the south. However, the reconstructed map in Figure 4 allows to better appreciate the finely defined areas where either NO2 columns are enhanced or diminished over the studied area. Lacustrian and mountainous areas inside the MCMA are, for instance, more clearly observed in Figure 4 than in Figure 3. Also, the industrial and more densely populated areas are better defined as well.
Unique distribution structures of NO2 columns above the MCMA were identified through the reconstruction of mean NO2 column distribution in Figure 5a. NO2 column enhancements can be observed on the northern and central parts of the city, which are characterized by intense and heavy traffic (central and northern parts, respectively) and industrial activities (northern part). On the southern, eastern and western sides, the NO2 columns tend to decrease with distance of the city center, while on the northern part, the enhancement of NO2 columns is even extended towards the Tula industrial area (Fig. 5b), from which a clear footprint was observed and differentiated from the MCMA.
NO2 total columns distribution over different areas of central Mexico reconstructed from OMI measurements conducted between June 2006 and December 2011. (a) The Mexico City metropolitan area; (b) the heavily industrialized region around Tula, Hidalgo; and (c) the Cuernavaca valley. Note that different scales have been used in the plots.
The region around Tula shown in Figure 5b, is located in the state of Hidalgo, ~50 km northwest of Mexico City. It is characterized by the presence of a power plant, a refinery, several cement works and other manufacturing industries. The OMI footprint clearly shows higher NO2 column values close to the industries. As in the case of the MCMA, the NO2 column enhancement gradually decreases with distance of the pollution sources with an increase of NO2 columns towards Mexico City. From Figure 5b, an empiric dispersion pattern of the plume of the Tula industrial complex towards the southeast can be inferred. This result is consistent with previous studies using dispersion models addressing the transport of the emissions from the Tula industrial complex towards the MCMA (CCA-UNAM, 2008; de Foy et al., 2009).
The third identified area in central Mexico where a clear, though less intense NO2 column enhancement is observed, is the Cuernavaca valley shown in Figure 5c. This result was somewhat unexpected when compared to the relatively small column enhancements observed above the cities of Puebla and Toluca, which are similar to Cuernavaca in terms of population and industrial activity. Cuernavaca valley, however, is much lower in altitude (< 1500 masl), host of a large number of sugar cane plantations. During the industrial processing and harvesting intense biomass burning is practiced, which could explain the increased NO2 column values observed with the OMI instrument. Another possible explanation is the impact of emissions by the surrounding states. Recently Salcedo et al. (2012) reported an efficient transport of ozone and its precursors to this area by the surrounding Mexico City and Puebla state, while Melamed et al. (2009) detected during a study in 2006, NO2 transport events southeast of Mexico City from and towards the Cuautla region.
3.2.2Seasonal variabilityAn additional approach was used analyzing OMI data into three different seasons: (1) dry and cold covering the months of November, December, January and February; (2) dry and warm for the months of March, April and May; and (3) the rainy season including the months of June, July, August, September and October. Figure 6 depicts NO2 total column distributions reconstructed from OMI measurements conducted between June 2006 and December 2011 for the dry and cold (Fig. 6a), dry and warm (Fig. 6b) and rainy seasons (Fig. 6c). The same scale is used in all three plots for comparison. Higher NO2 columns were measured during the dry and cold season, followed by the dry and warm period and finally the lowest NO2 columns were found during the raining season. For the dry-cold and dry-warm seasons the pattern is confirmed by the ground-based measurements at UNAM (Fig. 7). Wallace and Kanaroglou (2009) have reported sensitivity of OMI to subtle temperature changes, while Boersma et al. (2009) observed NO2 columns 1.6-2.7 times higher in winter than in summer in Israeli cities. The atmospheric stability characteristic of the dry and cold season promotes a larger number of molecules to be present in the studied area, while during the raining season, water is mainly responsible for the removal of pollutants in the atmosphere.
NO2 total columns distribution over central Mexico reconstructed from OMI measurements conducted between 2006-2011 on: (a) November, December. January and February (dry and cold); (b) March, April and May (dry and warm); (c) June, July, August, September and October (rainy season). Scale from 0 to 1.8×1016 molecules/cm2.
The NO2 total column distribution reconstructed during the rainy season (Fig. 6c) allows to better identify the presence of areas with NO2 column enhancements over the cities of Toluca and Puebla, whose contribution is somewhat lost by the large column densities measured over the MCMA, the Tula industrial area or even the Cuernavaca valley during other months of the year. Also, a transport between the MCMA and Cuernavaca is visible during the dry seasons both through the Cuautla and Cuernavaca sides, where the highway connecting both cities runs. This confirms the pollution transport highlighted in the previous section and addressed by previous studies (Melamed et. al, 2009; Salcedo et al, 2012).
3.3Comparison of ground-based and space-borne measurementsA comparison of direct coincident measurements was done for NO2 ground-based and OMI measured columns. All ground-based measurements conducted during the overpass time of OMI over UNAM (~14:21 LT) were used with a tolerance of ± 20 min. OMI measurements having a pixel center within a certain radius are used. The number of coincident days for the direct comparison is increasing when using larger radius. Ground-based and OMI measurements conducted on the coincident days and times were averaged and are reported in Table I along with its standard deviation, correlation coefficient, linear-fit offset and linear-fit slope. In general, ground-based measured NO2 columns were 2.7-6 times higher than OMI measured columns, which could be explained by the spatial discrepancy between the ground-based measurements (basically local) and the satellite averaged measurements. The highest correlation coefficient (0.807) was found when comparing OMI data with a pixel center falling within a 20 km radius around UNAM; however, good correlation was also observed for data with pixel centers of 15 to 50 km radius around the ground-based site.
Comparison of NO2 ground-based columns measured at UNAM and NO2 columns measured by OMI at overpass time 14:21 LT ± 20 min.
| Radius (km) | Number of coincident days | UNAM average column (molecules/cm2) | UNAM stdev column (molecules/cm2) | OMI average column (molecules/cm2) | OMI stdev column (molecules/cm2) | Correlation coefficient | Linear fit offset | Linear fit slope |
|---|---|---|---|---|---|---|---|---|
| 9 | 4 | 2.48E + 016 | 9.86E + 015 | 7.20E + 015 | 2.96E + 015 | −0.059 | 7.64E + 015 | −0.018 |
| 10 | 5 | 2.91E + 016 | 1.28E + 016 | 8.12E + 015 | 3.31E + 015 | 0.439 | 4.83E + 015 | 0.113 |
| 15 | 6 | 3.52E + 016 | 1.87E + 016 | 1.27E + 016 | 1.22E + 016 | 0.781 | 5.16E + 015 | 0.509 |
| 20 | 6 | 3.52E + 016 | 1.87E + 016 | 1.11E + 016 | 1.00E + 016 | 0.807 | 4.12E + 015 | 0.433 |
| 25 | 6 | 3.52E + 016 | 1.87E + 016 | 1.05E + 016 | 9.24E + 015 | 0.772 | 2.88E + 015 | 0.381 |
| 30 | 6 | 3.52E + 016 | 1.87E + 016 | 1.01E + 016 | 8.30E + 015 | 0.764 | 1.86E + 015 | 0.339 |
| 40 | 6 | 3.52E + 016 | 1.87E + 016 | 8.70E + 015 | 5.88E + 015 | 0.734 | 5.87E + 014 | 0.231 |
| 50 | 6 | 3.52E + 016 | 1.87E + 016 | 7.86E + 015 | 4.05E + 015 | 0.703 | 2.51E + 015 | 0.152 |
| 60 | 7 | 3.49E + 016 | 1.71E + 016 | 7.21E + 015 | 3.24E + 015 | 0.687 | 2.66E + 015 | 0.130 |
| 70 | 7 | 3.49E + 016 | 1.71E + 016 | 6.89E + 015 | 2.76E + 015 | 0.583 | 3.61E + 015 | 0.094 |
| 80 | 7 | 3.49E + 016 | 1.71E + 016 | 6.67E + 015 | 2.28E + 015 | 0.512 | 4.29E + 015 | 0.068 |
| 90 | 7 | 3.49E + 016 | 1.71E + 016 | 6.50E + 015 | 2.10E + 015 | 0.422 | 4.70E + 015 | 0.052 |
| 100 | 8 | 3.46E + 016 | 1.58E + 016 | 5.78E + 015 | 2.17E + 015 | 0.334 | 4.19E + 015 | 0.046 |
The averages from OMI NO2 columns shown in Table I, however, are underestimated since the instrument’s footprint of 13×24 km2 at nadir covers a large area and the UNAM site is within a region with large gradients in its horizontal distribution. A better representation of the fine structure is obtained by the reconstruction shown in Figure 4 according to the methodology described in section 2.1.1. The reconstructed value of the average over UNAM for the entire time period is 1.7×1016 molecules/cm2. This value is approximately three times lower than the average ground-based NO2 column measured over UNAM.
Furthermore, the satellite observation might not be sensitive enough close to the surface, where most of the pollution is present. A sensitivity of 0.5 would actually explain the difference by a factor of 2 between the mean values of the ground-based and space-borne DOAS measurements. The surface albedo of Mexico City in the 450 nm spectral region was measured by Castro et al. (2001) and recently during the MILAGRO field experiment by Coddington et al. (2008) and found to be in average around 0.05 and 0.07 respectively. Therefore, the estimated sensitivity of space-borne DOAS retrievals on surface near NO2 anomalies (local pollution in the mixing layer) should range between 0.2 and 0.7, which correspond to the two examples for surface albedo, 0.02 and 0.15, showed by Eskes and Boersma (2003).
Monthly average comparisons were also made between ground-based measurements, OMI and RAMA calculated NO2 columns (Fig. 7). Again, only ground-based measurements within ± 20 min of the OMI overpass time (~14:21 LT) above UNAM are kept. OMI averages were restricted to columns with pixel centers within ±0.1 degrees latitude/ longitude of UNAM. Ground-based reported RAMA NO2 columns were calculated from the average concentrations measured by five stations located near UNAM, multiplied by the measured mixing layer height (section 2.4) and correcting for pressure and temperature (section 2.5). Amonthly cycle was clearly identified, however, it would seem to differ slightly depending on the origin of the data. The calculated columns from surface RAMA data were the highest, followed by the ground-based DOAS measurements and finally the OMI data yielded the lowest values. This result was expected since the ground-based RAMA stations sample NO2 at a single sampling point very close to the sources and/ or when it has just been recently converted. The DOAS ground-based measurements are conducted close to the sources as well, however, its location and vegetation around the UNAM campus could partly dampen the detection of NO2 at the current sampling point. Furthermore, while the ground-based DOAS instrument is sampling at a single point, the observed columns by OMI are effectively averaged over the instrument’s field of view (13×24 km2 at nadir) yielding lower NO2 columns than the other two instruments probably because of the averaging of polluted and unpolluted areas.
3.4Spatial distribution of NO2 columns reconstructed from UNAM measurementsAn additional analysis was made in the interest of gaining more information regarding the possible origin of the NO2 columns measured over UNAM, especially due to the inhomogeneous distribution of NO2 columns in the city, attested by OMI measurements. The averaged surface wind speed and wind direction measured at five stations near UNAM were used to compute the statistical wind propagation pattern that indicates the probability that a column measured at UNAM is originated or transported from a certain area of the city. Wind directions were analyzed on 30° intervals, and a probability was assigned to each sector according to the relative frequencies. Each sector was differentiated into three regions according to the intervals of certain probability (15.85,68.3, and 15.85%) related to different intervals of confidence of a log-normal distribution. Accordingly, wind speeds were expected to follow a log-normal distribution as well, described by the median and the log-normal standard deviation (Fig. 8a).
The false color of each region in Figure 8a indicates the probability that an airmass will be measured at UNAM with the certain time delay of 2.25 h. This time delay used for the visualization (Fig. 8) is the estimated lifetime of NO2 in the mixing layer with a typical temperature of the Mexico City boundary layer and according to Dils (2008). The measured NO2 columns were then distributed to the corresponding sectors using wind direction and wind speed at the time when measurements were conducted. The distribution pattern of the averaged NO2 total columns reconstructed from the measurements at UNAM and wind measurements is presented in Figure 8b. The result of this approach indicates a tendency to measure higher NO2 columns when the wind direction is predominant from the north-northwest, coinciding with the most industrialized and populated area of the city.
4ConclusionsNO2 columns were measured at UNAM in Mexico City and compared to space-borne observations conducted by the OMI instrument. Space-borne NO2 columns were lower than ground-based measurements, which presented large day-to-day variations. The difference between space-borne and ground-based columns could be explained by the strong horizontal inhomogeneity of the lower layer of the measured NO2 columns and by the difference in horizontal sensitivity of the two techniques (local versus averaging over 13×24 km2). More instruments located in and around the MCMA would allow to better evaluate the spatial distribution obtained from satellite measurements. A network of several MAX-DOAS has already been installed starting the fall of 2012, which will result in better ground-based coverage and accuracy of the reported columns.
A reconstructed distribution of weighted NO2 columns over a fine grid was made over a large region in central Mexico, allowing to successfully define three main areas with enhanced pollution: the dominating metropolitan area of Mexico City, the heavily industrialized Tula region and the Cuernavaca valley. Other mayor urban areas such as Toluca and Puebla showed in average relatively lower column densities than Cuernavaca. The considerably higher altitude of Toluca probably favors the photolytic destruction of NO2, and this area does not appear to receive as much pollution from Mexico City, as the Cuernavaca valley. The importance of considering changes in albedo, aerosols load and their refractions indexes in NO2 photolysis rates should also be noted. The unanticipated enhancement of NO2 columns above the Cuernavaca valley, corroborates results from previous studies conducted by Melamed et al. (2009) and Salcedo et al. (2012), where a possible transport of pollutants from the MCMA have been discussed. A plausible source arising from the biomass burning activities in that region is suggested for the high NO2 columns detected.
From the reconstructed maps it is also possible to observe an important interconnection between the MCMA and the Tula industrial region, where emissions from one can evidently affect the other and the transport of pollutants between them most probably occurs. In addition, distinct seasonal variability of NO2 columns over central Mexico was found, concurring with previous studies, with higher NO2 columns measured during the dry and cold season, followed by the dry and warm period and finally the lowest NO2 columns were found during the rainy season. Furthermore, the spatial distribution of NO2 columns reconstructed from measurements conducted at UNAM indicate a tendency to measure higher NO2 columns when the wind direction is predominant from the north-northwest, coinciding with the most industrialized area of the city.
The authors gratefully acknowledge financial support provided by DGAPA-UNAM through project number IN119310 and a DGAPA-UNAM postdoctoral research scholarship. The OMI project is managed by the Netherlands Agency for Aerospace Programs (NIVR) and the Royal Netherlands Meteorological Institute (KNMI). The authors acknowledge the free use of NO2 total column data from the OMI sensor accessed via the NASA Aura Validation Data Center (AVDC) and Goddard Earth Science (GES) Data and Information Services Center (DISC). Caroline Fayt and Michel van Roozendael from the Belgian Institute for Space Aeronomie are acknowledged for the free use of the QDOAS software used for the evaluation of the ground-based data presented in this work. SIMAT and the SMA-GDF are thanked for providing the in situ measurements and the Pro-grama de Estaciones Meteorológicas del Bachillerato Universitario (PEMBU) network is acknowledged for the meteorological data. Alfredo Rodríguez from the Centro de Ciencias de la Atmósfera workshop is acknowledged for providing technical support on the installation of the instrumentation used to conduct the ground-based measurements.