APP下载

中国极端天气气候研究

2024-04-25陈海山张耀存张文君尹志聪华文剑况雪源陈国森马红云韩婷婷

大气科学学报 2024年1期
关键词:归因机理预测

陈海山 张耀存 张文君 尹志聪 华文剑 况雪源 陈国森 马红云 韩婷婷

摘要 全球变暖背景下,极端天气气候事件频发,并表现出群发性、持续性、复合性等特点,不可预测性增加;持续性强降水、极端低温、复合型极端高温干旱、群发性热浪和台风等极端天气气候事件对我国经济社会和可持续发展影响巨大。然而,上述极端天气气候事件的新特征、关键过程和机理尚不完全清楚,重大极端事件的预报预测水平亟待提升。文章首先简要介绍“地球系统与全球变化”重点专项项目“中国极端天气气候事件的形成机理及其预测和归因”的基本情况。项目拟在分析全球变化背景下对我国造成重大影响的极端天气气候事件新特征的基础上,深入研究多尺度海-陆-气耦合过程影响极端天气气候事件的机理,挖掘极端天气气候事件次季节-季节预测的前兆信号;发展动力与物理统计相结合的极端事件预测新方法,研制针对中国极端事件的新一代高分辨率数值预报与检测归因系统。文章重点总结了自2022年12月项目立项至今取得的最新研究成果和进展。

关键词极端天气气候事件;海气相互作用;陆面过程;海陆气耦合;机理;预测;归因

1 研究背景

全球变化背景下,极端天气气候事件频发,并表现出群发性、持续性、复合性等特点,不可预测性增加。持续性强降水、极端低温、复合型极端高温干旱、群发性热浪和台风等极端天气气候事件对我国经济社会和可持续发展影响巨大。然而,上述极端天气气候事件的新特征、关键过程和机理尚不完全清楚,重大极端事件的预报预测水平亟待提升。

首先,辨识全球变化背景下极端天气气候事件的新特征是备受关注的热点问题(Coumou and Rahmstorf,2012;Zscheischler et al.,2020)。IPCC第6次评估报告(AR6)指出:20世纪50年代以来,气候变暖导致全球水循环加剧,极端天气气候事件的变率增大。极端降水事件发生频次增加,强度增强,极端性愈发明显;极端高温事件增多,强度和持续性加强;冷事件减少,但极寒事件也时有发生(Seneviratne et al.,2021)。此外,随着全球变暖加剧,复合型极端高温干旱、群发性热浪变得更为频繁(Rahmstorf and Coumou,2011;Zscheischler et al.,2018)。近20年来,海洋热浪发生的频率更高、持续时间更长、范围更广、强度更大,强台风比例增加、登陆中国的台风平均强度增强(中国气象局气候变化中心,2022)。近年来,我国极端强降水和大范围极端高温事件显著增多,持续时间增长,破纪录的极端事件频繁发生,如2020年长江流域的超强梅雨、2021年7月河南郑州的极端强降水、2022年夏季持续的高温热浪(Yin et al.,2023)。在全球变暖背景下,全球不同区域和我国发生的极端天气气候事件表现出与以往显著不同的特征(AghaKouchak et al.,2020;Fischer et al.,2021),亟须加强极端天气气候事件新的变化特征、规律及其成因的研究。

其次,理解极端天气气候事件的形成机理是全球变化领域的研究难点。国内外学者在极端天气气候事件的形成机理及其归因方面取得了一些重要进展(Nie et al.,2020;Neelin et al.,2022;Barriopedro et al.,2023),强调了人类活动对极端天气气候事件变化的作用。值得注意的是,海-陆-气相互作用作为全球气候变化的重要驱动力,在极端天气气候事件的形成中也起到至关重要的作用(Buzan and Huber,2020;Thiery et al.,2020;Wouters et al.,2022)。一方面,极端天气气候事件的发生发展受到次季节到年际等多尺度海气耦合主要模态(如ENSO、IOD、MJO等)的调控,这些典型海气模态同时又受到全球增暖的显著影响(Jin et al.,2014;Cai et al.,2019;Hsu et al.,2021);另一方面,陆面过程与陆气相互作用通过改变地表能量和水分循环影响区域极端天气气候事件,且上述过程受到不同海洋热力状态的明显调制(Seneviratne et al.,2010;Schumacher et al.,2022)。然而,海气和陆气相互作用过程在全球气候变化的背景下呈现出多样性和复杂性,对其物理事实和动力机制的认识还有待提高。此外,陆面过程对极端天气气候的影响尚未受到足够的重视,陆面反馈在海洋影响极端天气气候事件过程中的作用的相关研究还十分薄弱。因此,深入理解全球变暖背景下海气相互作用和陆面过程变化的机理,探究多尺度海-陆-气耦合过程与我国极端天气气候事件的联系,明确其影响极端事件的关键物理过程,有望开辟我国极端天气气候事件变化研究的新思路,为我国极端天气气候事件的预测提供理论基础。

再者,提高极端事件的预测水平是当前面临的重大科学挑战。相对于平均气候,极端事件与气候灾害的联系更紧密,预测难度更大。目前,国内发展了一系列气候预测理论和方法,但针对极端天气气候事件尚缺乏成熟的预测理论和方法。在物理统计预测方面,年际增量、深度学习和最优子集回归等方法可以在一定程度上提取前期预测信号,进而预测暴雨、暴雪、干旱和高温等极端事件的发生频次等(张德宽等,2003;Fan and Tian,2013;Han et al.,2020;Tian and Fan,2022;Xie et al.,2023);但已有的物理统计模型多从海温、海冰、土壤、积雪等要素中直接提取预测信号,未能充分考虑气候系统多圈层相互作用过程,尤其是多尺度海-陸-气耦合过程中的有效信息。在动力预测方面,模式物理过程的完善、分辨率的提升、模式初值的改善和集合成员的增加能够提升动力系统对于极端事件的预测能力,且概率预测的结果优于确定性预测(Becker et al.,2013;Xie and Zhang,2017)。目前,欧美和我国都在发展更高分辨率的季节预测系统,或利用动力降尺度方法进一步提升对极端事件的分辨能力。此外,如何将物理统计模型和动力预测的结果有效集成?目前主要依靠主观经验,尚缺少客观化的集成方案,在一定程度上制约了我国极端事件的预测水平。因此,如何充分挖掘多尺度海-陆-气耦合过程中的预测信息,提升动力预测系统对极端事件的预测性能,并有效地将动力预测和物理统计模型的结果进行集成,是目前制约极端事件预测性能提升的重要问题,相关研究突破将为我国防灾减灾提供更直接、更有效的科学支撑。

综上所述,迫切需要深入分析全球变化背景下极端天气气候事件变化的新特征,从海-陆-气耦合的角度来探讨海洋和陆面影响极端事件的多尺度关键过程和机理,发展极端事件次季节-季节预测的新方法,构建中国新一代高分辨率的极端事件数值预报与检测归因系统,提升我国极端天气气候事件的预测水平,为防灾减灾和经济社会可持续发展提供有力支撑。

2 项目的关键科学问题和预期目标

科技部于2022年启动实施了“地球系统与全球变化”重点专项,2023年南京信息工程大学牵头,联合南京大学、中国科学院大气物理研究所、中国气象科学研究院、国家气候中心等高校、科研院所和业务部门申报的“中国极端天气气候事件的形成机理及其预测和归因”项目获科技部立项资助。

项目拟围绕“中国极端天气气候事件的形成机理及其预测和归因”这一核心主题开展多学科交叉集成研究,回答以下关键科学问题:1)在全球变化背景下,极端天气气候事件呈现出什么新的特点和演变规律?人类活动对这些极端天气气候事件有多大影响?2)海气相互作用、陆面过程、海-陆-气耦合影响重大极端天气气候事件的关键过程和机理是什么?3)如何获取极端天气气候事件预测的多尺度前兆信号,发展极端天气气候事件预测的新方法,建立基于动力降尺度与物理统计相结合的新一代高分辨率极端事件数值预报与检测归因系统?

围绕以上关键科学问题,项目拟在分析全球变化背景下对我国造成重大影响的极端天气气候事件新特征的基础上,深入研究多尺度海-陆-气耦合过程影响极端天气气候事件的机理,挖掘极端天气气候事件次季节-季节预测的前兆信号;发展动力与物理统计相结合的极端事件预测新方法,研制针对中国极端事件的新一代高分辨率数值预报与检测归因系统。

项目拟以变化事实、物理机制、预测方法和预测系统为主线,设立4个课题:第1课题,极端天气气候事件的变化事实、演变特征与归因;第2课题,多尺度海气相互作用影响极端天气气候事件的关键过程;第3课题,陆面过程与海-陆-气耦合对极端天气气候事件的影响及其机理;第4课题,极端事件的预测方法和高分辨率数值预报与检测归因系统(图1)。

项目研究旨在:1)揭示全球变化背景下对我国造成重大影响的极端天气气候事件的新特征;2)厘清多尺度海-陆-气耦合过程影响极端天气气候事件的机理,明确极端天气气候事件预测的前兆信号;3)发展基于动力与物理统计相结合的极端事件预测的新方法,研制针对中国极端事件的新一代高分辨率数值预报与检测归因系统,提高极端事件检测归因与预报预测业务水平,为极端天气气候预测业务提供科学支撑。同时,培养一批活跃在大气科学前沿,尤其是气候变化与极端天气气候研究领域的青年学术带头人。

3 项目最新研究进展

3.1 极端天气气候事件的检测识别

极端天气气候事件的检测识别是开展极端事件监测和研究的基础性工作,项目在群发性剧烈降温事件检测、基于逐日土壤湿度的干旱过程识别及中国群发性区域极端降水检测识别方面开展了研究。

1)中国近60 a群发性剧烈降温事件的检测识别:利用中国1961—2020年724站逐日平均气温资料,依据24 h(或48 h)降温超过10 ℃的标准,对中国区域冬季群发性剧烈降温天气事件进行检测,识别出388次群发性剧烈降温事件。表1为根据综合指标给出的1961—2020年冬季群发性剧烈降温前10名事件。不难发现,其中近20 a就出现了3次(2008年12月3—6日;2010年1月20—23日;2016年2月13—15日),这意味着在全球变暖背景下,群发性剧烈降温事件发生频繁且严重。从各年代不同等级群发性剧烈降温事件发生频数来看,21世纪以来群发性剧烈降温事件在持续天数及影响站数等指标上都高于上个世纪末,这种变暖背景下剧烈降温事件给社会发展造成了更大的危害,值得更多关注。

2)基于逐日土壤湿度的干旱起止时间的检测识别:在全球变暖背景下,干旱事件频发,如何在日尺度上识别干旱的持续时间是极端干旱研究的基础性工作。项目利用ERA5提供的逐日土壤湿度数据,计算每天的土壤湿度相对变率Di=SMi-SMSM,i=1,2,…,n,其中SMi为某天的土壤湿度,SM为同一天土壤湿度的多年平均值。之后利用日土壤湿度变率与每年土壤湿度的-1.5倍标准差进行比较,其中土壤湿度变率在1个月(30 d)内有至少20 d小于等于-1.5倍标准差的判定为1次干旱事件,第一次满足该条件的识别为干旱起始时间,最后一个不满足该条件的识别为干旱结束时间。图2为华中地区2017年3—10月逐日土壤湿度相对变率。由图2可较好识别出2017年6月27日—8月6日华中地区的持续干旱。

3)中国群发性区域极端降水的检测:在全球变暖背景下,极端降水事件频率和强度显著增强,区域性群发极端降水事件频繁发生。项目研究基于再分析格点资料并利用高相关范围的面积来表征极端降水过程的区域性强度。某格点极端降水的区域相关性是指该格点与其周围格点极端降水频数的相关情况,表示该格点发生极端降水时,其他格点同期也发生极端降水的可能性大小。高相关区面積越大,区域相关性相对越好。结果表明:中国夏季和秋季最易发生区域性极端降水,而春季区域性极端降水的强度自21世纪后有所减弱,冬季区域性极端降水主要发生在东南部地区。具体而言,东北极端降水区域性在春季和冬季不显著,而在夏季和秋季均显著且增强,即21世纪后,东北地区夏季和秋季更易发生大范围的极端降水;华北极端降水区域性春季减弱,夏季和秋季增强;江南一带全年极端降水均具有明显区域性特征,且冷季强于暖季。

3.2 极端天气气候事件的时空变化、环流型及动力学机理与归因研究

项目最近的工作围绕暖干复合事件、极端寒潮、区域性极端降水的时空变化、环流型及动力学机理开展了研究。此外,针对典型极端天气气候事件(2022年夏季长江中游极端热浪)开展了归因研究。主要进展如下:

1)揭示了我国东部地区暖干复合事件频数变化的主模态及其与异常环流型的可能联系(Huang D Q et al.,2023)。基于我国台站观测和再分析资料研究表明,我国东部地区的暖干复合事件频数变化的主模态主要为一致变化型和南北偶极子型;与之伴随的强度变化表现为北方强烈的偏暖和偏干,而南方地区则以偏暖为主。两个主模态分别与减弱的陆地上东亚副热带急流和南移的海洋上的东亚副热带急流相关(图3)。进一步分析发现,太平洋年代际振荡的负位相、大西洋高纬度地区的变暖和北极振荡的正位相共同维持急流的减弱和暖干复合事件频数的一致变化型;赤道太平洋和印度洋的变暖、北大西洋振荡的负位相和日本-太平洋遥相关型共同维持急流的南移和暖干复合事件频数的偶极型。

2)利用新的动力学变量-局地有限振幅波活动(Local finite-amplitude Wave Activity,LWA),探讨了中国东部极端寒潮的前兆信号、爆发条件和关键动力学机理(Shi et al.,2023)。结果表明,约40%的极端寒潮具有以下重要环流特征:在极端寒潮爆发的前4 d左右,巴尔喀什-贝加尔湖地区会出现显著的LWA正异常,这是一个明显的前兆信号(图4a),对应经向增强的环流特征,并伴随着冷空气在欧亚大陆中高纬地区堆积。上述LWA正异常逐渐增强并东移至贝加尔湖以东(图4b),形成局地极端LWA异常(超过95%百分位)。通过对比不同类型的寒潮发现,贝加尔湖以东出现极端LWA异常是爆发极端寒潮的必要条件。LWA可作为我国极端寒潮事件监测和预测的重要因子。

3)揭示了地形复杂地区(青藏高原中东部与四川盆地)夏季区域性极端降水主控环流天气型(Xu X K et al.,2023)。利用GPM卫星降水资料、ERA5大气再分析资料以及谱聚类方法研究发现,青藏高原中东部地区夏季西太平洋副热带高压的南退配合高层南亚高压的东进控制着高原夏季区域性极端降水主要在区域内分布偏南(Ⅰ型占比为55%),而西太平洋副热带高压的北抬与高层南亚高压的东进北抬则控制了高原夏季区域性极端降水在区域内分布偏北(Ⅱ型占比为45%)(图5)。对四川盆地而言,西太平洋副热带高压与南亚高压的相向北抬/南退配合低层西南涡的位置略微偏移控制了夏季区域极端降水分布偏西北(Ⅰ型占比为52.6%)/偏东南(Ⅱ型占比为47.4%)。

4)分析了江淮流域区域性极端降水事件的变化特征并揭示了其对应的主控环流(Zeng et al.,2023)。利用逐小时降水台站资料和ERA5再分析资料,研究发现江淮流域区域性极端降水事件对总的夏季极端降水量和极端降水频率的贡献都较大,大值集中在靠近长江的区域,且区域性事件的贡献随着气候变化呈现出了一定的增加趋势。进一步采用谱聚类方法,研究发现长江中下游区域性极端降水事件的主控环流为受东亚季风影响的梅雨锋控制的梅雨锋型(占比为92%)和登陆台风影响的台风型(占比为8%),其中梅雨锋型可进一步分为偏南梅雨锋型和偏北梅雨锋型。偏南梅雨锋型频发于6月,降水落区偏南,峰值多在早上8时,主要受夜间惯性振荡产生的低空急流的动力影响。偏北梅雨锋型频发于6—7月,降水落区偏北,峰值多出现在下午17時,主要受午后热力影响。台风型主要发生在8月,降水落区位于登陆台风东北侧(图6)。

5)开展了典型极端天气气候事件(2022年夏季长江中游极端热浪)的归因研究(Hua et al.,2023a):2022年夏季,北半球许多地区经历了极端高温天气。例如,欧洲遭遇了前所未有的热浪,英国气温首次超过40 ℃。在亚洲,中国经历了有记录以来范围最广、持续时间最长的热浪,其中长江中下游地区夏季热浪和干旱尤为严重。利用多源地面站点和格点观测数据、再分析资料(ERA5和JRA-55),借助先进地球系统模式(CESM2)的超级集合试验、偏差校正方法等,分析了2022年夏季长江中游极端热浪的特征,量化了气候系统内部变率和人为强迫对2022年夏季极端热浪的贡献,最后根据未来气候变化预估了此类极端事件发生的概率。结果表明,2022年长江中游地区有两个极热时期(集中在7月初和8月中旬),导致夏季极端温度超过历史最大值。在7 d到40 d的时间尺度内,破纪录值的区域约占研究区域的83%。长江中游地区最高气温比气候平均值(1981—2010年)高13.1 ℃,重现期长达数万年。针对2022年长江中游夏季高温热浪事件分析发现,异常的大气环流和持续的高压系统促成了2022年长江中游的高温热浪,使得此类事件发生的可能性增加了6~7倍(图7)。CESM2和其他CMIP6耦合模式的多模式集合模拟表明,这种异常环流主要是由大气和下垫面条件的内部变率造成的,而人为气候变化则增强了此类事件发生的可能性。随着气候变暖,在相对较高的SSP3-7.0排放情景下,预计到2070年至2100年,长江中游地区罕见的热浪将变得更加普遍,并且可能每8.5 a发生一次。发生此类事件的可能性大幅增加,主要是由于长江中游地区日最高温度的概率密度函数(PDF)平均态偏移(正偏移),其次是变率增加(PDF平坦化)。

3.3 热带太平洋海温变异机理、未来变化及其影响研究

海气相互作用与极端事件具有密切的联系,多尺度海气相互作用及其相关机理研究是从海气相互作用的角度深入理解极端事件形成机理的前提。项目在热带太平洋海温变异机理、未来变化及其影响方面取得了阶段性进展。具体如下:

1)揭示了连续性拉尼娜事件的独特海洋过程及前兆信号(Liu F Y et al.,2023)。ENSO具有多样性和复杂性,厄尔尼诺通常在成熟后迅速衰退,而拉尼娜往往在第二年持续,甚至在随后冬季再次加强,造成连续性拉尼娜(MLN)事件,其物理过程有待阐明。MLN相关的异常纬向风应力与海表温度呈现出一致的持续性特征,但是纬向流却显示出规律振荡,在发展年8月至第2年8月期间呈现出近一年的周期特征(图8)。基于纬向流的这种独特现象,研究定义了Uo-MLN指数(100°~170°W、5°S~5°N区域平均的正位相与负位相纬向流平均值的差值)。该指数刻画了两次纬向流异常反转的海洋过程,区分了多年拉尼娜和单年拉尼娜事件。MLN事件发展年4—10月太平洋海洋热力状况与第二个冬季Nio3.4指数呈现出高负相关关系(图9),表明赤道中东太平洋的前期热含量状态为MLN提供了一个潜在的前兆信号。

2)探讨了未来全球变暖背景下强El Nio的可能变化及其对东亚冬季降水的影响(Huang Y et al.,2023)。强El Nio事件会引起全球显著的天气气候异常,造成区域极端性干旱或者洪涝,尤其会引起中国冬季南方出现持续性降水异常。然而,在未来全球变暖背景下,强El Nio事件对中国冬季降水影响是否发生变化值得关注。以往研究表明,强El Nio事件通过罗斯贝波响应激发孟加拉湾至南海区域出现低空异常的反气旋,进而使得中国华南降水增多。项目基于CMIP模式发现,在未来高排放情景下强El Nio事件可造成孟加拉湾至南海区域反气旋出现北移现象,这一北移有利于东亚北部地区降水增多、南部地区降水减少(图10)。进一步研究发现,华南地区降水负异常随着排放增多出现减少的趋势,而冬季江淮流域以北地区降水增多的特征随着排放增多更加显著。

3)提出了热带海温变率与陆面过程的相互作用影响台风活动的可能機制(Cai et al.,2023)。研究表明,20世纪末以来,青藏高原冬春积雪异常引起的热力作用改变了局地位势高度异常,导致热带印度洋海盆海温冷异常。随后的夏秋季节,高原积雪导致的海陆热力差异较小,减弱了孟加拉湾爆发的西南季风,加强了印度洋东西海温梯度,驱动了印度尼西亚西海岸的异常东风,进一步增大了印度洋海温异常的纬向差异。此海洋-大气相互反馈的过程促使印度洋在台风快速增强活跃期加速印度洋类偶极子模态的演变。研究结果强调了青藏高原热力作用在印度洋海气相互作用中的重要性,也揭示了印度洋海温对快速增强台风的影响机制,为理解热带气旋强度的变化以及热带气旋强度的季节预测提供了新思路(图11)。

3.4 关键陆面因子的异常变化与机理研究

项目围绕积雪、地表蒸散发、植被等关键陆面因子的异常变化及相关机理开展了大量的研究。此外,近期研究还评估了土地利用/土地覆盖变化的气候效应。相关工作为后续研究陆面过程和海陆气耦合对极端事件的影响奠定了基础。主要进展如下:

1)揭示了欧亚大陆中高纬秋季积雪/冬季降雪异常的主要模态及其关键环流因子(Zhang et al.,2023a,b)。欧亚中高纬地区深秋(11月)积雪异常的东-西偶极子型模态能够通过调控对流层-平流层相互作用过程,对后期冬季北极涛动和北大西洋涛动等大范围大气环流系统以及区域性极端气候产生不可忽视的影响。研究表明,东-西偶极子型格局是欧亚深秋积雪年际变化的主导模态,其表现为欧亚中纬度东、西部积雪反向变化的特征(图12a、b)。研究发现,欧亚中高纬冬季降雪年际变化的第一、二模态分别表现为欧亚高纬地区、西西伯利亚地区降雪量的一致性变化特征(图12c—f)。秋季东-西偶极子型积雪异常格局受到同时期东大西洋型、东大西洋-西俄罗斯型和斯堪的纳维亚型三种大气遥相关型以及北大西洋经向偶极型异常海温和巴伦支海域海冰变化的共同作用,它们通过影响欧亚大陆上空异常大气环流及其导致的气温和降雪量的变化控制积雪异常模态的形成。诊断分析和数值试验表明,异常大气环流及其引起的大气动力条件和水汽输送状况的变化是导致降雪异常的主要原因,与冬季降雪第一模态相关的异常环流可以归因于极地-欧亚型和斯堪的纳维亚型遥相关的影响,二者解释了约2/3的变化;第二模态除了受到欧亚型遥相关和北大西洋涛动等大气内部变率的影响外,还与巴芬湾的海冰异常有关。

2)探讨了欧亚大陆中高纬度春季融雪的基本特征及其驱动机制(Sun and Chen,2023)。欧亚大陆中高纬度春季融雪能显著影响局地水文过程和后续夏季东北亚气候变异。研究发现,欧亚大陆中高纬度春季融雪主要表现为两种类型:第一类异常融雪型,前冬降雪异常积累伴随春季异常偏暖加剧融雪;第二类异常融雪型,虽然前冬积雪较少,但是春季强烈增暖更利于融雪。具体而言,第一类融雪事件中,前冬关键区受异常低压槽控制,利于水汽输送和局地辐合,伴随低层上升运动,促进前冬积雪,为融雪提供有利的初始条件。4月关键区受高压脊控制,异常暖平流和下沉运动导致地表升温,融雪增强。第二类融雪事件中,冬季环流异常不明显,但是4月受到更强的异常高压控制,导致了更剧烈的融雪过程。局地地表能量平衡的分析结果表明,大气向下的感热加热是影响第一类融雪事件的主要因素,短波辐射起次要作用,而第二类融雪事件中短波辐射的作用更强,感热加热次之(图13)。

3)基于Shuttleworth-Wallace模型构建了1982—2015全球5 km逐月潜在蒸散发数据集(Sun et al.,2023a)。作为蒸散发(evapotranspiration,ET)的理论上限,潜在蒸散发(potential evapotranspiration,PET)一直以来都是各个学科重点关注的一个关键变量。然而,由于非真实的下垫面设置和不合理的模型参数化,目前PET数据集及主要PET模型依然存在诸多不确定性。项目研究充分考虑植被参数(如植被气孔导度、CO2效应对植被气孔导度的影响)的种间差异,并应用经过严格筛选的全球96个通量塔观测的PET(即不受水分胁迫的ET)数据,对Shuttleworth-Wallace(SW)双源模型进行了率定和参数化。基于1982—2015年4套气象数据(MSWX-Past、CRU TS4.06、ERA-5和MERRA-2)、重制的全球植被冠层高度、土地利用/覆盖数据集(GLASS-GLC)、叶面积指数数据(GLASS AVHRR LAI)及格点化全球CO2浓度数据,生产了全球(180°~90°W~0°~90°E~180°,60°S~85°N)逐月5 km的PET、潜在蒸腾(potential transpiration,PT)和潜在土壤蒸发(PE)产品,并最终提供了PET、PT和PE的集合产品(https://doi.org/10.11888/Terre.tpdc.300193)。

4)揭示了青藏高原主要生长季植被的年际变化特征及其对印度夏季风活动的响应机制(Mao et al.,2023)。研究发现,印度夏季风不仅通过调节高原降水来影响高原植被,而且还能通过引起高原热力因子的变化来影响高原植被生长。在强印度夏季风的情况下,其相应的对流层上层和中下层环流结构有助于将更多的水汽输送到高原,并增强高原上空的对流。增强的对流和降水能够促进高原植被的生长,表现出一致的生长变化。与此同时,高原植被与热源之间存在相互作用,一方面植被会影响高原地表热力情况,而变化的热源又会促进高原植被生长。并且,多元回归分析表明,印度夏季风及其引起的局地气候因子的变化可占高原主要生长季植被年际变化的52%以上。研究结果有助于加深对高原植被与区域及全球气候变化联系的理解、提高未来对高原植被年际生长变化的预测。

5)评估了土地利用/土地覆盖变化的气候效应(Liu S Y et al.,2023)。人类活动导致的土地利用变化极大地改变了陆地表面(如毁林等),对极端气候也有重要影响。然而,毁林变化的局地和非局地影响仍存在很大的不确定性。项目通过模拟研究,首先强调了森林变化非局地气候效应的重要性,如果不考虑非局部效应,那么云量变化与森林砍伐之间的关系将被误解(Hua et al.,2023b);其次揭示了毁林的生物物理影响,量化了不同生物物理因素的相对贡献以及强调了不同模式对温度模拟的差异可能主要受非局地大气反馈的调制。

3.5 陆面过程与陆气耦合过程对极端事件的影响

项目在陆面过程与陆气耦合过程对极端事件的影响方面也取得了阶段性成果,具体包括:

1)阐明了地球三极热力异常强迫与2022年北半球群发性極端事件的关联(Zhang et al.,2024)。在全球变暖背景下,地球三极(北极、南极、青藏高原)的海冰、陆面冰川和积雪融化速率加快,水分相态的变化通过反照率效应和水文效应等影响海表和地表净辐射、地表能量分配、水循环等,进而通过下垫面反馈过程影响局地和全球天气气候。项目针对2022年夏季极端事件及其与三极热力异常的关联性开展了研究。图14给出了2022年7—8月北半球气温异常和极端降水的空间分布,可见北美西部、西欧-东欧、东西伯利亚和中国长江流域等区域出现了3 ℃以上的温度正异常,反映了这些区域的夏季极端高温事件。此外,在温度负异常的几个区域则出现了3次破纪录的极端强降水事件,包括巴基斯坦夏季洪涝、5月21日—6月21日华南龙舟水、8月7—11日韩国洪涝和中国东北夏季极端降水,上述极端事件具有明显群发性特征。进一步分析发现,上述群发性极端事件与北极、南极和青藏高原的水热异常之间存在密切关联;而且三极水热异常对北半球极端事件的群发性有重要影响。北极暖异常减小北极与中低纬度的温度梯度并引起NAO模态变化,进而造成西风急流、波列及副热带高压向北伸展,同时造成波动振幅加大,有助于阻塞活动以及极端温度和降水发生;另外,春季北极暖异常也影响春季陆面雪盖和土壤湿度异常,陆面反馈为夏季极端温度的发生提供了非绝热加热。南极暖异常减小南极与中低纬度的温度梯度,进而造成波列及副热带高压向南极伸展;作为亚洲季风的成员之一,偏南的马斯克林高压将拉大季风环流的空间范围,进而减弱南北半球季风环流系统之间的气压梯度,从而减弱亚洲季风。此外,南极暖异常也会使得整个季风环流系统偏南,有助于水汽在巴基斯坦及华南辐合,加强暴雨发生;南极增暖还通过影响垂直洋流变化影响AMOC,进而对北半球天气气候产生影响。偏高的青藏高原春季土壤湿度减少了高原热力作用,一方面促进了北非和西太平洋副热带高压的对接,另一方面使得亚洲季风减弱。此外,三者均可激发异常涡旋和波列直接影响极端天气气候(图15)。

2)揭示了中南半岛前期土壤湿度异常对2020年超强梅雨的影响及相关机制(Dong et al.,2023)。2020年夏季江淮流域发生的超强梅雨造成了巨大的经济损失和人员伤亡,最新的研究从大气环流异常、海温异常等角度探究了此次超强梅雨发生的可能原因,但陆面过程的影响并未受到关注。项目组基于观测分析和数值模拟探讨了陆面过程对此次超强梅雨的可能贡献,结果表明中南半岛前期的土壤湿度异常对此次超强梅雨有重要的影响。研究发现,5月中南半岛土壤湿度异常与6—7月长江中下游地区降水量显著相关(图16),中南半岛地区前期干的土壤有利于地表温度升高,土壤的“记忆性”使得这种异常可以持续到梅雨期。增暖的地表导致感热通量增加从而加热低层大气,有利于位势高度抬升,使中南半岛到西北太平洋上空出现异常反气旋,西太平洋副热带高压西伸加强,长江中下游地区上空的西南风和垂直运动增强,并导致水汽输送和水汽辐合加强,最终造成梅雨期降水显著增加。基于WRF模式的数值试验再现了观测的环流和降水异常(图17),结果表明中南半岛5月土壤湿度异常使6—7月长江中下游地区降水量增加10.6%。中南半岛土壤湿度减少有利于地表温度升高,有更多的能量以感热形式加热大气。尽管加热局限于中南半岛区域,但被中南半岛地表加热的暖空气在背景气流西南风的作用下使下游大气升温、位势高度增加,同时暖平流使南北温度梯度增加,有利于持续性锋面降水的发生。

3)发现春季中南半岛生物质燃烧可通过陆气反馈增加长江流域夏季降水异常(Wei et al.,2023)。中南半岛是全球生物质燃烧(BB)最严重的地区之一,每年春季排放大量的BB气溶胶。但上述区域BB气溶胶如何影响区域气候和陆气耦合尚不清楚。通过模式试验和观测分析发现,BB气溶胶使中南半岛大气更稳定,从而减少降水和土壤含水量(图18)。由于土壤湿度的记忆性及副高的响应和反馈,较干的土壤从春季持续到夏季,会进一步使副高西伸,并最终导致更多水汽向长江流域输送,从而增加长江流域降水(图19)。此外,BB气溶胶还会使中南半岛土壤湿度异常对长江流域降水的影响几乎加倍,相关结果证实BB气溶胶对区域陆气相互作用及气候具有十分重要的影响。

4)揭示了CO2浓度升高的植被生理效应对干旱影响的区域性差异及相关机制(Sun et al.,2023b)。基于两套模拟试验(即有无CO2效应),深入分析了CO2效应对中国1971—2010年PET(potential evapotranspiration)、SPEI(standardized precipitation-evapotranspiration index)及干旱特征的影响。结果发现,CO2效应可以导致PET下降、SPEI增大以及干旱强度和持续时间减小;但是这种影响存在明显的区域性差异,总体呈现为北方的影响最大(除中国最西部的PET)。通过公式推导分析发现,这种空间性差异不仅与CO2浓度存在一定的函数关系,而且与气候背景条件也存在一定的函数关系,如:就PET而言,为气候态的风速(U)和饱和水汽压斜率(s),而就SPEI而言,则为气候态的U和气候态的降水(P)。进一步的归因分析显示,气候态U的空间差异决定了CO2引起的PET变化的空间差异,而气候态U和P的空间差异决定了CO2引起的SPEI变化的空间差异。

3.6 极端事件预测相关研究

围绕既定目标,项目在极端事件的预测方面开展了初步研究,代表性成果包括:

1)中国东部极端温度预测研究(Pan et al.,2024)。夏季,长江中下游(YRB)极端高温天气频发,然而动力模式(例如美国环境预报中心第二代气候预报系统CFSv2等)在预测YRB夏季极端高温日数方面存在局限,无法很好地捕捉其年际变化、趋势和极值。因此,项目组基于年际增量方法,使用4个预测因子,即4—5月欧洲西北部融雪、3月中西伯利亚高原积雪深度、CFSv2预测的夏季海温和青藏高原上空200 hPa位势高度,建立了长江中下游极端高温日数(HD35)动力-统计相结合的季节预测模型(HMYRB)。HMYRB在1983—2015年的交叉检验结果中表现出色,预测值与观测值之间的相关系数为0.58,同号率(PSS)为76%(图20)。在独立试报中(2016—2022年),HMYRB仍保持高PSS技巧(为86%)。

2)极地-中高纬度信号“暖北极-冷欧亚”(WACE)模态和中国东部冬季气温的预测研究(Xu T B et al.,2023)。极地-中高纬度信号是中国东部冬季气候可预测性的重要来源,其中WACE模态与中国冬季气候关系密切。因此,准确预测WACE对改善中国冬季气候预测具有重要意义。然而,CFSv2对WACE的预测能力非常有限。

为实现对“暖北极-冷欧亚”模态指数(GradTAE)的准确预测,课题组基于年际增量方法,选择了前期观测的9月巴伦支-喀拉海海冰、10月北大西洋海温、9月北美土壤湿度和CFSv2预测的同期巴芬湾-拉布拉多海区域海冰作为预测因子,建立了动力-统计混合预测模型。该模型对GradTAE及其年际增量(GradTAE-DY)的年际变化、極值和3个不同时间段的趋势都有良好的预测效果(图21)。相比CFSv2,其预测技巧有了很大提升。在2016—2020年的独立预测中,预测模型保持了较高的准确性。同时对2021/2022年及2022/2023年冬季WACE的预测进行个例检验,其预测效果也保持稳定。CFSv2对中国东部气温的预测技巧较低,观测值与预测值之间的相关系数仅为0.02。基于建立的动力-统计模型对WACE预测,并将预测结果加入CFSv2,对中国东部气温的预测结果进行订正,这样可有效提高CFSv2对中国东部地表气温的预测效果(图22)。

4 结语

自2022年12月项目实施以来,项目紧密围绕既定研究内容,在极端事件的检测识别、极端天气气候事件基本特征与环流型及相关动力学机制、多尺度海气相互作用、陆面过程及其对极端事件的影响与极端事件的预测方面取得了一些进展,为项目的顺利实施奠定了坚实的基础。

然而,关于极端天气气候事件的研究是当前大气科学领域的重大挑战,如何聚焦中国极端天气事件变化的新特征,更多关注多尺度海-陆-气耦合过程影响极端天气气候事件的机理,挖掘极端天气气候事件次季节-季节预测的有效前兆信号,发展极端事件预测新方法,提升极端天气气候事件的预测水平,还任重道远。期待通过项目的研究,能在极端天气气候事件机理的科学认知上有新的突破,为极端天气气候事件预测水平提升、防灾减灾和经济社会可持续发展提供有力支撑。

参考文献(References)

AghaKouchak A,Chiang F,Huning L S,et al.,2020.Climate extremes and compound hazards in a warming world[J].Annu Rev Earth Planet Sci,48:519-548.doi:10.1146/annurev-earth-071719-055228.

Barriopedro D,García-Herrera R,Ordóez C,et al.,2023.Heat waves:physical understanding and scientific challenges[J].Rev Geophys,61(2):e2022RG000780.doi:10.1029/2022rg000780.

Becker E J,van den Dool H,Pea M,2013.Short-term climate extremes:prediction skill and predictability[J].J Climate,26(2):512-531.doi:10.1175/jcli-d-12-00177.1.

Buzan J R,Huber M,2020.Moist heat stress on a hotter earth[J].Annu Rev Earth Planet Sci,48:623-655.doi:10.1146/annurev-earth-053018-060100.

Cai W J,Wu L X,Lengaigne M,et al.,2019.Pantropical climate interactions[J].Science,363(6430):eaav4236.doi:10.1126/science.aav4236.

Cai Y H,Zhao H K,Klotzbach P J,et al.,2023.Can Tibetan Plateau snow depth influence the interannual association between tropical Indian Ocean sea surface temperatures and rapidly intensifying typhoons?[J].J Climate,36(11):3781-3800.doi:10.1175/jcli-d-22-0697.1.

Coumou D,Rahmstorf S,2012.A decade of weather extremes[J].Nature Climate Change,2:491-496.doi:10.1038/nclimate1452.

Dong Y S,Chen H S,Dong X,2023.Impact of antecedent soil moisture anomalies over the Indo-China Peninsula on the super Meiyu event in 2020[J].J Meteor Res,37(2):234-247.doi:10.1007/s13351-023-2144-4.

Fan K,Tian B Q,2013.Prediction of wintertime heavy snow activity in Northeast China[J].Chin Sci Bull,58(12):1420-1426.doi:10.1007/s11434-012-5502-7.

Fischer E M,Sippel S,Knutti R,2021.Increasing probability of record-shattering climate extremes[J].Nature Climate Change,11:689-695.doi:10.1038/s41558-021-01092-9.

Han T T,Li S F,Hao X,et al.,2020.A statistical prediction model for summer extreme precipitation days over the northern central China[J].Int J Climatol,40(9):4189-4202.doi:10.1002/joc.6451.

Hsu P C,Fu Z,Murakami H,et al.,2021.East Antarctic cooling induced by decadal changes in Madden-Julian oscillation during austral summer[J].Sci Adv,7(26):eabf9903.doi:10.1126/sciadv.abf9903.

Hua W J,Dai A G,Qin M H,et al.,2023a.How unexpected was the 2022 summertime heat extremes in the middle reaches of the Yangtze River?[J].Geophys Res Lett,50(16):e2023GL104269.doi:10.1029/2023gl104269.

Hua W J,Zhou L M,Dai A G,et al.,2023b.Important non-local effects of deforestation on cloud cover changes in CMIP6 models[J].Environ Res Lett,18(9):094047.doi:10.1088/1748-9326/acf232.

Huang D Q,Zheng Z L,Xiao X C,et al.,2023.The variation of compound warm dry extremes over eastern China and the associated East Asian subtropical jet in the boreal summer[J].J Geophys Res:Atmos,128(24):e2023jd039392.doi:10.1029/2023jd039392.

Huang Y,Ren H L,Kug J S,et al.,2023.Projected change of East-Asian winter precipitation related to strong El Nio under the future emission scenarios[J].Clim Change,176(7):1-21.doi:10.1007/s10584-023-03551-y.

Jin F F,Boucharel J,Lin I I,2014.Eastern Pacific tropical cyclones intensified by El Nio delivery of subsurface ocean heat[J].Nature,516(7529):82-85.doi:10.1038/nature13958.

Liu F Y,Zhang W J,Jin F,et al.,2023.New insights into multi-year La Nia dynamics from the perspective of a near-annual ocean process[J].J Climate,36(21):7435-7450.doi:10.1175/JCLI-D-22-0505.1.

Liu S Y,Hua W J,Zhou L M,et al.,2023.Local and non-local biophysical impacts of deforestation on global temperature during boreal summer:CMIP6-LUMIP multimodel analysis[J].J Geophys Res:Atmos,128(11):e2022jd038229.doi:10.1029/2022jd038229.

Mao X,Ren H L,Liu G,et al.,2023.Influence of the Indian summer monsoon on inter-annual variability of the Tibetan-plateau NDVI in its main growing season[J].Remote Sens,15(14):3612.doi:10.3390/rs15143612.

Neelin J D,Martinez-Villalobos C,Stechmann S N,et al.,2022.Precipitation extremes and water vapor[J].Curr Clim Change Rep,8(1):17-33.doi:10.1007/s40641-021-00177-z.

Nie J,Dai P X,Sobel A H,2020.Dry and moist dynamics shape regional patterns of extreme precipitation sensitivity[J].Proc Natl Acad Sci USA,117(16):8757-8763.doi:10.1073/pnas.1913584117.

Pan S F,Yin Z C,Duan M K,et al.,2024.Seasonal prediction of extreme high-temperature days over the Yangtze River Basin[J].Sci China:Earth Sci.doi:10.1007/s11430-023-1265-2.

Rahmstorf S,Coumou D,2011.Increase of extreme events in a warming world[J].Proc Natl Acad Sci USA,108(44):17905-17909.doi:10.1073/pnas.1101766108.

Schumacher D L,Keune J,Dirmeyer P,et al.,2022.Drought self-propagation in drylands due to land-atmosphere feedbacks[J].Nat Geosci,15(4):262-268.doi:10.1038/s41561-022-00912-7.

Seneviratne S I,Corti T,Davin E L,et al.,2010.Investigating soil moisture-climate interactions in a changing climate:a review[J].Earth Sci Rev,99(3/4):125-161.doi:10.1016/j.earscirev.2010.02.004.

Seneviratne S I,Zhang X,Adnan M,et al.,2021.Weather and climate extreme events in a changing climate[R]//Masson-Delmotte V,et al.Climate Change 2021:The Physical Science Basis:Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change.Cambridge:Cambridge University Press:1513-1766.doi:10.1017/9781009157896.013.

Shi N,Ekwacu S,Fu S Y,et al.,2023.Features of local finite-amplitude wave activity during extreme cold waves over eastern China[J].J Climate,36(21):7583-7597.doi:10.1175/jcli-d-22-0682.1.

Sun S L,Bi Z Y,Xiao J F,et al.,2023a.A global 5 km monthly potential evapotranspiration dataset (1982—2015) estimated by the Shuttleworth-Wallace model[J].Earth Syst Sci Data,15(11):4849-4876.doi:10.5194/essd-15-4849-2023.

Sun S L,Zhang Y F,Chai R F,et al.,2023b.Spatial differences in impacts of CO2 effect on Chinas meteorological droughts:analysis from surface resistance perspective[J].J Hydrol,621:129568.doi:10.1016/j.jhydrol.2023.129568.

Sun Y,Chen H S,2023.Atmospheric circulation anomalies and key physical processes behind two categories of anomalous Eurasian spring snowmelt[J].J Hydrometeorol,24(8):1349-1363.doi:10.1175/jhm-d-23-0010.1.

Thiery W,Visser A J,Fischer E M,et al.,2020.Warming of hot extremes alleviated by expanding irrigation[J].Nat Commun,11(1):290.doi:10.1038/s41467-019-14075-4.

Tian B Q,Fan K,2022.New downscaling prediction models for spring drought in China[J].Int J Climatol,42(13):6960-6975.doi:10.1002/joc.7623.

Wei J F,Mao Q Q,Shan Y F,et al.,2023.Spring biomass burning in Indochina enhances summer Yangtze River Valley rainfall through land-atmosphere interactions[J].NPJ Clim Atmos Sci,6:183.doi:10.1038/s41612-023-00514-z.

Wouters H,Keune J,Petrova I Y,et al.,2022.Soil drought can mitigate deadly heat stress thanks to a reduction of air humidity[J].Sci Adv,8(1):eabe6653.doi:10.1126/sciadv.abe6653.

Xie J B,Zhang M H,2017.Role of internal atmospheric variability in the 2015 extreme winter climate over the North American continent[J].Geophys Res Lett,44(5):2464-2471.doi:10.1002/2017gl072772.

Xie J H,Hsu P C,Hu Y M,et al.,2023.Skillful extended-range forecast of rainfall and extreme events in East China based on deep learning[J].Wea Forecasting,38(3):467-486.doi:10.1175/WAF-D-22-0132.1.

Xu T B,Yin Z C,Ma X Q,et al.,2023.Hybrid seasonal prediction of meridional temperature gradient associated with “Warm Arctic-Cold Eurasia”[J].Adv Atmos Sci,40(9):1649-1661.doi:10.1007/s00376-023-2226-3.

Xu X K,Huang A N,Huang D Q,et al.,2023.What are the dominant synoptic patterns leading to the summer regional hourly extreme precipitation events over central-eastern Tibetan Plateau and Sichuan Basin?[J].Geophys Res Lett,50(5):e2022GL102342.doi:10.1029/2022gl102342.

Yin Z C,Zhou B T,Duan M K,et al.,2023.Climate extremes become increasingly fierce in China[J].The Innovation,4(2):100406.doi:10.1016/j.xinn.2023.100406.

Zeng J W,Huang A N,Wu P L,et al.,2023.Typical synoptic patterns responsible for summer regional hourly extreme precipitation events over the middle and lower Yangtze River Basin,China[J].Geophys Res Lett,50(17):e2023GL104829.doi:10.1029/2023gl104829.

張德宽,杨贤为,邹旭恺,2003.均生函数-最优子集回归在高温极值预测中的应用[J].气象,29(4):44-47. Zhang D K,Yang X W,Zou X K,2003.Application of mean generating function-optimal subset regression to the prediction of high temperature extremes[J].Meteor Mon,29(4):44-47.doi:10.3969/j.issn.1000-0526.2003.04.010.(in Chinese).

Zhang J,Chen H S,Fang X Y,et al.,2024.Warming-induced hydrothermal anomaly over the Earths three Poles amplifies concurrent extremes in 2022[J].NPJ Clim Atmos Sci.doi:10.1038/s41612-023-00553-6.

Zhang T T,Feng Y Y,Chen H S,2023a.Revealing the formation of the dipole mode of Eurasian snow cover variability during late autumn[J].J Geophys Res:Atmos,128(6):e2022jd038233.doi:10.1029/2022jd038233.

Zhang T T,Zhu S G,Song Y M,et al.,2023b.The leading modes of northern Eurasian winter snowfall variability and the potential influencing factors[J].J Climate,36(22):7811-7826.doi:10.1175/jcli-d-22-0881.1.

中国气象局气候变化中心,2022.中国气候变化蓝皮书(2022)[M].北京:科学出版社. CMA Climate Change Centre,2022.Blue book on climate change in China(2022)[M].Beijing:Science Press.(in Chinese).

Zscheischler J,Westra S,van den Hurk B J J M,et al.,2018.Future climate risk from compound events[J].Nature Climate Change,8:469-477.doi:10.1038/s41558-018-0156-3.

Zscheischler J,Martius O,Westra S,et al.,2020.A typology of compound weather and climate events[J].Nat Rev Earth Environ,1:333-347.doi:10.1038/s43017-020-0060-z.

·ARTICLE·

Research on weather and climate extremes over China:brief introduction and recent progress of the National Key R&D Program of China for Earth System and Global Change

CHEN Haishan1,2,ZHANG Yaocun3,ZHANG Wenjun1,2,YIN Zhicong1,2,HUA Wenjian1,2,KUANG Xueyuan3,CHEN Guosen1,2,MA Hongyun1,2,HAN Tingting1,2

1Key Laboratory of Meteorological Disaster,Ministry of Education (KLME)/Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disaster (CIC-FEMD),Nanjing University of Information Science and Technology,Nanjing 210044,China;

2School of Atmospheric Sciences,Nanjing University of Information Science and Technology,Nanjing 210044,China;

3School of Atmospheric Sciences,Nanjing University,Nanjing 210093,China

Abstract Under the background of global warming,extreme weather and climate events may be unprecedented when they occur with increased frequency,evident temporal and spatial clustering (concurrence),stronger persistence and new combinations (compound).This increases the difficulty of predicting extreme events.Extreme weather and climate events,such as persistent heavy rainfall,extreme low temperature,extreme high temperature and drought compound events,heat waves and strong typhoon,have produced a significant influence on economic,social and sustainable development.However,new features of extreme weather and climate events,coupled with associated key processes and relevant mechanism,remain unclear,and the prediction of extreme weather and climate events needs further improvement.This paper first gives a brief introduction of the National Key R&D Program of China for Earth System and Global Change;supported by this program,the new features of extreme weather and climate events causing serious influence over China under changing climate are investigated.Next,the paper explores the physical mechanisms of ocean-land-atmosphere coupling processes that affect extreme weather and climate events,and attempts to seek the precursors for sub-seasonal to seasonal predictions of extreme events.This program aims to develop a new statistical-dynamical combined prediction method of extreme events and establish a new generation of high-resolution numerical prediction in China,along with a detection and attribution system of extreme events.This paper also highlights recent main achievements of this program.

Keywords extreme weather and climate events;air-sea interaction;land surface process;ocean-land-atmosphere coupling;mechanism;prediction;attribution

doi:10.13878/j.cnki.dqkxxb.20240110007

(責任编辑:倪东鸿)

猜你喜欢

归因机理预测
无可预测
错解归因 寻根溯源
选修2-2期中考试预测卷(A卷)
选修2-2期中考试预测卷(B卷)
隔热纤维材料的隔热机理及其应用
煤层气吸附-解吸机理再认识
积极归因方式:当船漏水了
初中数学学习中教师归因与学生自我归因的比较研究
不必预测未来,只需把握现在
雾霾机理之问