青藏高原臭氧谷的分布及其与太阳辐射的关系
文章快速检索     高级检索
  高原气象  2017, Vol. 36 Issue (5): 1201-1208  DOI: 10.7522/j.issn.1000-0534.2016.00106
0

引用本文 [复制中英文]

焦铂洋, 苏昱丞, 郭胜利, 等. 2017. 青藏高原臭氧谷的分布及其与太阳辐射的关系[J]. 高原气象, 36(5): 1201-1208. DOI: 10.7522/j.issn.1000-0534.2016.00106
[复制中文]
Jiao Boyang, Su Yucheng, Guo Shengli, et al. 2017. Distribution of Ozone Valley and Its Relationship with Solar Radiation over the Qinghai-Tibetan Plateau[J]. Plateau Meteorology, 36(5): 1201-1208. DOI: 10.7522/j.issn.1000-0534.2016.00106.
[复制英文]

资助项目

国家自然科学基金项目(61573193, 91537213, 41305039, 41675039);江苏高校优势学科建设工程项目(PAPD)

通讯作者

苏昱丞.E-mail:15951005455@163.com
郭胜利.E-mail:shlguo@nuist.edu.cn

作者简介

焦铂洋(1990-), 女, 河南省洛阳人, 硕士研究生, 主要从事平流层-对流层相互耦合的研究, E-mail:boyangjiao@sina.com

文章历史

收稿日期: 2016-06-23
定稿日期: 2016-10-21
青藏高原臭氧谷的分布及其与太阳辐射的关系
焦铂洋1,2, 苏昱丞1,3, 郭胜利1,3, 郭栋3, 施春华3, 李婧媛1, 苍中亚1,4, 傅帅1     
1. 空间天气研究所/南京信息工程大学数学与统计学院, 南京 210044;
2. 福建省气象培训中心, 福州 350007;
3. 南京信息工程大学气象灾害教育部重点实验室/气候与环境变化国际合作联合实验室/气象灾害预报预警与评估协同创新中心, 南京 210044;
4. 中国科学院国家天文台, 北京 100012
摘要: 使用2005-2015年夏季Aura卫星微波临边探测器(MLS)逐日臭氧观测资料, 讨论了夏季青藏高原臭氧低值区的三维分布。研究发现, 青藏高原臭氧不仅在对流层顶附近存在着臭氧低值区, 而且在平流层上层(20~1 hPa)也存在显著的低值区。高原区上空50%以上的臭氧存在于21.5~1.2 hPa的范围内, 因此平流层上层高原臭氧低值对高原臭氧谷来说也很重要。使用合成分析法对MLS夏季北半球昼夜臭氧进行研究, 结果表明该低值区仅存在于白天。根据高原区平流层上层臭氧模拟数据的集合经验模态分解, 得到IMF4的平均频率为0.09, 平均周期为11.1年, 正好对应太阳活动最强的周期。说明太阳辐射是影响高原平流层上层臭氧低值中心的一个因素。
关键词: 青藏高原    平流层中上层    臭氧低值中心    太阳辐射    
1 引言

平流层臭氧能够吸收太阳的紫外辐射, 因而被称为地球的保护伞, 平流层臭氧也可以改变平流层热力和动力结构, 对地球系统及气候变化起着重要作用, 因此臭氧成为平流层研究的热点应运而生(施春华等, 2010)。由于人类活动引起的臭氧层损耗, 尤其是极地臭氧洞的发现, 使得臭氧研究更加引起学术界的重视。

中国科学家周秀骥等(1995)通过对TOMS卫星资料分析发现, 青藏高原(下称高原)上空夏半年会形成一个臭氧总量低值区, 并推测夏半年高原上空上升运动的动力作用以及上升气流携带的对流层大气成分产生的化学作用可能是其形成的原因(周秀骥等, 2004)。Zou et al(1996, 1997)也利用TOMS、SAGE Ⅰ以及SAGE Ⅱ臭氧资料证实了高原上空存在臭氧亏损, 并给出了高原热源对该臭氧低值作用的证据。Liu et al(2003)Tian et al(2008)卞建春等(2011)以及Guo et al(2012)均强调了与大气环流有关的动力输运过程的重要作用。Liu et al(2003)Tian et al(2008)指出动力作用要比化学作用强得多。还有的讨论了该地区对流层和平流层输送对臭氧的作用(Zhang et al, 2014; 陈闯等, 2012; 陈斌等, 2011; 王卫国等, 2010; Zou et al, 2000), 郭栋等(2012)利用三维化学气候模式模拟了该现象, 并诊断其化学与动力机制, 都认为动力作用较强, 化学作用较弱。Bian et al(2011)还认为对流层顶附近臭氧低值的形成不但与动力输运有关, 还和大地形导致的整层大气减少有关, Ye and Xu(2003)指出夏季臭氧低谷的形成还与高原夏季地表的强感热加热有关, 周顺武等(2012)则认为夏季臭氧低谷的形成与对流层高度有关, Zhang et al(2017)讨论了臭氧最大加热率所在高度和臭氧最大浓度所在高度的不一致原因。这些研究从不同侧面证实了周秀骥等(2004, 1995)关于高原夏半年臭氧总量低值区动力机制的猜测。

尽管上述研究表明高原臭氧低值中心的形成因素中动力作用较强, 化学作用较弱, 然而仍有一些研究结果与此不同。付超和李维亮(1997)利用二维光化学模式模拟了高原臭氧低值区, 发现平流层中上层化学作用较强。丛春华等(2001)认为夏季高原上空存在物质输送通道, 可能会产生化学作用。刘煜等(2010)认为动力因素可能对平流层上层的高原臭氧低值作用不大。

郭栋等(2012)在用全大气气候通用模式(WACCM3, The Whole Atmosphere Community Climate Model Ⅲ)模拟研究高原臭氧谷和南亚高压的关系时意外发现, 模拟资料在高原上空的15 hPa处存在臭氧纬向偏差的低值, 该高度上的臭氧低值用上升气流、南亚高压辐散等理论无法解释。Zhang et al(2014)用WACCM3模式模拟发现地表N2O增加会使高原平流层上层臭氧减小。李书博等(2015)张芳等(2016)使用BCC-AGCM-Chem0模式对20世纪全球臭氧进行模拟, 发现对流层中高层与以下的对流层中低层的臭氧季节变化峰值存在明显的不一致性。万凌峰等(2017)使用WACCM3模式对夏季高原臭氧谷双心结构的模拟性能进行了评估。仕仁睿等(2017)研究了高原臭氧亏损变化及其对太阳活动的响应。因此, 本文利用2005-2012年夏季(6-8月)Aura卫星微波临边探测器逐日的臭氧观测资料, 分析夏季高原臭氧低值区的三维图像, 重点讨论高原臭氧谷在平流层上层的空间分布, 并讨论其形成的可能机制。

2 资料选取和方法介绍

选用的资料包括微波临边探测器(MLS, Microwave Limb Sounder)的2级轨道数据和全大气气候通用模式(WACCM 3, The Whole Atmosphere Community Climate Model Ⅲ)的输出资料(Marsh et al, 2007)。

MLS是搭载在地球观测系统(Earth Observing System, EOS)系列的第三颗卫星微风(Aura)上的一台仪器。该卫星于2004年7月15日发射成功, 飞行高度705 km, 每日绕地13~14圈, 计划运行寿命为6年。卫星探测数据可以用来研究中层大气的大气成分, 跟踪污染物的移动以及对平流层臭氧的恢复情况进行评估, 进而研究局地大气化学成分对气候变化的影响。

MLS的臭氧数据集时间范围是2004年8月8日至2016年3月22日, 并随时间自动追加, 空间范围是180°W-180°E, 82°S-82°N; 垂直分辨率215~0.2 hPa处是3 km, 0.1~0.02 hPa处是4~6 km; 臭氧分层资料的垂直有效范围是215~0.02 hPa。MLS的数据格式为.he5, 每天产生一个数据, 数据包括臭氧在各层的浓度、质量情况以及该次探测的经纬度信息。

所用的MLS臭氧数据是版本3.3(V3.3) 的2级产品, 再依据Livesey et al(2011)给出的质量控制法则对臭氧逐日轨道数据进行预处理, 剔除异常值。MLS数据水平方向沿着运行轨道每1.5°分布一条廓线, 每个轨道有240条廓线, 但各个纬度间并不是等间距的。图 1为一天内MLS对地球的探测廓线, 每一个点表示一次探测, 为便于分析, 将2005-2015年的MLS卫星的臭氧轨道廊线数据插值到2°×2°的规则格点上。数据的时间范围是2005-2015年每年夏季(6-8月), 垂直高度使用261~1 hPa的30层臭氧数据, 平均每个格点有约220个数据, 可以较准确的反应臭氧在该点的浓度。研究使用浓度的纬向偏差百分率作为基础数据进行分析。其中, 纬向偏差百分率=(格点值-纬圈平均值)/纬圈平均值。

图 1 MLS探测器2010年1月1日的对地探测分布 Figure 1 Distribution of MLS detector on 1 January 2010

在不同高度, 大气密度不同, 因而使用ppm并不能直接在垂直方向积分, 所以在需要研究臭氧垂直积分的时候, 将MLS的原始混合比数据使用公式(1)、(2)(汪自军等, 2010)转换为多布森单位进行研究。臭氧的多布森单位为臭氧在压力为标准大气压, 温度为273 K的条件下, 千分之一厘米(10 μm)的厚度。

$ \mathit{\Omega} = \sum {{n_i} \times \delta h \times {{3.7197}^{ - 12}}} \;\;, $ (1)
$ {n_i} = 7.244 \times {10^{10}}{r_i}P/T\;\;, $ (2)

式中: ri为体积混合比(单位: ppm); P为气压(单位: Pa); T为温度(单位: K); ni为数密度(单位: molecule·cm-3); δh为高度间隔(单位: km); Ω为柱量(单位: DU)。为了保持数据的统一性, 式子中计算ni所需的温度数据也使用MLS探测器的温度数据。

WACCM是基于CAM(Community Atmosphere Model)的大气化学气候模式(Collins et al, 2004)。WACCM3垂直分层66层, 范围由地面至145 km。其在100 hPa高度下使用地形追随坐标系, 100 hPa之上使用等压坐标系来简化运算。在不同的高度有不同的垂直分辨率, 65 km以上分辨率为3.5 km, 平流层顶附近为1.75 km, 低平流层区域为1.1~1.4 km; 对流层分辨率为1.1 km。为了统一坐标系便于分析数据, 使用NCAR Command Language将原来的66个垂直分层重新插值到所需要的水平分层。

本文使用的WACCM3臭氧数据为苏昱丞等(2016)对IPCC给定的三种未来情景下的全球臭氧模拟数据。WACCM3模式输出资料中A1B、A2和B1三种排放情景下全球2001-2099年来的全球臭氧总量均呈增长趋势, 即臭氧恢复, 且A2情景下柱总量增长最快, B1情景下增长最慢。这与IPCC4中给出的预测结果一致。高原模拟结果也与全球臭氧变化趋势类似, 所以模拟结果是可靠的。A1B情景为各种能源平衡发展的、经济高速增长的未来社会, 因此, 本文使用模拟结果中的A1B情景下臭氧数据作为下一步分析的源数据。

经验模态分解方法(Empirical Mode Decomposition, EMD) (姜琼妃等, 2015), 处理思路是:通过多次移动过程, 逐个分解本征模态函数(Intrinsic Mode Function, IMF)的一种新的时间序列分析方法。它可以将一个复杂的一维时间序列信号进行多尺度的傅立叶变换处理, 将不同尺度以及不同周期的波动分离出来, 得到不同周期的本征模态函数。通过多次移动, 不仅消除了信号的骑行波, 还对序列进行了平滑处理。最后变换本征模态函数, 得到其随时间变化的瞬时频率和振幅, 进而得到平均周期和平均振幅。而且因为它得到的是瞬时振幅和瞬时频率, 克服了在分辨率上不清晰的缺点。由于EMD是自适应的, 同时基于信号的局部变化特性, 因而非常适合客观处理非线性、非平稳过程。

集合经验模态分解方法(Ensemble empirical mode decomposition, EEMD)是在EMD基础上的一种白噪声辅助方法(王文等, 2014), 并对EMD的混频现象进行了改进。EEMD分解原理为在EMD分解的初试时间序列上添加随机的白噪声, 并重复进行多次该步骤。附加的白噪声, 本身均匀的分布在一维时间序列上, 在每一次EMD分解时, 白噪声就会被滤波器滤波到不同尺度的信号上, 当然需要分解的信号也就更细化平均的被滤波出来。这时, 每个样本中都含有相当的白噪声, 但是当进行分解的样本足够多时, 白噪声就会因相互抵消而消除。这样, 由于每个频率都有白噪声存在, 所以每个频率都会进行分解, 就可以解决混频问题。最后, 对全体样本取均值作为最终真正的结果。

3 夏季高原臭氧低值区的三维空间分布特征

图 2为利用MLS的臭氧浓度纬向偏差百分率数据做了高原区21.5~1 hPa和216~17 hPa的北半球臭氧的三维空间分布, 分别以-1%和-10%作为三维结构的边界, 图中地面的填色部分为地形高度, 黄色表示高出海平面的地形, 黑色的等值线为纬向偏差百分率在垂直方向的积分, 空间中灰色部分分别为纬向偏差百分率-1%(-10%)的范围。从图 2中可以看出, 高原夏季在70 hPa附近和15 hPa附近存在着不同程度的低值。臭氧在高原区150~70 hPa处有强烈的低值中心, 在北美落基山脉有较强的低值区, 这个低值区域的形成与地形高度及其产生的热力动力因素有关, 在太平洋地区也有强烈的低值区, 可能与大洋强烈的热力抬升作用相关, 低空臭氧浓度低的气团被热力抬升到高空并稀释了上空的空气, 形成了低值区。这与之前学者(周秀骥等, 1995; Liu et al, 2003; 卞建春等, 2011; 郭栋等, 2017; Li et al, 2017)对高原臭氧的研究相符合。但是在高原区平流层上层也出现了低值区域, 并且该臭氧低值区从20 hPa延伸到2 hPa, 其水平范围与高原海拔3 km以上的地形吻合很好。这个平流层上层的臭氧低值区在以上的研究中很少被提及。同时, 平流层上层的高原东北方向和大西洋上空也存在着相对较明显的低值中心, 太平洋上空的低值中心可能与夏季海面的热力抬升有关; 而高原东北部的低值中心则与氯氧化物和氮氧化物的分布有关(Zhang et al, 2014)。

图 2 2005-2015年夏季(6-8月)北半球臭氧在21.5~1.2 hPa (a)和216~17 hPa (b)上的三维廓线 Figure 2 3D profile of summer (from June to August) ozone under 21.5~1.2 hPa (a) and 216~17 hPa (b) during 2005-2015

为进一步研究高原臭氧谷在不同高度上的分布, 对20°E-180°E, 0°N-80°N范围内, 高度分别在216~17.8 hPa和21.5~1.2 hPa上的臭氧做垂直积分, 以DU为单位, 研究其不同高度内的臭氧柱总量纬向偏差在高原的分布。从图 3中可以看出, 高原在216~17.8 hPa的臭氧比纬向平均偏低约12 DU, 在21.5~1.2 hPa比纬向平均偏低约1.2 DU。不论在高原区对流层顶附近, 还是平流层上层, 臭氧的垂直积分都存在明显的低值区。但是, 平流层上层的臭氧低值中心是大气动力输送无法解释的。

图 3 2005-2015年夏季北半球臭氧纬向偏差的垂直积分(单位: DU) (a) 21.5~1.2 hPa, (b) 216~17.8 hPa.矩形代表高原的位置 Figure 3 The vertical integration of the zonal deviation of ozone in the northern hemisphere during the summer of 2005 2015. Unit: DU. The rectangle denote the location of Qinghai-Tibetan Plateau
4 太阳辐射对夏季高原平流层上层臭氧低值中心的作用

高原平流层上层的动力作用较弱, 垂直输送无法达到21.5~1.2 hPa的高度。但在这个高度上, 太阳辐射比较强, 因此, 高原平流层的臭氧低值中心很可能与较强的太阳辐射有关。为了验证高原平流层上层臭氧低值中心是否与太阳辐射有关, 使用合成分析的方法对MLS白天和夜间分布在73°E-105°E的夏季北半球臭氧数据进行比较。白天的数据使用2005-2015年6-8月每日06:00-18:00的臭氧数据, 夜间的数据使用2005-2015年夏季18:00到次日06:00的臭氧数据, 然后对白天夜间的臭氧差异进行T检验。图 4为2005-2015年北半球夏季73°E-105°E高原臭氧昼夜分布, 等值线为高原区白天臭氧的纬向偏差百分率, 填色部分为白天夜间臭氧差异的显著性检验, 通过99%显著性检验的区域使用灰色填色, 其余为白色。从图 4可以看出, 白天5.6~1 hPa的高度上有显著的臭氧低值中心, 但是夜间同一高度上该低值中心并不存在。并且高原区臭氧在10~1 hPa之间的昼夜差异极为明显。臭氧的这种昼夜不对称性, 表明太阳辐射确实对高原区平流层上层的臭氧低值有显著影响。

图 4 2005-2015年北半球夏季73°E-105°E臭氧昼(a)夜(b)分布 阴影区为通过了99%显著性水平检验的部分 Figure 4 The distribution Northern Hemisphere summer 73°E-105°E ozone day (a) night (b) from 2005 to 2015. The shaded area have passed the significant level at 99%

太阳活动最明显的年际变化为准11年的周期变化。为了研究太阳年际变化对高原区平流层上层臭氧浓度的影响, 本文使用集合经验模态分解法(EEMD)对高原平流层上层的臭氧浓度进行分析, 进一步探讨太阳活动对高原上层臭氧低值中心的作用。需要说明的是, 由于卫星观测数据的时间覆盖有限, 所以使用WACCM3模式模拟的A1B情景下臭氧作为分析的源数据, 其时间长度为2001-2099年。

通过对高原区模拟数据21.5~1.2 hPa的臭氧总量的集合经验模态分解, 观察其本征模函数, 可以检验出高原平流层上层的臭氧是否可以对应太阳活动的周期。在实际操作中, 为了避免白噪声的影响, 将集合分解重复进行了50次, EEMD结果见图 5。第一模态为平流层上层臭氧的标准化变化曲线, 模态二到模态五为分解出的本征函数, 最后一列代表的是残差。从图 5可以看出, 第二模态本征函数有准两年的年际振荡, 可能对应平流层准两年振荡(QBO); 第三模态本征函数有3~8年的年际振荡, 可能对应ENSO年际变化信号; 而第四模态本征函数则有准11年的变化周期, 对应太阳活动的11年周期。接着计算了第四模态本征函数的瞬时以及平均频率(图 6), 得到第四模态本征函数的平均频率为0.09, 平均周期为11.1年, 正好对应太阳活动最强的周期。

图 5 高原区21.5~1.2 hPa臭氧总量的EEMD分解 Figure 5 Total ozone of 21.5~1.2 hPa in plateau area based on EEMD decomposition
图 6 IMF4的年际变化曲线(a)以及IMF4的瞬时频率(b) Figure 6 The interannual variation curve (a) and the instantaneous frequency (b) of IMF4

综上所述, 太阳辐射对高原区平流层臭氧谷的形成起了一定的作用:白天高原区的平流层上层有显著的臭氧低值中心, 而夜间同一高度上并不存在低值中心; 此外, 高原区平流层上层的臭氧浓度存在准11年的周期变化。因此, 太阳辐射对高原区平流层上层臭氧低值中心的形成起了一定的作用。

5 结论

利用MLS卫星数据和WACCM3模式数据, 根据合成分析及集合经验模态分解方法, 对夏季青藏高原臭氧低值区的三维图像进行了分析, 重点讨论高原臭氧谷在平流层上层的空间分布形态, 并对其可能成因进行了研究。得出以下主要结论:

(1) 对夏季北半球臭氧的三维图像分析表明, 平流层上层20~1 hPa存在一个明显的臭氧低值中心。分析各层臭氧柱量对总柱量的相对贡献, 发现高原区上空50%以上的臭氧存在于21.5~1.2 hPa的范围内。因此, 平流层上层的高原臭氧低值对高原臭氧谷来说也很重要。

(2) 使用合成分析的方法对MLS夏季北半球昼夜臭氧进行比较可以发现, 在白天高原平流层上层有明显的臭氧低值中心, 而夜间在同一高度并不存在低值中心。这种昼夜不对称性表明太阳辐射确实对高原区平流层上层的臭氧低值有显著影响。

(3) 对2001-2099年高原区臭氧总量的模拟数据进行集合经验模态分解, 得到第四模态本征函数的平均频率为0.09, 即平均周期为11.1年, 正好对应太阳活动最强的周期。这也说明太阳活动是导致高原平流层上层臭氧低值中心形成变化的一个因素。

参考文献
Bian J C, Yan R C, Chen H B, et al. 2011. Formation of the summertime ozone valley over the Tibetan Plateau:The Asian summer monsoon and air column variations[J]. Adv Atmos Sci, 28(6): 1318–1325. DOI:10.1007/s00376-011-0174-9
Collins W D, Rasch P J, Boville B A, et al. 2004. Desciption of the NCAR community atmosphere model (CAM3. 0)[J]. National Center for Atmospheric Research Ncar Koha Opencat, tn-464+str(TN-485+STR). https://www.researchgate.net/publication/224017933_Description_of_the_NCAR_community_atmosphere_model_CAM_30
Guo D, Wang P X, Zhou X J, et al. 2012. Dynamic effects of the south Asian high on the ozone valley over the Tibetan Plateau[J]. Acta Meteor Sinica, 26(2): 216–228. DOI:10.1007/s13351-012-0207-2
Li Z K, Qin H, Guo D, et al. 2017. Impact of Ozone Valley over the Tibetan on South Asian High in CAM5[J]. Adv Meteorol, in press. https://www.hindawi.com/journals/amete/2017/9383495/ref/
Liu Y, Li W L, Zhou X J, et al. 2003. Mechanism of formation of the ozone valley over the Tibetan Plateau in summer-transport and chemical process of ozone[J]. Adv Atmos Sci, 20(1): 103–109. DOI:10.1007/BF03342054
Livesey N J, Read W G, Froidevaux L, et al. 2011. EOS MLS version 3. 3 level 2 data quality and description document report[R]. Jet Propulsion Laboratory.
Marsh D R, Garcia R R, Kinnison D E, et al. 2007. Modeling the whole atmosphere response to solar cycle changes in radiative and geomagnetic forcing[J]. J Geophys Res Atmos, 112(D23): 6033–6044. DOI:10.1029/2006JD008306
Tian W S, Chipperfield M, Huang Q. 2008. Effects of the Tibetan Plateau on total column ozone distribution[J]. Tellus B, 60(4): 622–635. DOI:10.1111/j.1600-0889.2008.00338.x
Ye Z J, Xu Y F. 2003. Climate characteristics of ozone over Tibetan Plateau[J]. J Geophys Res, 108(D20): 1–15. DOI:10.1029/2002JD003139
Zhang J K, Tian W S, Xie F, et al. 2014. Climate warming and decreasing total column ozone over the Tibetan Plateau during winter and spring[J]. Tellus Series B-Chemical and Physical Meteorology, 66(8): 136–140. DOI:10.3402/tellusb.v66.23415
Zhang F, Hou C, Li J N, et al. 2017. A simple parameterization for the height of maximum ozone heating rate[J]. Infrared Phys Technol, accept. https://www.researchgate.net/publication/319936007_A_simple_parameterization_for_the_height_of_maximum_ozone_heating_rate
Zou H, Gao Y Q. 1997. Vertical ozone profile over Tibet using sage Ⅰ and Ⅱ data[J]. Adv Atmos Sci, 14(4): 505–512. DOI:10.1007/s00376-997-0068-z
Zou H, Ji C P, Zhou L B. 2000. QBO signal in total ozone over Tibet[J]. Adv Atmos Sci, 17(4): 562–568. DOI:10.1007/s00376-000-0019-4
Zou H. 1996. Seasonal variation and trends of TOMS ozone over tibet[J]. Geophys Res Lett, 23(9): 1029–1032. DOI:10.1029/96GL00767
卞建春, 严仁嫦, 陈洪滨. 2011. 亚洲夏季风是低层污染物进入平流层的重要途径[J]. 大气科学, 35(5): 897–902. DOI:10.3878/j.issn.1006-9895.2011.05.09 Bian Jianchun, Yan Renchang, Chen Hongbin. 2011. Tropospheric pollutant transport to the stratosphere by Asian summer monsoon[J]. Chinese J Atmos Sci, 35(5): 897–902. DOI:10.3878/j.issn.1006-9895.2011.05.09
陈斌, 施晓晖, 徐祥德, 等. 2011. 利用卫星大气成分资料分析夏季亚洲季风区平流层-对流层输送特征[J]. 高原气象, 30(1): 65–73. Chen Bin, Shi Xiaohui, Xu Xiangde, et al. 2011. Characteristics of troposphere stratosphere transports over Asian monsoon region in summer based on Aura-MLS data[J]. Plateau Meteor, 30(1): 65–73.
陈闯, 田文寿, 田红瑛, 等. 2012. 青藏高原东北侧臭氧垂直分布与平流层-对流层物质交换[J]. 高原气象, 31(2): 295–303. Chen Chuang, Tian Wenshou, Tian Hongying, et al. 2012. Vertical distribution of ozone and stratosphere-troposphere exchanges on the northeastern side of Tibetan Plateau[J]. Plateau Meteor, 31(2): 295–303.
丛春华, 李维亮, 周秀骥. 2001. 青藏高原及其邻近地区上空平流层-对流层之间大气的质量交换[J]. 科学通报, 46(22): 1914–1918. DOI:10.3321/j.issn:0023-074X.2001.22.017 Cong Chunhua, Li Weiliang, Zhou Xiuji. 2001. Mass exchange between the stratosphere and troposphere over the Tibetan Plateau and its adjacent regions[J]. Chinese Sci Bull, 46(22): 1914–1918. DOI:10.3321/j.issn:0023-074X.2001.22.017
付超, 李维亮. 1997. 夏季青藏高原上空臭氧总量低值区形成的模拟试验[C]周秀骥编著. 中国地区大气臭氧变化及其对气候环境的影响(二). 北京: 气象出版社, 274-285. Fu Chao, Li Weiliang. 1997. Simulation experiment on the formation of ozone low value area over the Tibetan Plateau in summer[C] Zhou Xiuji. eds. Effects of atmospheric ozone variation on Chinese climate and environment (two). Beijing:China Meteorological Press, 274-285. http: //d. wanfangdata. com. cn/Conference/39037
郭栋, 周秀骥, 刘煜, 等. 2012. 南亚高压对青藏高原臭氧谷的动力作用[J]. 气象学报, 70(6): 1302–1311. DOI:10.11676/qxxb2012.109 Guo Dong, Zhou Xiuji, Liu Yu, et al. 2012. The dynamic effects of the south Asian high on the ozone valley over the Tibetan Plateau[J]. Acta Meteor Sinica, 70(6): 1302–1311. DOI:10.11676/qxxb2012.109
郭栋, 徐建军, 苏昱丞, 等. 2017. 青藏高原和北美夏季臭氧谷垂直结构和形成机制的比较[J]. 大气科学学报, 40(3): 412–417. Guo Dong, Xu Jianjun, Su Yuchen, et al. 2017. Comparison of vertical structure and formation mechanism of summer ozone valley over the Tibetan Plateau and North America[J]. Trans Atmos Sci, 40(3): 412–417.
姜琼妃, 郭胜利, 梁浩文, 等. 2015. 基于EMD的东亚夏季风年代际变化特征及与太阳活动的关系[J]. 空间科学学报, 35(1): 50–55. DOI:10.11728/cjss2015.01.050 Jiang Qiongfei, Guo Shengli, Liang Haowen, et al. 2015. Research of EASM decadal variation characteristics based on the EMD method and the relationship with solar activity[J]. Chinese J Space Sci, 35(1): 50–55. DOI:10.11728/cjss2015.01.050
李书博, 吴统文, 张洁, 等. 2015. BCC-AGCM-Chem0模式对20世纪全球O3气候平均态及季节变化特征的模拟研究[J]. 高原气象, 34(6): 1601–1615. DOI:10.7522/j.issn.1000-0534.2014.00119 Li Shubo, Wu Tongwen, Zhang Jie, et al. 2015. Simulation study about climatological basic state and seasonal variations of global O3 in the 20th century[J]. Plateau Meteor, 34(6): 1601–1615. DOI:10.7522/j.issn.1000-0534.2014.00119
刘煜, 李维亮, 周秀骥. 2010. 非均相化学过程在青藏高原臭氧低谷形成中的作用[J]. 气象学报, 68(6): 836–846. DOI:10.11676/qxxb2010.079 Liu Yu, Li Weiliang, Zhou Xiuji. 2010. A possible effect of heterogeneous reactions on the formation of the ozone valley over the Tibetan Plateau[J]. Acta Meteor Sinica, 68(6): 836–846. DOI:10.11676/qxxb2010.079
施春华, 陈月娟, 郑彬, 等. 2010. 平流层臭氧季节变化的动力和光化学作用之比较[J]. 大气科学, 34(2): 399–406. DOI:10.3878/j.issn.1006-9895.2010.02.13 Shi Chunhua, Chen Yuejuan, Zheng Bin, et al. 2010. A comparison with the contribution of dynamics and chemistry in ozone'seasonal variation in the stratosphere[J]. Chinese J Atmos Sci, 34(2): 399–406. DOI:10.3878/j.issn.1006-9895.2010.02.13
苏昱丞, 郭栋, 郭胜利, 等. 2016. 未来百年夏季青藏高原臭氧变化趋势及可能机制[J]. 大气科学学报, 39(3): 309–317. DOI:10.13878/j.cnki.dqkxxb.20140925002 Su Yucheng, Guo Dong, Guo Shengli, et al. 2016. Ozone trends over the Tibetan Plateau in the next 100 Years and their possible mechanism[J]. Trans Atmos Sci, 39(3): 309–317. DOI:10.13878/j.cnki.dqkxxb.20140925002
仕仁睿, 周顺武, 孙绩华, 等. 2017. 青藏高原臭氧亏损变化及其对太阳活动的响应[J]. 云南大学学报(自然科学版), 39(1): 78–87. Shi Renrui, Zhou Shunwu, Sun Jihua, et al. 2017. On variation features of total ozone depletion over the Tibetan Plateau and its responds to solar activity[J]. Journal of Yunnan University (Natural Science), 39(1): 78–87.
汪自军, 陈圣波, 吕航, 等. 2010. 大气臭氧不同量纲之间的转换[J]. 气象与环境学报, 26(2): 63–67. DOI:10.3969/j.issn.1673-503X.2010.02.012 Wang Zijun, Chen Shengbo, Lü Hang, et al. 2010. Conversion on different dimensions of atmospheric ozone[J]. J Meteor Environ, 26(2): 63–67. DOI:10.3969/j.issn.1673-503X.2010.02.012
万凌峰, 郭栋, 刘仁强, 等. 2017. WACCM3对夏季青藏高原臭氧谷的双心结构的模拟性能评估[J]. 高原气象, 36(1): 57–66. DOI:10.7522/j.issn.1000-0534.2016.00004 Wan Lingfeng, Guo Dong, Liu Renqiang, et al. 2017. Valuation of the WACCM3 performance on simulation of the double core of ozone valley over the Qinghai-Xizang Plateau in summer[J]. Plateau Meteor, 36(1): 57–66. DOI:10.7522/j.issn.1000-0534.2016.00004
王卫国, 梁俊平, 王颢樾, 等. 2010. 青藏高原及附近区域穿越对流层顶的质量和臭氧通量研究[J]. 高原气象, 29(3): 554–562. Wang Weiguo, Liang Junping, Wang Yingyue, et al. 2010. Comparative study on mass and ozone fluxes cross tropopause over Qinghai-Xizang Plateau and its surrounding areas[J]. Plateau Meteor, 29(3): 554–562.
王文, 任冉, 李耀辉. 2014. 基于EEMD的黄河中上游夏季降水预报方法的研究[J]. 气象科学, 34(3): 261–266. DOI:10.3969/2014jms.0004 Wang Wen, Ren Ran, Li Yaohui. 2014. Summer rainfall prediction in the upper and middle reaches of the Yellow River with EEMD method[J]. J Meteor Sci, 34(3): 261–266. DOI:10.3969/2014jms.0004
张芳, 吴统文, 张洁, 等. 2016. BCC-AGCM-Chem0模式对20世纪对流层臭氧变化趋势的模拟研究[J]. 高原气象, 35(1): 158–171. DOI:10.7552/j.issn.1000-0534,2014,00118 Zhang Fang, Wu Tongwen, Zhang Jie, et al. 2016. Variations of tropospheric ozone in the 20th century sinmulated by BCC-ACGM-Chem0 model[J]. Plateau Meteor, 35(1): 158–171. DOI:10.7552/j.issn.1000-0534,2014,00118
周顺武, 杨双艳, 张人禾, 等. 2012. 近30年青藏高原臭氧总量亏损的可能原因及其与对流层顶高度的联系[J]. 高原气象, 31(6): 1471–1478. Zhou Shunwu, Yang Shuangyan, Zhang Renhe, et al. 2012. Possible causes of total ozone depletion over the Qinghai-Xizang Plateau and its relation to troposphere height in recent 30 years[J]. Plateau Meteor, 31(6): 1471–1478.
周秀骥, 李维亮, 陈隆勋, 等. 2004. 青藏高原地区大气臭氧变化的研究[J]. 气象学报, 62(5): 513–527. DOI:10.11676/qxxb2004.053 Zhou Xiuji, Li Weiliang, Chen Longxun, et al. 2004. Study of ozone change over Tibetan Plateau[J]. Acta Meteor Sinica, 62(5): 513–527. DOI:10.11676/qxxb2004.053
周秀骥, 史久恩. 1995. 中国地区臭氧总量变化与青藏高原低值中心[J]. 科学通报, 40(15): 1396–1398. Zhou Xiuji, Shi Jiuen. 1995. Variation of total ozone over China and ozone low center over the Tibetan Platea[J]. Chinese Sci Bull, 40(15): 1396–1398.
Distribution of Ozone Valley and Its Relationship with Solar Radiation over the Qinghai-Tibetan Plateau
JIAO Boyang1,2 , SU Yucheng1,3 , GUO Shengli1,3 , GUO Dong3 , SHI Chunhua3 , LI Jingyuan1 , CANG Zhongya1,4 , FU Shuai1     
1. Institute of Space Weather/School of Mathematics and Statistics, Nanjing University of Information Science and Technology, Nanjing 210044, China;
2. Training Center of Fujian Meteorology, Fuzhou 350007, China;
3. Nanjing University of Information Science and Technology Key Laboratory of Meteorological Disaster, Ministry of Education(KLME)/Joint International Research Laboratory of Climate and Environment Change(ILCEC)/Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters(CIC-FEMD), Nanjing 210044, China;
4. National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, China
Abstract: Using Aura satellite 10 years (2005-2015) summer (from June to August) microwave limb sounder (MLS) daily ozone observation data, the three-dimensional distribution of the core of ozone valley over the Qinghai-Tibetan Plateau in summer is analyzed.The results show that the ozone valley in the Qinghai-Tibetan Plateau is not only in the tropopause, but also in the upper stratosphere (20~1 hPa).More than 50% of the ozone over the plateau is in the range of 21.5~1.2 hPa, so for ozone valley upper stratospheric ozone plateau low plateau is also very important.Then by using synthetic analysis method, MLS ozone in the Northern Hemisphere summer day and night were compared, and we can see there is a core of ozone valley on the plateau during the day.Next, by using Ensemble Empirical Mode Decomposition (EEMD) to process total ozone analog data in the Plateau, the average frequency of IMF4 is 0.09, that is to say the average period is 11.1 years, exactly correspond to the strongest solar activity cycle.So solar activity is a factor for the core of ozone valley in the plateau upper stratosphere.
Key Words: Qinghai-Tibetan plateau    The middle and upper stratosphere    The core of ozone valley    Solar radiation