大气甲烷探测进展与全球甲烷分布分析
1.
2.
3.
Development of Atmospheric Methane Observation and Distribution of Global Methane
1.
2.
3.
通讯作者:
收稿日期: 2021-05-17 修回日期: 2022-02-09
| 基金资助: |
|
Received: 2021-05-17 Revised: 2022-02-09
作者简介 About authors
刘双慧(1996-),女,山东德州人,硕士研究生,主要从事定量遥感应用研究E⁃mail:
关键词:
Keywords:
本文引用格式
刘双慧, 李小英, 曹西凤, 张新苑.
Liu Shuanghui, Li Xiaoying, Cao Xifeng, Zhang Xinyuan.
1 引 言
甲烷(CH4)是大气中具有辐射活性和化学活性的气体[1],它不仅影响着地—气系统的辐射平衡过程,而且在大气和平流层化学过程中也起着非常重要的作用。它可以吸收大气长波辐射,使大气升温并改变大气层顶的净辐射通量。CH4在平流层中与氧气(O2)发生氧化反应生成水汽(H2O),在对流层中与羟基反应生成臭氧(O3)。虽然CH4在大气中的含量不及二氧化碳(CO2),但单位质量CH4的全球变暖潜力却是CO2的28倍[2]。随着人类文明的发展,CH4的总量一直在增加。工业革命前大气中的CH4浓度为700 ppbv[3],到20世纪90年代它的浓度达到1 714 ppbv[4]。按照来源可以将大气中的CH4源划分为两类。一类是自然源,如沼泽和湿地的厌氧腐烂、极地冻土融化等;另一类是人为源,有煤炭开采、生物质燃烧、化石燃料燃烧、水稻种植、动物肠道发酵、垃圾填埋等[5]。
大气CH4探测主要有地基探测、空基探测和星基探测三种方式。CH4的系统观测开始于19世纪80年代世界气象组织(WMO, World Meteorological Organization)的全球大气观测(GAW,Global Atmosphere Watch)计划。1990年以来,卫星遥感在大气CH4等痕量气体成分监测中的应用越来越广泛。白文广[6]阐述了2010年前发射的CH4探测器TES(Tropospheric Emission Spectrometer)、ACE-FTS(Atmospheric Chemistry Experiment Fourier Transform Spectrometer)、MIPAS(Michelson interferometer for passive atmospheric sounding)、IASI(Infrared Atmospheric Sounding Interferometer)、MOPITT(Measurements of Pollution in the Troposphere)、SCIAMACHY(Scanning Imaging Absorption Spectrometer for Atmospheric Chartography)的仪器参数、指标性能和反演方法。近红外波段和中红外波段均可探测CH4浓度,Kamieniak等[7]介绍了这两个波段探测CH4浓度的优缺点。Jacob等[8-9]分析了CH4天底传感器的仪器参数,指出探测仪的空间分辨率和甲烷探测精度均不断提高。Wang等[10]对监测中国大气环境的卫星传感器的仪器参数进行了描述,包括可进行CH4探测的仪器如SCIAMACHY、AIRS(Atmospheric Infrared Sounder)、TANSO-FTS(Thermal and Near-infrared Sensor for carbon Observation Fourier-Transform Spectrometer)、TANSO-FTS2、GMI(Greenhouse gas Monitoring Instrument)。上述研究主要是对CH4天底观测模式进行介绍,而对新兴的临边/掩星观测模式及地基、空基和星基3种探测方式之间的对比研究较少。
首先介绍地基探测网络和空基探测器,同时对不同观测模式下(天底、临边/掩星)高光谱传感器仪器参数及在大气CH4探测方面进行论述,综合对比了地基探测、空基探测与星基探测3种探测方式的优缺点、适用性及精度;然后对剥洋葱算法、权重函数差分吸收算法、最优化算法等3种CH4浓度反演方法的原理、特点及研究进展进行了梳理;最后对300 hPa、150 hPa全球大气CH4浓度的分布特征和变化趋势进行分析,为后续CH4研究工作的开展提供借鉴。
2 大气甲烷探测方式分析
2.1 地基探测
CH4的系统观测开始于19世纪80年代世界气象组织(WMO)的全球大气观测(GAW)计划。国际上主要的CH4地面观测网络是总碳柱观测网络(TCCON, Total Carbon Column Observing Network)和大气成分变化探测网络(NDACC, Network for the Detection of Atmospheric Composition Change)。
2.1.1 TCCON
TCCON是地面傅里叶变换光谱仪全球监测网络,成立于2004年。它可以测量CO2、CO、CH4、N2O等大气成分柱浓度,可为星基传感器(SCIAMA-CHY、TES、AIRS、OCO-2、GOSAT)提供验证数据。TCCON所有观测站点均采用Burker IFS 125HR光谱仪进行测量,仪器的光谱分辨率为0.02 cm-1,用于探测CH4的中心波长为5 938 cm-1、6 002 cm-1、6 076 cm-1。TCCON基于光谱拟合和逐行计算反演算法(GFIT, spectral fitting and line-by-line retrieval algorithm)对CH4柱平均干空气摩尔分数进行反演。该算法首先是计算CH4柱浓度和同一光谱中的O2柱浓度,然后计算CH4柱平均干空气摩尔分数,它的优点是大大减小了表面压力对柱浓度的影响,提高了TCCON的反演精度[11]。
TCCON所有站点的经纬度位置和海拔高度如表1所示。TCCON在不同区域的覆盖率差异较大,特别是在南美、非洲和亚洲区域的覆盖率低。站点数据可在https:∥tccondata.org/进行下载。
| 站点名称 | 纬度/°N | 经度/°E | 高度/m | 站点名称 | 纬度/°N | 经度/°E | 高度/m | ||
|---|---|---|---|---|---|---|---|---|---|
| 1 | Eureka | 80.05 | -86.42 | 610 | 19 | Anmyeondo | 36.54 | 126.33 | 30 |
| 2 | Ny-Alesund | 78.9 | 11.9 | 20 | 20 | Tsukuba | 36.05 | 140.12 | 30 |
| 3 | Sodankyla | 67.37 | 26.63 | 188 | 21 | Hefei | 31.90 | 118.67 | 29 |
| 4 | Los Alamos | 35.87 | -106.32 | 2 130 | 22 | Nicosia | 35.14 | 33.38 | 185 |
| 5 | East Trout Lake | 54.35 | -104.99 | 501.8 | 23 | Dryden | 34.96 | -117.88 | 699 |
| 6 | Sites/Bialystok | 53.23 | 23.03 | 180 | 24 | JPL, Pasadena | 34.2 | -118.18 | 390 |
| 7 | Bremen | 53.1 | 8.85 | 27 | 25 | Caltech | 34.14 | -118.13 | 230 |
| 8 | Harwell,Oxfordshire | 51.57 | -1.31 | 123 | 26 | Saga | 33.24 | 130.29 | 7 |
| 9 | Karlsruhe | 49.1 | 8.44 | 116 | 27 | Izana | 28.3 | -16.5 | 2 370 |
| 10 | Paris | 48.85 | 2.36 | 60 | 28 | Burgos | 18.53 | 120.65 | 35 |
| 11 | Orleans | 47.97 | 2.11 | 130 | 29 | Manaus | -3.21 | -60.6 | 50 |
| 12 | Garmisch | 47.48 | 11.06 | 740 | 30 | Ascension Island | -7.92 | -14.33 | 10 |
| 13 | Zugspitze | 47.42 | 10.98 | 2 960 | 31 | Darwin | -12.42 | 130.89 | 30 |
| 14 | Park Falls | 45.95 | -90.27 | 440 | 32 | Reunion Island | -20.9 | 55.49 | 87 |
| 15 | Rikubetsu | 43.46 | 143.77 | 380 | 33 | Wollongong | -34.41 | 150.88 | 30 |
| 16 | Indianapolis | 39.86 | -86 | 270 | 34 | Lauder | -45.04 | 169.68 | 370 |
| 17 | Four Corners | 36.8 | -108.48 | 1 643 | 35 | Arrival Heights | -77.83 | 166.66 | 250 |
| 18 | Lamont | 36.6 | -97.49 | 320 |
2.1.2 NDACC
NDACC是由90多个观测站点组成的全球监测网络,于1991年正式运行。NDACC采用不同的地面观测技术(探空仪,激光雷达,微波辐射计,傅里叶变化红外光谱仪,差分吸收光谱仪以及紫外光谱辐射计)对O3、H2O、CH4等大气成分进行测量,并为星基遥感探测提供验证数据。NDACC可在近红外和中红外波段进行工作,可获取多种气体的柱浓度和廓线浓度。它采用傅里叶红外变换光谱仪(FTIR, Fourier-Transform Infrared Spectrometer)对大气CH4廓线进行探测,探测范围是整个大气层。
NDACC观测网络主要采用Bruker 125系列、Bruker 120系列、JPL MkIV、Bomem DA8、Bomem DA3、Bomem DA2、Eocom等仪器进行测量。其中Bomem DA8傅里叶变换光谱仪系统以0.004 cm-1的分辨率记录太阳透射光谱[13];Bruker IFS 125HR(henceforth referred to as the 125HR)是扫描式迈克尔逊干涉仪,它的光谱分辨率高达0.002 4 cm-1[14];JPL MkIV干涉仪是一种高分辨率傅里叶变换红外光谱仪,用于远程探测大气成分,在光学上,它与大气痕量分子仪(ATMOS, Atmospheric Trace Molecule Spectrometer)类似。NDACC观测网络利用SFIT2对站点测量到的光谱数据进行分析,SFIT2是一种基于最优估计技术的辐射传输和廓线反演算法,通过调整气体廓线和辅助仪器参数,将计算光谱与观测光谱进行拟合[14]。
NDACC可探测CH4的23个站点经纬度、高度、采用的仪器等信息如表2所示,表中的C、P分别表示柱浓度、廓线浓度。FTIR测量数据可以在NDACC官网上获取(http:∥www.ndaccdemo.org/)。NDACC站点主要分布在欧洲和北美洲,非洲、亚洲、大洋洲和南美洲站点比较缺乏。
表2 NDACC可探测CH4浓度的地基站点信息表[15]
Table 2
| 站点名称 | 纬度/°N | 经度/°E | 高度/m | 仪器名称 | 分辨率/cm-1 | 柱(C)/廓线(P) | |
|---|---|---|---|---|---|---|---|
| 1 | Eureka | 80.05 | -86.42 | 610 | Bruker 125HR | 0.002 7 | P |
| 2 | Ny Alesund | 78.92 | 11.93 | 15 | Bruker 120M Bruker 120HR | 0.004 0.002 8 | C / P |
| 3 | Thule | 76.53 | -68.74 | 30~220 | Bruker 120M Bruker 125HR | 0.004 | C |
| 4 | Esrange | 67.9 | 21.1 | 341 | JPL MkIV | 0.006 | C |
| 5 | Kiruna | 67.84 | 20.41 | 419 | Bruker 120HR Bruker 125 | 0.002 | |
| 6 | Fairbanks | 64.82 | -147.87 | 135 | JPL MkIV | 0.006 | C |
| 7 | Harestua | 60.2 | 10.8 | 596 | Bruker 120M Bruker 125M | 0.003 5 | C |
| 8 | St. Petersburg | 59.9 | 29.8 | 20 | Bruker 125HR | 0.002 | P |
| 9 | Bremen | 53.1 | 8.8 | 27 | Bruker 120HR Bruker 125HR | 0.002 8 | C |
| 10 | Zugspitze | 47.42 | 10.98 | 2 964 | Bruker 120HR Bruker 125HR | 0.002 8 | C |
| 11 | Jungfraujoch | 46.55 | 7.98 | 3 580 | home-built FTS Bruker 120HR | P | |
| 12 | Toronto | 43.66 | -79.4 | 174 | Bomem DA8 | 0.004 | C |
| 13 | Rikubetsu | 43.46 | 143.77 | 380 | Bruker 120M Bruker 120/5HR | 0.004 0.003 5 | C |
| 14 | Mt. Barcroft | 37.58 | -118.24 | 3 800 | JPL MkIV | 0.006 | C |
| 15 | Table Mountain Facility | 34.4 | -117.7 | 2 300 | JPL MkIV | 0.006 | C |
| 16 | Izana | 28.3 | -16.48 | 2 367 | Bruker 120M Bruker 125HR | 0.003 5 0.003 5 | C |
| 17 | Mauna Loa | 19.54 | -155.58 | 3 397 | Bruker 120HR Bruker 125HR Bomem DA3 | 0.003 5 0.002 | C |
| 18 | Altzomoni | 19.12 | 98.66 | 3 985 | Bruker 120/5HR | 0.003 5 | P |
| 19 | Paramaribo | 5.75 | -55.2 | 23 | Bruker 120M Bruker 125M | P | |
| 20 | Reunion Island | -20.9 | 55.5 | 85 | Bruker 120M Bruker 125HR | P | |
| 21 | Wollongong | -34.41 | 150.88 | 30 | Bomem DA3 Bomem DA8 Bruker HR 125 | 0.02 0.004 0.004 | C |
| 22 | Lauder | -45.04 | -169.68 | 370 | Bruker 120M Bruker 120HR Bruker 125HR Bruker 120M Bomem DA2 | 0.003 5 0.003 5 0.003 5 0.004 0.02 | C |
| 23 | Arrival Heights | -77.83 | 166.67 | 184 | Bruker 120M Bruker 125HR Bruker 120M Bomem DA2 Eocom | 0.003 5 0.003 5 0.02 0.06 | C |
2.2 空基探测
机载激光红外吸收光谱仪ALIAS(Aircraft Laser Infrared Absorption Spectrometer)是光谱分辨率为0.000 3 cm-1的扫描可调二极管激光光谱仪,用于测量HCl、NO2、HNO3、CH4和N2O[16]。机载大气痕量气体色谱仪ACATS-IV(the Airborne Chromatograph for Atmospheric Trace Species IV)是测量对流层上部和平流层下部CH4的新型大气痕量气体气相色谱仪,它在NASA ER-2高空飞机上全自动工作,飞行时间长达8 h[17]。ALIAS-II是ALIAS的后继仪器,它的垂直分辨率为15 m,可以在气球或飞机上对N2O、CH4、CO、HCl和NO2等气体进行原位测量[18]。俄罗斯科学院普通物理研究所研制的机载可调谐二极管激光光谱仪TDLS(Tunable diode laser spectrometer)可以测量CO2、CH4、H2O等,TDLS测量CH4的灵敏度为20~25 ppbv [13]。
2.3 星基探测
1990年以来,卫星遥感在大气成分监测中的应用越来越广泛。许多国家的航天部门,如美国航空航天局(NASA, National Aeronautics and Space Administration)、欧空局等,陆续发射了许多颗大气探测卫星。通过星基探测获取CH4浓度的方式主要有天底观测、临边/掩星观测。天底观测是指传感器对星下点进行扫描成像,获取目标的后向散射及辐射信息;临边/掩星观测是通过切向观测方式获取大气的垂直分布信息[19]。本文按照卫星不同的观测模式和发射顺序进行介绍。
2.3.1 天底探测传感器发展
根据传感器用于探测CH4的波长范围,将天底模式传感器划分为短波红外天底传感器、中红外和热红外天底传感器。短波红外传感器采用1.3~3 μm范围内的波段来探测CH4,中红外和热红外天底传感器主要采用3~14 μm范围内的波段来探测CH4。
短波红外天底探测器:
表3 短波红外天底探测器主要参数
Table 3
| 卫星/传感器 | 发射时间 | 光谱范围/μm | CH4探测窗口/μm | 光谱分辨率 |
|---|---|---|---|---|
| Terra/MOPITT | 1999.12 | 2.2、2.3、4.7 | 2.17~2.34 | — |
| Envisat/SCIAMACHY | 2002.3 | 0.24~2.38 | 1.63~1.67 | 0.2~1.5 nm |
| Sentinel-5P/TROPOMI | 2017.10 | 0.27~0.775 2.305~2.385 | 2.305~2.385 | 0.25~1 nm |
| GF-5/ GMI | 2018.5 | 0.759~0.769 1.568~2.058 | 1.642~1.658 | 0.27~0.6 cm-1 |
(2)Envisat/SCIAMACHY。SCIAMACHY(Scan-ning Imaging Absorption Spectrometer for Atmospheric Chartography)是由欧空局研制,于2002年3月1日发射升空,自6月20日起开始以天底模式探测。SCIAMACHY的光谱范围是0.24~2.38 μm,光谱分辨率为0.2~1.5 nm,它可以测量大气中O3、O2、SO2、NO、CH4、CO2等大气成分的浓度和分布。采用SCIAMACHY 1.63~1.67 μm探测窗口,基于权重函数修正差分吸收算法(WFM-DOAS, Weighting Function Modified Differential Optical Absorption Spectroscopy)进行CH4柱浓度反演,结果表明反演的柱浓度误差小于6%[22-24]。
(3)Sentinel-5P/TROPOMI。TROPOMI(Tro-pospheric Monitoring Instrument)搭载于Sentinel-5P卫星上,由欧空局于2017年10月13日发射升空。TROPOMI光谱波段覆盖紫外可见光和短波红外,光谱分辨率是0.25~1 nm,空间分辨率为7 km×7 km,它可以测量O3、NO2、SO2、NO、CH4等大气成分。与OMI相比,它扩大了近红外和短波红外的波长范围,并且空间分辨率也有很大提高。TROPOMI反演的CH4柱浓度与地基数据的偏差为14 ppbv (0.8%)[25]。
(4)GF-5/GMI。GF-5/GMI(Greenhouse gas Monitoring Instrument)由中国于2018年5月发射升空。GMI采用空间外差光谱技术进行高光谱分光,是国际上首台基于该体制的星载温室气体遥感设备;它的光谱范围是0.759~0.769 um、1.568~2.058 um,光谱分辨率为0.27~0.6 cm-1,信噪比为300@
中红外与热红外天底探测器:
(1)ADEOS/IMG。ADEOS上的IMG(Interferometric Monitor for Greenhouse gases)由日本国家空间局(NASDA, National Space Development Agency of Japan)于1996年8月发射升空。IMG是第一个高分辨率天底红外仪器,可以探测H2O、CO2、O3、CH4等痕量气体柱浓度和垂直廓线,于1997年6月30日停止工作。IMG的光谱范围是3.3~15 μm,光谱分辨率为0.1~0.4 cm-1。IMG的3.3~4.31 μm可探测CH4浓度[27]。
(2)Aqua/AIRS。Aqua/AIRS(Atmospheric Infrared Sounder)由NASA于2002年5月4日发射升空。它可以探测大气中的CO2、CO、CH4、O3、SO2等痕量气体。AIRS光谱范围是3.8~15.4 μm,光谱分辨率为0.5~2 cm-1。Xiong等[28]基于7.6 μm附近的通道反演了CH4廓线,结果表明AIRS反演的CH4廓线偏差为-1.4%~0.1%。
表4 中红外与热红外天底探测器主要参数
Table 4
| 卫星/传感器 | 发射时间 | 光谱范围/μm | CH4探测窗口/μm | 光谱分辨率/cm-1 |
|---|---|---|---|---|
| ADEOS/IMG | 1996.8 | 3.3~15 | 3.3~4.31 | 0.1~0.4 |
| Aqua/AIRS | 2002.5 | 3.8~15.4 | 7.6及附近波长 | 0.5~2 |
| Aura/TES | 2004.7 | 3.3~15.4 | 7.52~9.09 | 0.06 |
| MetOp/IASI | 2006.10 | 3.62~15.5 | 7.75~8.06 | 0.3~0.5 |
| GOSAT/TANSO-FTS | 2009.1 | 0.757~14.33 | 1.64~1.68 7.63~8.2 | 0.27~0.6 |
| Suomi-NPP/CrIS | 2011.10 | 3.9~15.4 | 7.37~8.2 | 0.625~2.5 |
(3)Aura/TES。TES(Tropospheric Emission Spectrometer)是由喷气推进实验室(JPL, the Jet Propulsion Laboratory)研制,于2004年7月15日发射进入太阳同步轨道。它是地球观测系统Aura飞船上的4个仪器之一,TES可测量从地面到平流层低层中的O3、CH4、H2O、NO、NO2和HNO3的垂直廓线。它的光谱范围是3.3~15.4 μm,TES官方网站获取的一级产品光谱分辨率为0.06 cm-1。TES利用8 μm附近的窗口基于最优化算法反演CH4廓线[29]。TES反演的对流层上层CH4浓度与HIPPO数据的偏差为25.7ppbv[30]。
(4)MetOp-A/IASI。IASI(Infrared Atmospheric Sounding Interferometer)于2006年10月19日搭载在MetOp-A卫星平台上发射升空。IASI光谱范围是3.62~15.5 μm,光谱分辨率为0.3~0.5 cm-1,它可以探测大气中CO、CH4、O3、HNO3、NH3和SO2的分布。它基于最优化算法在7.75~8.06 μm(1 240~1 290 cm-1)波段反演CH4柱浓度,结果表明地面至12 km高度范围内的总反演误差为1%~2.5%[31]。
(5)GOSAT/TANSO-FTS。GOSAT卫星是日本的地球观测卫星,于2009年1月23日发射升空,主要用于监测地球上不同地点的CO2和CH4柱浓度。星上的TANSO-FTS(Thermal and Near-infrared Sensor for carbon Observation Fourier-Transform Spectrometer)是由高分辨率傅里叶变换光谱仪(FTS)、云和气溶胶成像仪(TANSO-CAI)组成。TANSO-FTS光谱范围是0.757~14.33 μm,光谱分辨率为0.27~0.6 cm-1。TANSO-FTS基于最优化算法利用7.63~8.20 μm(1 220~1 310 cm-1)和1.64~1.68 μm(5 937~6 110 cm-1)两个波段反演CH4廓线,结果表明同时利用热红外和短波红外光谱可提高CH4廓线的垂直分辨率[32]。基于TANSO-FTS数据反演的CH4廓线与机载测量数据之间的RMS误差为36.4 ppbv (1.9%),偏差为1.4 ppbv (0.07%)[33]。
(6)Suomi-NPP/CrIS。CrIS(Cross-track Infrared Sounder)搭载于Suomi-NPP卫星上,由美国于2011年10月28日发射升空。它可以探测大气温度、湿度和温室气体廓线。CrIS的光谱范围是3.9~15.4 μm,光谱分辨率为0.625~2.5 cm-1。CrIS用于探测CH4的通道为7.37~8.2 μm(1 220~1 356 cm-1)[34]。
2.3.2 临边/掩星探测传感器发展
表5 临边/掩星探测器主要参数
Table 5
| 卫星/传感器 | 发射时间 | 探测模式 | 光谱范围/μm | CH4探测窗口/μm | 光谱分辨率 |
|---|---|---|---|---|---|
| Nimbus-7/SAMS | 1978.10 | 临边 | 2.7~100 | 7.7及附近波长 | — |
| UARS/HALOE | 1991.9 | 掩星 | 2.45~10.04 | 3.3及附近波长 | 2~4 cm-1 |
| ADEOS-I/ILAS-I | 1996.8 | 掩星 | 0.753~0.784 6.21~11.76 | 7.4/7.8及附近波长 | 0.15 nm 0.129 μm |
| Envisat/SIAMACHY | 2002.3 | 掩星/临边 | 0.24~2.38 | 1.559~1.671 | 0.2~1.5 nm |
| Envisat/MIPAS | 2002.3 | 临边 | 4.1~14.6 | 6.67~8.23 | 0.035 cm-1 |
| ADEOS-II/ILAS-II | 2002.12 | 掩星 | 0.753~0.784 3~12.85 | 7.4 / 7.8 | 0.15 nm 0.003 2~0.129 μm |
| SCISAT-1/ACE-FTS | 2003.8 | 掩星 | 2.2~13.3 | 3.46~3.83 6.84~8.03 | 0.02 cm-1 |
| Aura/TES | 2004.7 | 临边 | 3.3~15.4 | 3.28~3.57 6.45~9.1 | 0.06 cm-1 |
| Aura/HIRDLS | 2004.7 | 临边 | 6~18 | 7.44及附近波长 | — |
| GF-5/AIUS | 2018.5 | 掩星 | 2.4~13.3 | — | 0.02 cm-1 |
(2)UARS/HALOE。UARS/HALOE(Halogen Occultation Experiment)由NASA于1991年9月12日发射升空。HALOE光谱范围是2.45~10.04 μm,光谱分辨率为2~4 cm-1,它的探测范围是80° S~80° N,垂直探测范围是15 km至60~130 km。它可以探测温度、HCl、HF、CH4、H2O、NO、NO2、O3和气溶胶垂直廓线浓度。HALOE采用3.3 μm及附近波长反演CH4廓线,结果表明0.3~50 mbar压强范围内的CH4廓线总误差小于15%,精度优于7%[37]。
(3)ADEOS-I/ILAS-I。ILAS(Improved Limb Atmospheric Spectrometer)是搭载在ADEOS卫星上的传感器,于1996年8月17日发射升空,它采用掩星方式对大气中的O3、HNO3、NO2、N2O、CH4、H2O垂直廓线浓度进行探测。ILAS的光谱范围是0.753~0.784 μm、6.21~11.76 μm,对应的光谱分辨率分别为0.15 nm、0.129 μm。ILAS-I探测CH4的波段与ILAS-II类似,为7.4 μm/ 7.8 μm及附近波长。ILAS-I数据反演的CH4浓度与HALOE CH4浓度相比,除了冬季CH4浓度偏差超过20%(15~30 km),总体上CH4浓度的偏差优于10%(10~50 km)[38]。
(4)Envisat/SIAMACHY。SCIAMACHY于2002年4月3日开始以临边模式进行探测,在4月15日开启辐射冷却器后,大于1.59 μm的光谱信噪比显著提高,结合临边-天底观测模式,SCIAMACHY可以获取对流层中O3、CH4、CO等大气成分柱浓度。掩星模式下,Noel等[39]将剥洋葱算法与WFM-DOAS算法结合,采用1.559~1.671 μm波段进行平流层CH4廓线反演,结果表明平流层的CH4廓线反演精度为5%~10%。
(6)ADEOS-II/ILAS-II。ILAS-II(Improved Limb Atmospheric Spectrometer-II)是由日本环境部研制的掩星仪器,它搭载在ADEOS-II卫星上于2002年12月14日发射升空,但在2003年2月因AD-EOS-II出现故障而退役。ILAS-II是ILAS-I的后继传感器,光谱范围是0.753~0.784 μm、3~12.85 μm,光谱分辨率分别为0.15 nm、0.003 2~0.129 μm,可以对大气中的O3、HNO3、NO2、N2O、CH4、H2O、CFC-11、CFC-12、ClONO2、N2O5等垂直廓线进行测量。ILAS-II采用7.4 μm/7.8 μm及附近波长基于剥洋葱算法反演CH4垂直廓线,结果表明与空基测量数据的偏差为-5%~3%(14~18 km)[42-43]。
(7)SCISAT-1/ACE-FTS。ACE-FTS(Atmospheric Chemistry Experiment Fourier Transform Spectrometer)搭载于SCISAT-1卫星上,于2003年8月12日发射升空。它采用掩星模式对大气痕量气体、温度、压力和气溶胶廓线进行探测。ACE-FTS光谱范围是2.2~13.3 μm,光谱分辨率为0.02 cm-1。ACE-FTS利用3.46~3.83 μm、6.84~8.03 μm(1 245.14~1 463 cm-1、2 610.2~2 888.48 cm-1)波段的数据反演CH4廓线,结果表明对流层至平流层下部CH4廓线精度在10%以内,在平流层中高层至中间层下部(<60 km)精度在25%以内[44]。
(8)Aura/TES。TES主要任务是对对流层的O3和对其含量产生影响的化学物质进行全球三维测量。临边模式下,TES的垂直分辨率为2.3 km,光谱分辨率为0.06 cm-1。它可探测0~35 km范围内的O3、H2O、CO、CH4、NO的廓线浓度。TES于2005年5月停止临边探测,以此来提高仪器和轴承的寿命。它的3.28~3.57 μm、6.45~9.1 μm(1 100~1 550 cm-1、2 800~3 050 cm-1)波段可以探测CH4[45]。
(9)Aura/HIRDLS。Aura/HIRDLS(High Resolution Dynamics Limb Sounder)于2004年7月15日发射升空,它可探测大气中的温度、O3、H2O、CH4、N2O、NO2、HNO3、N2O5等大气成分廓线浓度。HIRDLS具有21个通道,可探测的高度范围是8~80 km。它的光谱范围是6~18 μm,7.44 μm及附近波长可以探测CH4廓线[46]。
(10)GF-5/AIUS。AIUS(Atmospheric Infrared Ultraspectral Soun-der)搭载于我国高分五号卫星上,是我国自主研制的第一颗甚高分辨率红外掩星探测仪,于2018年5月8日发射升空。AIUS主要任务是测量南极上空对流层和平流层上层O3、CH4等痕量气体的化学过程,它的光谱范围为2.4~13.3 um,光谱分辨率为0.02 cm-1,轨道倾角是98.218°,轨道高度为705 km,视场角是1.25 mrad,信噪比为100~350[47]。基于最优化算法(OEM, Optimal Estimation Method)可对O3、H2O、HCl、NO2、N2O、HF等大气成分进行反演[48]。
2.4 不同探测方式和观测模式的比较
地基探测、空基探测和星基探测3种方式均可探测CH4。地基探测和空基探测可为星基探测提供验证资料,星基探测可以获取全球尺度的CH4廓线或柱浓度。3者各有特点,互相补充。
空基探测可以实现实时、实地的探测,也可以在短时间内实现远距离飞行来获取空间气象要素资料。它可以对疑似CH4泄露问题做出快速响应,不需要人员进出高危CH4浓度地区,保证了人员安全。航空器的巡航速度快,这要求仪器具有几秒钟响应的能力,若飞行速度与仪器的响应速度存在差异,这将导致探测数据不能真实反应大气成分状况;空基传感器还需具有重量轻、紧凑和便携等特点[51]。
虽然地基和空基探测可以准确获取局部尺度上的大气CH4浓度,但这两种探测方式的空间覆盖率低[25]。随着科学研究的不断深入,大尺度卫星遥感监测逐渐成为一种必要的研究手段。1979~1983年间,Nimbus-7卫星上的SAMS探测仪首次实现对平流层CH4廓线的探测[42, 52],为Aqua、Terra、Aura、Landsat和SeaWiFS等地球观测系统奠定了基础。随着卫星探测技术的发展和反演算法的改进,卫星传感器反演CH4精度逐渐提高。天底模式下TROPOMI传感器反演的CH4浓度与地面站点数据的偏差为14 ppbv(0.8%)[25];临边/掩星模式下ACE-FTS反演的对流层至平流层下部CH4廓线精度在10%以内[44]。星基探测具有天底、临边/掩星这3种观测模式,临边/掩星传感器的光谱分辨率高于天底传感器,前者的光谱分辨率高达0.02 cm-1(如ACE-FTS、AIUS);天底传感器具有较高的空间分辨率,如TROPOMI的空间分辨率为7 km×7 km;3种观测模式传感器相互补充,共同用于CH4探测。
综上所述,星基探测弥补了地基探测和空基探测空间覆盖率低的问题,地基探测和空基探测又为星基探测精度验证提供数据资料;地基、空基和星基探测3种探测方式互为补充,为大气CH4探测提供数据基础。
3 大气甲烷反演算法研究现状
目前,有许多天底、临边/掩星传感器对大气CH4等温室气体进行监测,并采用不同的方法和波段对CH4柱浓度和廓线浓度进行反演。针对上述内容中提及到的天底、临边/掩星传感器反演CH4柱浓度或廓线浓度的方法,本研究主要对剥洋葱算法、WFM-DOAS算法和最优化算法进行介绍。
3.1 剥洋葱算法
3.2 WFM-DOAS算法
传统的DOAS算法反演痕量气体柱浓度主要包括两个步骤,首先是基于线性最小二乘拟合算法对痕量气体吸收截面、一个低阶多项式和一个Ring参考光谱进行拟合,获取痕量气体斜柱密度;然后利用空气质量因子(AMF)和痕量气体斜柱密度计算垂直柱密度[56]。由于痕量气体在近红外光谱波段的吸收截面与压力、温度和波长有关,传统的DOAS算法不能用于近红外波段的痕量气体柱浓度反演。
其中:
3.3 最优化算法
其中:下标
4 全球大气甲烷分布现状
随着人类文明的发展,CH4的总量一直在增加。42万年前,间冰期温度较低时CH4的最高浓度只有600 ppbv[70]。到工业革命前CH4浓度为700 ppbv[3]。第一次工业革命期间大气CH4浓度增加了62 ppbv,到20世纪90年代大气CH4浓度达到1714 ppbv[4]。不同的研究学者对全球大气CH4浓度的分布特征和变化趋势进行了分析研究。Straume等[28,71-72]先后利用不同的卫星数据对不同时期的全球大气CH4柱浓度和廓线进行了反演,结果表明全球大气CH4浓度存在明显的纬度差异,并且北半球CH4浓度明显高于南半球。刘宝锋等[73]基于GOSAT数据对全球近地面CH4浓度时空分布特征进行了分析,发现亚欧非大陆CH4浓度大于美洲大陆,陆地CH4浓度明显大于海洋。Park等[74]发现北半球的稻田和南半球的常绿阔叶林区存在大量的CH4排放浓度。
4.1 全球甲烷长期变化趋势
图1展示了300 hPa、150 hPa 压强处2010~2020年全球CH4浓度的变化趋势。2010~2020年全球CH4浓度增长了约50 ppbv,年均增长率约为0.29%。300 hPa处,2010年全球CH4浓度为1767 ppbv,到2020年全球CH4浓度增长了55 ppbv;2016~2018年、2019~2020年增长较快。150 hPa处,2010年全球CH4浓度小于1 700 ppbv,到2020年全球CH4浓度增长到1 745.6 ppbv;除2019年外,全球CH4浓度逐年增加。从图中还可以看出,300 hPa处的CH4浓度明显高于150 hPa处的CH4浓度。
图1
图1
2010~2020年全球甲烷浓度变化趋势图
Fig.1
The change trend of global methane concentration from 2010 to 2020
4.2 全球甲烷分布特征
基于AIRS L3的CH4月平均浓度数据,分析了2020年全球CH4在300 hPa与150 hPa压强处的CH4浓度分布特征,如图2所示。
图2
图2
2020年全球大气CH4浓度分布图
Fig.2
Distribution map of global atmospheric CH4 concentration in 2020
5 结 语
大气CH4在地—气系统的辐射平衡过程和平流层化学过程中起着重要作用。目前探测CH4浓度的方式主要有地基探测、空基探测和星基探测。整体而言,星基探测弥补了地基探测和空基探测空间覆盖率低的问题,地基探测和空基探测又为星基探测精度验证提供数据资料;地基、空基和星基探测3种观测方式互为补充,为大气CH4探测提供数据基础。
首先主要对星载大气CH4探测仪的发展进行分析。星载CH4探测器主要在红外波段采用天底或临边模式进行CH4探测。临边/掩星传感器的光谱分辨率高于天底传感器,前者的光谱分辨率高达0.02 cm-1(如ACE-FTS、AIUS);天底传感器具有较高的空间分辨率,如TROPOMI的空间分辨率为7 km×7 km。随着卫星探测技术的发展和反演算法的改进,卫星传感器反演CH4精度逐渐提高。天底模式下TROPOMI传感器反演的CH4浓度与地面站点数据的偏差为14 ppbv(0.8%);临边/掩星模式下ACE-FTS反演的对流层至平流层下部CH4廓线精度在10%以内。
同时,对3种用于反演大气CH4廓线、柱浓度的算法原理、特点及研究进展进行了介绍。剥洋葱算法计算速度快,但它的计算方式容易导致误差累积,致使最低层的反演结果精度较低;权重函数差分吸收算法用于近红外和短波红外波段CH4柱浓度反演,它的优点在于不需要先验信息来约束要反演的柱浓度,有关大气的先验信息用于为无约束线性最小二乘法提供合理的线性化点;最优化算法具有普适性和定量估计误差的优点,它已广泛应用于天底、临边/掩星探测模式传感器大气成分的反演,但最优化算法需要输入先验信息,先验信息的准确与否直接关系到反演结果的精度。
最后基于AIRS L3数据对300 hPa和150 hPa全球CH4浓度分布特征进行分析,发现2010~2020年全球CH4浓度增长了约50 ppbv,年均增长率约为0.29%;全球CH4浓度北高南低,高值区主要分布在大西洋中部、非洲北部、中东地区和中国西部;低值区主要分布在南极地区。
受观测方式和运行轨迹的影响,卫星探测数据空间覆盖范围有限,希望随着科技的发展,可以使更多的传感器实现全球数据的获取。目前有关国产卫星GF5-AIUS传感器甲烷廓线及柱浓度的反演研究比较少,有望以后增进对国产卫星数据的应用。
参考文献
Implication of strongly increased atmospheric methane concentrations for chemistry-climate connections
[J].
Atmospheric CO and CH4 time series and seasonal variations on Reunion Island from ground-based in situ and FTIR (NDACC and TCCON) measurements
[J].
ATMOSPHERIC SCIENCE Ancient ice and the global methane cycle
[J].
Atmospheric methane between 1000 AD and present: Evidence of anthropogenic emissions and climatic variability
[J].
The “source” and “collection” of Greenhouse Gases and its control measure
[J],
王晓峰, 温室气体的“源”与“汇”及其防治对策
[J].
Preliminary study of satellite remote sensing of Greenhouse Gases methane
[D].
The latest developments in the analytical sensing of methane
[J].
Satellite observa-tions of atmospheric methane and their value for quantifying me-thane emissions
[J].
Emissions estimation from satellite retrievals: A review of current capability
[J].
Systematics of atmospheric environment monitoring in China via satellite remote sensing
[J].
The total carbon column observing network
[J].
Calibration of TCCON column-averaged CO2: The first aircraft campaign over European TCCON sites
[J].
Preliminary results of an aircraft system based on near-IR diode lasers for continuous measurements of the concentration of methane, carbon dioxide, water and its isotopes
[J].
A new bruker IFS 125HR FTIR spectrometer for the polar environment atmospheric research laboratory at eureka, nunavut, canada: Measurements and comparison with the existing bomem DA8 spectrometer
[J].
The Network for the Detection of Atmospheric Composition Change (NDACC): History, status and perspectives
[J].
Aircraft(ER-2) Laser Infrared-Absorption Spectrometer(ALIAS) for In-Situ Stra to Spheric Measurements of HCL,N2O,CH4,NO2,and HNO3
[J].
Airborne gas chromatograph for in situ measurements of long-lived species in the upper troposphere and lower stratosphere
[J].
Airborne Laser Infrared Absorption Spectrometer (ALIAS-II) for in situ atmospheric measurements of N2O,CH4,CO,HCl,and NO2 from balloon or remotely piloted aircraft platforms
[J].
Retrieval and monitoring of atmospheric trace gas concentrations in nadir and limb geometry using the space-borne SCIAMACHY instrument
[J].
Radiative transfer modeling for the EOS Terra satellite Measurement of Pollution in the Troposphere (MOPITT) instrument
[J].
Effects of a spectral surface reflectance on measurements of backscattered solar radiation: Application to the MOPITT methane retrieval
[J].
Tropical methane emissions: A revised view from SCIAMACHY onboard ENVISAT
[J].
Retrieval of CH4
Atmospheric carbon gases retrieved from SCIAMACHY by WFM-DOAS: version 0.5 CO and CH4 and impact of calibration improvements on CO2 retrieval
[J].
A scientific algorithm to simultaneously retrieve carbon monoxide and methane from TROPOMI onboard Sentinel-5 Precursor
[J].
Greenhouse gases Monitoring Instrument(GMI) on GF-5 satellite (invited)
[J].
Interferometric Monitor for Greenhouse Gases(IMG)
[C]∥ Scholl M S, eds.
Characterization and validation of methane products from the Atmospheric Infrared Sounder (AIRS)
[J].
Profiles of CH4, HDO, H2O, and N2O with improved lower tropospheric vertical resolution from Aura TES radiances
[J].
Validation of TES methane with HIPPO aircraft observations: implications for inverse modeling of methane sources
[J].
Characterization of methane retrievals from the IASI space-borne sounder
[J].
Methane vertical profiles retrieval from IASI/METOP and TANSO-FTS/GOSAT data
[C]∥ Matvienko G G, Romanovskii O A, eds.
Characterization and validation of CO2 and CH4 products derived from the GOSAT thermal infrared band
[C]∥ Shimoda H, Xiong X, eds.
Reprocessing of Suomi NPP Cris sensor data records and impacts on padiometric and spectral long-term accuracy and stability
[C]∥ IEEE eds.
The stratospheric and mesospheric sounder on Nimbus-7
[J].
J J. Retrieval of temperature and composition from Nimbus-7 Sams measurements
[J].
Validation of halogen occultation experiment CH4 measurements from the UARS
[J].
Validation and data characteristics of nitrous oxide and methane profiles observed by the Improved Limb Atmospheric Spectrometer (ILAS) and processed with the Version 5.20 algorithm
[J].
Stratospheric CH4 and CO2 profiles derived from SCIAMACHY solar occultation measurements
[J].
Retrieval of temperature, H2O, O3, HNO3, CH4, N2O, ClONO2 and ClO from MIPAS reduced resolution nominal mode limb emission measurements
[J].
Validation of MIPAS O3,NO2,H2O and CH4 Profiles(V4.61) with Colloca-ted Measurements of HALOE and SAGE II
[C]∥Danesy,eds.
Validation of the improved limb atmospheric spectrometer-II (ILAS-II) version 1.4 nitrous oxide and methane profiles
[J].
Improved Limb Atmospheric Spectrometer (ILAS) data retrieval algorithm for Version 5.20 gas profile products
[J].
Validation of ACE-FTS v2.2 methane profiles from the upper troposphere to the lower mesosphere
[J].
Tropospheric emission spectrometer for the Earth Observing System's Aura Satellite
[J].
Infrared filters and coatings for the high resolution dynamics limb sounder (6~18 μm)
[J].
Monitoring trace gases over the antarctic using atmospheric infrared ultraspectral sounder onboard GaoFen-5: algorithm description and first retrieval results of O3, H2O,and HCl
[J].
Assessment of retrieved N₂O, NO₂, and HF profiles from the atmospheric infrared ultraspectral sounder based on simulated spectra
[J].
Analysis of the annual variation trend of atmospheric methane over China
[C]∥ Gao W, Jackson T J, Wang J, et al.
Ground-based Station Network in Arctic and Subarctic Eurasia: An overview
[J].
A cryogenically operated laser diode spectrometer for airborne measurement of stratospheric trace gases
[J].
Observations of CH4 and N2O by the NIMBUS-7 SAMS-a Comparison wutg insitu data and TWO-Dunebsuibak Numerical-MODEL cakcykatuibs
[J].
Inversion algorithms for recovering minor species densities from limb scatter measurements at UV-visible wavelengths
[J].
Equivalence method for retrieving stratospheric constituent profiles from infrared solar occultation data
[J].
Evaluation of the impact of ionospheric asymmetry on GNSS radio occultation inversion products using NeQuick and IRI
[C]∥
A near-infrared optimized DOAS method for the fast global retrieval of atmospheric CH4, CO, CO2, H2O, and N2O total column amounts from SCIAMACHY Envisat-1 nadir radiances
[J].
WFM-DOAS algorithm research and application on the greenhouse gases column density retrieval
[D].
基于WFM-DOAS温室气体柱浓度反演方法及应用研究
[D].
Atmospheric methane and carbon dioxide from SCIAMACHY satellite data: initial comparison with chemistry and transport models
[J].
Carbon monoxide, methane and carbon dioxide over China retrieved from SCIAMACHY/ENVISAT by WFM-DOAS
[C]∥
Three years of greenhouse gas column-averaged dry air mole fractions retrieved from satellite-Part 2:Methane
[J].
Evaluation of single-footprint AIRS CH4 profile retrieval uncertainties using aircraft profile measurements
[J].
Atmospheric temperature, moisture and trace gas retrieval algorithm for the NPOESS/NPP Cross-track Infrared Microwave Sounder Suite (CrIMSS)
[C]∥ Shen S S, Lewis P E, eds.
Retrievals for the atmospheric chemistry experiment Fourier-transform spectrometer
[J].
ILAS (improved limb atmospheric spectrometer) ADEOA data retrieval algorithms
[C]∥ Susskind J, Singhrovy V, Tanaka S, eds.
The joint methane profiles retrieval approach from GOSAT TIR and SWIR spectra
[C]∥ Matvienko G G, Romanovskii O A,eds.
Retrieval algorithm for CO2 and CH4 column abundances from short-wavelength infrared spectral observations by the Greenhouse gases observing satellite
[J].
A fast radiative transfer model for simulating high-resolution absorption bands
[J].
An Improved constraint method in optimal estimation of CH4 from gosat swir observations
[C]∥ IEEE eds.
Atmospheric methane and nitrous oxide of the late Pleistocene from Antarctic ice cores
[J].
The global variation of CH4 and CO as seen by SCIAMACHY
[C]∥ Burrows J P, Eichmann K U.
Toward blobal mapping of methane with TROPOMI: First results and intersatellite comparison to GOSAT
[J].
Analysis of spatio-temporal distribution of methane concentrations near the ground gased on GOSAT satellite
[J].
基于GOSAT卫星的全球近地面甲烷浓度时空分布特征分析
[J].
The extraction of methane sources over the global area using satellite data
[C]∥
Spatial and temporal distribution of atmospheric methane in middle-low reaches of Yangtze River based on satellite observations
[J].
长江中下游六省大气甲烷柱浓度时空分布
[J].
The variation characteristics of methane concentration in the Three Gorges reservoir area were observed by AIRS
[J].
AIRS遥感观察三峡库区甲烷浓度变化特征
[J].
Seven years' observation of mid-upper tropospheric methane from atmospheric infrared sounder
[J].
The Vertical variations of Atmospheric Methane(CH4) concentrations over selected cities in Iraq based on AIRS data
[J].
Distribution of methane concentration in the world and East Asia
[C]∥
/
| 〈 |
|
〉 |

