环境科学  2019, Vol. 40 Issue (10): 4412-4422   PDF    
引用本文
王闰芳, 马晓丹, 赵天良, 王红磊, 丁涵, 郑小波. 利用地面观测资料和AIRS卫星资料评估MACC再分析臭氧数据中国地区的适用性[J]. 环境科学, 2019, 40(10): 4412-4422.
WANG Run-fang, MA Xiao-dan, ZHAO Tian-liang, WANG Hong-lei, DING Han, ZHENG Xiao-bo. Evaluation of MACC Reanalysis Ozone Data over China Using Ground-based and AIRS Satellite Observations[J]. Environmental Science, 2019, 40(10): 4412-4422.
利用地面观测资料和AIRS卫星资料评估MACC再分析臭氧数据中国地区的适用性
王闰芳1, 马晓丹1, 赵天良1, 王红磊1, 丁涵1, 郑小波2     
1. 南京信息工程大学气象灾害预报预警与评估协同创新中心, 中国气象局气溶胶-云-降水重点开放实验室, 南京 210044;
2. 贵州省山地环境气候研究所, 贵阳 550000
收稿日期: 2019-04-02; 修订日期: 2019-05-05
基金项目: 国家重点研发计划项目（2016YFC0203304）；国家自然科学基金项目（41830965，91744209）；江苏省研究生科研与实践创新计划项目（KYCX18_1027）
作者简介: 王闰芳(1994~), 女, 硕士研究生, 主要研究方向为大气环境变化, E-mail:wangrunfang1995@163.com
通信作者: 赵天良, E-mail:tlzhao@nuist.edu.cn
摘要: 用10 a（2003~2012年）的MACC（monitoring atmospheric composition and climate）再分析臭氧资料与同期中国大陆6个地面观测站的臭氧数据，以及AIRS（atmospheric infrared sounder）卫星资料反演的大气臭氧数据进行对比分析.结果表明，MACC臭氧柱浓度与所有地面站臭氧柱浓度相对偏差基本控制在17%之内，相关系数在0.79~0.97之间，MACC臭氧柱浓度与地面站臭氧柱浓度具有很好的一致性.对于多年平均臭氧柱浓度空间分布，MACC再分析与AIRS卫星反演臭氧柱浓度的相对偏差在-3%~5%之间，MACC臭氧柱浓度相比AIRS柱浓度在青藏高原、南部沿海地区偏高，在东北地区偏低.各季节空间变化，MACC臭氧柱浓度与AIRS基本一致.在瓦里关站，MACC近地面臭氧与地面观测臭氧浓度月平均变化趋势存在一致性.MACC再分析近地面臭氧数据可以反映春季、夏季和秋季地面臭氧浓度的变化趋势，但冬季MACC近地层臭氧资料出现了较大的偏差.
关键词: 臭氧      MACC产品      AIRS卫星      中国     
Evaluation of MACC Reanalysis Ozone Data over China Using Ground-based and AIRS Satellite Observations
WANG Run-fang1 , MA Xiao-dan1 , ZHAO Tian-liang1 , WANG Hong-lei1 , DING Han1 , ZHENG Xiao-bo2     
1. Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing University of Information Science and Technology, Nanjing 210044, China;
2. Guizhou Institute of Mountainous Environment and Climate, Guiyang 550000, China
Abstract: A comparative analysis was conducted using total ozone products derived from monitoring atmospheric composition and climate (MACC) reanalysis data validated with ozone data from the atmospheric infrared sounder (AIRS) satellite and ground-based ozone measurements. The results indicate that the relative deviation of total ozone from the MACC reanalysis data and the ground-based ozone total data is controlled within 17%, and all of the correlation coefficients were between 0.79 and 0.97. The total ozone values from the MACC reanalysis data showed good consistency with the ground-based ozone measurements. With respect to the spatial distribution of multi-year averages, the relative deviation of total ozone values in the MACC reanalysis data and the AIRS satellite data was between -3% and 5%. The total ozone values in the MACC reanalysis data were higher than those from AIRS measurements for the Qinghai-Tibet Plateau and the coastal areas of South China, and were lower for northeast China. Furthermore, the seasonal variations in total ozone values in the MACC reanalysis data were consistent with AIRS data. At Mt. Waliguan station, the monthly averaged trends for near-surface ozone in the MACC reanalysis data were also consistent with surface ozone concentrations; the MACC reanalysis data reflect the observed trends for surface-based ozone measurements in spring, summer, and autumn, but show a large deviation in winter.
Key words: ozone      MACC reanalysis      AIRS satellite      China     

臭氧是大气中的重要痕量气体, 大部分集中在10~30 km的大气平流层, 其中20~25 km处浓度最高[1], 能够大量吸收紫外线有效地保护地球生物, 其在全球的辐射平衡过程和大气加热率计算中都扮演了重要的角色, 并且是主要的温室气体[2].在对流层中, 臭氧具有活泼的光化学性质, 能够发生多种光化学反应, 特别是在近地层, 臭氧是光化学烟雾的重要组成部分[3], 过量的臭氧对人体健康[4]和植物生长[5]都不利.近些年来, 人类排放导致近地面大气臭氧浓度不断升高[6~8], 光化学烟雾频繁发生[9, 10].当前, 我国臭氧污染形势日益严峻, 尤其是夏日高温时臭氧污染最为严重[11~13], 且大有赶超PM2.5污染的势头.

欧洲中期天气预报中心ECMWF (European Center of Medium-range Weather Forecasts)的MACC再分析数据集提供了包括臭氧(O3)在内的多种大气成分的分析场和气溶胶[14]. Katragkou等[15]使用MACC再分析资料对2003~2010年欧洲地面臭氧进行了研究, 结果表明近地表臭氧在冬季呈现负偏差, 在暖季呈现正偏差. Gaudel等[16]使用MOZAIC-IAGOS数据评估了2003~2010年欧洲地区MACC数据的适用性, 并探究了臭氧和一氧化碳的分布情况. Georgoulias[17]等使用CALIOP/CALIPSO卫星观测对欧洲、北非和中东地区的MACC再分析沙尘产品进行三维评估.目前为止用MACC再分析臭氧产品的研究基本都是针对欧洲区域.然而, 由于中国区域缺乏较长序列臭氧地面观测数据, 对于MACC臭氧产品在中国地区评估仍然缺失, 因此, 进行这项工作很有必要.本文将10a(2003~2012年)中国区域的MACC再分析臭氧数据与同期中国6个地面站的观测臭氧数据, 以及AIRS卫星反演的臭氧数据进行对比, 比较分析大气臭氧柱浓度的时空变化, 及近地面臭氧浓度的时间变化, 探讨MACC资料的在中国区域的适用性, 以期为中国区域地面臭氧分布和变化研究中, 合理选用MACC再分析臭氧资料提供借鉴和参考.

1 材料与方法 1.1 MACC再分析

MACC再分析数据集是由CAMS (the Copernicus Atmosphere Monitoring Service)提供的全球大气成分数据集. MACC再分析数据集是通过将卫星数据同化到一个基于ECMWF综合预报系统(IFS)的全球模式数据同化系统里建立的.其中, IFS已扩展到包括化学反应气体[18, 19]、气溶胶[20]和温室气体[21], 以便其具有的四维变分数据同化(4D-Var)系统可以用来同化全球尺度上大气成分的卫星资料.将化学传输模型MOZART与IFS耦合, 为反应气体提供初始场、化学生成和损失率[18]. MACC再分析数据同化系统包括气象场, 大气成分的分析场[臭氧(O3)、一氧化碳(CO)、氮氧化物(NOx)、二氧化硫(SO2)、甲醛(HCHO)和气溶胶[14]]以及一个具有12 h分析窗的4维变分分析(4D-Var). MACC数据集提供全球范围数据, 水平分辨率为~80 km(T159光谱分辨率)和垂直方向分为60个混合西格玛压力层, 并延伸至0.1 hPa.

表 1为MACC再分析资料同化的臭氧卫星资料, MIPAS和MLS分别在光谱的中红外和微波部分进行测量, 是MACC再分析中唯一不受光照条件影响的臭氧数据.连同GOME的O3廓线数据, 它们也是具有较高垂直分辨率的臭氧廓线同化数据, 这些卫星数据在MACC再分析同化系统中已被证明对于获得真实的臭氧垂直分布至关重要[22, 23]. MACC再分析同化系统利用OASIS4耦合器将化学输送模式MOZART CTM耦合到ECMWF-IFS之中, 并将GOME臭氧廓线数据、MIPAS臭氧廓线数据、MLS臭氧廓线数据、SBUV/2臭氧廓线数据、SCIAMACHY柱浓度臭氧数据及OMI臭氧柱浓度数据同化到MACC再分析同化系统之中, 同化实验能够显著改善臭氧积分场的模拟效果, 特别是能够改善控制实验模拟的南极臭氧洞不够“深”的现象; 但在对MLS臭氧廓线数据使用变分偏差校正方案时导致对流层臭氧增加, 但并没有影响臭氧柱浓度场. 2008年1月1日MLS的偏差校正关闭后, 这个数据漂移问题不再出现. MACC再分析数据集目前发布了2003~2012年期间的臭氧(O3)、一氧化碳(CO)、氮氧化物(NOx)、二氧化硫(SO2)、甲醛(HCHO)、气溶胶和气象场数据, 以及1979~2014年期间的二氧化碳(CO2)、甲烷(CH4)和一氧化二氮(N2O)的大气成分分析场.本文主要是对中国大陆地区2003~2012年期间MACC再分析数据集的臭氧资料进行评估.本文选取的臭氧资料分辨率为0.5°×0.5°.

表 1 MACC再分析资料同化的臭氧卫星资料 Table 1 Ozone satellite data included in the MACC reanalysis

1.2 AIRS卫星数据

大气红外探测仪(AIRS)卫星资料提供了全球温度、水汽、臭氧等分布及变化.目前, AIRS发布的数据包含V5和V6版本, 由于基本的反演方法变化较少, V6版本只做了稍许改进, 如不确定性减小[24]. AIRS的level-2和level-3产品中都有臭氧数据, 前者的分辨率约为45 km. level-3数据是level-2的数据处理后生成, 空间分辨率为1°×1°. Bian等[25]的研究表明AIRS用(41个)通道测量臭氧柱浓度, 不确定度估计在热带地区为5%, 在极地地区为5%~40%; Rajab等[26]使用AIRS反演的资料对马来西亚半岛臭氧特征进行了分析.潘亮等[27]的研究表明AIRS在中国地区有很大的应用潜力.从表 1可看出, MACC未使用AIRS资料进行同化, 因此, 可以使用AIRS资料对MACC再分析资料进行评估.本文选择的是NASA GES DISC提供的2003~2012年V6版本level-3月臭氧产品, 利用NCL的插值函数interp2将AIRS数据插值到0.5°×0.5°的格点上.

1.3 地面臭氧观测

根据中国地理特征及数据获取的便利性, 选取了龙凤山站、香河站、临安站、瓦里关站、拉萨站以及昆明站这6个地面站观测的臭氧柱浓度数据同MACC臭氧柱浓度进行对比, 站点分布具体信息如图 1和表 2所示.龙凤山站位于黑龙江省五常市, 可代表东北区域; 香河地处华北平原北部廊坊市, 可作为华北平原代表站; 临安站位于浙江杭州, 可代表长三角区域; 青海的瓦里关和西藏的拉萨可代表青藏高原; 昆明地处云贵高原中部, 可代表中国西南.其中香河站和昆明站的地基臭氧柱浓度数据分别取自于Dobson#075和Dobson#003光谱仪的观测数据.龙凤山站、瓦里关站和临安站分别是来自于Brewer#054、Brewer#076和Brewer#077光谱仪的观测臭氧柱浓度.这些仪器测得的大气臭氧柱浓度单位是DU(Dobson unit), 为垂直单位气柱中的臭氧分子总数换算到标准温压状况下的厚度, 也称为大气厘米, 1DU=10-3 atm-cm=2.69×1016臭氧分子.另外, 还选取2003~2012年瓦里关站地面臭氧浓度数据同MACC再分析近地层臭氧资料进行比较.

图 1 臭氧地面观测站点分布示意 Fig. 1 Distribution of ground-based monitoring stations for total ozone

表 2 臭氧地面观测站点的信息 Table 2 Summary of ground-based monitoring stations for total ozone

2 结果与讨论 2.1 大气臭氧柱浓度时间变化

用这6个地面观测站的逐日臭氧柱浓度数据, 与同时期的MACC再分析臭氧柱浓度进行对比, 图 2为两种逐日臭氧柱浓度随时间的变化曲线, 图 3为各个地面站两种逐日臭氧柱浓度的相对偏差(即MACC柱浓度减去地面观测柱浓度, 除以地面观测柱浓度).各站点的臭氧柱浓度平均值如表 3所示.由图 2可看出臭氧柱浓度呈周期性变化, 春季柱浓度较高, 秋季柱浓度低.香河站地面观测资料的连续性最好, 观测臭氧柱浓度变化趋势与MACC的变化相当一致, 地面观测的臭氧柱浓度平均值为343.92DU, MACC的臭氧柱浓度平均值为339.08DU, 相对偏差基本控制在15%之内.在龙凤山站, 观测臭氧柱浓度的变化趋势与MACC臭氧柱浓度变化一致性较好, 地面观测的臭氧柱浓度平均值为362.34DU, MACC的臭氧柱浓度平均值为366.84DU, 相对偏差小于17%.在昆明站, 观测臭氧柱浓度的变化趋势与MACC臭氧柱浓度变化一致, 地面观测臭氧柱浓度平均值为264.51DU, MACC臭氧柱浓度平均值为267.86DU, 整体相对偏差基本在8%之内.在瓦里关站, MACC臭氧柱浓度整体相比于观测资料偏大, 但相对偏差基本控制在12%, 地面观测的臭氧柱浓度平均值为288.75DU, MACC的臭氧柱浓度平均值为297.48DU.临安站在2004年至2006年8月期间MACC臭氧柱浓度相比于地面观测的偏大, 但相对偏差基本控制在17%之内.在拉萨地面站, 2008~2009年期间MACC臭氧柱浓度相比于地面观测偏小, 2011年后MACC偏大, 但偏差都在10%之内.总的来说, MACC臭氧柱浓度与地面站观测的臭氧柱浓度具有很好的一致性.

图 2 各地面站观测的逐日臭氧柱浓度和MACC再分析臭氧柱浓度比较 Fig. 2 Daily variations in total column ozone concentrations retrieved from the MACC reanalysis and ground-based observations

图 3 各地面站观测的逐日臭氧柱浓度和MACC再分析臭氧柱浓度的相对偏差 Fig. 3 Relative differences in total column ozone concentrations retrieved from the MACC reanalysis and ground-based observations

表 3 MACC再分析臭氧柱浓度与地面观测站臭氧柱浓度的平均值及其间的线性回归系数 Table 3 Mean ozone concentrations and coefficients of linear regression for total ozone from the MACC reanalysis and ground-based observations

图 4为各个地面站观测的臭氧柱浓度和相应MACC臭氧柱浓度的线性回归, 回归系数如表 3所示.各站点MACC臭氧柱浓度与地面观测柱浓度的线性相关性都很高, 相关系数在0.79~0.97之间.不同站点MACC臭氧柱浓度与观测柱浓度的相关系数都大于0.7, 特别是香河、昆明、瓦里关和龙凤山的相关性更为明显, 达到0.9以上.其中龙凤山站整体的相关系数最高, 达到0.97, 斜率为0.97.相关系数最小的是拉萨站, 为0.79, 斜率为0.93.所有站点均通过了置信区间99.9%的显著性检验, 表明MACC再分析臭氧柱浓度数据可靠性较高.

图 4 各地面站点逐日臭氧柱浓度和MACC再分析臭氧柱浓度相关性分析 Fig. 4 Correlation analyses of daily variations in total ozone from the MACC reanalysis and ground-based observations

2.2 大气臭氧柱浓度空间变化

从图 5中可知, 在空间尺度上, MACC和AIRS的臭氧柱浓度分布基本一致, 呈现较明显的纬向分布, 从低纬度向高纬度逐渐增高, 在青藏高原出现明显的低值中心, 这与周秀骥等[34]的研究结果一致.中国区域臭氧柱浓度值基本在260DU以上.导致臭氧柱浓度空间分布变化的主要原因是太阳辐射和大气环流影响.臭氧柱浓度在青藏高原出现最低值, 主要是由于高海拔引起的[35, 36].

图 5 2003~2012年MACC再分析和AIRS卫星的大气臭氧柱浓度空间分布 Fig. 5 Spatial distribution of total ozone derived from the MACC reanalysis and ground-based observations during 2003-2012

图 6为2003~2012年AIRS臭氧柱浓度与MACC臭氧柱浓度的相对偏差(即MACC柱浓度减去AIRS柱浓度, 除以AIRS柱浓度).从中可知, 两种数据空间分布基本一致, MACC与AIRS的臭氧柱浓度的相对偏差在-3%~5%之间, 其中, 在东北地区, MACC的臭氧柱浓度小于AIRS的臭氧柱浓度, 相对偏差大约在1%~3%之间.在南部沿海, MACC的臭氧柱浓度大于AIRS的臭氧柱浓度, 相对偏差大约在1%~2%之间.另外, 在青藏高原, MACC的臭氧柱浓度大于AIRS的臭氧柱浓度, 相对偏差最大可达5%.这可能是青藏高原不适当的地表反照率输入到了MACC资料同化的卫星反演算法中引起, 而具体原因有待进一步探究.总体看来, MACC臭氧柱浓度与AIRS的臭氧柱浓度的空间分布基本一致, 偏差基本控制在6%之内, 说明MACC再分析臭氧柱浓度能够揭示出臭氧柱浓度的空间差异.

图 6 2003~2012年MACC和AIRS卫星的大气臭氧柱浓度空间相对差异 Fig. 6 Difference values for total ozone between MACC and AIRS datasets during 2003-2012

2.3 大气臭氧柱浓度季节分布

为探究MACC臭氧柱浓度资料在各季节的适用性, 将2003~2012年春(3~5月)、夏(6~8月)、秋(9~11月)和冬(12月~次年2月)四季的MACC臭氧柱浓度和AIRS臭氧柱浓度进行了对比, 见图 7和8.臭氧柱浓度存在明显的季节性变化, 而且高纬度地区臭氧柱浓度季节变化更为明显.其主要原因是太阳辐射的季节变化, 尤其是太阳紫外辐射的季节变化引起, 它是大气中产生臭氧光解作用的关键性因素.高纬度地区太阳辐射的季节变化较大, 因此, 由O2光解产生O3的变化也较大[37].对于各季节, MACC臭氧柱浓度和AIRS臭氧柱浓度空间分布基本一致.结合图 9也可以看出, MACC与AIRS臭氧柱浓度在一些地区存在差异, 春季, 在青藏高原及南部沿海地区MACC臭氧柱浓度相比于AIRS的偏高, 相对偏差在2%~6%之间; 而在东北地区略偏低, 相对偏差在2%~3%之间; 夏季, 在青藏高原MACC臭氧柱浓度比AIRS的高, 相对偏差为2%~5%, 而在东北地区MACC臭氧柱浓度偏低, 相对偏差最高达到6%;秋季, 在青藏高原及其周边地区MACC臭氧柱浓度偏高, 最大相对偏差为6%, 相对偏差高值中心出现在塔里木盆地地区; 冬季, 在青藏高原及其周边地区和南部沿海地区MACC臭氧柱浓度偏高, 相对偏差大约为2%~5%.虽然MACC和AIRS的臭氧数据存在一定的误差, 但是最大相对偏差仅为6%, 说明MACC臭氧柱浓度可以揭示出臭氧柱浓度季节变化特征.

图 7 2003~2012年MACC再分析的大气臭氧柱浓度各季节空间分布 Fig. 7 Spatial distribution of total ozone derived from the MACC reanalysis during 2003-2012

图 8 2003~2012年AIRS卫星的大气臭氧柱浓度(DU)各季节空间分布 Fig. 8 Spatial distribution of total ozone derived from ground-based observations during 2003-2012

图 9 2003~2012年MACC和AIRS的大气臭氧柱浓度(DU)空间相对差异 Fig. 9 Difference values for total ozone between the MACC reanalysis and AIRS datasets during 2003-2012

2.4 地面臭氧浓度分析

图 10为2003~2012年瓦里关站地面臭氧浓度与MACC近地面臭氧数据对比分析, 从中可知, 瓦里关站近地面臭氧浓度具有准周期季节波动. MACC近地面臭氧浓度可以大致反映出这种变化趋势, 但部分时间段内MACC臭氧数据偏差较大, 其中冬季整体偏大, 相对偏差最大超过50%.表明MACC近地面臭氧可以反映出近地面臭氧浓度月际变化的大致趋势.

图 10 MACC近地层月平均臭氧浓度和地面观测的月平均地面臭氧浓度对比分析 Fig. 10 Comparative analysis of monthly averaged ozone concentrations from the MACC reanalysis and ground-based observations

图 11为瓦里关站地面观测臭氧浓度与MACC近地面臭氧浓度的季节变化对比.春季, MACC近地面臭氧浓度与地面观测值具有相同的变化趋势, 相关系数达到0.75, 但具体数值存在一些差异, 相对偏差控制在±20%之内.夏季, MACC近地面臭氧浓度与地面观测的一致性较好, 但MACC比地面观测的值整体偏高, 相对偏差在-10%~22%.秋季, MACC近地面臭氧浓度可以反映出地面臭氧浓度的变化趋势, 但是具体数值差异偏大, 最大差异超过30%.冬季, MACC近地面臭氧浓度与地面观测的一致性较差, MACC臭氧值整体比地面观测臭氧值大, 相对偏差最大超过50%.根据以上分析, MACC近地面臭氧浓度基本可以反映地面臭氧浓度季节变化特征, 其中, 在春夏秋三季较好地反映了地面的变化趋势, 但冬季MACC近地面臭氧浓度与地面资料相比出现了较大的偏差. 图 12为瓦里关站观测的年平均地面臭氧浓度和同期MACC近地面臭氧浓度年平均值对比分析.在2004~2010年MACC与地面观测的臭氧浓度值变化趋势一致, MACC整体偏高, 其余年份MACC偏低, 但总的来说, MACC近地面臭氧浓度与地面观测臭氧浓度数值上差异较小, 相对偏差控制在12%之内.表明MACC近地面臭氧数据可以较好地反映出地面臭氧浓度年际变化.

图 11 MACC近地层臭氧浓度和地面观测臭氧浓度的季节变化对比 Fig. 11 Seasonal comparative analysis of ozone concentrations from the MACC reanalysis and ground-based observations

图 12 MACC近地层年平均臭氧浓度和地面观测的年平均地面臭氧浓度对比分析 Fig. 12 Comparative analysis of annually averaged ozone concentrations from the MACC reanalysis and ground-based observations

综上所述, 从瓦里关站的比较看, MACC反演局地短时间尺度的地面臭氧浓度变化尚存在局限性.出现这种现象可能是瓦里关站地处青藏高原东北坡, 地形复杂, 地表状态变化较为频繁[38], 不准确的地表反照率输入卫星遥感反演算法中引起卫星资料不准确; 此外, 选取的MACC资料分辨率较粗(约55 km), 用“面上”的平均数据与地面观测点进行比较难免会出现一些匹配差距.因此, 在使用冬季的MACC近地面臭氧浓度资料时必须非常谨慎.

3 结论

(1) MACC臭氧柱浓度与6个地面站观测的臭氧柱浓度相对偏差基本控制在17%之内, 相关系数均在0.79~0.79之间, 两种臭氧柱浓度资料具有很好的一致性.

(2) MACC臭氧柱浓度与AIRS卫星资料臭氧柱浓度的相对偏差在-3%~5%之间, 主要出现在青藏高原、南部沿海和东北地区.总体看来, MACC臭氧柱浓度与AIRS臭氧柱浓度值的空间分布基本一致.

(3) 在瓦里关站, MACC近地层臭氧数据可以反映出地面臭氧浓度月际变化的大致变化趋势.从季节看, MACC再分析月平均近地层臭氧数据基本可以反映地面臭氧浓度的季节变化特征, 但在冬季出现了较大的偏差.此外, MACC近地面臭氧可以较好地反映出地面臭氧浓度年际变化.

(4) MACC对于研究本世纪以来中国地面臭氧柱浓度时空变化有较好的参考意义.从瓦里关站比较的情况看, 由于我国地面臭氧观测站点相对较少、分布不均, 获取大范围、长序列的臭氧数据比较困难, 在国内大部分地区使用MACC资料进行地面月或者年尺度臭氧分析研究不乏为一种替代方法.但是, MACC反演局地短时间尺度的地面臭氧浓度变化尚存在局限性.今后还需要应用更多低海拔地区地面测站的资料与MACC进行比较验证后, 才能对MACC资料在反演地面臭氧的适应性进行较为全面的评价.

参考文献
[1] 唐孝炎, 张远航, 邵敏. 大气环境化学[M]. 第二版. 北京: 高等教育出版社, 2006.
[2] IPCC. Climate Change 2007: The physical science basis. contribution of working group Ⅰ to the fourth assessment report of the intergovernmental panel on climate change[R]. Cambridge: Cambridge University Press, 2007.
[3] Berntsen T, Isaksen I S A, Wang W C, et al. Impacts of increased anthropogenic emissions in Asia on tropospheric ozone and climate[J]. Tellus B:Chemical and Physical Meteorology, 1996, 48(1): 13-32. DOI:10.3402/tellusb.v48i1.15662
[4] 孔琴心, 刘广仁, 李桂忱. 近地面臭氧浓度变化及其对人体健康的可能影响[J]. 气候与环境研究, 1999, 4(1): 61-66.
Kong Q X, Liu G R, Li G C. Surface ozone concentration variation and possible influences on human health[J]. Climatic and Environmental Research, 1999, 4(1): 61-66. DOI:10.3878/j.issn.1006-9585.1999.01.14
[5] 王春乙, 关福来. O3浓度变化对我国主要农作物产量的可能影响[J]. 应用气象学报, 1995, 6(S1): 69-74.
Wang C Y, Guan F L. The potential effects of O3 concentration change on yield of main crops in China[J]. Quarterly Journal of Applied Meteorology, 1995, 6(S1): 69-74.
[6] Reinsel G C, Tiao G C, Wuebbles D J, et al. Seasonal trend analysis of published ground-based and TOMS total ozone data through 1991[J]. Journal of Geophysical Research:Atmospheres, 1994, 99(D3): 5449-5464. DOI:10.1029/93JD03517
[7] Schneider N, Lezeaux O, de La Noё J, et al. Validation of ground-based observations of stratomesospheric ozone[J]. Journal of Geophysical Research:Atmospheres, 2003, 108(D17): 4540.
[8] 张倩倩, 张兴赢. 基于卫星和地面观测的2013年以来我国臭氧时空分布及变化特征[J]. 环境科学, 2019, 40(3): 1132-1142.
Zhang Q Q, Zhang X Y. Ozone spatial-temporal distribution and trend over China since 2013:insight from satellite and surface observation[J]. Environmental Science, 2019, 40(3): 1132-1142. DOI:10.3969/j.issn.1000-6923.2019.03.029
[9] Mazzeo N A, Venegas L E, Choren H. Analysis of NO, NO2, O3 and NOx concentrations measured at a green area of Buenos Aires city during wintertime[J]. Atmospheric Environment, 2005, 39(17): 3055-3068. DOI:10.1016/j.atmosenv.2005.01.029
[10] Clapp L J, Jenkin M E. Analysis of the relationship between ambient levels of O3, NO2 and NO as a function of NOx in the UK[J]. Atmospheric Environment, 2001, 35(36): 6391-6405. DOI:10.1016/S1352-2310(01)00378-8
[11] 单文坡, 殷永泉, 杜世勇, 等. 夏季城市大气O3浓度影响因素及其相关关系[J]. 环境科学, 2006, 27(7): 1276-1281.
Shan W P, Yin Y Q, Du S Y, et al. Ozone pollution, influence factor and their correlation at urban area in summer[J]. Environmental Science, 2006, 27(7): 1276-1281. DOI:10.3321/j.issn:0250-3301.2006.07.004
[12] 陆克定, 张远航, 苏杭, 等. 珠江三角洲夏季臭氧区域污染及其控制因素分析[J]. 中国科学:化学, 2010, 40(4): 407-420.
Lu K D, Zhang Y H, Su H, et al. Regional ozone pollution and key controlling factors of photochemical ozone production in Pearl River Delta during summer time[J]. SCIENTIA SINICA Chimica, 2010, 40(4): 407-420.
[13] 沈劲, 黄晓波, 汪宇, 等. 广东省臭氧污染特征及其来源解析研究[J]. 环境科学学报, 2017, 37(12): 4449-4457.
Shen J, Huang X B, Wang Y, et al. Study on ozone pollution characteristics and source apportionment in Guangdong Province[J]. Acta Scientiae Circumstantiae, 2017, 37(12): 4449-4457.
[14] Inness A, Baier F, Benedetti A, et al. The MACC reanalysis:an 8 yr data set of atmospheric composition[J]. Atmospheric Chemistry and Physics, 2013, 13(8): 4073-4109. DOI:10.5194/acp-13-4073-2013
[15] Katragkou E, Zanis P, Tsikerdekis A, et al. Evaluation of near-surface ozone over Europe from the MACC reanalysis[J]. Geoscientific Model Development, 2015, 8(7): 2299-2314. DOI:10.5194/gmd-8-2299-2015
[16] Gaudel A, Clark H, Thouret V, et al. On the use of MOZAIC-IAGOS data to assess the ability of the MACC reanalysis to reproduce the distribution of ozone and CO in the UTLS over Europe[J]. Tellus B:Chemical and Physical Meteorology, 2015, 67(1): 27955. DOI:10.3402/tellusb.v67.27955
[17] Georgoulias A K, Tsikerdekis A, Amiridis V, et al. A 3-D evaluation of the MACC reanalysis dust product over Europe, northern Africa and Middle East using CALIOP/CALIPSO dust satellite observations[J]. Atmospheric Chemistry and Physics, 2018, 18(12): 8601-8620. DOI:10.5194/acp-18-8601-2018
[18] Flemming J, Inness A, Flentje H, et al. Coupling global chemistry transport models to ECMWF's integrated forecast system[J]. Geoscientific Model Development, 2009, 2(2): 253-265. DOI:10.5194/gmd-2-253-2009
[19] Inness A, Flemming J, Suttie M, et al. GEMS data assimilation system for chemically reactive gases[R]. Shinfield Park, Reading: ECMWF, 2009.
[20] Benedetti A, Morcrette J J, Boucher O, et al. Aerosol analysis and forecast in the European centre for medium-range weather forecasts integrated forecast system:2. Data assimilation[J]. Journal of Geophysical Research:Atmospheres, 2009, 114(D13): D13205. DOI:10.1029/2008JD011115
[21] Engelen R J, Serrar S, Chevallier F. Four-dimensional data assimilation of atmospheric CO2 using AIRS observations[J]. Journal of Geophysical Research:Atmospheres, 2009, 114: D03303.
[22] Dethof A. Assimilation of ozone retrievals from the MIPAS instrument on board ENVISAT[R]. Shinfield Park, Reading: ECMWF, 2003.
[23] Flemming J, Inness A, Jones L, et al. Forecasts and assimilation experiments of the Antarctic ozone hole 2008[J]. Atmospheric Chemistry and Physics, 2011, 11(5): 1961-1977. DOI:10.5194/acp-11-1961-2011
[24] Aumann H H, Chahine M T, Gautier C, et al. AIRS/AMSU/HSB on the aqua mission:design, science objectives, data products, and processing systems[J]. IEEE Transactions on Geoscience and Remote Sensing, 2003, 41(2): 253-264. DOI:10.1109/TGRS.2002.808356
[25] Bian J C, Gettelman A, Chen H B, et al. Validation of satellite ozone profile retrievals using Beijing Ozonesonde data[J]. Journal of Geophysical Research:Atmospheres, 2007, 112: D06305.
[26] Rajab J M, Matjafri M Z, Tan F, et al. Analysis of ozone column burden in Peninsular Malaysia retrieved from atmosphere infrared sounder (AIRS) data: 2003-2009[A]. Proceedings of 2011 IEEE International Conference on Imaging Systems and Techniques[C]. Penang, Malaysia: IEEE, 2011. 29-33.
[27] 潘亮, 牛生杰. AIRS/TOVS/TOMS反演的臭氧总量的对比分析[J]. 遥感学报, 2008, 12(1): 54-63.
Pan L, Niu S J. A Comparison of total column ozone values derived from AIRS, TOVS and TOMS[J]. Journal of Remote Sensing, 2008, 12(1): 54-63.
[28] Siddans R, Reburn W J, Kerridge B J, et al. Height-resolved ozone information in the troposphere and lower stratosphere from global ozone monitoring experiment (GOME)[R]. British: British Atmospheric Data Centre, 2007.
[29] Carli B, Alpaslan D, Carlotti M, et al. First results of MIPAS/ENVISAT with operational Level 2 code[J]. Advances in Space Research, 2004, 33(7): 1012-1019. DOI:10.1016/S0273-1177(03)00584-2
[30] Waters J W, Froidevaux L, Harwood R S, et al. The earth observing system microwave limb sounder (EOS MLS) on the aura satellite[J]. IEEE Transactions on Geoscience and Remote Sensing, 2006, 44(5): 1075-1092. DOI:10.1109/TGRS.2006.873771
[31] Bhartia P K, Wellemeyer C. TOMS-V8 total O3 algorithm, OMI algorithm theoretical basis document, Volume Ⅱ, OMI Ozone Products[R]. Greenbelt, MD: NASA Goddard Space Flight Center, 2002. 15-31.
[32] Bhartia P K, McPeters R D, Mateer C L, et al. Algorithm for the estimation of vertical ozone profiles from the backscattered ultraviolet technique[J]. Journal of Geophysical Research:Atmospheres, 1996, 101(D13): 18793-18806. DOI:10.1029/96JD01165
[33] Eskes H J, Van Der A R J, Brinksma E J, et al. Retrieval and validation of ozone columns derived from measurements of SCIAMACHY on Envisat[J]. Atmospheric Chemistry and Physics Discussions, 2005, 5(4): 4429-4475. DOI:10.5194/acpd-5-4429-2005
[34] 周秀骥, 罗超, 李维亮, 等. 中国地区臭氧总量变化与青藏高原低值中心[J]. 科学通报, 1995, 40(15): 1396-1398. DOI:10.3321/j.issn:0023-074X.1995.15.016
[35] 肖钟湧, 江洪. 利用遥感监测青藏高原上空臭氧总量30a的变化[J]. 环境科学, 2010, 31(11): 2569-2574.
Xiao Z Y, Jiang H. Variation of total ozone over Tibetain Plateau during 30 years monitored by remotely sensed data[J]. Environmental Science, 2010, 31(11): 2569-2574.
[36] 乜虹, 牛生杰, 王治邦, 等. 青藏高原清洁地区近地面层臭氧的特征分析[J]. 干旱气象, 2004, 22(1): 1-7.
Nie H, Niu S J, Wang Z B, et al. Characteristic analysis of surface ozone over clean area in Qinghai-Xizang Plateau[J]. Arid Meteorology, 2004, 22(1): 1-7. DOI:10.3969/j.issn.1006-7639.2004.01.001
[37] Coldewey-Egbers M, Weber M, Lamsal L N, et al. Total ozone retrieval from GOME UV spectral data using the weighting function DOAS approach[J]. Atmospheric Chemistry and Physics, 2005, 5(4): 1015-1025. DOI:10.5194/acp-5-1015-2005
[38] 陈爱军, 曹晓云, 韩琛惠, 等. 2000~2016年青藏高原地表反照率时空分布及动态变化[J]. 气候与环境研究, 2018, 23(3): 355-365.
Chen A J, Cao X Y, Han C H, et al. Spatial-temporal distribution and variation of land surface albedo over the Tibetan Plateau during 2000-2016[J]. Climatic and Environmental Research, 2018, 23(3): 355-365.