A combined surface and tropospheric ozone climatology and interannual variability study was performed for the first time using co-located ozone photometer measurements (2013–2015) and tropospheric ozone differential absorption lidar measurements (2000–2015) at the Jet Propulsion Laboratory Table Mountain Facility (TMF; elev. 2285 m), in California.
The surface time series were investigated both in terms of seasonal and diurnal variability. The observed surface ozone is typical of high-elevation remote sites, with small amplitude of the seasonal and diurnal cycles, and high ozone values, compared to neighboring lower altitude stations representative of urban boundary layer conditions. The ozone mixing ratio ranges from 45 ppbv in the winter morning hours to 65 ppbv in the spring and summer afternoon hours. At the time of the lidar measurements (early night), the seasonal cycle observed at the surface is similar to that observed by lidar between 3.5 and 9 km.
Above 9 km, the local tropopause height variation with time and season impacts significantly the ozone lidar observations. The frequent tropopause folds found in the vicinity of TMF (27 % of the time, mostly in winter and spring) produce a dual-peak vertical structure in ozone within the fold layer, characterized by higher-than-average values in the bottom half of the fold (12–14 km), and lower-than-averaged values in the top half of the fold (14–18 km). This structure is consistent with the expected origin of the air parcels within the fold, i.e., mid-latitude stratospheric air folding down below the upper tropospheric sub-tropical air. The influence of the tropopause folds extends down to 5 km, increasing the ozone content in the troposphere.
No significant signature of interannual variability could be observed on the
2000–2015 de-seasonalized lidar time series, with only a statistically
non-significant positive anomaly during the years 2003–2007. Our trend
analysis reveals however an overall statistically significant positive trend
of 0.3 ppbv year
A classification of the air parcels sampled by lidar was made at 1 km intervals between 5 and 14 km altitude, using 12-day backward trajectories (HYSPLIT, Hybrid Single Particle Lagrangian Integrated Trajectory Model). Our classification revealed the influence of the Pacific Ocean, with air parcels of low ozone content (43–60 ppbv below 9 km), and significant influence of the stratosphere leading to ozone values of 57–83 ppbv down to 8–9 km. In summer, enhanced ozone values (76 ppbv at 9 km) were found in air parcels originating from Central America, probably due to the enhanced thunderstorm activity during the North American Monsoon. Influence from Asia was observed throughout the year, with more frequent episodes during spring, associated with ozone values from 53 to 63 ppbv at 9 km.
Ozone is an important constituent in the troposphere, impacting climate,
chemistry and air quality (The Royal Society, 2008). As a greenhouse gas
(Forster et al., 2007), it contributes to the Earth's global warming with an
estimated radiative forcing of 0.40
Several studies show that background ozone levels have increased significantly since preindustrial times (Mickley et al., 2001; Parrish et al., 2012; Staehelin et al., 1994; Volz and Kley, 1988) and these levels continued rising in the last decades in both Hemispheres (Derwent et al., 2007; Jaffe et al., 2004; Lee et al., 1998; Naja and Akimoto, 2004; Oltmans et al., 2006; Parrish et al., 2012; Simmonds et al., 2004; Tanimoto et al., 2009; Zbinden et al., 2006; Lelieveld et al., 2004). Nevertheless, after air quality regulations were implemented in the 1970s, the increasing trend has slowed down or even reversed in regions such as the eastern USA and Europe (Cooper et al., 2012, 2014; Granier et al., 2011). The situation is not the same for emerging economies such as Asia, where emissions are increasing with a corresponding increase in ozone levels (Dufour et al., 2010; Gao et al., 2005; Strode et al., 2015; Tie et al., 2009; Wang et al., 2006).
In most cases, variability and trend studies have revealed very large ozone variability with time, location and altitude (Cooper et al., 2014). This variability is mostly due to the large heterogeneity and variability of the ozone sources themselves, the different chemical processes affecting the formation and depletion of tropospheric ozone and its variable lifetime in the troposphere. Ozone atmospheric lifetime goes from a few hours in the polluted boundary layer to several weeks in the free troposphere, allowing it to travel over distances of intercontinental scale (Stevenson et al., 2006; Young et al., 2013). Additional factors that have been observed to influence tropospheric ozone variability are climate variability and related global circulation patterns such as ENSO (El Niño-Southern Oscillation) or PDO (Pacific Decadal Oscillation) (e.g., Lin et al., 2014, 2015a; Neu et al., 2014). Tropopause folds also play a key role on tropospheric ozone interannual variability, as they influence the ozone budget in the troposphere and can even affect air quality near the surface (e.g., Lin et al., 2015a; Brown-Steiner and Hess, 2011; Langford et al., 2012). In order to obtain statistically significant results and be able to assess tropospheric ozone interannual variability and trends, a large long-term monitoring data set with global coverage is required. In the last decades, efforts have been made in this respect and the number of tropospheric ozone measurements has considerably increased throughout the globe. However, it is still necessary to increase the current observation capabilities to characterize tropospheric ozone variability more accurately.
Long-term records of tropospheric ozone have been available since the 1950s
(Feister and Warmbt, 1987; Parrish et al.,
2012), but it is not until the 1970s that the number of ozone monitoring
stations became significant (Cooper et al., 2014 and references therein).
Currently, a considerable number of ozone monitoring sites are operating as
part of regional networks or international programs (e.g., World
Meteorological Observation Global Atmosphere Watch WMO/GAW, Acid Deposition
Monitoring Network in East Asia EANET, Clean Air Status and Trends Network
CASTNET). In addition to these ground-based networks, tropospheric
ozone measurements from satellite (TOMS, TES, OMI, etc.) or aircraft
(MOZAIC/IAGOS) platforms have been successfully implemented. Nevertheless, a
large fraction of the tropospheric ozone measurements are still only surface
or column-integrated measurements whilst the number of them with information
on the vertical coordinate is very small. Until recently, mainly ozonesonde
profiles have been used to provide altitude-resolved ozone variability
information in the troposphere (Logan,
1994; Logan et al., 1999; Naja and Akimoto, 2004; Oltmans et al., 1998,
2006, 2013; Newchurch et al., 2003), but the cost of an ozonesonde launch
has kept the sampling interval to one profile per week (or less) for a given
location. Ozone vertical profiles have also been obtained from aircraft
platforms through programs such as MOZAIC and IAGOS, available since 1995
(e.g., Zbinden et al., 2013; Logan et al., 2012). However, aircraft data are
limited to air traffic routes and the temporal resolution depends on the
frequency of the commercial flights. Differential absorption lidar (DIAL)
systems, which started to be used to measure tropospheric ozone in the late
1970s (Bufton et al., 1979; Proffit and Langford, 1997),
complement the ozonesonde and aircraft records, providing higher temporal
resolution thanks to their inherent operational configuration (from minutes
to days of continuous measurements). Currently, tropospheric ozone lidars
are still very scarce, but the implementation of observation networks such
as the international Network for the Detection of Atmospheric Composition
Change (NDACC;
As part of NDACC and TOLNet, a tropospheric ozone DIAL system located at TMF has been operating since 1999. In this study, an analysis of 16 years of lidar profiles measured at the station is presented together with the analysis of the surface ozone measurements that have been available at the site since 2013. The objective is to provide the first-ever published study of tropospheric ozone variability above TMF using both the surface and lidar data sets. The work presented here is particularly valuable due to the rising interest in the detection of long-term trends in the western USA and the scarcity of long-term measurements of ozone vertical profiles in this region. The high-terrain elevation and the deep planetary boundary layer of the intermountain western USA region facilitate inflow of polluted air masses originating in the Asian boundary layer and ozone-rich stratospheric air down to the surface, thus highly influencing air quality in the region (Brown-Steiner and Hess, 2011; Cooper et al., 2004; Langford et al., 2012; Liang et al., 2004; Lin et al., 2012a, b; Stohl, 2002). After a brief description of the instrumentation and data sets (Sect. 2), an analysis of the seasonal and interannual variability of tropospheric ozone above TMF for the period 2000–2015 will be presented in Sect. 3. The study includes a characterization of the air parcels sampled by lidar by identification of the source regions based on backward trajectories analysis. Concluding remarks are provided in Sect. 4.
TMF is located in the San Gabriel Mountains, in southern California
(34.4
The instrument temporal sampling can be set to any value from a few seconds to several hours and the vertical sampling can be set to any multiple of 7.5 m, depending on the science or validation need. For the routine measurements contributing to NDACC over the period 1999–2015 and used for the present work, the standard settings have typically ranged between 5 and 20 min for temporal sampling, and between 7.5 and 75 m for the vertical sampling. Profiles routinely archived at NDACC are averaged over 2 h, with an effective vertical resolution varying from 150 m to 3 km, decreasing with altitude. These temporal and vertical resolution settings yield a standard uncertainty of 7–14 % throughout the profile. The system operates routinely at nighttime, but daytime measurements with reduced signal-to-noise ratio are occasionally performed in special circumstances such as process studies, and aircraft or satellite validation. The total number of routine 2 h ozone profiles used in this study and archived at NDACC for the period 2000–2015 is included in Table 1.
Number of measurements, by month and years, performed at TMF with the tropospheric ozone DIAL system. NA indicates data not available at the time of the study.
The TMF ozone lidar measurements have been regularly validated using
simultaneous and co-located electrochemical concentration cell (ECC) sonde
measurements (Komhyr, 1969; Smit et al., 2007). In the troposphere the
precision of the ozonesonde measurement is approximately 3–5 % with
accuracy of 5–10 % below 30 km. TMF has had ozonesonde launch
capability since 2005 and 32 coincident profiles were obtained over the
period 2005–2013. Results from the lidar and the ECC comparison are included
in Fig. 1. Figure 1a shows the averaged relative difference between the lidar
and ECC ozone number density profiles for the 32 cases. The lidar and sonde
measurements are found to be in good agreement, with an average difference of
7 % in the bulk of the troposphere and most of the values under 10 %
(Fig. 1b), which is within the combined uncertainty computed from both the
lidar and sonde measurements. Note that a non-negligible fraction of the
differences is due to the different measurement geometry of the lidar and
ozonesonde: 2 h averaged, single location for lidar, and horizontally
drifting 1 sec measurements for the ozonesonde usually rising at
5 m s
Continuous surface ozone measurements have been performed at TMF since 2013 using the UV (ultraviolet) photometry technique (Huntzicker and Johnson, 1979) with a UV photometric ozone analyzer (Model 49i from Thermo Fisher Scientific, US). The operation principle is based on the absorption of UV light at 254 nm by the ozone molecules (Sinha et al., 2014). The instrument collects in situ air samples at 2 m above ground taken from an undisturbed forested environment adjacent to the lidar building. It provides ozone mixing ratio values at 1 min time intervals with a lower detection limit of 1 ppbv. Uncertainty has been reported to be below 6 % in previous studies (Sinha et al., 2014).
Figure 2a shows the surface ozone seasonal cycle at TMF and nearby stations from the California Air Resources Board (ARB) air quality network for the period 2013–2015. The seasonal cycle at TMF comprises a maximum in spring and summer and a minimum in winter, consistent with the ARB stations shown, as well as other stations in the US west coast (e.g., Schnell et al., 2015). Nonetheless, the seasonal cycle obtained at TMF from the hourly samples (left plot) presents larger ozone values and lower variability throughout the year compared to the other ARB stations, all of which are at lower altitudes. The mean surface value for the complete period at TMF is 55 ppbv, whereas the seasonal values are 57, 57, 52 and 45 ppbv in spring (March–April–May), summer (June–July–August), fall (September–October–November) and winter (December–January–February), respectively. These values are in good agreement with those obtained from surface measurements at high elevation sites in the Northern Hemisphere and reported in the review by Cooper et al. (2014). When using the 8hMDA (8 h maximum daily average; right plot), larger seasonal-cycle amplitudes occur, especially at stations affected by anthropogenic pollution such as Crestline or San Bernardino. These polluted stations present larger values in summer than those recorded at high-elevation remote stations like Joshua Tree or TMF. The mean 8hMDA at TMF is 58 ppbv and the seasonal averages are 62, 66, 57 and 49 for spring, summer, fall and winter, respectively. The observed low seasonal variability is typical of high-elevation remote sites with low urban influence (Brodin et al., 2010). A similar behavior can be observed at the Phelan, Joshua Tree or the Mojave National Preserve stations, all sites being at high elevation with low or negligible urban influence. In Fig. 2a a secondary minimum is observed at TMF and most of the ARB nearby stations in July–August, followed by a secondary maximum in fall.
In Fig. 2a a clear combined effect of the altitude and proximity to
anthropogenic pollution sources on the ozone levels is observed. In general,
higher ozone levels and lower variability are observed at higher altitudes.
The lowest altitude Pico Rivera instrument measures the lowest ozone levels,
and the highest-altitude TMF instrument measures the highest ozone levels
throughout the year when considering the hourly sampled data set. A mean
difference of
The difference between the seasonal cycle retrieved from the 1 h averaged data and the 8hMDA can be easily explained from the differences in the daily cycles at the different stations. The mean surface ozone diurnal cycle at TMF and nearby ARB stations is shown in Fig. 2b for the four seasons. Minimum values are observed at nighttime, whereas maxima appear in late afternoon. As for the seasonal cycle, the daily cycle at TMF, Joshua Tree, Mojave National Preserve and Phelan stations exhibit low variability compared to the other stations located at lower altitude and more affected by urban pollution. On average, daily values are larger at high-elevation remote sites such as TMF or Joshua Tree. However, the afternoon maximum is larger at polluted stations such as Crestline, especially in the summer season. In addition, the maximum at TMF and the ARB stations of Joshua Tree and Mojave National Preserve occurs later than at the other stations. The difference in timing is likely due to the different chemical species involved in the ozone formation and depletion processes due to the low influence of anthropogenic pollution (Brodin et al., 2010; Gallardo et al., 2000; Naja et al., 2003). In winter, a minimum is observed at TMF in the afternoon instead of the maximum observed at the other stations. This difference in diurnal pattern has been observed at other remote or high-elevation sites and has been attributed to the shorter day length and the lack of ozone precursors compared to urban sites. The resulting daytime photochemical ozone formation is insufficient to produce an ozone diurnal variation maximizing in the afternoon (Brodin et al., 2010; Gallardo et al., 2000; Naja et al., 2003; Oltmans and Komhyr, 1986; Pochanart et al., 1999; Tsutsumi and Matsueda, 2000).
The red curve in Fig. 3a (left plot) shows the average ozone profile in the
troposphere and the UTLS (upper troposphere–lower stratosphere) region
obtained by the TMF lidar for the period 2000–2015. The cyan horizontal bars
show the corresponding standard deviation at
The seasonally averaged profiles are shown in Fig. 3b. These averages represent larger values in spring and summer in the troposphere, whereas in the stratosphere maximum values occur in winter and spring. Within the troposphere, below 9 km, the seasonally averaged profiles show average values of 62, 60, 51 and 50 ppbv in spring, summer, fall and winter, respectively. These values are in good agreement with the average ozone concentrations (50–70 ppbv) obtained in previous studies (Thompson et al., 2007; Zhang et al., 2010) above the western USA. In the altitude range 9–16 km (UTLS) a much larger variability in ozone is observed, as indicated by the large standard deviation (left plot) and the differences between the seasonally averaged profiles (right plot). This large variability results from the horizontal and vertical displacement of the tropopause above the site, causing the lidar to sound either the ozone-rich lowermost stratosphere or the ozone-poor sub-tropical upper troposphere for a given altitude.
Composite monthly mean ozone mixing ratio (2000–2015) computed from the TMF lidar measurements. The dashed line indicates the climatological tropopause above the site (WMO definition). Bottom strip: composite monthly mean ozone mixing ratio (2000–2015) from the surface measurements.
The two-dimensional (2-D) color contours of Fig. 4 show the composite (2000–2015) monthly mean ozone climatology measured by lidar (main panel, 4–20 km). A similar 2-D color contour representation was used just below the main panel to represent the composite (2013–2015) monthly mean surface ozone. The climatological tropopause height at TMF is also included in the main panel (blue dotted line), with mean values ranging between 12 and 15 km. As discussed previously in this paper, the tropopause height variability is the main cause of the larger standard deviation observed in Fig. 3a in this region. Between the surface and 9 km, a very consistent seasonal pattern occurs, with maximum values in April–May and minimum values in winter. The spring–summer maximum in the free-troposphere has been consistently observed at other stations in Europe and North America and is commonly attributed to photochemical production (Law et al., 2000; Petetin et al., 2015; Zbinden et al., 2006). The maximum values in the western USA are also usually related to the influence of Asian emissions reaching the US west coast (Jaffe et al., 2003; Parrish et al., 2004; Cooper et al., 2005b; Neuman et al., 2012; Zbinden et al., 2013). Above 9 km, the seasonal maximum occurs earlier (i.e., in March and April between 10 and 12 km and February and March at higher altitudes) consistent with the transition towards a dynamically driven lower-stratospheric regime. At these altitudes, the ozone minimum is also displaced earlier in the year (August–October), which is consistent with the findings of Rao et al. (2003) above Europe.
The TMF surface and lidar data are found to be very consistent, both in terms of seasonal-cycle phase and amplitude, and in term of absolute mixing ratio values. The mean value obtained from the lidar measurements in the troposphere is very similar to the mean value obtained from the surface measurements (around 55 ppbv). This consistency points out that the TMF surface measurements are representative of the lower part of the free troposphere (i.e., below 7 km), at least during the nighttime lidar measurements. This is mostly due to the fact that the station is not affected by the boundary layer during most of the time because of its high-elevation. Additional daytime lidar measurements will be performed in 2016 to assess whether such consistency also exists at other times of the day, especially in the afternoon.
The 2000–2015 time series of the de-seasonalized ozone mixing ratio is shown in Fig. 5. Anomalies, expressed in percent, resulted from subtracting the climatological ozone monthly mean profiles computed for the period 2000–2015 to the measured lidar profiles. Large ozone variability with time is clearly observed, highlighting the difficulty in identifying trends and patterns. No clear mode of interannual variability is observed for the analyzed period here. However, positive anomalies seem to predominate throughout the troposphere during the period 2003–2007, especially below 7 km. On average, ozone mixing ratio values in the lower troposphere were 5 ppbv larger in 2003–2007 than during the entire period 2000–2015.
De-seasonalized ozone mixing ratio above TMF. Anomalies (in %) were computed with respect to the climatological (2000–2015) monthly means.
Time series of the median (blue), 5th (orange) and 95th (yellow) percentile ozone values at different altitude layers for the full year (top) and for selected seasons and altitude layers (bottom) obtained from the TMF lidar measurements. Dashed lines represent the linear fit for each time series.
Following a procedure similar to that described in Cooper et al. (2012), a
trend analysis was performed at different altitude levels (Tables 2 and 3
and Fig. 6). Figure 6 shows the time series of the median, 95th and
5th percentile values, obtained every year between 2000 and 2015 for
different layers and different seasons using the lidar profiles measured at
TMF. In order to obtain the trends, linear fits (shown in Fig. 6) of the
median, 95th and 5th percentiles were performed independently
using the least squares method. The ozone rate of change in ppbv year
The calculated trends depend on altitude and season. Table 2 contains the
ozone rate change expressed in ppbv year
Analyzing each season separately, a significant positive trend occurs in the
upper troposphere (7–10 km) for both spring and summer, with an ozone
increasing rate of 0.71
Ozone mixing ratio trends for the median, 5th and 95th percentiles
(P. is percentiles in table) over the period 2000–2015 as shown in Fig. 6
(see text for details) in ppbv year
The positive trend at TMF in spring for the median values is larger than the
trend obtained by Cooper et al. (2012) for the free troposphere in 1995–2011
(0.41
The springtime positive trend estimates reported in the western USA oppose ozone decrease in the eastern part. These results indicate that the 2-decade-long efforts to implement regulations to control air quality and anthropogenic emissions in the USA have led to a clear decrease in ozone levels in the eastern USA, but not in the western USA (e.g., Copper et al., 2012, 2014). This different regional behavior has been attributed to the inflow of elevated ozone, mainly from East Asia, and to the increasing contribution of stratospheric intrusions (Cooper et al., 2010; Jacob et al., 1999; Parrish et al., 2009; Reidmiller et al., 2009; Lin et al., 2012a, 2015a: Lefohn et al., 2011, 2012). But again, differences in sampling can impact significantly the interpretation of our trend estimates. As pointed out by Lin et al. (2015b), further coordination efforts at both global and regional scales are necessary in order to reduce biases introduced by inhomogeneity in sampling.
Standard errors in ppbv year
In an attempt to characterize the air parcels sounded by lidar above TMF
based on their travel history and analyze the influence of the different
source regions on the ozone profiles, 12-day backward trajectories ending at
TMF between 5 and 14 km altitude were computed using the HYSPLIT4 model
(Draxler and Rolph, 2003;
Our trajectory analysis comprises two steps. First, the 12-day backward trajectories computed by HYSPLIT and ending at different altitude levels were grouped using the HYSPLIT clustering tool (Draxler et al., 2009) in order to identify the most significant paths followed by the air masses arriving over the station. Based on the results of this preliminary analysis, five main regions were identified: the stratosphere, the Asian boundary layer (ABL), the free-troposphere above Asia (AFT), Central America and the Pacific Ocean. Once these geographical areas were identified, we performed a classification of the air parcels according to the criteria described next.
An air parcel was classified as “Stratospheric” if the 12-day backward
trajectory intercepted the tropopause and resided at least 12 h above
the local tropopause. The tropopause height information comes from the
global tropopause height data derived once a day by the NOAA Physical
Sciences Division (
Next, the air parcels that were not classified as “stratosphere” were then classified as “Central America” for trajectories comprising a minimum residence time of 4 days within the area labeled “Central America” in Fig. 7. According to the sensitivity test, a 4-day residence time period is long enough to avoid the influence of additional source regions and short enough to avoid an underestimation of the “Central America” cases.
Geographical boundaries used to characterize the air parcels associated with the 12-day backward trajectories ending at TMF during the lidar measurements over the period 2000–2015.
The air parcels not classified as “stratosphere”, or “Central America” were then classified as Asian if they comprised a minimum residence time of 6 h within the area labeled as “Asia” in Fig. 7. The Asian trajectories are subdivided in ABL if they come from an altitude below 3 km and AFT if they come from altitudes above 3 km. According to the sensitivity test, a residence time of 6 h is enough to clearly identify the signature of Asian emissions on the ozone profiles observed at TMF.
Examples of HYSPLIT 12-day backward trajectories arriving at TMF at 7 km altitude for four selected seasons and categories (see text for details).
Number of air parcels ending at TMF during lidar measurements over the period 2000–2015, classified as “stratosphere”, “Central America”, “ABL” (Asian boundary layer), “AFT” (Asian free troposphere), “Pacific Ocean” and “RT” (residual trajectories) (see text for details).
The air parcels not classified in any of the previous categories were classified as “Pacific Ocean” if a minimum residence time of 276 h (11.5 days) within the area labeled “Pacific” in Fig. 7 was reached. A residence time of 276 h guarantees that no influence from additional sources affects the air masses reaching TMF and the “Pacific” region can be considered as a background region.
Trajectories that did not match any of the previous categories were grouped as “residual trajectories” (RT). They will be considered for statistical purposes, but not for the analysis of the ozone mixing ratio values.
Distribution of the five categories identified for each trajectory ending at TMF during the lidar measurements over the period 2000–2015. The number of occurrences is given in percentage for each month of the year, and for four different altitude layers.
Box plot of the ozone mixing ratios measured within the air masses arriving at TMF at 9 km for the five identified categories (see text for details) and the four seasons. The black dot represents the mean value, the red line is the median and the box limits correspond to the 25th and 75th percentiles. The numbers between parentheses indicate the number of associated trajectories.
The classification of the air parcels took place sequentially, which means that each category is exclusive from the others. The classification was made for each of the four seasons separately in order to account for the seasonal changes in synoptic circulation. Examples of the corresponding classified back trajectories are shown in Fig. 8. The number and frequency of occurrences of each air parcel category for all seasons is compiled in Table 4. A monthly distribution of these occurrences is shown in Fig. 9. With the selection criteria we have set, air masses are predominantly associated with the “AFT” region below 11 km, ranging between 32 and 42 % from 5 to 11 k with maximum number of cases in spring. A very low number of parcels classified as “ABL” are found (between 0 and 9 %). Increasing influence of the stratosphere is observed at upper levels, with values increasing from 5 % at 5 km to 80 % at 14 km. Higher influence is observed during winter and spring, which agrees well with previous studies in the western USA (Sprenger, 2003; Stohl, 2003). A statistically significant Central American influence was identified in summer with a frequency of occurrence varying between 12 and 3 %, decreasing with altitude. The Central America influence coincides with the establishment of the North American Monsoon circulation from July to September which affects Central America and the southern USA.
Composite ozone profiles and statistical parameters were estimated for each category of air parcel and for altitudes between 5 and 14 km at 1 km altitude intervals. Figure 10 shows the ozone mixing ratio mean (open circles), median (red bars), 25th and 75th percentiles (blue bars) at 9 km altitude for each of the identified categories and season. The number of occurrences for each category is mentioned between parentheses. The ozone statistics obtained when a low number of occurrences was found should be ignored (e.g., Central America except for summer, or ABL for winter and fall). Figure 11a shows, for each season, the composite ozone profiles constructed from the ozone mixing ratio median values found for a particular category at a given altitude. The same profiles, but focused on the troposphere (5–10 km), are shown in Fig. 11b. In order to keep the most statistically significant results, composite values computed using less than 5 % of the total number of samples for a given season were not plotted, leaving out certain sections of the composite profiles.
Not surprisingly, the analysis reveals that the largest ozone mixing ratio values were mostly observed when the air masses were classified as “stratospheric” regardless of the season (median values between 17 and 35 ppbv larger than for the Pacific Ocean at 9 km). In spring and winter, the influence of the stratosphere goes down to 5 km, with ozone values ranging from 3 to 13 ppbv larger than for the Pacific category below 9 km. For this category, large ozone variability was found, as indicated by the 25th and 75th percentiles in Fig. 10. As altitude increases, the influence of the stratosphere is more important, exceeding 40 % above 12 km, resulting in higher ozone mixing ratio values (red curves in Fig. 11).
Conversely, low ozone mixing ratio values (40–61 ppbv below 9 km) were consistently associated with the air parcels classified as “Pacific Ocean” (cyan curves). This region can be considered as a source of “background ozone”, since no anthropogenic source is expected to affect the local ozone budget.
Higher ozone content (from 2 to 13 ppbv higher than for the Pacific region)
is systematically found for air parcels classified as “AFT”. Values are
especially larger in summer, when differences of at least 8 ppbv with the
Pacific region are found for altitudes between 5 and 13 km. In general, the
number of occurrences for air parcels classified as “ABL” remains very
small to provide any meaningful interpretation. Nonetheless, values in the
lower part of the troposphere during spring and winter, when the number of
occurrences is higher, are similar to those observed for the AFT. The
occurrence of the Asian air masses is mostly observed in spring (Figs. 9 and
10), and ozone associated with Asian emissions has been frequently detected
in the western USA during this season in previous studies (e.g., Cooper et
al., 2005b; Zhang et al., 2008; Lin et al., 2012b). Even though less
frequent, our results indicate that Asian pollution episodes observed during
summer are associated with larger ozone values than in spring. These larger
values are due to more active photochemical ozone production observed over
China in summer (Verstraeten et al., 2015), associated with larger ozone
values than those in spring. The influence of the air parcels classified as
Central America is mainly observed during summer, with ozone median values
5–28 ppbv larger than those observed for the Pacific region between 5 and
9 km (yellow curve in Fig. 11). Ozone mean values of 72 ppbv were found at
9 km altitude for the 74 air parcels classified as ”Central America”
(Fig. 10). The corresponding values for the 115 air parcels classified as
“Pacific Ocean” are about 52 ppbv, which is 20 ppbv lower. The larger
ozone values associated with the “Central America” category possibly points
to the lightning-induced enhancement of ozone within the more frequent
occurrence of thunderstorms during the North American summer monsoon.
Previous studies (Cooper et al., 2009), have observed enhanced ozone values
associated with the North American Monsoon, mainly due to ozone production
associated with lightning (Choi et al., 2009; Cooper et al., 2009). However,
this feature was observed in the eastern USA. Because of the synoptic
conditions during the monsoon, the western USA is not as much influenced and
no significant regional ozone increase was reported (Barth et al., 2012;
Cooper et al., 2009). Nevertheless, Cooper et al. (2009) reported higher
modeled lightning-induced NO
In the previous section, a large variability in the composite ozone content
was found for the air parcels classified as “Stratospheric”. In the
current section, we provide at least one clear explanation for this large
variability. Tropopause folds are found primarily in the vicinity of the
subtropical jets, in the 20–50
Double tropopauses are usually expected to result from tropopause folds in
the layer between the two identified tropopauses. Therefore, a common method
used in the literature to identify tropopause folds is to detect the presence
of double tropopauses based on temperature profiles (e.g., Chen et al.,
2011). The MERRA (Modern-Era Retrospective analysis for Research and
Applications; Rienecker et al., 2011) reanalysis data (1 km vertical
resolution, 1
Monthly distribution of occurrences (in %) of double tropopauses above TMF. The number of days with tropopause folds is normalized to the total number of measurements per month compiled in Table 1.
Using this methodology, we found that 27 % of the TMF tropospheric ozone lidar profiles were measured in the presence of double tropopauses. This high frequency of double tropopause occurrences was expected considering the latitude of TMF, i.e., near the subtropical jet, where frequent tropopause folds occur. Figure 12 shows the number of cases with double tropopauses above TMF distributed per months, with the number of days with double tropopause being normalized to the total number of measurements every month (compiled in Table 1). As we can see, the presence of double tropopauses was especially frequent during winter and spring, which coincides with the higher frequency of stratospheric air masses arriving at TMF estimated by the backward trajectories analysis (Fig. 9) and is in agreement with previous studies (Randel et al., 2007). The altitude of detected single tropopauses is found around 13 km in winter and spring, and 16–17 km in summer and fall (Fig. 13a–d). When a double tropopause is identified, the altitude of the lower tropopause ranges between 8 and 15 km, with the distribution peak centered around 12–13 km (Fig. 13e–h), and the second tropopause is detected typically around 17–18 km (Fig. 13i–l).
Figure 14a shows an example of an ozone profile measured on 8 January 2013, when a double tropopause was detected above TMF. The average of all tropospheric ozone lidar profiles measured in winter in cases of single tropopause is plot as reference. In Fig. 14b, the average of all tropospheric ozone lidar profiles measured in winter (blue curves) and spring (red curves) in the presence of a double tropopause (solid curves), and in the presence of a single tropopause (dashed curves) are included. The right panel (Fig. 14c) is simply a lower tropospheric-zoomed version of the middle panel (Fig. 14b). Only winter and spring are shown because they are the seasons most affected by double tropopause cases as previously stated. In the presence of double tropopauses a clear dual vertical structure in ozone is observed. For the specific case on 8 January 2013 (Fig. 14a), the lower tropopause was located at 9 km and stratospheric air reached down to approximately 6 km, considerably increasing the ozone content in the troposphere. On the other hand, ozone values were lower than the winter average in the lower stratosphere (11–19 km). In the case of the average profiles (Fig. 14b), the dual vertical structure presents higher ozone values between 12 and 14 km and lower mixing ratio values between 14 and 18 km. The dual ozone structure observed by lidar coincides with the expected location of the fold, and consists of systematically higher-than-average mixing ratios in the lower half of the fold (12–14 km), and lower-than-average mixing ratios in the upper half of the fold (14–18 km). This dual structure is consistent with the expected origin of the air masses within a tropopause fold. Stratospheric air, richer in ozone, is measured within the lower half of the fold, while tropospheric ozone-poor air is measured within the upper half of the fold.
In the case of deep stratospheric intrusions, ozone-rich stratospheric air masses embedded in the lower half of the fold can reach lower altitudes, and occasionally the planetary boundary layer mixing down to the surface (Chung and Dann, 1985; Langford et al., 2012, 2015; Lefohn et al., 2012; Lin et al., 2012a), leading to an ozone increase in the lower troposphere (Fig. 14b). In our case, the mean increase is around 2 ppbv below 6 km for both spring and winter. This increase is consistent with previous reports of the importance of the stratosphere as an ozone source in the lower troposphere (Cooper et al., 2005a; Langford et al., 2012; Lefohn et al., 2011; Trickl et al., 2011), with a 25 to 50 % contribution to the tropospheric budget (Davies and Schuepbach, 1994; Ladstätter-Weißenmayer et al., 2004; Roelofs and Lelieveld, 1997; Stevenson et al., 2006).
Combined ozone photometer surface measurements (2013–2015) and tropospheric ozone DIAL profiles (2000–2015) at the JPL-Table Mountain Facility were presented for the first time. The high ozone values and low interannual and diurnal variability measured at the surface, typical of high elevation remote sites with no influence of urban pollution, constitute a good indicator of background ozone conditions over the southwestern USA.
The 16-year tropospheric ozone lidar time series is one of the longest lidar records available and is a valuable data set for trend analysis in the western USA, where the number of long-term observations with high vertical resolution in the troposphere is very scarce. A statistically significant positive trend was observed in the upper troposphere, in agreement with previous studies. This ozone increase points out to the influence of long-range transport and/or a change in stratospheric influence, since ozone precursor emissions have been decreasing in the USA over the past 2 decades.
Influence of five main regions (stratosphere, Central America, Asian
boundary layer, Asian free troposphere and Pacific Ocean) on the ozone
profiles sampled above TMF was detected using 12-day backward trajectories.
This trajectories analysis revealed the large influence of the stratosphere,
especially in the UTLS and the upper troposphere, leading to high ozone
values. The influence of the stratosphere reached down to 5 km in spring and
winter, with ozone values ranging from 3 to 13 ppbv larger than for the
Pacific category, considered as a background region. In summer, enhanced
ozone values (5–28 ppbv larger than for the Pacific region) were found in
air parcels originating from Central America, probably due to the enhanced
thunderstorm activity during the North American Monsoon. Frequent air masses
coming from Asia were also observed, mainly in spring, associated with ozone
values 2 to 13 ppbv larger than those from the background region. Ozone
vertical distribution above TMF is also affected by the frequent occurrence
of tropopause folds. A dual vertical structure in ozone within the fold
layer was clearly observed, characterized by above-average values in the
bottom half of the fold (12–14 km), and below-averaged values in the top
half of the fold (14–18 km). Above-average ozone values were also observed
near the surface (
Part of the data used in this publication were obtained as part of the
Network for the Detection of Atmospheric Composition Change (NDACC) and are
publicly available (see
The work described in this paper was carried out at the Jet Propulsion
Laboratory, California Institute of Technology, under a Caltech Postdoctoral
Fellowship sponsored by the NASA Tropospheric Chemistry Program. Support for
the lidar, surface and ozonesonde measurements was provided by the NASA Upper
Atmosphere Research Program. The authors would like to thank M. Brewer,
T. Grigsby, J. Howe and members of the JPL lidar team, who assisted in the
collection of the data used here. The authors gratefully acknowledge the NOAA
Air Resources Laboratory (ARL) for the provision of the HYSPLIT transport and
dispersion model and/or READY website (