Surface ozone concentrations in Mexico City frequently exceed the Mexican standard and have proven difficult to forecast due to changes in meteorological conditions at its tropical location. The Madden–Julian Oscillation (MJO) is largely responsible for intraseasonal variability in the tropics. Circulation patterns in the lower and upper troposphere and precipitation are associated with the oscillation as it progresses eastward around the planet. It is typically described by phases (labeled 1 through 8), which correspond to the broad longitudinal location of the active component of the oscillation with enhanced precipitation. In this study we evaluate the intraseasonal variability of winter and summer surface ozone concentrations in Mexico City, which was investigated over the period 1986–2014 to determine if there is a modulation by the MJO that would aid in the forecast of high-pollution episodes.
Over 1 000 000 hourly observations of surface ozone from five stations around the metropolitan area were standardized and then binned by active phase of the MJO, with phase determined using the real-time multivariate MJO index. Highest winter ozone concentrations were found in Mexico City on days when the MJO was active and in phase 2 (over the Indian Ocean), and highest summer ozone concentrations were found on days when the MJO was active and in phase 6 (over the western Pacific Ocean). Lowest winter ozone concentrations were found during active MJO phase 8 (over the eastern Pacific Ocean), and lowest summer ozone concentrations were found during active MJO phase 1 (over the Atlantic Ocean). Anomalies of reanalysis-based cloud cover and UV-B radiation supported the observed variability in surface ozone in both summer and winter: MJO phases with highest ozone concentration had largest positive UV-B radiation anomalies and lowest cloud-cover fraction, while phases with lowest ozone concentration had largest negative UV-B radiation anomalies and highest cloud-cover fraction. Furthermore, geopotential height anomalies at 250 hPa favoring reduced cloudiness, and thus elevated surface ozone, were found in both seasons during MJO phases with above-normal ozone concentrations. Similar height anomalies at 250 hPa favoring enhanced cloudiness, and thus reduced surface ozone, were found in both seasons during MJO phases with below-normal ozone concentrations. These anomalies confirm a physical pathway for MJO modulation of surface ozone via modulation of the upper troposphere.
Ozone is hazardous to human health (WHO, 2008) and is an ubiquitous problem
in many megacities around the world. Tropospheric ozone is a secondary
pollutant produced by complex photochemistry from anthropogenic emissions,
and high ozone events typically affect mid-latitude urban areas during
summer, while in the tropics, such events can be observed throughout the
year. The problem of the incidence of high surface ozone events is
exacerbated in Mexico City, a megacity with 21 million inhabitants, because
of the intense solar radiation received at its relatively high elevation
(more than 2200 m above sea level) and tropical latitude (19.4
Seasonal variability in maximum surface ozone concentrations is not large in
Mexico City due to its geographical location (Raga and LeMoyne, 1996). Both
in the dry winter (December–February) and wet summer (June–August) months,
clear skies and strong insolation in the morning hours promote rapid
generation of surface ozone via photochemical conversions from anthropogenic
precursor emissions near the surface. In both seasons, as the day
progresses, the boundary layer becomes unstable from solar radiation and
deepens, diluting pollutant concentrations near the surface. The growth of
the boundary layer in Mexico City occurs over the course of a few hours,
with typical heights reaching at least 1.2 km above the surface (Nickerson
et al., 1992; Perez Vidal and Raga, 1998), even during the winter months when
insolation is reduced at this latitude. Highest ozone concentrations during
the winter months are often seen on days with strong insolation and light or
no surface wind (Lei et al., 2007). In summer months, clouds and
precipitation generally reduce the number of days with extremely elevated
surface ozone concentrations. However, when large-scale atmospheric
conditions are favorable, such as when a high-pressure regime and associated
clear skies affect the Mexico City basin, elevated concentrations of surface
ozone are also recorded in summer (Raga and Le Moyne, 1996). Hourly surface
ozone concentrations routinely exceed the national standard, set at 110 ppb
in 1993 (by law NOM-020-SSA1-1993) and modified in 2014 to 95 ppb (by law
NOM-020-SSA1-2014). In 2015, hourly maximum O
The problem of air quality in Mexico City has been the subject of numerous field programs over the years, typically limited in time but more comprehensive in terms of the number of parameters measured. One such campaign was MILAGRO (Megacity Initiative: Local and Global Research Observations), a very large international field campaign that took place in March 2006. The results of the large number of publications from that project are summarized by Molina et al. (2010). These results provided new insight into several processes related with pollutant transformations and chemical pathways, emerging from the analysis of the data collected with the large suite of sophisticated instrumentation deployed and the modeling performed. However, intensive field campaigns limited to 1 month, cannot address the seasonal and intraseasonal variability of the high surface ozone within the city. Past studies have examined the variability of surface ozone in Mexico City at different timescales, e.g., hourly (Raga and Le Moyne, 1996; Huerta et al., 2004; Lei et al., 2007), daily (Fast and Zhong, 1998), weekly (Stephens et al., 2008), monthly (Rodríguez et al., 2016), and seasonal (Thompson et al., 2008). All of these studies noted a primary relationship between ozone concentration in Mexico City and ultraviolet (UV) radiation, where days with more UV radiation were associated with elevated surface ozone concentrations. Furthermore, UV radiation received at the surface is strongly modulated by cloud cover (El-Nouby Adam and Ahmed, 2016). However, as yet, no study has explored surface ozone variability in Mexico City on the intraseasonal (30–60 day) timescale, despite known relationships between the leading mode of atmospheric intraseasonal variability, the Madden–Julian Oscillation (MJO; Madden and Julian, 1971), and tropical cloud cover (Riley et al., 2011) and circulation (Madden and Julian, 1972; Zhang, 2005). The MJO is largely responsible for intraseasonal variability in the tropics. Circulation patterns in the lower and upper troposphere and precipitation are associated with the oscillation as it progresses eastward around the planet. It is typically described by phases (labeled 1 through 8), which correspond to the broad longitudinal location of the active component of the oscillation with enhanced precipitation.
In this study we evaluate the intraseasonal variability of winter and summer
surface ozone concentrations in Mexico City over the period 1986–2014 to
determine if there is a modulation by the MJO that would aid in the forecast
of high-pollution episodes. Based on the relationships between surface ozone
and UV radiation, UV radiation and cloud cover, and cloud cover and the MJO,
the primary hypothesis tested in this study was the following:
The physical pathway hypothesized to support this intraseasonal variability
was as follows:
Station names, locations, period of record, and number and type of observations.
The government monitoring network, Red Automática de Monitoreo Atmosférico (Automated Atmospheric Monitoring Network, RAMA) has been operational since January 1986 measuring all criteria pollutants, with instrumentation certified by the US Environmental Protection Agency (EPA). In particular, the instrument to measure ozone is produced by Thermo Environmental Instrument Model 49, by UV absorbance. The RAMA currently has 33 stations within the Mexico City basin, but only a few have records dating back to 1986.
We selected five stations with the longest periods of record (Table 1), one station from each of the five geographic regions in the metropolitan area identified by several previous studies and summarized by Raga et al. (2001). Hourly observations from Tlalnepantla (TLA, in the northwest sector of the city, NW), Xalostoc (XAL, in the northeast sector, NE), Merced (MER, in the Center), Pedregal (PED, in the southwest sector, SW), and Universidad Autónoma Metropolitana-Iztapalapa (UIZ, in the southeast sector, SE) were available beginning in January 1986 and up to December 2014. See Fig. 1 for station locations and Table 1 for numbers of observations and elevations of each station. Since the ozone time series were non-stationary, standard anomalies (also called normalized anomalies) were calculated by subtracting a mean value from each observation and then dividing that result by a standard deviation (Wilks, 2011). Those mean values and standard deviations for each hour were estimated applying a 30-day (approximately monthly) running window, and the 30-day period was selected to avoid influence from both seasonal variability and also the long-term trend. We did not stratify by day of the week based on Stephens et al. (2008), who found that ozone in Mexico City exhibited relatively little variability by day of the week. Furthermore, we defined a “low” ozone concentration day as one with mean afternoon (12:00 to 16:00 local time) ozone standard anomalies (averaged across the five observing stations) below the 10th percentile. Percentiles were determined separately for each season using standard anomalies on all days in that season from 1986 to 2014. Similarly, we defined a “high” ozone concentration day as one where mean afternoon ozone standard anomalies exceeded the 90th percentile, again calculating winter and summer percentiles separately.
Locations of RAMA surface ozone stations used in this study (colored dots; abbreviations defined in Table 1) and topographic height (shaded, in m) of the Mexico City metropolitan region. State boundaries shown as black contours. Surface meteorology station at Tacubaya (TCBY) also indicated. The inset in the upper right corner shows the location of Mexico City within Mexico.
The MJO phase was determined using the Real-time Multivariate MJO (RMM) index
(Wheeler and Hendon, 2004). The daily RMM is based on time series of two
principal components derived from empirical orthogonal functions of
equatorially (5
Values of geopotential height (in m) and
We selected the winter (December–February; DJF) and summer (June–August; JJA) seasons for this study because of the homogeneity in synoptic-scale weather patterns in those seasons. More details on the climatological variability of ozone in Mexico City can be found in Klaus et al. (2001).
Finally, daily values of surface wind at the Tacubaya station (TCBY in Fig. 1) were taken from the NOAA National Centers for Environmental Information (NCEI) Integrated Surface Database (ISD; Smith et al., 2011). Anomalies of those values, calculated with respect to seasonal means, were binned by MJO phase to give composite anomalies for each season. For UV and total cloud cover in Mexico City itself, the gridded ERA-Interim value at the point closest to the mean latitude and longitude of the five RAMA stations was selected.
The diurnal cycle of ozone concentrations at each of the stations exhibited a daily minimum around 07:00 LT just prior to sunrise and a peak between 12:00 and 15:00 LT, with highest concentrations at the southernmost stations (PED and UIZ) and lowest in the northernmost station (XAL) (Fig. 2a). Additionally, highest ozone concentrations occurred 1–2 h earlier in spring (March–May; MAM) than in winter (December–February; DJF) at both PED and XAL (Fig. 2b), and peak ozone at PED in the south occurred one to two hours after peak ozone in XAL in the north, as a result of weak northeasterly surface winds transporting ozone and photochemical precursors southward during the day (Bossert, 1997).
Mean ozone concentrations in spring were nearly 30 % higher at all
stations than the rest of the year (Fig. 3, with observations smoothed by a
30-day running mean), and the effects of increased UV radiation during the
“mid-summer drought” (
Annual cycles of surface ozone concentrations (ppb) for five observing stations for hours 12:00–16:00 LT (local time) from the RAMA network, 1986–2014. Observations are smoothed using a 30-day running mean.
One of the challenges in examining intraseasonal variability of ozone is the need for a stationary record over a long period. In Mexico City, ozone concentrations have steadily decreased from the early 1990s to the 2010s (Fig. 4a; also Rodríguez et al., 2016) as a result of pollution-control measures (Molina and Molina, 2004). In order to remove the long-term trend, while keeping the intraseasonal variability at hourly resolution, hourly observations were converted to standard anomalies as described in Sect. 2. Results of this transformation of hourly observations to standard anomalies for station PED are shown in Fig. 4a (original hourly observations) and Fig. 4c (hourly standard anomalies). Standard anomalies for the other four stations show very similar results.
We note that overnight minimum observations from 1991 to 1993 were probably overestimated in the observational record (Fig. 4a), an artifact also seen in the other four stations (not shown). However, because in this study we focused on afternoon values (from 12:00 to 16:00 LT), that potential overestimation did not materially impact our results.
By transforming each hourly observation into a standard anomaly, the
distribution of relative frequencies shifted from highly non-Gaussian, with
peaks near zero and very long right tails (Fig. 4b), to more Gaussian, with
peaks near
Before examining ozone variability by MJO phase, it was important to first establish the synoptic-scale patterns associated with days of low and high ozone concentrations (defined in Sect. 2) in each season.
In winter (DJF), the synoptic pattern on days with low afternoon surface
ozone concentration featured a 250 hPa ridge over northwest Mexico and the
southwest US (height anomalies up to
As in Fig. 5, but for summer (JJA) days.
The synoptic pattern for winter days with high surface ozone concentration
was opposite that for the low ozone days. Over northwest Mexico and the
southwest US, a trough was seen at 250 hPa (anomalies
Summer days with low surface ozone concentration featured a slight anomalous
ridge (height anomalies of
On an annual basis, afternoon (12:00 to 16:00 LT) surface ozone concentrations in Mexico City were found to vary by MJO phase. Highest ozone concentrations were noted on days when MJO was active and in phases 3, 4, and 5, while lowest ozone concentrations were noted on days when the MJO was active and in phases 1 and 2 (Fig. 7a). This variability was seen at all five stations, regardless of geographic position within the basin. Normalized anomalies of surface UV radiation and total cloud fraction from ERA-Interim reanalysis strongly supported the observed surface ozone variability: MJO phases with highest ozone concentrations also had highest UV anomalies and lowest total cloud fraction anomalies, while MJO phases with lowest ozone concentrations had the most negative UV and the most positive cloud fraction anomalies (Fig. 7d). We found this agreement remarkable, particularly because the two data sets independently presented the same intraseasonal pattern.
Mean standard anomalies of midday (hours 12:00–16:00 LT) surface
ozone concentrations by active MJO phase for
On a seasonal basis, surface ozone concentrations in Mexico City were also found to vary by MJO phase. However, the dependence on phase was found to change between winter and summer, meaning a phase associated with higher ozone concentrations in winter would not necessarily be associated with higher ozone concentrations in summer. We attribute these differences to seasonality in both the convective properties of the MJO itself (e.g., Zhang and Dong, 2004; Wu et al., 2006) and in the extratropical atmosphere, whose circulation the MJO modulates (Gloeckler and Roundy, 2013). Despite the phase-to-phase variability in maximum and minimum ozone concentrations throughout the year, in all seasons, there remained good agreement between phases with highest (lowest) ozone concentrations and phases with highest (lowest) UV and lowest (highest) total cloud fraction. That is, the sunnier phases were consistently associated with the highest ozone concentrations.
Relative frequency of extreme ozone days in winter (top two rows)
and summer (bottom two rows). A high ozone day was defined as one with a mean
afternoon (12:00 to 16:00 LT) ozone anomaly across the five observing
stations greater than the long-term (1986–2014) 90th percentile. Similarly,
a low ozone day was defined as one with a mean afternoon anomaly across the
5 observing stations less than the long-term 10th percentile. Values in bold
(winter phase 2; summer phase 6) indicate phases with highest mean ozone
concentrations in those seasons; values in italics (winter phase 8; summer
phase 1) indicate phases with lowest mean ozone concentrations in those
seasons. The number of days (
In winter months (DJF), highest ozone concentrations were found on days when the MJO was in phase 2, and lowest ozone concentrations were found on days when the MJO was in phase 8 (Fig. 7b). Highest UV radiation, and lowest total cloud fraction, were seen on days when the MJO was in phase 2, and lowest UV radiation and second-highest cloud fraction were seen on days when the MJO was in phase 8 (Fig. 7e). In summer months (JJA), highest ozone concentrations were found on days when the MJO was in phases 5, 6, and 7, and lowest ozone concentrations were found on days when the MJO was in phases 1 and 8 (Fig. 7c). Highest UV radiation, and lowest total cloud fraction, was seen on days when the MJO was in phase 6, and lowest UV radiation and highest total cloud fraction were seen on days when the MJO was in phase 1. In both winter and summer, UV radiation and cloud-cover anomalies strongly supported observed surface ozone anomalies, whereby the cloudiest MJO phases featured lowest ozone and the sunniest phases featured highest ozone. We again consider this agreement remarkable, given the independence of the ozone and reanalysis data sets. Summer months (JJA) featured the greatest range in mean ozone concentrations by MJO phase: a difference in 0.25 standard anomaly units between the phases with the highest ozone concentrations (phases 5 and 6) and the phases with the lowest ozone concentrations (phases 1 and 8) (Fig. 7c). Summer months also featured the largest spread in both UV and total cloud fraction standard anomalies (Fig. 7f).
An examination of the frequency of “extreme” ozone days in each MJO phase (here a day with an “extreme” ozone value was defined for each season as an afternoon standard anomaly either above the 90th percentile value or below the 10th percentile value) provides additional insight into the character of the MJO modulation of ozone. In both winter and summer, the phases associated with highest ozone concentrations (phase 2 in winter and phase 6 in summer) featured the fewest occurrences of days with extremely low ozone (days with concentrations below the 10th percentile; Table 2). Those phases also featured either the highest (in summer) or near-highest (in winter) occurrences of days with concentrations above the 90th percentile (Table 2). Furthermore, the phases associated with lowest ozone concentration (phase 8 in winter and phase 1 in summer) featured the highest occurrences of days with low ozone (Table 2) and below-normal occurrence of days with high ozone. These results confirm that one manner in which the MJO modulates ozone concentration in Mexico City is to reduce (or augment) the frequency of days with afternoon ozone concentrations either below the 10th or above the 90th percentiles.
Composites of 250 hPa height (in m), height anomaly (in m), and
mean wind
To examine physical mechanisms for the observed variability in ozone
concentration and cloud cover by MJO phase, composite anomalies of 250 hPa
height and
As in Fig. 8, but for summer days in active MJO
phase 1
Maximum winter ozone concentrations occurred on days when the MJO was active
and in phase 2, and on those days, a synoptic-scale pattern opposite to that
of phase 8 was seen: anomalous 250 hPa troughing was seen over northern
Mexico and the south-central US (height anomalies of
In JJA, minimum ozone concentrations occurred on days when the MJO was in
phase 1. In that phase, anomalous 250 hPa ridging was seen over northwest
Mexico and the southwest US (anomalies up to
The final physical variable examined for intraseasonal variability by MJO phase was the surface wind vector at 18:00 UTC (12:00 LT) at Tacubaya (TCBY in Fig. 1) in the center-west portion of the metropolitan area (Fig. 1). In winter, days in phase 8 (lowest ozone concentrations) featured anomalous northeasterly surface winds (blue vectors; Fig. 10), resulting in observed wind speeds up to 40 % stronger than climatology (red vectors in Fig. 10). Days in phase 2 (highest ozone concentrations) featured anomalous westerly winds, resulting in winds up to 50 % weaker in magnitude (Fig. 10) than climatology. In summer, days in phases 8 and 1 (lowest ozone concentrations) featured surface winds very similar to climatology in both magnitude and direction. In summer, the wind direction on days in phase 8 was more from the north-northwest, while climatology was from the north-northeast, resulting in a very small westerly anomaly. Days in phase 6 (highest ozone concentrations) also featured winds with similar direction as the seasonal mean, but with speeds up to 30 % faster (Fig. 10). Despite these variations by MJO phase across all seasons, we do not consider the surface wind anomalies to be physically consistent or representative of a large-scale pattern, for two reasons. First, because Mexico City is located in a basin, surface flow fields do not normally respond to synoptic-scale pattern variability (Stephens et al., 2008). Indeed, the majority of the day-to-day variability in surface wind speed and direction is controlled by mesoscale, thermally driven mountain–valley circulations (Doran et al., 1998). With the exception of “cold surge” events in winter that have been associated with cloudy days, the two dominant ozone patterns identified by de Foy et al. (2005) only served to identify whether the ozone maximum would be in the southern or northern parts of the metropolitan area. Second, the wind anomalies by MJO phase resulted in only subtle changes in either direction, or speed, or both (Fig. 10). Moreover, none of the wind anomalies identified in DJF would meet the northerly “cold surge” of de Foy et al. (2005), suggesting that the “cold surge” events can occur during different MJO phases unrelated to modulation from the MJO. Finally, the smallness of the surface wind variability by MJO phase supports our argument that variability in surface ozone concentrations by MJO phase are primarily driven by variability in total cloud cover and surface UV radiation, which in turn are related to anomalies in upper-tropospheric circulation.
Mean 10 m winds at Tacubaya station (TCBY in Fig. 1) at 18:00 UTC (12:00 LT). Mean surface wind vectors for each season, DJF and JJA, are on row one and indicated by red arrows. Mean (black arrows) and anomaly (blue arrows) vectors for the MJO phases associated with lowest surface ozone (phase 8 in DJF and phase 1 in JJA) are on the middle row. Mean (black arrows) and anomaly (blue arrows) vectors for the MJO phases associated with highest surface ozone (phase 2 in DJF and phase 6 in JJA) are on the bottom row. Note that the mean winds for low ozone in DJF and high ozone in JJA are very similar to the seasonal mean winds, so the anomaly (blue) vector is very small. All wind data are from NOAA National Centers from Environmental Information, 1986–2014.
In this study, we investigated the intraseasonal variability of winter (DJF) and summer (JJA) surface ozone concentrations in Mexico City. After standardizing over 1 000 000 hourly observations of surface ozone from five stations around the metropolitan area, we binned them by phase of the active MJO. We found that highest winter ozone concentrations occurred on days when the MJO was active and in phase 2 (in the Indian Ocean), and highest summer ozone concentrations occurred on days when the MJO was active and in phase 6 (in the western Pacific Ocean) in summer. Lowest ozone concentrations were found on winter days in MJO phase 8 (in the eastern Pacific Ocean) and summer phase 1 (in the Atlantic Ocean). This intraseasonal variability in surface ozone concentrations agreed well with anomalies in cloud cover and UV-B radiation: phases with highest ozone concentration had highest UV-B radiation and lowest cloud cover, while phases with lowest ozone concentration had lowest UV-B radiation and highest cloud cover. This agreement was found for both winter and summer. Circulation anomalies at 250 hPa were found to support the observed variability in ozone and cloud cover. In winter, height and circulation anomalies favoring reduced cloudiness, and thus elevated surface ozone, were found on days when the MJO was in phase 2, and height and circulation anomalies favoring enhanced cloudiness, and thus reduced surface ozone, were found on days when the MJO was in phase 8. In summer, monsoon-like 250 hPa circulation patterns that favor enhanced cloudiness, and thus reduced surface ozone, were found on days when the MJO was in phase 1, and 250 hPa circulation patterns opposite to the monsoon, favoring reduced cloudiness and thus elevated surface ozone, were found on days when the MJO was in phase 6. We did not find physically meaningful variability in surface wind direction by MJO phase, despite earlier studies suggesting a relationship between surface wind and surface ozone in Mexico City. This suggests that the intraseasonal variability in both summer and winter surface ozone by MJO phase is driven primarily by variability in cloud cover via modulation of upper-troposphere circulation.
Hourly ozone concentrations for Mexico City
are available from the governmental monitoring network and can be downloaded
from
The ERA-Interim reanalysis is available to the community courtesy of ECMWF
from
Partial funding for Bradford S. Barrett was provided by the Fulbright Scholar program of the US State Department and the Programa de Estancias de Investigación, Dirección General de Personal Académico, Universidad Nacional Autónoma de México (DGAPA-UNAM). The air quality data were obtained from the databases of the Mexico City's Air Quality Monitoring Network, operated by the Ministry of Environment of Mexico City. ERA-Interim data were provided courtesy of ECMWF. Edited by: N. Harris Reviewed by: two anonymous referees