气候统计方法和应用：第一课 课程引论+气候跃变.pptx

第一课 课程引论+气候跃变; 通过课堂教学和课外实习，引导学生 理解、运用气候统计分析中的基本概念和常用方法 加深认识统计分析在气候学特别是气候变化研究中的意义;课程大纲;参考书;考核方式; 理论背景 定义/实例 检测方法 研究案例 讨论 ;Climatic Jump – 理论背景;Climatic Jump – 理论背景;混沌、多平衡态、复杂系统 上世纪七、八十年代后上述新概念和理论越来越多地在各种学科领域得到应用和验证。人们也开始重新审视气候变化的基本特征 Lorenz (1963) 混沌吸引子 Charney (1975, 1979) 多平衡态; 气候系统无疑是复杂的，其变量的时间序列则必然是非平稳的、非周期的、时而跃变的; 气候跃变(climatic jump) - 泛指气候从一种状态到另一种状态的较迅速（跳跃性）转变的现象 Yamamoto et al (1986) 时序中的气候跃变 Berger and Labeyrie (1987eds) Abrupt Climate Change，主要反映长期地质气候过程中的有关现象 严中伟等 (1990) 全球要素分布场中的气候跃变 Yan et al (1992) 历史气候跃变的检测 符淙斌、王强 (1992) 突变定义和检测 动力学结构跃变（He, Feng, Wu et al 2008 NPG) ;气候跃变(climatic jump)和突变(abrupt change) 英文文献里多数情况下，abrupt change特指全球变暖过程中，北大西洋温盐环流停顿而导致的相对突然的变冷事件。近年来有人用rapid change似更贴切。;Climatic Jump – 实例;冰后期中的跃变事件 1. Younger Dryas：距今最近的一个冰期于15ky前结束，全球变暖，行将进入全新世；然而在11-12ky之间却突然发生了一次急剧变冷事件（右上图）。这就是著名的“新仙女木” 2. 温暖的全新世期间也发生过类似事件，如距今8.2ky;冰后期跃变事件 机制解释 距今8.2ky前后，全球温暖期期间突然变冷。 右图给出该事件的一种可能机制，即淡水脉冲。模拟过程中第550-570年间在Labrador海域注入适量淡水，结果导致北海经圈翻转流急剧变弱达百年之久，而后缓慢恢复。;;热带季风爆发 - 季节性跃变 5月下旬TBB跃变性地降低，表征大规模的对流云发展; 滑动T-test（均值检验） 滑动F-test（方差检验） 滑动符号检验法 Mann-Kendall’s test 小波分析法 动力学结构分析法;滑动T-test（均值检验） T(i) = (M1 - M2) / SQRT(C s) 其中C = (L1 + L2) / (L1 L2) s = (L1 S12 + L2 S22) / (L1 + L2 -2) 若|T(i)| Tc，说明i时间点附近变化较剧烈； 对所有i计算可得T(i)序列，其中的显著极值点，被定义为跃变点;滑动T-test（均值检验）应用问题 优点：物理意义简明，即两段子序列均值差越大，变率越小，就越可能是一个跃变点 缺点： L1/L2不能自然确定 Note - 气候变化分析中常用距平序列，即为去掉季节cycle，否则长期气候变化很难“显著”;滑动F-test（方差检验） F(i) = S12 / S22 （子序列长L） 若F(i) Fc1 或F(i) Fc2，说明序列的变率大小在 i 时间点附近发生显著的改变 对所有i计算可得F(i)序列，其中显著的极大（小）值点被定义为跃变点 Note - 滑动检验有偶然性，即一个随机序列也可能通过滑动检验找到“跃变点”。如何避免这种偶然性？有统计方法 ;滑动符号检验 设序列X(i)在 i 时间点前后各有子序列（长度L） 比较X(i+k)和X(i-k), p为前者大于后者的情形数，q为前者小于后者的情形数，n = p+q，则 S(i) = min{p,q} 近似服从正态分布（n自由度） 若S(i) Sc，说明 i 时间点前后有显著变化。对所有i计算可得S(i)序列，其极小值点被定义为跃变点。;Mann-Kendall’s test（趋势检验） 设有序列X(i), 定义Kendall系数T = S Ni 其中Ni是X(j)的个数，满足X(i)X(j), ij 在序列无趋势假设下，T渐进逼近正态分布，其均值和方差分别为E = N(N-1)/4; V = N(N-1)(2N+5)/72，