Linearization of the effect of slit function changes 1 for improving OMI ozone profile retrievals

9 We introduce a method that accounts for errors caused by the slit function in an optimal estimation 10 based spectral fitting process to improve ozone profile retrievals from the Ozone Monitoring Instrument 11 (OMI) ultraviolet measurements (270-330 nm). Previously, a slit function was parameterized as a standard 12 Gaussian by fitting the Full Width at Half Maximum (FWHM) of the slit function from climatological OMI 13 solar irradiances. This cannot account for the temporal variation of slit function in irradiance, the intra-orbit 14 changes due to thermally-induced change and scene inhomogeneity, and potential differences in the slit 15 functions of irradiance and radiance measurements. As a result, radiance simulation errors may be induced 16 due to convolving reference spectra with incorrect slit functions. To better represent the shape of the slit 17 functions, we implement a more generic super Gaussian slit function with two free parameters (slit width 18 and shape factor); it becomes standard Gaussian when the shape factor is fixed to be 2. The effects of errors 19 in slit function parameters on radiance spectra, referred as “Pseudo Absorbers (PAs)”, are linearized by 20 convolving high-resolution cross sections or simulated radiances with the partial derivatives of the slit 21 function with respect to the slit parameters. The PAs are included in the spectral fitting scaled by fitting 22 coefficients that are iteratively adjusted as elements of the state vector along with ozone and other fitting 23 parameters. The fitting coefficients vary with cross-track and along-track pixels and show sensitivity to 24 heterogeneous scenes. The PA spectrum is quite similar in the Hartley band below 310 nm for both standard 25 and super Gaussians, but is more distinctly structured in the Huggins band above 310 nm with the use of 26 super Gaussian slit functions. Finally, we demonstrate that some spikes of fitting residuals are slightly 27 smoothed by accounting for the slit function errors. Comparisons with ozonesondes demonstrate noticeable 28 improvements when using PAs for both standard and super Gaussians, especially for reducing the 29 systematic biases in the tropics and mid-latitudes (mean biases of tropospheric column ozone reduced from 30 -1.4 ~ 0.7 DU to 0.0 ~0.4 DU) and reducing the standard deviations of tropospheric ozone column 31 differences at high-latitudes (by 1 DU for the super Gaussian). Including PAs also makes the retrievals 32 consistent between standard and super Gaussians. This study corroborates the slit function differences 33 between radiance and irradiance demonstrating that it is important to account for such differences in the 34 ozone profile retrievals. 35


Introduction
The fitting of measured spectra to simulated spectra is the most basic concept for analysis of the Earth's atmospheric constituents from satellite measurements.Therefore, accurate calibration and simulation of measurements are essential for the successful retrieval of atmospheric constituents.The knowledge of the instrumental spectral response function (ISRF) or slit function could affect the accuracies of both calibration and simulation, as it is required for the convolution of a high-resolution reference spectrum to instrument's spectral resolution in the wavelength calibration and for the convolution of high-resolution absorption cross section spectra or simulated radiance spectra in the calculation of radiance at instrumental resolution.
Compared to other trace gases, the retrieval of ozone profiles can be more susceptible to the accuracy of ISRFs due to the large spectral range, where the radiance spans a few orders of magnitude and to the fact that the spectral fingerprint for the tropospheric ozone is primarily provided by the 310-330 nm absorption features residing in the temperature-dependent Huggins bands.Therefore, the efforts to characterize and verif the ISRFs have preceded the analyses of ozone profiles from satellite and aircraft measurements (Liu et al., 2005(Liu et al., , 2010;;Cai et al., 2012;Liu et al., 2015;Sun et al. 2017;Bak et al., 2017).
For space-borne instruments, ISRFs are typically characterized as a function of the detector dimensions using a tunable laser source prior to the launch (Dirksen et al., 2006;Liu et al., 2015;van Hees et al., 2018) and directly used in ozone profile retrievals (e.g., Kroon et al., 2011;Mielonen et al., 2015;Fu et al., 2013;2018).However, the preflight measured ISRFs could be inconsistent with those after launch due to the orbital movement and the instrument temperature change (Beirle et al., 2017;Sun et al., 2017).Therefore, the post-launch ISRFs have been fitted from the preflight ones (e.g., Bak et al., 2017;Sun et al., 2017) or parameterized through a cross-correlation of the measured solar irradiance to a high-resolution solar spectrum (Caspar and Chance, 1997), assuming Gaussian-like shapes (e.g., Liu et al. 2005;2010).The direct retrieval of the ISRFs from radiances has not typically been done due to the complication of taking the atmospheric trace gas absorption and Ring effect into account in the cross-correlation procedure and the slow-down of the fitting process.However, slit function differences between radiance and irradiance could exist due to scene heterogeneity, differences in stray light between radiance and irradiance, and intraorbit instrumental changes (such as instrument temperature change) (Beirle et al., 2017;Sun et al., 2017).
In addition, using temporally invariant slit functions derived from climatological solar spectra in the retrievals could cause the long-term trend errors if instrument degradation occurs.Therefore, there is room for improving our trace gas retrievals by accounting for the effects of the different ISRFs between radiance and irradiance on the spectral fitting on a pixel-to-pixel basis.The "Pseudo Absorber (PA)" is a common concept in spectral fitting to account for the effect of physical phenomena that are difficult or computationally demanding to be simulated in radiative transfer calculations, like spectral misalignments (shift and stretch) between radiance and irradiance, Ring effect, spectral undersampling, and additive straylight offsets.The pseudo absorption spectrum can be derived from a finite-different scheme (e.g.Azam and Richter, 2015) or a linearization scheme via a Taylor expansion (e.g.Beirle et al., 2013;2017); the latter approach is more efficient than the former one, but less accurate because only the first term of the Taylor series is typically taken into account for simplicity.Beirle et al. (2013) introduced a linearization scheme to account for spectral misalignments between radiance and irradiance and then included them as a pseudoabsorber in DOAS-based NO 2 and BrO fittings.Similarly, Beirle et al. (2017) linearized the effect of the change of the ISRF parameterized as a super Gaussian on GOME-2 solar irradiance spectra to characterize the slit function change over time and wavelength.Sun et al. (2017) derived on-orbit slit functions from solar irradiance spectra measured by the Ozone Monitoring Instrument (OMI) (Levelt et al., 2006) assuming standard Gaussian, super Gaussian, and preflight ISRFs with adjusted widths.The derived on-orbit slit functions, showing significant cross-track dependence that cannot be represented by preflight ISRFs, substantially improve the retrievals by the Smithsonian Astrophysical Observatory (SAO) ozone profile algorithm.However, it is not fully understood why the use of super Gaussian or stretched preflight functions, which are supposed to better model the OMI spectra as indicated by smaller mean fitting residuals, does not improve the retrievals over the use of standard Gaussian especially in the standard deviations of the differences with relative to ozonesonde observations.This study suggests that the slit functions derived from solar spectra might not fully represent those in radiance spectra.
As such, the objective of this paper is to expand the slit function linearization proposed by Beirle et al. (2017) into the optimal estimation based spectral fitting of the SAO ozone profile algorithm.The slit function linearization is used to account for the radiative transfer calculation errors caused by the slit functions differences between radiance and irradiance on a pixel-by-pixel basis, and ultimately to improve OMI ozone profile retrievals.This paper is organized as follows: after a mathematical description of the linearization of slit function changes using the generic super Gaussian function, we introduce their practical application in an optimal estimation based spectral fit procedure (Section 2).This linearization scheme is implemented differently, depending on the simulation scheme of measured spectra using high resolution or effective cross section data, respectively.Section 3 characterizes the derived pseudo absorber spectra, along with evaluations of ozone profile retrievals using independent ozonesonde observations as a reference dataset.Finally, the summary of this study is given in Section 4.

Super Gaussian linearization
The slit function parameterization and linearization are briefly summarized as in Beirle et al. (2017), focusing on what we need to derive the pseudo absorbers in the terms of the optimal estimation based fitting process.The slit function can be parameterized with the slit width , and shape factor  assuming the supper Gaussian, S as: where A(, ) is  (2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015) of the fitted slit width and shape factor in 310-330 nm are displayed in Figure 1.a.The FWHM and shape factor of the super Gaussian function is on average 0.44 nm and 2.9, respectively, while the FWHM of the standard Gaussian is 0.395 nm.The sharp change and random-noise of these derived slit function parameters might be influenced by the decreasing signal-to-noise ratio (SNR) of solar spectra later in the OMI mission and radiometric errors in solar irradiance due to the row anomaly (Sun et al., 2017).
and thus the effect of changes of S on the convolved high-resolution spectrum can be parameterized as where the convolved spectrum is I = S ⨂  ℎ .Consequently, the partial derivatives of I with respect to slit parameters p are defined as while the response of the unit change of the slit width to the convolved spectrum is dominant against that of the shape factor.

Implementation of the slit function linearization in the SAO ozone profile algorithm
In Beirle et al. (2017) a slit function linearization was implemented only to fit solar irradiances from GOME-2.We implement the slit function linearization to fit radiances in the SAO ozone profile algorithm (Liu et al. 2010), which is routinely being performed to produce the OMI PROFOZ product (https://avdc.gsfc.nasa.gov/index.php?site=1389025893&id=74).Two spectral windows (270-309 nm in the UV1 band and 312-330 nm in the UV2 band) are employed to retrieve ozone profiles from OMI BUV measurements.To match the different spatial resolutions between UV1 and UV2 bands, every two crosstrack pixels are averaged for UV2 band, resulting into 30 positions with the spatial resolution of 48 km (across-track) × 13 km (along-track) at nadir position.Partial ozone columns at 24 layers between the surface and 60 km are iteratively estimated toward minimizing the fitting residuals between measured and simulated radiances and simultaneously between a priori and estimated ozone values using the optimal estimation inversion method.A priori ozone information is taken from a tropopause-based (TB) ozone profile climatology (Bak et al., 2013).The Vector Linearized Discrete Ordinate Radiative Transfer model (VLIDORT; Spurr, 2008) is used to simulate the radiances and their derivatives with respect to geophysical parameters.The radiance calculation is made for the Rayleigh atmosphere, where the incoming sunlight is simply absorbed by ozone and other trace gases, scattered by air molecules, and reflected by surfaces/clouds assumed as a Lambertian surface.Besides these, other physical phenomena are treated as PAs to the spectral response such as Ring effect, additive offset, and spectral shifts due to misalignments of radiance relative to irradiance and ozone cross sections.In the SAO algorithm, these PAs are derived using finite differences of the radiances with and without perturbation to a phenomenon, except for the Ring spectrum that is calculated using a first-order single scattering rotational Raman scattering model (Sioris and Evans, 2000).
In this paper, we introduce new PAs to account for the radiance simulation errors caused by the slit function errors.The OMI ISRFs have been parameterized as a standard Gaussian from climatological OMI solar irradiances for each UV1 and UV2 band and thereby these PAs could take into account the spectral fitting responses caused by temporal variations of the slit function.This ozone fitting procedure uses ISRFs to convolve high resolution absorption spectra, taken from Brion et al. (1993) for ozone absorption cross sections and Wilmouth et al. (1999) for BrO absorption cross sections.In DOAS analysis, the pseudo absorber is defined as ∂S  ⨂ ℎ (  ℎ is a high-resolution absorption cross section), which could be calculated at a computationally low-cost.In our optimal estimation based ozone profile retrievals, it is conceptually defined as ∂S  ⨂ ℎ (  ℎ is a high-resolution simulated radiance), which is computationaly very expensive because of on-line radiative calculation for a ~ 60 nm wide fit window on the spatial pixel-to-pixel basis.
We now introduce how to implement the slit function linearization to derive the derivatives of the OMI radiances with respect to slit function changes in two different radiative transfer approaches used in the SAO ozone profile algorithm, i.e., the effective cross section approach in Liu et al (2010) and the updated high-resolution convolution approach described in Kim et al. ( 2013), respectively.
In Liu et al (2010), VLIDORT simulates the radiances at OMI spectral grids (λ omi ) using effective cross sections that are produced by convolving high-resolution cross sections with the OMI ISRFs.
Therefore, we apply a similar convolution process of matching the high-resolution cross section spectra with OMI spectra to derive the partial derivative of   with respect to slit parameter, p as follows: where  ,ℎ is a high-resolution absorption spectrum for ozone or BrO.Due to the dominant absorption of O 3 over BrO, the derivative of the BrO cross section with respect to p is neglected here.This partial derivative of ozone is then converted to the partial derivative of radiance through the chain rule with the analytical ozone weighting function ( ), calculated from VLIDORT, as follows: This simulation process is hereafter referred to as "effective resolution cross section (ER) simulation".
As described in Kim et al. (2013), the radiative transfer calculation in the SAO ozone profile algorithm has been performed using high-resolution extinction spectra at the optimized sampling intervals for resolving the ozone absorption features, which are a ~1.0 nm below 300 nm and ~0.4 nm above 300 nm.
These sampling intervals are coarser than actual OMI sampling grids with approximately half the number of wavelengths.The coarser sampled simulated radiances are then interpolated to a fine grid of 0.05 nm assisted by the weighting functions with respect to absorption and Rayleigh optical depth: where ∆   and ∆   are the optical thickness (the product of cross section and layer column density) at each layer for trace gas absorption and Rayleigh scattering, respectively.The convolution is then applied to these simulated high-resolution radiances, I(λ h ) with assumed slit functions and derivatives, respectively, and thereby I(λ omi ) and Furthermore, this linearization process can be formulated with n-order polynomial fitting parameters (∆ i ) to account for the wavelength-dependent change of the slit parameters around a central wavelength λ ̅ , which is expressed as . ( 9)

Results and Discussion
We characterize the effect of including the PA ( ∂lnI ∂ • △ ) on ozone profile retrievals using both Super Gaussian and standard Gaussian slit functions.Hereafter, the correction spectrum ( ) is derived using the HR simulation.The PA coefficient (∆ i ) (one for each channel and for each order) is included as part of the state vector to be iteratively and simultaneously retrieved with ozone.The a priori value is set to be zero for all fitting coefficients, while the a priori error is set to be 0.1, empirically.We should note that the empirical "soft calibration" is applied to OMI radiances before the spectral fitting, in order to eliminate the wavelength and cross-track dependent systematic biases, due to the interference of the PA coefficients with systematic measurement errors during the fitting process.Gaussian due to the interference between shape factor and slit width.

Characterization of the pseudo absorbers in ozone fitting procedure
Examples of the PAs (eq.9) are illustrated in Figure 5 when (a) zero and (b) first-order polynomial coefficients are fitted, respectively.In the UV1 range, the sum of PAs multiplied by corresponding coefficients, regardless of which Gaussian is assumed as slit function, is very similar because the spectral structure caused by the slit width change is dominant.It implies that OMI ISRFs in the UV1 band are similar to the standard Gaussian, for both radiance and irradiance measurements, consistent with the prelaunch characterization (Dirksen et al., 2006).However, in the UV2 range, the spectral structures are generated by the shape factor change rather than the slit width change and therefore PAs show noticeable discrepancies for different Gaussian assumptions.Our results indicate that the PA for the shape factor change is required to adjust the spectral structures due to the differences in the slit functions between radiance and irradiance over the UV2 band.In the case of the wavelength dependent PA coefficient fit, the impact of first-order PAs on OMI radiances is relatively visible in the wavelength range of 300-310 nm.
This result is physically consistent with the wavelength dependent property shown in the slit parameters derived from OMI irradiances as shown in Figure 6 where slit parameters are characterized in 10-pixel increments assuming the super Gaussian slit function.In UV1, the slit widths plotted as FWHM slightly decrease by ~ 0.1 nm at shorter wavelengths than 288 nm, but vary more sharply by up to ~ 0.2 nm at longer wavelengths.Compared to slit widths, the wavelength dependences of the shape factors are less noticeable, except at boundaries of the window.In the UV2 window, both slit width and shape factor are highly invariant.

Impact of including pseudo absorbers on ozone profile retrievals
Figures 7 to 9 evaluate the impact of including zero-order PAs on ozone profile retrievals.Figure 7 illustrates how different assumptions in the slit functions affect the ozone profile retrievals with respect to the retrieval sensitivity and the fitting accuracy from the case shown in The fitting accuracy is assessed as the root mean square (RMS) of the relative differences (%) between measured and calculated radiances over the UV1 and UV2 ranges, respectively.Including the PAs makes little difference in the UV1 fitting residuals for most of individual pixels (1-5 %), but significantly reduces residuals in the UV2 range.The adjusted amounts of the residuals with PAs are generally larger when assuming super Gaussian slit functions.This comes from different assumptions for slit functions in deriving soft calibration spectra, where slit functions were parameterized as standard Gaussians.Therefore, applying soft calibration to OMI spectra entails somewhat artificial spectral structures if ISRFs are assumed as super Gaussian in ozone retrievals, and hence the impact of PAs on the spectral fitting becomes more considerable.
Figure 8 compares how the spectral residuals are adjusted with PAs when soft calibration is turned on and off, respectively.Using super Gaussians causes larger amplitudes of the spectral fitting residuals than using standard Gaussians, if soft calibration is turned on and PAs are excluded.On the other hand, some residuals are reduced and more broadly structured if soft calibration is turned off.Including PAs eliminates or reduces some spikes of fitting residuals as well as improves the consistency of the fitting accuracy between using standard and super Gaussians at wavelengths above 300 nm.
The benefit of this implementation on ozone retrievals is further assessed through comparison with Electrochemical Concentration Cell (ECC) ozonesondes collected from the WOUDC (https://woudc.org/) and SHADOZ (https://tropo.gsfc.nasa.gov/shadoz/)networks.This evaluation is limited to the period of 2005 through 2008 to avoid interferences with row-anomaly effects appearing in 2007 and becoming serious in early 2009 (Schenkeveld, et al 2017).We select 13 SHADOZ sites in the tropics and 38 WOUDC sites in the northern mid/high latitudes.The collocation criteria is within +/-1 ° in latitude and longitude and within 12 hours in time.For comparison, high-vertical resolution (~100 nm) profiles of ozonesondes are interpolated onto OMI retrieval grids (~2.5 km thick).We limit OMI/ozonesonde comparisons to OMI solar zenith angle < 85°, effective cloud fraction < 0.4, surface albedo < 20 % (100 %) in tropics and midlatitudes (high latitude), top altitude of ozonesondes > 30 km, ozonesonde correction factors ranging from 0.85 to 1.15 if they exist, and data gaps for each ozonesonde no greater than 3km.Comparisons between OMI and ozonesondes are performed for the tropospheric ozone columns (TCOs) over 3 different latitude bands and for ozone profiles including all the sites, with and without PAs (zero-order) for standard and super Gaussian slit function changes, respectively.
In Table 1, the comparison statistics of tropospheric ozone columns between OMI and ozonesonde are summarized as a function of latitude bands.Without using PAs, the comparison results show a noticeable discrepancy in mean biases (1.3-2.1 DU or 3.9-6.4%)due to different assumptions on the slit function shape, with positive biases of 0.3-0.7 DU for super Gaussians and negative biases of 1.0-1.4DU for standard Gaussians.Overall, OMI retrievals are in a better agreement with ozonesonde measurements using super Gaussians.The correlations and standard deviations are very similar in the tropics and mid-latitudes, but the retrievals with standard Gaussians show better correlation and smaller standard deviations at highlatitudes.As in Sun et al. (2017), the retrievals show significant differences between using standard and super Gaussians, although there are some inconsistencies in comparing OMI and ozonesondes; the main inconsistent factors are: In this study, soft calibration is turned on and a priori information is taken from the TB climatology to perform OMI ozone profile retrievals, whereas soft calibration is turned off and a priori information is taken from the LLM climatology in Sun et al. (2017).OMI/ozonesonde data filtering criteria are quite similar to each other, except that the criteria of the solar zenith angle and cloud fraction are relaxed from 75° and 0.3 to 85° and 0.4, respectively, and the adjustment of ozonesondes with correction factors given for the WOUDC dataset is turned on in this study.Comparison is performed by latitudes here whereas global comparison is analyzed in Sun et al. (2017).After accounting for the slit differences between radiances and irradiances using PAs, the retrievals are improved for both standard and super Gaussians and these two retrievals become consistent except for the use of super Gaussians in the tropics.The mean biases in the tropics and mid-latitudes are almost eliminated, but the standard deviations and correlation do not change much.In the high-latitudes, the standard deviations and correlation are significantly improved due to applying PAs with super Gaussian ISRFs.The lack of improvement with PAs in the tropics with super Gaussians illustrates that ISRFs of radiances are quite similar to those of irradiances in the tropics, while super Gaussians better parameterize OMI ISRFs than standard Gaussians.This is consistent with the comparison of the fitting accuracy of the UV2 band as shown in Figure 7, where the fitting residuals are slightly reduced in the tropics when super Gaussians are linearized, but the standard Gaussian linearization significantly improves the fitting accuracy.The mean biases of the profile comparison as shown in Figure 9 clearly show that including PAs to account for ISRF differences reduces mean biases by up to ~ 5 % below 10 km and their general altitude dependence, and improves the consistency between using standard and super Gaussians; in addition, the standard deviations are slightly improved in the 10-20 km altitude range for both Gaussians.The improvement at all latitudes corroborates the change of ISRFs between radiance and irradiance along the orbit as conjectured by Sun et al. (2017).The consistency between using standard and super Gaussians after using PAs is mainly because there is strong anti-correlation between the slit width and shape partial derivatives as shown in Figure 2, so the adjustment of slit width only in the use of standard Gaussians can achieve almost the same effect as the adjustment of both parameters in the use of super Gaussians.Accounting for the wavelength dependent change of the ISRFs with first-order PAs makes insignificant differences to both fit residuals and ozone retrievals (not shown here).This could be mainly explained by the fact of negligible wavelength dependence of OMI ISRFs especially in UV2 as shown in Figure 5, where the PA spectrum (  ln • ∆) shows almost no variance except at the upper boundary of UV1, as well as in Figure 6 where the UV2 slit parameters derived from irradiances in the subfit windows vary within 0.05 nm for FWHM and 0.2 for shape factor.

Summary
The knowledge of the Instrument Spectral Response Functions (ISRFs) or slit functions is important for ozone profile retrievals from the Hartley and Huggins bands.ISRFs can be measured in the laboratory prior to launch, but they have been typically derived from solar irradiance measurements assuming Gaussian-like functions in order to account for the effect of the ISRF changes after launch.However, the parameterization of the ISRFs from solar irradiances could be inadequate for achieving a high accuracy of the fitting residuals as ISRFs in radiances could significantly deviate from those in solar radiances (Beirle et al., 2017) and might affect ozone profile retrievals as suggested in Sun et al. (2017).Therefore, this study implements a linearization scheme to account for the spectral errors caused by the ISRF changes as Pseudo Absorbers (PAs) in an optimal estimation based fitting procedure for retrieving ozone profiles from OMI BUV measurements using the SAO ozone profile algorithm.The ISRFs are assumed to be the generic super Gaussian that can be used as standard Gaussian when fixing the shape factor to 2. This linearization was originally introduced in Beirle et al. (2017) for DOAS analysis, but this study extends this application and more detail how to implement in practice using two different approaches to derive radiance errors from slit function partial derivatives with respect to slit parameters.These two approaches correspond to the two methods of simulating radiances at instrument spectral resolution, one using effective cross sections which were previously used in the SAO ozone profile algorithm and are still used in most of the trace gas retrievals from the UV and visible, and the other calculating radiances at high resolution before convolution, which is the preferred method in the SAO ozone profile algorithm.Consistent PAs are derived with these two approaches, as expected.
The fitting coefficients (△ ) to the PAs, representing the difference of slit parameters between radiance and irradiance, are iteratively fitted as part of the state vector along with ozone and other parameters.The UV1 slit parameters show distinct cross-track-dependent differences, especially in high latitudes.In Insignificant wavelength dependence on OMI slit functions is demonstrated from slit function parameters derived from irradiances in the sub-fit window, which leads to little difference in ozone profile retrievals when zero and first-order wavelength dependent PA coefficients are implemented to fit the spectral structures caused by slit function errors, respectively.Therefore we evaluate the benefit of including the zero-order PAs fit on both the accuracy of the fitting residuals and the quality of retrieved ozone profiles through validation against ozonesonde observations.Some spikes in the fitting residuals are reduced or eliminated.Commonly, including PAs makes little change on both fit residuals and ozone retrievals in the tropics if super Gaussians are assumed as ISRFs but this is not the case for the standard Gaussian assumption.In the TCO comparison between OMI and ozonesonde, the mean biases are reduced by 0.2 (0.6) DU and 0.6 (1.4) DU in the tropics (mid-latitude) when super and standard Gaussians are linearized, respectively.In particular, applying PA improves the standard deviations at high latitudes by 1.0 DU for super Gaussian and 0.5 DU for standard Gaussian.The profile comparison generally shows improvements in mean biases (~ 5% in the lower troposphere) as well as in standard deviation, slightly in the altitude range 10-20 km by applying PAs.More importantly, using these PAs make the retrieval consistent between standard and super Gaussians.Such consistency is due to the anti-correlation between slit width and shape PAs.This study demonstrates the slit function differences between radiance and irradiance and their usefulness to account for such differences on a pixel-to-pixel basis.In this experiment, the soft spectrum, derived with the standard Gaussian assumption, is applied to remove systematic measurement errors before spectral fitting, indicating that the evaluation of ozone retrievals might be unfairly performed for the super Gaussian function implementation.Nonetheless, OMI ozone profile retrievals show better agreement with ozonesonde observations when the super Gaussian is linearized.
Actually, the fitting residuals are slightly more broadly structured with super Gaussians than with standard Gaussians if the soft-calibration and PAs are turned off, indicating the benefit of deriving a soft calibration with the super Gaussians.Therefore, there is still room for achieving better benefits when using the PAs on ozone profile retrievals by applying the soft calibration derived with super Gaussians.
Author contributions.JB and XL designed the research; XL provided oversight and guidance; KC and JB developed the methodology together; JB conducted the research and wrote the paper; XL, SK, KC, and JHK contributed to the analysis and writing.
Competing interests.The authors declare that they have no conflict of interest.
representing the gamma function.This equation allows many forms of distributions by varying : the top-peaked function (k<2), the standard Gaussian function (k=2), and the flat-topped function (k>2). is converted to the Full Width at Half Maximum (FWHM) via the relationship of FWHM = 2 √2   .We investigate the impact of including one more slit parameter k on the OMI ISRF fitting results over the standard Gaussian using OMI daily solar measurements.As an example, time-series

Figure 1 .
Figure 1.b illustrates the high wavelength stability (0.003 nm) in the OMI mission, verifying that better calibration stability is performed with super Gaussian slit functions as abnormal deviations of wavelength shifts are derived with standard Gaussian slit functions.
∂lnI ∂ is calculated.This simulation process is hereafter referred to as "highresolution cross section (HR) simulation."The ER simulation is more commonly implemented in trace gas retrievals in the UV and visible, but the HR simulation allows for more accurate fitting residuals, to better than 0.1 %(Kim et al., 2013)  as well as shorter computation time.∂lnI ∂ is scaled by the fitting coefficients, ∆, to account for the actual size of the spectral structures caused by the slit function differences between radiance and irradiance spectra.The "pseudo absorber (PA)" for the super Gaussian slit function linearization is expressed as: the form of the logarithm of normalized radiance, PA is physically related to the optical depth change ∆τ. Figure 3 compares the partial derivatives of radiances to slit parameters,   in HR and ER simulations.Little difference is found even though convolution error for ozone cross sections is only accounted for in the ER simulation due to the overwhelming impact of ozone cross section convolution errors over other cross section data.The amplitude of   varies with different satellite pixels (e.g., ozone profile shape, geometry, and cloud/surface property), but the spectral peak positions do not change because they arise from the errors due to the convolution process of high-resolution absorption cross-sections dominated by ozone.It should be noted that these spectral structures are weakly correlated with the partial derivatives of radiances with respect to other state vectors (ozone, BrO, cloud fraction, surface albedo, radiance/irradiance shift, radiance/ozone cross section shift, Ring, mean fitting residual scaling factor) within ± 0.3 and ± 0.1 in the UV 1 and UV 2, respectively.

Figure 4
Figure4displays how the zero-order PA coefficients (∆) vary within one orbit when slit functions are assumed as standard and Super Gaussians, respectively, along with variation of cloud fraction, surface albedo, and cloud pressure from the retrievals.These retrieved coefficients physically represent the deviation of ISRFs in radiances from those in solar measurements.We normalize them with the slit parameters derived from OMI solar irradiances for a better interpretation.Cross-track dependent features are shown in slit width.The relative change of the slit width is more distinct in the UV1 band than in the UV2 band, whereas the change of the shape factor is more distinct in the UV2 band.The UV2 slit widths increase typically within 5 % over the given spatial domain.However, the UV1 slit widths increase from 10 % at most pixels up to 50 % at off-nadir positions in the high latitudes, which might be caused by stray light differences between radiance and irradiance and intra-orbit instrumental changes.An abnormal change of the UV1 slit parameters due to the scene heterogeneity is detected at the along-track scan positions of ~300 and 900, respectively, where upper-level clouds are present.The UV2 shape factor changes show a coherent sensitivity to bright surfaces under clear-sky condition over the northern high latitudes.Fitting coefficients for the standard Gaussian show a quite similar spatial variation for the UV1 slit width

Figure 4 .
In this figure, the Degrees of Freedom for Signal (DFS) represents the independent pieces of ozone information available from measurements, which typically decreases as ozone retrievals are further constrained by other fitting variables.The reduced DFS values (< 5 %) imply that the ozone retrievals are correlated slightly with PAs.
addition, an abnormal △ p caused by scene heterogeneity is observed around bright surfaces and cloudy scenes.The PA spectrum ( ∂I ∂ •△ ) illustrates that the slit width change causes most of the spectral structures in the UV1 band because the OMI ISRFs are close to Gaussian.Otherwise, the ISRF change results into different spectral responses in the UV2 band with different Gaussian functions because the adjustment of the shape factor becomes more important in accounting for the convolution error when using super Gaussians.

Figure 1 .
Figure 1.Time series of (a) slit parameters and (b) wavelength shifts for OMI daily irradiance measurements (310-330 nm) at nadir cross track position when super Gaussians (solid line) and standard Gaussians (dotted line) are parameterized as slit function shapes, respectively.

Figure 3 .
Figure 3. Derivatives of an OMI radiance spectrum simulated using high-resolution (HR) and effective resolution (ER) cross section spectra with respect to slit parameters assuming a super Gaussian function./ is multiplied by a factor of 10 to visually match / on the same y-axis.

Figure 4 .Figure
Figure 4. Pseudo absorption coefficients (∆, ∆) for fitting OMI radiances to account for slit function changes assuming (a) standard Gaussian and (b-c) super Gaussian, for the first orbit of measurements on 1 July 2006, with (d-f) the corresponding geophysical parameters.∆ and ∆ are displayed after being normalized with   , and   , the slit parameters derived from OMI solar irradiance measurements.

Figure 6 .
Figure 6.OMI ISRF FWHM (nm) and shape factor () as functions of the center wavelength, as derived from OMI solar irradiances assuming Super Gaussian functions over a range of 31 spectral pixels in 10pixel increments.Different colors represent different cross-track positions from 1 (blue) to 30 (red).

Figure 7 .
Figure 7. Same as Figure 4, but for comparisons of the Degrees of Freedom for Signal (DFS) and the Root Mean Square (RMS) of spectral fitting residuals in UV 1 and UV2 with and without zero-order pseudo absorber.Positive values indicate that both fitting residuals and DFSs are reduced due to the pseudo absorber.

Figure 8 .
Figure 8.Average differences (%) between measured (OMI) and simulated (VLIDORT) radiances at the nadir cross-track pixel in the tropics (30°S-30°S) from measurements used in Figure 4, without (a) and with (b) zero-order pseudo absorbers (PA) when the standard Gaussian (black line) and the super Gaussian (red line) are assumed as ISRFs, respectively.Upper/lower panels represent the fitting results with soft calibration being turned on/off.The residuals in the UV1 (< 310 nm) are scaled by a factor of 2 to fit in the given y-axis.In the legend, the RMS of residuals (%) are given for UV1 and UV2 wavelength ranges, respectively.