Research article 10 Aug 2021
Research article  10 Aug 2021
Slit homogenizer introduced performance gain analysis based on the Sentinel5/UVNS spectrometer
 ^{1}Meteorological Institute LMU Munich, Theresienstraße 37, Munich, Germany
 ^{2}Airbus Defence and Space GmbH, WillyMesserschmittStr. 1, 82024 Taufkirchen, Germany
 ^{3}Airbus Defence and Space GmbH, ClaudeDornierStr., 88090 Immenstaad, Germany
 ^{1}Meteorological Institute LMU Munich, Theresienstraße 37, Munich, Germany
 ^{2}Airbus Defence and Space GmbH, WillyMesserschmittStr. 1, 82024 Taufkirchen, Germany
 ^{3}Airbus Defence and Space GmbH, ClaudeDornierStr., 88090 Immenstaad, Germany
Correspondence: Timon Hummel (timon.hummel@airbus.com)
Hide author detailsCorrespondence: Timon Hummel (timon.hummel@airbus.com)
Spatially heterogeneous Earth radiance scenes affect the atmospheric composition measurements of highresolution Earth observation spectrometer missions. The scene heterogeneity creates a pseudorandom deformation of the instrument spectral response function (ISRF). The ISRF is the direct link between the forward radiative transfer model, used to retrieve the atmospheric state, and the spectra measured by the instrument. Hence, distortions of the ISRF owing to radiometric inhomogeneity of the imaged Earth scene will degrade the precision of the Level2 retrievals. Therefore, the spectral requirements of an instrument are often parameterized in the knowledge of the ISRF over nonuniform scenes in terms of shape, centroid position of the spectral channel and the full width at half maximum (FWHM).
The Sentinel5/UVNS instrument is the first pushbroom spectrometer that makes use of a concept referred to as a slit homogenizer (SH) for the mitigation of spatially nonuniform scenes. This is done by employing a spectrometer slit formed by two parallel mirrors scrambling the scene in the along track direction (ALT) and hence averaging the scene contrast only in the spectral direction. The flat mirrors do not affect imaging in the across track direction (ACT) and thus preserve the spatial information in that direction. The multiple reflections inside the SH act as coherent virtual light sources and the resulting interference pattern at the SH exit plane can be described by simulations using scalar diffraction theory. By homogenizing the slit illumination, the SH strongly modifies the spectrograph pupil illumination as a function of the input scene. In this work we investigate the impact and strength of the variations of the spectrograph pupil illumination for different scene cases and quantify the impact on the ISRF stability for different types of aberration present in the spectrograph optics.
The Ozone Monitoring Instrument (OMI) was the first instrument identifying the issue arising from nonuniform Earth scenes on the shape and maximum position of the spectral response of the instrument (Voors et al., 2006). In gratingbased imaging spectrometers, the Earth ground scene is imaged by the telescope onto the instrument entrance slit plane. The scanning over the ground area is achieved by either a scanning mirror or a pushbroom configuration, where different areas of the surface are imaged as the satellite flies forward. In the subsequent spectrograph, the slit illumination gets spectrally resolved by a dispersive element and reimaged on the focal plane array (FPA) by an imaging system. The limited spectral resolving power of the instrument arising from diffraction and aberration is described by a convolution of the slit image with the spectrometer and detector point spread function (PSF). In this study, we interpret the resulting intensity pattern on the FPA in the spectral direction as the instrument spectral response function (ISRF). In fact, there exist other definitions of the ISRF. The differentiation of the definitions become particularly important in the presence of spectrometer smile effects (Caron et al., 2017). As we neglect such effects, we will continue with the previously described definition of the ISRF.
Depending on the observed scene heterogeneity, the entrance slit will be inhomogeneously illuminated. In the case of a classical slit, this will alter the shape of the ISRF (see Fig. 1). Moreover, a scene dependency in the PSF will also affect the ISRF, which will be particularly discussed in this paper. As the ISRF is the direct link between the radiative transfer model and the spectrum measured by the instrument, a scenedependent shape of the ISRF will have an immediate impact on the accuracy of the Level2 retrieval products. Figure 2 depicts a representative topofatmosphere spectrum (solar zenith angle 10^{∘}, albedo 0.05) for the Sentinel5/UVNS (ultraviolet/visible/nearinfrared/SWIR) SWIR3 spectrometer, incident on the instrument's entrance aperture. The monochromatic spectrum will be smeared by means of a convolution with an exemplary ISRF, which depends on the imaging properties of the instrument for any given wavelength. In general, the ISRF is a wavelength and fieldofviewdependent instrument characteristic and hence varies over the FPA position. It is experimentally determined prior to launch in onground characterization campaigns. Whenever the inorbit ISRF shape deviates from the onground characterized shape due to, for example, heterogeneous scenes, it will affect the measured spectrum from which the Level2 products are retrieved (e.g., CH_{4} and CO in the SWIR3 channel of Sentinel5/UVNS).
This effect is particularly prominent for instruments with a high spatial resolution. The along track motion of the satellite during the integration times results in a temporal averaging of the ISRF variation, which reduces the impact of scene heterogeneity. The impact of, e.g., albedo variations depends on the instantaneous field of view (IFOV) and the sampling distance in ALT (for Sentinel5/UVNS: IFOV = 2.5 km, ALT SSD = 7 km). Spectrometers with a large scan area like GOME (Burrows et al., 1999) or SCIAMACHY (Bovensmann et al., 1999; Burrows et al., 1995) are less vulnerable to contrast in the Earth scene due to the small ratio between the slit footprint and the smear distance. In contrast, recent highresolution hyperspectral imaging spectrometers with an IFOV comparable to the sampling distance (or scan area) are more strongly affected and therefore demand a set of stringent requirements on the inflight knowledge and stability of the ISRF. This is necessary, as distortions in the ISRF due to nonuniform scenes will introduce biases and pseudorandom noise in the Level2 data and therefore in the precision of atmospheric composition products. For the Sentinel5 Precursor (S5P) satellite (Veefkind et al., 2012), launched in 2017 with the Tropospheric Monitoring Instrument (TROPOMI) being the single payload, Hu et al. (2016) showed that the stability and knowledge of the ISRF is the main driver of all instrument calibration errors for the retrieval accuracy. Landgraf et al. (2016) estimate the error of the retrieved CO data product due to nonuniform slit illumination to be on the order of 2 % with quasirandom characteristics. Noël et al. (2012) quantify the retrieval error for the upcoming Sentinel4 UVN imaging spectrometer for tropospheric O_{3}, NO_{2}, SO_{2} and HCHO. They identify a difference in the retrieval error depending on the trace gas under observation. The largest error occurs for NO_{2} with a mean error of 5 % and a maximum error of 50 %. They propose a software correction algorithm based on a wavelength calibration scheme individually applied to all Earth radiance spectra. As discussed by Caron et al. (2019), this type of software correction can only be applied to dedicated bands (UV, VIS, NIR) but particularly fails in the SWIR bands due to the strong absorption lines of highly variable atmospheric components.
Sentinel5/UVNS (Irizar et al., 2019) is the first pushbroom spectrometer that employs an onboard concept to mitigate the effect of nonuniform scenes in the alongtrack direction. A hardware solution called a slit homogenizer (SH) is implemented, which reduces the scene contrast of the Earth radiance in the along track direction (ALT) of the satellite flight motion by replacing the classical slit with a pair of two parallel extended mirrors (Fig. 3). The two parallel rectangular mirrors composing the entrance slit have a distance of b=248 µm, side lengths of 65 mm in ACT and a length of 9.91 mm (SWIR3) along the optical axis. Thereby, the light focused by the telescope optics onto the slit entrance plane is scrambled by multiple reflections in the ALT direction, whereas in ACT the light passes the SH without any reflection. Heterogeneous scenes in the ACT direction may also affect the ISRF stability in the presence of spectrometer smile. This effect will not be covered in this study and instead we refer the reader to Gerilowski et al. (2011) and Caron et al. (2017). For a realistic reference Earth scene of the Sentinel5/UVNS mission provided by the ESA (Fig. 5), the ISRF shall meet the requirements of <2 % ISRF shape knowledge error, <1 % relative full width at half maximum (FWHM) knowledge error and 0.0125 nm centroid error in the SWIR3. Meister et al. (2017) and Caron et al. (2019) presented simulation results providing a firstorder prediction of the performance of the SH principle, which are relevant to achieve the performance requirements above. However, so far several secondorder effects haven't been quantitatively addressed in the prediction of the homogenizing performance. This paper extends the existing firstorder models and provides a more elaborated and comprehensive description of the SH and its impact on performance and instrument layout. We present an endtoend model of the Sentinel5/UVNS SWIR3 channel (2312 nm). In particular, we determine the spectrograph pupil illumination, which is altered by the multiple reflections inside the SH. This effect changes the weighting of the aberration present in the spectrograph optics and consequently results in a scene dependency in the optical PSF. As the ISRF is not only a function of the slit illumination but also of the spectrograph PSF, a variation in the intensity distribution across the spectrograph pupil will ultimately put an uncertainty and error contribution to the ISRF. The severity of the spectrograph illumination distortion highly depends on the slit input illumination and the strength and type of aberration present in the spectrograph. In order to quantify the achievable ISRF stability, we simulate several input scenes and different types of aberration.
The outline of this paper is as follows: Sect. 2 describes the model we deployed to propagate the light through the SH by Huygens–Fresnel diffraction formula. Applying Fourier optics, we formulate the propagation of the complex electric field from the SH exit plane up to the grating position, representing the reference plane for the evaluation of the spectrograph pupil intensity distribution. In Sect. 3 we quantify the spectrograph pupil intensity distribution for several Earth scene cases. The scenedependent weighting of the aberration in the spectrograph and its impact on the ISRF properties is discussed and quantified in Sect. 4. Finally, we summarize our results in Sect. 5.
This section describes the underlying models and the working principle of the SH. The first part briefly summarizes the model developed by Meister et al. (2017), which propagates the field through the Sentinel5/UVNS instrument up to the SH exit plane by using a scalardiffraction approach. In the second part a novel modeling technique of the spectrograph optics is introduced. We put a particular focus on the scene dependency of the spectrograph illumination while using an SH.
2.1 Near field
The light from objects on the Earth that are imaged at one spatial position (along slit) within the homogenizer entrance slit arrive at the Sentinel5/UVNS telescope entrance pupil as plane waves, where the incidence angle θ is between $\pm \mathrm{0.1}{}^{\circ}$. The extent of the wavefront is limited by the size and shape of the telescope aperture. By neglecting geometrical optical aberration, the telescope would create a diffractionlimited point spread function in the telescope image where the SH entrance plane is positioned. Depending on the angle of incidence, the PSF centroid will be located at a dedicated position within the SH entrance plane. The electric field of the diffraction pattern in the SH entrance plane is given as the Fourier transform of the complex electric field over the telescope pupil. For a square entrance pupil, the diffraction pattern is calculated as (Goodman, 2005, p. 103)
where (x_{t},y_{t}) are the coordinate positions in the telescope entrance pupil and (u_{a},v_{a}) are the respective coordinates in the SH entrance plane. Ω denotes the twodimensional entrance pupil area, f is the focal length of the telescope, A the amplitude of the plane wavefront at the telescope entrance pupil, D the full side length of the quadratic telescope entrance pupil and $k=\frac{\mathrm{2}\mathit{\pi}}{\mathit{\lambda}}$ is the wavenumber. Further, the relation ${\int}_{a}^{a}{e}^{ixc}=\mathrm{2}a\mathrm{sin}c\left(ca\right)$ and a Fresnel approximation was applied in Eq. (2). The propagation of ${\stackrel{\mathrm{\u0303}}{U}}_{f}$ through the subsequent SH is described by the Huygens–Fresnel principle (Goodman, 2005, p. 66). The reflections at the two mirrors are accounted for by inverting the propagation component in ALT upon every reflection n as
where R is the reflectivity, b is the slit width and e^{inπ} describes a phase jump upon every reflection n. Inserting Eq. (2) into (3) and applying the Huygens–Fresnel diffraction principle yields the expression for the intensity distribution at the SH exit plane for a given incidence angle θ, SH length l and position $r({u}_{a},{v}_{a})=\sqrt{{l}^{\mathrm{2}}+({u}_{b}{u}_{a}{)}^{\mathrm{2}}+({v}_{b}{v}_{a}{)}^{\mathrm{2}}}$ as
where u_{b},v_{b} are the coordinates of the position at the SH exit plane. Evaluating Eq. (4) for every incidence angle of the Sentinel5/UVNS field of view (FoV) results in the socalled SH transfer function (Fig. 3b), which maps any field point originating from Earth to an intensity distribution at the SH exit plane. In a purely geometric theory and a perfect SH configuration in terms of length, every point source would be distributed homogeneously in the ALT direction (Fig. 3a). However, as is quantified in Eq. (4), the field distribution at the SH output plane highly depends on interference effects due to path differences of the reflected light inside the SH, resulting in a nonuniform transfer function as shown in Fig. 3b.
A full experimental validation of the propagation model through the SH is still missing. An initial approach to validate the model in a breadboard activity was conducted by ITO Stuttgart and published in Irizar et al. (2019).
2.2 Far field
In a spacebased imaging spectrometer equipped with a classical slit acting as a field stop, a point source on the Earth surface enters the instrument as a plane wavefront with a uniform intensity over the telescope pupil. As this principle applies for every point source in a spatial sample on the Earth, the telescope pupil intensity homogeneity is independent of the radiance variation among the point sources in a spatial sample. Besides some diffraction edge effects in the slit plane, the telescope pupil intensity distribution gets retrieved in the spectrograph pupil. This is not the case when introducing a mirrorbased SH. Existing SH models (Meister et al., 2017; Caron et al., 2019) implement the spectrometer as a simple scaling factor and the ISRF on the FPA is obtained via the convolution of the SH output intensity distribution, the pixel response implemented as a characteristic function and the spectrograph PSF. In this contribution we model the propagation through the spectrograph more accurately by including the spectrograph optics, such as the collimator, a dispersive element and the imaging optics. In particular, the inclusion of these optical parts becomes important because the SH not only homogenizes the scene contrast in the slit, but it also significantly modifies the spectrograph pupil illumination. A schematic diagram of the SH behavior and the instrument setup is shown in Fig. 4. A plane wavefront with incidence angle Θ is focused by a telescope on the SH entrance plane. In the ACT direction, the light is not affected by the SH. After a distance l, corresponding to the SH length, the diffraction limited PSF at the SH entrance plane is converted to the farfield pattern of the diffraction pattern. Independent of the applied scene in ACT, the telescope pupil intensity distribution in ACT is mostly retrieved again at the spectrograph pupil. The exact distribution of the spectrograph pupil illumination is affected by magnification factor and a truncation of the electric field at the SH entrance plane, which leads to a slight broadening and small intensity variations with a high frequency in angular space (Berlich and Harnisch, 2017). In ALT the diffraction pattern in the SH entrance plane undergoes multiple reflections on the mirrors, so that eventually the whole exit plane of the SH is illuminated. To preserve the full image information along the swath, the entrance plane of the SH must be imaged; to homogenize the scene in ALT the exit plane of the SH must be imaged. This is achieved by an astigmatism in the collimator optics. Moreover, the multiple reflections inside the SH lead to a modification of the system exit pupil illumination. In other words, the SH output plane (near field) and the spectrograph pupil intensity variation (far field) strongly depend on the initial position of the incoming plane wave and therefore on the Earth scene radiance in ALT direction. Following a first simple geometrical argument as discussed by (Caron et al., 2019), we consider a point source at the SH entrance. The rays inside the cone emerging from this source will undergo a number of reflections depending on the position of the point source and the angle of the specific ray inside the cone. The maximum angle is given by the telescope fnumber. With this geometrical reasoning it becomes obvious that the number of reflections differs among the rays inside the cone. If the number of reflections is even, a ray keeps its nominal pupil position, whereas if the number is odd, its pupil coordinate will be inverted. From this argument we deduce that the spectrograph pupil illumination will be altered with respect to the telescope pupil illumination. Note that the reallocation of the angular distribution of the light has a different origin than the remaining inhomogeneity at the SH exit plane. The achieved nearfield homogenization is dependent on the remaining interference fluctuations in the SH transfer function. In contrast, the variations in the spectrograph illumination are based on a geometrical reallocation of the angular distribution of the light exiting the SH in combination with interference effects in the spectrograph pupil plane.
In the following we make the geometrical argument rigorous using diffraction theory. A general case for the connection between slit exit plane and spectrograph pupil plane is considered by Goodman (2005, p. 104). In the scenario discussed there, a collimated input field U_{l}(x_{s},y_{s}) propagates through a perfect thin lens at a distance d. The field in the focal plane of the lens is then given by
where x_{s},y_{s} are the position in the spectrometer pupil plane and u_{b},v_{b} are the coordinates in the image plane at the SH exit. Indeed, the field at the lens focal plane is proportional to the twodimensional Fourier transform. In contrast, our situation is inverted as we are interested in U_{l}(x_{s},y_{s}), i.e., the collimated field distribution at the spectrometer pupil originating from the SH output plane. Further, we need to incorporate the astigmatism in the collimation optics and the diffraction grating. These steps are covered in the following two sections.
2.3 Collimator astigmatism
In order to keep the full image information in ACT while imaging the homogenized SH output image, the collimator needs an astigmatism. In our model, this is implemented via Zernike polynomial terms on the collimation lens. We follow the OSA/ANSI convention for the definitions of the Zernike polynomials and the indexing of the Zernike modes (Thibos et al., 2000). The focal length of the collimator in ALT images the SH exit plane, while in ACT the SH entrance plane is imaged. In the simulation this is realized with three terms: a focal length term where the focal length is that of the collimator in ALT, a defocus term to shift the object plane and an astigmatism term to separate the ALT (tangential) and ACT (sagittal) object planes.
The Zernike polynomials are given by
and
where c_{nm} are the Zernike coefficients defining the strength of the aberration and ${Z}_{n}^{m}$ are the Zernike polynomials. Due to the elegant and orthonormal definition of the Zernike polynomials, a perfect matching of defocus and astigmatism amplitude is straightforward, as the difference between the sagitta and tangential plane of the astigmatism is solely dependent on the radial term of the Zernike polynomial. Therefore, in order to match the corresponding difference given by the SH length, the weighting of the astigmatism has to be larger than the defocus term by a factor of $\sqrt{\mathrm{2}}$. Hence, the combined Zernike term will be
Including the astigmatism of the collimation optics, applying $d={f}_{\mathrm{col},\mathrm{ALT}}$ and solving Eq. (5) for U_{l,θ} by using the coordinate transformation ${x}_{s}^{\prime}=\frac{k}{f}{x}_{s}$ and ${y}_{s}^{\prime}=\frac{k}{f}{y}_{s}$, we get the field distribution at the diffraction grating as
Equation (9) yields the field distribution incident on the diffraction grating. The implementation of the diffraction grating, which is responsible for the wavelength dispersion, will be introduced in the next section.
2.4 Diffraction grating
The primary goal of the spectrometer is to distinguish the intensity of the light as a function of the wavelength and spatial position. In order to separate the wavelengths, a diffractive element is placed in the spectrograph pupil and disperses the light in the ALT direction. For our analysis, we place the diffraction grating at a distance $d={f}_{\mathrm{col},\mathrm{ALT}}$ after the collimator and on the optical axes. Further, we model the dispersive element as a onedimensional binary phase diffraction grating. Such gratings induce a π phase variation by thickness changes of the grating medium. Three design parameters are used to describe the grating and are unique for every spectrometer channel: the period of the grating Λ, the phase difference Φ between the ridge (of width d) and the groove regions of the grating, and the fill factor $d/\mathrm{\Lambda}$. Physically, the phase difference itself is induced by two parameters: the height or thickness t of the ridge and the refractive index of the material of which the grating is made. In most cases, the refractive index of the used material is fixed and the thickness of the material is the primary parameter. The phase profile with a fill factor of 0.5, which provides the maximum efficiency in the ALT direction, is given by
The complex electric field of the spectrograph pupil wavefront after the diffraction grating is then given by
The intensity distribution after the grating is given by inserting Eq. (9) in (11) and applying the absolute square:
The implementation of the diffraction grating is a simplified model, which is an approximation of the real, more complex case. In Sentinel5/UVNS, the SWIR spectrograph is equipped with a silicon immersed grating. The simplified approach is also valid for this case, as the SH does not affect the general behavior of the grating.
The farfield intensity distribution is dependent on the contrast of the Earth scene in ALT and therefore on the SH entrance plane illumination. We characterize the amplitude of the variations of the spectrograph pupil illumination by introducing two types of heterogeneous scenes. The first is an applicable Earth scene as defined by the ESA for the Sentinel5/UVNS mission, which aims at representing a realistic Earth scene case. The onground albedo variations of this scene can be parameterized as a linear interpolation between two spectra representing the same atmospheric state but obtained with either a dark or bright albedo (Caron et al., 2017). The spatial variation of the scene heterogeneity is described by introducing interpolation weights w_{k}. The resulting spectrum for a given ALT subsample k is then calculated as
where the reference spectra correspond to a tropical bright scene (L_{bright}: albedo = 0.65) and a tropical dark scene (L_{dark}: albedo = 0.05). The weighting factors that were used for this study have been derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) surface reflectance products with 500 m spatial resolution and total coverage of 25 km for relevant conditions of Sentinel5/UVNS (EOP PIO, 2011). The slit smearing due to platform movement is accounted for by convolving the onground scene with the motion boxcar of the spatial sampling distance (SSD). The platform movement acts like a lowpass filter and averages out short albedo variations with respect to the SSD and the instrument's FoV. However, without an SH, remaining inhomogeneities are present in the slit, which yield up to 20 % slit illumination variations in ALT directions. Figure 5 depicts the onground albedo contrast given in terms of weighting factors w_{k}, the scene after smearing due to the motion of the platform and the location of the SH entrance plane. We assume the scene to be homogeneous in the ACT direction. In fact, heterogeneous scenes in the ACT direction may also affect the ISRF stability in the presence of spectrometer smile (see Gerilowski et al., 2011; Caron et al., 2017). The second scene considered represents an artificial calibration (CAL) scene where 50 % of the slit is illuminated and 50 % is dark. These kind of instantaneous transitions are impossible to be observed by a pushbroom instrument with finite FoV and integration time. However, they are convenient to be applied in experimental measurements and will serve as reference to experimentally validate the SH performance models.
Figure 6 depicts the simulation results for the pupil intensity distribution in the SWIR3 (2312 nm) for the applied test scenes as well as a homogeneous slit illumination. As expected, the uniformity of the input telescope pupil illumination is completely conserved in the ACT direction due to the absence of interaction, i.e., reflection, with the SH. Therefore the tophat intensity distribution of the telescope is, besides the diffraction edge effects, completely preserved. On the contrary, the intensity distribution in ALT is dependent on the contrast of the applied scene. Even for a homogeneous scene the SH modifies the pupil intensity (Fig. 6a) and consists of symmetrical variations. The intensity pattern just varies slightly for the applicable Earth scene (Fig. 6b) due to the moderate gradient of the slit illumination variation. The CAL scenes (Fig. 6c, d) highlight the previously made geometrical argument for the nonuniform pupil illumination as parts of the pupil are left with only a fraction of the light. For illustration, we show a case where the upper 50 % of the slit are illuminated and another case where the lower 50 % of the slit are illuminated (representing the ALT illumination).
In the next section we will investigate the impact of nonuniform pupil illumination in combination with spectrograph aberration on the ISRF stability.
The main impact of the above described variations in the spectrometer pupil illumination is the scenedependent weighting of the aberration inherent to the spectrograph optics. In the case of a classical slit, it is valid to calculate the ISRF of an imaging spectrometer as the convolution of the slit illumination, the pixel response on the FPA and the optical PSF of the spectrograph optics. When using an SH, a scene dependency of the spectrograph pupil illumination will weight the aberration of the system accordingly and thereby create a variation in the PSF, which will ultimately also change the ISRF properties. Therefore, it is necessary to keep the complex phase of the electric field during the propagation through the instrument.
Instead of a convolution, we propagate the spectrograph pupil illumination through the imaging optics by diffraction integrals. For the description of the aberration present in the Sentinel5/UVNS instrument we use again the formulation of Zernike theory. We know the expected PSF size on the FPA of the Sentinel5/UVNS SWIR3 channel, which in the case of a classical slit can be approximated by the standard deviation of a normal distribution. In order to assess the impact of aberration, we impinge different types of aberration on the spectrograph imaging optics and match the PSF size to the instrument prediction. As the shape of the PSF for an arbitrary aberration is not given by a normal distribution, we define the PSF size as the area where 80 % of the encircled energy (EE) is contained. Then we tune the strength of the aberration coefficients in such a way that the size of the aberrated PSF matches that of the normal distributed PSF. For the transformation of the spectrograph pupil illumination to the FPA including aberration, we apply the thin lens formula and expand it by adding the phase term for the Zernike aberration (Goodman, 2005, p. 145). Our starting point for the propagation is the grating position where, for the case of Sentinel5/UVNS, the distance d matches the focal length of the imaging optics. In that case, the formulation simplifies again and is given by a relation that has the form of a Fourier transform:
where s,t are the coordinates at the FPA, f_{im} is the focal length of the imager, U_{g,θ} is the field distribution at the grating and H(r,ϕ), with $r=r\left({x}_{s},{y}_{s}\right)$ and $\mathit{\varphi}=\mathit{\varphi}\left({x}_{s},{y}_{s}\right)$, is the respective Zernike aberration that we apply. Any spatially incoherent monochromatic input scene can be distributed in plane wavefronts with amplitude A(Θ). Each such wavefront leads to an intensity $I={I}_{\mathrm{\Theta}}(s,t)={U}_{\mathrm{FPA},\mathit{\theta}}{}^{\mathrm{2}}$ on the FPA. As we have no SH impact in the ACT direction, we collapse this dimension and sum along it. This yields the onedimensional ISRF intensity distribution on the FPA as a function of the incidence angle Θ as I_{Θ}(t). The respective scene will weight the intensities on the FPA depending of their strength and is therefore the linear operator:
Note that for a homogeneous scene, A(Θ)=1 for every incidence angle. Finally, the normalized ISRF on the FPA is given by:
where χ is the characteristic function, which is 1 inside a pixel area and 0 elsewhere, α is a scaling factor to give the ISRF in units of wavelength (λ) and N_{σ} is the density function of a normal distribution with zero mean value and standard deviation σ. The latter factor accounts for the modulation transfer function (MTF) of the detector (not the MTF of the whole optical system).
In order to assess the stability of the ISRF we define three merit functions:

shape error: the maximum difference of the ISRF calculated for a homogeneous and heterogeneous scene, respectively:
$$\begin{array}{}\text{(18)}& \text{Shape error}:=\underset{\mathit{\lambda}}{max}{\displaystyle \frac{{\mathrm{ISRF}}_{\mathrm{hom}}\left(\mathit{\lambda}\right){\mathrm{ISRF}}_{\mathrm{het}}\left(\mathit{\lambda}\right)}{{max}_{\stackrel{\mathrm{\u0303}}{\mathit{\lambda}}}{\mathrm{ISRF}}_{\mathrm{hom}}\left(\stackrel{\mathrm{\u0303}}{\mathit{\lambda}}\right)}};\end{array}$$ 
centroid error: the shift of the position of the spectral channel centroid, where the centroid is defined as
$$\begin{array}{}\text{(19)}& \text{Centroid}:={\displaystyle \frac{\underset{\mathrm{FPA}}{\int}\mathrm{ISRF}\left(\mathit{\lambda}\right)\mathit{\lambda}\mathrm{d}\mathit{\lambda}}{\underset{\mathrm{FPA}}{\int}\mathrm{ISRF}\left(\mathit{\lambda}\right)\mathrm{d}\mathit{\lambda}}};\end{array}$$ 
the spectral resolution of the ISRF given by the FWHM.
We consider two cases for the assessment of the induced impact on the ISRF stability. In the first case, we neglect any variation of the spectrograph illumination and use the PSF as a convolution kernel of the ISRF given as a constant and sceneindependent normal distribution defined as
where σ is the standard deviation representing the size of the PSF. The spot size value for a representative field point in the SWIR3 spectrometer of Sentinel5/UVNS is about 6.85 µm. When convolving with a gaussian PSF, we neglect the nonuniformity in the pupil and the spectrometer aberration, and the ISRF errors are only driven by the slit exit illumination (near field). For the second case, we impinge a certain amount of aberration on the imaging optics to get the same spot size for the PSF as in the first case. In this case, the ISRF errors are a combination of the remaining inhomogeneities at the SH exit plane (near field) as well as effects due to nonuniform spectrograph illumination (far field). The aberration present in the Sentinel5/UVNS spectrograph are dependent on the position on the FPA in spectral and spatial direction. In the upcoming characterization and calibration campaign, the specific types of aberration of the final instrument will not be determined, but only the size of the spots. Therefore, although it is not a realistic case, we impinge pure aberration of a single type in order to determine critical Zernike terms for the ISRF stability. We also test two mixtures of different types of aberration, which represent more realistic field points of the Sentinel5/UVNS. The ISRF for a homogeneous scene including aberration will be extensively characterized on ground. We want to investigate how the ISRF based on several Zernike terms behave under the condition of nonuniform scenes and how the ISRF deviation evolves with respect to each aberrationtypespecific, homogeneous ISRF. Therefore, in the next paragraph, we calculate the relative change in the ISRF figures of merit functions.
In the following, we present the ISRF figures of merit resulting from the simulation of several Zernike polynomials for the Sentinel5/UVNS applicable heterogeneous Earth scene and a 50 % stationary calibration scene. Further, we compare the results to the case of a classical slit without scene homogenization. Tables 1 and 2 summarize the results for the ISRF figures of merit. Note that the errors for the calibration scene are much larger than the errors for a realistic Earth scene. The calibration scene can be used in a laboratory to characterize the SH performance and compare it with the prediction. All Zernike polynomials increase the error in the ISRF knowledge compared to the case where the ISRF is calculated as the convolution with a constant gaussian PSF. The error magnitude variation ranges from only small increasing errors (defocus, vertical astigmatism) to a notable increase of the error (oblique quadrafoil, horizontal coma). The aberration changes both the maximum amplitude of the errors and the specific shape of the ISRF. Figure 7b depicts the ISRF assuming pure vertical coma, pure spherical aberration and pure oblique trefoil for a heterogeneous 50 % calibration scene. The lower part of the plot shows the ISRF shape difference for each specific homogeneous reference scene. Note that the shape error is defined as the maximum amplitude of the difference plot. As none of the field points in the real Sentinel5/UVNS instrument will contain a pure singular type of aberration, we tested two sets of aberration mixtures, which is more representative of a real field point in the Sentinel5/UVNS instrument. Although our study does not provide a rigorous mathematical argument, the results indicate that the error of the combined Zernike polynomials lies within the errors of the individual contributors.
This argument is supported by Fig. 8, where we plotted the ISRF shape error going from a pure oblique quadrafoil aberration to a pure defocus aberration. In each step we reduced the fraction of the quadrafoil aberration by 20 % and tuned the defocus aberration coefficient in such a way that we ended up with the same PSF size of 6.85 µm (80 % EE). The ISRF errors always remain in the corridor between the case of pure oblique quadrafoil and pure defocus aberration. This behavior was tested for several other Zernike combinations and we conclude that the errors given in Tables 1 and 2 for the respective Zernike polynomials span the error space within which mixtures of aberration lie.
Although the phenomena of the variations of the pupil illumination in combination with spectrometer aberration increases the errors, the SH still homogenizes the scene well and significantly improves the stability of the ISRF compared to a classical slit. In Fig. 7a we compare the ISRF shape difference for a 50 % stationary calibration scene for a case with a classical slit and a case with SH. The SH improves the ISRF stability by almost an order of magnitude. Considering the applicable Earth scene and including the farfield variations, the SH still provides sufficiently stable ISRF stability with respect to the mission requirements for moderate heterogeneous scenes of Sentinel5/UVNS. This would not be the case for an instrument employed with a classical slit.
In certain scenarios, Sentinel5/UVNS will fly over Earth scenes with higher contrasts than specified in the applicable Earth scene. This will be the case when flying over cloud fields, water bodies or city to vegetation transitions. However, these scenes are excluded from the mission requirements in terms of scene homogenization. Although sufficient for the purposes of Sentinel5/UVNS, the capability of the SH to homogenize the scene is not perfect. This imperfection is particularly prominent when considering the calibration scenes. The imperfections originate from the remaining interference fluctuations in the SH transfer function and are dependent on the wavelength. Higher wavelengths show smaller frequencies and larger peaktovalley amplitudes of the maxima in the SH transfer function, which leads to reduced homogenization efficiency. Therefore, the SWIR3 wavelength channel is the most challenging in terms of scene homogenization.
We observe that increasing the number of reflections inside the SH will increase the number of stripes in the spectrometer pupil illumination (see Fig. 6c, d) and reduce the peak to valley amplitude. This would lead to a more homogeneous pupil illumination. More reflection in the SH can be achieved by either increasing the length of the SH or adapting the telescope F_{#}. However, it is advantageous to keep the SH length small to reduce the collimator astigmatism requirements. Note that a longer SH would not increase the nearfield homogenization performance. In addition, more reflections in the SH lead to greater transmission losses at the mirrors. As the errors due to the pupil illumination are small compared to achieved nearfield homogenization, it seems favorable to prioritize the firstorder design rule given in Caron et al. (2019) and Meister et al. (2017). The SH shows the best nearfield homogenization performance if ${F}_{\mathit{\#}\mathrm{tel}}=l/\left(\mathrm{2}bn\right)$, where F_{#tel} is the telescope Fnumber, l is the SH length, b is the SH width and n is the number of reflections. For Sentinel5/UVNS, the optimal parameters for SWIR3 are a telescope F_{#} of 9.95, a slit length of 9.91 mm and a slit width of 248 µm.
The simulation results of this study still require experimental validation. An initial approach to validate the SH transfer functions was published in Irizar et al. (2019), where they showed good agreement between the simulation and the experimental result for a single SH incidence angle. The verification of the full transfer function including the full FoV range is pending. The SH farfield effects investigated in this study could be determined by measuring the pupil intensity distribution at the grating position by means of an appropriate test bench. The test bench would need to be capable of illuminating the SH entrance plane through a telescope with angles representing the Sentinel5/UVNS FoV. Further, the astigmatism of the SH needs to be compensated for, which could be done by introducing a cylindrical lens in the collimator system.
Apart from the mirrorbased SH discussed in this study, future remote sensing instruments investigate the technology of another slit homogenizer technology, which is based on rectangular multimode fiber bundles. These devices are based on the same principle as the mirrorbased SH but enable one to homogenize the scene in the ACT and ALT direction (Amann et al., 2019) and provide enhanced performance over extreme albedo variations.
The presented study continues the investigation by Caron et al. (2019) and Meister et al. (2017) on mirrorbased slit homogenizer technology. While the preceding studies considered the homogenization of the SH exit plane, here we extend the models by including the electric field propagation through the subsequent spectrograph. The slit homogenizer not only homogenizes the slit illumination but also modifies the spectrograph illumination dependent on the input scene. The variations in the spectrograph pupil illumination will lead to a scenedependent weighting of the geometrical aberration in the optical system, which causes an additional distortion source of the ISRF. The phenomena is particularly prominent in the presence of extreme onground albedo contrasts. This will be the case when the instrument flies over clouds or water bodies. However, in the context of the Sentinel5/UVNS instrument, these scenes are excluded from the mission requirements.
We observe that the impact of spectrograph pupil illumination variations is small compared to the error due to nonuniform slit illumination, and the ISRF distortion is primarily driven by the remaining nearfield variations after the SH. The inhomogeneity remnants arise from the fluctuations of the interference pattern at the SH exit plane. The strength of the variations increases with wavelength. Therefore, this study was conducted in the SWIR3 channel in order to cover the worst case.
We quantify the ISRF in terms of shape error, FWHM error and centroid error at 2312 nm by an endtoend propagation through the SH and the subsequent spectrograph optics. With regard to these figures of merit, our simulation results suggest an increase of the errors depending on the specific type of aberration impinged on the optics. ISRF errors of combined Zernike polynomials are always within the maximum errors of the individual Zernike constituent. Although the SH changes the spectrometer illumination, it still has significant performance advantages in stabilizing the ISRF compared to a classical slit. For an applicable heterogeneous Earth scene, the SH improves the ISRF shape stability by a factor of 5–10. The remaining residual errors are well below the Sentinel5/UVNS system requirements, which are 2 % shape error, <1 % relative FWHM error and <0.0125 nm (SWIR3) centroid error.
The datasets generated and/or analyzed for this work are available from the corresponding author on reasonable request, subject to confirmation of Airbus Defence and Space GmbH.
TH developed, implemented, applied and evaluated the methods for the endtoend modelling of the Sentinel5/UVNS instrument with input from CM. TH performed the performance gain analysis which was supported by CM and CK and revised by JK and MW. TH prepared the paper, with contributions and critical revision from all coauthors.
The authors declare that they have no conflict of interest.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
We thank Tobias Lamour, Jess Köhler and Markus Melf (all at Airbus Defence and Space) for the helpful comments on a previous version of the manuscript.
This paper was edited by Ulrich Platt and reviewed by Bernd Sierk and Ruediger Lang.
Amann, S., DuongEderer, Q., Haist, T., Sierk, B., Guldimann, B., and Osten, W.: Characterization of fiberbased slit homogenizer devices in the NIR and SWIR, in: International Conference on Space Optics – ICSO 2018, 7–10 October 2014, La Caleta, Tenerife, Canary Islands, edited by: Sodnik, Z., Karafolas, N., and Cugny, B., vol. 11180, 2276–2286, International Society for Optics and Photonics, SPIE, https://doi.org/10.1117/12.2536147, 2019. a
Berlich, R. and Harnisch, B.: Radiometric assessment method for diffraction effects in hyperspectral imagers applied to the earth explorer 8 mission candidate flex, in: International Conference on Space Optics – ICSO 2014, 7–10 October 2014, La Caleta, Tenerife, Canary Islands, edited by: Sodnik, Z., Cugny, B., and Karafolas, N., vol. 10563, 1475–1483, International Society for Optics and Photonics, SPIE, https://doi.org/10.1117/12.2304079, 2017. a
Bovensmann, H., Burrows, J. P., Buchwitz, M., Frerick, J., Noël, S., Rozanov, V. V., Chance, K. V., and Goede, A. P. H.: SCIAMACHY: Mission Objectives and Measurement Modes, J. Atmos. Sci., 56, 127–150, https://doi.org/10.1175/15200469(1999)056<0127:SMOAMM>2.0.CO;2, 1999. a
Burrows, J., Hölzle, E., Goede, A., Visser, H., and Fricke, W.: SCIAMACHY–scanning imaging absorption spectrometer for atmospheric chartography, Acta Astronaut., 35, 445–451, https://doi.org/10.1016/00945765(94)00278T, 1995. a
Burrows, J. P., Weber, M., Buchwitz, M., Rozanov, V., LadstätterWeißenmayer, A., Richter, A., DeBeek, R., Hoogen, R., Bramstedt, K., Eichmann, K.U., Eisinger, M., and Perner, D.: The Global Ozone Monitoring Experiment (GOME): Mission Concept and First Scientific Results, J. Atmos. Sci., 56, 151–175, https://doi.org/10.1175/15200469(1999)056<0151:TGOMEG>2.0.CO;2, 1999. a
Caron, J., Sierk, B., Bezy, J.L., Loescher, A., and Meijer, Y.: The CarbonSat candidate mission: radiometric and spectral performances over spatially heterogeneous scenes, in: International Conference on Space Optics – ICSO 2014, 7–10 October 2014, La Caleta, Tenerife, Canary Islands, edited by: Sodnik, Z., Cugny, B., and Karafolas, N., vol. 10563, 1019–1027, International Society for Optics and Photonics, SPIE, https://doi.org/10.1117/12.2304186, 2017. a, b, c, d
Caron, J., Kruizinga, B., and Vink, R.: Slit homogenizers for Earth observation spectrometers: overview on performance, present and future designs, in: International Conference on Space Optics – ICSO 2018, 9–12 October 2018, Chania, Greece, edited by: Sodnik, Z., Karafolas, N., and Cugny, B., vol. 11180, 402–417, International Society for Optics and Photonics, SPIE, https://doi.org/10.1117/12.2535957, 2019. a, b, c, d, e, f
EOPPIO: IPDRSESA18, Issue 1, Revision 2, 10 February 2011, Sentinel5UVNS Instrument Phase A/B1 Reference Spectra, ESA, Noordwijk, the Netherlands, 2011. a
Gerilowski, K., Tretner, A., Krings, T., Buchwitz, M., Bertagnolio, P. P., Belemezov, F., Erzinger, J., Burrows, J. P., and Bovensmann, H.: MAMAP – a new spectrometer system for columnaveraged methane and carbon dioxide observations from aircraft: instrument description and performance analysis, Atmos. Meas. Tech., 4, 215–243, https://doi.org/10.5194/amt42152011, 2011. a, b
Goodman, J. W.: Introduction to Fourier optics, vol. 1, McGrawHill, San Francisco, CA, USA, 2005. a, b, c, d
Hu, H., Hasekamp, O., Butz, A., Galli, A., Landgraf, J., Aan de Brugh, J., Borsdorff, T., Scheepmaker, R., and Aben, I.: The operational methane retrieval algorithm for TROPOMI, Atmos. Meas. Tech., 9, 5423–5440, https://doi.org/10.5194/amt954232016, 2016. a
Irizar, J., Melf, M., Bartsch, P., Koehler, J., Weiss, S., Greinacher, R., Erdmann, M., Kirschner, V., Albinana, A. P., and Martin, D.: Sentinel5/UVNS, in: International Conference on Space Optics – ICSO 2018, 9–12 October 2018, Chania, Greece, edited by: Sodnik, Z., Karafolas, N., and Cugny, B., vol. 11180, 41–58, International Society for Optics and Photonics, SPIE, https://doi.org/10.1117/12.2535923, 2019. a, b, c
Landgraf, J., aan de Brugh, J., Scheepmaker, R., Borsdorff, T., Hu, H., Houweling, S., Butz, A., Aben, I., and Hasekamp, O.: Carbon monoxide total column retrievals from TROPOMI shortwave infrared measurements, Atmos. Meas. Tech., 9, 4955–4975, https://doi.org/10.5194/amt949552016, 2016. a
Meister, C., Bauer, M., Keim, C., and Irizar, J.: Sentinel5/UVNS instrument: the principle ability of a slit homogenizer to reduce scene contrast for earth observation spectrometer, in: SPIE Proceedings, Sensors, Systems, and NextGeneration Satellites XXI, Vol. 10423, 104231E, https://doi.org/10.1117/12.2278619, 2017. a, b, c, d, e
Noël, S., Bramstedt, K., Bovensmann, H., Gerilowski, K., Burrows, J. P., Standfuss, C., Dufour, E., and Veihelmann, B.: Quantification and mitigation of the impact of scene inhomogeneity on Sentinel4 UVN UVVIS retrievals, Atmos. Meas. Tech., 5, 1319–1331, https://doi.org/10.5194/amt513192012, 2012. a
Thibos, L. N., Applegate, R. A., Schwiegerling, J. T., and Webb, R.: Standards for Reporting the Optical Aberrations of Eyes, in: Vision Science and its Applications, Optical Society of America, SuC1, https://doi.org/10.1364/VSIA.2000.SuC1, 2000. a
Veefkind, J., Aben, I., McMullan, K., Förster, H., de Vries, J., Otter, G., Claas, J., Eskes, H., de Haan, J., Kleipool, Q., van Weele, M., Hasekamp, O., Hoogeveen, R., Landgraf, J., Snel, R., Tol, P., Ingmann, P., Voors, R., Kruizinga, B., Vink, R., Visser, H., and Levelt, P.: TROPOMI on the ESA Sentinel5 Precursor: A GMES mission for global observations of the atmospheric composition for climate, air quality and ozone layer applications, Remote Sens. Environ., 120, 70–83, https://doi.org/10.1016/j.rse.2011.09.027, 2012. a
Voors, R., Dobber, M., Dirksen, R., and Levelt, P.: Method of calibration to correct for cloudinduced wavelength shifts in the Aura satellite's Ozone Monitoring Instrument, Appl. Opt., 45, 3652–3658, https://doi.org/10.1364/AO.45.003652, 2006. a