the Creative Commons Attribution 4.0 License.

the Creative Commons Attribution 4.0 License.

# Satellite data validation: a parametrization of the natural variability of atmospheric mixing ratios

### Alexandra Laeng

### Thomas von Clarmann

### Quentin Errera

### Udo Grabowski

### Shawn Honomichl

High-resolution model data are used to estimate the statistically typical mixing ratio variabilities of trace species as a function of distance and time separation. These estimates can be used to explain the fact that some of the differences between observations made with different observing systems are due to the less-than-perfect co-location of the measurements. The variability function is approximated by a two-parameter regression function, and lookup tables of the natural variability values as a function of distance separation and time separation are provided. In addition, a reparametrization of the variability values as a function of latitudinal gradients is proposed, and the seasonal independence of the linear approximation of such a function is demonstrated.

This paper tackles a problem that typically arises when remotely sensed data
from different instruments are compared within the framework of validation
studies. In quantitative validation, the common approach is to calculate
differences between pairs of measurements of the same air mass from the two instruments
under comparison. With the aid of *χ*^{2} statistics, it is tested if
the observed differences can be explained by the estimated error of the
differences (Rodgers and Connor, 2003). The estimated error of the differences includes measurement noise as well as parameter and model errors,
as far as those coming from different instruments are uncorrelated (von Clarmann, 2006). Furthermore, the different impact of prior
information on the result has to be considered. However,
the differences will often be too large to be explained by the combined error budget
of the measurements under comparison. The reason for this is that the instruments
typically do not sound exactly the same air mass. Spatial and temporal mismatch of the
measurements as well as natural variability of the measured state variable
contribute to the observed differences. This source of differences
is only quantified in validation papers in a few exceptional cases
(see, e.g., Sheese et al., 2021, for an example where models were used to quantify the related effect of ozone variability).
Instead, natural variability is often used as a universal excuse to defend measurements if validation studies suggest that the related retrieval errors are underestimated.

In this study, we present a user-friendly tool to provide quantitative estimates of the component of the observation differences that can be attributed to the spatial and temporal mismatch and natural variability. The underlying method is based on highly resolved model fields of temperature and mixing ratios of trace species, as described in Sect. 2. These model fields were smoothed according to the horizontal resolution of the instruments whose precision was to be validated; in this study, we have chosen the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS; Fischer et al., 2008) as an example. From these smoothed fields, the typical variabilities are evaluated as a function of spatial and temporal mismatch (Sect. 3). In order to avoid unnecessarily high data traffic and to reduce the impact of model imperfections on the calculated fields as much as possible, parametrizations of these dependencies (for different trace gases, altitudes and latitude bins) are developed by prescribing a particular shape to the natural variability function, which arises from the general theory of random functions with stationary increments (Sect. 4) and is confirmed by calculations using model data. A reparametrization is developed to improve the validity of the inductive generalization towards other gases and seasons (Sect. 5). Instructions on how to use this reparametrization are given in Sect. 6. The adequacy of our suggested method is critically discussed in Sect. 7, and final recommendations are given in Sect. 8.

The model fields used in this study came from two models.

## 2.1 BASCOE (Belgian Assimilation System for Chemical ObsErvations)

The main set of model fields has been produced by the Belgian Assimilation System for Chemical ObsErvations (BASCOE; Errera et al., 2008).
While generally used in the context of stratospheric chemical data assimilation (e.g., Errera et al., 2019), the chemistry transport model (CTM)
of the system has also been used to study the evolution of the stratospheric composition (Chabrillat et al., 2018; Prignon et al., 2019; Minganti et al., 2020).
Here, BASCOE was run for the period from 25 September to 1 October 2008, and wind and temperature were taken from the ERA-Interim reanalysis (Dee et al., 2011).
The model was run on a $\mathrm{1}{}^{\circ}\times \mathrm{1}{}^{\circ}$ horizontal grid using the native 60 vertical levels of ERA-Interim (from the surface to 0.1 hPa) and a time step of 30 min. Hourly global fields of 28 relevant trace gases^{1} were used for this study.

## 2.2 WACCM6 (Whole Atmosphere Community Climate Model 6)

The auxiliary set of model fields used came from the Whole Atmosphere Community Climate Model 6 (WACCM6), which is the atmospheric component
of the Community Earth System Model,
version 2 (CESM2; Danabasoglu et al., 2019; Emmons et al., 2019; Gettelman et al., 2019; Tilmes et al., 2019), run at the Atmospheric Chemistry Observations and Modelling Laboratory of the National Center of Atmospheric Research (UCAR/NCAR/ACOM).
The model has a horizontal resolution of $\mathrm{0.9}{}^{\circ}\times \mathrm{1.25}{}^{\circ}$ with 88 vertical hybrid sigma-pressure levels and is run using specified dynamics, with nudging of temperature, and $\mathit{U}/\mathit{V}$ winds from the NASA Goddard Earth Observing System, version 5 (GEOS-5), forecast model.
The model simulations contain the fields of three species (O_{3}, H_{2}O and NO) for 4 weeks in the year 2020, one in each season. The data were regridded on the same fixed height grid with a 1 km step using the fields of geopotential
height and temperature provided with the data.

## 3.1 Structure functions

Let *X* be a random variable defined by the amount of the target trace gas in a given infinitely
small air parcel, centered around a point in the atmosphere at a given moment in time and
reported as a volume mixing ratio (VMR).
The amount of the trace gas at any point in the atmosphere at a fixed time
(or at any moment in time at a fixed point in the atmosphere) can be viewed as a state of
a one-dimensional random process *X*(*t*), where *t* parameterizes the distance
from an initial point (or the time elapsed from an initial moment). This random process is nonstationary,
as its statistical characteristics can change with *t*. The increments $X(t+\mathit{\tau})-X\left(t\right)$ of the process *X*(*t*)
represent the change in the amount of the trace gas with distance (or over time).
Based on the literature, we assume that, in a given sufficiently narrow latitude band, at a fixed altitude
and in a given season, the distribution of the differences $X(t+\mathit{\tau})-X\left(t\right)$ does not depend on *t*,
which means that *X*(*t*) is a process with stationary increments.
The basic characteristics of real-valued random process with stationary increments are the mean value of the increment $E\left[X(t+\mathit{\tau})-X\left(t\right)\right]$ and the correlation function of the increment

also called the “structure function” of the process *X*(*t*) (Yaglom, 1986, chap. 23). In our case, as $E\left[X(t+\mathit{\tau})-X\left(t\right)\right]=\mathrm{0}$,

The natural variability of a trace gas is the square root of the structure function of the process *X*(*t*):

This provides the formal link between the intuitive definition of the natural variability, as the variability of differences, and the mathematical machinery of random processes with stationary increments, which is widely used in studying the processes at smaller spatiotemporal scales, for example, in the theory of atmospheric turbulence. This link will allow us to draw conclusions about the nature of the process *X*(*t*) based on the shape of the obtained statistics and will justify the choice of the form of the regression function for the natural variability of trace gases. The next section explains how the estimation of $\sqrt{D}$ is done out of model fields.

## 3.2 Estimation of variabilities

In order to obtain a statistic of the variability of differences, the model fields were first transformed from their native hybrid sigma-pressure vertical grid to a fixed 1 km step geometrical height grid. For this, the geopotential height for each model knot was restored and transformed into geometrical height using the temperature values from the model, which allowed the interpolation of the profiles on a fixed altitude grid. Second, the model fields were smoothed according to the horizontal resolution of the instruments whose precision was to be validated. We have chosen the MIPAS instrument as an example. Its cross-track resolution roughly corresponds to an east–west resolution and is defined by the width of the instantaneous field of view at the tangent point, which is 30 km. The along-track smearing, which roughly corresponds to a north–south horizontal resolution, is roughly 200 km on average (von Clarmann et al., 2009). This smoothing operation does not influence the shape of the obtained curves, but it does reduce the obtained variability values by around 0.05 %. No vertical smoothing is applied, as vertical smoothing is typically considered in an explicit manner in validation via the averaging kernels. These smoothed fields are the basis for the statistics of the horizontal and temporal variability of the atmospheric state.

### 3.2.1 Horizontal variability

We take the model data within a fixed 10^{∘} latitude bin and at a fixed height; as each of the five model data sets used (one from BASCOE and four from WACCM6) covers only 1 week, there is no need to fix a season.
For all possible pairs of points in the obtained subset, the normalized differences of the VMR of the target trace species within a predefined radius
of 1500 km are calculated:

with VMR_{mean} being the mean VMR values of the target trace gas in the chosen latitude band at the chosen height. The constant time index *t* indicates that only differences
where the subtractor and the subtrahend refer to the same time are considered.
These differences are binned according to their horizontal separation distance.
The following bins were used: 0–100 km between
the two points, 100–200 km between
the two points, etc. We calculate the standard deviation of the sample of these normalized differences, which provides us with an estimator of the natural variability of the target trace gas as a function of distance separation.
The obtained estimation of the natural variability of ozone at 35 km altitude as a function of distance is shown in Fig. 1a. The fast growth of the variability for separation distances over 1000 km at high northern latitudes reflects that, in many pairs of corresponding samples, one point lies inside the polar vortex whereas the other lies outside. Note also that the calculated variability values for the distances from 0 to 100 km are zero for the tropical latitudes (yellow curves in Fig. 1) and present a peak at subtropical latitudes (clear orange and clear green curves in Fig. 1). This is due to the models' resolutions: the samples for 100 km distance separation
are empty or very small at low latitudes. Therefore, the 100 km point will not be taken into account in the calculation of regression coefficients at these latitudes.

### 3.2.2 Temporal variability

In a similar fashion to that described above, for all possible pairs of data points for the entire data set, the differences of the VMR of the target trace gas within a predefined time period of 72 h are calculated for each latitude band and height:

The constant location index “location” indicates that
differences are only considered when the subtractor and the subtrahend refer to the same location. These differences are sorted according to their time lag.
Similar to the horizontal variability, for each time lag, the differences are normalized by the mean VMR within given latitude band at given altitude, and the standard deviation of the sample of normalized differences is then calculated. This quantity is the estimator of the natural variability of the target species
as a function of time separation, and we note it as *σ*_{rel,time}; its values for ozone at 35 km altitude are shown in Fig. 1b. As for most satellite validation exercises, time separation within co-location criteria stays within 5 h, and we made the choice to restrain our analysis to a maximum time separation lag of 5 h.

## 3.3 Combination of horizontal and temporal variability

Despite the fact that advection can admittedly cause correlations between horizontal and temporal components of the variability, at the scales considered here, we assume that horizontal and temporal variations of the atmospheric state are uncorrelated. Tests using a statistic of combined horizontal and temporal differences of the type

have shown that, at the scales considered here, the error due to the neglect of correlations is below 0.1 % and, thus, not usually worth an additional effort. In our analysis, we offer independent parametrizations for each of them, which are recommended to be combined with their quadratic sum; we also provide a software that performs this summation for a reparametrization of these quantities on latitudinal gradients, and this software is described in the Sect. 6.

We do not consider the vertical variability in the present work, and the fields are calculated independently for each altitude level. Note that atmospheric variability as a function of distance and time separation could be approximately represented by a two-dimensional random process; however, this is outside the scope of this paper: our choice is to treat the distance and time mismatch dependences separately, as this is what a typical validation exercise does.

## 4.1 Motivation

The goal of the present work is to provide the community with information on the natural variability of trace gas mixing ratios as a function of distance and time separation. This information is meant to be used in the context of validation studies. Instead of providing the entire variability data set, we consider the use of a simple and easy-to-use parametrization as more adequate. The reasons for this are as follows: (1) the use of parametrizations avoids a considerable amount of data traffic, (2) the fine structure of the fields reflects the actual conditions of the days actually covered by the model run rather than the general behavior of the atmosphere, and (3) our parametrization using continuous regression functions allows for easy interpolation.

## 4.2 The regression function

In view of the shape of the curves produced from model data (Sect. 3.2.2 and 3.2.1), the natural variability function can be parameterized in the form *D*(*τ*)=*A**τ*^{γ}, where *A*>0, and $\mathrm{0}<\mathit{\gamma}<\mathrm{1}$. An interesting side conclusion that can be made from the obtained shape of the structure function of *X*(*t*) is that the process of the atmospheric variability (horizontal and temporal variability) of mixing ratios is self-similar: the form of its structure function *D*(*τ*)=*A**τ*^{γ} is invariant under a group of similarity transformations,
*t*→*h**t* and *X*→*a*(*h*)*X* (Kolmogorov, 1940; Yaglom, 1986); in other words, no characteristic scale can be associated with their structure function.

As pointed out in the Sect. 3.2.1, variability grows very rapidly at high latitudes and distances
over 1000 km. Moreover, at low latitudes, the values of
variability as a function of the distance mismatch is meaningless at 100 km: it is calculated on the samples from too small to empty, due to the model horizontal resolution. The choice of a 5 h upper limit for the time-dependent variability is driven by the typical values of the time mismatch occurring in satellite validation studies and on the shape of the obtained curves. We calculate the regression coefficients *A* and *γ* by minimizing the quantity

via the Sequential Least SQuares Programming (SLSQP) optimizer (Kraft, 1988) and by giving the value at 100 km as a first constraint for *A* and 0.5 as a first constraint for *γ*. The obtained regression curves and the initially calculated model curves for some species are shown in
Fig. 2 for distance separation and in
Fig. 3 for time separation. Note that the assumed parameterization does not apply for a high-latitude case at large distances. The initial and smoothed (regressed) natural variability
surfaces as a function of latitude and time separation are shown in Fig. 4 for ozone at 35 km altitude. In a range between 100 and 1000 km, the parametrizations fit the data very well.

An obvious deficiency of our approach is that the variability fields are calculated from only 1 week of data. Operational constraints did not allow us to generate a better coverage for that many species at the required resolution. The variability, which we have calculated as a function of latitude and distance (time) mismatch, is also dependent on season, as different seasons correspond to different inclinations of the Earth's axis and, intuitively, should all be shifted in latitude as the season changes. There is, however, a way to tackle the problem: if the user is willing to calculate one additional quantity from this data, namely the latitudinal gradients of the gas under validation, then the workaround consists of the reparametrization of the variability of the latitudinal gradients of the species. Latitudinal gradients of a gas are defined as follows:

where ${\text{VMR}}_{{l}_{\mathrm{1}}}$ is the mean VMR of the gas in a latitude band, ${\text{VMR}}_{{l}_{\mathrm{2}}}$ is the mean VMR in the northward neighboring latitude band, and *l*_{1}−*l*_{2} is the width of the latitude band (in our case, this is 10^{∘}). In the relative version of the latitudinal gradient, this quantity is normalized with respect to the mean VMR in both bands and is multiplied by 100. For a fixed distance (or time) mismatch, as a first approximation,
because the variability is calculated as the square root of the variance of differences, and the latitudinal gradients are calculated as differences,
we expect a linear dependence of the variability from the latitudinal gradients.

To test the abovementioned theoretical considerations in practice, especially to see if the linear dependence of the natural variability from the latitudinal gradients is the same in different seasons, the link between the natural variability and the latitudinal gradients was tested on the data of a worse resolved model with less species, although using a model for which the data from different seasons were available, namely data from WACCM6.

Figure 5 shows the ozone distributions in each season of WACCM data. The distributions are as expected, and the latitudinal shift in different season is visible; therefore, the data are suitable for our test.

Figures 6, 7 and 8 show the natural variability as a function of latitudinal gradients for a particular distance separation of 400 km in four seasons of WACCM for H_{2}O at 30 km, O_{3} at 35 km and NO at 40 km. Each point in these figures corresponds to one 10^{∘} latitude band; for the sake of completeness, we include the points from all of the latitudes. One can observe that the regression lines are similar in all cases, with some deviation in summer and winter
that stems from the points corresponding to high (>70^{∘}) latitudes, which is expected. Similar linear dependencies, which are comparable for different seasons, are observed when changing the 400 km separation distance. This hints at a seasonally independent linear approximation of the natural variability as a function of latitudinal gradients. Hence, the natural variability calculated in just one season and reparameterized on latitudinal gradients of the species provides all of the information needed for the validation exercise in any season. Finally, we would like to remark that the variability fields of well-correlated tracers are close, which provides additional confidence in the method.

The variabilities as functions of time and distance mismatch are added quadratically, and this provides the final variability value for the
chosen co-location criteria. In practice, the users will have to calculate just one additional quantity from the data under validation,
namely the latitudinal gradients of the species under validation.
This quantity should be calculated on the whole sample in order to increase the significance of the statistics.
In addition to the regression coefficient values, we provide a software that uses the species' name
(among the 30 available), the distance and time mismatch chosen, the latitude band, the value of the latitudinal gradients of these data and the height as input and then outputs the value of the natural variability of the gas at a given altitude, latitude band and mismatch.
If the validation study is performed in a latitude domain larger then 10^{∘}, the values in the corresponding 10^{∘} bands should be added quadratically.

In the sense of playing devil's advocate, we try to raise possible objections to our method and rebut them. As we do not use model data directly, instead using only differences between model data, the additive model biases cancel out. Critical minds might plead that there could still be multiplicative biases in the model data that would affect our statistics of differences. However, these are not harmful if the gradient-related parametrization is used instead of the latitude–month-related parametrization. The rough reason for this is as follows: a model bias affects the horizontal gradients in the same way as the differences used for our statistics. Thus, the effect of a multiplicative model bias also cancels out.

An obvious objection to our method is that the model data used cover only a short time period and might not be inductively generalizable towards other time periods. We agree that, due to the annual cycle, the typical meteorological regimes are shifted in latitude over the year. However, once again, when the gradient-related parametrization is used, it is likely that the statistics of the correct meteorological regime is chosen, even if it is found at different latitudes to that of the validation experiment. The explanation for this is that the natural variability of mixing ratios of most trace species is predominantly driven by the latitudinal gradients. It goes without saying that this does not hold for fast-reacting species, particularly those that are in photochemical equilibrium.

These parametrizations should not be used where polar vortex dynamics may play a role or for spatial mismatches beyond 1000 km and temporal mismatches beyond 5 h; however, these situations are not the preferred validation scenarios anyway.

In validation exercises, a universal excuse used to explain the residual discrepancy between the data is the natural atmospheric variability due to imperfect co-locations. This work is the first attempt to quantify this atmospheric variability for a large sample of atmospheric constituents and to provide the user with a tool to subtract the natural atmospheric variability portion from the residual variability. The fields of natural atmospheric variability as a function of distance and time mismatch were calculated from highly resolved BASCOE model data. The variability data were described with an easy-to-use regression function, and the regression coefficients as well as the software that calculates the value of the corresponding natural variability for given gas, latitudinal gradient, height, and co-location criteria are provided to the community. The independence of the linear approximation of the natural variability as a function of latitudinal gradients from season was demonstrated on WACCM model data. An application of the method to ozone and temperature fields will be provided in upcoming validation papers of version 8 of the MIPAS data.

The regression coefficients of the parametrization on latitudes for autumn and the software for calculating the variability as a function of latitudinal gradients are available from https://doi.org/10.5445/IR/1000137514 (Laeng, 2021).

Whole Atmosphere Community Climate Model (WACCM) data can be downloaded from the ACOM website: https://www2.acom.ucar.edu/gcm/waccm (Atmospheric Chemistry Observations & Modeling/National Center for Atmospheric Research/University Corporation for Atmospheric Research, 2020c). WACCM forecast maps are available from https://www.acom.ucar.edu/waccm/forecast/ (Atmospheric Chemistry Observations & Modeling/National Center for Atmospheric Research/University Corporation for Atmospheric Research, 2020a).

The WACCM model can be downloaded from the NCAR/UCAR Research Data Archive: https://rda.ucar.edu/datasets/ds313.6/ (Atmospheric Chemistry Observations & Modeling/National Center for Atmospheric Research/University Corporation for Atmospheric Research, 2020b).

AL carried out the statistical analysis and developed the parametrization and reparametrization. TvC identified the problem, suggested the solution, supervised the work and contributed to the text. UG ensured the implementation of statistical calculations out of model fields. QE and SH provided the model fields and guidance to their use for the BASCOE and WACCM models, respectively.

At least one of the (co-)authors is a member of the editorial board of *Atmospheric Measurement Techniques*. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.

Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

This article is part of the special issue “Towards Unified Error Reporting (TUNER)”. It is not associated with a conference.

This work was performed as part of the European Space Agency (ESA) VACUUM'R project. The authors are thankful to Claus Zehner for very effective supervision of the project. The WACCM data were provided by the WACCM forecast team, including Simone Tilmes and Douglas Kinnisoni, and the website/technical team, including Carl Drews and Garth DAttilo. ISSI financed and hosted the second SPARC TUNER (Toward UNified Error Reporting) meeting, where the intermediate results of this work were presented and critically accessed by the members of the TUNER consortium. The first author is particularly thankful to Viktoria Sofieva for numerous helpful discussions.

This research has been supported by the ESA Validation And Calibration of Uncertainty estimates Using Model data Revisited
(VACUUM'R) project.

The article processing charges for this open-access publication were covered by the Karlsruhe Institute of Technology (KIT).

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

Atmospheric Chemistry Observations & Modeling/National Center for Atmospheric Research/University Corporation for Atmospheric Research: WACCM forecast maps, updated daily, NCAR/UCAR [data set], https://www.acom.ucar.edu/waccm/forecast (last access: 11 April 2022), 2020a. a

Atmospheric Chemistry Observations & Modeling/National Center for Atmospheric Research/University Corporation for Atmospheric Research: Whole Atmosphere Community Climate Model (WACCM) Model Output, updated daily, Research Data Archive at the National Center for Atmospheric Research, Computational and Information Systems Laboratory [data set], https://rda.ucar.edu/datasets/ds313.6/ (last access: 11 April 2022), 2020b. a

Atmospheric Chemistry Observations & Modeling/National Center for Atmospheric Research/University Corporation for Atmospheric Research: Whole Atmosphere Community Climate Model (WACCM), NCAR/UCAR [code], https://www2.acom.ucar.edu/gcm/waccm (last access: 11 April 2022), 2020c. a

Chabrillat, S., Vigouroux, C., Christophe, Y., Engel, A., Errera, Q., Minganti, D., Monge-Sanz, B. M., Segers, A., and Mahieu, E.: Comparison of mean age of air in five reanalyses using the BASCOE transport model, Atmos. Chem. Phys., 18, 14715–14735, https://doi.org/10.5194/acp-18-14715-2018, 2018. a

Danabasoglu, G., Lamarque, J.-F., Bacmeister, J., Bailey, D. A., DuVivier, A. K., Edwards, J., Emmons, L. K., Fasullo, J., Garcia, R., Gettelman, A., Hannay, C., Holland, M. M., Large, W. G., Lauritzen, P. H., Lawrence, D. M., Lenaerts, J. T. M., Lindsay, K., Lipscomb, W. H., Mills, M. J., Neale, R., Oleson, K. W., Otto-Bliesner, B., Phillips, A. S., Sacks, W., Tilmes, S., van Kampenhout, L., Vertenstein, M., Bertini, A., Dennis, J., Deser, C., Fischer, C., Fox-Kemper, B., Kay, J. E., Kinnison, D., Kushner, P. J., Larson, V. E., Long, M. C., Mickelson, S., Moore, J. K., Nienhouse, E., Polvani, L., Rasch, P. J., and Strand, W. G.: The Community Earth System Model Version 2 (CESM2), J. Adv. Model. Earth Sy., 12, e2019MS001916, https://doi.org/10.1029/2019MS001916, 2019. a

Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Hólm, E. V., Isaksen, L., Kållberg, P., Köhler, M., Matricardi, M., McNally, A. P., Monge‐Sanz, B. M., Morcrette, J., Park, B., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J., and Vitart, F.: The ERA‐Interim reanalysis: configuration and performance of the data assimilation system, Q. J. Roy. Meteor. Soc., 137, 553–597, https://doi.org/10.1002/qj.828, 2011. a

Emmons, L. K., Schwantes, R. H., Orlando, J. J., Tyndall, G., Kinnison, D., Lamarque, J.-F., Marsh, D., Mills, M. J., Tilmes, S., Bardeen, C., Buchholz, R. R., Conley, A., Gettelman, A., Garcia, R., Simpson, I., Blake, D. R., Meinardi, S., and Pétron, G.: The Chemistry Mechanism in the Community Earth System Model version 2 (CESM2), J. Adv. Model. Earth Sy., 12, e2019MS001882, https://doi.org/10.1029/2019MS001882, 2019. a

Errera, Q., Daerden, F., Chabrillat, S., Lambert, J. C., Lahoz, W. A., Viscardy, S., Bonjean, S., and Fonteyn, D.: 4D-Var assimilation of MIPAS chemical observations: ozone and nitrogen dioxide analyses, Atmos. Chem. Phys., 8, 6169–6187, https://doi.org/10.5194/acp-8-6169-2008, 2008. a

Errera, Q., Chabrillat, S., Christophe, Y., Debosscher, J., Hubert, D., Lahoz, W., Santee, M. L., Shiotani, M., Skachko, S., von Clarmann, T., and Walker, K.: Technical note: Reanalysis of Aura MLS chemical observations, Atmos. Chem. Phys., 19, 13647–13679, https://doi.org/10.5194/acp-19-13647-2019, 2019. a

Fischer, H., Birk, M., Blom, C., Carli, B., Carlotti, M., von Clarmann, T., Delbouille, L., Dudhia, A., Ehhalt, D., Endemann, M., Flaud, J. M., Gessner, R., Kleinert, A., Koopman, R., Langen, J., López-Puertas, M., Mosner, P., Nett, H., Oelhaf, H., Perron, G., Remedios, J., Ridolfi, M., Stiller, G., and Zander, R.: MIPAS: an instrument for atmospheric and climate research, Atmos. Chem. Phys., 8, 2151–2188, https://doi.org/10.5194/acp-8-2151-2008, 2008. a

Gettelman, A., Mills, M. J., Kinnison, D. E., Garcia, R. R., Smith, A. K., Marsh, D. R., Tilmes, S., Vitt, F., Bardeen, C. G., McInerny, J., Liu, H.-L., Solomon, S. C., Polvani, L. M., Emmons, L. K., Lamarque, J.-F., Richter, J. H., Glanville, A. S., Bacmeister, J. T., Phillips, A. S., Neale, R. B., Simpson, I. R., DuVivier, A. K., Hodzic, A., and Randel, W. J.: The whole atmosphere community climate model version 6 (WACCM6), J. Geophys. Res.-Atmos., 124, 12380–12403, https://doi.org/10.1029/2019JD030943, 2019. a

Kolmogorov, A. N.: Wiener's spiral and some other interesting curves in Hilbert space, Dokl. Akad. Nauk SSSR, 26, 115–118, 1940. a

Kraft, D.: A Software Package for Sequential Quadratic Programming, DFVLR-FB 88–28, 1988. a

Laeng, A.: Regression Coefficients of parametrisation of function of natural atmospheric variability and reparametrization on latitudinal gradients software, Institut für Meteorologie und Klimaforschung – Atmosphärische Spurenstoffe und Fernerkundung (IMK-ASF), KITopen [data set], https://doi.org/10.5445/IR/1000137514, 2021. a

Minganti, D., Chabrillat, S., Christophe, Y., Errera, Q., Abalos, M., Prignon, M., Kinnison, D. E., and Mahieu, E.: Climatological impact of the Brewer–Dobson circulation on the N_{2}O budget in WACCM, a chemical reanalysis and a CTM driven by four dynamical reanalyses, Atmos. Chem. Phys., 20, 12609–12631, https://doi.org/10.5194/acp-20-12609-2020, 2020. a

Prignon, M., Chabrillat, S., Minganti, D., O'Doherty, S., Servais, C., Stiller, G., Toon, G. C., Vollmer, M. K., and Mahieu, E.: Improved FTIR retrieval strategy for HCFC-22 (CHClF_{2}), comparisons with in situ and satellite datasets with the support of models, and determination of its long-term trend above Jungfraujoch, Atmos. Chem. Phys., 19, 12309–12324, https://doi.org/10.5194/acp-19-12309-2019, 2019. a

Rodgers, C. D. and Connor, B. J.: Intercomparison of remote sounding instruments, J. Geophys. Res., 108, 4116, https://doi.org/10.1029/2002JD002299, 2003. a

Sheese, P. E., Walker, K. A., Boone, C. D., Degenstein, D. A., Kolonjari, F., Plummer, D., Kinnison, D. E., Jöckel, P., and von Clarmann, T.: Model estimations of geophysical variability between satellite measurements of ozone profiles, Atmos. Meas. Tech., 14, 1425–1438, https://doi.org/10.5194/amt-14-1425-2021, 2021. a

Tilmes, S., Hodzic, A., Emmons, L. K., Mills, M. J., Gettelman, A., Kinnison, D. E., Park, M., Lamarque, J.-F., Vitt, F., Shrivastava, M., Campuzano-Jost, P., Jimenez, J. L., and Liu, X.: Climate forcing and trends of organic aerosols in the Community Earth System Model (CESM2), J. Adv. Model. Earth Sy., 11, 4323–4351, https://doi.org/10.1029/2019MS001827, 2019. a

von Clarmann, T.: Validation of remotely sensed profiles of atmospheric state variables: strategies and terminology, Atmos. Chem. Phys., 6, 4311–4320, https://doi.org/10.5194/acp-6-4311-2006, 2006. a

von Clarmann, T., De Clercq, C., Ridolfi, M., Höpfner, M., and Lambert, J.-C.: The horizontal resolution of MIPAS, Atmos. Meas. Tech., 2, 47–54, https://doi.org/10.5194/amt-2-47-2009, 2009. a

Yaglom, A.: Correlation Theory of Stationary and Related Random Functions, Vol. II, Springer-Verlag, ISBN 978-1461290902, 1986. a, b

^{1}

BrNO, BrO, CCl_{4}, CFC_{11}, CFC_{12}, CFC_{113}, CH_{3}Cl, CH_{4}, ClO, ClONO_{2}, CO, CO_{2}, H_{2}O, HBr, HCl, HF, HNO_{3}, HNO_{4}, HO_{2}, HOBr, HOCl, N_{2}O, N_{2}O_{5}, NO, NO_{2}, NO_{3}, O_{3}, OH and temperature

- Abstract
- Introduction
- Model data
- Variabilities
- Parametrization
- Reparametrization on latitudinal gradients
- What to do in practice: the software
- Discussion
- Conclusions
- Code and data availability
- Author contributions
- Competing interests
- Disclaimer
- Special issue statement
- Acknowledgements
- Financial support
- Review statement
- References

- Abstract
- Introduction
- Model data
- Variabilities
- Parametrization
- Reparametrization on latitudinal gradients
- What to do in practice: the software
- Discussion
- Conclusions
- Code and data availability
- Author contributions
- Competing interests
- Disclaimer
- Special issue statement
- Acknowledgements
- Financial support
- Review statement
- References