Articles | Volume 12, issue 9
Research article
06 Sep 2019
Research article |  | 06 Sep 2019

Ozone profile climatology for remote sensing retrieval algorithms

Kai Yang and Xiong Liu

New ozone (O3) profile climatologies are created from the Modern-Era Retrospective Analysis for Research and Applications version 2 (MERRA-2) O3 record between 2005 and 2016, within the period of Aura Microwave Limb Sounder (MLS) and Aura Ozone Monitoring Instrument (OMI) assimilation. These two climatologies consist of monthly mean O3 profiles and the corresponding covariances dependent on the local solar time, longitude (15), and latitude (10), which are parameterized by tropopause pressure and total O3 column. They are validated through comparisons, which show good agreements with previous O3 profile climatologies. Compared to a monthly zonal mean climatology, both tropopause- and column-dependent climatologies provide improved a priori information for profile and total O3 retrievals from remote sensing measurements. Furthermore, parameterization of the O3 profile with total column O3 usually reduces the natural variability of the resulting climatological profile in the upper stratosphere further than the tropopause parameterization, which usually performs better in the upper troposphere and lower stratosphere (UTLS). Therefore tropopause-dependent climatology is more appropriate for profile O3 retrieval for complementing the vertical resolution of backscattered ultraviolet (UV) spectra, while the column-dependent climatology is more suited for use in total O3 retrieval algorithms, with an advantage of complete profile specification without requiring ancillary information. Compared to previous column-dependent climatologies, the new MERRA-2 column-dependent climatology better captures the diurnal, seasonal, and spatial variations and dynamical changes in O3 profiles with higher resolutions in O3, latitude, longitude, and season. The new MERRA-2 climatologies contain the first quantitative characterization of O3 profile covariances, which facilitate a new approach to improve O3 profiles using the most probable patterns of profile adjustments represented by the empirical orthogonal functions (EOFs) of the covariance matrices. The MERRA-2 daytime column-dependent climatology is used in the combo O3 and SO2 algorithm for retrieval from the Earth Polychromatic Imaging Camera (EPIC) on board the Deep Space Climate Observatory (DSCOVR) satellite, the Ozone Mapping and Profiler Suite Nadir Mapper (OMPS-NM) on the Suomi National Polar Partnership (SNPP), and the Ozone Monitoring Instrument (OMI) on the Aura spacecraft.

1 Introduction

Remote sensing instruments measure spectral radiances, from which information about light absorbers such as ozone (O3) and other trace gases may be inferred using retrieval algorithms. Absorption signals carried in the measured radiances come from the interaction between photons and these absorbers, which are naturally distributed throughout the atmosphere. Consequently, for a band within the spectral range of significant atmospheric absorptions, its measured radiance is sensitive to the profiles, i.e., vertical distributions of light absorbers, and depends on other atmospheric state and surface variables. However, multispectral or even hyperspectral radiance measurements from nadir-viewing instruments do not have sufficient vertical resolution to fully disentangle the absorption signals for the determination of an absorber profile. Hence, the retrieval of quantitative information about an absorber requires some knowledge of its vertical distribution. Frequently, prescribed or a priori profiles are used to fill the knowledge gap on the altitudes from which absorption signals are originated but not differentiated by the measurements.

For retrieval of O3 from remote sensing measurements, the a priori knowledge is usually taken from an O3 profile climatology, which provides average O3 profiles and their variances. Most, if not all, total O3 algorithms (e.g., Mateer et al.1971; Klenk et al.1982; Bhartia and Wellemeyer2002; Coldewey-Egbers et al.2005; Eskes et al.2005; Veefkind et al.2006; Lerot et al.2010, 2014; Loyola et al.2011; Van Roozendael et al.2006, 2012; Wassmann et al.2015) rely on an O3 profile climatology to uniquely and completely specify the vertical distribution of a retrieved total column. Profile O3 algorithms (e.g., Hoogen et al.1999; Liu et al.2005; Wei et al.2010; Bhartia et al.2013; Miles et al.2015) based on the widely adopted optimal estimation (OE) inversion technique (Rodgers2000) require the a priori O3 profiles and their covariances to constrain the retrievals from deviating too far from the a priori O3 distributions. An OE-retrieved O3 profile is a combination, i.e., weighted average, of the real and the a priori profiles. Therefore, the accuracy of a total column or a profile retrieval is significantly affected by the selection of the a priori profile. A closer match between the a priori and the actual vertical distributions allows a higher accuracy retrieval to be achieved.

In this paper, we outline the characteristics of O3 vertical distribution and existing climatological data sets that are commonly used in many O3 retrieval algorithms. To improve a priori knowledge of O3 vertical distributions and its covariances and to address the deficiencies in existing O3 profile climatologies, we construct new O3 profile climatologies from the long-term global O3 field provided by a modern reanalysis system. We present comparisons to validate the new climatologies and summarize the present work in the last section.

2 Characteristics and climatologies of O3 vertical distribution

Figure 1Annual zonal mean O3 mixing ratio in parts per million by volume as a function of latitude and pressure altitude Z*. Panel (a) is from the new O3 profile climatology described in this paper, and panel (b) is from the ML climatology (McPeters and Labow2012). The dotted black line in each panel indicates the pressure altitude of the annual zonal mean tropopause. The pressure altitude Z* is defined as 16log10[psp], where p is pressure level (hPa) and ps=1013.25 hPa.


O3 profile climatologies (e.g., Wellemeyer et al.1997; Fortuin and Kelder1998; Bhartia and Wellemeyer2002; Lamsal et al.2004; McPeters et al.2007; McPeters and Labow2012; Wei et al.2010; Bak et al.2013; Sofieva et al.2014; Labow et al.2015) usually provide the a priori profiles for O3 retrieval from nadir-viewing satellite observations. Most of these climatologies are constructed by merging lower-atmosphere ozonesonde data with upper-atmosphere O3 profile measurements from one or more satellite instruments, such as the Microwave Limb Sounder (MLS) and the Stratospheric Aerosol and Gas Experiment II (SAGE-II) on the Upper Atmosphere Research Satellite (UARS), the MLS on Aura, or the Solar Backscatter Ultraviolet (SBUV) and SBUV/2 on NASA and NOAA satellites. All these climatologies capture the main characteristic of O3 vertical distribution, which is determined by the balance of the chemical processes of O3 production and destruction, as well as by atmospheric motions. Specifically, the vertical O3 distribution varies strongly with latitude, as shown in Fig. 1. The highest O3 mixing ratio is found at an altitude of 30–40 km in the equatorial region (see Fig. 1), in which most atmospheric O3 production induced by strong solar ultraviolet (UV) radiation takes place. The peak values of O3 profiles in the mixing ratio decrease with higher latitudes, as the large-scale meridional Brewer–Dobson circulation carries the O3 in the tropical stratosphere towards the poles, and slowly transfers the O3 to the lower stratosphere at middle and high latitudes. As a result of this atmospheric-circulation-driven process of transport, descent, and accumulation, O3 profiles in partial pressure or number density exhibit higher maximum values occurring at lower altitudes as latitude increases towards the poles (see Fig. 13). In addition to latitude dependence, some climatologies (e.g., Fortuin and Kelder1998; McPeters et al.2007; McPeters and Labow2012; Wei et al.2010; Bak et al.2013; Sofieva et al.2014) include monthly zonal mean O3 profiles to describe systematic profile changes caused by the seasonal variations in O3 photochemistry and atmospheric circulation, as well as significant hemispheric profile asymmetry resulting from hemispheric differences in orography, atmospheric temperature, and circulation transport (Maeda and Heath1983; Perliski and London1989; Cariolle et al.1992). While the seasonal dependence is accounted for in these O3 profile climatologies, they have not included the diurnal cycle of O3, which is significant in the upper stratosphere and mesosphere (Haefele et al.2008; Huang et al.2008; Sakazaki et al.2013; Schranz et al.2018).

Tropopause marks the interface between the stratosphere and the troposphere, across which there is a steep vertical gradient in O3 mixing ratio, which is higher in the stratosphere than in the troposphere (see Fig. 1). In the lower stratosphere (between the tropopause and ∼30 km altitude) where O3 lifetime is quite long (on the order of weeks or longer), the deviation of O3 vertical distribution from its climatological mean is mainly controlled by atmospheric dynamics. Consequently, O3 profile and column amount have large daily variations that are associated with meteorological conditions (Reed1950). In particular, the rise and fall of tropopause affect O3 columns and their vertical profiles directly through shifting the O3 mixing ratio gradient in the upper troposphere and the lower stratosphere (UTLS). This dynamical connection between total O3 and tropopause pressure (or altitude) has been investigated and documented in a number of studies (e.g., Ohring and Muench1960; Salby and Callaghan1993; Steinbrecht et al.1998; Krzyścin et al.1998; Weiss et al.2001; Varotsos et al.2004), which reveal a strong positive (negative) correlation between total O3 and tropopause pressure (altitude). This correlation implies that tropopause height (or pressure) and total O3 column amount are excellent indicators for selecting O3 profile shape in the UTLS region, where O3 profiles have the highest dynamical variability. To capture dynamical profile variations, O3 profile climatologies have included tropopause-sensitive zonal mean O3 profiles (Wei et al.2010; Bak et al.2013; Sofieva et al.2014), as well as column classifications for which zonal mean O3 profiles are compiled for a range of possible total O3 columns (Wellemeyer et al.1997; Bhartia and Wellemeyer2002; Lamsal et al.2004; Labow et al.2015).

The accuracy of the global O3 measurements from the series of Total Ozone Mapping Spectrometer (TOMS) is owed in large part to the a priori knowledge of O3 vertical distribution provided by the total-ozone-column-classified O3 climatology (Bhartia and Wellemeyer2002; McPeters et al.2007) created for the version 8 (V8) total O3 algorithm. For mapping between the O3 column and profile, many recent total O3 algorithms (e.g., Bhartia and Wellemeyer2002; Eskes et al.2005; Veefkind et al.2006; Van Roozendael et al.2006, 2012; Lerot et al.2010, 2014; Loyola et al.2011; Wassmann et al.2015) use the TOMS-V8 climatology, which is a combination of latitude and total O3-dependent profiles, also known as the standard profiles (Bhartia and Wellemeyer2002), and the latitude- and month-dependent Labow–Logan–McPeters (LLM) climatology (McPeters et al.2007). The standard profiles are 21 annual mean O3 profiles, covering the possible total O3 ranges in steps of 50 DU for three latitude zones: 225–325 DU at low latitudes, 225–575 DU at mid-latitudes, and 125–575 DU at high latitudes, and they are used to expand the latitude and month dependent LLM climatology to include the variation with total O3. While the LLM climatology provides latitude-dependent O3 profiles that capture the north–south asymmetry, the profile dependence on total O3 taken from the standard profiles is independent of season and makes no distinction between the Northern Hemisphere and Southern Hemisphere, therefore ignoring the hemispheric asymmetry in O3 profile deviation from the climatological mean. Furthermore, profile changes associated with total column variations represented by the standard profiles exhibit a large morphology difference between two adjacent latitude zones, as a result of these standard profiles being binned over wide (30) latitude zones. Consequently, total O3 algorithms relying on TOMS-V8 O3 profile climatology for profile specification usually have systematic latitudinal biases in retrieved O3 columns due to the use of hemispherically symmetric and time-independent profile adjustments, and O3 column discontinuities across the latitude zone boundaries due to large morphology differences. The deficiency of symmetric profile changes may be addressed with the newer column-dependent O3 climatologies with hemispheric separation (Lamsal et al.2004; Labow et al.2015), but their latitude dependencies are compiled on the same 30 wide latitude zones as the standard profiles, with a coarse (semi-annual) seasonal variation (Lamsal et al.2004) or without seasonal distinction (Labow et al.2015).

Since the TOMS-V8 climatology is widely used today by many different total O3 algorithms, there is a need to eliminate its deficiencies by creating a new O3 profile climatology. The primary objective of the investigation presented in this paper is to develop a new column-classified O3 climatology using recent data to provide realistic a priori O3 profile specifications in total O3 retrievals from the ultraviolet (UV) observations from the Earth Polychromatic Imaging Camera (EPIC) on board the Deep Space Climate Observatory (DSCOVR) satellite, the Ozone Mapping and Profiler Suite Nadir Mapper (OMPS-NM) on the Suomi National Polar Partnership (SNPP) and NOAA-20 satellites, the Ozone Monitoring Instrument (OMI) on Aura, or any other instruments on current and future satellites.

Profile O3 retrieval algorithms based on the optimal inversion method (Rodgers2000) need not only a priori O3 profiles but also the associated profile covariance matrices to constrain the retrieved profiles. However, O3 climatologies (e.g., McPeters et al.2007; McPeters and Labow2012; Bak et al.2013; Sofieva et al.2014) that are frequently employed in profile retrievals do not contain information about O3 profile covariance; instead the needed covariance matrices are constructed by assuming positive covariance between different atmospheric layers (e.g., Hoogen et al.1999; Liu et al.2005; Bhartia et al.2013; Miles et al.2015). This assumption is often unrealistic, as O3 profile changes resulting from atmospheric vertical motions usually have negative correlations among layers in the UTLS. Improved knowledge about the variability of O3 vertical distributions, especially how changes among different layers are related, benefit both profile and total column retrievals. Therefore, it is important to develop a new O3 profile climatology that provides quantitative information about O3 profile covariance.

3 MERRA-2 O3 field

The Modern-Era Retrospective Analysis for Research and Applications version 2 (MERRA-2, Bosilovich et al.2015; Gelaro et al.2017) reanalysis project produces assimilated data products, including meteorological fields (such as atmospheric temperature and tropopause pressure and altitude) and an O3 field. This global O3 field is driven by atmospheric dynamics and constrained by satellite O3 measurements and is continuous in time and three-dimensional space. Beginning in October 2004, Aura MLS provides constraints on the O3 profile from the lower mesosphere to the upper troposphere. Though MLS does not have information from the lower troposphere, Aura OMI and MLS jointly provide constraints on the tropospheric O3 column, with its vertical distribution controlled by atmospheric transport and simplified chemistry with parameterized O3 loss implemented in the MERRA-2 assimilation system (Wargan et al.2015; Bosilovich et al.2015). Evaluation of MERRA-2 O3 shows excellent agreement between the assimilated O3 profile in the stratosphere and upper troposphere (between 1 and 500 hPa) with independent satellite and ozonesonde measurements. It also shows that MERRA-2 assimilation realistically reproduces the stratospheric and upper-tropospheric O3 variability (Wargan et al.2015, 2017). In the lower troposphere, the correlations between MERRA-2 O3 profiles and the coincident ozonesonde measurements are lower than the correlations in the UTLS (Wargan et al.2015, 2017). But the tropospheric O3 columns (which are O3 amounts resulting from surface-to-tropopause profile integration) agree very well and have a high degree of correlation with the corresponding columns from ozonesonde measurements (Ziemke et al.2014), validating the integrated tropospheric profiles of MERRA-2. In short, a MERRA-2 assimilated O3 field provides a realistic representation of atmospheric O3 that faithfully captures both short-term (daily or shorter) variations and seasonal changes in vertical and horizontal distributions and thus contains vast information to be harvested for improved characterization of O3 vertical distributions and variations.

Figure 2MERRA-2 monthly mean O3 profiles for four months (January, April, July, and October) and seven tiles that cover the same longitude range (0–15 E) and seven different latitude zones (90–80 S, 50–40 S, 30–20 S, 00–10 N, 20–30 N, 40–50 N, and 80–90 N). Each profile (plotted as a colored line) is calculated from samples with the same UTC time, with different colors representing different UTC times shown in the legend. These mean profiles are averages over 15 longitude, which is equivalent to 1 h of local solar time (LST).


4 MERRA-2 O3 and temperature profile climatologies

The MERRA-2 record since the Aura MLS and OMI assimilation, owing to its realistic representation of atmospheric O3, is suitable for creating climatologies to provide improved knowledge of O3 and temperature vertical distributions.

The long-term MERRA-2 record is stored as a four-dimensional (latitude, longitude, atmospheric pressure, and time) data set that covers the globe uniformly at high spatial and temporal resolutions. More precisely, the global coverage is provided at 3-hourly intervals with a horizontal resolution of 0.5 latitude by 0.625 longitude and a vertical grid covering 72 layers between the surface and 0.01 hPa. The MERRA-2 O3 profiles from 2005 to 2016, within the period of Aura MLS and OMI assimilation, are analyzed to create a new O3 profile climatology. The enormous number of data from this period facilitate reliable statistical representation of mean O3 profiles and their variations on more finely resolved dependent variables, including tropopause pressure, total column amount, latitude, longitude, and time.

MERRA-2 data fields (such as O3 and temperature profiles) are represented on a hybrid terrain-following pressure coordinate (Wargan et al.2017), and thus have the lowest altitude level at the surface. These data fields are interpolated or extrapolated (for grids below the surface down to the sea level) to the uniform pressure altitude grid before statistical computation to create sea-level-based MERRA-2 climatologies. Here the pressure altitude, Z*=16log10[psp], is a dimensionless quantity, but is assigned units of kilometers because it is quite close to the altitude (km) at pressure level p above the surface at ps=1013.25 hPa, with a difference usually less than 1 km when Z*30 km.

4.1 Local-solar-time-dependent O3 profile climatology

The photochemical O3 production and destruction depend strongly on solar illumination, which changes systematically for a location during the course of a day, resulting in diurnal variation in O3 vertical profiles (Haefele et al.2008; Huang et al.2008; Sakazaki et al.2013; Schranz et al.2018), but existing O3 profile climatologies have not included O3 profile dependence on local solar time. Remote sensing O3 measurements are collected at different local solar times, even from the same instrument, such as the DSCOVR EPIC, which observes the Earth from sunrise to sunset simultaneously (Herman et al.2017). Proper accounting for the diurnal variations in O3 vertical distribution would enable more accurate O3 retrieval from remote sensing observations. Therefore we construct a local-solar-time-dependent O3 profile climatology from the MERRA-2 O3 field.

We divide the globe equally into 24×18 rectangle tiles, and each has the size of 15 longitude × 10 latitude. Monthly mean profiles and their variances are calculated from samples that fall within a tile at the eight UTC times of the MERRA-2 O3 field. Note that data samples are weighted by the cosine of latitude in performing the statistics to ensure equal weighting by area. Since each climatological profile of a tile is compiled from O3 profiles at the same UTC, it represents an hourly mean because a tile's 15 longitude range is equivalent to 1 h of local solar time (LST).

Figure 2 shows the sample O3 profiles from the MERRA-2 local-solar-time-dependent climatology, illustrating the diurnal, seasonal, and latitudinal variations in O3 vertical distributions. Each panel in Fig. 2 contains plots of eight climatological O3 mixing ratio profiles for a tile next to and east of the prime meridian, for eight different 1 h LST periods. The differences among the curves in a panel reveal the diurnal O3 mixing ratio changes, which in general increase significantly with altitude from the upper (Z*30 km) stratosphere into the mesosphere (Z*50 km). The O3 mixing ratio level in the mesosphere reaches its minimum after sunrise but then increases after sunset, recovering its nighttime value quickly, while in the stratosphere where the peak value of O3 mixing ratio is located (Z* between ∼35 and 40 km), O3 diurnal cycle has its maximum in the afternoon and a lower value close to its minimum in the morning. The diurnal variations, which are mainly driven by photochemistry in the upper atmosphere, exhibit a strong latitudinal (panel rows in Fig. 2) and seasonal (panel columns) dependence, as expected from the variations in solar insolation with the latitude and the season. The mean nighttime O3 mixing ratio can be an order of magnitude higher than the daytime value in the mesosphere at Z*=70 km. In general, the magnitude of diurnal difference reduces with the altitude, with a maximum of ∼7 % near the stratospheric O3 mixing ratio peak. In the troposphere, the cycle of diurnal variation changes with season and location, with a peak-to-trough difference usually less than 5 % in the lower troposphere. These significant diurnal O3 mixing ratio variations correspond to typically 0.5 % and a maximum of ∼4 % peak-to-trough differences in diurnal variation in total O3 vertical column. The variability of the tropospheric O3 profile, including the diurnal cycle, is mostly subdued in the MERRA-2 O3 field because the assimilation captures only the average behavior of tropospheric O3 but not its variation realistically, a consequence of the limited tropospheric O3 information contained in the source data (OMI columns and MLS profiles). The results for the mesosphere and stratosphere are in general consistent with previous findings of O3 diurnal variations (Haefele et al.2008; Huang et al.2008; Sakazaki et al.2013; Schranz et al.2018, and references therein). Thus, this local-solar-time-dependent climatology correctly captures diurnal variation in upper-atmospheric O3 vertical distributions.

Figure 3MERRA-2 daytime monthly zonal mean tropopause altitude and its standard deviation as a function of latitude for April (M04) and October (M10).


Figure 4Profile comparisons between M2TPO3 and TPO3 of Sofieva et al. (2014) for four months (January, April, July, and October) and seven latitude zones (90–80 S, 50–40 S, 30–20 S, 00–10 N, 20–30 N, 40–50 N, and 80–90 N). Colored solid lines represent M2TPO3 profiles, while the dotted ones represent TPO3 profiles. The color of a solid line indicates the percentage occurrence of the profile, which is calculated as the percentage of profiles in the month–latitude class that fall within the Z* tropopause pressure bin. The line legends display the average tropopause altitude and the average total O3 column for the corresponding climatological profile. The solid gray line represents the downgraded M2TPO3 profile, i.e., the monthly zonal mean profile.


Figure 5Similar to Fig. 4, but for standard deviation (SD) comparisons between M2TPO3 and TPO3 of Sofieva et al. (2014).


4.2 Tropopause-pressure-classified O3 profile (M2TPO3) climatology

Tropopause altitude varies with month and latitude systematically (Hoinka1998) and is usually highest in the tropics and drops toward the poles with steep declines in the midlatitudes but exhibits hemispheric asymmetry (e.g., see Fig. 3). This figure also illustrates that the dynamic variability of tropopause altitude is characterized by daily standard deviations of approximately 1–2 km, which agrees well with the characterization by Seidel and Randel (2006). Since O3 concentration is highest in the lower stratosphere above the tropopause and decreases rapidly in the troposphere, the large short-term fluctuation of tropopause altitude results in high variability in the total O3 column, accounting for ∼60 % of its daily variation (Wei et al.2010). Thus grouping O3 profiles according to the tropopause pressure or altitude generates climatological profiles that reflect the dynamical influences on O3 vertical distributions, as demonstrated by the tropopause-altitude-classified O3 profile climatology (named the TPO3 climatology) created by Sofieva et al. (2014) using ozonesonde and SAGE-II data.

Similarly, we create a tropopause-pressure-classified O3 profile climatology (referred to as M2TPO3 climatology hereafter) using the MERRA-2 O3 record, which provides a large number of data covering a diverse range of conditions, likely including all possible tropopause pressures. These mean profiles and their variances are calculated by statistically analyzing the sets of MERRA-2 O3 profiles at the same UTC time, binned by tropopause pressure in 1 km pressure altitude (Z*) steps, calendar month, and 15× 10 longitude–latitude tiles. To compare the M2TPO3 with the TPO3 (Sofieva et al.2014), we create a daytime zonal mean climatology by combining M2TPO3 tiles with the same latitude zone and LST from 09:00 through 17:00. The resulting daytime climatology contains 2154 mean profiles and standard deviations (shown in Appendix A), spanning the range of tropopause pressure (altitude) from 605 hPa (3.56 km) to 62 hPa (19.3 km), distributed in 12 calendar months and 18 10 latitude zones that cover the latitude range from −90 to 90.

Figures 4 and 5 display a subset of daytime M2TPO3 profiles and standard deviations and their comparisons with those of TPO3 (Sofieva et al.2014). The results in Fig. 4 show that both M2TPO3 and TPO3 have similar O3 profile shapes for similar tropopause altitudes in each month–latitude class. Especially in the tropics, M2TPO3 and TPO3 show very close profiles that are nearly independent of tropopause altitude and season. For profiles at higher latitudes, Fig. 4 illustrates that the profile change associated with tropopause altitude variation occurs mostly below the O3 concentration peak, exhibiting increasing UTLS O3 with lower tropopause altitude. The O3 column amounts (displayed in the figure legends) of TPO3 profiles are generally (about 1 % to 4 %) higher than the corresponding M2TPO3 profiles. This comparison indicates that there is likely a slight positive bias in the total columns in the TPO3 climatology since MERRA-2 total O3 columns have a latitude-dependent (up to 2 %) low biases (Wargan et al.2017). Comparing with the TPO3, there are more tropopause pressure bins with sufficient samples for reliable statistics, and hence more M2TPO3 profiles corresponding to a broader range of tropopause altitudes in each month–latitude class. Under the O3 hole condition, strong O3 spatial variations sampled differently by ozonesondes likely contribute to significant profile differences between M2TPO3 and TPO3 (see Fig. 4 lower left panel). Note that TpO3 is binned in 1 km tropopause altitude steps, and the corresponding legend shows the midpoint of the altitude bin. Consequently, the tropopause altitude in the legends of Fig. 5 changes regularly (in 1 km steps). On the other hand, M2TOP3 is binned in 1 km Z* tropopause pressure altitude steps, and the average tropopause altitude of the binned profiles is displayed in the corresponding legend. This average value does not change regularly from one bin to the next because the distribution of tropopause altitude is not necessarily symmetric with respect to the bin center.

The standard deviations in Fig. 5 show similar magnitudes and vertical structures between the M2TPO3 and the TPO3 profiles with similar tropopause altitudes, demonstrating that the M2TPO3 climatology represents O3 variability realistically. The main differences occur in the lower troposphere of tropical latitude zones, where TPO3 shows higher variability than the M2TPO3, due to ozonesondes better capturing the influence from tropospheric O3 productions, which are not explicitly included in the MERRA-2 reanalysis. Significant differences are also observed over polar latitude zones (90–80 S and 80 N–90 S), likely due to differences in sampling over highly inhomogeneous spatial distributions in these zones.

Results in Fig. 5 show that the standard deviations of the high-occurrence profiles are in general smaller than those for the downgraded (monthly zonal mean) profile, indicating that tropopause pressure is a good predictor of O3 profile shape, leading to a more accurate profile specification than simply taking the monthly zonal mean profile.

The comparison in this section shows overall good agreement of the M2TPO3 means and standard deviations versus those of the TPO3, with an expected discrepancy in the O3 variability in the tropical lower troposphere, validating the realism of MERRA-2 O3 record and its suitability for climatology construction.

Figure 6Profile comparisons between M2TCO3 and TOMS-V8 for four months and seven latitude zones (same as those in Fig. 4). Colored solid lines represent M2TCO3 profiles, while the dotted ones represent TOMS-V8 profiles. The color of a solid line indicates the percentage occurrence of the climatological profile, and its line legend displays the average tropopause altitude and the average total column O3. The solid gray line represents the downgraded M2TCO3 (monthly zonal mean) profile, which is the same as the downgraded M2TPO3 profile shown in Fig. 4. Note that since the total O3 distribution is usually not symmetric with the bin center, the average total O3 column does not change regularly from one bin to the next (see O3 columns in the legends). On the other hand, TOMS-V8 profiles are created to cover preset (latitude-zone-dependent) O3 column ranges in 50 DU steps.


Figure 7Standard deviation (SD) comparisons between M2TPO3 and M2TCO3 for four months and seven latitude zones (same as those in Fig. 4). Colored solid lines represent M2TCO3 SDs, while colored dotted lines represent M2TPO3 SDs. The color of a solid or dotted line indicates the percentage occurrence of the climatological profile, and its line legend displays the average tropopause altitude and the average total column O3. Only SDs for climatological profiles with greater than 5 % occurrence are plotted. The solid gray line represents the downgraded M2TCO3 profile SD, i.e., the SD for the monthly zonal mean profile, and it is identical to the SD of the downgraded M2TPO3 profile. The SDs quantify the natural variability of O3, and a smaller one signifies that the climatological mean provides a more reliable O3 profile specification.


4.3 Total-column-classified O3 profile (M2TCO3) climatology

The tropopause-pressure-dependent O3 profiles contained in the M2TPO3 climatology capture profile changes resulting from short-term meteorological disturbances in the UTLS region. The near-linear relationship between mean total column O3 and tropopause altitude (e.g., see legend tables in Fig. 4) quantified for each month–latitude class implies that the total column O3 may serve as a good predictor of O3 profile shape as well. Grouping profiles by total column O3 may actually better capture O3 distribution and variability driven by both dynamical and chemical processes; we therefore create a total-column-classified O3 profile climatology (referred to as the M2TCO3 climatology hereafter) from the long-term MERRA-2 O3 record.

The M2TCO3 climatology is constructed by binning the MERRA-2 O3 field at the same UTC time by total column in 25 DU steps, calendar month, and 15× 10 longitude–latitude tiles. To compare the M2TCO3 with the TOMS-V8 (Bhartia and Wellemeyer2002; McPeters et al.2007), we create a daytime zonal mean climatology by combining M2TCO3 tiles with the same latitude zone and LST from 00:09 through 17:00. The resulting climatology contains 1644 mean profiles and standard deviations (shown in Appendix A), spanning the range of total O3 from 94 to 584 DU, distributed in 12 calendar months and 18 latitude zones that cover the latitude range from −90 to 90. Figure 6 displays a subset of daytime M2TCO3 profiles and the corresponding TOMS-V8 profiles, illustrating the distinct characteristics of O3 vertical distributions, as well as significant differences between the two climatologies. Both M2TCO3 and TOMS-V8 profiles exhibit higher altitudes of O3 peak with lower O3 columns for a month–latitude class, lower peak altitudes with higher latitudes for a total O3 column, narrow dynamical O3 ranges and nearly identical profile shapes for different seasons in the tropics but higher seasonal and dynamical variations in midlatitude and high-latitude regions, and hemispherical asymmetry.

The differences between M2TCO3 and TOMS-V8 (shown in each panel of Fig. 6) reflect the improvements in the realism of O3 profile representation with the M2TCO3 climatology. The column-averaged TOMS-V8 profiles are represented by the LLM climatology (McPeters et al.2007), which is updated to the (McPeters and Labow2012) (ML) climatology, and are very close to the downgraded M2TCO3 profiles (see comparisons in Figs. 1 and 13). Therefore the differences between M2TCO3 and TOMS-V8 profiles shown in Fig. 6 are attributed to the profiles deviations applied to the mean derived from the standard column-dependent profiles of TOMS-V8 (Bhartia and Wellemeyer2002), which are independent of season, coarse in latitude resolution, and symmetric with respect to the Equator. These deficiencies limit the realism of TOMS-V8 climatological profiles, and are eliminated with the new M2TCO3 climatology, which provides improved O3 profile representation that better captures profile changes associated with column variations and their dependence on season and latitude.

The total O3 range of the TOMS-V8 climatology is the same as the set of TOMS-V8 standard (annual mean) profiles. These preset (latitude-zone-dependent) ranges may be different (narrower or broader) than those of M2TCO3 at some months and latitude zones (see Fig. 6), but they are likely inconsistent with actual O3 range due to the imposition of annual O3 range on the monthly variation.

Figure 7 shows the standard deviations of a subset of M2TCO3 profiles with high (>5 %) occurrence percentage and comparisons with those of M2TPO3. A vast majority of high-occurrence climatological profiles of M2TCO3 and M2TPO3 shown in Fig. 7 exhibit significantly reduced standard deviations compared to those of monthly zonal mean (i.e., downgraded) profiles, illustrating that both tropopause pressure and total column O3 provide information for more precise specification of O3 profiles. Comparisons between M2TPO3 and M2TCO3 show that column classification usually leads to greater reductions in standard deviations in the upper stratosphere but smaller ones in the UTLS than tropopause pressure classification. Since there is much more O3 in the upper stratosphere than below, an overall more realistic specification of the O3 profile is achieved using the M2TCO3 climatology based on the total column O3 than using the M2TPO3 based on the tropopause pressure. Hence the column-dependent climatology is most appropriate in mapping between the O3 column and vertical profile needed in total O3 column retrieval algorithms. However, with the knowledge of tropopause pressure, the tropopause-dependent climatology usually provides closer matches to actual profiles and stronger constraints in the UTLS; therefore its usage significantly improves the accuracy of O3 profile retrieval (Bak et al.2013) from backscattered UV measurements, which have lower vertical resolutions in the troposphere than in the stratosphere (Liu et al.2010).

Figure 8Correlation matrices from M2TPO3 (a) and M2TCO3 (b) for four months and seven latitude zones (same as those in Fig. 4). Here the matrix for the highest-occurrence climatological profile is shown for each month–latitude class.


Figure 9 


4.4O3 profile covariance climatology

Previous O3 climatologies (e.g., Fortuin and Kelder1998; McPeters et al.2007; McPeters and Labow2012; Bak et al.2013; Sofieva et al.2014) used in O3 profile retrievals include profile standard deviations but not information about O3 profile covariance because their sources (ozonesonde and satellite measurements) have limited coincident samples to quantify the joint O3 variability between different altitudes throughout the atmosphere. In contrast, the MERRA-2 assimilation provides complete profiles simultaneously, allowing direct statistical computation of profile covariance matrices. We expand the classification analysis to quantify the O3 profile covariance. Specifically, we construct an O3 profile covariance matrix from each bin used for the calculation of an M2TPO3 or M2TCO3 climatological profile. The resulting covariance climatology first provides quantification of O3 vertical distribution variabilities, the correlations among different levels, and their dependence on tropopause pressure or total column O3 for different months and longitude–latitude tiles.

Figure 8 shows a subset of correlation matrices, which are standardized, i.e, diagonal element normalized to covariance matrices, from the daytime M2TPO3 and daytime M2TCO3 climatologies. These density plots highlight the varying degree of level-to-level correlations and their contrasts with the diagonal-constant matrices typically used in O3 profile retrievals (e.g., Hoogen et al.1999; Liu et al.2005; Miles et al.2015), which assume positive layer-to-layer correlation that decreases monotonically and exponentially with distance between layers. Contrary to this typical assumption, Fig. 8 illustrates that the degree of correlation between two levels fluctuates with the level separation and negative correlations are quite common.

The quantification of O3 profile covariance offers a new way to represent O3 vertical distribution realistically and efficiently. Figure  shows the first three leading empirical orthogonal functions (EOFs), which are the ordered set of eigenvectors for each covariance matrix shown in Fig. 8. Typically the first three leading EOFs combined explain between 65 % and 85 % of the class variance, and the first 15 leading EOFs account for over 95 % of the class variance. Much like a climatological profile describes the likely O3 vertical distribution, the EOFs describe the most probable patterns of profile deviations from the climatological mean. In general, an O3 profile X can be expressed as

(1) X = X m + k = 1 n ω k e k ,

where Xm is a climatological mean profile, ek the kth EOF, ωk the kth coefficient, and n the total number of EOFs with its maximum limited to the number of levels used to represent the O3 profile. Since a high fraction of class variance may be explained with just a few leading EOFs, they provide the most efficient adjustments to improve the representation of O3 profile when it deviates from the climatological mean. When interpreted physically, the leading EOFs correspond to processes of O3 accumulation, reduction, or redistribution. For instance, the first three EOFs for an M2TPO3 covariance matrix (see Fig. a) describe (1) column increase or decrease, represented by the EOF-1 (monopole) pattern, (2) up or down shift of a profile, represented by the EOF-2 (dipole) pattern, showing O3 increase at one altitude and decrease at another, and (3) shrink or stretch of a profile, represented by the EOF-3 (tripole) pattern, showing O3 decrease (increase) at one altitude and increase (decrease) at adjacent (above and below) altitudes. The results in Fig. a indicate that the highest or the second highest contribution to the variance of a tropopause class originates from the total column variation. However, the highest contribution to the variance of a total column class comes from O3 profile shift, represented by EOF-1 (dipole) in Fig. b, usually linked with tropopause variation, since the dominant profile change resulting from column variation is accounted for with different O3 column classes. The second EOF (tripole) of an M2TCO3 covariance matrix describes the shrink or stretch of a profile, similar to the third EOF of an M2TPO3 matrix, while subsequent EOFs describe more complex rearrangement of the O3 profile.

Figure 10Climatological temperature profiles corresponding to the M2TCO3 O3 profiles shown in Fig. 6. The color of a solid line indicates the percentage occurrence of the climatological profile, and its line legend displays the average tropopause altitude and the average total column O3. The solid gray line represents the monthly zonal mean temperature profile, and the dashed dark blue line shows the coefficient of correlation between O3 partial pressure and temperature as a function of pressure altitude Z*.


4.5 Temperature profile climatology

Using the same binning schemes employed in the O3 profile climatologies, we create temperature profile climatologies from the MERRA-2 assimilated atmospheric temperature field. The resulting climatologies contain the mean and covariance of temperature profiles, paired with the respective M2TPO3 and M2TCO3 climatological O3 profiles.

Figure  shows sample MERRA-2 climatological temperature profiles corresponding to the sample M2TCO3 climatological O3 profiles shown in Fig. 6, and the O3–temperature correlation profiles for the sample month–latitude classes. These results illustrate the systematic behavior among O3 column amount, tropopause altitude, and atmospheric temperature: higher O3 column amounts occur with a warmer lower stratosphere and a colder troposphere, the condition for lower tropopause altitudes. This relationship, as well as its dependence on season and latitude, is captured in the MERRA-2 O3 and temperature climatologies and is consistent with the findings of long-term O3 and temperature profile measurements (Steinbrecht et al.1998).

Results in Fig.  also reveal a significant north–south asymmetry in climatological temperature profiles. For instance, the stratospheric temperatures in the austral spring (e.g., see the October panel) in the southern polar region (90–80 S) are significantly lower than those in the boreal spring (e.g., see the April panel) in the northern polar region (80–90 N), and these temperature differences are associated with lower total O3 columns in the Antarctic than those in the Arctic.

The profile of correlation coefficient in each panel of Fig.  illustrates a mutual relationship between O3 concentration and atmospheric temperature. From the upper (Z*35 km) stratosphere to the lower mesosphere (Z*60 km), O3 and temperature are mainly negatively correlated because in this region O3 concentration is mostly governed by photochemical reactions, for which higher temperature speeds up the rate of O3 destruction. From the lower stratosphere down to the upper troposphere (10 km Z* 35 km) in which O3 concentration is controlled primarily by atmospheric motions (Brasseur and Solomon2005), O3 and temperature are positively correlated since higher O3 concentrations likely occur from adiabatic air parcel compression, which increases the parcel temperature as well, and additionally higher O3 concentrations absorb more UV radiation, thus raising the atmospheric temperature. The positive correlation quickly becomes negative as the altitude descends into the lower troposphere (Z* 10 km), but swings back and may become positive as the altitude falls further. In general, the degree of O3–temperature correlation is smaller in the lower troposphere than in the atmosphere above, indicating a weaker connection between them.

The pattern of negative correlation above the upper stratosphere and the positive correlation in the upper troposphere and lower stratosphere between O3 and temperature fields were elucidated with dynamical chemical models (Rood and Douglass1985; Froidevaux et al.1989; Smith1995). The O3–temperature correlation profiles in Fig.  are consistent with those from numerical modeling (Rood and Douglass1985) and observational data analysis (Fortuin and Kelder1996), indicating that the MERRA-2 O3 and temperature climatologies represent O3 and temperature distributions and their interrelationship realistically.

4.6 Spatial distribution of O3

Figure 11Climatological stratospheric column O3 (SCO) maps from MERRA-2 and Ziemke et al. (2011) for the 12 months (January–June, upper panels; July–December, lower panels) of the year and their differences.

Figure 12Similar to Fig. 10, except for tropospheric column O3 (TCO).

Figure 13MERRA-2 climatological O3 concentration (i.e., mixing ratio in parts per million by volume) at Z*=4.5 km for the 12 months (January–June, upper panels; July–December, lower panels) of the year and comparison with the tropospheric O3 profile climatology of G. Liu et al. (2013). Note that the mixing ratio of G. Liu et al. (2013) is at 4.5 km altitude.

The spatial distribution of O3 is controlled by the various chemical and dynamical processes that drive the production, destruction, and transport of atmospheric O3. Since the distributions of O3 sources and surface topography are inhomogeneous over the globe, systematic patterns with significant longitudinal variations are present in horizontal O3 distribution. The MERRA-2 baseline climatologies presented in previous sections focus on the latitudinal and altitudinal distribution of O3, with a coarse longitudinal resolution retained through the binning of the globe with equal size (15 longitude × 10 latitude) tiles. To better capture the spatial variation in O3 and to compare directly with previous global O3 climatologies (Ziemke et al.2011; Liu et al.2013, referred to respectively as Ziemke2011 and GLiu2013 hereafter), we create a higher-spatial-resolution O3 profile climatology from MERRA-2 by reducing the tile bin size to 5 longitude × 5 latitude, which is the bin size for both Ziemke2011 and GLiu2013 climatologies. This higher-spatial-resolution MERRA-2 climatology consists of monthly statistics of O3 volume mixing ratio profile, tropopause pressure, stratospheric column O3 (SCO), and tropospheric column O3 (TCO) for 2592 (=36×72) 5× 5 tiles. Here the SCO is an integration of an O3 profile from 0.01 hPa down to the tropopause and the TCO from the tropopause down to the surface.

For better comparisons with the Ziemke2011 climatology, we partition the MERRA-2 total O3 column into SCO and TCO using the climatological tropopause pressure provided in Ziemke2011, which is compiled from the National Centers for Environmental Prediction (NCEP) tropopause pressures. Figures 10 and 11 show the spatial distributions of the NCEP-based MERRA-2 SCO and TCO, respectively, for 12 months of the year and their comparisons with Ziemke2011. Results in Figs. 10 and 11 reveal excellent agreement between MERRA-2 and Ziemke2011 in the low-latitude zone (within ±30), with ΔSCO ±5 DU and ΔTCO ±6 DU for most months and tiles. This agreement becomes worse in the higher-latitude regions, mostly showing larger differences in SCO and TCO between MERRA-2 and Ziemke2011, likely due to higher longitudinal variability in stratospheric O3, degrading the accuracy of SCO from gap-filling interpolation of Aura MLS data (Ziemke et al.2011). From September to May, MERRA-2 SCO in the polar latitudes is higher, while TOC is lower than the corresponding columns from Ziemke2011, with the broadest spread occurring in February–April. But more significantly, Figs. 10 and 11 show that both spatial distributions of MERRA-2 SCO and TCO and their seasonal cycles closely resemble those of Ziemke2011.

MERRA-2 SCO displays a strong latitude dependence (see Fig. 10) that is shaped primarily by the Brewer–Dobson circulation, showing low values in the tropics (∼220 DU) and elevated values in the middle and high latitudes, with the highest column amount exceeding 400 DU in the Northern Hemisphere (NH) (between 70 and 80 N) during February–April and exceeding 350 DU in the Southern Hemisphere (SH) (between 45 and 60 S) in September–October. The SCO longitudinal variability is low in the tropics but increases significantly along with higher columns at midlatitudes and high latitudes. In the midlatitude to high-latitude regions centered around 60 N, the longitudinal variation exhibits an oscillatory pattern with high SCO over North America that stretches across the Pacific Ocean to over eastern Asia. This pattern, which is influenced by the semipermanent atmospheric pressure systems resulting from the unique orography distribution in the NH, persists over the year with amplitude varying with the season and reaching its maximum during February. In the SH, a wavelike high-SCO pattern with longitude center near the dateline (the 180th meridian) also appears in the middle–high-latitude band (around 60 S), which evolves with the progression of the Antarctic O3 hole from formation to breakup, exhibiting SCO buildup and decay outside the stratospheric polar vortex and reaching its maximum in September–October when the O3 hole drops to its minimum (∼125 DU).

As illustrated in Figs. 10 and 11, TCO behaves differently from SCO, reflecting different sources and dynamical processes that affect its spatiotemporal distribution. The NH mean SCO rises in January, reaches its maximum during February–April, drops in May, and reaches its minimum in September. By comparison, the NH mean TCO starts to increase in March, maximizes in the summer months, and begins to drop in the fall, arriving at its minimum in January–February. These different seasonal cycles indicate a weak relationship between TCO and SCO in NH.

In the tropics, TCO is small (15 to 25 DU) in the Pacific but high (35 to 45 DU) in the Atlantic, exhibiting a stable wavelike pattern that detaches from the SCO distribution, which is almost longitudinally invariant in the same latitude zone. The small TCO in the tropical Pacific is likely due to atmospheric deep convection that uplifts marine boundary air, which is O3 poor, into the middle and upper troposphere. The elevated TCO in the tropical South Atlantic and the connected high-TCO band along the edge of the southern tropics (30 S) between 40 W and 80 E are contributed by regional O3 production, including lightning, biomass burning, fossil fuel combustion, and soil emissions, as well as influx from the stratosphere (Thompson et al.2000; Lelieveld and Dentener2000; Moxim and Levy2000; Martin et al.2002; Edwards et al.2003; Sauvage et al.2007). This enhancement persists throughout the year, but it is strongest in September–November and weakest in March–May. Its longitudinal variation and seasonal dependence are influenced by large-scale atmospheric circulations (Wang et al.2006). In the NH extratropics (between 20 and 40 N), TCO displays a zonal band elevation with strong longitudinal variability and the highest TCO occurs over the Mediterranean and eastern Asia in June–July. The TCO enhancement in this industrialized zonal band is significantly contributed from anthropogenic emissions (Li et al.2001; Liu et al.2011; Jiang et al.2016).

Aura MLS and OMI O3 measurements are combined to create the Ziemke2011 climatology, and they are assimilated to generate the MERRA-2 O3 field. Consequently, the close resemblance described in this section between MERRA-2 and Zimeke2011 climatologies is expected, even though MERRA-2 uses a period that is twice as long as that of Zimeke2011. To further evaluate the MERRA-2 climatology, we compare it with the GLiu2013 climatology, which is created from trajectory mapping of the long-term ozonesonde record (Liu et al.2013). Figure 12 shows monthly maps of O3 concentrations (i.e., mixing ratio in parts per million by volume) at 4.5 km altitude from GLiu2013 and at 4.5 km pressure altitude from MERRA-2 for 12 months of the year. These maps in Fig. 12 show similar O3 spatial distributions and temporal evolution between the two climatologies. For instance both exhibit persistent O3 concentration enhancements in the NH extratropics and in the tropical and subtropical South Atlantic. However, while GLiu2013 enhancements are located in roughly the same area as MERRA-2, they differ in shape, likely due to the spatial gaps of ozonesonde data. The seasonal cycles agree well with each other: both show the strongest NH enhancements in June–August and the weakest in January–March and the strongest SH enhancements in September–November and the weakest in March–May.

The close resemblances of O3 columns (both TCO and SCO) between MERRA-2 and Ziemke2011 and general similarities of O3 concentrations between MERRA-2 and GLiu2013 in terms of their spatial distributions and seasonal cycles demonstrate that the MERRA-2 climatology captures the spatiotemporal O3 distribution realistically.

5 Validation

Figure 14Annual zonal mean profiles and standard deviations from M2TCO3 and ML climatologies and their differences. (a, b, c) Latitude-dependent O3 partial pressure (mPa) profiles: M2TCO3 (a), ML (b), and their relative differences (c). (d, e, f) Latitude-dependent O3 profile standard deviations (mPa): M2TCO3 (a), ML (b), and their percent differences (c).


The baseline climatologies constructed from MERRA-2 data describe the systematic behavior of the O3 profile and variance and their spatial and temporal dependence on tropopause pressure and total O3 column. Since there is no other O3 profile climatology with a similar resolution and dependency, we validated the MERRA-2 climatologies by downgrading and then comparing them with independent climatological data sets. In Sect. 4.2, we validated the daytime tropopause-dependent (downgraded in LST and longitude) M2TPO3 climatology, which shows good agreement with the TPO3 climatology (Sofieva et al.2014) compiled from independent O3 profile data. In Sect. 4.6, we validated the spatiotemporal O3 variations represented in the downgraded (in LST and O3 columns) daytime M2TCO3 climatology, which shows good agreement with the Ziemke2011 and GLiu2013 climatologies. In this section, we present further comparisons of vertical profiles and integrated quantities (i.e., vertical columns) to demonstrate the validity of the M2TCO3 climatology.

We compare the daytime annual zonal mean (i.e., downgraded in LST, season, longitude, and O3 column) M2TCO3 climatology with the annual zonal mean (downgraded in season) ML climatology (McPeters and Labow2012), which is an improved version of the LLM climatology (McPeters et al.2007). Figure 13 shows the annual mean O3 profiles (panels a and b) and their standard deviations (panels d and e) as functions of altitude and latitude from the daytime M2TCO3 and ML climatologies. This figure also includes the plots of the relative (percent) difference between the mean profiles (panel c) and between the standard deviations (panel f) from these two climatologies. In the upper stratosphere (Z*>20 km), annual mean profiles have an excellent agreement between these two climatologies, with relative differences mostly within ±3 % (see Fig. 13a–c). However, relative differences are larger but within ±30 % in the lower stratosphere and troposphere (Z*<20 km), mainly due to sampling differences and due to the lack of tropospheric O3 production in the MERRA-2 reanalysis (Wargan et al.2015; Bosilovich et al.2015). These differences between M2TCO3 and ML climatologies are within the uncertainties estimated from the differences between two climatologies compiled from different sources, including the ML and the LLM in McPeters and Labow (2012), the TPO3 and the ML, and the TPO3 and the LLM in Sofieva et al. (2014).

As shown in Fig. 13d–f, the standard deviations of the downgraded M2TCO3 climatology are similar to those of the ML climatology: both capture the high variability in the upper troposphere and lower stratosphere (UTLS) and the low variability in the upper stratosphere, with a vast majority of differences within ±30 %. The sampling differences contribute to the standard deviation differences exhibited in Fig. 13f.

Figure 15Comparison of monthly zonal mean total column O3 from the daytime M2TCO3 (a) and the ML (b) climatologies and their relative differences (c). Both climatologies capture the strong annual cycle of total O3, which exhibits a maximum in late winter–early spring in the NH and a minimum in late summer–early fall in the SH.


Figure 14 shows the total O3 columns versus month and latitude from the ML and the downgraded M2TCO3 climatologies and their percent differences, illustrating the good agreements between the two for a vast majority of months and latitude zones, with O3 absolute differences less than 4 %. Absolute differences exceeding 4 % occur in a few areas only, notably in the northern polar zone for several months during which M2TCO3 total columns are lower than those of ML by more than 4 %, and in the SH above 70 latitude for the months from August to December during which M2TCO3 total columns are higher than those of ML by more than 4 %. These large biases likely result from spatially inhomogeneous O3 distributions, such as the longitudinally and seasonally dependent tropospheric O3 distributions (Ziemke et al.2011; Liu et al.2013) and the strong O3 gradient across the Antarctic ozone hole, which are sampled differently in creating the two climatologies. While the MERRA-2 assimilation provides a spatially and temporally uniform O3 field, the ozonesonde stations distributed unevenly around the globe provide intermittent O3 profile measurements, and the MLS on the polar-orbiting Aura platform samples more densely at higher latitudes. As a result, the ML climatology, which relies on the ozonesonde data and the Aura MLS O3 profile measurements, samples each latitude zone unevenly, thus contributing to the differences in Fig. 14.

In summary, comparisons of climatological profiles and standard deviations between M2TPO3 and TPO3 (see Figs. 4 and 5), annual zonal mean profiles and standard deviations between M2TCO3 and ML (see Fig. 13), and integrated monthly zonal mean profiles (i.e., the total column O3) between M2TCO3 and ML (see Fig. 14) show good overall agreement, with differences similar to those between two climatologies constructed from different data sources, thus validating the daytime M2TPO3 and daytime M2TCO3 climatologies.

6 Conclusions

Tropopause-pressure-classified (M2TPO3) and total-ozone-column-classified (M2TCO3) climatologies are created from the MERRA-2 O3 profile record between 2005 and 2016, within the period of Aura MLS and OMI assimilation. The enormous number of MERRA-2 O3 profile data that cover the globe uniformly and continuously enable precise and accurate representations of systematic behaviors and dynamical variations in O3 vertical distributions. The resulting set of O3 profiles and covariance matrices captures their dependence on longitude, latitude, local solar time, and season, as well as on tropopause pressure or O3 abundance, more accurately over a broader range and at a higher resolution than other O3 profile climatologies. Parameterization of O3 profile with tropopause pressure or total column O3 reduces the variability in stratosphere and troposphere compared to the month–latitude-dependent climatology, therefore providing improved a priori knowledge of O3 vertical distribution. Both M2TPO3 and M2TCO3 climatologies contain quantitative information about O3 profile covariances, which is not included in previous O3 profile climatologies. The profile covariances provide more realistic constraints on O3 profile retrievals based on the OE inversion technique. Moreover, the EOFs of the climatological covariance matrices facilitate a new scheme to represent the O3 profile and guide a retrieval algorithm to successively improve O3 profile user information contained in spectral measurements.

For profile retrieval algorithms, a closer match between actual and a priori O3 profiles, especially in the region where spectral measurements have low vertical resolution, improves the retrieval accuracy. Thus tropopause-dependent O3 climatology, which reduces the variability further in the UTLS region, is more appropriate for use with O3 profile algorithms. However, the variability reduction is overall higher with column O3 parameterization, indicating more realistic O3 profile assignment based on total column abundance. Therefore the column-dependent O3 climatology is uniquely suited for use in total O3 retrieval algorithms, as the retrieved column determines the likely O3 vertical distribution without needing additional information.

The M2TCO3 climatology provides improved O3 profile representation, capturing systematic profile changes resulting from column variations and their dependence on season and spatial location, which are missing from or insufficiently represented by previous column-classified climatologies (e.g., Wellemeyer et al.1997; Bhartia and Wellemeyer2002; Lamsal et al.2004; Labow et al.2015) with coarse latitude and time resolutions. The smooth profile change between adjacent dependent variables (illustrated in figures in Appendix A) implies that merging profiles from different columns, months, and tiles can provide spatially and temporally continuous representation of the O3 profile. The MERRA-2 temperature climatology and the M2TCO3 climatology are used in the O3 and SO2 combination algorithm applied to retrievals from DSCOVR EPIC (EPIC Science Team2018), SNPP OMPS-NM (Yang2017), and Aura OMI. The description and validation of these O3 and SO2 products will be presented in separate papers.

Data availability

The MERRA-2 climatologies are available by contacting the author (

Appendix A: Figures of the MERRA-2 O3 and temperature profile climatologies

Figure A1The daytime M2TPO3 climatology contains 2154 O3 profiles that are distributed among the 12×18 month–latitude classes. The color of a solid line indicates the percentage occurrence of the profile. The line legends display the average tropopause altitudes and the average total O3 columns. The solid gray line represents the downgraded M2TPO3 profile, i.e., the monthly zonal mean profile.


Figure A2Similar to Fig. A1, except for profile standard deviations in daytime M2TPO3 climatology. This illustrates that tropopause pressure classification in general reduces the variability of the climatological profile.


Figure A3The daytime M2TCO3 climatology contains 1644 O3 profiles that are distributed among the 12×18 month–latitude classes. The color of a solid line indicates the percentage occurrence of the profile. The line legends display the average tropopause altitudes and the average total O3 columns. The solid gray line represents the downgraded M2TCO3 profile, i.e., the monthly zonal mean profile.


Figure A4Similar to Fig. A3, except for profile standard deviations in daytime M2TCO3 climatology. This illustrates that total O3 column classification in general reduces the variability of the climatological profile.


Figure A5O3 profile correlation matrices corresponding to daytime monthly zonal mean profiles. The correlation matrices are standardized or normalized covariance matrices (with 1 s in the main diagonal).


Figure A6Climatological temperature profiles corresponding to the daytime M2TCO3 O3 profiles shown in Fig. A3. The color scheme is the same as that in Fig. A3. The dashed dark blue line represents the coefficient of correlation between O3 partial pressure and temperature as a function of pressure altitude Z*.


The baseline climatologies consist of statistics of MERRA-2 O3 and temperature profiles collected from 24×18 rectangular tiles at eight different UTC times each separated by 3 h. Since each tile has the size of 15 longitude by 10 latitude, the tile statistics at a UTC represent the statistics of an hour of local solar time. Combining the statistics of different tiles within a range of local solar time downgrades the baseline climatologies to month–latitude (i.e., no distinction in longitude) climatologies, as illustrated in this paper with the example of daytime (09:00–17:00) climatologies. Similarly morning, afternoon, or nighttime climatologies may be created by selecting the proper local solar time range in combining tiles with the same latitude zone.

In this appendix, each O3 profile in the daytime M2TPO3 and the daytime M2TCO3 climatologies and the corresponding profile variance are displayed in Figs. A1 to A4, which show O3 partial pressure or standard deviation (mPa) as a function of pressure altitude Z* from 0 to 71 km (or 1013.25 to 0.04 hPa). Colored solid lines represent the climatological profiles in Figs. A1 and A3 or the corresponding standard deviations in Figs. A2 and A4. The color of a solid line indicates the percentage occurrence of the climatological profile, and the line legend displays the average tropopause altitude and the average total column O3. The solid gray line represents the downgraded (monthly zonal mean) profile in Figs. A1 and A3 or the corresponding standard deviations in Figs. A2 and A4. Panels in each row show change with latitude, while those in each column reflect seasonal variation.

The MERRA-2 climatologies contain quantitative characterizations of O3 profile covariances. Correlation matrices, which are normalized covariance matrices, associated with daytime monthly zonal means are shown in Fig. A5.

Accompanied by the O3 profile climatologies, temperature profile climatologies are created. Figure A6 shows the climatological temperature profiles associated with the daytime M2TCO3 climatology. Figure A6 includes the monthly zonal O3–temperature correlation profiles to illustrate the consistent correlation pattern for all latitude zones and seasons.

Appendix B: Climatology usage


The MERRA-2 baseline climatologies provide local-solar-time-dependent coverage over the globe with equally sized (15 longitude × 10 latitude) tiles, which are quite coarse compared to the spatial resolutions of most spaceborne observations. For retrieval applications, the a priori information for a specific time and location may be determined from climatological data via spatial and temporal interpolation to ensure continuous representation in time and space. Frequently, the baseline climatologies are downgraded in longitude to generate the latitude-zone-dependent climatologies (e.g., see Figs. A1 and A3). A retrieval application using these latitude-zone-dependent climatologies usually performs temporal and latitudinal interpolation to provide time- and location-dependent a priori information.

Setting a climatological profile to a user-specified total column

Frequently, an O3 retrieval algorithm needs to set a climatological profile X to a specific total column Ω. For the M2TPO3 climatology, this is accomplished by finding the coefficient ω1, such that the profile integration of

(B1) X = X m + ω 1 e 1

is equal to Ω. For the M2TCO3 climatology, this can be done by linear interpolation or extrapolation of the column-dependent climatological profile Xm1) and Xm2),

(B2) X = w 1 X m ( Ω 1 ) + w 2 X m ( Ω 2 ) ,

where Ω1 and Ω2 are the two climatological columns that are closest to Ω, and the coefficients w1=(Ω2-Ω)/(Ω2-Ω1) and w2=(Ω-Ω1)/(Ω2-Ω1). When Ω is too far (−15 DU > or >15 DU) outside the valid M2TCO3 range, Eq. (B1) with Xm and e1 from the M2TCO3 climatology may be used to set the profile X to the desired total column Ω.

The O3 column used in selecting a profile from the M2TCO3 climatology refers to the total column from profile integration down to the sea level (Z*=0 km). When the bottom level of a user-defined vertical grid is significantly above the sea level, such as the cases over the Himalayan Plateau or Antarctica, the column from the integration of the user grid needs to add the climatological column from the bottom level to the sea level for profile selection from the M2TCO3 climatology.

Mapping of a climatology to a user-defined vertical grid

The baseline climatologies contain O3 concentration (i.e., mixing ratio) profiles specified at m=72 equally spaced pressure altitude levels (Z*=0 to 71 km by 1 km) and their covariance matrices. Frequently an O3 retrieval application needs to recast the concentration profiles and their covariance matrices into those of column amount of the vertical layers. Next, we describe the equations for mapping a level-based climatology to user-defined atmospheric layers.

Given an atmospheric layering scheme, it is straightforward to convert a mixing ratio profile X into a column density profile Ψ of n layers {Ψj=XjAj,j=1,,n}, where Xj is the mean mixing ratio of the jth layer, and Aj is the layer air column density, which is determined by the difference between pressures at the layer boundaries. Usually, Xj is determined via interpolation or extrapolation (if a user grid is outside the grid range of the climatology) from X.

An m×m covariance matrix Sm, which is real and symmetric by definition, may be decomposed exactly as

(B3) S m = Q Λ Q T ,

where Q=[e1,e2,,em] is an orthogonal matrix of m columns, with its ith column ei being the ith eigenvector of the covariance matrix Sm, and Λ=diag(λ1,λ2,,λm) is a diagonal matrix, with its ith element λi being the ith eigenvalue of Sm. Similar to converting a mixing ratio profile into a layer column density profile, the eigenvectors {ei,i=1,,m}, also known as the empirical orthogonal functions (EOFs), are converted into layer column density vectors, {ei,i=1,,m}, and the elements of the ith column vector are {eij=eijAj,j=1,,n}, where eij is the mean value of ei at layer j, usually obtained through interpolation or extrapolation from ei. The n×n layer-to-layer covariance matrix Sn is constructed as

(B4) S n = Q Λ Q T ,

where Q=[e1,e2,,em] is a matrix of m column vectors, each with the length of n.

A covariance matrix example

Figure B1(a) The covariance matrix for October and the midlatitude (50–40 S) zone, (b) eigenvalues of the covariance matrix in (a), and (c) the first nine leading EOFs of the matrix in (a).


We show an example of a covariance matrix Sm in Fig. B1a, its eigenvalues in Fig. B1b, and the nine leading EOFs in Fig. B1c. The eigenvalues drop rapidly with higher index i (see Fig. B1b) and are typically between 15 and 20 leading EOFs to account for 99 % of the total variances. The Sn may be reconstructed with a smaller number (<m) of leading EOFs in Eq. (B2). Doing so may reduce the total variance represented by Sn, thus improving the numerical stability of O3 profile retrievals that use Sn as an a priori constraint.

Author contributions

KY (the lead author) developed the overall design of this research, constructed the MERRA-2 climatologies, conducted analyses, performed validations, and prepared the paper. XL (coauthor) contributed to the development and presentation of diurnal variation in the O3 profile, tropopause-dependent climatology, ozone–temperature correlation, and temperature profile climatology. XL also contributed to the design and usage of the Appendix section: mapping of a climatology to a user-defined vertical grid.

Competing interests

The authors declare that they have no conflict of interest.


The MERRA-2 data used in this study are provided by the Global Modeling and Assimilation Office (GMAO) at NASA Goddard Space Flight Center and are available at the NASA Goddard Earth Sciences (GES) Data and Information Services Center (DISC). This work is supported by NASA.

Financial support

This research has been supported by NASA (grant nos. NNX15AB10G, NNX14AR20A, and NNX17AF56G).

Review statement

This paper was edited by Mark Weber and reviewed by two anonymous referees.


Bak, J., Liu, X., Wei, J. C., Pan, L. L., Chance, K., and Kim, J. H.: Improvement of OMI ozone profile retrievals in the upper troposphere and lower stratosphere by the use of a tropopause-based ozone profile climatology, Atmos. Meas. Tech., 6, 2239–2254,, 2013. a, b, c, d, e, f

Bhartia, P. K. and Wellemeyer, C. G.: TOMS-V8 Total O3 Algorithm, in: OMI Algorithm Theoretical Basis Document, edited by: Bhartia, P. K., vol. II, chap. 2, pp. 15–32, NASA Goddard Space Flight Center, Greenbelt, Maryland, USA, 2 edn., available at: (last access: 18 June 2019), 2002. a, b, c, d, e, f, g, h, i

Bhartia, P. K., McPeters, R. D., Flynn, L. E., Taylor, S., Kramarova, N. A., Frith, S., Fisher, B., and DeLand, M.: Solar Backscatter UV (SBUV) total ozone and profile algorithm, Atmos. Meas. Tech., 6, 2533–2548,, 2013. a, b

Bosilovich, M., Akella, S., Coy, L., Cullather, R., Draper, C., Gelaro, R., Kovach, R., Liu, Q., Molod, A., Norris, P., Wargan, K., Chao, W., Reichle, R., Takacs, L., Vikhliaev, Y., Bloom, S., Collow, A., Firth, S., Labow, G., Partyka, G., Pawson, S., Reale, O., Schubert, S. D., and Suarez, M.: MERRA-2 : Initial Evaluation of the Climate, in: NASA Technical Report Series on Global Modeling and Data Assimilation NASA/TM-2015-104606, edited by: Randal D. Koster, vol. 43, p. 139, Goddard Space Flight Center, Greenbelt, Maryland, USA, available at: (last access: 2 April 2019), 2015. a, b, c

Brasseur, G. P. and Solomon, S.: Aeronomy of the Middle Atmosphere, vol. 32 of Atmospheric and Oceanographic Sciences Library, Springer-Verlag, Berlin/Heidelberg,, 2005. a

Cariolle, D., Amodei, M., and Simon, P.: Dynamics and the ozone distribution in the winter stratosphere: modelling the inter-hemispheric differences, J. Atmos. Terr. Phys., 54, 627–640,, 1992. a

Coldewey-Egbers, M., Weber, M., Lamsal, L. N., de Beek, R., Buchwitz, M., and Burrows, J. P.: Total ozone retrieval from GOME UV spectral data using the weighting function DOAS approach, Atmos. Chem. Phys., 5, 1015–1025,, 2005. a

Edwards, D. P., Lamarque, J., Attié, J., Emmons, L. K., Richter, A., Cammas, J., Gille, J. C., Francis, G. L., Deeter, M. N., Warner, J., Ziskin, D. C., Lyjak, L. V., Drummond, J. R., and Burrows, J. P.: Tropospheric ozone over the tropical Atlantic: A satellite perspective, J. Geophys. Res., 108, 4237,, 2003. a

EPIC Science Team: DSCOVR/EPIC L2_SO2O3AI HDF-5 File – Version 1,, 2018. a

Eskes, H. J., van der A, R. J., Brinksma, E. J., Veefkind, J. P., de Haan, J. F., and Valks, P. J. M.: Retrieval and validation of ozone columns derived from measurements of SCIAMACHY on Envisat, Atmos. Chem. Phys. Discuss., 5, 4429–4475,, 2005. a, b

Fortuin, J. P. F. and Kelder, H.: Possible links between ozone and temperature profiles, Geophys. Res. Lett., 23, 1517–1520,, 1996. a

Fortuin, J. P. F. and Kelder, H.: An ozone climatology based on ozonesonde and satellite measurements, J. Geophys. Res.-Atmos., 103, 31709–31734,, 1998. a, b, c

Froidevaux, L., Allen, M., Berman, S., and Daughton, A.: The mean ozone profile and its temperature sensitivity in the upper stratosphere and lower mesosphere: An analysis of LIMS observations, J. Geophys. Res., 94, 6389,, 1989. a

Gelaro, R., McCarty, W., Suárez, M. J., Todling, R., Molod, A., Takacs, L., Randles, C. A., Darmenov, A., Bosilovich, M. G., Reichle, R., Wargan, K., Coy, L., Cullather, R., Draper, C., Akella, S., Buchard, V., Conaty, A., da Silva, A. M., Gu, W., Kim, G. K., Koster, R., Lucchesi, R., Merkova, D., Nielsen, J. E., Partyka, G., Pawson, S., Putman, W., Rienecker, M., Schubert, S. D., Sienkiewicz, M., and Zhao, B.: The modern-era retrospective analysis for research and applications, version 2 (MERRA-2), J. Climate, 30, 5419–5454,, 2017. a

Haefele, A., Hocke, K., Kämpfer, N., Keckhut, P., Marchand, M., Bekki, S., Morel, B., Egorova, T., and Rozanov, E.: Diurnal changes in middle atmospheric H2O and O3: Observations in the Alpine region and climate models, J. Geophys. Res., 113, D17303,, 2008. a, b, c

Herman, J., Huang, L., McPeters, R., Ziemke, J., Cede, A., and Blank, K.: Synoptic ozone, cloud reflectivity, and erythemal irradiance from sunrise to sunset for the whole earth as viewed by the DSCOVR spacecraft from the earth–sun Lagrange 1 orbit, Atmos. Meas. Tech., 11, 177–194,, 2018. a

Hoinka, K. P.: Statistics of the Global Tropopause Pressure, Mon. Weather Rev., 126, 3303–3325,<3303:SOTGTP>2.0.CO;2, 1998. a

Hoogen, R., Rozanov, V. V., and Burrows, J. P.: Ozone profiles from GOME satellite data: Algorithm description and first validation, J. Geophys. Res.-Atmos., 104, 8263–8280,, 1999. a, b, c

Huang, F. T., Mayr, H. G., Russell, J. M., Mlynczak, M. G., and Reber, C. A.: Ozone diurnal variations and mean profiles in the mesosphere, lower thermosphere, and stratosphere, based on measurements from SABER on TIMED, J. Geophys. Res.-Space Phys., 113, A04307,, 2008. a, b, c

Jiang, Z., Miyazaki, K., Worden, J. R., Liu, J. J., Jones, D. B. A., and Henze, D. K.: Impacts of anthropogenic and natural sources on free tropospheric ozone over the Middle East, Atmos. Chem. Phys., 16, 6537–6546,, 2016. a

Klenk, K. F., Bhartia, P. K., Kaveeshwar, V. G., McPeters, R. D., Smith, P. M., and Fleig, A. J.: Total Ozone Determination from the Backscattered Ultraviolet (BUV) Experiment, J. Appl. Meteorol., 21, 1672–1684,<1672:TODFTB>2.0.CO;2, 1982. a

Krzyścin, J. W., Degórska, M., and Rajewska-Wȩch, B.: Seasonal acceleration of the rate of total ozone decreases over Central Europe: impact of tropopause height changes, J. Atmos. Sol.-Terr. Phys., 60, 1755–1762,, 1998. a

Labow, G. J., Ziemke, J. R., McPeters, R. D., Haffner, D. P., and Bhartia, P. K.: A total ozone-dependent ozone profile climatology based on ozonesondes and Aura MLS data, J. Geophys. Res.-Atmos., 120, 2537–2545,, 2015. a, b, c, d, e

Lamsal, L. N., Weber, M., Tellmann, S., and Burrows, J. P.: Ozone column classified climatology of ozone and temperature profiles based on ozonesonde and satellite data, J. Geophys. Res.-Atmos., 109, 1–15,, 2004. a, b, c, d, e

Lelieveld, J. and Dentener, F. J.: What controls tropospheric ozone?, J. Geophys. Res.-Atmos., 105, 3531–3551,, 2000. a

Lerot, C., Van Roozendael, M., Lambert, J.-C., Granville, J., van Gent, J., Loyola, D., and Spurr, R. J. D.: The GODFIT algorithm: a direct fitting approach to improve the accuracy of total ozone measurements from GOME, Int. J. Remote Sens., 31, 543–550,, 2010. a, b

Lerot, C., Van Roozendael, M., Spurr, R. J. D., Loyola, D., Coldewey-Egbers, M. R., Kochenova, S., Van Gent, J., Koukouli, M., Balis, D., Lambert, J. C., Granville, J., and Zehner, C.: Homogenized total ozone data records from the European sensors GOME/ERS-2, SCIAMACHY/Envisat, and GOME-2/MetOp-A, J. Geophys. Res.-Atmos., 119, 1639–1662,, 2014. a, b

Li, Q., Jacob, D. J., Logan, J. A., Bey, I., Yantosca, R. M., Liu, H., Martin, R. V., Fiore, A. M., Field, B. D., Duncan, B. N., and Thouret, V.: A tropospheric ozone maximum over the Middle East, Geophys. Res. Lett., 28, 3235–3238,, 2001. a

Liu, G., Liu, J., Tarasick, D. W., Fioletov, V. E., Jin, J. J., Moeini, O., Liu, X., Sioris, C. E., and Osman, M.: A global tropospheric ozone climatology from trajectory-mapped ozone soundings, Atmos. Chem. Phys., 13, 10659–10675,, 2013. a, b, c, d, e

Liu, J. J., Jones, D. B. A., Zhang, S., and Kar, J.: Influence of interannual variations in transport on summertime abundances of ozone over the Middle East, J. Geophys. Res., 116, D20310,, 2011. a

Liu, X., Chance, K. V., Sioris, C. E., Spurr, R. J. D., Kurosu, T. P., Martin, R. V., and Newchurch, M. J.: Ozone profile and tropospheric ozone retrievals from the Global Ozone Monitoring Experiment: Algorithm description and validation, J. Geophys. Res., 110, D20307,, 2005. a, b, c

Liu, X., Bhartia, P. K., Chance, K., Spurr, R. J. D., and Kurosu, T. P.: Ozone profile retrievals from the Ozone Monitoring Instrument, Atmos. Chem. Phys., 10, 2521–2537,, 2010. a

Loyola, D. G., Koukouli, M. E., Valks, P., Balis, D. S., Hao, N., Van Roozendael, M., Spurr, R. J. D., Zimmer, W., Kiemle, S., Lerot, C., and Lambert, J.-C.: The GOME-2 total column ozone product: Retrieval algorithm and ground-based validation, J. Geophys. Res., 116, D07302,, 2011. a, b

Maeda, K. and Heath, D. F.: Asymmetries of the Upper Stratospheric Ozone Distribution Between Two Hemispheres, J. Atmos. Sci., 40, 1353–1359,<1353:AOTUSO>2.0.CO;2, 1983. a

Martin, R. V., Jacob, D. J., Logan, J. A., Bey, I., Yantosca, R. M., Staudt, A. C., Li, Q., Fiore, A. M., Duncan, B. N., Liu, H., Ginoux, P., and Thouret, V.: Interpretation of TOMS observations of tropical tropospheric ozone with a global model and in situ observations, J. Geophys. Res., 107, 4351,, 2002. a

Mateer, C. L., Heath, D. F., and Krueger, A. J.: Estimation of Total Ozone from Satellite Measurements of Backscattered Ultraviolet Earth Radiance, J. Atmos. Sci., 28, 1307–1311,<1307:EOTOFS>2.0.CO;2, 1971. a

McPeters, R. D. and Labow, G. J.: Climatology 2011: An MLS and sonde derived ozone climatology for satellite retrieval algorithms, J. Geophys. Res. Atmos., 117, 1–8,, 2012. a, b, c, d, e, f, g, h

McPeters, R. D., Labow, G. J., and Logan, J. A.: Ozone climatological profiles for satellite retrieval algorithms, J. Geophys. Res. Atmos., 112, 1–9,, 2007. a, b, c, d, e, f, g, h, i

Miles, G. M., Siddans, R., Kerridge, B. J., Latter, B. G., and Richards, N. A. D.: Tropospheric ozone and ozone profiles retrieved from GOME-2 and their validation, Atmos. Meas. Tech., 8, 385–398,, 2015. a, b, c

Moxim, W. J. and Levy, H.: A model analysis of the tropical South Atlantic Ocean tropospheric ozone maximum: The interaction of transport and chemistry, J. Geophys. Res.-Atmos., 105, 17393–17415,, 2000. a

Ohring, G. and Muench, H. S.: Relationships between ozone and meteorological parameters in the lower stratosphere, J. Meteorol., 17, 195–206,<0195:RBOAMP>2.0.CO;2, 1960. a

Perliski, L. M. and London, J.: Satellite observed long-term averaged seasonal and spatial ozone variations in the stratosphere, Planet. Space Sci., 37, 1509–1525,, 1989. a

Reed, R. J.: The role of vertical motions in ozone-weather relationships, J. Meteorol., 7, 263–267,<0263:TROVMI>2.0.CO;2, 1950. a

Rodgers, C. D.: Inverse Methods for Atmospheric Sounding, World Scientific, Singapore,, 2000. a, b

Rood, R. B. and Douglass, A. R.: Interpretation of ozone temperature correlations: 1. Theory, J. Geophys. Res., 90, 5733,, 1985. a, b

Sakazaki, T., Fujiwara, M., Mitsuda, C., Imai, K., Manago, N., Naito, Y., Nakamura, T., Akiyoshi, H., Kinnison, D., Sano, T., Suzuki, M., and Shiotani, M.: Diurnal ozone variations in the stratosphere revealed in observations from the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES) on board the International Space Station (ISS), J. Geophys. Res.-Atmos., 118, 2991–3006,, 2013. a, b, c

Salby, M. L. and Callaghan, P. F.: Fluctuations of total ozone and their relationship to stratospheric air motions, J. Geophys. Res.-Atmos., 98, 2715–2727,, 1993. a

Sauvage, B., Martin, R. V., van Donkelaar, A., and Ziemke, J. R.: Quantification of the factors controlling tropical tropospheric ozone and the South Atlantic maximum, J. Geophys. Res., 112, D11309,, 2007. a

Schranz, F., Fernandez, S., Kämpfer, N., and Palm, M.: Diurnal variation in middle-atmospheric ozone observed by ground-based microwave radiometry at Ny-Ålesund over 1 year, Atmos. Chem. Phys., 18, 4113–4130,, 2018. a, b, c

Seidel, D. J. and Randel, W. J.: Variability and trends in the global tropopause estimated from radiosonde data, J. Geophys. Res., 111, D21101,, 2006. a

Smith, A. K.: Numerical simulation of global variations of temperature, ozone, and trace species in the stratosphere, J. Geophys. Res.-Atmos., 100, 1253–1269,, 1995. a

Sofieva, V. F., Tamminen, J., Kyrölä, E., Mielonen, T., Veefkind, P., Hassler, B., and Bodeker, G. E.: A novel tropopause-related climatology of ozone profiles, Atmos. Chem. Phys., 14, 283–299,, 2014. a, b, c, d, e, f, g, h, i, j, k, l

Steinbrecht, W., Claude, H., Köhler, U., and Hoinka, K. P.: Correlations between tropopause height and total ozone: Implications for long-term changes, J. Geophys. Res.-Atmos., 103, 19183–19192,, 1998. a, b

Thompson, A. M., Doddridge, B. G., Witte, J. C., Hudson, R. D., Luke, W. T., Johnson, J. E., Johnson, B. J., Oltmans, S. J., and Weller, R.: A tropical Atlantic Paradox: Shipboard and satellite views of a tropospheric ozone maximum and wave-one in January–February 1999, Geophys. Res. Lett., 27, 3317–3320,, 2000. a

Van Roozendael, M., Loyola, D., Spurr, R., Balis, D., Lambert, J.-C., Livschitz, Y., Valks, P., Ruppert, T., Kenter, P., Fayt, C., and Zehner, C.: Ten years of GOME/ERS-2 total ozone data – The new GOME data processor (GDP) version 4: 1. Algorithm description, J. Geophys. Res., 111, D14311,, 2006. a, b

Van Roozendael, M., Spurr, R. J. D., Loyola, D., Lerot, C., Balis, D., Lambert, J. C., Zimmer, W., Van Gent, J., Van Geffen, J., Koukouli, M., Granville, J., Doicu, A., Fayt, C., and Zehner, C.: Sixteen years of GOME/ERS-2 total ozone data: The new direct-fitting GOME Data Processor (GDP) version 5Algorithm description, J. Geophys. Res. Atmos., 117, 1–18,, 2012. a, b

Varotsos, C., Cartalis, C., Vlamakis, A., Tzanis, C., and Keramitsoglou, I.: The Long-Term Coupling between Column Ozone and Tropopause Properties, J. Climate, 17, 3843–3854,<3843:TLCBCO>2.0.CO;2, 2004. a

Veefkind, J., de Haan, J., Brinksma, E., Kroon, M., and Levelt, P.: Total ozone from the ozone monitoring instrument (OMI) using the DOAS technique, IEEE T. Geoscience Remote, 44, 1239–1244,, 2006. a, b

Wang, P.-H., Fishman, J., Harvey, V. L., and Hitchman, M. H.: Southern tropical upper tropospheric zonal ozone wave-1 from SAGE II observations (1985–2002), J. Geophys. Res., 111, D08305,, 2006. a

Wargan, K., Pawson, S., Olsen, M. A., Witte, J. C., Douglass, A. R., Ziemke, J. R., Strahan, S. E., and Nielsen, J. E.: The global structure of upper troposphere-lower stratosphere ozone in GEOS-5: A multiyear assimilation of EOS Aura data, J. Geophys. Res.-Atmos., 120, 2013–2036,, 2015. a, b, c, d

Wargan, K., Labow, G., Frith, S., Pawson, S., Livesey, N., and Partyka, G.: Evaluation of the Ozone Fields in NASA's MERRA-2 Reanalysis, J. Climate, 30, 2961–2988,, 2017. a, b, c, d

Wassmann, A., Borsdorff, T., aan de Brugh, J. M. J., Hasekamp, O. P., Aben, I., and Landgraf, J.: The direct fitting approach for total ozone column retrievals: a sensitivity study on GOME-2/MetOp-A measurements, Atmos. Meas. Tech., 8, 4429–4451,, 2015. a, b

Wei, J. C., Pan, L. L., Maddy, E., Pittman, J. V., Divarkarla, M., Xiong, X., and Barnet, C.: Ozone profile retrieval from an advanced infrared sounder: experiments with tropopause-based climatology and optimal estimation approach, J. Atmos. Ocean. Technol., 27, 1123–1139,, 2010. a, b, c, d, e

Weiss, A. K., Staehelin, J., Appenzeller, C., and Harris, N. R. P.: Chemical and dynamical contributions to ozone profile trends of the Payerne (Switzerland) balloon soundings, J. Geophys. Res.-Atmos., 106, 22685–22694,, 2001. a

Wellemeyer, C. G., Taylor, S. L., Seftor, C. J., Mcpeters, R. D., and Bhartia, P. K.: A correction for total ozone mapping spectrometer profile shape errors at high latitude, J. Geophys. Res., 102, 9029,, 1997. a, b, c

Yang, K.: OMPS-NPP L2 NM Sulfur Dioxide (SO2) Total and Tropospheric Column swath orbital V2, Greenbelt, MD, USA, Goddard Earth Sciences Data and Information Services Center (GES DISC),, 2017.  a

Ziemke, J. R., Chandra, S., Labow, G. J., Bhartia, P. K., Froidevaux, L., and Witte, J. C.: A global climatology of tropospheric and stratospheric ozone derived from Aura OMI and MLS measurements, Atmos. Chem. Phys., 11, 9237–9251,, 2011. a, b, c, d

Ziemke, J. R., Olsen, M. A., Witte, J. C., Douglass, A. R., Strahan, S. E., Wargan, K., Liu, X., Schoeberl, M. R., Yang, K., Kaplan, T. B., Pawson, S., Duncan, B. N., Newman, P. A., Bhartia, P. K., and Heney, M. K.: Assessment and applications of NASA ozone data products derived from Aura OMI/MLS satellite measurements in context of the GMI chemical transport model, J. Geophys. Res.-Atmos., 119, 5671–5699,, 2014. a

Short summary
We constructed total-ozone-dependent and tropopause-dependent climatologies from MERRA-2 ozone data to describe the dynamic variations in the ozone profile in response to changing meteorological conditions. The new climatologies contain the first quantitative characterization of ozone profile covariances, which facilitate a new approach to improve ozone profiles using the most probable patterns of profile adjustments represented by the empirical orthogonal functions of the covariance matrices.