留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

中国夏季降水预测因子潜在技巧分布图及应用

刘伯奇 祝从文

引用本文:
Citation:

中国夏季降水预测因子潜在技巧分布图及应用

Potential Skill Map of Predictors Applied to the Seasonal Forecast of Summer Rainfall in China

  • 摘要: 影响我国夏季汛期降水异常的因子繁多,不同因子之间复杂的相互作用制约我国夏季降水季节预测水平。目前动力模式对降水预测技巧水平较低,如何开发客观统计预报方法,提高我国夏季降水预报技巧依然存在挑战。该文基于最小二乘法拟合和交叉检验方法,提出一种搜索预测因子潜在预测技巧的方法(潜在技巧分布图),并基于该方法开发预测因子自动选择器,建立中国夏季降水异常自动统计预测模型。与传统线性相关分析相比,潜在技巧分布图不受极端气候事件影响,可直观展现具有显著预测技巧的前兆信号,而预测因子自动选择器则能从潜在技巧分布图中自动筛选最优预测因子,获得逐年不同的预测因子,更符合中国夏季降水异常影响因子多样性的客观事实。在完全剔除预测当年信息的回报试验中,该预测模型对1999—2019年中国夏季汛期降水异常的历史回报技巧明显高于动力模式。通过方差订正,历史回报降水的PS评分从71.00分提高到82.10分,显示了该模型的潜在预报潜力。
  • 图 1  基于1981—2019年EOF主模态和多元回归系数重构降水异常场和观测降水异常场(参考态)的泰勒图

    (不同颜色的点表示基于不同数量EOF主模态和多元回归系数的重构结果)

    Fig.1  Taylor diagram of the reconstructed rainfall anomaly field based on EOF modes and multiple regression coefficients referring to the observation during 1981-2019

    (dots in different colors denote reconstructed results using different numbers of EOF modes and principle components)

    图 2  基于1981—2019年前12个EOF主模态和多元回归系数重构的降水异常场与观测逐站降水异常序列时间相关系数空间分布

    (斜线和打点区分别表示达到0.05和0.01显著性水平)

    Fig.2  Correlation coefficient of reconstructed rainfall anomaly using the first 12 EOF modes and multiple regression coefficients to observed rainfall anomaly at each station during 1981-2019

    (areas with significance exceeding 0.05 and 0.01 levels are slashed and stippled, respectively)

    图 3  基于1981—2019年前12个EOF主模态和多元回归系数重构的降水异常场与观测降水异常场空间相关系数逐年时间序列

    Fig.3  Anomaly correlation coefficient of reconstructed rainfall anomaly using the first 12 EOF modes and multiple regression coefficients to observed rainfall anomaly during 1981-2019

    图 4  1999—2019年R2与不同季节预测因子的相关分布

    (打点区表示达到0.05显著水平,预测因子包括30°S~30°N地区降水和南半球、北半球中高纬度地区200 hPa位势高度场)

    Fig.4  Spatial distribution of temporal correlation coefficient between the second regression coefficient(R2) and the predictors in different seasons during 1999-2019

    (the stippled denotes passing the test of 0.05 level, predictors include the rainfall in 30°S-30°N and 200 hPa geopotential height in the mid-high latitude)

    图 5  1999—2019年不同季节预测因子对R2的潜在技巧分布图

    (打点区表示达到0.05显著性水平,预测因子包含30°S~30°N地区降水和南半球、北半球中高纬度地区200 hPa位势高度场)

    Fig.5  Potential skill map of the second regression coefficient(R2) referring to the predictors in different seasons during 1999-2019

    (the stippled denotes passing the test of 0.05 level, predictors include rainfall in 30°S-30°N and 200 hPa geopotential height in the mid-high latitude)

    图 6  预测因子自动选择器提取的1999—2019年前4个多元回归系数的预测因子热度图

    (预测因子包含30°S~30°N地区降水和南半球、北半球中高纬度地区200 hPa位势高度场)

    Fig.6  Heat map of predictors of the first 4 multiple regression coefficients during 1999-2019 obtained by the predictor automatic selection scheme

    (predictors include the rainfall in 30°S-30°N and 200 hPa geopotential height in the mid-high latitude)

    图 7  1999—2019年前4个多元回归系数的回报检验

    (阴影区表示回归系数范围)

    Fig.7  Reforecast test of the first 4 multiple regression coefficients during 1999-2019

    (the shaded denotes range of regression coefficients)

    图 8  采用前12个多元回归系数和EOF模态的回报与观测降水站点相关系数空间分布

    (打点区表示达到0.1显著性水平)

    Fig.8  Correlation coefficient of reforecast anomalous rainfall using the first 12 multiple regression coefficients and EOF modes to observed rainfall anomaly

    (the stippled denotes passing the test of 0.1 level)

    图 9  预测模型对1999—2019年中国夏季汛期平均降水异常的回报检验

    (a)采用前12个多元回归系数和EOF模态的回报结果相对于观测的同号率,(b)基于不同多元回归系数的回报结果相对于观测降水异常的空间相关系数的逐年变化(m表示前1~12个多元回归系数和EOF主模态的回报结果)

    Fig.9  Reforecast test of Chinese summer rainfall anomaly during 1999-2019 using new predicting method

    (a)the same sign rate between reforecast and observation using the first 12 multiple regression coefficients and EOF modes, (b)anomaly correlation coefficients between observation and reforecast based on different numbers of multiple regression coefficients(m, ranging from 1 to 12, indicates the reforecast generated by different numbers of multiple regression coefficients and EOF modes)

    图 10  基于1999—2019年不同多元回归系数的回报结果相对于观测降水异常的空间相关系数多年平均值

    Fig.10  Mean value of anomaly correlation coefficients between observation and reforecast based on different numbers of multiple regression coefficients during 1999-2019

    图 11  回报降水距平百分率的PS评分对方差订正系数的响应函数(a)和方差订正前后的PS评分(b)

    Fig.11  Response PS score curve of the reforecast percentage of rainfall anomaly to the variance correction parameter(a) and PS score before and after variance corrected(b)

  • [1] 李维京.现代气候业务.北京:气象出版社, 2012:202-203.
    [2] Wang B, Xiang B Q, Lee J Y.Subtropical High predictability establishes a promising way for monsoon and tropical storm predictions.PNAS, 2013, 110:2718-2722. doi: 10.1073/pnas.1214626110
    [3] Li T, Wang B, Wu B, et al.Theories on formation of an anomalous anticyclone in western North Pacific during El Nio:A review.J Meteor Res, 2017, 31:987-1006. doi: 10.1007/s13351-017-7147-6
    [4] Wang B, Wu R G, Fu X H.Pacific-East Asian teleconnection:How does ENSO affect East Asian climate?J Climate, 2000, 13:1517-1536. doi: 10.1175/1520-0442(2000)013<1517:PEATHD>2.0.CO;2
    [5] Xie S P, Kosaka Y, Du Y, et al.Indo-Western Pacific Ocean capacitor and coherent climate anomalies in post-ENSO summer:A Review.Adv Atmos Sci, 2016, 33:411-432. doi: 10.1007/s00376-015-5192-6
    [6] Yang J L, Liu Q U, Xie S P, et al.Impact of the Indian Ocean SST basin mode on the Asian summer monsoon.Geophys Res Lett, 2007, 34, DOI:10.1029/2006GL028571.
    [7] 张人禾, 闵庆烨, 苏京志.厄尔尼诺对东亚大气环流和中国降水年际变异的影响:西北太平洋异常反气旋的作用.中国科学(地球科学), 2017, 47(5):544-553. 
    [8] 陆端军, 张先恭.中国降水和温度对ENSO响应的特征.应用气象学报, 1995, 6(1):118-123. 
    [9] 刘永强, 丁一汇.ENSO事件对我国天气气候的影响.应用气象学报, 1992, 3(4):473-481. 
    [10] Feng J, Chen W, Tam C Y, et al.Different impacts of El Nio and El Nio Modoki on China rainfall in the decaying phases.Int J Climatol, 2011, 31:2091-2101. doi: 10.1002/joc.2217
    [11] 王钦, 李双林, 付建建, 等.1998和2010年夏季降水异常成因的对比分析:兼论两类不同厄尔尼诺事件的影响.气象学报, 2012, 70(6):1207-1222. 
    [12] 陈文, 丁硕毅, 冯娟, 等.不同类型ENSO对东亚季风的影响和机理研究进展.大气科学, 2018, 42(3):640-655. 
    [13] Yuan Y, Yang S, Zhang Z.Different evolutions of the Philippine Sea anticyclone between the eastern and central Pacific El Nio:Possible effects of Indian Ocean SST.J Climate, 2012, 25:7867-7883. doi: 10.1175/JCLI-D-12-00004.1
    [14] 袁媛, 杨辉, 李崇银.不同分布型厄尔尼诺事件及对中国次年夏季降水的可能影响.气象学报, 2012, 70(3):467-478. 
    [15] Srinivas G, Chowdary J S, Gnanaseelan C, et al.Impact of differences in the decaying phase of El Nio on South and East Asia summer monsoon in CMIP5 models.Int J Climatol, 2019, 39:5503-5521. doi: 10.1002/joc.6168
    [16] Zhou X, Liu F, Wang B, et al.Different responses of East Asian summer rainfall to El Nio decays.Climate Dyn, 2019, 53:1497-1515. doi: 10.1007/s00382-019-04684-6
    [17] 薛峰, 刘长征.中等强度ENSO对中国东部夏季降水的影响及其与强ENSO的对比分析.科学通报, 2007, 52(23):2798-2805. 
    [18] 袁媛, 高辉, 贾小龙, 等.2014-2016年超强厄尔尼诺事件的气候影响.气象, 2016, 42(5):532-539. 
    [19] 陈桂英.El Nio和La Nia冬季增强型和减弱型及其对中国夏季旱涝的影响.应用气象学报, 2000, 11(2):154-164. 
    [20] Yuan Y, Yang H, Zhou W, et al.Influences of the Indian Ocean dipole on the Asian summer monsoon in the following year.Int J Climatol, 2008, 28:1849-1859. doi: 10.1002/joc.1678
    [21] 肖子牛, 晏红明, 李崇银.印度洋地区异常海温的偶极振荡与中国降水及温度的关系.热带气象学报, 2002, 18(4):335-344. 
    [22] Qiu Y, Cai W, Guo X, et al.The asymmetric influence of the positive and negative IOD events on China's rainfall.Sci Rep, 2014, 4, DOI:10.1038/srep04943.
    [23] Xue F, Wang H, He J H.Interannual variability of Mascarene High and Australian High and their Influences on East Asian summer monsoon.J Meteor Soc Japan, 2004, 82:1173-1186. doi: 10.2151/jmsj.2004.1173
    [24] 贾小龙, 李崇银.南印度洋海温偶极子型振荡及其气候影响.地球物理学报, 2005, 48(6):1238-1249. 
    [25] 杨明珠, 丁一汇.中国夏季降水对南印度洋偶极子的响应研究.大气科学, 2007, 31(4):685-694. 
    [26] Wu Z W, Wang B, Li J P, et al.An empirical seasonal prediction model of the east Asian summer monsoon using ENSO and NAO.J Geophys Res Atmos, 2009, 114, DOI:10.1029/2009JD011733.
    [27] Zuo J Q, Li W J, Sun C H, et al.Impact of the North Atlantic sea surface temperature tripole on the East Asian summer monsoon.Adv Atmos Sci, 2013, 30:1173-1186. doi: 10.1007/s00376-012-2125-5
    [28] 任宏昌, 左金清, 李维京.1998年和2016年北大西洋海温异常对中国夏季降水影响的数值模拟研究.气象学报, 2017, 75(6):877-893. 
    [29] Li J P, Zheng F, Sun C, et al.Pathways of influence of the Northern Hemisphere mid-high latitudes on East Asian climate:A review.Adv Atmos Sci, 2019, 36:902-921. doi: 10.1007/s00376-019-8236-5
    [30] Guo D, Gao Y, Bethke I, et al.Mechanism on how the spring Arctic sea ice impacts the East Asian summer monsoon.Theor Appl Climatol, 2014, 115:107-119. doi: 10.1007/s00704-013-0872-6
    [31] Wu B Y, Zhang R H, Wang B, et al.On the association between spring Arctic sea ice concentration and Chinese summer rainfall.Geophys Res Lett, 2009, 36, DOI:10.1029/2009GL037299.
    [32] Liu X, Yanai M.Influence of Eurasian spring snow cover on Asian summer rainfall.Int J Climatol, 2002, 22:1075-1089. doi: 10.1002/joc.784
    [33] Yim S Y, Jhun J G, Lu R, et al.Two distinct patterns of spring Eurasian snow cover anomaly and their impacts on the East Asian summer monsoon.J Geophys Res Atmos, 2010, 115, DOI:10.1029/2010JD013996.
    [34] 张人禾, 张若楠, 左志燕.中国冬季积雪特征及欧亚大陆积雪对中国气候影响.应用气象学报, 2016, 27(5):513-526. 
    [35] Xiao Z X, Duan A M.Impacts of Tibetan Plateau snow cover on the interannual variability of the East Asian summer monsoon.J Climate, 2016, 29:8495-8514. doi: 10.1175/JCLI-D-16-0029.1
    [36] Kim J E, Hong S-Y.Impact of soil moisture anomalies on summer rainfall over East Asia:A regional climate model study.J Climate, 2007, 20:5732-5743. doi: 10.1175/2006JCLI1358.1
    [37] Zuo Z Y, Zhang R H.The spring soil moisture and the summer rainfall in eastern China.Chinese Sci Bull, 2007, 52:3310-3312. doi: 10.1007/s11434-007-0442-3
    [38] Zuo Z Y, Zhang R H.Influence of soil moisture in eastern China on the East Asian summer monsoon.Adv Atmos Sci, 2016, 33:151-163. doi: 10.1007/s00376-015-5024-8
    [39] Wang B, Lee J Y, Kang I S, et al.How accurately do coupled climate models predict the leading modes of Asian-Australian monsoon interannual variability?Climate Dyn, 2008, 30:605-619. doi: 10.1007/s00382-007-0310-5
    [40] Wang B, Lee J Y, Kang I S, et al.Advance and prospectus of seasonal prediction:Assessment of the APCC/CliPAS 14-model ensemble retrospective seasonal prediction (1980-2004).Climate Dyn, 2008, 33:93-117. 
    [41] Zhang W J, Jin F F, Stuecker M F, et al.Unraveling El Nio's impact on the East Asian monsoon and Yangtze River summer flooding.Geophys Res Lett, 2016, 43:11375-11382. doi: 10.1002/2016GL071190
    [42] 高辉, 高晶.黑潮冬季海温对我国东北地区夏季降水预测信号的增强.海洋学报, 2014, 36(7):27-33. 
    [43] 冯蕾, 魏凤英, 朱艳峰.基于前春对流层温度和北大西洋涛动的中国夏季降水统计预测模型.大气科学, 2011, 35(5):963-976. 
    [44] 刘颖, 任宏利, 张培群, 等.利用高原积雪信号改进我国南方夏季降水预测的新方法及其在2014年降水预测中的应用试验.大气科学, 2017, 41(2):313-320. 
    [45] 刘婷婷, 陈海山, 蒋薇, 等.基于土壤湿度和年际增量方法的我国夏季降水预测试验.大气科学, 2016, 40(3):591-603. 
    [46] Xing W, Wang B, Yim S Y.Long-lead seasonal prediction of China summer rainfall using an EOF-PLS regression-based methodology.J Climate, 2016, 29:1783-1796. doi: 10.1175/JCLI-D-15-0016.1
    [47] 陈丽娟, 顾伟宗, 伯忠凯, 等.黄淮地区夏季降水的统计降尺度预测.应用气象学报, 2017, 28(2):129-141. 
    [48] 魏凤英.气候统计诊断与预测方法研究进展——纪念中国气象科学研究院成立50周年.应用气象学报, 2006, 17(6):736-742. 
    [49] 高辉, 王永光.ENSO对中国夏季降水可预测性变化的研究.气象学报, 2007, 65(1):131-137. 
    [50] Xie P, Arkin P A.Global precipitation:A 17-year monthly analysis based on gauge observations, satellite estimates, and numerical model outputs.Bull Amer Meteor Soc, 1997, 78:2539-2558. doi: 10.1175/1520-0477(1997)078<2539:GPAYMA>2.0.CO;2
    [51] Kobayashi S, Ota Y, Harada Y, et al.The JRA-55 Reanalysis:General specifications and basic characteristics.J Meteor Soc Japan, 2015, 93:5-48. doi: 10.2151/jmsj.2015-001
    [52] 庞轶舒, 祝从文, 刘凯.中国夏季降水异常EOF模态的时间稳定性分析.大气科学, 2014, 38(6):1137-1146. 
    [53] Michaelsen J.Cross-validation in statistical climate forecast models.J Climate Appl Meteor, 1987, 26:1589-1600. doi: 10.1175/1520-0450(1987)026<1589:CVISCF>2.0.CO;2
    [54] Yim S Y, Wang B, Xing W.Peak-summer East Asian rainfall predictability and prediction part Ⅱ:Extratropical East Asia.Climate Dyn, 2016, 47:15-30. doi: 10.1007/s00382-015-2849-x
    [55] Xing W, Wang B, Yim S Y, et al.Predictable patterns of the May-June rainfall anomaly over East Asia.J Geophys Res Atmos, 2017, 122:2203-2217. doi: 10.1002/2016JD025856
    [56] Trenberth K E, Shea D J.Atlantic hurricanes and natural variability in 2005.Geophys Res Lett, 2006, 33, DOI:10.1029/2006GL026894.
    [57] Newman M, Alexander M A, Ault T R, et al.The Pacific decadal oscillation, revisited.J Climate.2016, 29:4399-4427. doi: 10.1175/JCLI-D-15-0508.1
    [58] Zhu C W, Park C K, Lee W S, et al.Statistical downscaling for multi-model ensemble prediction of summer monsoon rainfall in the Asia-Pacific region using geopotential height field.Adv Atmos Sci, 2008, 25:867-884. doi: 10.1007/s00376-008-0867-x
    [59] 刘长征, 杜良敏, 柯宗建, 等.国家气候中心多模式解释应用集成预测.应用气象学报, 2013, 24(6):677-685. 
    [60] 任宏利.动力季节预测中预报误差与物理因子的关系.应用气象学报, 2008, 19(3):276-286. 
    [61] 封国林, 赵俊虎, 支蓉, 等.动力-统计客观定量化汛期降水预测研究新进展.应用气象学报, 2013, 24(6):656-665. 
  • [1] 谭桂容,  段浩,  任宏利.  中高纬度地区500 hPa高度场动力预测统计订正 . 应用气象学报, 2012, 23(3): 304-311.
    [2] 李春晖,  梁建茵.  基于Shepard和OI方法对雨量计逐时资料的分析 . 应用气象学报, 2010, 21(4): 416-422.
    [3] 谢炯光.  广东省前汛期降水分布趋势的分析与预测 . 应用气象学报, 1997, 8(4): 452-459.
    [4] 范可,  王会军.  南极涛动异常与2006年我国东部夏季降水形势预测 . 应用气象学报, 2006, 17(3): 383-384.
    [5] 赵俊虎,  杨杰,  封国林,  张世轩.  2011年我国夏季降水动力统计预测与异常成因 . 应用气象学报, 2013, 24(1): 43-54.
    [6] 刘娜,  李双林.  基于时间尺度分离的中国东部夏季降水预测 . 应用气象学报, 2015, 26(3): 328-337. DOI: 10.11898/1001-7313.20150308
    [7] 何敏,  李小泉.  热带环流异常与我国夏季降水分布的关系 . 应用气象学报, 1992, 3(2): 181-189.
    [8] 尤凤春,  丁裕国,  周煜,  史印山.  奇异值分解和奇异交叉谱分析方法在华北夏季降水诊断中的应用 . 应用气象学报, 2003, 14(2): 176-187.
    [9] 孙颖,  丁一汇.  1997年东亚夏季风异常活动在汛期降水中的作用 . 应用气象学报, 2002, 13(3): 277-287.
    [10] 曾晓青,  邵明轩,  王式功,  刘还珠.  基于交叉验证技术的KNN方法在降水预报中的试验 . 应用气象学报, 2008, 19(4): 471-478.
    [11] 封国林,  赵俊虎,  支蓉,  龚志强,  郑志海,  杨杰,  熊开国.  动力-统计客观定量化汛期降水预测研究新进展 . 应用气象学报, 2013, 24(6): 656-665.
    [12] 任宏利.  动力季节预测中预报误差与物理因子的关系 . 应用气象学报, 2008, 19(3): 276-286.
    [13] 陈丽娟,  顾伟宗,  伯忠凯,  刘向文.  黄淮地区夏季降水的统计降尺度预测 . 应用气象学报, 2017, 28(2): 129-141. DOI: 10.11898/1001-7313.20170201
    [14] 陈官军,  魏凤英,  姚文清,  周璇.  基于南海夏季风季节内振荡的降水延伸预报试验 . 应用气象学报, 2016, 27(3): 273-284. DOI: 10.11898/1001-7313.20160302
    [15] 覃志年,  陈丽娟,  唐红玉,  黄颖.  月尺度动力模式产品解释应用系统及预测技巧 . 应用气象学报, 2010, 21(5): 614-620.
    [16] 祝从文,  刘伯奇,  左志燕,  袁乃明,  刘舸.  东亚夏季风次季节变化研究进展 . 应用气象学报, 2019, 30(4): 401-415. DOI: 10.11898/1001-7313.20190402
    [17] 毕硕本,  徐寅,  覃志年,  陈最,  王必强.  EMD在广西季节降水预报中的应用 . 应用气象学报, 2010, 21(3): 366-371.
    [18] 张人禾,  殷永红,  李清泉,  刘益民,  牛涛.  利用ARGO资料改进ENSO和我国夏季降水气候预测 . 应用气象学报, 2006, 17(5): 538-547.
    [19] 肖科丽,  赵国令,  方建刚,  孙娴.  影响陕西夏季降水主要因子及增量预测方法 . 应用气象学报, 2017, 28(4): 458-469. DOI: 10.11898/1001-7313.20170407
    [20] 吴洪宝,  王盘兴,  林开平.  广西夏季降水量潜在可预报性估计 . 应用气象学报, 2005, 16(4): 445-452.
  • 加载中
图(11)
计量
  • 文章访问数:  1012
  • HTML全文浏览量:  146
  • PDF下载量:  17
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-04-25
  • 修回日期:  2020-06-03
  • 刊出日期:  2020-09-30

中国夏季降水预测因子潜在技巧分布图及应用

  • 中国气象科学研究院, 北京 100081

摘要: 影响我国夏季汛期降水异常的因子繁多,不同因子之间复杂的相互作用制约我国夏季降水季节预测水平。目前动力模式对降水预测技巧水平较低,如何开发客观统计预报方法,提高我国夏季降水预报技巧依然存在挑战。该文基于最小二乘法拟合和交叉检验方法,提出一种搜索预测因子潜在预测技巧的方法(潜在技巧分布图),并基于该方法开发预测因子自动选择器,建立中国夏季降水异常自动统计预测模型。与传统线性相关分析相比,潜在技巧分布图不受极端气候事件影响,可直观展现具有显著预测技巧的前兆信号,而预测因子自动选择器则能从潜在技巧分布图中自动筛选最优预测因子,获得逐年不同的预测因子,更符合中国夏季降水异常影响因子多样性的客观事实。在完全剔除预测当年信息的回报试验中,该预测模型对1999—2019年中国夏季汛期降水异常的历史回报技巧明显高于动力模式。通过方差订正,历史回报降水的PS评分从71.00分提高到82.10分,显示了该模型的潜在预报潜力。

English Abstract

    • 东亚夏季风异常是造成我国夏季汛期降水异常的关键因素,中国夏季汛期降水异常的季节预测既是中国气象局的核心业务之一,也是国家防灾减灾工作的第一道防线。东亚夏季风包含南海—西北太平洋热带季风系统和东亚副热带季风系统,其对应的降水期包括江南春雨、江淮梅雨、华北—东北雨季和华南后汛期,由于降水期延续时间长、空间跨度广、影响因素多、预测难度大,夏季汛期降水异常的季节预测一直是东亚季风研究的重点和难点。已有研究表明:季节-年际预测的可预报性来源于大气下垫面外强迫[1],包括海温、海冰、积雪和土壤湿度等陆面过程。国内外学者针对下垫面外强迫前兆信号如何影响东亚夏季风降水已开展大量研究,取得了丰富成果,如热带ENSO(El Niño-Southern Oscillation)事件能够通过西北太平洋局地海-气相互作用[2-4]或热带印度洋全区一致型海温异常的“电容器”效应[5-6]影响夏季对流层低层的西北太平洋副热带高压[7](以下简称副高),进而影响我国夏季汛期降水的空间分布[8-9]。而ENSO事件的多样性,如中部型、东部型和混合型事件[10-14]、ENSO事件衰减快慢[15-16]、超强El Niño事件[17-18]等,使其对我国夏季汛期降水异常的影响变得更加复杂[19]。同时,热带印度洋偶极型海温异常[20-22]和南印度洋偶极子型海温异常[23-25]也能通过调控Walker环流和越赤道气流,改变东亚夏季风水汽输送和副高位置及强度,进而引起我国夏季汛期降水异常。此外,北大西洋的三极子型海温异常及其相关的北大西洋涛动(NAO)可激发出跨越欧亚大陆的准正压纬向波列影响我国夏季汛期降水的分布特征[26-29]。就海冰强迫而言,春季北极海冰通过改变北太平洋海温异常[30]或欧亚大陆积雪[31],影响东亚夏季风环流和我国夏季汛期降水异常。就陆面过程而言,冬季、春季欧亚大陆[32-34]和青藏高原西部积雪异常[35]是我国夏季汛期降水异常的重要前兆信号,而东亚地区春季土壤湿度通过其记忆性影响东亚地区纬向海陆热力对比,从而影响东亚夏季风环流和降水[36-38]

      目前,数值模式已从单独大气环流模式发展成为海-陆-气-冰耦合的气候系统模式,基于模式产品的动力预测逐渐成为提高汛期降水季节预测能力的重要手段。尽管动力模式对热带大气环流和气温降水的预测技巧大幅提高,但由于模式尚无法完全准确模拟ENSO演变及其遥相关和海-陆-气相互作用[39],因此模式对热带外地区,尤其是对我国夏季汛期降水的预测技巧仍十分有限[40]。考虑到动力预测能力的不足,基于大气外强迫影响东亚夏季风的机理研究,很多学者从不同角度建立了我国夏季汛期降水异常的统计预测模型。这些模型有的基于ENSO对副高的调控[41],有的基于黑潮延伸体海温异常的气候效应[42],有的基于北大西洋NAO的下游影响[43],有的基于欧亚大陆积雪[44]或土壤湿度[45],有的则基于上述多种因子[46]

      较动力预测,这些经验统计预测模型具有更高的预测技巧[47-48]。传统的统计预测建模大多基于线性相关、回归或奇异值分解(SVD)方法确定预测因子,通过独立样本检验反映模型的预测技巧。但这种建模流程在确定预测因子时不可避免地受到极端气候事件的影响(如超强El Niño事件),而独立样本检验又假定某一因子能长期稳定地影响预测对象。实际业务工作证明,每年影响我国夏季汛期降水异常的关键前兆信号不同,且同一预测因子(如ENSO)的气候效应也随时间变化[49]。为了克服传统统计预测模型的不足,本文提出基于交叉检验寻找预测因子的方案——潜在技巧分布图,设计用于确定逐年最优预测因子的因子自动选择器,建立中国夏季汛期降水异常的客观预测方法。

    • 本文预测对象是中国160站夏季(6—8月)平均降水异常,参考气候态为1981—2010年平均。大气环流的异常演变包含对下垫面强迫响应的重要信息,是大气外强迫因子与预测对象之间的重要桥梁。为兼顾热带和热带外地区对流和环流对外强迫响应的差异性,本文选取的预测因子包括低纬度地区(30°S~30°N)美国CPC的CMAP(Climate Prediction Center Merged Analysis of Precipitation)逐月降水[50]和南半球、北半球中高纬度地区(30°S以南和30°N以北)日本JRA-55再分析产品的逐月200 hPa位势高度场[51]。预测对象和预测因子的研究时段均为1981—2019年,从前期秋季(9—11月)开始至前期冬季(12月—次年2月)结束,针对预测因子进行3个月滑动平均,分别从前期4个不同阶段(9—11月、10—12月、11月—次年1月、12月—次年2月)对夏季降水预测。

    • 假设有N年降水异常资料Y,预测当年为第i年,其降水异常场为Yi,首先用剔除预测当年信息的N-1年资料进行EOF展开,得到m个EOF主模态(Vki,1≤km),再将第i年的观测资料投影到这些主模态上,获得第i年的多元回归系数(R1, i, R2, i, R3, i, …, Rm, i), 则Yi可以展开为

      (1)

      其中,Rki表示Yi对第k个EOF主模态(Vki)的回归系数(即预测对象)。与传统的EOF分析相比,本方法中逐年的EOF模态(Vki)相对独立,虽然Vki会随预测年份变化,但由于本方法针对每个预测年独立计算其多元回归系数Rki,因此无论Vki如何变化,Yi总能通过RkiVki正确地重构。本方法将EOF模态视为时变函数,并使Rki随之逐年调整,可避免将预测当年信息通过主成分带入交叉检验的回报模型中。

      包含预测年信息时,取前4个EOF主模态和主成分能够较全面刻画我国夏季汛期降水异常的整体时空变化特征[52]。剔除预测年信息时,则需要更多模态和多元回归系数方能反映我国夏季汛期降水异常的年际变化。分析发现,随着模态数的增加,重构场和观测场的空间相关系数逐渐升高,同时两者的空间方差比(重构场的空间方差除以观测场的空间方差)也逐渐加大(图 1)。在该过程中,空间相关系数和空间方差比的增幅随模态为2个增至12个的过程中迅速加大,随后逐渐减小。这种非线性增幅说明,截取前12个EOF模态和多元回归系数(R1, R2, R3, …, R12)不会对预报模型和预报结果造成不稳定影响。基于前12个EOF模态和多元回归系数的重构场能很好地重现我国东部黄河以南地区夏季汛期降水异常的年际波动(图 2),重构场和观测场的平均空间相关系数为0.601,其年际变化范围为0.35~0.85(图 3)。值得注意的是,基于交叉检验和空间投影的EOF主模态和多元回归系数重构场对我国华北—东北和西部地区夏季汛期降水年际变率的刻画能力较弱,这是因为上述地区夏季降水以短时强降水过程为主,故在季节平均降水异常场上信号较弱,当采用空间投影法获取多元回归系数时,很可能无法准确捕获这些地区的夏季降水异常。这与我国西部和北方夏季降水的季节可预报性低于东部和南方的传统认知相符。此外,重构场和观测场相比空间方差较小(图 1),这与EOF分析方法的空间平滑特性有关,也说明有必要对预测结果进行合理的方差订正。

      图  1  基于1981—2019年EOF主模态和多元回归系数重构降水异常场和观测降水异常场(参考态)的泰勒图

      Figure 1.  Taylor diagram of the reconstructed rainfall anomaly field based on EOF modes and multiple regression coefficients referring to the observation during 1981-2019

      图  2  基于1981—2019年前12个EOF主模态和多元回归系数重构的降水异常场与观测逐站降水异常序列时间相关系数空间分布

      Figure 2.  Correlation coefficient of reconstructed rainfall anomaly using the first 12 EOF modes and multiple regression coefficients to observed rainfall anomaly at each station during 1981-2019

      图  3  基于1981—2019年前12个EOF主模态和多元回归系数重构的降水异常场与观测降水异常场空间相关系数逐年时间序列

      Figure 3.  Anomaly correlation coefficient of reconstructed rainfall anomaly using the first 12 EOF modes and multiple regression coefficients to observed rainfall anomaly during 1981-2019

    • 通过EOF分析,将预测对象转换为12个多元回归系数序列,将预测的多元回归系数和剔除预测当年信息的EOF空间模态进行重构,即可得到预测当年降水场。因此,如何准确预测每一个回归系数成为关键问题。本章将结合降水场多元回归系数的预测问题,介绍搜索预测因子潜在预测技巧的方案和客观自动选择预测因子的思路。

    • 传统的经验统计预测往往采用相关系数分布或SVD分析方法确定预测模型,但这与实际的气候预测存在较大差异。首先,相关系数反映的预测因子是预测对象的必要条件,而实际预测却需要寻找预测对象的充分条件;其次,相关系数易受极端气候事件的影响,且预测对象和预测因子之间的相关性还易受气候变化和年代际变率等缓变过程影响。因此,基于相关分析得到的预测因子随时间变化不稳定,直接导致预测模型的预测技巧不稳定,表现出强烈的年际和年代际波动。因此,基于交叉检验和最小二乘法线性拟合设计预测因子的潜在技巧分布图。

      具体计算步骤如下:假设预测对象多元回归系数R(mN)和预测因子场X(x, y, N)的总样本量为N,预测当年为第i年,其中m表示多元回归系数个数,xy分别表示预测因子场的经向和纬向格点数,那么,①从R(mN)和X(x, y, N)中剔除掉第i-1, ii+1的数据,形成训练数据集。连续剔除预测当年及其相邻两年的数据是为了在年际尺度上完全去除预测当年信息,以尽量减小过度拟合的影响。该方法被广泛应用于季节预测回报试验的交叉检验中[53-55]。②基于训练数据集,利用最小二乘法拟合,建立预测因子场X(x, y, N)每个格点数据序列和预测对象的线性函数关系,得到回归参数K和截距B;再将第i年的前1年秋季至冬季的预测因子格点数据代入以上函数关系,得到第i年夏季降水场回归系数的预测值;遍历X的所有格点(x, y),得到基于不同格点数据的第i年预测回归系数场。③将预测年i循环N次,得到N年基于预测因子场X每个格点数据序列的预测回归系数场;再计算观测回归系数序列与预测回归系数场的相关系数,将其中具有显著相关性的格点视为存在潜在预测技巧,故将预测回归系数场和观测回归系数序列之间的相关系数场定义为潜在技巧分布图(potential skill map,简称PSM)。

      与传统相关分析相比,PSM采用实际统计预测建模的计算流程,能够准确提取预测对象的充分条件,同时由于采用交叉检验思路,其结果可消除极端气候事件影响,因此得到的统计关系具有更好的时间稳定性。需要说明的是,在使用PSM时,需要人为确定总样本N的分析时段。本文采用的分析时段是1999—2019年,这是因为气候系统自然变率在20世纪90年代末期发生了显著的年代际调整,表现为北大西洋多年代振荡(AMO)[56]和太平洋年代际振荡(PDO)相继发生位相转换[57],受其影响,东亚夏季风的影响因子发生明显改变,这说明对统计预测而言,数据长度并非越长越好,应根据气候变化的观测事实挑选合理的训练期。

      R2为例,图 4图 5展示传统相关系数分布图和PSM的区别。以R2为例,是因为它的前兆信号范围相对较广,更容易体现二者差别。事实上,其他回归系数也有类似的差别(图略)。虽然与R2显著相关的低纬度对流出现在前期秋季、冬季的热带西南印度洋(图 4),但PSM显示,上述低纬度对流对R2并无潜在预测技巧(图 5)。同时,前期秋季、冬季与R2显著相关的中高纬度环流异常位于欧亚大陆、北美大陆和南极大陆上空(图 4),这些中高纬度环流信号在PSM上也显示出对R2的显著预测技巧(图 5)。在当年春季和夏季,与R2显著相关的大范围低纬度对流位于热带东太平洋(图 4),但它们在PSM上却并不明显(图 5)。夏季同期具有潜在技巧的预测因子则是欧洲北部的高空环流(图 4图 5)。同期预测因子为建立动力-统计相结合预测模型提供基础。尽管动力模式对中国夏季降水异常的直接预测技巧偏低,但其对低纬度降水和中高纬度高空环流却具有较好预测技巧。由于同期预测因子影响预测对象的物理过程更直接、不确定性更小,因此,可应用具有动力预测技巧的夏季热带外高空环流与观测多元回归系数的PSM,提高动力模式对中国夏季降水异常的预测技巧。以上对比说明,PSM比传统相关系数分布图更加严格,具有潜在技巧的前兆信号一定和预测因子具有显著相关,但和预测因子显著相关的前兆信号却未必具有潜在预测技巧。

      图  4  1999—2019年R2与不同季节预测因子的相关分布

      Figure 4.  Spatial distribution of temporal correlation coefficient between the second regression coefficient(R2) and the predictors in different seasons during 1999-2019

      图  5  1999—2019年不同季节预测因子对R2的潜在技巧分布图

      Figure 5.  Potential skill map of the second regression coefficient(R2) referring to the predictors in different seasons during 1999-2019

    • 尽管PSM能过滤掉和预测对象显著相关却没有预测技巧的前兆信号,但仍需从中确定每年的关键预测因子,以便准确开展统计预测。传统的统计预测中,挑选关键因子依赖于预报员或科研人员的经验,具有很强的主观性。为克服这一不足,在PSM的基础上,借鉴集合预报思想,设计预测因子自动选择器,其计算方法分为两步:①从具有潜在预测技巧的预测因子格点中,挑选出预测和观测回归系数符号一致率达到一定阈值的格点,进入预测因子自动选择器。对1999—2019年共21年数据,选择的符号一致率样本量阈值为17,即21年的预测中,有17年预测结果和观测值同号。②将达到符号一致率阈值的所有预测回归系数值的中位数定义为回归系数的预测值。因此,在最终建模时,仅基于单一格点上的预测因子给出回归系数的预测值,这样可有效地避免多因子建模时的过度拟合问题。

      预测因子自动选择器的优点之一是完全客观化,挑选的预测因子不会因人而异;它的另一优点是能针对不同年份挑选出不同预测因子,更符合实际情况。热度图能够反映这些多元回归系数的关键因子,且不同回归系数的关键因子存在明显年际差异(图 6)。每幅热度图上都有21个点,对应为在21年的回报试验中,每个回归系数逐年不同的关键因子,若某些区域内点的聚集程度越高,则表示这些区域内的预测因子在回报试验中所起作用越大。以前4个多元回归系数为例,R1的关键因子主要集中在10—12月热带西太平洋暖池对流和12月—次年1月南太平洋的高空环流R2的关键因子则更多地分布在前期秋季、冬季南半球中高纬度高空环流中,对R3而言,前期秋季、冬季的大西洋低纬度对流成为大多数年份的关键因子,而R4关键因子的空间分布则较为分散,说明其预测难度相对较大。对其余多元回归系数而言,R5R7的关键因子主要是前期秋季和初冬欧亚大陆上空的高空环流,R9R10的关键因子主要是11月—次年1月和12月—次年2月南极洲上空的高空环流,R12的关键因子则包括前年10—12月北太平洋低纬度对流和12月—次年1月南太平洋高空环流,而R8和R11关键因子的空间分布较分散(图略)。

      图  6  预测因子自动选择器提取的1999—2019年前4个多元回归系数的预测因子热度图

      Figure 6.  Heat map of predictors of the first 4 multiple regression coefficients during 1999-2019 obtained by the predictor automatic selection scheme

      预测因子自动选择器挑选的关键因子对预测对象的影响应具有较清晰的物理过程。需要指出的是,预测对象和预测因子之间时空尺度的一致性是考察预测模型物理属性的重要前提。本方法将大范围的预测对象和小范围的预测因子建立统计关系,后者对前者的影响可能有两种途径:①直接影响,即小范围预测因子通过某种“升尺度”物理过程直接影响中国大范围汛期降水异常型。以R1的关键影响因子为例,10—12月热带西太平洋暖池对流反映前冬ENSO事件通过暖池区热带对流影响我国南方地区夏季降水异常的过程,在这一过程中,海-气相互作用和热带-热带外遥相关型将西太平洋暖池的局地对流异常信号和影响整个东亚地区的东亚—太平洋遥相关波列相联系。②间接影响,即小范围预测因子代表某种大尺度信号的统计投影,在这种情况下,需要分析预测因子对应的大尺度信号,并进一步分析其影响东亚夏季降水异常的物理过程。图 6揭示的某些关键因子尚无明确的机理解释,这说明对东亚夏季风和我国夏季汛期降水年际变率的理解尚不全面。因此,潜在技巧图也为深入研究东亚夏季风年际变率机理提供新的切入点。

    • PSM和预测因子自动选择器可产生逐年降水场多元回归系数的历史回报值及其预测技巧。初筛后的每个预测因子均产生1个回归系数预测值,对应图 7的阴影区。以前4个多元回归系数为例,通过预测因子自动选择器第1步初筛的预测因子个数依次为115,312,157和410。回报的多元回归系数值基本覆盖观测回归系数的变化范围,说明自动挑选的因子能反映回归系数变化的多样性。多元回归系数预测值的中位数集合就是各年回报的回归系数,对应图 7的红线,且每个回归系数回报所用的预测因子具有明显年际差异(图 6)。同时,前4个多元回归系数回报结果和观测的相关系数依次为0.70,0.57,0.82和0.87,均达到0.01的显著性水平。对第5至第12个多元回归系数而言,回报结果和观测数据的相关系数依次为0.8,0.7,0.74,0.85,0.47,0.59,0.8和0.76,均达到0.01显著性水平(图略)。以上结果说明该方法对多元回归系数具有显著回报技巧。

      图  7  1999—2019年前4个多元回归系数的回报检验

      Figure 7.  Reforecast test of the first 4 multiple regression coefficients during 1999-2019

      基于每年预测的多元回归系数和剔除该年信息的EOF主模态,能重构出预测年夏季汛期降水异常值,通过和观测降水异常值对比反映该方法的整体预测技巧。当采用前12个多元回归系数和EOF模态进行回报时,该预测模型对我国东部大部分地区的夏季汛期降水异常具有较高的季节预测技巧(提前3个月,基于前期秋季至冬季的预测因子),大值中心位于长江流域、黄河流域和珠江流域等人口密集区(图 8),这与预测技巧上限(图 2)的分布特征一致。同时,该模型对我国华北—东北地区和西部地区夏季降水的季节预测能力相对较弱。从降水型的整体分布看,该模型预测的我国夏季汛期降水异常场和观测结果的符号一致率变化范围为40%~70%,平均符号一致率达到60%(图 9a),两者空间相关系数的变化范围为0.25~0.65,平均空间相关系数为0.436(图 9b),均远高于现有主流的气候模式动力预测技巧(空间相关系数小于0.1)[40, 46, 58]。当采用前12个多元回归系数和EOF空间模态进行回报试验时,回报结果和观测的符号一致率及空间相关系数的年际波动相对较小,说明该模型的预测技巧具有良好的时间平稳性。

      图  8  采用前12个多元回归系数和EOF模态的回报与观测降水站点相关系数空间分布

      Figure 8.  Correlation coefficient of reforecast anomalous rainfall using the first 12 multiple regression coefficients and EOF modes to observed rainfall anomaly

      图  9  预测模型对1999—2019年中国夏季汛期平均降水异常的回报检验

      Figure 9.  Reforecast test of Chinese summer rainfall anomaly during 1999-2019 using new predicting method

      需要指出的是,多元回归系数的物理含义是降水场在由某组EOF向量作为基底所支撑的相空间中各个分量轴上的投影,因此,需要考虑EOF模态作为相空间基底的时间稳定性对预测性能的影响。为定量分析这种影响,可分别取不同的多元回归系数和EOF模态(m=1, 2, 3, …, 12)进行回报检验(图 9b图 10)。结果表明:除了2001年,其余年份回报结果相对于观测降水的空间相关系数均随着模态数的增加而增大(图 9b)。多年平均结果也证明,高阶模态的引入能够进一步提高回报技巧(图 10)。这说明本文中的高阶模态仍具有潜在的物理意义,它们通过增加回报结果的多样性提高预测技巧。此外,空间相关系数随m值的增长并非线性,在m为1~4阶段,空间相关系数增长最快,随后增长变缓,这说明前4个多元回归系数和主模态对降水的整体分布起决定性作用,这与庞轶舒等[52]的结果一致。若基于PSM中具有显著技巧的某一区域大范围平均的前兆信号进行我国夏季降水异常的回报试验,回报效果(1999—2019年平均空间相关系数为0.15~0.23,依赖于选取的预测因子)远低于自动因子选择器,这说明变因子预测更符合实际情况,即每年影响我国夏季降水异常的关键因子可能不同。

      图  10  基于1999—2019年不同多元回归系数的回报结果相对于观测降水异常的空间相关系数多年平均值

      Figure 10.  Mean value of anomaly correlation coefficients between observation and reforecast based on different numbers of multiple regression coefficients during 1999-2019

      图 11是将历史回报的降水异常转换为降水距平百分率后,基于国家气候中心预测技巧评分标准(PS评分)的回报技巧。结果表明:在回报时段内,该模型预测的我国夏季汛期降水距平百分率平均PS评分为71.00分,除2009年PS评分偏低以外,其余年份的PS评分均在70分左右,回报结果PS评分的时间平稳性良好(图 11b)。如前文所述,该模型在采用EOF分析和统计建模的过程中,不可避免地使预测降水的方差衰减。为合理定义方差订正系数,需要首先计算预测结果PS评分对方差订正系数A(线性放大倍数)的响应曲线(图 11a)。如图 11a所示,在A=5.5后,PS评分不再随A的增大而增加,说明A=5.5是适用于该模型的相对合理的方差订正系数。订正后(A=5.5),该模型回报结果的平均PS评分技巧从71.00分进一步提升至82.10分,其变化范围为70分~90分,展现出稳定的高回报技巧(图 11b)。

      图  11  回报降水距平百分率的PS评分对方差订正系数的响应函数(a)和方差订正前后的PS评分(b)

      Figure 11.  Response PS score curve of the reforecast percentage of rainfall anomaly to the variance correction parameter(a) and PS score before and after variance corrected(b)

    • 由于传统的相关分析和主观挑选预测因子进行季节预测存在较大不确定性,本文提出PSM和预测因子自动选择器,并基于二者建立我国夏季汛期降水异常的季节预测模型。该模型对我国夏季汛期降水异常具有稳定的较好历史回报技巧。具体结论如下:

      1) PSM采用交叉检验思路,能够反映预测因子对预测对象的潜在预测技巧,且不受极端样本影响。PSM完全基于实际的统计预测建模流程,在逻辑上体现寻找预测对象的充分条件的特点,是对传统相关系数分布图(寻找必要条件)的重要补充。结合两者结果,可揭示符合预测对象充要条件的预测因子。

      2) 预测因子自动选择器借鉴集合预报思想,从PSM中挑选出具有最显著潜在预测技巧的预测因子,再通过大量预测结果的集合产生最终预测产品,实现逐年自动挑选预测因子的功能。该方案不仅克服了依赖于预报员主观经验选择预测因子的不足,也为进一步深入东亚夏季风年际变率可预报性研究提供新切入点。

      3) 基于PSM和预测因子自动选择器的统计预测模型对我国夏季汛期降水异常的回报技巧较高。在基于前期秋季、冬季预测因子的21年回报试验中,预测结果和观测的平均符号一致率为60%,平均空间相关系数为0.436,平均PS评分为71.00分,经方差订正后,平均PS评分可达82.10分,远高于现有动力模式的预测技巧。

      目前,动力模式直接输出的降水预测技巧偏低,但其对东亚夏季风主要环流系统(西北太平洋副热带高压、东亚大槽、高空副热带西风急流等)的预测技巧正在逐步提高。因此,通过统计方法利用动力模式输出的环流信息提高我国夏季汛期降水异常预测能力成为气候预测研究的新热点,且已有一些动力-统计相结合预测汛期降水异常的成功案例[59-61]。本文采用的预测因子既包含低纬度降水,也包含中高纬度高空环流,动力模式对二者具有一定的预测能力,如何将本文提出的统计预测模型和动力预测结果相结合,进一步发展全新的动力-统计相结合预测模型,将成为未来的工作重点。

参考文献 (61)

目录

    /

    返回文章
    返回