Figure 7 summarizes the application of the k-means clustering method (with k = 6) to the seasonal subsets,
Hwinter,
Hspring,
Hsummer, and
Hautumn, of the data objects (events) of wind circulation and pollutant-flows. This figure shows an arrangement with four columns (one per seasonal period) and five rows (one for wind circulation, and four for the pollutant-flows of NO
2, O
3, PM
10, and SO
2). Each graph shows six plots (one per cluster) differentiated by the line color. The algorithm enumerates the clusters according to their sizes, from the larger to the smaller, so colors mean no more. Each plot presents the
hourly population of one cluster, that is to say, the number of data objects belonging to the cluster as a function of the hour of the day (Mexico City local time, UTC-6 h).
For wind circulation or any given pollutant-flow, a qualitative comparison of graphs in the respective row allowed recognizing some seasonal similarities and differences between the clusters. For example, for the wind circulation (first row), if we pay attention to the clusters represented by a blue-line plot in the four seasons, we can observe an evident regularity: the hourly populations have similar trends and indicate that the majority of their events occurred during the night, mainly from midnight until dawn. We interpreted this regularity as the possible existence of a wind pattern. The plots of the pollutant-flows also show a similar regularity.
Eight of the sixteen parameters of the 4-cell model are the components of the flow vectors in the quadrants (NW, NE, SW, and SE) of the MCMA. This complementary information and the hourly population plots allowed recognizing the wind circulation patterns shown in
Figure 8, and the pollution-flow patterns shown in
Figure 9,
Figure 10,
Figure 11 and
Figure 12. We found seven patterns for wind circulation and nine patterns for each of the pollutant-flows (NO
2, O
3, PM
10, and SO
2). In
Figure 8, we show the hourly population of each wind pattern and a graphic representation of the wind velocity at the quadrants.
Figure 9,
Figure 10,
Figure 11 and
Figure 12 show the hourly population of each pattern, a graphic representation of the flow vector at the quadrants, the surface distribution of the pollutant concentration, and the mean spatial pollutant concentration (expressed in ppb for NO
2, O
3, and SO
2, and in µgm
−3 for PM
10). In these figures, the hourly population of a pattern is an average over the seasons the pattern occurred; the vectors (the wind velocity and the pollutant-flow vectors) are averages over the elements of the pattern and the seasons; the spatial distribution of the pollutant concentrations is the average over the pattern elements and seasons; and the mean spatial concentration of a pollutant is the average over the positions of the average pollutant distribution. The lengths of the arrows that represent the wind and flow vectors were scaled to get the larger one fitted in the square that represents the quadrant, allowing a qualitative comparison of the vector magnitudes among the patterns of wind circulation or the patterns of any of the pollutant flows.
3.1. Wind Circulation Patterns
The proposed methodology for examining the database of the wind events that occurred in Mexico City from 2001 to 2010 allowed discovering there the presence of seven wind circulation patterns and the estimation of their seasonal and annual occurrence frequencies, which we summarized in
Table 2. We denoted these patterns as WIND-Pn, with n = 1, 2,…, 7, and enumerated them trying to follow the order of their appearance during the day. In
Figure 8, we presented the hourly population of each wind pattern, including a sketch out of the corresponding wind velocity in the city quadrants.
In the following paragraphs, we provide brief descriptions of the wind circulation patterns. We used the Mexico City local time (UTC-6 h) in all figures and descriptions.
WIND-P1: early morning downslope winds. WIND-P1 was the wind pattern with the highest annual frequency (30%). It was observed systematically throughout the year during the night, showing a large peak around hour 6 (
Figure 8) and a small peak around hour 15. The mean wind velocities at the quadrants suggest downslope winds from the surrounding mountains, converging towards Mexico City basin. The converging winds represented by the small peak reflect the presence of convective upwards winds due to the urban heat island (UHI) in the Mexico City region. Jauregui [
7] and Klaus et al. [
33] pointed out that the Mexico City downslope winds are reinforced by UHI.
WIND-P2: northeasterly and easterly winds. The annual average frequency of this wind pattern was 21%. It is observed systematically throughout the year, especially during winter (25%) and summer (24%). The mean velocities correspond to winds blowing from NE in almost all the quadrants of the city (trade winds). The wind events of this pattern occurred from the sunrise (hour 7) up to the hour 17. Its hourly population presents a high peak between hours 10 and 11.
WIND-P3: midday northerly winds. This pattern had an annual frequency of 9%. It was observed during spring (12%), summer (9%), and autumn (14%) from hour 7 to hour 21, with a peak around hour 13. The mean wind velocities in the quadrants indicate winds blowing from the north and northeast sectors.
WIND-P4: westerly and southerly winds. This wind pattern was observed only winter season with a seasonal frequency of 20%; its events occurred throughout the day, but mainly from noon to midnight, reaching its maximum population around hour 17. The mean wind velocities represent winds with a westerly component in the west quadrants and southerly winds in the east quadrants. These velocities correspond to winds blowing from the south and southeast city sectors through a gap between the mountains of the Sierra Ajusco-Chichinautzin and the Sierra Nevada, following the ventilation channel S–N located at the east side of the city (
Figure 1). The local winds of this pattern could be related to the subtropical jet stream of winter [
34] or with the westerlies that are permanently occurring in subtropical and middle latitudes, but coupled with the winds driven by the gap between the mountains located at the southeast. Doran and Zhong [
35] described the main characteristics of a gap wind system in the southeastern corner of Mexico City that produces low-level jets occurring regularly during the late winter.
WIND-P5: afternoon northerly winds. This pattern was observed throughout the year with an annual frequency of 12%. It occurred from noon to midnight, with its maximum around hour 19. The mean wind velocities at the quadrants indicate winds blowing from the north sectors.
WIND-P6: midnight downslope winds. This pattern was observed systematically throughout the year, mainly during the nighttime hours, from the sunset (hour 18) to the sunrise (hour 6 of the next day), with its maximum around midnight. It got an annual average frequency of 18%. The population plot of this pattern has a small peak close to hour 15, which could be related to the UHI effect of the city. The wind velocities in the western quadrants of the city represent downslope flows from Sierra de las Cruces (at west) and Sierra del Ajusco-Chichinautzin (located at SW and S) toward the town. In the eastern quadrants, otherwise, late afternoon northerly gap winds are observed guided by the S–N ventilation channel of the city.
WIND-P7: UHI-driven winds. This wind pattern was observed only during the spring (9%) and autumn (10%) seasons, with a low annual occurrence frequency of 5%. The wind events of this pattern occurred from the hour 14 to the hour 22, and its population reached a maximum value close to the hour 17. The mean wind velocities of the city quadrants indicate cyclonic winds converging towards Mexico City downtown, which seems to be related to an upwards flow driven by the UHI.
Summarizing, the wind patterns with the most significant frequencies were WIND-P1 and WIND-P2, with annual values of 30% and 21%, respectively. These two patterns, however, correspond to small velocities in comparison with the other five wind patterns. An outstanding feature of WIND-P1 is that, although its events occurred mainly during the night with converging winds, it also comprises daylight events of winds with the same characteristic of convergence towards the city downtown. We cannot interpret these last events as katabatic winds produced by the mountain-valley system; instead, we must understand them as winds associated with convective upwards flows driven by the urban heat island (UHI) effect. The pattern WIND-P6 also occurred during the night, but it peaked at midnight. This wind pattern comprises the downslope winds of the first hours of the night and also reflects the effect of the late afternoon northerly gap winds. The patterns WIND-P3 and WIND-P5 reflect the well-known prevailing winds from N and NE of the Mexico City region, and WIND-P2 evidence the occurrence of the trade winds.
The wind patterns with the most significant frequencies were WIND-P1 and WIND-P2, with annual values of 30% and 21%, respectively. These two patterns, however, correspond to small velocities in comparison with the other five wind patterns. An outstanding feature of WIND-P1 is that, although its events occurred mainly during the night with converging winds, it also comprises daylight events of winds with the same characteristic of convergence towards the city downtown. We cannot interpret these last events as katabatic winds produced by the mountain–valley system; instead, we must understand them as winds associated with convective upwards flows driven by the urban heat island (UHI) effect. The pattern WIND-P6 also occurred during the night, but it peaked at midnight. This wind pattern comprises the downslope winds of the first hours of the night and also reflects the effect of the late afternoon northerly gap winds. The patterns WIND-P3 and WIND-P5 reflect the well-known prevailing winds from N and NE of the Mexico City region, and WIND-P2 evidence the occurrence of the trade winds.
Our wind patterns agree quite well with the Mexico City wind patterns obtained in [
25] for the period 2001–2006 with a slightly different methodology. In
Table 3, we described the correspondence between these two sets of patterns. The first column contains the wind patterns WP1-WP7 reported in [
25]; the second column presents the wind patterns obtained in the present work; and, the last column provides a brief description of the main characteristics of these patterns.
The main differences between both sets of patterns were:
The pattern WIND-P4 includes the wind patterns WP4 (southerly winds) and WP5 (westerly winds) reported in [
25]. However, we detected the pattern WIND-P4 only during the winter season, while the pattern WP4 occurred throughout the year (although with tiny frequencies in spring, summer, and autumn), and the pattern WP5 was only observed during the first semester of the year (winter and spring seasons).
We recognized a pattern (WIND-P7) associated with cyclonic winds strongly converging towards Mexico City downtown during the daylight hours. The winds in this pattern indicate an upwards flow driven by the UHI. No report exists in [
25] about a similar pattern.
3.2. Pollutant Flow Patterns
In
Table 4,
Table 5,
Table 6 and
Table 7, we enumerated the pollutant-flow patterns (NO
2, O
3, PM
10, and SO
2) that we identified for the Mexico City region, their seasonal and annual frequencies and flow intensities (the magnitudes of the flow vectors), and the associated wind circulation patterns. In
Figure 9,
Figure 10,
Figure 11 and
Figure 12, for each pollutant-flow pattern, we presented the average population plot, the mean pollutant-flow vectors at the quadrants NW, NE, SW, and SE of the city, the pollutant surface distribution averaged over the pattern, and its mean spatial value. We denoted the pollution-flow patterns as POL-Pn, where POL indicates the pollutant (POL = NO2, O3, PM10, and SO2), and n is a consecutive integer number from 1 to 9. The cells with the value 0.00 in the tables indicate the seasonal periods where some flow patterns were not detected.
In general, although the pollutant concentrations modulate the magnitudes of the flow vectors, the pollutant-flow patterns reflected the directional characteristics inherited from the wind patterns.
3.2.1. The NO2 Flow Patterns
We identified nine NO
2 flow patterns for the Mexico City region.
Table 4 shows the seasonal and annual average frequencies and average flow intensities of the patterns, expressed in percent (%) and µgm
−2s
−1, respectively. As shown in
Figure 9, these flow patterns closely resemble the wind circulation patterns, and
Table 4 summarized the relationship between the pollutant-flow and wind circulation patterns.
The NO2 flow patterns with the most significant annual frequencies were NO2-P1 (30%), NO2-P4 (20%), and NO2-P7 (12%), but the patterns with the most substantial annual flow intensities were NO2-P5, NO2-P7, and NO2-P9 with 25.4, 25.3, and 24.5 µgm−2s−1, respectively.
The flow patterns NO2-P4, NO2-P7, NO2-P8, and NO2- P9 carry nitrogen dioxide from the north to the south quadrants of the city, particularly the pattern NO2-P9, although it occurred only during the second semester of the year. The events of these patterns take the ozone precursor to the SW and SE quadrants, contributing actively to the ozone formation in this area throughout the year.
The patterns NO2-P2 and NO2-P5 carry the pollutant from the eastern to the western quadrants of the city following the trade winds. The events of the NO2-P3 and NO2-P6 patterns, differently, take nitrogen dioxide from west to east on the west side of the city, but from south to north on the eastern side. However, while the pattern NO2-P3 reflects a coupling of the westerly winds and the afternoon southerly gap winds, the pattern NO2-P6 reflects cyclonic transport driven by the UHI.
During the night, the events of the pattern NO2-P1, follow the downslope winds from the surrounding mountains and carry NO
2 to the city downtown. In comparison with the patterns of other pollutants (O
3, PM
10, and SO
2), the NO2-P1 flow pattern is the only one that revealed considerable nocturnal transport due to the downslope winds. The main reason is that NO
2 accumulates during the night since there are no photochemical reactions that consume it, and some of the mobile sources, which are responsible for approximately 86% of nitrogen oxide emissions in the Mexico City region [
3], remain active during the night.
The NO2 surface distributions, which are averages over the events of the NO2 flow patterns, reveal a North–South part of the city as the zone with the higher NO2 concentrations. This zone extends around the N–S axis that separates the west from the east quadrants (although slightly shifted to the eastern quadrants). The highest NO2 levels occurred, in general, in the south sector, close to the Huizachtepetl (Cerro de la Estrella).
The surface NO2 distributions with the highest concentrations were in connection with the NO2-P5 and NO2-P6 flow patterns, with mean spatial levels of 46 and 41 ppb. These high levels of NO2 were detected mainly during the winter and autumn. We note that, in these two cases, the high levels of NO2 extended spatially, covering almost completely the city. In the first case, it was a consequence of blocked transport by the mountains of the Sierra las Cruces, while in the second case, the high levels were a consequence of the cyclonic wind convergence driven by the UHI.
3.2.2. The O3 Flow Patterns
In
Table 5 and
Figure 10, we summarized the main characteristics of the O
3 flow patterns that we identified. In
Table 5, we included the average seasonal and annual occurrence frequencies (%) and flow intensities (µg m
−2s
−1) of the flow patterns, and the wind patterns from which they inherited their circulation characteristics.
The flow vectors of the diurnal flow patterns O3-P6, O3-P8, and O3-P9 have considerable northerly components and convey ozone from the NW and NE to the SW and SE city quadrants, providing substantial contributions to the high ozone concentrations frequently observed in the south sectors (especially in SW) of the city throughout the year.
The flow patterns O3-P6, O3-P8, and O3-P9 revealed the highest flow intensities (43, 82, and 74 µgm−2s−1, on annual average). Seasonally, we observed the highest O3 flow intensities during the spring, summer, and autumn.
The flow patterns O3-P3 and O3-P5 bring ozone from East to West in the city following the trade winds (pattern WIND-P2). These patterns reveal flow vectors from NE in the west quadrants of the city, and because of the mountains barrier (Sierra Las Cruces and Sierra Ajusco Chichinautzin), they contribute to the ozone accumulation at SW sector of the Mexico City region. The O3 surface distributions, averaged over the events of these flow patterns, reveal, in fact, an evident and significant ozone accumulation (ranging from 70 to 100 ppb, approximately) at the SW quadrant.
The flow patterns O3-P1 and O3-P7 comprised nighttime flow events and had the first and second-largest occurrence frequencies. However, their flow intensities are tiny because of the low ozone concentrations (12 ppb and 22 ppb on average, respectively) in the nocturnal atmosphere (ozone production is of photochemical nature). It is interesting to note that the patterns O3-P1 and O3-P7 and the patterns NO2-P1 and NO2-P4 reveal opposite behaviors (during the night, one observes NO2 accumulation while O3 decreases to its lowest levels) due to its relationship through the atmospheric photochemistry.
The O3-P2 is a flow pattern detected during the evening, driven by the UHI winds. It brings ozone from the surrounding parts of the city towards downtown, where an upwards convective flow takes ozone out there again, keeping the ozone levels relatively small in the city.
During winter, the flow events of the O3-P4 pattern bring ozone from west to east in the west side quadrants with small flow intensities, but the same flow pattern carries a considerable amount of ozone from south to north in the east side quadrants of the city, following the ventilation channel located at the west side of Sierra Nevada.
In general, the O3 surface distributions show relative small concentrations where the NO2 surface distributions show large emissions and vice versa. It is an evident behavior because of the NOx-O3 photochemical interactions.
3.2.3. The PM10 Flow Patterns
We recognized nine PM
10 flow patterns, which we enumerated in
Table 6 and sketched out in
Figure 11. In
Table 6, we presented the seasonal and annual averages of the occurrence frequencies and the flow intensities of the PM
10 flow patterns, expressed in % and µg/m
2s, respectively. The directional characteristics of wind inherited by the PM
10 patterns are also briefly indicated in this table.
The PM10 flow patterns with the highest annual flow intensities were PM10-P7 (43 µgm−2s−1), PM10-P6 (39 µgm−2s−1), PM10-P5 (35 µgm−2s−1), and PM10-P8 (26 µgm−2s−1); but their annual frequencies were of the smaller (3%, 4%, 2%, and 5%, respectively). All these four patterns reveal a remarkable flow intensity at the NE city quadrant throughout the year, although with different flow directions. The patterns PM10-P6, PM10-P7, and PM10-P8 exhibited also flow vectors with intense northerly components in the NW quadrant, which convey particle PM10 from the north to the south in the west side of the city. The events of these patterns also carried particulate matter from the northeast to the south sectors of the city, particularly in the NE quadrant. The flow intensities of the flow pattern PM10-P6 were 71, 56, and 27 µgm−2s−1 during the winter, spring, and summer seasons; for PM10-P7 were 110 and 63 µgm−2s−1 during the spring and summer, and for PM10-P8 were 43 and 55 µgm−2s−1 during the summer and autumn, respectively.
The events of the PM10-P4 and PM10-P5 flow patterns carried particles from south to northwest in the quadrants of the east side of the city: PM10-P4 with flow intensities of 29 and 27 µgm−2s−1 during the winter and spring, and PM10-P5 with flow intensities of 67 and 71 µgm−2s−1 during the winter and autumn, respectively.
The events of the flow patterns PM10-P1 and PM10-P2 occurred throughout the year with the highest frequencies of occurrence (44% and 22%, on an annual average). The flow intensity of the PM10-P1 (a nocturnal pattern) was too small all year, but the flow intensity of the PM10-P2 (a diurnal pattern) was similar to those of the other flow patterns. The patterns PM10-P1 and PM10-P3 reflect the influence of the downslope winds driven by the mountain-valley system. These patterns, furthermore, reflect the effect of the urban heat island, which is revealed in the hourly population by the small peak around the hour 15 (daylight winds converging to the downtown).
For all the PM10 flow patterns, the average surface distributions of the pollutant reveal the NE quadrant as the zone of the city with the highest PM10 concentrations, and the SW quadrant as the zone with the lowest levels of PM10. It suggests that, at the NE of the town, in the surrounding area of Xalostoc on the east side of the Sierra de Guadalupe, there exists considerable sources of particulate matter, which release PM10 to the atmosphere all year long.
It is interesting to observe that the events of the pattern PM10-P5 revealed flow vectors at the NE and SE quadrants with intense southerly components, which produced, on average, the highest mean spatial concentration of PM10 (81 µgm−3) in the city. This result displays the east side quadrants of the city as the most polluted areas by PM10 during the winter and autumn.
3.2.4. The SO2 Flow Patterns
Table 7 and
Figure 12 enumerate and sketch out the SO
2 flow patterns that we identified.
Table 7 presents the average seasonal and annual occurrence frequencies and flow intensities of these patterns, including their relationship with the wind circulation patterns.
Figure 12 shows the hourly population of the patterns, the corresponding flow vectors, and the surface distribution of the SO
2 emissions produced by the events of these flow patterns.
We observed, in general, that all the SO2 flow patterns presented small flow intensities and produced small surface concentrations (we underline that the Mexican air quality standard for SO2 is 110 ppb on a 24 h average, and 25 ppb on an annual average). Nevertheless, it is particularly interesting to observe that the flow patterns SO2-P3 (winter, spring, and autumn), SO2-P4 (winter, summer, and autumn), SO2-P5 (winter and autumn), and SO2-P9 (spring and summer), which were detected mainly during the night hours, exhibit the flow vector at the NW quadrant with a flow intensity relatively larger than in the other quadrants. However, the corresponding wind patterns show winds of similar magnitudes in all the quadrants.
The events of the flow patterns SO2-P7 (spring, summer, and autumn) and SO2-P8 (spring and summer) occurred during daylight hours and show the NW-quadrant flow vector larger than in the other quadrants. Consistently with these observations, the mean spatial SO2 surface distributions reveal the NW quadrant of the city as the zone with the highest levels of concentration, particularly in the case of the patterns SO2-P4 and SO2-P5. This situation seems to indicate that in the NW quadrant of the city, there are intense activities that release sulfur dioxide to the atmosphere during the night.
The SO2-P1 and SO2-P2 flow patterns, as was the case with the other pollutants, were the most frequent patterns (33% and 19%, respectively, on annual average) and the only ones detected all over the year; however, their flow intensities were tiny.