中国首个地基高分辨率傅里叶变换光谱观测站的大气成分长期观测研究

摘要

长期观测主要大气成分的体积混合比（VMR）廓线和柱总量对掌握中国气候变化和碳收支具有重要意义。本文对我国首个全球碳总柱观测网（TCCON）高分辨率傅里叶变换光谱仪（FTS）地基观测站，即合肥观测站开展了系统性研究。该观测站可观测30多种大气成分的柱总量和VMR廓线。本文公布了2014年以来合肥观测站观测到的部分关键大气成分的时间序列，总结了迄今为止该站点取得的主要研究成果，包括光谱反演表征和归一化、关键大气成分的总体演变特征、排放估计、卫星和化学传输模式（CTM）校验以及对大气污染来源和传输的相关研究。同时，本文还对合肥观测站的观测、科学研究和未来研究计划进行了展望。中国明确提出2030年二氧化碳排放达到峰值，2060年实现碳中和。合肥观测站将为中国政府制定绿色经济政策、实现碳中和及《巴黎协定》目标提供科学支撑。

正文

《1、引言》

1、引言

近十年来，空气污染对我国社会经济发展和人民健康带来了巨大挑战[1‒2]。特别是中国东部的特大城市群，由于人类活动密集，空气污染最严重[3‒9]。近年来，由于采取了一系列减排措施，中国各地的空气污染已显著减少[10‒11]。然而，在工业化程度高、人口密集的中国东部地区，大气污染仍然较严重[12‒13]。该地区大气污染的严重程度、范围和复杂性仍远超世界其他地区[11,13‒15]。长期观测主要大气成分的体积混合比（VMR）廓线和柱总量，对掌握中国空气污染和碳排放的来源、化学机制和输送过程，以及潜在的调控措施都具有重要意义。此外，精确探测关键对流层物种[如CO、CO₂、CH₄、N₂O、C₂H₆、SF₆、碳酰硫（OCS）、H₂CO、HCN和过氧乙酰硝酸盐（PAN）]和平流层物种（如O₃、HNO₃、HF、HCl、ClONO₂）的变化趋势，并深入探索它们对气候变化的影响，对于建立气候变化与大气成分之间的科学联系，验证卫星、机载或其他地面观测平台，支撑聚焦大气传输过程的外场观测计划，改进理论化学传输模式（CTM）都具有重要意义[16‒19]。

地基高分辨率直接太阳光傅里叶变换光谱仪（FTS）已成为探测大气成分VMR廓线和柱总量的有力工具[16‒17,20‒23]。全球碳总柱观测网（TCCON）和大气成分变化探测网（NDACC）均利用高分辨率FTS观测大气成分变化。TCCON和NDACC分别成立于2004年和1992年，可以探测30多种大气成分的柱总量和VMR廓线。这些探测结果已广泛应用于碳循环[24‒30]、大气污染与输送[8,10,13,31‒34]、平流层O₃[18]研究，光谱线模型开发[35‒45]，卫星观测验证和卫星反演算法开发[46‒56]，CTM评估等[57‒61]。然而，大多数TCCON/NDACC站都位于北美和欧洲，世界上其他地区的观测站仍然非常稀少。目前，中国只有合肥站（117.2°E, 32.0°N，海拔高度30.0 m）和香河站（116.96°E, 39.75°N）两个合格的TCCON运行站点[8,15]。

本研究综述了中国第一个地基高分辨率TCCON观测站——合肥观测站的系统性研究成果。合肥观测站靠近中国东部的长江三角洲（YRD），该地区是中国工业化程度最高、人口最密集的地区之一，污染和碳排放源复杂。本文总结了迄今为止该站点取得的主要科学成果，包括光谱反演表征和归一化、关键大气成分的总体演变特征、排放估计、CTM校验以及对大气污染来源和传输的研究结果。同时，本文还对合肥站的观测、科学研究和未来研究计划进行了展望。

《2、站点描述和遥感系统介绍》

2、站点描述和遥感系统介绍

《2.1 站点描述》

2.1 站点描述

FTS观测站位于中国东部高度工业化和人口密集的合肥郊区——科学岛[62‒63]。合肥站自2014年7月开始进行TCCON观测，2018年正式成为TCCON标准站（https://tccondata.org/）。合肥站点地理位置及气象场如图1所示。该站点尚被未纳入NDACC观测网，但自2015年起，其日常观测遵循NDACC标准公约[8,17]。合肥FTS观测站是一个跨学科的观测平台，由中国科学院安徽光学精密机械研究所（AIOFM-CAS）与中国科学技术大学（USTC）合作开发。该观测站是国家重大科研项目大气环境观测与模拟（AEOS）基础设施的一部分，包括多种环境监测仪器和模拟器，用于密集的大气科学研究。该观测站是国际公认的中国首个大气成分长期观测站，对评估中国东部地区的卫星观测、CTM模拟、空气质量以及人为和自然排放产生的大气污染物的长距离传输具有重要意义[63]。

《图1》

图1 中国合肥高分辨率FTS观测系统的主要部分。（a）IFS125HR光谱仪；（b）圆顶和太阳跟踪器[62‒63]。

合肥站具有独特的地理位置，是研究不同来源气团的物理过程及气团与大气化学相互作用的理想场所。合肥是长三角地区合肥-上海轴线城市群的起点。夏季为东南风，冬季为东北风，处于区域人为源和自然源的下风位置。因此，合肥站是研究经过充分混合和化学反应的大气成分的物理和化学过程的理想观测站[64]。与此同时，合肥地处中国另一个高度工业化和人口密集的京津冀（Beijing-Tianjin-Hebei, BTH）地区的南端。偏北冬季风可以将BTH地区密集的化石燃料和生物质燃烧源的区域性羽流带到合肥观测站[13]。此外，受夏季风影响，合肥观测站处于生物源排放较强的东南亚热带雨林区和西南亚热带雨林区的下风位置[15,64]。这些地区的生物质燃烧排放也会传输到合肥观测站。

《2.2 遥感系统介绍》

2.2 遥感系统介绍

合肥观测站包括一台高分辨率FTS（IFS125HR, Bruker, Germany）、一个太阳跟踪器和一个气象站（图1）。光谱仪放置在室内，太阳跟踪器和气象站直接安装在屋顶。两台空调、一台除湿机和一台空气净化器不间断运行，以获得恒温、干燥和清洁的室内条件[62‒63]。IFS125HR光谱仪具有卓越的精度和稳定性，已广泛应用于TCCON/NDACC观测网[22,65]。合肥观测站的IFS125HR有6个探测器、4个分束器，最大光程差为937.5 cm，对应的最大分辨率为0.001 cm^-1。所有探测器和分束器类型及其工作范围详见附录A中的表S1和表S2。光谱仪的光谱范围覆盖420 cm^-1到50 000 cm^-1。太阳跟踪器安装在一个由两个球形罩组成的圆顶内。穹顶可以保护太阳跟踪器不受雨、雪和大风等恶劣天气的影响。穹顶的旋转状态和两个球形盖的开关状态由三个特定的电机控制。太阳跟踪器通过调节两个镀铝折叠镜来跟踪太阳，并将阳光引导到FTS中。太阳跟踪器的太阳跟踪精度为±0.1 mrad [1]（1 rad = 57.3°）。安装在圆顶附近的气象站包括空气温度（±0.30 °C）、相对湿度（±3.0%）、气压（±0.1 hpa）、风速（±0.50 m·s^-1）、风向（±5.0°）、太阳辐射（±5.0%）和降雨的传感器。这些气象参数用于改进光谱反演。IFS125HR光谱仪配有多级涡旋真空泵，使干涉仪处于真空状态。真空条件提高了仪器的稳定性，最大程度减少了水蒸气的吸收干扰。

在常规观测中，按照附录A中表S3的顺序交替采集近红外（NIR）和中红外（MIR）太阳光谱。NIR光谱和MIR光谱的光谱范围分别为4000~11 000 cm^-1和500~8500 cm^-1，时间分辨率分别为150 s和288 s，光谱分辨率分别为0.020 cm^-1和0.005 cm^-1。根据天气条件不同，每个观测日的NIR光谱和MIR光谱数量位于11到103之间。采用InGaAs探测器和CaF₂分束器进行NIR光谱采集，采用InSb或碲镉汞（MCT）探测器、KBr分束器和7个滤光片进行MIR光谱采集。NIR光谱用于反演CO₂、CH₄、N₂O、HF、CO、H₂O、HCl和HDO的柱总量[17,62]。MIR光谱用于反演O₃、CO₂、NO₂、N₂、NO、HNO₃、CO、HCl、HF、N₂O、CH₄、C₂H₆、HCN、ClONO₂、H₂O、HDO、ClO、OCS、H₂CO、C₂H₂、C₂H₄、CF₄、SF₆、CH₃D、COF₂、CCl₃F、CCl₂F₂、CHClF₂、C₃H₈、HCOOH、CH₃OH、PAN和NH₃的VMR廓线和柱总量[18,66]。表1总结了这些气体的主要研究相关性和意义，包括地基验证、全球变暖、臭氧化学和区域污染。

《表1》

表1 合肥站所观测的大气成分的主要研究意义

Focus area	Significance/effect	Key atmospheric constituents
Validation	Satellite or CTM	CO₂, CH₄, N₂O, HF, CO, H₂O, HCl, HDO, O₃, NO₂, HNO₃, HCl, NO, C₂H₆, HCN, ClONO₂, H₂CO, C₂H₂, C₂H₄, SF₆, CCl₃F, CCl₂F₂, CHClF₂, HCOOH, CH₃OH, and NH₃
Global warming (radiative forcing)	Climate/environmental change	H₂O, HDO, CO₂, CH₄, CH₃D, N₂O, SF₆, O₃, C₂H₆, CCl₃F, CCl₂F₂, and CHClF₂
Ozone chemistry	Biological effects of ultraviolet-visible (UV-vis) exposure	O₃, ClO, HCl, ClONO₂, COF₂, NO, NO₂, and HNO₃
Regional pollution	Human health	CO, C₃H₈, HCN, H₂CO, COF₂, and O₃
	Photochemical smog	NO, HNO₃, NO₂, O₃, C₂H₂, H₂CO, C₂H₄, C₂H₆, PAN, and CO
	Acid rain	SO₂, NH₃, NO, HNO₃, and NO₂
	Oxidation efficiency of atmosphere	HCOOH, CH₃OH, CO, H₂CO, CH₄, NO, NO₂, and O₃
	Source of aerosols and precursors	SO₂, HNO₃, CF₄, NH₃, and OCS (for stratospheric aerosols)

《3、反演算法和表征》

3、反演算法和表征

《3.1 近红外光谱的反演算法及表征》

3.1 近红外光谱的反演算法及表征

采用TCCON标准分析软件GGG从近红外太阳光谱中反演大气成分的柱总量[17]。首先，GGG处理原始干涉图、校正相位误差和太阳强度变化，并将其转换为光谱。然后，生成气象参数和痕量气体的先验廓线，其中温度、湿度和压力先验廓线基于美国国家环境预报中心（National Centers for Environmental Prediction）和美国国家大气研究中心（NCEP/NCAR）的再分析结果，CO₂、CO、CH₄和N₂O先验廓线基于MkIV球载观测值、卫星大气化学实验的傅里叶变换光谱仪（ACE-FTS）数据和GLOBALVIEW数据。其次，GGG逐次迭代先验廓线生成模拟光谱，当测量光谱与模拟光谱相吻合时，停止迭代，推算出大气成分柱总量并进一步转换为柱平均干空气摩尔分数（DMF）。附录A中的表S4概括了基于NIR光谱反演所有气体的相关参数设置。

采用Wunch等[22]的方法来计算所有气体的反演误差，得出 X_CO $_{2}$ 、X_CH $_{4}$ 、X_N $_{2}$ _O、X_CO、X_H $_{2}$ _O、X_HDO和X_HF的反演误差分别小于0.25%（约1 ppm）、0.50%（约5 ppb）、1.00%（约3 ppb）、4.00%（约10 ppb）、1.30%、2.00%和8.00%（约4 ppb）。合肥观测站的柱总量平均核函数值（TAK）如图2所示。由于平流层光谱线的压强展宽比对流层窄，导致平流层的TAK随着空气质量的增加而减少。基于合肥观测站的NIR光谱反演的X_CO $_{2}$ 、X_CH $_{4}$ 、X_N $_{2}$ _O、X_CO、X_H $_{2}$ _O、X_HDO、X_HCl、X_HF时间序列如图3所示。作为一个正式的TCCON观测站，合肥观测站的所有数据集都存档在TCCON数据库（https://tccondata.org/）中，供全球开放访问。

《图2》

图2 基于合肥观测站的NIR光谱反演大气柱总量的平均核函数（TAK）。

《图3》

图3 基于合肥观测站的NIR光谱反演的气体时间序列。

《3.2 中红外光谱的反演算法及表征》

3.2 中红外光谱的反演算法及表征

采用NDACC最新的标准分析软件SFIT4反演MIR光谱[16,68]。如附录A中表S3所示，合肥站的MIR光谱能够反演30多种大气成分的柱总量和VMR廓线，其中CH₄、N₂O、CO、HCN、C₂H₆、O₃、ClONO₂、HNO₃、HF和HCl为NDACC强制观测大气成分。所有NDACC气体的反演参数设置如附录A中表S5 [69‒71]所示。除H₂O以外的所有气体的先验廓线来自全大气气候模型（WACCM）1980年至2020年模拟结果的统计平均值[8]。压力、温度和H₂O的先验廓线基于NCEP/NCAR的6 h再分析结果。根据气体种类不同，光谱参数来自于HITRAN 2008或HITRAN 2012数据库[8,10,67,71]。先验协方差矩阵_a的对角线值设置为WACCM v6模拟值的标准偏差，噪声协方差矩阵_ε的对角线值设置为光谱信噪比平方的倒数。_a和_ε的非对角线值均设为零，反演中的仪器函数（ILS）来自于实际测量值[72‒73]。

基于合肥观测站的MIR光谱反演的部分气体的平均核函数值（AK）及其面积如图4所示。HCN、H₂CO、SF₆、CO、C₂H₆的高灵敏范围位于对流层，HF、NO₂、HCl的高灵敏范围位于平流层，O₃和CH₄的高灵敏范围同时包括平流层和对流层[74]。CH₄、NO₂、CO、HCN、C₂H₆、O₃、H₂CO、SF₆、HF和HCl的整层反演自由度（DOFS）分别为3.5、1.0、3.0、1.1、1.5、4.8、1.1、1.0、1.3和1.5。

《图4》

图4 利用合肥站的MIR光谱反演的对流层气体（C₂H₆、CH₄、CO、H₂CO、HCN和SF₆）和平流层（HCl、HF、NO₂和O₃）气体的平均核函数（AK）（彩色细线）。每种气体的AK面积用黑色虚线表示，并按0.1倍缩放。

采用的方法估算合肥站点的MIR光谱反演误差，CH₄、NO₂、CO、HCN、C₂H₆、O₃、SF₆、H₂CO、HF和HCl的反演误差分别为4.5%、12.4%、4.7%、11.3%、6.2%、5.7%、5.8%、12.3%、5.1%和5.2%。基于合肥站点的MIR光谱反演的部分气体柱总量时间序列如图5所示。这些气体都具有明显的日变化和季节性。

《图5》

图5 基于合肥站点的MIR光谱反演的部分气体柱总量时间序列，垂直误差条表示反演不确定度。

《4、合肥观测站取得的主要研究成果》

4、合肥观测站取得的主要研究成果

合肥观测站已经连续运行7年，反演了很多种大气成分的时间序列。国内外研究学者利用这些数据集开展了一系列研究，包括高分辨率光谱反演与归一化、关键大气成分的总体演变特征、排放估计、卫星和大气化学模式评估、空气污染和化学反应机制等。表2 [8,10,13,15,62‒63,66‒67,71,73,75‒81]概括了基于2014—2020年合肥站的观测数据取得的研究成果清单，按主要发现、数据时间段、主题、波段、大气层范围和主题进行了分类。合肥站的主要成果如下：

《表2》

表2 基于2014—2020年合肥站的观测数据的研究论文清单，按主要发现、数据时间段、主题（主要气体）、波段、大气层范围和主题分类

Authors	Key findings	Period	Theme	Waveband	Domain	Topic
Wang et al. [62]	X_CO $_{2}$ at the Hefei station showed an obvious seasonal cycle but X_CO displayed no clear seasonality. X_CO $_{2}$ at the Hefei station are in good agreement with both the Greenhouse gases observing satellite (GOSAT) and Orbiting Carbon Observatory-2 (OCO-2) observations	2014.06‒2016.04	CO₂, CO	IR	Troposphere and stratosphere	Overall characteristics and satellite evaluation
Tian et al. [63]	Diurnal X_CH $_{4}$ exhibits an increasing change rate in summer, and a decreasing change rate or remains constant in other seasons. X_CH $_{4}$ shows a seasonal maximum in summer and a minimum in winter. X_CH $_{4}$ from 2014 to 2017 showed an increasing change rate of (0.56% ± 0.15%) a^-1	2014.06‒2017.07	CH₄	IR	Troposphere and stratosphere	Retrieval harmonization, satellite evaluation, and model evaluation
Sun et al. [75]	Attenuators cause negligible ILS drifts and can be used to adapt the received intensity of the detector. The usage of different sizes of field stop may cause dramatically ILS drifts due to inconsistent mechanical errors among the selected field stops	2015	ILS	IR & MIR	Troposphere and stratosphere	ILS characterization
Sun et al. [73]	Tropospheric gases are less sensitive to ILS drift than the stratospheric gases, and, for the same level of drift, the negative ILS influence is smaller than the positive ILS. To limit the difference in total column of ClONO₂ within 10% and other gases within 1%, the upper limits of ILS for CH₄, N₂O, CO, HCN, C₂H₆, O₃, ClONO₂, HNO₃, HF, and HCl are suggested	2015.08‒2016.08	CH₄, N₂O, CO, HCN, C₂H₆, O₃, ClONO₂, HNO₃, HF, HCl	MIR	Troposphere and stratosphere	ILS characterization and retrieval harmonization
Sun et al. [66]	The effect of ILS drift on partial column is gas and altitude dependent. To suppress the influence on partial column within 2%, the upper limits of ILS for CH₄, N₂O, O₃, and CO are suggested	2015.08‒2016.08	CH₄, N₂O, O₃, CO	MIR	Troposphere and stratosphere	ILS characterization and retrieval harmonization
Sun et al. [8]	The level and variability of tropospheric O₃in spring and summer (March-April-May/June-July-August; MAM/JJA) are higher than those in autumn and winter (September-October-November/December-January-February; SON/DJF). Reductions in VOC and NO_x could effectively mitigate O₃ pollution in SON/DJF and MAM/JJA seasons, respectively	2014.09‒2017.09	O₃	MIR	Troposphere	Overall characteristics and source attribution
Shan et al. [76]	CO emissions over the megacity Hefei are calculated by utilizing the enhancement ΔCO/ΔCO₂ ratios derived from the NIR measurements at the Hefei station. The discrepancy between the CO emissions inferred from emission inventories and those from the measurements is also investigated	2015.09‒2017.08	CO	IR	Troposphere	Emission estimates
Yin et al. [71]	Stratospheric NO₂ column at the Hefei station shows a seasonal maximum in June and a minimum in January. GEOS-Chem model can simulate the burden, seasonal cycle, and interannual trend of stratospheric NO₂ over the polluted eastern China	2015.07‒2019.01	NO₂	MIR	Stratosphere	Overall characteristics, satellite evaluation, and model evaluation
Hedelius et al. [77]	FTS measurements at the Hefei station can be used as one of the reference datasets for evaluation of measurements of pollution in the troposphere (MOPITT) X_CO version 7 retrievals	2015.09‒2016.12	CO	NIR	Troposphere and stratosphere	Satellite evaluation
Yin et al. [78]	HCl total column at the Hefei station shows a seasonal maximum in January and a minimum in September. GEOS-Chem model can simulate the burden, seasonal cycle, and interannual trend of stratospheric HCl over the polluted eastern China	2015.07‒2019.04	HCl	MIR	Troposphere and stratosphere	Overall characteristics and model evaluation
Oshio et al. [79]	FTS measurements at the Hefei station can be used as one of the reference datasets for bias correction of the CH₄/CO₂total column ratio retrieved from GOSAT observations	2015.09‒2016.12	CO₂, CH₄,	NIR	Troposphere and stratosphere	Satellite evaluation
Sun et al. [10]	Tropospheric columns of HCN at the Hefei station show pronounce seasonal variations with three seasonal peaks in May, September, and December, respectively. Elevated BB emissions in Oceania dominated the HCN enhancements in the second half of 2015. Elevated BB emissions in SEAS dominated the HCN enhancements in the first half of 2016	2015‒2018	HCN, CO	MIR	Troposphere	Overall characteristics, model evaluation, and source attribution
Shan et al. [80]	GEOS-Chem model can simulate the burden, seasonal cycle, and interannual trend of CO₂ over the polluted eastern China	2015.07‒2019.12	CO₂	NIR	Troposphere	Overall characteristics and model evaluation
Sun et al. [13]	H₂CO photolysis plays a significant role in OH generation over eastern China. The NMVOCs related summertime H₂CO enhancements were largely caused by the NMVOCs emissions within eastern China. Due to the increase in photochemical H₂CO resulting from the increases in both NMVOCs and CH₄, H₂CO from 2015 to 2019 showed an increasing trend	2015‒2019	H₂CO	MIR	Troposphere and Stratosphere	Overall characteristics, model evaluation, and source attribution
Sun et al. [15]	The observed C₂H₆ variability at the Hefei station was mainly driven by the variation of C₂H₆ emissions within China. Decrease in C₂H₆ from 2015 to 2020 points to an improvement of air quality in China, which is due to the reduction in transported and local C₂H₆ emissions in recent years	2015‒2020	C₂H₆	MIR	Troposphere	Overall characteristics, model evaluation, and source attribution
Yin et al. [67]	The X_HF time series at the Hefei station reached a seasonal maximum in March and a minimum in September. HF columns showed a decreasing change rates of (‒0.38% ± 0.22%) a^-1between 2015 and 2020	2015‒2020	HF	MIR	Stratosphere	Overall characteristics and retrieval harmonization
Shan et al. [81]	The stratospheric HNO₃ column at the Hefei station reached a seasonal maximum in March and a minimum in September, and showed a decreasing change rate of (‒9.45% ± 1.20%) a^-1 between 2017 and 2019	2017‒2029	HNO₃	MIR	Stratosphere	Overall characteristics and satellite evaluation

《4.1 仪器函数表征及反演归一化》

4.1 仪器函数表征及反演归一化

TCCON/NDACC观测网的大多数FTS观测站都使用不同尺寸的光阑或衰减器来调整探测器的接收光强。这些行为可能会改变高分辨率FTS的仪器函数，增大反演偏差。Sun等[75]首次量化了高分辨率FTS在不同衰减器下的仪器函数漂移。Sun等[75]在干涉仪的不同位置插入了5个衰减器，并进行了仪器函数测量，发现这5种衰减器造成的仪器函数漂移可以忽略不计，都可用于调整探测器的接收光强。然而，使用不同尺寸的光阑可能会产生严重的仪器函数漂移。

Sun等[66,73]进一步模拟了不同类型的仪器函数漂移，量化了它们对所有NDACC气体反演的影响。Sun等[66,73]发现，当仪器函数漂移为10%时，CH₄、N₂O、CO、HCN、C₂H₆、O₃、ClONO₂、HNO₃、HF和HCl的柱总量分别变化了0.04%、0.20%、2.10%、0.75%、1.10%、1.90%、23.00%、0.70%、3.00%和4.00% [73]。为了使ClONO₂柱总量偏差控制在10%以内，其他气体偏差控制在1%以内，O₃、ClONO₂、HNO₃、HF、HCl、CO、HCN和C₂H₆的仪器函数正向漂移上限分别为6%、5%、15%、5%、5%、5%、13%和9%；O₃、HF和HCl的仪器函数负漂移上限分别为6%、12%和12%；对于CH₄和N₂O，仪器函数漂移的影响可以忽略不计[73]。

Sun等[66]研究得出仪器函数漂移对气体偏柱的影响取决于气体种类和高度。仪器函数漂移对CH₄和N₂O的影响小于对O₃和CO的影响，为了将偏柱的影响抑制在2%以内，建议反演CO和O₃时使用精确的仪器函数漂移值。CH₄在0~7 km高度和16~37 km高度的仪器函数正向漂移上限分别为13.0%和1.5%。在海拔7~16 km和16~37 km范围内，仪器函数负向漂移上限分别为12.0%和4.5%。仪器函数负向漂移对0~7 km处CH₄偏柱的影响和正向漂移对7~16 km处CH₄偏柱的影响可以忽略不计[66]；N₂O在0~5.0 km、5.0~11.5 km、11.5~20.0 km和20.0~35.0 km高度范围内的仪器函数正向漂移上限分别为6.0%、4.0%、1.5%和1.5%。在11.5~20.0 km和20.0~35.0 km高度范围内，仪器函数负向漂移上限分别为2.0%和1.5%。然而，仪器函数负向漂移对0~5.0 km和5.0~1.5 km N₂O偏柱的影响可以忽略不计[66]。

合肥观测站采集的NIR和MIR光谱基于不同的滤光片、分束器和探测器，针对它们的反演采用了不同的反演波段、光谱建模方案和迭代方法，这可能会导致反演偏差。为了诊断NIR和MIR光谱反演的HF数据集是否可以合并为一个可靠的数据集，Yin等[67]对这两个数据集进行了比较，发现NIR和MIR X_HF数据具有一致的季节性，但NIR X_HF日平均值普遍低于MIR数据，平均差值为(6.90 ± 1.07) pptv（按体积计算，万亿分之一）[67]。因此，Yin等[67]在结合这两个数据集进行研究之前，对NIR数据集增加了6.90 pptv的固定偏差。

在合肥观测站的早期阶段，NIR光谱使用了不同的光阑和探测器。同时，这些早期的NIR光谱也缺乏同步气象数据，并且由于观测协议不完整而受到计时误差的影响。所有这些观测缺陷都有可能产生系统偏差。为了获得一致的X_CH $_{4}$ 长期时间序列，Tian等[63]进行了一系列修正以消除这些潜在的系统偏差，最终修正了2015年7月之前X_CH $_{4}$ 的11.0 ppb系统偏差[63]。

《4.2 关键大气成分的演变特征》

4.2 关键大气成分的演变特征

合肥观测站的连续观测数据已经用于研究中国东部地区主要大气成分的总体演变特征。Wang等[62]基于两年的NIR测量结果研究了CO₂和CO的季节变化规律，发现X_CO $_{2}$ 具有明显的季节周期，而X_CO没有明显的季节性。Wang等[62]进一步分析了X_CO和X_CO $_{2}$ 之间的关系，发现这两种气体仅在个别观测时间段具有相关性，且随盛行风条件而变化[62]。Tian等[63]基于合肥站2014至2017年的NIR观测数据推导了X_CH $_{4}$ 的年变化率为(0.56 ± 0.15%) a^-1。X_CH $_{4}$ 在夏季具有较强的日增长趋势，而其他季节日变化率呈下降趋势或保持不变。X_CH $_{4}$ 时间序列的季节最大值出现在夏季，最小值出现在冬季。1月份的X_CH $_{4}$ 年平均值比8月份低30~55 ppb。

Sun等[8,10,13,15]利用合肥站不同时间尺度的MIR观测数据研究了中国东部地区O₃、HCN、H₂CO和C₂H₆的季节性和年际变化趋势。对流层O₃（0~12 km）具有6月最大、12月最小的季节性特征。6月对流层O₃柱浓度比12月高50%，春夏季节的对流层O₃日变化一般大于秋冬季节。这与Vigouroux等[82]的观点不同，后者认为相对干净地区的对流层O₃通常在春季达到季节性最大值。

合肥站HCN对流层柱浓度具有明显的季节变化特征，5月、9月和12月分别出现三个季节峰值[10]。总体而言，合肥站HCN的对流层柱浓度高于NDACC其他站点的观测值。2015年9月至2016年7月HCN对流层柱浓度比其他年份同期值大5%~46%。合肥站X_H $_{2}$ _CO时间序列在7月达到季节最大值，在1月达到季节最小值，1月份的X_H $_{2}$ _CO值比7月份低42.9%。2015至2019年合肥站X_H $_{2}$ _CO时间序列呈(2.37% ± 0.70%) a^-1的正变化趋势。

Sun等[15]发现C₂H₆对流层柱平均干空气摩尔分数（troDMF）在12月达到季节最大值，7月达到最小值，2015—2020年呈负增长趋势，变化率为(-2.6% ± 1.33%) a^-1。由于夏季对流层OH氧化能力高于冬季[15]，C₂H₆ troDMF具有冬季最大值和夏季最小值的季节性特征。排放和气象因素对夏秋季节的C₂H₆ troDMF具有负向影响，对春冬季节的C₂H₆ troDMF具有正向影响。

Yin等[71]研究发现，合肥站平流层NO₂在上半年呈递增趋势，下半年呈递减趋势，6月出现季节性最大值，1月出现季节性最小值。6月份的NO₂平流层柱浓度比1月份平均高39.20% ± 8.95%。Yin等[67‒78]和Shan等[81]进一步研究得出平流层HCl、HF和HNO₃的季节变化趋势与对流层顶高度的季节性波动有关。平流层HCl、HF和HNO₃柱浓度在1月份出现季节性最大值，在9月份出现季节性最小值[67,78,81]。平流层NO₂、HCl、HF和HNO₃柱浓度变化率分别为(-0.34% ± 0.05%) a^-1（2015—2018年）、(-1.83% ± 0.13%) a^-1（2015—2019年）、(-0.38% ± 0.22%) a^-1（2015—2020年）和(-9.45% ± 1.2%) a^-1（2017—2019年）。

《4.3 排放评估》

4.3 排放评估

对于CO和CO₂有共同来源的区域，可以通过研究CO和CO₂的相关性来推断CO的排放量。Shan等[76]利用合肥站NIR光谱反演的ΔCO/ΔCO₂比值计算了超大城市合肥的CO排放。Shan等[76]还深入研究了从排放清单推断的CO排放量与从测量结果推断的CO排放量之间的差异。为此，Shan等[76]比较了合肥站2015年9月至2017年8月通过FTS、卫星、现场测量和排放清单获得的ΔCO/ΔCO₂比值。在Shan等[76]的研究中使用了全球大气研究排放数据库（EDGAR）和北京大学（PKU）的排放清单。研究得出，基于两种排放清单的ΔCO/ΔCO₂比值远远大于基于地基FTS、原位仪器和卫星测量推算的ΔCO/ΔCO₂比值。

Shan等[76]进一步利用FTS观测值推算的ΔCO/ΔCO₂比值，以及从EDGAR和PKU清单中得出的CO₂排放量来推断CO排放量。基于EDGAR清单的CO₂排放和基于FTS观测的ΔCO/ΔCO₂比值计算的CO排放，2015—2016年约为(11.27 ± 0.91) Tg·a^-1，2016—2017年约为(12.35 ± 0.74) Tg·a ^-1。基于PKU清单的CO₂排放和基于FTS观测的ΔCO/ΔCO₂比值估算的同期CO排放量分别约为(10.96 ± 0.88) Tg·a^-1和(11.95 ± 0.71) Tg·a^-1。由于PKU和EDGAR清单可以很好地再现实际的CO₂排放量，基于FTS ΔCO/ΔCO₂比值与两份清单计算的CO排放量吻合较好。

《4.4 卫星地基校验》

4.4 卫星地基校验

合肥站的FTS连续观测结果为卫星遥感地基验证提供了一系列有价值的数据集。为了比较卫星和地面FTS的测量结果，考虑了它们在垂直分辨率和先验廓线上的差异。首先，以FTS先验廓线的垂直间隔对卫星数据进行重新网格化，减少两者先验廓线的相关差异[56]。然后通过采用卫星数据的平均核函数对FTS数据进行平滑处理，最小化垂直分辨率的相关差异。最后，再将重新网格化和平滑处理的两种观测数据进行对比。

Wang等[62]发现合肥站的地基FTS X_CO $_{2}$ 时间序列与日本温室气体观测卫星（GOSAT）和美国碳观测卫星（OCO-2）观测结果吻合较好[62]。与地面FTS观测数据相比，GOSAT和OCO-2 X_CO $_{2}$ 时间序列的日均值偏差分别为(-0.52 ± 1.63) ppm和(0.81 ± 1.73) ppm。GOSAT和OCO-2观测可以再现FTS X_CO $_{2}$ 时间序列的日中位数，相关系数（r）分别为0.79和0.83。Yuan等[63]发现合肥站GOSAT与FTS X_CH $_{4}$ 数据的差异（GOSAT与FTS观测值之差）均在两个数据集的误差范围内，且与季节无关，其平均值为(1.60 ± 13.0) ppb（0.09% ± 0.70%）[50,83]。GOSAT和FTS观测到的X_CH $_{4}$ 的季节周期符合较好，相关系数（r）为0.77。

Yin等[71]发现臭氧监测仪（OMI）卫星观测到的2015—2018年之间的平流层NO₂时间序列与合肥站同期的地基FTS数据吻合较好，相关系数（r）为0.841 [71]。OMI与FTS的年平均差值（OMI减去FTS）为(1.48 ± 5.33) × 10¹⁴ molecules∙cm^-2（4.82% ± 17.37%）[71]，在两个数据集的不确定度范围内。OMI与地基FTS观测的2015—2018年平流层NO₂具有一致的负增长趋势，约为(-0.91% ± 0.09%) a^-1。此外，Oshio等[79]将合肥站的FTS测量数据作为参考数据集之一，对GOSAT观测的CH₄/CO₂总柱比进行了偏差校正，Hedelius等[77]将FTS测量数据用于评估对流层污染监测仪（MOPITT）X_CO v7的反演结果。

Shan等[81]将合肥站观测的平流层HNO₃和HCl柱总量与MLS卫星数据进行了对比，发现MLS卫星数据与FTS数据具有相似的季节变化和年际变化率，两个数据集的平流层HNO₃和HCl柱浓度的相关系数（r）分别为0.87和0.88。卫星观测的平流层HNO₃柱和HCl柱与傅里叶变换红外光谱（FTIR）数据的平均偏差分别为-8.58% ± 12.22%和4.58% ± 13.09%。

Wang等[62]、Tian等[63]、Yin等[71]、Hedelius等[77]、Oshio等[79]、Shan等[81]的结果验证了合肥观测站能够探测对流层和平流层关键大气成分的日、季、年际变化，并验证了卫星观测的有效性。

《4.5 大气化学模式评估》

4.5 大气化学模式评估

合肥站的连续FTS观测数据已经用于评估大气CTM在中国东部地区的模拟性能。由于CTM的垂直分辨率与地基FTS存在差异，评估之前，对CTM廓线进行了平滑修正[84]。首先，将廓线重新网格化到FTS垂直高度。然后，再以FTS的先验廓线和AK对重新网格化的廓线进行平滑[68,84]。最后，再将重新网格化和平滑后的CTM数据与地基FTS数据进行比较。

利用合肥站获得的CH₄、NO₂、H₂CO、C₂H₆、CO₂等多种气体的FTS时间序列对GEOS-Chem模型进行验证。GEOS-Chem模式总体上可以再现地基FTS观测到的这些气体的季节周期和年际趋势。对于X_CH $_{4}$ 、NO₂平流层柱、X_H $_{2}$ _CO、C₂H₆ troDMF和CO₂ troDMF，GEOS-Chem与FTS的相关系数（r）分别为0.860、0.858、0.780、0.880和0.890 [13,15,63,71,80]。然而，GEOS-Chem模式很难再现FTS观测到的上述气体的日变化[13]。例如，夏季GEOS-Chem和FTS之间的X_CH $_{4}$ 差异比冬季大35~55 ppb [63]。对于NO₂平流层柱，秋季GEOS-Chem和FTS之间的差异要明显小于其他季节[71]。对于C₂H₆ troDMF，GEOS-Chem模式[15]在7月高估FTS月平均观测值35.6%，在12月低估FTS月平均观测值17.4%。这些差异主要是因为GEOS-Chem模式在粗空间分辨率下模拟的垂直混合和水平传输方案存在一定不确定性，难以与单点柱状观测相匹配，也与排放清单和化学机制的不确定性有关。

GEOS-Chem与FTS的X_CH $_{4}$ 、NO₂平流层柱、X_H $_{2}$ _CO、C₂H₆ troDMF和CO₂ troDMF差异分别为(-8.09 ± 17.83) ppb（0.45% ± 0.95%）、(2.36 ± 2.33) × 10¹⁴ 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%）和(-1.68 ± 1.23) ppbv（-0.42% ± 0.31%），均在各自气体的FTS反演误差范围内。结果验证了GEOS-Chem能够模拟中国东部污染地区的CH₄、NO₂、H₂CO、C₂H₆和CO₂浓度、季节循环和年际变化趋势[13,15,63,71,80]。

《4.6 大气污染、来源与传输》

4.6 大气污染、来源与传输

Sun等[8]利用合肥站观测的对流层O₃、H₂CO和CO柱浓度，OMI提供的对流层NO₂柱浓度，结合气团轨迹，研究了合肥地区的O₃光化学生成机制[8]。Sun等[8]得出，中国北部的山东和河北，东部的江苏和安徽，中部的陕西、湖北和湖南，西北部的河南和山西等地的排放都对合肥地区的O₃生成有影响[8]。因为春夏季节的气团更多来自高度工业化和人口密集地区，春夏季节的对流层O₃浓度要明显高于秋冬季节。Sun等[8]进一步以H₂CO/NO₂对流层柱比值为指示因子，探索了O₃（PO₃）的光化学生成敏感性。Sun等[8]发现春夏季节PO₃主要受NO_x控制，秋冬季节主要受VOC控制或VOC-NO_x混合控制。

GEOS-Chem模式捕捉到了合肥站CO、H₂CO、C₂H₆的浓度变化，因此Sun等将其用于这些气体的溯源研究[10,13,15]。Sun等[10]发现对流层HCN（ΔHCN）和CO（ΔCO）的增强具有很强的相关性，表明对流层CO和HCN的增强具有相同的来源[10]。利用这一优势，Sun等[10]利用GEOS-Chem标记的CO模拟确定了HCN的来源，得出HCN在5月的季节最大值主要归因于东亚和东南亚（SEAS）（41.0% ± 13.1%）、非洲（AF）（22.0% ± 4.7%）及欧洲和北亚（EUBA）（21.0% ± 9.2%）的生物质燃烧排放。HCN在9月的季节最大值主要来自EUBA （38.0% ± 11.3%）、SEAS（14.0% ± 3.3%）、AF（26.0% ± 6.7%）和北美（NA）（13.7% ± 8.3%）的生物质燃烧排放。HCN在12月的季节最大值主要受AF（36.0% ± 7.0%）、NA（18.7% ± 5.2%）和EUBA（21.0% ± 5.2%）影响。

Sun等[13]利用GEOS-Chem模型量化了不同地理区域和排放类别对合肥站夏季H₂CO增强的贡献[13]。结果表明，非甲烷挥发性有机化合物（NMVOC）和CH₄氧化对合肥市夏季H₂CO增强的贡献分别为56.73%和43.27%；与夏季H₂CO增强有关的NMVOC主要是由中国东部NMVOC排放引起的。由于NMVOC和CH₄的增加导致光化学H₂CO的增加，因此2015—2019年H₂CO呈现增加趋势[13]。基于Sun等[13]早期相似的灵敏度研究方法，Sun等[15]得出中国东部地区人为（生物燃料+化石燃料）排放占C₂H₆浓度的49.2%，自然（生物质燃烧+生物排放）排放占37.0%。合肥站观测到的C₂H₆变化主要归因于中国境内C₂H₆排放变化（74.0%），其中华北、华东和华中地区贡献最大（57.6%）。2015—2020年C₂H₆呈下降趋势，表明中国的空气质量有所改善，这是由于近年来外源传输和本地排放的C₂H₆量有所减少[15]。

《5、结论及未来展望》

5、结论及未来展望

本研究总结了中国首个TCCON地基高分辨率FTS观测站——合肥观测站，2014年以来取得的主要研究成果。这些研究提高了在高分辨率光谱反演和表征、气候变化、排放估计、卫星和大气化学模式验证、大气污染的来源和传输方面的认识。关于物理过程和化学机制的新发现将会被应用到CTM中。反过来，经验证的CTM将被用于进一步改善对大气产生和消耗过程的认识。根据合肥站所取得的主要研究成果和总体构想，未来合肥站预计将开展以下方面科学研究。

（1）利用地基柱总量或廓线反演结果开展不同尺度的大气研究，并在不同观测平台之间提供标准转换，将地基FTS数据集转换为与世界气象组织（WMO）尺度相同的验证尺度。目前还未将合肥观测站的地基FTS数据进行基于机载观测的原位廓线测量结果的验证。这种校验将提高地基FTS观测的精度水平。本研究简单地假设其他TCCON站点的有效校准值也适用于合肥站，并根据Wunch等[17]的观测偏差来定标本文NIR观测结果。因此，未来的研究期待使用机载廓线测量来进一步提高合肥站点FTS观测的准确性。

（2）为了理解大气成分与气候变化之间的相互作用，需要对很多种大气成分进行长期测量。因此，合肥站将继续高质量、稳定测量所有大气成分。未来更多的物种将被研究，如PAN、OCS、NH₃、含氯氟烃（CFC），随着反演技术和新光谱数据的发展，这些物种数量将随着时间的推移而不断增加。与TCCON/NDACC观测网的其他光谱仪相比，合肥站的IFS125HR光谱仪覆盖了紫外光谱范围，但在常规观测中没有尝试采集紫外光谱。计划在未来采集紫外光谱，这将进一步增加可观测大气成分的数量。合肥站将为中国政府制定绿色经济政策、实现碳中和和《巴黎协定》目标提供科学支持。

（3）作为大气环境观测与模拟基础设施的一部分，高精度FTS数据集将被广泛用于验证大气环境观测与模拟框架内的其他仪器。尽管IFS125HR光谱仪在精度和稳定性等方面表现出色，但也有局限性。IFS125HR是一种昂贵而笨重的光谱仪，它的运行需要大量的基础设施。定期维护其光学准直度既繁琐又受限于操作人员的专业技能。为了克服这些问题，国际上开发了更小、更便宜且更容易运输的光谱仪，如EM27/SUN和VERTEX-80/SUN。已有研究表明EM27/SUN和VERTEX-80/SUN对于某些气体（如CO₂、CO和H₂CO）具有与IFS125HR相似的探测特性，但具有更好的流动性[70,85]。在使用EM27/SUN和VERTEX-80/SUN仪器计算城市或工业设施的排放通量之前，将利用高分辨率FTS的数据集对其进行验证。经验证的EM27/SUN和VERTEX-80/SUN观测数据可以补充中国现有的高分辨率FTS观测数据集。此外，合肥站点的FTS数据集还将用于验证未来的中国环境卫星。

（4）合肥站积累了近7年的高质量大气成分观测数据。这些数据将与其他平台或CTM模拟数据一起使用，并将继续提供给国内外科研团体，以更好地了解中国和全球的气候变化、空气污染、来源、传输及其对农作物和人类健康的影响。

（5）合肥站取得的研究成果得益于合肥站点的研究团队与TCCON/NDACC观测网的密切合作。合肥站未来将更多地参与国际TCCON/NDACC的观测和科学活动，为深入研究全球范围内的碳循环和气候变化做出贡献。

（6）利用合肥站的高光谱分辨率和高精度地基FTS数据集研究新的光谱特性。光谱特性的深入研究将有助于提高大气成分的反演精度和准确度。本文研究团队将与光谱学研究团体合作，改进大气吸收和消光的解释。

最后，合肥站的研究成果表明，新建一个观测站可以在区域和全球范围内提高人们对气候变化、空气污染和传输的认识。考虑人类正面临着前所未有的气候变化挑战，有必要在全球不同地区，特别是在只有少量地基观测的非洲、南美洲和亚洲建立更多这样的观测站。只有充分的测量和各站点之间的密切合作，人们才能更好地了解大气，并利用气候模式预测“未来气候”。

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.

展示更多