2018, Vol. 42
2 上海市气候中心, 上海 200030
3 上海市公共气象服务中心, 上海 200030
, GUO Dong1, SHI Chunhua1
, LI Zhenkun2, ZHOU Shunwu1, HUANG Ying1, SU Yucheng1, WANG Linwei3 2 Shanghai Climate Center, Shanghai 200030
3 Shanghai Public Meteorological Service Center, Shanghai 200030
北半球夏半年，南亚高压因亚洲季风区的对流加热作用形成并维持在上对流层下平流层(UTLS)区域(Hoskins and Rodwell, 1995；Highwood and Hoskins, 2010；Jin et al., 2013)。同时，其发展和变化也受到青藏高原热力作用的影响(Ye, 1981；Yeh, 1982；黄荣辉，1985；吴国雄等，1997；刘屹岷等，1999；吴国雄等，2002；Duan and Wu, 2005；刘伯奇等，2009；周秀骥等，2009)。在南亚高压反气旋环流影响下，其邻近UTLS区的大气成分分布与变化具有独特的区域特征(陈洪滨等，2006)。周秀骥等(1995)利用1979~1991年平均的TOMS卫星资料发现在夏半年(4~9月)青藏高原上空的臭氧总量比同纬度的中国东部地区要低5 DU(多布森)以上，称其为青藏高原臭氧谷。Zou(1996)同样利用TOMS资料，证实了青藏高原臭氧谷的存在。随着观测资料的丰富，Guo et al., (2015)利用卫星资料和再分析资料发现了青藏高原臭氧谷的双心结构，最强中心位于UTLS区，较弱的中心位于平流层上层，并指出该弱中心主要由化学作用导致(郭栋等，2017)。关于青藏高原臭氧谷的形成与变化特征，前人已经做了许多研究(刘煜和李维亮，2001；苏昱丞等，2016；郭栋等，2017；Shi et al., 2017；仕仁睿等，2017；万凌峰等，2017)，并且大多数工作认为与南亚高压有关的动力过程对UTLS区的臭氧低值形成起到主要作用，而化学作用较弱(郭世昌和徐裕华，1986；Liu et al., 2003；苏绍基和王卫国，2004；周任君和陈月娟，2005；Tian et al., 2008；Liu et al., 2009；Bian et al., 2011；Guo et al., 2012；苏昱丞等，2016；Guo et al., 2017)。其中，对流层顶高度变化与大尺度环流变化导致的UTLS区臭氧异常在南亚高压影响其邻近区臭氧低值的变化中起着关键作用(刘煜和李维亮，2001；周任君和陈月娟，2005；Zhou et al., 2013；Zhang et al., 2014, 2015)。另外，高原排开空气导致臭氧较同纬度少也是青藏高原臭氧谷形成原因之一(Tian et al., 2008；Bian et al., 2011)。
然而，以上研究仅关注了南亚高压对其邻近UTLS区臭氧分布的影响，并没有讨论南亚高压邻近区域臭氧对南亚高压的反作用。事实上，臭氧可以吸收短波和长波辐射并释放长波辐射，是平流层最重要的热源。310~400 nm的太阳辐射称为Huggins臭氧吸收带，主要由UTLS区的臭氧吸收(Brasseur and Solomon, 2005)，而且在下平流层区，由于臭氧变化所引起的冷却效应与水汽增加所引起的长波辐射冷却效应是相当的(de F. Forster and Shine, 1999)。因此UTLS区臭氧可以通过影响辐射平衡，来影响大气温度结构和大气环流(Ramaswamy et al., 1996；Bengtsson et al., 1999；Hu et al., 2015)。而UTLS区臭氧的纬向非对称性会在南亚高压区产生纬向非对称的辐射强迫，进而可能会引起南亚高压的变化(Li et al., 2017)。因此，本文研究南亚高压与其邻近区域臭氧变化的相互作用是有意义的。
另外，皮尔逊相关系数只反映了两变量间的相关性，无法判断两者间的因果联系。在以往气候研究中往往采用时间滞后相关分析的方法判断时间序列间的因果关系，但由于相关分析不区分方向性，在应用上无法处理许多问题，例如重复性出现的过程。而Liang-Kleeman的信息流方法(Liang, 2008, 2013, 2014)，则可以用给定时间序列间单位时间内传递的信息来表征两者间的因果关系。因此，本文利用Liang-Kleeman的信息流方法研究了南亚高压与其邻近区域臭氧变化的因果联系。
2 资料与方法本文采用欧洲中期天气预报中心(ECMWF)提供的1979~2015年夏半年(5~9月)ERA-interim再分析月平均资料，包括100 hPa高度场和臭氧质量混合比(水平分辨率为0.5°×0.5°)，以及散度、垂直速度、长波加热率和短波加热率(水平分辨率为1°×1°)。
本文利用的统计分析方法有皮尔逊相关系数、t检验、Liang-Kleeman信息流方法，其中，因果关系定量地由两时间序列间单位时间内传递的信息来表征。一个事件到另一事件存在非零的信息流，或者说信息传递，逻辑上表明了前者与后者因果联系的强弱，若不存在因果联系则信息流为零。Liang-Kleeman信息流与因果关系的联系由一个重要的性质保证，该性质也是经过证明在因果分析中必须遵循的定理：一个动态事件A的发展独立于另一事件B，则从B到A的信息流是零。给出一个二维系统：
| $ \left\{ \begin{array}{l} \frac{{{\rm{d}}{x_1}}}{{{\rm{d}}t}} = {F_1}\left({{x_1}, {x_2}, t} \right) + {b_{11}}{{\dot w}_1} + {b_{12}}{{\dot w}_2}, \\ \frac{{{\rm{d}}{x_2}}}{{{\rm{d}}t}} = {F_2}\left({{x_1}, {x_2}, t} \right) + {b_{21}}{{\dot w}_1} + {b_{22}}{{\dot w}_2}, \end{array} \right. $ | (1) |
其中，
| $ {T_{2 \to 1}} = - E\left({\frac{1}{{{\rho _1}}}\frac{{\partial \left({{F_1}{\rho _1}} \right)}}{{\partial {x_1}}}} \right) + \frac{1}{2}E\left({\frac{1}{{{\rho _1}}}\frac{{{\partial ^2}\left({b_{11}^2 + b_{12}^2} \right){\rho _1}}}{{\partial x_1^2}}} \right), $ | (2) |
其中，ρ1为x1的边际密度，E为数学期望(单位：nats/单位时间，nats是自然信息量单位)。显然，方程(2)包含了上述因果联系的严格定理，即若x1的演变独立于x2，则T2→1=0。
当只给出两时间序列X1和X2，我们首先需要一个计算信息流的模型，对于线性系统，方程(2)的极大似然估计在形式上将变得非常简洁(Liang，2014)：
| $ {T_{2 \to 1}} = \frac{{{C_{11}}{C_{12}}{C_{2, d1}} - C_{12}^2{C_{1, d1}}}}{{C_{11}^2{C_{22}} - {C_{11}}C_{12}^2}}, $ | (3) |
式中，Cij为Xi与Xj的样本协方差，Ci, dj为Xi与
| $ {\dot X_{j, n}} = \frac{{{X_{j, n + k}} - {X_{j, n}}}}{{k\Delta t}}, $ | (4) |
Δt为时间步长。若T2→1，则X2不会引起X1的变化，即X2不是X1变化的原因。若T2→1＞0，则X2是X1变化的原因，且X2使X1变化的不确定性增强，可预报性降低。若T2→1＜0，则X2是X1变化的原因，且X2使X1变化趋于稳定，可预报性增强。相应地，Liang(2014)给出了检验信息流显著性的方法。另外，该方法已经在研究El Niño与印度洋偶极子(IOD)间因果关系时得到验证(Liang，2014)。
本文采用强度指数作为南亚高压的特征指数。参考张琼等(2000)对南亚高压强度指数的定义：100 hPa等压面上，(0°~40°N，40°E~120°E)范围内16600 gpm以上所有格点的位势高度减去16600 gpm差值的总和。如上定义的强度指数具有明显的年际和年代际变化，在以往的研究中被多次应用来代表南亚高压的变化特征(张琼和吴国雄，2001；徐忠峰和钱永甫，2003；陈延聪等，2009；周利敏等，2016)。同时，本文利用均方差来表征100 hPa高度上位势高度与臭氧质量混合比的纬向偏差(以下简称臭氧纬向偏差)的变化幅度。
3 南亚高压与邻近地区臭氧的时空分布特征南亚高压在夏半年一直存在，除5月偏南以外，6~9月均位于伊朗和青藏高原上空。5~9月南亚高压特征等值线分别取为16620、16720、16760、16760、16690 gpm，代表南亚高压的范围。5月份南亚高压自西向东横跨印度半岛、孟加拉湾、中南半岛直至西太平洋，中心位于(21.5°N，95.5°E)，强度为16674 gpm(图 1a)。到了6月，南亚高压西北抬，高压中心位于青藏高原上空，中心位于(29°N，88°E)，高压强度有所增强，强度增大至16776 gpm(图 1b)。7月份南亚高压达到最西北位置，中心位于伊朗高原上空(33°N，58°E)，强度达到最强，为16855 gpm(图 1c)。8月份南亚高压大体位置与强度稍变，中心位于(32°N，80.5°E)，强度为16827 gpm(图 1d)。到了9月份，南亚高压南撤，中心南撤至(29°N，88.5°E)，强度减弱为16721 gpm(图 1e)。
|
图 1 (a) 5月、(b)6月、(c)7月、(d)8月、(e)9月月平均的100 hPa位势高度(填色，单位：gpm)、位势高度的均方差(等值线，单位：gpm)，以及相应的(f–j)臭氧质量混合比纬向偏差(填色，单位：10-7 kg kg-1)、臭氧质量混合比纬向偏差均方差(等值线，单位：10-7 kg kg-1)。图a–e中粗实线代表南亚高压的范围，图f–j中粗虚线代表臭氧低值的范围，下同 Figure 1 Distributions of monthly mean 100-hPa geopotential height (shaded, units: gpm) and its standard deviation (contours, units: gpm) in (a) May, (b) June, (c) July, (d) August, (e) September, (f–j) corresponding zonal deviation of ozone mass mixing ratio (shaded, units: 10-7 kg kg-1) and its standard deviation (contours, units: 10-7 kg kg-1). In Figs. a–e, solid lines indicate the region of South Asia high, in Figs. f–j, thick dashed lines indicate the region of ozone low, the same below |
同时，南亚高压区臭氧低值(以下简称臭氧低值)也一直存在，均位于伊朗和青藏高原上空。取-1.5×10-7 kg kg-1为臭氧纬向偏差的特征等值，代表 100 hPa上臭氧低值的范围。对比图 1a–e可以看出，夏半年各月份在南亚高压及邻近区总存在臭氧纬向偏差负值区。5月份负值中心位于青藏高原北部上空(图 1f)，而到了6月份负值中心西进且有一分为二，一个中心位于青藏高原西北侧而另一个中心位于伊朗高原北部(图 1g)。7月份负值中心合二为一且强度减弱，其位置达到最西，位于里海与伊朗高原交汇处(图 1h)。8月份负值中心减弱并且东退(图 1i)。9月份负值中心东退至了青藏高原西北侧(图 1j)。
由100 hPa上位势高度的均方差分布可以看出，南亚高压的变化情况总体上是高纬地区变化大而低纬地区变化较小，并且5月份在南亚高压的西北部存在变化大值中心(图 1a)，7月份在南亚高压的北部及东北部也存在变化大值中心(图 1c)。而臭氧纬向偏差的均方差分布的总体情况与南亚高压的类似(图 1a–e)。高纬地区臭氧纬向偏差变化大而低纬地区变化较小，并且在南亚高压北部，位势高度场与臭氧纬向偏差的均方差具有相似的分布。
4 南亚高压与邻近UTLS区臭氧低值的联系为进一步研究南亚高压强度变化与同期臭氧低值变化的相关及因果联系，即探究南亚高压强度(臭氧低值)的变化是否能引起臭氧低值(南亚高压强度)的变化，或者是否两者相互影响，计算了1979~2015年夏半年(5~9月)各月份100 hPa南亚高压强度指数与其邻近地区臭氧纬向偏差变化的相关系数(图 2)与信息流(图 3、4)。
|
图 2 (a) 5月、(b)6月、(c)7月、(d)8月、(e)9月南亚高压强度与同期臭氧质量混合比纬向偏差相关系数分布。阴影由浅到深分别为通过80%、90%、95%信度水平的显著性检验区域 Figure 2 Distributions of correlation coefficient between the South Asia high intensity and the zonal deviation of ozone mass mixing ratio in (a) May, (b) June, (c) July, (d) August, (e) September. Shaded areas from light to dark indicate the values exceed the 80%, 90%, 95% confidence levels, respectively |
|
图 3 (a) 5月、(b)6月、(c)7月、(d)8月、(e)9月南亚高压强度对臭氧质量混合比纬向偏差的信息流分布。彩色阴影部分通过80%信度水平的显著性检验 Figure 3 Distributions of information flow from intensity of South Asia high to the zonal deviation of ozone mass mixing ratio in (a) May, (b) June, (c) July, (d) August, (e) September. Shaded areas indicate the values exceed the 80% confidence level |
|
图 4 (a) 5月、(b)6月、(c)7月、(d)8月、(e)9月臭氧质量混合比纬向偏差对南亚高压强度的信息流分布。彩色阴影部分通过80%信度水平的显著性检验 Figure 4 Distributions of information flow from zonal deviation of ozone mass mixing ratio to the intensity of South Asia high in (a) May, (b) June, (c) July, (d) August, (e) September. Shaded areas indicate the values exceed the 80% confidence level |
5月份南亚高压较臭氧低值区域偏南。在阿拉伯海、印度半岛以及中南半岛地区，即南亚高压西部和臭氧低值南部，南亚高压强度与臭氧纬向偏差呈显著负相关，相关系数达到-0.35，并通过95%信度水平的显著性检验(图 2a)。而从信息流分布可以看到，5月份在南亚高压西部和臭氧低值区域南部，南亚高压强度变化对该地区臭氧纬向偏差产生影响，并使其趋于不稳定(图 3a)；而该地区臭氧纬向偏差的变化反过来也影响南亚高压强度变化，使南亚高压强度趋于稳定，可预报性增强(图 4a)。即5月，在南亚高压西部和臭氧低值区域南部，臭氧低值增强(减弱)可能部分导致了南亚高压的增强(减弱)。
到了6月，由图 2b可以看到南亚高压和臭氧低值的位置范围较一致，青藏高原西部、孟加拉湾北部以及中南半岛地区为显著负相关，青藏高原西部和孟加拉湾北部地区相关系数达到-0.35，中南半岛地区相关系数达-0.4；在南亚高压西侧的西亚地区向北延伸至里海存在一个负相关中心，相关系数达到-0.45，以上区域均通过95%信度水平的显著性检验。并且可以看到信息流高绝对值中心与强相关中心位置对应，在青藏高原、孟加拉湾、中南半岛以及南亚高压西侧的西亚地区向北延伸至里海地区，南亚高压强度变化是该地区臭氧纬向偏差变化的部分原因，并使其可预报性增强(图 3b)；而只有伊朗高原与青藏高原交界地区的臭氧纬向偏差对南亚高压有影响，使其强度趋于稳定(图 4b)。即6月，在南亚高压和臭氧低值的中部和西部边缘，臭氧低值增强(减弱)可能和南亚高压增强(减弱)互为因果，相互影响，而在南亚高压和臭氧低值的东南侧，南亚高压增强(减弱)可能部分导致了臭氧低值增强(减弱)。
7月份南亚高压最强，且南亚高压的位置与臭氧低值较为一致。然而，南亚高压和臭氧低值的关系并不密切，南亚高压强度和臭氧纬向偏差相关系数和信息流通过检验的区域均位于南亚高压和臭氧低值边缘与外部(图 2c、3c、4c)。因此，认为7月南亚高压和臭氧低值年际变化的因果关系不明显。
8月份南亚高压强度与南亚高压和臭氧低值西部的伊朗高原和东部的河西走廊地区臭氧纬向偏差为负相关，相关系数分别达到-0.35和-0.25，分别通过95%和80%信度水平的显著性检验。信息流显示，在伊朗高原、河西走廊地区，南亚高压强度变化是该地区臭氧纬向偏差变化的部分原因，并使其可预报性增强(图 3d)；而在南亚高压和臭氧低值主体上，臭氧纬向偏差变化对南亚高压的影响并不显著(图 4d)。因此，8月在南亚高压和臭氧低值的西部和东部，南亚高压增强(减弱)可能部分导致了臭氧低值增强(减弱)。
9月份，臭氧低值中心位于南亚高压北部。南亚高压强度和臭氧纬向偏差在臭氧低值中心的塔里木盆地和南亚高压西部的伊朗高原地区均存在显著负相关，相关系数达到-0.4，且均通过95%信度水平的显著性检验(图 2e)。在伊朗高原地区，南亚高压强度变化可能是该地区臭氧纬向偏差变化的部分原因，并使其趋于稳定，与其为单向因果关系(图 3e)；在臭氧低值中心的塔里木盆地地区，南亚高压与臭氧纬向偏差的变化互为因果，均使对方可预报性增强(图 3e、4e)。即9月，在南亚高压北部和臭氧低值中心区，臭氧低值增强(减弱)可能和南亚高压增强(减弱)互为因果，而在南亚高压西部南亚高压增强(减弱)可能部分导致了臭氧低值增强(减弱)。
为进一步分析在垂直方向上南亚高压与臭氧低值的相互作用，分别取上述分析中相互作用较为明显的5月、6月以及9月的关键区取垂直剖面(图 5)。
|
图 5 (a) 5月、(b)6月、(c)9月关键区垂直剖面上(5月、6月和9月分别沿75°E，80°E和55°E做垂直剖面)南亚高压强度对臭氧质量混合比纬向偏差的信息流分布以及(d)5月、(e)6月、(f)9关键区垂直剖面上(5月、6月和9月分别沿70°E、72°E和80°E做垂直剖面)臭氧质量混合比纬向偏差对南亚高压强度的信息流分布。彩色阴影部分通过80%信度水平的显著性检验 Figure 5 Distributions of information flow from intensity of South Asia high to the zonal deviation of ozone mass mixing ratio on the vertical sections (along 75°E, 80°E, and 55°E in May, June, and September) in the key region in (a) May, (b) June, (c) September. Distributions of information flow from zonal deviation of ozone to the intensity of South Asia high on the vertical sections (along 70°E, 72°E, and 80°E in May, June, and September) in the key region in (d) May, (e) June, (f) September. Shaded areas indicate the values exceed the 80% confidence level |
在南亚高压强度对臭氧低值信息流方向，5月、6月和9月分别沿75°E，80°E和55°E做垂直剖面。可以看到，在5月份，南亚高压强度变化对125~100 hPa以及200 hPa以下臭氧低值具有影响，且使臭氧低值的变化趋于不稳定(图 5a)。而6月份，南亚高压强度变化主要影响了100 hPa附近的臭氧低值变化，且使其变化趋于稳定，可预报性增强(图 5b)。到了9月份，南亚高压强度变化主要影响了关键区内300~100 hPa的臭氧低值(图 5c)。
而臭氧低值对南亚高压强度的信息流，5月、6月和9月分别沿70°E、72°E和80°E做垂直剖面。可以看到在5月和6月，100 hPa以下臭氧低值的变化对南亚高压强度变化具有影响，且使其趋于稳定，可预报性增强(图 5d、5e)，而9月份，臭氧低值对南亚高压强度的影响在125 hPa以上为使其变化趋于稳定，125 hPa以下为使其变化趋于不稳定(图 5f)。
5 南亚高压与邻近UTLS区臭氧低值相互影响机制关于南亚高压对臭氧低值的影响机制，前人已经有了丰富的研究。多数研究认为，动力过程对臭氧低值的形成及变化具有重要影响，而南亚高压是其中的主要影响因素(卞建春等，1997；Liu et al., 2003；Tian et al., 2008；Guo et al., 2012)。二维全球平流层—对流层动力、辐射和光化学耦合模式(付超等，1997)以及三维化学模式OSLO CT M2(Liu et al., 2003)的模拟结果指出，5月份南亚高压从低纬度带来的低浓度臭氧是臭氧低值中心初期建立的主要原因，之后随着南亚高压移上高原和强度增强，对流层低层低浓度的臭氧被输送到UTLS区，从而使得臭氧低值得以加强和维持。同样的，周任君和陈月娟(2005)指出，南亚高压增强时，120~50 hPa的等位温面上升，对应上升气流，从而推测此上升气流带来下层低臭氧浓度的空气从而使该层臭氧浓度减少。周秀骥等(2004)也对南亚高压影响臭氧低值的可能机制进行了详细的讨论，因此本文不再详述，以下着重分析臭氧低值影响南亚高压的可能机制。
为分析臭氧低值影响南亚高压的可能机制，我们选取5月、6月以及9月臭氧低值影响南亚高压的关键区域进行分析。根据上述分析中臭氧低值影响南亚高压的区域，5月选取区域为(11°~23°N，57°~85°E)，6月为(17°~37°N，60°~80°E)，9月为(30°~42°N，70°~95°E)。我们将37年每年各月关键区域内臭氧纬向偏差求和，以表征该区域臭氧浓度的变化，之后与辐射场以及动力场进行相关分析，并将5月、6月和9月分别取沿70°E、72°E、80°E的垂直剖面分析以探究可能的影响机制。
由图 6可以看到，5月关键区内，100 hPa高度以上(以下)臭氧浓度与长波加热率为显著负(正)相关(图 6a)，相关系数绝对值均超过0.5。而与短波加热率在100 hPa以下为显著负相关(图 6b)，相关系数达到-0.5。但臭氧浓度与总加热率(长波加热率与短波加热率之和)的相关分布与长波加热率的分布类似，即100 hPa高度以上(以下)臭氧浓度与总加热率为显著负(正)相关(图 6c)，强相关中心从低层到高层向北倾斜，相关系数绝对值达到0.5，并通过99%信度水平的显著性检验。因此认为该区域臭氧浓度变化对辐射变化的影响主要集中在长波波段。类似的臭氧浓度与辐射场的相关分布也存在于6月以及9月的关键区(图 7a–c、图 8a–c)。
|
图 6 5月沿70°E垂直剖面上臭氧浓度与(a)长波加热率、(b)短波加热率、(c)总加热率、(d)散度以及(e)垂直速度的相关系数分布。阴影由浅到深分别为通过90%、95%、99%信度水平的显著性检验区域 Figure 6 Distributions of correlation coefficients between ozone concentration and (a) long wave heating rate, (b) short wave heating rate, (c) total heating rate, (d) divergence, and (e) vertical velocity on the vertical section along 70°E in May. Shaded areas from light to dark indicate the values exceed the 90%, 95%, 99% confidence levels, respectively |
|
图 7 同图 6，但为6月沿72°E的相关系数分布 Figure 7 As in Fig. 6, but for distributions of correlation coefficients along 72°E in June |
|
图 8 同图 6，但为9月沿80°E的相关系数分布 Figure 8 As in Fig. 6, but for distributions of correlation coefficients along 80°E in September |
同时，从臭氧浓度的变化与动力场的相关分布看，关键区内5月份臭氧浓度与散度在300~50 hPa上均为负相关(图 6d)，而6月和9月在150 hPa以上(以下)为负(正)相关(图 7d、8d)。至于臭氧浓度与垂直速度的相关关系，3个月份的分布类似，基本均为正相关(图 6e、7e、8e)。
根据以上分析，我们得出臭氧低值影响南亚高压的可能机制：当关键区臭氧浓度为正异常时，辐射加热在100 hPa以上(以下)为负异常(正异常)。在加热条件的变化下，高层(低层)辐散为负异常(正异常)，从而使得上升运动出现负异常。高层辐散以及上升运动的负异常最终使得南亚高压呈现负异常。而臭氧浓度负异常导致南亚高压正异常的机制与之类似。
6 结论与讨论通过分析1979~2015年夏半年5~9月份100 hPa上南亚高压与臭氧纬向偏差的气候特征以及南亚高压强度与其邻近区域臭氧纬向偏差的相互作用，得出如下结论：
除7月外，夏半年南亚高压与臭氧低值的强度年际变化存在因果联系(表 1)。在部分地区，6月和9月，南亚高压和臭氧低值强度变化相互影响，然而在5月和8月二者的因果关系仅仅是单向的。在6月南亚高压和臭氧低值的中部和西部边缘，以及9月南亚高压北部和臭氧低值中心区，臭氧低值增强(减弱)可能和南亚高压增强(减弱)互为因果，相互影响。在6月南亚高压和臭氧低值的东南侧、8月南亚高压和臭氧低值的西部和东部，以及9月南亚高压的西部，南亚高压增强(减弱)可能部分导致了臭氧低值增强(减弱)。在5月南亚高压西部和臭氧低值南部，臭氧低值增强(减弱)可能部分导致了南亚高压增强(减弱)。
|
表 1 5~9月南亚高压与邻近UTLS区臭氧低值相互影响的区域分布特征 Table 1 Distribution characteristics of interaction between South Asia high and ozone low in the adjacent UTLS (Upper Troposphere and Lower Stratosphere) region from May to September |
根据相关分析，推测当南亚高压区臭氧浓度出现正异常，辐射加热在其上部(下部)为负异常(正异常)，导致高层(低层)辐散出现负异常(正异常)，从而使得上升运动出现负异常。高层辐散以及上升运动的负异常最终使南亚高压呈现负异常。而臭氧浓度负异常导致南亚高压出现正异常的机制与之类似。
此外，上升运动异常可能会导致水汽等温室气体异常，其涉及相变的潜热释放过程较为复杂，可能会对结果具有影响，我们会在之后工作中加入水汽的影响，并通过偏相关分析以及数值模拟来进一步验证该机制。另外，臭氧变化对南亚高压强度变化的反馈作用的位置差异等问题还不明确。如6月份，伊朗高原与青藏高原交界地区，南亚高压强度变化与臭氧低值相互作用大致以80°E为界，以东主要为南亚高压对臭氧低值的作用，以西则相反(图 2b、图 3b、图 4b)。以及不同月份南亚高压和臭氧低值因果关系存在很大差异的原因，这些问题都需要进行进一步研究。
致谢: 感谢欧洲中期天气预报中心提供的资料，感谢两位审稿老师提出的建设性意见，感谢编辑老师的辛劳。
Bengtsson L, Roeckner E, Stendel M. 1999. Why is the global warming proceeding much slower than expected? [J]. Geophys. Res., 104(D4): 3865-3876. DOI:10.1029/1998JD200046
|
卞建春, 李维亮, 周秀骥. 1997. 青藏高原及其邻近地区流场结构季节性变化的特征分析[M]//周秀骥. 中国地区大气臭氧变化及其对气候环境的影响(二). 北京: 气象出版社, 257-272. Bian Jianchun, Li Weiliang, Zhou Xiuji. 1997. Analysis of the seasonal variation feature of the wind structure over Tibetan Plateau and its surroundings[M]//Zhou Xiuji. Atmospheric Ozone and Its Impact on Climate and Environment in China (Ⅱ) (in Chinese). Beijing: China Meteorological Press, 257-272.
|
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
|
Brasseur G P, Solomon S. 2005. Aeronomy of the Middle Atmosphere [M]. Netherlands: Springer. DOI:10.1007/1-4020-3824-0
|
陈洪滨, 卞建春, 吕达仁. 2006. 上对流层-下平流层交换过程研究的进展与展望[J]. 大气科学, 30(5): 813-820. Chen Hongbin, Bian Jianchun, Lü Daren. 2006. Advances and prospects in the study of stratosphere-troposphere exchange (in Chinese)[J]. Chinese Journal of Atmospheric Sciences, 30(5): 813-820. DOI:10.3878/j.issn.1006-9895.2006.05.10
|
陈延聪, 王盘兴, 周国华, 等. 2009. 夏季南亚高压的一组环流指数及其初步分析[J]. 大气科学学报, 32(6): 832-838. Chen Yancong, Wang Panxing, Zhou Guohua, et al. 2009. A set of circulation indices for summer South Asian high and their preliminary analyses (in Chinese)[J]. Transactions of Atmospheric Sciences, 32(6): 832-838. DOI:10.3969/j.issn.1674-7097.2009.06.013
|
Forster P M D F, Shine K P. 1999. Stratospheric water vapour changes as a possible contributor to observed stratospheric cooling [J]. Geophys. Res. Lett., 26(21): 3309-3312. DOI:10.1029/1999GL010487
|
Duan A M, Wu G X. 2005. Role of the Tibetan Plateau thermal forcing in the summer climate patterns over subtropical Asia [J]. Climate Dyn., 24(7-8): 793-807. DOI:10.1007/s00382-004-0488-8
|
付超, 李维亮, 周秀骥. 1997. 夏季青藏高原上空臭氧总量低值区形成的模拟实验[M]//周秀骥. 中国地区大气臭氧变化及其对气候环境的影响会议(二). 北京: 气象出版社, 274-285. Fu Chao, Li Weiliang, Zhou Xiuji. 1997. Numerical simulation of the formation of ozone valley over Tibetan Plateau in summer[M]//Zhou Xiuji. Atmospheric Ozone and Its Impact on Climate and Environment in China (Ⅱ) (in Chinese). Beijing: China Meteorological Press, 274-285.
|
郭世昌, 徐裕华. 1986. 大地形加热对亚洲上空臭氧经向分布不连续现象的可能贡献[J]. 云南大学学报, 8(4): 433-439. Guo Shichang, Xu Yuhua. 1986. The possible contribution of the large scale topographic heating to the discontinuous feature of ozonosphere meridionals effect over Asia (in Chinese)[J]. Journal of Yunnan University, 8(4): 433-439.
|
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
|
Guo D, Su Y C, Shi C H, et al. 2015. Double core of ozone valley over the Tibetan Plateau and its possible mechanisms [J]. Journal of Atmospheric and Solar-Terrestrial Physics: 130-131. DOI:10.1016/j.jastp.2015.05.018
|
郭栋, 徐建军, 苏昱丞, 等. 2017. 青藏高原和北美夏季臭氧谷垂直结构和形成机制的比较[J]. 大气科学学报, 40(3): 412-417. Guo Dong, Xu Jianjun, Su Yucheng, et al. 2017. Comparison of vertical structure and formation mechanism of summer ozone valley over the Tibetan Plateau and North America (in Chinese)[J]. Trans. Atmos. Sci., 40(3): 412-417. DOI:10.13878/j.cnki.dqkxxb.20160315001
|
Guo D, Su Y C, Zhou X J, et al. 2017. Evaluation of the trend uncertainty in summer ozone valley over the Tibetan Plateau in three reanalysis datasets [J]. J. Meteor. Res., 31(2): 431-437. DOI:10.1007/s13351-017-6058-x
|
Highwood E J, Hoskins B J. 2010. The tropical tropopause [J]. Quart. J. Roy. Meteor. Soc., 124(549): 1579-1604. DOI:10.1002/qj.49712454911
|
Hoskins B J, Rodwell M J. 1995. A model of the Asian summer monsoon. Part Ⅰ:The global scale [J]. Journal of Atmospheric Sciences, 52(9): 1329-1340. DOI:10.1175/1520-0469(1995)052<1329:AMOTAS>2.0.CO;2
|
Hu D Z, Tian W S, Xie F, et al. 2015. Impacts of stratospheric ozone depletion and recovery on wave propagation in the boreal winter stratosphere [J]. J. Geophys. Res., 120(16): 8299-8317. DOI:10.1002/2014JD022855
|
黄荣辉. 1985. 夏季青藏高原上空热源异常对北半球大气环流异常的作用[J]. 气象学报, 43(2): 208-220. Huang Ronghui. 1985. The influence of the heat source anomaly over Tibetan Plateau on the Northern Hemispheric circulation anomalies (in Chinese)[J]. Acta Meteorologica Sinica, 43(2): 208-220. DOI:10.11676/qxxb1985.026
|
Jin Q, Yang X Q, Sun X G, et al. 2013. East Asian summer monsoon circulation structure controlled by feedback of condensational heating [J]. Climate Dyn., 41(7-8): 1885-1897. DOI:10.1007/s00382-012-1620-9
|
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. Meteor., 2017: 9383495. DOI:10.1155/2017/9383495
|
Liang X S. 2008. Information flow within stochastic dynamical systems [J]. Phys. Rev. E, 78(3 Pt 1): 031113. DOI:10.1103/PhysRevE.78.031113
|
Liang X S. 2013. The Liang-Kleeman information flow:Theory and applications [J]. Entropy, 15(1): 327-360. DOI:10.3390/e15010327
|
Liang X S. 2014. Unraveling the cause-effect relation between time series [J]. Phys. Rev. E, 90(5-1): 052150. DOI:10.1103/PhysRevE.90.052150
|
刘煜, 李维亮. 2001. 青藏高原臭氧低谷的加深及其可能的影响[J]. 气象学报, 59(1): 97-106. Liu Yu, Li Weiliang. 2001. Deepening of ozone valley over Tibetan Plateau and its possible influences (in Chinese)[J]. Acta Meteorologica Sinica, 59(1): 97-106. DOI:10.11676/qxxb2001.011
|
刘屹岷, 吴国雄, 刘辉, 等. 1999. 空间非均匀加热对副热带高压形成和变异的影响——Ⅲ:凝结潜热加热与南亚高压及西太平洋副高[J]. 气象学报, 57(5): 525-538. Liu Yimin, Wu Guoxiong, Liu Hui, et al. 1999. The effect of spatially nonuniform heating on the formation and variation of subtropical high. Part Ⅲ:Condensation heating and South Asia high and western Pacific subtropical high (in Chinese)[J]. Acta Meteorologica Sinica, 57(5): 525-538. DOI:10.11676/qxxb1999.051
|
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]. Advances in Atmospheric Sciences, 20(1): 103-109. DOI:10.1007/BF03342054
|
刘伯奇, 何金海, 王黎娟. 2009. 4~5月南亚高压在中南半岛上空建立过程特征及其可能机制[J]. 大气科学, 33(6): 1319-1332. Liu Boqi, He Jinhai, Wang Lijuan. 2009. Characteristics of the South Asia high establishment processes above the Indo-China Peninsula from April to May and their possible mechanism (in Chinese)[J]. Chinese Journal of Atmospheric Sciences, 33(6): 1319-1332. DOI:10.3878/j.issn.1006-9895.2009.06.17
|
Liu Y, Wang Y, Liu X, et al. 2009. Tibetan middle tropospheric ozone minimum in June discovered from GOME observations [J]. Geophys. Res. Lett., 36(5): L05814. DOI:10.1029/2008GL037056
|
Ramaswamy V, Schwarzkopf M D, Randel W J. 1996. Fingerprint of ozone depletion in the spatial and temporal pattern of recent lower-stratospheric cooling [J]. Nature, 382(6592): 616-618. DOI:10.1038/382616a0
|
Shi C H, Zhang C X, Guo D. 2017. Comparison of electrochemical concentration cell ozonesonde and microwave limb sounder satellite remote sensing ozone profiles for the center of the South Asian high [J]. Remote Sens., 9(10): 1012. DOI:10.3390/rs9101012
|
仕仁睿, 周顺武, 孙绩华, 等. 2017. 青藏高原臭氧亏损变化及其对太阳活动的响应[J]. 云南大学学报(自然科学版), 39(1): 78-87. Shi Renrui, Zhou Shunwu, Sun Jihua, et al. 2017. Variation features of total ozone depletion over the Tibetan Plateau and its responds to solar activity (in Chinese)[J]. Journal of Yunnan University, 39(1): 78-87. DOI:10.7540/j.ynu.20160536
|
苏绍基, 王卫国. 2004. 大气臭氧层变化与南亚高压活动关系的探讨[J]. 云南大学学报(自然科学版), 22(4): 293-296. Su Shaoji, Wang Weiguo. 2004. The relation between the total ozone variation in Asia and the South Asian high (in Chinese)[J]. Journal of Yunnan University, 22(4): 293-296. DOI:10.3321/j.issn:0258-7971.2000.04.013
|
苏昱丞, 郭栋, 郭胜利, 等. 2016. 未来百年夏季青藏高原臭氧变化趋势及可能机制[J]. 大气科学学报, 39(3): 309-317. Su Yucheng, Guo Dong, Guo Shengli, et al. 2016. Ozone trends over the Tibetan Plateau in the next 100 years and their possible mechanism (in Chinese)[J]. Transactions of Atmospheric Sciences, 39(3): 309-317. DOI:10.13878/j.cnki.dqkxxb.20140925002
|
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
|
万凌峰, 郭栋, 刘仁强, 等. 2017. WACCM3对夏季青藏高原臭氧谷的双心结构的模拟性能评估[J]. 高原气象, 36(1): 57-66. Wan Lingfeng, Guo Dong, Liu Renqiang, et al. 2017. Evaluation of WACCM3 performance on simulation of the double core of ozone valley over the Qinghai−Xizang in summer (in Chinese)[J]. Plateau Meteorology, 36(1): 57-66. DOI:10.7522/j.issn.1000-0534.2016.00004
|
吴国雄, 李伟平, 郭华, 等. 1997. 青藏高原感热气泵和亚洲夏季风[M]//叶笃正. 赵九章诞辰九十周年纪念文集. 北京: 科学出版社, 116-120. Wu Guoxiong, Li Weiping, Guo Hua, et al. 1997. Thermal air pumps and the Asian summer monsoon on the Tibetan Plateau[M]//Ye Duzheng. Zhao Jiuzhang's 90th Anniversary Collection (in Chinese). Beijing: Science Press, 116-120.
|
吴国雄, 刘新, 张琼, 等. 2002. 青藏高原抬升加热气候效应研究的新进展[J]. 气候与环境研究, 7(2): 184-201. Wu Guoxiong, Liu Xin, Zhang Qiong, et al. 2002. Progresses in the study of the climate impacts of the elevated heating over the Tibetan Plateau (in Chinese)[J]. Climatic and Environmental Research, 7(2): 184-201. DOI:10.3878/j.issn.1006-9585.2002.02.06
|
徐忠峰, 钱永甫. 2003. 100 hPa高压环流和东风气流的季节、年际和年代际变化[J]. 热带气象学报, 19(3): 225-233. Xu Zhongfeng, Qian Yongfu. 2003. The seasonal, interannual and interdecadal variations of the anticyclone and the easterly current at 100 hPa (in Chinese)[J]. Journal of Tropical Meteorology, 19(3): 225-233. DOI:10.3969/j.issn.1004-4965.2003.03.001
|
Ye D Z. 1981. Some characteristics of the summer circulation over the Qinghai−Xizang (Tibet) Plateau and its neighborhood [J]. Bull. Amer. Meteor. Soc., 62(1): 14-19. DOI:10.1175/1520-0477(1981)062<0014:SCOTSC>2.0.CO;2
|
Yeh T C. 1982. Some aspects of the thermal influences of the Qinghai-Tibetan Plateau on the atmospheric circulation [J]. Archives for Meteorology, Geophysics, and Bioclimatology, Series A, 31(3): 205-220. DOI:10.1007/BF02258032
|
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 B, 66(1): 23415. DOI:10.3402/tellusb.v66.23415
|
Zhang J K, Tian W S, Wang Z W, et al. 2015. The influence of ENSO on northern midlatitude ozone during the winter to spring transition [J]. J. Climate, 28(12): 4774-4793. DOI:10.1175/JCLI-D-14-00615.1
|
张琼, 吴国雄. 2001. 长江流域大范围旱涝与南亚高压的关系[J]. 气象学报, 59(5): 569-577. Zhang Qiong, Wu Guoxiong. 2001. The large area flood and drought over Yangtze River valley and its relation to the South Asia high (in Chinese)[J]. Acta Meteorologica Sinica, 59(5): 569-577. DOI:10.11676/qxxb2001.061
|
张琼, 钱永甫, 张学洪. 2000. 南亚高压的年际和年代际变化[J]. 大气科学, 24(1): 67-78. Zhang Qiong, Qian Yongfu, Zhang Xuehong. 2000. Interannual and interdecadal variations of the South Asia high (in Chinese)[J]. Chinese Journal of Atmospheric Sciences, 24(1): 67-78. DOI:10.3878/j.issn.1006-9895.2000.01.07
|
周任君, 陈月娟. 2005. 青藏高原和伊朗高原上空臭氧变化特征及其与南亚高压的关系[J]. 中国科学技术大学学报, 35(6): 899-908. Zhou Renjun, Chen Yuejuan. 2005. Ozone variations over the Tibetan and Iranian Plateaus and their relationship with the South Asia high (in Chinese)[J]. Journal of University of Science and Technology of China, 35(6): 899-908. DOI:10.3969/j.issn.0253-2778.2005.06.027
|
周利敏, 陈海山, 彭丽霞, 等. 2016. 青藏高原冬春雪深年代际变化与南亚高压可能联系[J]. 高原气象, 35(1): 13-23. Zhou Limin, Chen Haishan, Peng Lixia, et al. 2016. Possible connection between interdecadal variations of snow depth in winter and spring over Qinghai− Xizang Plateau and South Asia high in summer (in Chinese)[J]. Plateau Meteorology, 35(1): 13-23. DOI:10.7522/j.issn.1000-0534.2014.00152
|
Zhou L B, Zou H, Ma S P, et al. 2013. The Tibetan ozone low and its long-term variation during 1979-2010 [J]. Acta Meteorologica Sinica, 27(1): 75-86. DOI:10.1007/s13351-013-0108-9
|
周秀骥, 罗超, 李维亮, 等. 1995. 中国地区臭氧总量变化与青藏高原低值中心[J]. 科学通报, 40(15): 1396-1398. Zhou Xiuji, Luo Chao, Li Weiliang, et al. 1995. The change of total ozone in China and the low value center of the Tibetan Plateau (in Chinese)[J]. Chinese Science Bulletin, 40(15): 1396-1398. DOI:10.3321/j.issn:0023-074X.1995.15.016
|
周秀骥, 李维亮, 陈隆勋, 等. 2004. 青藏高原地区大气臭氧变化的研究[J]. 气象学报, 62(5): 513-527. Zhou Xiuji, Li Weiliang, Chen Longxun, et al. 2004. Study of ozone change over Tibetan Plateau (in Chinese)[J]. Acta Meteorologica Sinica, 62(5): 513-527. DOI:10.11676/qxxb2004.053
|
周秀骥, 赵平, 陈军明, 等. 2009. 青藏高原热力作用对北半球气候影响的研究[J]. 中国科学D辑:地球科学, 52(11): 1473-1486. Zhou Xiuji, Zhao Ping, Chen Junming, et al. 2009. Impacts of thermodynamic processes over the Tibetan Plateau on the Northern Hemispheric climate (in Chinese)[J]. Science in China Series D:Earth Sciences, 52(11): 1679-1693. DOI:10.1007/s11430-009-0194-9
|
Zou H. 1996. Seasonal variation and trends of Toms ozone over Tibet [J]. Geophysical Research Letters, 23(9): 1029-1032. DOI:10.1029/96GL00767
|