《1. Introduction》

1. Introduction

Air pollution has been one of China’s greatest challenges in recent decades [1,2]. In particular, megacity clusters in east China have endured heavy air pollution due to intensive anthropogenic activities [3–9]. In recent years, however, air pollution across China has dramatically decreased due to a series of measures taken to reduce emissions [10,11]. Nevertheless, the atmospheric pollution over highly industrialized and densely populated east China is still severe [12,13]. The severity, expanse, and complexity of the atmospheric pollution in this region remain unrivaled in comparison with the rest of the world [11,13–15]. Long-term observations of the volume mixing ratio (VMR) profiles and total columns of key atmospheric constituents are necessary to understand the sources, chemical mechanisms, and transport processes of air pollution and carbon emissions in China, as well as for potential regulatory and control purposes [15]. In addition, precise detection of the burdens and variabilities of key tropospheric (e.g., CO, CO2, CH4, N2O, C2H6, SF6, carbonyl sulfide (OCS), H2CO, HCN, and peroxyacetyl nitrate (PAN)) and stratospheric (e.g., O3, HNO3, HF, HCl, and ClONO2) species and an understanding of their influence on climate change are important for establishing scientific links between climate change and atmospheric constituents; validating satellite, airborne, or other ground-based observations; supporting processfocused field campaigns; and improving theoretical chemical transport models (CTMs) [16–19]. 

Ground-based high-resolution direct solar Fourier-transform spectrometry (FTS) has been established as a powerful tool to derive VMR profiles and total columns of atmospheric constituents [16,17,20–23]. Both the Total Carbon Column Observing Network (TCCON) and the Network for Detection of Atmospheric Composition Change (NDACC) utilize high-resolution FTS spectrometers to observe the variabilities of atmospheric constituents. The TCCON and NDACC networks have been operating since 2004 and 1992, respectively, and provide time series of total columns and VMR profiles of more than 30 atmospheric constituents. These results have been widely used in scientific investigations of the carbon cycle [24–30], atmospheric pollution and transport [8,10,13,31–34], and stratospheric ozone (O3) [18]; the development of spectroscopic line lists and models [35–45]; the validation of satellite observations and satellite retrieval algorithm development [46–56]; and the evaluation of atmospheric CTMs [57–61]. However, most TCCON/NDACC stations are operated in North America and Europe, while operational stations in the rest of the world remain sparse. For now, there is only one qualified operational site within China: Hefei station (117.2°E, 32.0°N, 30.0 m above sea level (a.s.l.)) [8,15]. 

In this study, we provide an overview of the first ground-based high-resolution FTS observation station in China. The site is close to the Yangtze River Delta (YRD) in east China—one of the most highly industrialized and densely populated regions in the country, with complex pollution and carbon emission sources. We introduce the overall envisaged concept of the observation station, describe the site and instrumentation, and summarize the major achievements that have been made at this station, which include spectral retrieval characterization and harmonization, investigation of the overall characteristics of key atmospheric constituents, emission estimates, satellite and CTM evaluations, and results on the sources and transport patterns of the studied atmospheric constituents. An outlook of the envisaged plan for observations, scientific studies, and data usage at the Hefei station is also presented.

《2. Site description and instrumentation》

2. Site description and instrumentation

《2.1. Envisaged concept and site description》

2.1. Envisaged concept and site description

The FTS station is located on Science Island in the suburbs of Hefei in highly industrialized and densely populated east China [62,63]. The Hefei station has run TCCON observations since July 2014 and became a formal TCCON station in 2018 . This station is not yet affiliated with the NDACC network, but its observation routine has followed the NDACC standard conventions since 2015 [8,17]. The FTS observatory is a cross-disciplinary observation platform developed by the Anhui Institute of Optics and Fine Mechanics, Chinese Academy of Sciences (AIOFM-CAS), in collaboration with the University of Science and Technology of China (USTC), and other TCCON/NDACC partners such as the University of Wollongong, Australia, and the University of Bremen, Germany. The observatory is part of the Atmospheric Environmental Observation and Simulation (AEOS) infrastructure, which is a National Mega-Project of Science Research of China; it includes many kinds of environmental instruments, simulators, and platforms for intensive atmospheric science research. 

 https://tccondata.org/.

Hefei is populated by 7 million people, and the city center is located to the southeast of the station. In the other directions, the observation station is surrounded by cultivated lands or wetlands. The regional landscape around the station is mostly flat, with a few hills. The observatory is the only qualified station in China that has been performing long-term FTS observations of key atmospheric constituents; thus, it plays a crucial role in evaluating satellite observations, CTM simulations, air quality, and the long-range transport of atmospheric pollutants arising from anthropogenic and natural emissions in this important region [63].

Benefiting from its unique geographical location, the Hefei station is an ideal site for investigating the physical processes of air masses from different origins and the atmospheric chemistry interactions along the transport paths. Hefei is the starting point of a series of city clusters along the Hefei-to-Shanghai axis in the YRD region. It is generally downwind of the regional anthropogenic and natural sources under the southeasterly wind in summer and northeasterly wind in winter. As a result, its location is ideal for investigating the physical and chemical processes of atmospheric constituents after sufficient mixing and chemical reactions [64]. Hefei is also located at the southern tip of the North China Plain, which contains another highly industrialized and densely populated region in China: the Beijing–Tianjin–Hebei (BTH) region. The northerly winter monsoon carries regional plumes from intensive fossil fuel and biomass burning (BB) sources in the BTH region to the Hefei station [13]. Furthermore, due to the influence of the summer monsoon in the warm season, the Hefei station is basically downwind of the tropical forest region in Southeast Asia (SEAS) and the subtropical forest region in southwest China, which have strong biogenic emissions. However, from a sub-synoptic perspective, the air masses arriving at Hefei in all seasons are influenced by cyclones or anti-cyclones originating from different regions [15,64].

《2.2. Instrumentation》

2.2. Instrumentation

The observatory includes a high-resolution FTS spectrometer (IFS125HR, Bruker, Germany), a solar tracker, and a meteorological station (Fig. 1) [62,63]. The spectrometer is inside a room, whereas the solar tracker and weather station are installed directly on the building roof. Two air conditioners, a dehumidifier, and an air cleaner are operated uninterruptedly to obtain a constant temperature with dry and clean indoor conditions [62,63]. IFS125HR spectrometers exhibit outstanding capabilities in precision and stability during normal operations and have been widely used within the TCCON/NDACC networks [22,65]. The IFS125HR spectrometer at the Hefei station has six detectors, four beam splitters, and a maximum optical path difference (OPD) of 937.5 cm, corresponding to a maximum possible resolution of 0.001 cm–1 . All detector and beam splitter types, along with their work ranges, are summarized in Tables S1 and S2 in Appendix A, respectively. The spectrometer covers a wide spectral range from 420 to 50 000 cm–1 . The solar tracker is mounted inside a dome that consists of two spherical covers, which protects the solar tracker from bad weather such as rain, snow, and high wind speed. The rotation of the dome and the on/off state of the two spherical covers are controlled by three specific motors. The solar tracker tracks the Sun by adjusting two aluminized folding mirrors and guides the sunlight into the FTS spectrometer. Thanks to its Camtracker mode, the solar tracker can track the Sun with a precision of ±0.1 mrad [1] (1 rad = 57.3). The meteorological station mounted near the dome includes sensors for air temperature (±0.30 °C), relative humidity (±3.0%), air pressure (±0.1 hPa), wind speed (±0.50 m·s–1 ), wind direction (±5.0°), solar radiation (±5.0% during daylight), and the presence of rain. These meteorological parameters are used to improve the spectral retrievals. The IFS125HR instrument is equipped with a multi-stage scroll pump to keep the interferometer under a vacuum, which enhances the stability of the instrument and minimizes the absorption interference of water vapor.

《Fig. 1》

Fig. 1. The main components of the high-resolution FTS observatory at Hefei, China. (a) An IFS125HR spectrometer; (b) a dome and a solar tracker [62,63].

We have operated the FTS observatory almost continuously since 2014. Before July 2015, we operated the FTS instrument to follow TCCON requirements and only collected near-infrared (NIR) solar spectra. We then modified the FTS instrument to collect NDACC observations in the middle-infrared (MIR) spectral range as well. Thereafter, NIR and MIR solar spectra were alternately collected in routine observations following the order shown in Table S3 in Appendix A. NIR and MIR spectra are collected with spectral ranges of 4000–11 000 cm–1 and 500–8500 cm–1 , temporal resolutions of 150 and 288 s, and spectral resolutions of 0.020 and 0.005 cm–1 , respectively. Both spectral resolutions can resolve the atmospheric absorption structures of all gases. Depending on the weather conditions, the number of both NIR and MIR spectra on each observation day vary from 11 to 103. An InGaAs detector and CaF2 beam splitter are used for NIR spectra acquisition, but no filters are needed. For MIR spectra acquisition, an InSb or cadmium telluride mercury (MCT) detector, KBr beam splitter, and seven filters with different transmission wavenumbers are used. NIR spectra are used to retrieve the total columns of CO2, CH4, N2O, HF, CO, H2O, HCl, and HDO [17,62], while MIR spectra are used to retrieve VMR profiles and the total columns of O3, CO2, NO2, N2, NO, HNO3, CO, HCl, HF, N2O, CH4, C2H6, HCN, ClONO2, H2O, HDO, ClO, OCS, H2CO, C2H2, C2H4, CF4, SF6, CH3D, COF2, CCl3F, CCl2F2, CHClF2, C3H8, HCOOH, CH3OH, PAN, and NH3 [18,66]. The major research relevance and significance of these gases are summarized in Table 1, including validation, global warming, ozone chemistry, and regional pollution.

《Table 1》

Table 1 Major research relevance and significance of the atmospheric constituents observed at the Hefei station.

《3. Methodology and characterization》

3. Methodology and characterization

《3.1. Retrieval and characterization of the NIR spectra suite》

3.1. Retrieval and characterization of the NIR spectra suite

We use the TCCON standard analysis software GGG to retrieve the total columns of atmospheric gases from the NIR solar spectra [17]. The GGG software was recently updated to version 2020, which has improved the a priori profiles for primary gases. Although GGG2020 will be used at the Hefei station in the near future, the NIR solar spectra that are reported in the present work were processed with the well-established GGG2014. The first step in GGG is to process the raw interferograms, correct them for phase errors and solar intensity variations, and convert them into spectra. The a priori profiles of meteorological parameters and trace gases are then generated, with the temperature, humidity, and pressure being based on analyses by the National Centers for Environmental Protection and the National Center for Atmospheric Research (NCEP/NCAR), and with the profiles of CO2, CO, CH4, and N2O being based on the Mark IV Balloon Interferometer (MkIV) balloon flights, satellite atmospheric chemistry experiment-FTS (ACEFTS), and GLOBALVIEW data. A series of influences such as the seasonal cycle, a secular increase, interhemispheric gradient, and stratospheric decay are considered. GGG then scales the a priori profiles iteratively to generate modeled spectra until an optimal fit to the measured spectra is achieved. Finally, the scaled profiles are integrated to compute the total column and are further converted into column-averaged dry air mole fractions (DMFs; referred to as XG for gas G), with the total column of oxygen (O2) being retrieved from the same spectra. This conversion process reduces the influence associated with surface pressure and H2O variations, and improves the precision. Setups for retrieving all gases with the NIR spectra suite are summarized in Table S4 in Appendix A. 

We follow the TCCON convention and use the formalism of Wunch et al. [22] to calculate the retrieval error. We conduct a series of sensitivity studies and perturb every source of error in the GGG forward model by a realistic range. We then calculate the fractional difference for each error source relative to the unperturbed case [17,62,63,67]. These sensitivities are calculated for spectra from a clear-sky day, over a large range of surface temperatures, solar zenith angles (SZAs), and water vapor amounts. In the sensitivity study, a priori profiles of and XHDO are modified by reducing the surface concentrations by 50%. Other gases are modified by shifting the profiles down by 1 km in altitude. The a priori profiles of temperature and pressure are increased by 1 K and 1 hPa at all altitudes, respectively. The continuum curvature is modified by fitting a third-order continuum curvature to all windows wider than 20 cm–1 . Zero level offsets, pointing offsets, surface pressure, field of view (FOV), and observer–Sun Doppler stretch (OSDS) are increased by 0.1%, 0.05°, 1 hPa, 7%, and 2 parts per million by volume (ppmv), respectively. The alignment of the instrument is perturbed by ±5% of the ideal condition. The total error is calculated in quadrature for each individual error. The results show that the largest contributors to the uncertainty budget for all gases are shear misalignments, OSDS, a priori profile shape errors, a priori temperature profile errors, continuum curvature, and zero-level offsets. For SZAs less than 80°, the total errors of  XHDO, and XHF are less than 0.25% (~1 parts per million (ppm)), 0.50% (~5 parts per billion (ppb)), 1.00% (~3 ppb), 4.00% (~10 ppb), 1.30%, 2.00%, and 8.00% (~4 ppb), respectively. 

《Fig. 2》

Fig. 2. TAKs of typical gases retrieved with the NIR spectra suite at the Hefei station AK: averaging kernel.

Total column averaging kernels (TAKs) describe the heightdependent sensitivity in the retrieval to errors in the profile shape. For an ideal column observation, the TAKs should be 1.0 at all heights; in reality, however, greater or lesser sensitivities to some heights in comparison with others also exist. The TAKs at all TCCON stations are similar and are a function of the fitted parameters and the GGG profile scaling algorithm. The TAKs are also slightly dependent on the SZA during the measurement time. Fig. 2 shows the TAKs for the Hefei station. Since the pressure broadening of the spectroscopic lines in the stratosphere is narrower than that in the troposphere, gases with saturated atmospheric absorption lines are insensitive to variations in the stratosphere. As a result, the stratospheric TAKs decrease as the air mass increases.

Time series of XHDO, and XHF retrieved with the NIR spectra suite at the Hefei station are shown in Fig. 3. As a formal TCCON station, all these datasets are archived in the TCCON database under open access. The main achievements with these datasets are summarized in Section 4.

↑  https://tccondata.org/.

《Fig. 3》

Fig. 3. Time series of the DMFs of typical gases retrieved with the NIR spectra suite at the Hefei station. Green and blue dots represent individual and daily mean results, while vertical error bars represent retrieval uncertainties ppt: parts per trillion.

《3.2. Retrieval and characterization of the MIR spectra suite》

3.2. Retrieval and characterization of the MIR spectra suite

We use the latest NDACC standard analysis software SFIT4 (version 0.9.4.4) with an optimal estimation (OE) technique to process the MIR spectra suite [16,68]. As listed in Table S3, the MIR spectra suite at the Hefei station is capable of retrieving the total columns and VMR profiles of more than 30 atmospheric constituents, including ten NDACC-mandatory species—namely, CH4, N2O, CO, HCN, C2H6, O3, ClONO2, HNO3, HF, and HCl—and many researchoriented species. Thus far, we have retrieved, analyzed, and investigated all the NDACC species and four research-oriented species: H2CO, NO2, CO2, and SF6. The retrieval settings for all NDACC gases follow the NDACC recommendations, as tabulated in Table S5 in Appendix A. For gases outside of the NDACC network, the retrieval settings are based on recommendations in Refs. [69–71]. The a priori profiles of all gases except water vapor (H2O) are from the statistical averages of the Whole Atmosphere Community Climate Model (WACCM) v6 simulations from 1980 to 2020 [8]. The a priori profiles of pressure, temperature, and H2O are extracted from the NCEP/ NCAR 6 h reanalysis data. Depending on the gas, the spectroscopic parameters are provided by the HITRAN 2008 or HITRAN 2012 databases [8,10,67,71]. We deal with the absorption interference of H2O differently for different gases. For all gases except ClONO2 and HNO3, the H2O profile is retrieved together with the profile of the target gas. For ClONO2 and HNO3, we treat H2O as another interfering species—that is, we only scale the a priori profile of H2O. Since the abundance and interference are gas dependent, different deweighting schemes are used for different gases. We use a deweighting signal-to-noise ratio (SNR) of 500 for CO and 300 for HCl retrievals, but no de-weighting SNRs are used for other gases [72,73]. The empirical usage of these de-weighting SNRs can reduce the retrieval residual and improve the accuracy of the spectral retrieval for each gas. The diagonal values of the a priori covariance matrices Sa are set to the standard deviation of the WACCM v6 simulations. The abundances of some gases (e.g., CO2, CH4, and N2O) have changed significantly since 1980, so the standard deviation of the WACCM simulations may overestimate the natural variabilities of these gases. In order to minimize the influence, Gaussian or exponential inter-layer correlation shapes are used for these gases to constrain the shapes of the retrieved profiles. The diagonal values of the noise covariance matrices are set to the inverse square of the spectral SNR. The non-diagonal values of Sa and are set to zero. The instrument line shape (ILS) available from optical alignment diagnosis is adopted in the retrieval [72,73].

For each retrieval, the averaging kernels (AKs) describe the height-dependent sensitivity. The area of the AKs at a selected height is calculated by summing up all the elements of the corresponding AKs [74]. This represents the fraction of the retrieval at that height, which comes from the measurement [68]. A value close to zero indicates that the retrieval is nearly independent of the measurement but is determined by the a priori information. The degrees of freedom for the signal (DOFS) are calculated as the trace of the AKs matrix and quantify the number of independent partial columns in the retrieval. The AKs and their areas for some MIR gases at the Hefei station are shown in Fig. 4. Highsensitivity ranges (i.e., a height range with a sensitivity of more than 0.5) for HCN, H2CO, SF6, CO, and C2H6 are in the troposphere; for HF, NO2, and HCl, these are in the stratosphere, while for O3 and CH4, they are in both the stratosphere and troposphere. The typical DOFS over the total atmosphere for CH4, NO2, CO, HCN, C2H6, O3, H2CO, SF6, HF, and HCl are 3.5, 1.0, 3.0, 1.1, 1.5, 4.8, 1.1, 1.0, 1.3, and 1.5, respectively. The potential independent altitude information to be retrieved from the measurements for CH4, O3, and CO at Hefei can be separated into three, four, and three independent partial layers, respectively. However, we can roughly obtain one degree of freedom on the profiles of HCl, NO2, H2CO, HF, SF6, C2H6, and HCN. 

《Fig. 4》

Fig. 4. AKs (colored fine lines) of some tropospheric (C2H6, CH4, CO, H2CO, HCN, and SF6) and stratospheric (HCl, HF, NO2, and O3) gases retrieved with the MIR spectra suite at the Hefei station. The area of AKs for each gas is plotted in a black dashed line and scaled by a factor of 0.1. Demonstrations for all gases are prescribed from randomly selected retrievals.

The error budget for MIR spectra retrievals at Hefei is calculated using the formalism of Rodgers [68], and all error components are separated into random and systematic errors, as listed in Table S5. The covariance matrix of each interfering gas is prescribed from the WACCM v6 climatology. The covariance matrix of temperature is set to 3–7 K in the stratosphere and 2–5 K in the troposphere, based on the differences between the sonde and the NCEP temperature profiles. The covariance matrix of the measurement error is set to the inverse square of the spectral SNR. We regularly diagnose the optical alignment of the instrument and realign the interferometer when indicated. The FTS instrument at Hefei is assumed to be close to the ideal state, and the uncertainties for the FOV, OPD, interferogram phase, zero level, and background curvature are all assumed to be 1%. The covariance matrix for the retrieval parameter and smoothing error are prescribed from the SFIT4 retrieval outputs. Depending on the gas, we use 5% or 10% for uncertainties of line intensity, air-broadening, and self-broadening coefficients. Total retrieval error is calculated in quadrature for all systematic and random errors. Typical retrieval errors over the total atmosphere at the Hefei station for CH4, NO2, CO, HCN, C2H6, O3, SF6, H2CO, HF, and HCl are estimated to be 4.5%, 12.4%, 4.7%, 11.3%, 6.2%, 5.7%, 5.8%, 12.3%, 5.1%, and 5.2%, respectively.

The total column time series of some typical gases retrieved with the MIR spectra suite at the Hefei station are shown in Fig. 5. Pronounced day-to-day variabilities and seasonality of these gases are observed. The main achievements that were made with the results retrieved from the MIR spectra suite are summarized in Section 4.

《Fig. 5》

Fig. 5. Total column time series of some typical gases retrieved with the MIR spectra suite at the Hefei station. The seasonality of each gas (represented by blue dots) is fitted using a bootstrap resampling model with a third Fourier series. Vertical error bars represent retrieval uncertainties. The units of the y-coordinate are moleculars·m–2 .

《4. Main achievements obtained at the Hefei station》

4. Main achievements obtained at the Hefei station

The Hefei station has archived almost seven years of high-quality data for many atmospheric constituents. We have conducted a series of studies with these datasets to understand high-resolution spectral retrieval and harmonization, investigate the overall characteristics of key atmospheric constituents, perform emission estimates and satellite and CTM evaluations, and determine seasonal levels of air pollution, along with its sources, transport, and chemical mechanisms. Table 2 [8,10,13,15,62,63,66,67,71,73,75–81] provides a summarized inventory of these studies based on FTS observations at the Hefei station during 2014–2020, categorized by key findings, data period, theme (i.e., the main gases considered), waveband, atmospheric domain, and topic. The main achievements at the Hefei station are as follows.

《Table 2 》

Table 2 Inventory of papers based on FTS observations at the Hefei station during 2014–2020, categorized by key findings, data period, theme (i.e., main gases considered), waveband, atmospheric domain, and topic.

《4.1. ILS characterization and retrieval harmonization 》

4.1. ILS characterization and retrieval harmonization 

Most FTS stations within the TCCON/NDACC networks utilize different sizes of field stop or attenuators to adapt the received optical intensity of the detector. These behaviors may modify the ILS of a high-resolution FTS spectrometer, which can cause retrieval biases if the resulting ILS drifts are not properly addressed. Sun et al. [75] first quantified the ILS drifts of a high-resolution FTS spectrometer with respect to different kinds of attenuators. They inserted five attenuators in front of or behind the interferometer and conducted ILS measurements to monitor the ILS of the high-resolution FTS instrument. Sun et al. [75] concluded that any of the five attenuators caused negligible ILS drifts and could be used to adapt the received optical intensity of the detector. However, the usage of different sizes of field stop may cause dramatic ILS drifts due to inconsistent mechanical errors among the selected field stops.

Sun et al. [66,73] further simulated different kinds of ILS drifts and quantified their effect on the retrieval of all NDACC gases. The sensitivities of the total column and partial column to different levels of ILS drift were investigated. Sun et al. found that the tropospheric gases are less sensitive to ILS drift than the stratospheric gases; furthermore, for the same level of drift, the negative ILS influence is smaller than the positive ILS. For an ILS drift of 10%, the total columns of CH4, N2O, CO, HCN, C2H6, O3, ClONO2, HNO3 , HF, and HCl changed by 0.04%, 0.20%, 2.10%, 0.75%, 1.10%, 1.90%, 23.00%, 0.70%, 3.00%, and 4.00%, respectively [73]. To limit the difference in the total column of ClONO2 within 10% and those of the other gases within 1%, the upper limits of positive ILS drifts for O3 , ClONO2, HNO3, HF, HCl, CO, HCN, and C2H6 are 6%, 5%, 15%, 5%, 5%, 5%, 13%, and 9%, respectively; the upper limits of negative ILS drifts for O3, HF, and HCl are 6%, 12%, and 12%, respectively; and the influence of ILS drift for CH4 and N2O is negligible [73] .

Sun et al. [66] concluded that the effect of ILS drift on a partial column is gas and altitude dependent. The influences of ILS drift on CH4 and N2O are smaller than those on O3 and CO. In order to suppress the influence on a partial column within 2%, it is suggested that a precise ILS be used for CO and O3. For CH4, the upper limits of the positive ILS drift for altitude ranges of 0–7 and 16–37 km are 13.0% and 1.5%, respectively. The upper limits of the negative ILS drift for altitude ranges of 7–16 and 16–37 km are 12.0% and 4.5%, respectively. The influence of the negative ILS drift on the partial column of CH4 at 0–7 km and that of the positive ILS drift at 7–16 km are negligible [66]; for N2O, the upper limits of the positive ILS drift for altitude ranges of 0–5.0, 5.0–11.5, 11.5–20.0, and 20.0–35.0 km are 6.0%, 4.0%, 1.5%, and 1.5%, respectively. The upper limits of the negative ILS drift for altitude ranges of 11.5–20.0 and 20.0–35.0 km are 2.0% and 1.5%, respectively. However, the influence of the negative ILS drift on partial columns of N2O at 0–5.0 and 5.0–11.5 km are negligible [66] .

The NIR and MIR spectra at the Hefei station were collected with different optical filters, beam splitters, and detectors and were retrieved with different Micro-windows (MWs), spectral modeling schemes, and iterative methods, which could result in retrieval biases. In order to determine whether the NIR and MIR HF datasets at the Hefei station can be combined into a single reliable dataset, Yin et al. [67] compared these two datasets and identified their discrepancies. They found that the NIR and MIR XHF data exhibit consistent seasonality, but the daily mean values of NIR XHF are generally lower than the MIR data, with an average difference of (6.90 ± 1.07) parts per trillion by volume (pptv) [67]. As a result, Yin et al. [67] added a bias of 6.90 pptv to the NIR dataset before combining these two datasets for investigation.

Inconsistent field stops and detectors were used for the NIR spectra in the early stage of operation. These NIR spectra also lacked synchronous meteorological data and were subject to timTable 2 ing errors due to an incomplete observation protocol. All these observation deficiencies would cause biases between individual measurements if not properly treated. In order to obtain a consistent long-term  time series, Tian et al. [63] conducted a series of corrections to eliminate these potential biases and applied a bias of 11.0 ppb to the dataset derived before July 2015 [63].

《4.2. Overall characteristics of key atmospheric constituents》

4.2. Overall characteristics of key atmospheric constituents

We have conducted several works using the continuous FTS observations at the Hefei station to investigate the overall characteristics of the key atmospheric constituents over east China. Wang et al. [62] first investigated the seasonal cycles of CO2 and CO with two year NIR measurements at the Hefei station. Wang et al. [62] found that the FTS time series of showed an obvious seasonal cycle, while the time series of displayed no clear seasonality. The researchers further analyzed the relationship between and , and found that these two species were only correlated for individual days and varied according to prevailing wind conditions. The correlation slopes of these two species in Hefei are usually smaller than those in Beijing, indicating that the two megacities have different emission characteristics. Tian et al. [63] reported an increasing change rate of (0.56% ± 0.15%) a–1 in based on three years’ worth of NIR observations at the Hefei station (2014–2017). Diurnal generally exhibited an increasing change rate in summer (June–July–August; JJA) and showed a decreasing change rate or remained constant in other seasons. The time series show seasonal maxima in JJA and minima in winter (December–January–February; DJF). The annual means in January were 30–55 ppb lower than those in August.

Sun et al. [8,10,13,15] investigated the seasonality and interannual trends of O3, HCN, H2CO, and C2H6 over east China using different time scales of MIR observations at the Hefei station. Tropospheric O3 (0–12 km) was shown to be characterized by a seasonal maximum in June and minimum in December [8]. Tropospheric columns of O3 in June were 50% higher than those in December. The day-to-day variabilities of tropospheric O3 in spring (March–April–May; MAM) and JJA were generally larger than those in autumn (September–October–November; SON) and DJF, and a broad high-value range within MAM/JJA was observed. These findings differ from those of Vigouroux et al. [82], who showed that tropospheric columns of O3 at relatively clean sites usually reach a seasonal maximum in springtime.

The tropospheric columns (0–15 km) of HCN at the Hefei station showed pronounced seasonal variations, with three seasonal peaks in May, September, and December, respectively [10]. The tropospheric column of HCN was characterized by a seasonal maximum in May and minimum in November. In general, the tropospheric columns of HCN at the Hefei station are higher than the FTS observations at the NDACC sites arrival heights (78°S), Lauder (45°S), Reunion Maido (21°S), Mauna Loa (20°N), Izana (28°N), Rikubetsu (43°N), Toronto (44°N), Jungfraujoch (47°N), Bremen (53°N), Kiruna (68°N), and Ny Alesund (79°N). The tropospheric columns of HCN between September 2015 and July 2016 were 5%–46% larger than those during the same period in other years.

The time series at the Hefei station reached a seasonal maximum in July and a minimum in January [13]. The FTS values in January were 42.9% lower than those in July. The variabilities of at the Hefei station are in good agreement with those of surface H2CO measurements; thus, the FTS column observations were applied as a representation of near-surface conditions by Sun et al. [13] for source separation. The time series at the Hefei station showed a positive change rate of (2.37% ± 0.70%) a–1 from 2015 to 2019. In the studied years, Sun et al. [13] found that hydroxyl (OH) production from H2CO photolysis was comparable with that from O3 photolysis, indicating that H2CO photolysis plays an important role in OH generation over east China.

Sun et al. [15] found that the tropospheric (0–12 km) columnaveraged dry air mole fraction (troDMF) of C2H6 reached a seasonal maximum in December and a minimum in July, and showed a negative change rate of (–2.60% ± 1.33%) a–1 from 2015 to 2020. The C2H6 troDMF seasonality is characterized by a winter maximum and a summer minimum because the tropospheric OH oxidation capability in summer is higher than that in winter [15]. Both emission and meteorological factors have negative effects on C2H6 troDMF in JJA/SON and positive effects in DJF/MAM. Regional and long-range transport of C2H6 over east China is largely driven by the Asian monsoon and mid-latitude westerlies [15].

In addition to tropospheric gases, the seasonal and interannual variabilities of a few stratospheric gases have been investigated. Yin et al. [71] found that the stratospheric NO2 columns at the Hefei station incrementally increase in the first half of the year and decrease in the second half of the year, with a seasonal maximum in June and a minimum in January. The stratospheric columns of NO2 in June are on average 39.20% ± 8.95% higher than those in January. High values of summertime stratospheric NO2 columns are attributed to a faster photolysis of N2O5 and HNO3 by higher solar intensity in summer [71]. Further studies by Yin et al. [67,78] and Shan et al. [81] concluded that the seasonal cycles of HCl, HF, and HNO3 in the stratosphere are associated with the variability of tropopause height. Stratospheric HCl, HF, and HNO3 columns showed a seasonal maximum in January and a minimum in September [67,78,81]. Stratospheric NO2, HCl, HF, and HNO3 columns showed decreasing change rates of (–0.34% ± 0.05%) a–1 between 2015 and 2018, (–1.83% ± 0.13%) a–1 between 2015 and 2019, (–0.38% ± 0.22%) a–1 between 2015 and 2020, and (–9.45% ± 1.20%) a–1 between 2017 and 2019, respectively.

《4.3. Emission estimates》

4.3. Emission estimates

For regions where CO and CO2 share common sources, CO emission can be inferred by studying the correlation between CO and CO2. Shan et al. [76] utilized the enhancement ratios derived from the NIR measurements at the Hefei station to calculate the CO emission over the megacity Hefei. They also investigated the discrepancy between the CO emissions inferred from emission inventories and those from the measurements. In order to do so, the researchers compared the enhancement ratios at the Hefei station derived from FTS, satellite, in situ measurement, and emission inventory from September 2015 to August 2017. The Emission Database for Global Atmospheric Research (EDGAR) emission inventories and the Peking University (PKU) emission inventories were used in the study. Shan et al. [76] concluded that the two emission inventory-based ratios are substantially larger than those based on measurements from ground-based FTS, in situ instruments, and satellite. Furthermore, the error bars of the emission inventory-based ratios are higher than those based on measurements, indicating that the ratios derived from the two emission inventories have high uncertainties [76].

Shan et al. [76] further used the FTS-based ratio and the CO2 emissions derived from the EDGAR and PKU inventories to infer CO emissions. The CO emissions calculated with the EDGAR inventory-based CO2 emissions and the ground-based FTS ratios were about (11.27 ± 0.91) Tg·a–1 for 2015–2016 and (12.35 ± 0.74) Tg·a–1 for 2016–2017. The CO emissions during the same periods estimated from the PKU inventory-based CO2 emissions and the ground-based FTS ratios were about (10.96 ± 0.88) and (11.95 ± 0.71) Tg·a–1 , respectively. As the PKU and EDGAR inventories can well reproduce actual CO2 emissions, the CO emissions calculated with the two inventory-based CO2 emissions and the FTS ratios are in good agreement. However, Shan et al. [76] observed a large discrepancy between the CO emissions derived from the emission inventories and those derived from ground-based FTS observations. The CO emissions estimated from the EDGAR inventory are 19.12 and 17.88 Tg·a–1 , while those from the PKU inventory are 38.07 and 35.60 Tg·a–1 for 2015–2016 and 2016–2017, respectively. This finding indicates that both emission inventories substantially overestimated the real CO emissions near Hefei.

《4.4. Evaluating satellite observations》

4.4. Evaluating satellite observations

The continuous FTS observations at the Hefei station provide a series of valuable datasets that can be used to validate satellite observations. To compare the measurements obtained from satellites and ground-based FTS, their differences in vertical resolution and the a priori profile are considered. First, the satellite data are re-gridded to the FTS a priori profile to minimize the a priori profile-associated difference [56]. The FTS data are then smoothed by the concurrent satellite AKs to minimize the vertical resolutionassociated difference. Finally, the re-gridded and smoothed datasets of these two observations are compared.

Wang et al. [62] found that the ground-based FTS time series at the Hefei station are in good agreement with both the greenhouse gases observing satellite (GOSAT) and Orbiting Carbon Observatory-2 (OCO-2) observations [62]. Compared with the ground-based FTS observations, the daily median of the GOSAT and OCO-2 time series from 2015 to 2017 have biases of (–0.52 ± 1.63) and (0.81 ± 1.73) ppm, respectively. The GOSAT and OCO-2 observations can reproduce the daily median of FTS time series with a correlation coefficient (r) of 0.79 and 0.83, respectively. Tian et al. [63] observed that the differences between the GOSAT and FTS data at the Hefei station (GOSAT–FTS) were within the error budgets of both datasets and were season independent, with a mean value of (1.6 ± 13.0) ppb (0.09% ± 0.70%). This mean difference is comparable to similar comparisons at other TCCON sites [50,83]. The seasonal cycles of observed by GOSAT and FTS are in good agreement with a correlation coefficient (r) of 0.77.

Yin et al. [71] found that the stratospheric NO2 time series from 2015 to 2018 observed by the ozone monitoring instrument (OMI) satellite are in good agreement with the concurrent ground-based FTS data at the Hefei station. The seasonal cycles of stratospheric NO2 observed by OMI and FTS are correlated with a correlation coefficient (r) of 0.841 [71]. The annual mean difference between OMI and FTS (OMI – FTS) is (1.48 ± 5.33) × 1014 molecules∙cm–2 (4.82% ± 17.37%) [71], which is within the uncertainties of both datasets. The OMI versus FTS differences in MAM are higher than those in the other seasons. Yin et al. [71] also found that the negative stratospheric NO2 trend between 2015 and 2018 obtained by OMI observations (–0.91% ± 0.09%) a–1 is consistent with that obtained by ground-based FTS. The FTS measurements at Hefei station have also been used by Oshio et al. [79] as one of the reference datasets for the bias correction of the CH4/CO2 total column ratio retrieved from GOSAT observations and by Hedelius et al. [77] for the evaluation of measurements of pollution in the troposphere (MOPITT) XCO v7 retrievals.

The stratospheric HNO3 and HCl columns at the Hefei station were compared with Microwave Limb Sounder (MLS) satellite data. Shan et al. [81] found that MLS satellite data showed similar seasonal variations and annual rates as the Fourier-transform infrared (FTIR) data, and that the stratospheric HNO3 and HCl columns of the two datasets have correlation coefficients (r) of 0.87 and 0.88, respectively. The mean bias between the satellite and FTIR data of the stratospheric HNO3 and HCl columns is –8.58% ± 12.22% and 4.58% ± 13.09%, respectively.

The results from Wang et al. [62], Tian et al. [63], Yin et al. [71], Hedelius et al. [77], Oshio et al. [79], and Shan et al. [81] have verified the ability of the ground-based FTS observations at the Hefei station to detect diurnal, seasonal, and interannual variabilities of key atmospheric constituents in both the troposphere and stratosphere and to validate satellite observations.

《4.5. Evaluating CTM simulations》

4.5. Evaluating CTM simulations

Several studies have been conducted with the continuous FTS observations at the Hefei station to evaluate the performance of CTM in the polluted region of east China. Since the vertical resolution of CTM differs from the ground-based FTS observations, a smoothing correction is applied to the model profiles [84]. First, the model profiles are re-gridded to the FTS height grid to reconcile a common height grid. The re-gridded profiles are then smoothed by a priori profiles and the concurrent AKs of FTS [68,84]. Subsequently, the smoothed model profiles are converted to a total column or DMF using the concurrent re-gridded air density profiles from the CTM. Finally, the converted CTM data are compared with the ground-based FTS data.

The FTS time series of many gases, such as CH4, NO2, H2CO, C2H6, and CO2, obtained at the Hefei station have been used to validate the GEOS-Chem model. In general, GEOS-Chem simulations can reproduce the seasonal cycles and interannual trends of the gases observed by ground-based FTS. The correlation coefficients (r) between GEOS-Chem and the FTS for , the stratospheric column of NO2, , the C2H6 troDMF, and the CO2 troDMF are 0.860, 0.8580, 0.780, 0.880, and 0.890, respectively [13,15,63,71,80]. However, the observed day-to-day variabilities of these gases cannot always be reproduced by the GEOS-Chem model [13]. For example, the GEOS-Chem versus FTS differences for in JJA are 35–55 ppb larger than those in DJF [63]. For the stratospheric column of NO2, the differences in SON are smaller than those in other seasons [71]. For the C2H6 troDMF, the observed monthly mean values are overestimated by 35.6% in July and underestimated by 17.4% in December by the GEOS-Chem model [15]. These discrepancies can mostly be attributed to uncertainties in the vertical mixing and horizontal transport schemes simulated by the GEOS-Chem model at a coarse spatial resolution, which make it difficult to match column observation over a single point. In particular, the Hefei station is located in a highly industrialized and densely populated region in east China. The regional differences in the atmospheric burdens of these gases could affect the comparison with the FTS observations. In addition, the differences between simulation and measurement could be associated with the uncertainties in emission inventories and chemical mechanisms. For example, many oxidation pathways of volatile organic compounds (VOCs) may not be optimally implemented or may not even be considered in the GEOS-Chem model. Furthermore, the geographical distributions of various emission sources and how these sources affect the air masses over the polluted Hefei station may not always be properly addressed in the model. 

Nevertheless, the GEOS-Chem versus FTS differences for , the stratospheric column of NO2, , the C2H6 troDMF, and the CO2 troDMF are (–8.09 ± 17.83) ppb (0.45% ± 0.95%), (2.36 ± 2.33) × 1014 molecules∙cm–2 (7.66% ± 7.49%), (–0.05 ± 0.20) ppbv (–2.60% ± 10.40%), (–0.02 ± 0.05) ppbv (–1.60% ± 4.20%), and (–1.68 ± 1.23) ppbv (–0.42% ± 0.31%), respectively, which are within the FTS uncertainty budget of the respective gases. These results verify that GEOS-Chem can simulate the burdens, seasonal cycles, and interannual trends of CH4, NO2, H2CO, C2H6, and CO2 over the polluted east China [13,15,63,71,80]. 

《4.6. Air pollution, sources, and transport》

4.6. Air pollution, sources, and transport

With tropospheric columns of O3, H2CO, and CO provided by FTS, the tropospheric column of NO2 provided by OMI, and the air mass trajectories, Sun et al. [8] investigated the O3 photochemical regime at the Hefei station. Sun et al. concluded that the emissions in the Shandong and Hebei Provinces in north China, the Jiangsu and Anhui Provinces in east China, the Shaanxi, Hubei, and Hunan Provinces in central China, and the Henan and Shanxi Provinces in northwest China contribute to the O3 burdens at the Hefei station [8]. The level and variability of tropospheric O3 in MAM/JJA are higher than those in SON/DJF because air masses in MAM/JJA have a stronger tendency to originate from highly industrialized and densely populated regions. Sun et al. [8] further used the H2CO/NO2 tropospheric column ratio as a proxy to explore the photochemical production sensitivity of O3 (PO3), and found that the temporal PO3 during the measurement time is mainly limited to NOx in MAM/JJA and is limited to VOCs or mixed VOCs–NOx in SON/DJF. Emissions control for VOC and NOx could effectively mitigate O3 pollution in the SON/DJF and MAM/JJA seasons, respectively [8].

The GEOS-Chem model captured the observed variabilities of CO, H2CO, and C2H6 at the Hefei station and was thus used for source attributions of these gases by Sun et al. [10,13,15]. Sun et al. [10] found that the enhancements of tropospheric HCN () and CO () are correlated, which indicates that the CO and HCN enhancements in the troposphere can be attributed to the same sources [10]. Taking this advantage, Sun et al. [10] used GEOSChem-tagged CO simulation to determine the source of HCN and concluded that the seasonal maximum of HCN in May can be mainly attributed to BB emissions in SEAS (41.0% ± 13.1%), Africa (AF) (22.0% ± 4.7%), and Europe and boreal Asia (EUBA) (21.0% ± 9.2%). The seasonal maximum of HCN in September can be mainly attributed to BB emissions in EUBA (38.0% ± 11.3%), SEAS (14.0% ± 3.3%), AF (26.0% ± 6.7%), and North America (NA) (13.7% ± 8.3%). With respect to the seasonal maximum of HCN in December, dominant influences are from AF (36.0% ± 7.0%), NA (18.7% ± 5.2%), and EUBA (21.0% ± 5.2%). The enhanced tropospheric columns of HCN between 2015 and 2016 were due to elevated BB emissions in SEAS, Oceania (OCE), and EUBA during that period. More specifically, the elevated BB emissions in OCE dominated the HCN enhancements in the second half of 2015, while elevated BB emissions in SEAS dominated the HCN enhancements in the first half of 2016 [10].

Sun et al. [13] used the GEOS-Chem model to quantify the contributions of different geographical regions and emission categories to the summertime H2CO enhancements at the Hefei station [13]. The results showed that the oxidation of nonmethane volatile organic compound (NMVOCs) and CH4 caused the summertime H2CO enhancements over Hefei, with contributions of 56.73% and 43.27%, respectively. The NMVOC-related summertime H2CO enhancements were largely caused by NMVOC emissions in east China. Due to the increase in photochemical H2CO resulting from increases in both NMVOCs and CH4, H2CO emissions from 2015 to 2019 showed an increasing trend [13]. With similar sensitivity simulations as those used in their earlier work [13], Sun et al. [15] concluded that anthropogenic (biofuel plus fossil fuel) emissions accounted for 49.2% and natural (BB plus biogenic) emissions accounted for 37.0% of the C2H6 burden over east China. The observed C2H6 variability at the Hefei station was mainly driven by the variation in C2H6 emissions within China (74.0%), where north, east, and central China accounted for most of the contribution (57.6%). A decrease in C2H6 from 2015 to 2020 points to an improvement in the air quality in China, which can be attributed to a reduction in transported and local C2H6 emissions in recent years [15].

《5. Conclusion and future outlook》

5. Conclusion and future outlook

This study summarized the major achievements obtained from the long-term observations of atmospheric constituents at the first ground-based high-resolution FTS station in China, which has been located in Hefei, east China, since 2014. These results have improved our current knowledge on the retrieval and characterization of high-resolution spectra, climate change, emission estimates, satellite and CTM validations, and the sources and transport of atmospheric pollution. New findings on physical processes and chemical mechanisms will be implemented into CTMs in the near future. In turn, the validated CTMs and the measurement data will be used to further understand the impact of atmospheric production and depletion processes on a larger scale. Based on these achievements and the overall envisaged concept of this station, the following scientific studies are expected to be conducted at the Hefei station in future.

(1) In order to use ground-based total column or profile retrievals to investigate different scales of atmospheric research and provide a standard transfer among different observation platforms, the ground-based FTS dataset must be converted to the same validation scale as the World Meteorological Organization (WMO) scale. We have not yet validated the ground-based FTS data at Hefei to in situ profile measurements based on airborne instruments. Such a comparison will enhance the accuracy level of the ground-based FTS observations. At present, we simply assume that the calibration value that is valid for other TCCON sites is also valid for Hefei, and scale our NIR observations by the observed biases reported by Wunch et al. [17]. Therefore, we look forward to using airborne profile measurements to further improve the accuracy of the observations.

(2) To understand the interaction between atmospheric constituents and climate change, long-term measurements of a sufficient number of atmospheric constituents are needed. Therefore, trend-quality stable measurements of all atmospheric constituents at the Hefei station will be continued. A few more species will be investigated, such as PAN, OCS, NH3, and chloro–fluoro–carbons (CFCs); this list of species will continue to increase over time, in keeping with the evolution of the retrieval techniques and new spectral data. Unlike other spectrometers within the TCCON/ NDACC networks, the IFS125HR used at the Hefei station covers the UV spectral range; however, we have not attempted to collect UV spectra in routine observations. We plan to collect UV spectra in future, which will further increase the number of atmospheric constituents. In order to contribute to the mitigation of climate change, China has explicitly proposed reaching a peak in its CO2 emissions by 2030 and realizing carbon neutrality by 2060. The Hefei station will provide scientific assistance for the Chinese Government to develop green economy policies and to achieve carbon neutrality and the goals of the Paris Agreement.

(3) As a part of the AEOS infrastructure, the high accuracy of the FTS dataset will allow it to be used extensively to validate other instruments within the AEOS framework and above. Despite their outstanding capabilities in precision and stability, IFS125HR spectrometers have limitations. The IFS125HR is an expensive and ponderous spectrometer, and its operation requires a significant amount of infrastructure. Furthermore, it is difficult and time consuming to regularly maintain its optical alignment. To overcome these problems, the usage of smaller, cheaper, andmore easily transportable spectrometers such as the EM27/SUN and the VERTEX-80/ SUN was recently investigated. Both the EM27/SUN and the VERTEX-80/SUN have been verified to have similar features as the IFS125HR for certain gases (e.g., CO2, CO, and H2CO) but have better mobility [70,85]. The high accuracy of the FTS dataset will be used to validate the flexible instruments EM27/SUN and VERTEX-80/SUN before they are used to calculate the emission flux of a city or industrial facility. The validated EM27/SUN and VERTEX-80/SUN observations can then complement the existing high-resolution FTS observation dataset in China. The FTS dataset will be also used to validate the envisaged Chinese environmental satellites in future.

(4) The Hefei station has archived almost seven years of highquality data for many atmospheric constituents. These data will be used together with the measurements from other platforms or CTM simulations and will continue to be available to the community in order to help us better understand climate change, air pollution, pollution sources and transport, and their impacts on crops and human health in China and around the world.

(5) The achievements made at the Hefei station have benefited from the close cooperation between AIOFM-CAS and the TCCON/ NDACC networks. It is planned for the Hefei station to be more involved in international TCCON/NDACC observations and scientific activities, which will contribute to an understanding of the carbon cycle and climate change on a global scale.

(6) New spectroscopic properties will be studied with the high spectral resolution and precision of the ground-based FTS datasets at the Hefei station. These analyses of the spectral lines will improve the retrieval precision and accuracy of all kinds of observations at the Hefei station. We will collaborate with the spectroscopic community to improve the interpretation of atmospheric absorption and extinction.

Finally, the achievements based on our work at the Hefei station exemplify how a new station can improve the current understanding of climate change, air pollution, and pollution transport both on a regional scale and globally. Considering that humankind is facing the unprecedented challenge of climate change, it is necessary to build more stations of this kind in different regions around the globe—especially in Africa, South America, and Asia, where only sparse ground-based observations are available. It is only with sufficient measurements and close cooperation between all sites that we can better understand the atmosphere and predict the ‘‘future climate” using climate models.

《Acknowledgments》

Acknowledgments

This work is supported by the Youth Innovation Promotion Association, Chinese Academy of Sciences (CAS) (2019434) and the Sino–German Mobility programme (M-0036). We thank the NDACC/TCCON networks for providing the state-of-the-art software. We thank David Griffith and Nicholas Jones from university of Wollongong, James Hannigan from National Center for Atmospheric Research (NCAR), Thomas Blumenstock and Frank Hase from Karlsruhe Institute of Technology (KIT), Corinne Vigouroux from Royal Belgian Institute for Space Aeronomy, Emmanuel Mahieu from University of Liège, and many other colleagues within NDACC/TCCON networks for their help in setting up/operating the FTS spectrometer or analyzing the spectra at the Hefei station.

《Compliance with ethics guidelines》

Compliance with ethics guidelines

Cheng Liu, Youwen Sun, Changgong Shan, Wei Wang, Justus Notholt, Mathias Palm, Hao Yin, Yuan Tian, Jixi Gao, and Huiqin Mao declare that they have no conflict of interest or financial conflicts to disclose.

《Authors’ contribution》

Authors’ contribution

Cheng Liu and Youwen Sun designed and wrote the paper with inputs from all coauthors. Changgong Shan is in charge of TCCON retrieval, analysis, and plotting. The rest of the authors contributed to this work by providing constructive comments.

《Appendix A. Supplementary data》

Appendix A. Supplementary data

Supplementary data to this article can be found online at https://doi.org/10.1016/j.eng.2021.11.022.

《Data availability》

Data availability

Dataset retrieved from the NIR spectra suite is archived in the TCCON database https://tccondata.org/. Dataset retrieved from the MIR spectra suite can be available from https://doi.org/10. 5281/zenodo.5205257. Data format and usage of the NIR and MIR datasets at the Hefei station follow the TCCON and NDACC conventions, respectively.