Retrievals of X CO 2 , X CH 4 and X CO from portable, near-infrared Fourier transform spectrometer solar observations in Antarctica.

. The Collaborative Carbon Column Observing Network (COCCON), uses low-resolution, portable EM27/SUN Fourier Transform Spectrometers (FTSs) to make retrievals of column-averaged dry-air mole fractions (DMFs, represented as X gas ) of CO 2 , CH 4 , CO and H 2 O from near infrared solar absorption spectra. The COCCON has developed rapidly over recent years and complements the Total Carbon Column Observing Network (TCCON). In this work we provide details of the first seasonal timeseries of near infrared X CO 2 , X CH 4 and X CO retrievals from 5 measurements made in Antarctica during the deployment of an EM27/SUN to the Arrival Heights laboratory on Ross Island over the austral summer of 2019/20 under the auspices of the COCCON. The DMFs of all three species were lower in Antarctica than at mid-latitude and for X CO 2 and X CO the retrieved values were less variable. For X CH 4 however, the variability was significantly greater and it was found that this was strongly correlated to the proximity of the polar vortex. 10

The Bruker EM27/SUN is a portable, low-resolution Fourier transform spectrometer (FTS) with a built-in solar tracker that can measure near-infrared, solar absorption spectra. From these spectra, it is possible to retrieve columnaveraged DMFs (represented as X gas ) of CO 2 , CH 4 , CO and H 2 O with a precision and accuracy similar to or better than the TCCON (Gisi et al., 2012 andHase et al., 2016). A network based on the EM27/SUNs is being developed, known as the Collaborative Carbon Column Observing Network (COCCON, Frey et al., 2019).
The comparatively low cost, portability and relative easeof-use of these instruments have meant that they can be utilised in greater numbers or to cover specific targets in locations where the more complex instruments cannot be deployed. For example, Velazco et al. (2019) used an EM27/SUN in a semiarid region of Australia to validate retrievals from GOSAT, Knapp et al. (2021) have conducted observations from a ship transiting the Pacific Ocean and Frey et al. (2021) have established a COCCON site in Namibia while Hase et al. (2015) and Dietrich et al. (2021) have used dense networks of these instruments to estimate the carbon fluxes of Berlin and Munich, respectively. Tu et al. (2022) have gone on to demonstrate the utility of these instruments to quantify fluxes at the facility level.
In this work, we present the data gathered during the deployment of an EM27/SUN to the Arrival Heights laboratory on Ross Island, Antarctica (77.83 • S, 166.66 • E, 205 m a.m.s.l.) over the austral summer of 2019-2020. The retrieved time series of X gases is available from the COCCON central facility hosted by the ESA Atmospheric Validation Data Centre (EVDC) https://doi.org/10.48477/coccon.pf10.arrivalheights.R02 (Pollard, 2021). This dataset represents the first seasonal time series of X gases retrieved from near-infrared solar spectra in Antarctica and will provide a useful source of validation data at these high, southern latitudes which are not well covered by existing networks (the southernmost TCCON station is Lauder at 45 • South).
In the next section, we will introduce the EM27/SUN instrument, briefly describe the retrieval scheme used to derive X gases from the measured solar spectra, and describe the other datasets that we will compare the EM27/SUN results with. The following sections will describe the measurements made at Arrival Heights, discuss the results and demonstrate the robustness and utility of the data. Conclusions will be drawn in Sect. 4.

Instrumentation and data processing
In this section, we will briefly describe the EM27/SUN instrument, followed by the retrieval scheme used to infer the column-averaged dry-air mole fractions of trace gases from the measured solar absorption spectra.
We will also provide a high-level overview of the TC-CON data from the Lauder site, which were compared to the EM27/SUN retrievals before and after the period of deployment, and Sentinel-5 precursor (S5P) data that were compared to the data collected at Arrival Heights.

EM27/SUN Fourier transform spectrometer
The Bruker EM27/SUN is a portable, low-resolution Fourier transform spectrometer (FTS) with a built-in solar tracker using a camera trained on the detector aperture to provide active feedback. The instruments measure DC-coupled interferograms with a spectral resolution of 0.5 cm −1 using two indium gallium arsenide (InGaAs) detectors at room temperature, one with a spectral range of 5500-11 000 cm −1 and a second, wavelength-extended detector measuring in the 4000-5500 cm −1 range (Gisi et al., 2012 andHase et al., 2016).

EM27/SUN data processing
The raw interferograms are processed and X gas retrievals made using the PROFFAST software of the COCCON-PROCEEDS framework which was developed on behalf of the ESA; it is open-source and freely available. It has previously been described by Sha et al. (2020) and will be summarised herein. The DC-coupled, double-sided interferograms are processed into spectra by the PROFFAST-PREPROCESS code which includes a DC correction, phase correction and quality checks of the resulting spectra. Separately, cross sections of the molecular absorption are calculated for each day by the PCXS module based on the meteorological and trace gas priors generated using the TCCON method and for a representative surface pressure. The trace gas retrievals are then computed by the INVERS module which uses a least squares fitting algorithm to scale the prior profiles, with the option of adjusting the surface pressure to that which was actually measured at the time of the observation.
A further correction is applied to the retrieved X CO 2 values. This is to account for a compounded correction factor applied by PROFFAST, which had been derived for an earlier version of the code. The earlier version included a bug that caused a pointer offset of magnitude of 1 in the spectrum handling. This bug was fixed in the PROFFAST version 2020-08-10, and while it was found that the correction factors applied to all other retrieved species remained valid, it was necessary to apply a further, solar zenith angledependent correction to X CO 2 of the following form: where X CO 2 corr and X CO 2 are the corrected and uncorrected values and SZA is the solar zenith angle in degrees. This correction has been applied to the R02 version of the dataset described in this work. A previous dataset version (R01) was published before the need to apply the correction was known. Further details of this correction and the background to it are available in "Technical note on X CO 2 bias in current PROFFAST distribution" (https://www.imk-asf.kit.edu/ english/3225.php, last access: 27 January 2022).

TCCON
The Lauder TCCON station (45.04 • S, 169.68 • E, 370 m a.m.s.l.) is the southernmost in the network and has one of the longest continuous records. The station has been previously described in Pollard et al. (2017Pollard et al. ( , 2021. In this work, we have used the standard output from the GGG2014 data-processing version (Wunch et al., 2015) which is publicly available from the TCCON data archive (Pollard et al., 2019).

Sentinel-5 precursor
The ESA's Sentinel-5 precursor (S5P) satellite, which carries the TROPOMI instrument as its sole payload, is orbiting in a Sun-synchronous, low-Earth polar orbit and has an observation swath of 2600 km wide across track, resulting in daily global coverage and a pixel size of 5.5 × 7 km for CH 4 and CO.
The S5P operational X CH 4 data are retrieved using the RemoteTeC-S5P algorithm , which produces retrievals of X CH 4 only under cloud-free conditions. In this work, we compared both the standard product and a bias-corrected version of the dataset. The details of the bias correction can be found in the algorithm theoretical baseline document (ATBD) for S5P methane retrievals (Hasekamp et al., 2021). The S5P operational total column density of carbon monoxide (CO) is retrieved using a short-wave infrared carbon monoxide retrieval (SICOR) algorithm . The retrievals of CO are performed simultaneously with interfering trace gases and effective cloud parameters, such as cloud height and optical thickness. The offline (OFFL) operational data version 01.03.02 was used in this paper.

Methods and results
In this section, we will first demonstrate the instrumental stability of both the EM27/SUN and the TCCON instrument at the Lauder site through the instrument line shape (ILS) retrieved for both.
In the subsequent two subsections, the experimental setups used for the measurements made at Arrival Heights and at the Lauder site before and after the deployment will be described along with a discussion of the retrievals resulting from those measurements.
Comparisons to the S5P retrievals will also be examined to demonstrate the application of this dataset to satellite validation activities.

Instrument stability
In order to monitor the stability of the alignment of an FTS, the ILS can be retrieved (Hase et al., 2013).
For the Lauder TCCON instrument, regular (approximately monthly) measurements are made of an internal cell, containing a calibrated quantity of HCl, illuminated by a lamp source. From these spectra, the ILS is retrieved using version 14.5 of the LINEFIT software described in Hase et al. (2013).
The EM27/SUN does not have the capability to measure HCl cells, but an ILS retrieval can be achieved by taking long (4 m) path measurements of a lamp source and making use of water vapour adsorption as described in Frey et al. (2015) and updated by Alberti et al. (2022). For the EM27/SUN instrument used in this study (serial number 53), ILS measurements were conducted at the Karlsruhe Institute of Technology before and after the instrument was shipped to New Zealand.
A time series of the results of these ILS retrievals, parameterised in terms of the modulation efficiency (ME) at maximum optical path difference (OPD) and the maximum value of the phase error, are shown in Fig. 1 for both instruments. A noticeable feature of Fig. 1 is a step change in the ME at max OPD for the TCCON instrument in May 2019. This is related to a change of the instrument's metrology laser on 28 May. Before this change, the mean ME at max OPD was 0.9991 (with a standard deviation of 0.0006), after the change, it was 0.9902 (0.0010). This 1 % change in ME at max OPD will not have a significant impact on X gas retrievals since Hase et al. (2013) estimated that a 4 % change would result in only a 0.035 % error in X CO 2 . The TCCON max phase error remained virtually the same across this change at −0.0019 (0.0004) before and −0.0022 (0.0007) after. The EM27/SUN had an ME at max OPD and max phase error of 0.9891 and 0.0006, respectively before leaving Karlsruhe and 0.9875 and 0.0010 when it returned. We therefore conclude that both instruments maintained their alignment throughout the presented dataset.

EM27/SUN observations and retrievals at Arrival Heights
At Arrival Heights, the EM27/SUN was positioned on a bench inside the laboratory. A small amount of the beam from the solar tracker (Robinson et al., 2020) that is used for a Bruker IFS125HR as part of the Network for the Detection of Atmospheric Composition Change (NDACC, De Mazière Under normal operations, 10 interferogram scans over a period of about 70 s are co-averaged into a single EM27/SUN measurement. To halve the amount of data needing to be transferred, it was decided to co-average 20 interferograms at Arrival Heights. Since the rate-of-change in solar zenith angle is relatively small at high latitudes, this change is not expected to have any impact on retrievals. The pressure data used in the retrievals at Arrival Heights is taken from the weather station at the nearby Scott Base (NIWA electronic weather station, EWS, Scott Base, ID:12740, 77.85 • S, 166.76 • E, 20 m a.m.s.l.) and corrected for the altitude difference of 185 m between the two locations. Figure 2 shows the time series of each of the X gases retrieved from the EM27/SUN and the Lauder TCCON station before, during and after the measurement campaign at Arrival Heights.
The X CO 2 measured at Arrival Heights is systematically lower than at Lauder. This is an expected result as highlatitude southern hemispheric CO 2 concentrations are less than at mid-latitudes (Stephens et al., 2013). Also, the sea-sonality that can be seen at the lower latitude is not obvious in these measurements. However, over the period of deployment, the growth rate of X CO 2 is also not seen at Arrival Heights. This is potentially because the seasonal drawdown is of a similar magnitude and obscures the growth over this short period.
The X CH 4 shows considerable variability over the period of deployment and a negative trend which contradicts the positive trend seen at lower latitudes. To investigate this structure further, we examined the isentropic modified potential vorticity (MPV) (Lait, 1994) over Arrival Heights derived from the Modern Era Retrospective-Analysis for Research and Applications reanalysis product (MERRA2) (Gelaro et al., 2017). Although not an absolute diagnostic of the position of the polar vortex relative to Arrival Heights, the MPV value will be negatively correlated with the influence of the vortex . The top panel of Fig. 3 shows the MPV at the 460 K isentropic level (corresponding to the lower stratosphere) over the course of the instrument campaign, while the lower panel shows the daily averaged X CH 4 retrievals for the same period. There is clear correlation between the two (r = 0.82, 95 % CI: 0.68-0.90). The correlation of X CH 4 with MPV confirms a weak barrier effect of the polar vortex (Choi et al., 2002).
The lower panel of Fig. 2 shows X CO . The data collected at Arrival Heights are clearly measuring baseline X CO concentrations whilst at lower latitudes, the Lauder TCCON data Figure 2. Time series of daily averaged X CO 2 (a), X CH 4 (b) and X CO (c) retrieved from spectra measured using the Lauder TCCON instrument (blue circles), the EM27/SUN whilst it was at Lauder (red squares) and when it was at the Arrival Heights laboratory (green triangles).
show spikes caused by stratospheric transport of air masses affected by biomass burning in the tropics.

Comparison to TCCON
Before and after the EM27/SUN was deployed to Arrival Heights, it was operated for a period at Lauder, alongside the TCCON station there. Between 8 March and 22 July 2019, the instrument was in a workshop where it could be moved outside on fine days to make observations using the built-in solar tracker. An automatic scheduling application (Geddes et al., 2018) was used to run the measurements and this was able to interrogate the output from the EM27/SUN's Camtracker software and pause measurements when clouds obscured the Sun. On 22 July, the EM27/SUN was relocated to the same laboratory as the TCCON instrument; here, the built-in solar tracker was parked and illuminated by a small amount (approximately 5 %) of the parallel solar beam from a solar tracker coupled to another Bruker 125HR. The difference in altitude between these locations is 10 m and corrections have been applied to account for this in the surface pressure used for trace gas retrievals. This pressure measurement is from the Vaisala PTB100A sensor that is part of the Lauder climate station.
After returning to New Zealand from Antarctica, the instrument was returned to Lauder where it was installed un-der another, permanently installed solar tracker and further measurements were collected from 5 June until 3 September 2020 when it was returned to the Karlsruhe Institute of Technology (KIT).
In total, measurements were collected on 72 d alongside the Lauder TCCON station. To make a meaningful comparison between the EM27/SUN and TCCON, we first average the retrievals from both instruments into bins of 10 min. The window of 10 min. was chosen to include sufficient measurements to reduce the effects of random uncertainties whilst not aliasing in slowly varying signals due to e.g. air mass dependence. Figure 4 shows a time series of the differences between the 10 min averages before and after the deployment for the three species, with box and whisker plots of these in the inset panel. A summary of the before-and-after statistics is also given in Table 1. There is a small, but statistically significant, difference for all three species between the beforeand-after comparisons. As described previously, retrievals of the instrument line shape (ILS) for both instruments spanning the period in question show no large drifts in the alignment of either. Therefore, it is reasonable to expect that this is due to a seasonally dependent bias between the two instruments. While it is difficult to draw conclusions from this limited time span plotted in Fig. 4, it is not unreasonable to suggest that this is the case. Such an effect was seen previ-  ously by Sha et al. (2020) when comparing low-and highresolution instruments. This effect is likely caused by the averaging kernels of the low-and high-resolution instruments aliasing different errors in the common priors into the respective retrievals. In any case, the change in offset with respect to TCCON is 0.07 % for X CO 2 (0.16 % for X CH 4 and 1.72 % for X CO ), which is below the 0.1 % precision value generally accepted as the requirement for carbon cycle studies and satellite validation (Wunch et al., 2015).

Comparison to S5P
The spatial and temporal coincidence criteria used to select S5P measurements corresponding to EM27/SUN observations are the same as Sha et al. (2021), i.e. within a radius of 100 km (50 km) around the site for methane (carbon monox-ide) validation and with a maximal time difference of 1 h for EM27/SUN observations around the S5P overpass time. Because of the high latitude of the Arrival Heights laboratory and the inclination of the S5P orbit, this often results in coincidences on more than one S5P orbit for a particular day's EM27/SUN observing period. The average S5P pixel values are compared to EM27/SUN retrievals which have had the a priori alignment applied to compensate for or correct its contribution to the smoothing equation (Rodgers and Connor, 2003). The co-located pairs are selected only if a minimum of five S5P pixels were found in applying the co-incidence criteria. Figure 5 shows the comparison between the ground-based EM27/SUN and the S5P retrievals for bias-corrected X CH 4 and standard X CO (left and right plots, respectively) with all values retrieved by each instrument (pale symbols) as well as co-located overpass means (bold) (upper panels) and the relative differences (lower panels).
The bias-corrected S5P X CH 4 product shows a bias of 2.14 %, which exceeds the S5P bias requirement (1.5 %). The bias of the standard S5P X CH 4 product without the albedodependent correction is 1.05 %. This shows that the S5P X CH 4 product at this location is strongly dependent on the surface albedo and the respective correction applied to the product. The standard deviation of the relative bias, which is a measure of the random error, is below 0.3 % for both standard and bias-corrected S5P X CH 4 products, which is well below the requirement of 1 %. This high bias at Arrival Heights was also seen by Sha et al. (2021) when they per- Figure 4. Time series of the difference between 10 min averages of (a) X CO 2 , (b) X CH 4 and (c) X CO retrieved by the TCCON and EM27/SUN instruments at Lauder before and after the deployment to Arrival Heights. The inset panel shows box and whisker plots of these differences. formed a comparison with the NDACC station there. This was attributed to the highly variable topography in the region. Similar effects have also been noted at high northern latitudes (Lambert et al., 2021, p. 129), again attributed to variability on surface albedo and topography. These results further highlight the value of having reliable, high-quality, ground-based measurements which are independent of surface effects to validate satellite retrievals, particularly in challenging regions such as these.
For S5P X CO , the mean bias is 3.77 %, which is well within the S5P bias requirement (15 %). The standard deviation of the relative bias is 4.73 %, which is well below the requirement of 10 %. Applying a cone co-location criterion to the S5P data following the ground-based EM27/SUN line of sight, as recommended in Sha et al. (2021) for highlatitude sites, we find a mean bias of 5.89 % compared to the EM27/SUN. The bias with respect to the EM27/SUN is less than the bias of 11.99 % that Sha et al. (2021) found when comparing to the NDACC data using a common a priori for a period of about 3 years. The discrepancy is mostly because of the different time periods between the NDACC and EM27/SUN comparisons, a priori difference or spectroscopy limitations. The S5P bias shows a seasonality with respect to NDACC data with high values in September and October, a period when the EM27/SUN did not measure. Furthermore, a bias change of about +3.7 % has been reported for the S5P validation with NDACC when comparing directly versus using the S5P a priori as the common prior, whereas this change is about −0.26 % for EM27/SUN comparisons.

Conclusions
We have described the first seasonal time series of X CO 2 , X CH 4 and X CO retrieved from near-infrared solar spectra in Antarctica, gathered during the deployment of an EM27/SUN to the Arrival Heights laboratory on Ross Island, Antarctica over the austral summer of 2019-2020 under the auspices of the COCCON.
Through monitoring of the ILS and a comparison to the Lauder TCCON station before and after the deployment, we have demonstrated that the precision of the X CO 2 retrievals was 0.07 % (0.16 % for X CH 4 and 1.72 % for X CO ) over the duration of the measurement campaign, and that the instrument was unaffected by being shipped to and from Antarctica.
The retrieved column-averaged abundances of all three measured species (X CO 2 , X CH 4 and X CO ) were lower in Antarctica than at the mid-latitude Lauder TCCON station as expected. However, the range of X CH 4 values observed at Arrival Heights was larger than at Lauder and are well correlated with the proximity of the polar vortex edge.
When comparing the EM27/SUN retrievals at Arrival Heights to S5P, it was found that the S5P, bias-corrected X CH 4 product had a mean difference of 2.14 % which exceeded the mission bias requirement of 1.5 %. However, the product without the albedo-dependent bias correction only differed by 1.05 %, suggesting that the albedo-dependent bias correction is not valid under these surface conditions. This finding is consistent with previous studies and highlights the value of these high-quality, ground-based measurements which are independent of surface effects to validate satellite retrievals. For S5P X CO , the mean bias is 3.77 %.
It is expected that further deployments of EM27/SUN instruments to Arrival Heights will be undertaken in the future and the data added to the COCCON archive.
Author contributions. DFP operated the EM27/SUN at Lauder, installed it at Arrival Heights, conducted the data processing and analysis and wrote this manuscript. MS performed the comparison with S5P, FH supplied the EM27/SUN and developed the PROF-FAST code. DS provided expertise on instrument deployment at Arrival Heights, derivation of the MPV time series and assisted in manuscript preparation. CA and DD performed EM27/SUN calibrations and ILS retrievals at Karlsruhe Institute of Technology (KIT). Review statement. This paper was edited by David Carlson and reviewed by two anonymous referees.