第２４卷第５期 ２００８年１０月 气象与环境学报 ＪＯＵＲＮＡＬ ＯＦ Ｍ哐ＴＥＯＲＯＬＯＧＹ Ａ ＮＤ ＥＮＶＩＲｏＮＭ旺ＮＴ Ｖｏ１．２４ Ｎｏ．５ Ｏｃｔｏｂｅｒ ２００８ １９８１－－－２００５年中国东北榆树物候变化特征及模拟研究 李荣平 ’２’ 周广胜 ， 郭春明４ 孙守军 （１．中国科学院植物研究所植被与环境变化重点实验室，北京，１０００９３；２．中国科学院研究生院，北京，１０００３９； ３．中国气象局沈阳大气环境研究所，辽宁沈阳，１１００１６；４．吉林省气象科学研究所，吉林长春，１３００６１； ５．黑龙江省气象科学研究所，黑龙江哈尔滨，１５００３０） 摘要：探讨了１９８１－２００５年中国东北榆树物候的时空变化与气候因子的关系。结果表明：近２５ ａ，东北地区１３个站点 榆树的展叶期和枯黄期变化趋势较大，１３个物候时间系列，展叶初期显著提前趋势占５４％，枯黄初期显著推后趋势占３８％；榆 纣物候与气温呈显著相关关系：展叶期与３—４月平均气温呈显著的负相关关系，与平均展叶期的积温呈显著的负相关关系；各 站点枯黄期与８—９月平均气温关系不显著，整个东北地区呈显著的正相关关系。研究表明，榆树物候与降水、日照关系不显 著。以日平均气温最高日为界，枯黄期积温值可分为呈线性关系的两段，从日平均气温最高日开始，连续１０ ｄ积温值低于固定 值时，植物开始枯黄。 关键词：中国东北；榆树物候；生长期；降水；气温 中图分类号：Ｑ９４８．１１２ 文献标识码：Ａ 文章编号：１６７３—５０３Ｘ（２００８）０５—００２０—０５ １引言 气温最高日为界限，枯黄期的积温值分为２部分，且 这２部分存在线性关系；（２）从１３平均气温最高日开 始，连续１０ ｄ的积温值小于一个定值，则植物开始枯 还探讨了榆树物候的时空变化特征及其 植物物候是植物发育（如萌芽、展叶、开花和落 黄。同时，叶期）的季节性发生规律。众所周知，植物物候对短 与气候因子的关系，为陆地生态系统模型提供物候 期和长期的气候变化很敏感，是气候变化的敏感指 参数。示器【卜２ｌ。近几十年来，全球气候变化显著，有关植 ２资料与方法 ２．１数据 物物候与气候因子的关系研究日益受到重视。葛全 胜等｜３Ｊ对中国１９６１－－１９９９年气候和物候的关系研 究发现，物候与气温显著相关；物候期的开始日期与 物候数据由辽宁、吉林、黑龙江省气象局提供。 根据观测序列长、连续性较好、覆盖面广的原则，分 别选取１９８１－－２００５年１３个东北地区农业气象站 （表１）的榆树物候观测资料，物候期为展叶初期和枯 黄初期。气候数据为与物候观测年份段相对应的ｌ３ 其前２—３个月的气温有显著的相关关系＿４ ；气温的 高低决定了中纬度植物的春季物候＿５ ；桃花的始花 期与开花前２—３旬的平均积温和平均最高积温关 系密切ｌ ；３月和４月平均气温显著影响芦苇的展 叶［７ ；大豆生长期主要受积温和降水的影响，出苗主 要受积温和日照的影响ｌ＿８ ；２月平均气温是影响银杏 花期的主要因子，花期至成熟期的积温又是影响银 个站点的日平均气温、日降水量和日照时数。 ２．２研究方法 采用儒略日（Ｊｕｌｉａｎ ｄａｙｓ）换算方法，将物候出现 日期转化为距１月１日的实际Ｅｔ数，即年序歹ｆＪ累积日 数，得到各物候期的时间序列。利用ＳＰＳＳ统计软件 杏果实的主要气候因子Ｅ９ ；地中海针叶林物候主要 受影响叶面积发育的水分驱动 。然而，在区域尺 进行统计分析，利用ＦＯＲＴＲＡＮ语言编程进行数据 度上，针对某一特定植物种的整个物候期对气候因 处理。 子响应的研究较少，特别是对植物枯黄期的模拟很 枯黄期模拟的方法。 本文以榆树（Ｕｌｍｕｓ ｐｕｍｉｌａ）整个物候期作为研 究对象，试图证明以下２个假设：（１）１ ａ中的日平均 模拟模型采用积温模型模拟榆树的展叶初期， （１）榆树展叶初期模拟模型 Ｓ＝ （Ｔ，一Ｔｔ，） Ｔ ＞Ｔｈ （１） 少，在国际物候学研究中还未出现普遍认同的植物 利用分段积温来模拟榆树的枯黄初期。 收稿日期：２００８—０５—１１；修订日期：２００８—０７—１７。 基金项目：中国气象局沈阳大气环境研究所生态环境重点实验室项目“东北主要植物物候监测预报产品与模式研发”资助。 作者简介：李荣平，男，１９７５年生，助理研究员，主要从事气候变化方面的研究，Ｅ—ｍａｉｌ：ｒｏｎｇｐｉｎｇｌｉ＠１６３．ｃｏｒｎ。 通信作者：周广胜，Ｅ－ｍａｉｌ：ｇｓｚｈｏｕ＠ｉｅａｓ．ａｂＣ．Ｃｒｌ 第５期 李荣平等：１９８１－－２００５年中国东北榆树物候变化特征及模拟研究 ２１ ———一———一——————一式（１）中，Ｓ为积温阈值（℃・ｄ）；７＂ｉ为日平均气温 表１物候观测站点和观测年份段 站点 所辖省份纬度尸Ｎ经度尸】 海拔／ｍ观测年份段 朝阳 辽宁 ４１．５５ １２０．４５ １７６．０ １９８１－－２００５ 开原 辽宁 ４２．５３ １２４．０５ ９９．０ １９８１－－２００５ 建平 辽宁 ４１．８７ １１９．６３ ６６１．６ １９８１－－２００５ 新宾 辽宁 ４１．７３ １２５．０５ ３２８．７ １９８２－－２００５ （℃）； 为基础温度（℃）；Ｙ为展叶期（ｄ）。 （２）榆树枯黄初期模拟模型 假设将植物的积温过程分成２段，第１段为积温 速率由小到大，第２段为积温速率由大到小，以榆树 的枯黄初期为终点，以积温速率最大时为分界点，这 ２段积温存在极显著的线性关系，因为从积温速率最 大时开始，环境温度是逐渐降低的。采用窗口移动 方法，从积温速率最大时开始移，当窗口时间范围内 的积温小于一个阈值的时候，植物开始枯黄初期。 因此，选取窗口时间尺度为１０ ｄ，则： Ｓ ＝ ：∑ Ｔ Ｔ ＞０℃ Ｔ＞０℃ （（）２） ５ ｒ １０ 式（２）中，Ｓ为连续１０ ｄ大于０℃的积温（℃・ｄ）；ｖ 为枯黄期（ｄ）。 ３结果分析 ３．１榆树的物候期特征 ３．１．１时间特征 分析表明，１９８７－－２００４年，中国东北地区榆树平 均展叶初期显著提前，每１０ ａ提前４．１ ｄ。１９８１— ２００５年，东北地区ｌ３个榆树物候观测站中展叶期提 前的站点占６９％，其中显著提前占５４％。朝阳、建平 和黑河表现为极显著的提前趋势，每１０ ａ分别提前 ５．７、９．１ ｄ和６．２ ｄ；开原、梅河口、松原和泰来表现 为显著提前趋势，每１０ ａ分别提前２．４、４．０、２．５ ｄ和 ３．６ ｄ；其余的站点表现不显著（表１）。这与欧洲物 候特征相似，１９７１－－２０００年欧洲５４２种植物７８％的 展叶、开花和结实提前，其中３０％显著提前，３％表现 为显著推后ｌｌ１一。在德国，１９８４－－２００４年植物春季物 候的２３个系列中２２个为物候显著提前，其中６个表 现为每１０ ａ提前ｉ０．０ ｄ Ｌｌ２一。 一——————————— ———————～————————————一 １９８７－－２００５年，东北地区榆树平均枯黄期变化 趋势不显著。１９８１－－２００５年，东北地区１３个榆树物 候站点中，桦甸、榆树和穆棱表现为极显著推后，每 １０ ａ分别推后７．２、７．２ ｄ和６．９ ｄ；朝阳和建平表现 为显著推后，每１０ ａ分别推后４．９ ｄ和５．８ ｄ；其余表 现为变化趋势不显著（表２）。 表２ １９８１－－２００５年东北地区榆树物候期变化 注：Ｋ为变化斜率；尺为物候期与时问的相关系数；Ｎ为 样本数；＊为Ｐ＜０．０５；＊＊为Ｐ＜０．０１。 以上分析结果表明，东北地区榆树展叶初期存 在显著提前趋势，枯黄初期存在显著退后趋势的现 象，整个生长季延长。 ３．１．２空间特征 榆树物候在空间上呈现明显的变化趋势。从中 国辽宁至吉林至黑龙江省，随着纬度的增加，榆树展 叶期是逐渐推后的，这也再一次证明榆树的展叶期 与气温呈负相关关系。同时，榆树的枯黄期随着纬 度的升高呈显著的提前趋势，温度的降低导致榆树 枯黄期的提前。生长期与枯黄期在空问上的变化相 似，越向北，榆树的生长期越短，生长期随着纬度的 增加而呈缩短趋势（图１）。Ａ．ｈａｓ等 Ｆ Ｊ指出，在欧洲 ８个地区纬度跨度４７．２２。一５７．８８。，桦树在欧洲的南 部４月４日就开始展叶，在欧洲的北部５月３０日才 展叶。 ３．２榆树物候期与气候因子的关系 分析榆树展叶初期与３—４月的日平气温、日平 均日照时数、降水量和１—１２０ ｄ的积温（大于０℃） 的相关关系，榆树展叶初期仅仅与３—４月的日平均 气温和１—１２０ ｄ的积温存在显著的负相关关系，其 他关系不显著。分析榆树枯黄初期与８—９月的日 平均气温、日平均日照时数、降水量和年积温（大于 第５期 李荣平等：【９８１一１００５年中国东北榆树物候变化特征及模拟研究 个站点榆树展叶初期趋势不同 但主要为提前趋 表－ｌ东北地区榆树枯黄期的第１段积温（ ）与 第２段积温（Ｙ）的关系 站点 朝阳 桦甸 Ｙ一Ｙ＝ｌ 一Ｙ一势，研究表明＿Ｉ 植物展叶是与展叶初期的温度具有 显著的负相关关系，本研究也发现了同样的规律（Ｒ ＝一０．６１ ～一０．８１一），因此，榆树展叶初期的特 征反映了展叶初期前的温度变化特征 榆树的枯黄 初期在各站点也变化不同，但是榆树的枯黄期对温 度的响喧不是很清楚。研究发现，植物枯黄期的变 化是模糊的，在ｌ４个欧洲国家，秋天叶变色和落叶， 植物物候与气温趋势的关系不是很清楚ｌ】 ，本研究 ∞ ∞ 如 ∞ 蛐 ∞ ∞ 加 方程式 ０．８７９ｌ ＋１７ｌ５．４ ０．８４９２ ＋ｌ９７２．６ ０．９２６７ 十ｌ７０５．６ ０．８９ｌ７ ＋２ｌ５１．Ｏ 相关系数（Ｒ） ０．９５３８“ ０．８０６４… —０．９５７８ ０．９５ｌ８“ 黑河 开原 建平 Ｙ＝０．９２１６ ＋２４７５．０ ０．９３８８… 梅河口 Ｙ一０．９１８５ ＋２１４６．０ ０．９２４７一 结果与之相同。 穆棱 Ｙ一０．９７０２ ＋ｌ９６９．７ ０．８７６６… （２）在空间上，东北榆树的物候在纬度上存在明 庆安 ｙ—１．０７１５Ｘ＋２１（１５．８ ０．９３３０一 显的梯度趋势，展叶初期随纬度的增加而展叶推后， 青冈 、 ：一１．１３０９ ＋２ｌ９６．８ —０．９：７１ 枯黄初期随着纬度的增加而展叶提前。这主要是由 松原 Ｙ＝ １．０１【】８ ＋２４ｊ３．６ ０．９４５３… 于在纬度上，气温存在一个明显的梯度关系，随着纬 泰来 ＝一０．９９｛）ｌ ＋２３４１．６ —０．９３ｌ８ ’一 度的增加，气温显著降低，而榆树的物候期显著受气 新宾 ｙ＝０．９７５８ ‘十！ｌ６７．７ ０．８８８９ 榆树 ０ ９５１９．１。 ２２２６．５ —０ ９４５９ 温的影响。研究发现，榆树枯黄初期在单个站点上 ＝东北 Ｙ ＝０．６９２６ 十Ｉ９ｌ３．３ ０．６８４３ 一 并不表现为显著的与８～９月平均气温的相关性，而 注：＊＊＊为Ｐ＜０．００１：第１段积温为１月１日至日平均 在整个东北区域却呈现显著的相关关系，这也说明 气温最高日积温，第 段积温位勾Ｅｔ平均气温最高臼至叶枯 榆树枯黄初期在时间上不显著的推后趋势，而在空 黄初日积温 间上表现为提前趋势。偷树枯黄初期主要受一段时 ４ ２７ ４０ ５３ ６６ ７９ ９２ １（１５ ｌ１ ８ １ ３１ Ｉ４４ １ ５ ７ １ ７０ １ ８３ １９６ ２０９ ２２２ ２３ ５ ２４８ ２６ｌ ２７４ ２８ ７ 样本数／个 －……观 则随：—— 测值 图３ ｌ９８ｌ一２ＯＯ５年东北地区榆树枯黄初期的观测值与预测值变化曲线 间内的温度持续偏低而引起一当一段时间温度累加 低于某一个值的时候，榆树枯黄初期也就开始，而与 参考文献 某段固定时间的平均气温的高低无必然联系 １ Ｂａｄｅｃｋ Ｆ Ｗ，Ｉ３０ＩｌｄｅａＨ Ａ，Ｂ６ｔｔｃｈｅｒ Ｋ，ｅｔ ａ１．Ｒｅｓｔｘ）ｎｓｍ ｏｆ （３）积温模型可以较好地模拟榆树展叶初期，对 ｓｐｒｉｎｇ ｐｈｏｎｏｌｏｇｙ ｔｏ ｃｌｉｍａｔｅ ｃｈａｎｇｅ［］ ．Ｎｅｗ Ｐｈｙｔｏｌｏｇｉｓｔ， ２００４，１６２（２）：２９５—３０９． ２９０个观测数据的验证，解释率达６２％，关系极显 ２］Ｃｈｕｉｎｅ Ｉ，Ｙｉｏｕ Ｉ ，Ｖｉｏｖｙ Ｎ，ｅｔ ａ１．Ｇｒａｐｅ ｒｉｐｅｎｉｎｇ ａｓ ａ ｐａｓｔ 著 植物枯黄期的模拟至今仍然不存在一种令科学 ｄｉｎｌａｔｅ ｉｎｄｉｃａｔｏｒ －】Ｊ．Ｎａｔｕｒｅ，２００４，４３２（７０１５）：２８９— 家普遍接受的方程，在过去的研究中，出现丁不同的 ２９０． 研究方法，但是结果都较差。本研究所建物候模型 ３］葛全胜，郑景云，张学霞．近４０年中国气候与物候的 的解释率达８３％，尽管偏差为１２．７ ｄ，但是与其他模 变化研究［Ｊ］．自然科学进展，２００３，１３（１０）：１０４８ １０５３． 拟枯黄期方法相比，具有简单和相对准确的特点。 ４ Ｊ Ａｈａｓ Ｒ，Ｊａａｇｕｓ Ｊ，Ａａ Ａ． Ｆｈｅ ｐｈｅｎｏｌｏｇｉｃａｌ ｃａｌｅｎｄａｒ ｏｆ （４）对榆树展叶初期的３８％和枯黄初期的ｌ７ ｔ圮 Ｆ￣ｓｔｏｎｉａ ａｎｄ ｉｔｓ ｃｏｒｒｅｌａｔｉｏｎ ｗｉｔｈ ｍｅａｎ ａｉｒ ｔｅｍｐｅｒａｍｒｅ 的变化模型不能解释，主要原因：（１）数据来源于不 Ｊ ．Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｏｆ Ｂｉｏｍｅｔｅｏｒｏｌｇｙ，２０００．４４ 同的站点，观测记录的人不同，存在较大的人为误 （４）：１５９—１６６． ５ Ｊ Ｃｈｍｉｅｌｅｗｓｋ Ｆ Ｍ．Ｒｏｔｚｅｒ Ｔ．Ｒｅｓｐｏｎｓｅ ｏｆ ｔｒｅｅ ｐｈｅｎｏｌｏｇｙ 差；（２）在不同地点的榆树物种可能存在基因变化， ｔｏ ｃｌｉｍａｔｅ ｃｈａｎｇｅ ａｃｒｏｓｓ Ｅｕｒｏｐｅ Ｊ］．Ａｇｒｉｃｕｌｔｕｒａｌ ａｎｄ 造成一部分不能解释的变化。 Ｆｏｒｅｓｔ Ｍｅｌｅｏｒｏｌｏｇ５￣，２００１，１０８（２）：ｌ０ｌ １ １２ ２４ 气象与环境学报 第２４卷 张秀英，胡东燕．桃花花期预报的探讨［Ｊ］．北京林业大 Ｂｉｏｍｅｔｅｏｒｏｌｏｇｙ，２０００，４４（２）：７６—８１． 学学报，１９９５，１７（４）：８８—９２． Ｍｅｎｚｅｌ Ａ，Ｓｐａｒｋｓ Ｔ，Ｅｓｔｒｅｌｌａ Ｎ，ｅｔ ａ１．Ｅｕｒｏｐｅａｎ ｐｈｅｎｏ— 李荣平，刘晓梅，周广胜．盘锦湿地芦苇物候特征及其 ｌｏｇｉｃａｌ ｒｅｓｐｏｎｓｅ ｔｏ ｃｌｉｍａｔｅ ｃｈａｎｇｅ ｍａｔｃｈｅｓ ｔｈｅ ｗａｒｍｉｎｇ 对气候变化的响应［Ｊ］．气象与环境学报，２００６，２２（４）： ｐａｔｔｅｒｎ［Ｊ］．Ｇｌｏｂｄ Ｃｈａｎｇｅ Ｂｉｏｌｏｇｙ，２００６（１２）：１９６９— ３０—３４． １９７６． 刘景利，杨扬，史奎桥，等．１９８５－－２００５年锦州地区大 Ｓｃｈａｂｅｒ Ｊ．Ｐｈｅｎｏｌｏｇｙ ｉｎ ｇｅｒｍａｎｙ ｉｎ ｔｈｅ ２０ｔｈ ｃｅｎｔｕｒｙ： 豆物候期变化及其气候响应［Ｊ］．气象与环境学报， ｍｅｒｔｈｏｄｓ，ａｎａｌｙｓｅｓ ａｎｄ ｍｏｄｅｌｓ［Ｄ］．Ｇｅｒｍａｎｙ：Ｕｎｉｖｅｒｓｉｔｙ ２０Ｏ７，２３（４）：２９—３２． ｏｆ Ｐｏｓｔｓｄａｒｎ． 冯颖竹，梁红．广东省南雄、和平银杏的物候特征及其 Ａｈａｓ Ｒ，Ａａｓａ，Ａ，Ｍｅｎｚｅｌ Ａ，ｅｔ ａ１．Ｃｈａｎｇｅｓ ｉｎ ｅｕｒｏｐｅａｎ 与气候因子的关系［Ｊ］．植物资源与环境，１９９９，８（３）： ｓｐｒｉｎｇ ｐｈｅｎｏｌｏｇｙ［Ｊ］．Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｏｆ Ｃｌｉｍａｔｏｌｏ— ３４—３８． ｇｙ，２００２（２２）：１７２７—１７３８． Ｋａｍｅｒ Ｋ，Ｌｅｉｎｏｎｅｎ Ｉ，Ｌｏｕｓｔａｕ Ｄ．Ｔｈｅ ｉｍｐｏｒｔａｎｃｅ ｏｆ Ｃｈｕｉｎｅ Ｉ，Ｃａｍｂｏｎ Ｇ，Ｃｏｍｐｔｏｉｓ Ｐ．Ｓｃａｌｉｎｇ ｐｈｅｎｏｌｇｏｙ ｐｈｅｎｏｌｏｇｙ ｆｏｒ ｔｈｅ ｅｖａｌｕａｔｉｏｎ ｏｆ ｉｍｐａｃｔ ｏｆ ｃｌｉｍａｔｅ ｃｈａｎｇｅ ｆｒｏｍ ｔｈｅ ｌｏｃａｌ ｔｏ ｔｈｅ ｒｅｇｉｏｎａｌ ｌｅｖｅｌ：ａｄｖａｎｃｅｓ ｆｒｏｍ ｓｐｅｃｉｅｓ ｏｎ ｇｒｏｗｔｈ ｏｆ ｂｏｒｅａｌ，ｔｅｍｐｅｒａｔｕｒｅ ａｎｄ Ｍｅｄｉｔｅｒｒａｎｅｍ￣ｆｏｒ— —ｓｐｅｃｉｆｉｃ ｐｈｅｎｏｌｇｏｉｃＮ ｍｏｄｅｌｓ［Ｊ］．Ｇｌｏｂａｌ Ｃｈａｎｇｅ Ｂｉｏｌｏ— ｅｓｔ ｅｃｏｓｙｓｔｅｍｓ：Ａｎ ｏｖｅｒｖｉｅｗ［Ｊ］．Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｏｆ ｇＹ，２０００（６）：９４３—９５２． Ｐｈｅｎｏｌｏｇｙ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｕｌｍｕｓ ｐｕｍｉｌａ ａｎｄ ｉｔｓ ｓｉｍｕｌａｔｉｏｎ ｆｒｏｍ １９８１ ｔｏ ２００５ ｉｎ Ｎｏｒｔｈｅａｓｔ Ｃｈｉｎａ ＬＩ Ｒｏｎｇ—ｐｉｎｇ ， ， ＺＨＯＵ Ｇｕａｎｇ ｓｈｅｎｇ ＇。 ＧＵＯ Ｃｈｕｎ—ｍｉｎｇ４ ＳＵＮ Ｓｈｏｕ—ｊｕｎｓ （１．Ｋｅｙ Ｌａｂｏｒａｔｏｒｙ ｏｆ Ｖｅｇｅｔａｔｉｏｎ ａｎｄ Ｅｎｖｉｒｏｎｍｅｎｔａｌ Ｃｈａｎｇｅ，Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｂｏｔａｎｙ，ｔｈｅ Ｃｈｉｎｅｓｅ Ａｃａｄｅｍｙ ｏｆ Ｓｃｉｅｎｃｅｓ，Ｂｅｒｉｎｇ １０００９３，Ｃｈｉｎａ；２．Ｇｒａｄｕａｔｅ Ｓｃｈｏｏｌ ｏｆ Ｃｈｉｎｅｓｅ Ａｃａｄｅｍｙ ｏｆ Ｓｃｉｅｎｃｅｓ，Ｂｅｉｊｉｎｇ １０００３９， Ｃｈｉｎａ；３．Ｉｎｓｔｉｔｕｔｅ ｏｆ Ａｔｍｏｓｐｈｅｒｉｃ Ｅｎｖｉｒｏｎｍｅｎｔ，Ｃｈｉｎａ Ｍｅｔｅｏｒｏｌｏｇｉｃａｌ Ａｄｍｉｎｉｓｔｒａｔｉｏｎ，Ｓｈｅｎｙａｎｇ １１００１６，Ｃｈｉｎａ；４．Ｊｉｌｉｎ Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｍｅｔｏｅｒｏｌｏｇｉｃａｌ Ｓｃｉｅｎｃｅｓ，Ｃｈａｎｇｃｈｕｎ １３００６１，Ｃｈｉｎａ； ５．Ｈｅｉｌｏｎｇｊｉａｎｇ Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｍｅｔｏｅｒｏｌｏｇｉｃａｌ Ｓｃｉｅｎｃｅｓ，Ｈａｒｒｂｉｎ １５００３０，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ｔｈｅ ｔｅｍｐｏｒａｌ ａｎｄ ｓｐａｔｉａｌ ｃｈａｎｇｅｓ ｏｆ ＵＺＩＥｌＵＳ ｐｕｍｉｌａ ｐｈｅｎｏｌｏｇｙ ａｎｄ ｔｈｅ ｒｅｌａｔｉｏｎｓｈｉｐ ｗｉｔｈ ｃｌｉｍａｔｅ ｆａｃｔｏｒｓ ｉｎ Ｎｏｒｔｈｅａｓｔ Ｃｈｉｎａ ｆｒｏｍ １９８１ ｔｏ ２００５ ｗｅｒｅ ｄｉｓｃｕｓｓｅｄ．Ｔｈｅ ｒｅｓｕｌｔｓ ｓｈｏｗ ｔｈａｔ ｔｈｅ ｃｈａｎｇｅ ｔｒｅｎｄｓ ｏｆ ｌｃａｖｅ ｕｎｆｏｌｄｉｎｇ ｐｈａｓｅ ａｎｄ ｌｃａｖｅ ｃｏｌｏｒｉｎｇ ｐｈａｓｅ ｏｆ ｕｌｌｎｕｓ ｐｕｍｉｌａ ａｒｅ ｏｂｖｉｏｕｓ ｄｕｒｉｎｇ ｒｅｃｅｎｔ ２５ ｙｅａｒｓ ｂａｓｅｄ ｏｎ ｔｈｅ ｄａｔａ ｆｒｏｍ １３ ｏｂ— ｓｅｒｖａｔｉｏｎ ｓｔａｔｉｏｎｓ ｉｎ Ｎｏｒｔｈｅａｓｔ Ｃｈｉｎａ．Ｔｉｍｅ ｓｅｑｕｅｎｃｅ ｎｕｍｂｅｒｓ ｔｈａｔ ｌｅａｖｅ ｕｎｆｏｌｄｉｎｇ ｐｈａｓｅ ｉｓ ｓｉｇｎｉｆｉｃａｎｔｌｙ ａｄｖａｎｃｅｄ ａｃｃｏｕｎｔ ｆｏｒ ５４％ｏｆ ｔｈｅ ｔｏｔａｌ ｌｅａｖｅ ｕｎｆｏｌｄｉｎｇ ｔｉｍｅ ｓｅｑｕｅｎｃｅ ｎｕｍｂｅｒｓ，ａｎｄ ｔｈｏｓｅ ｔｈａｔ ｌｃａｖｅ ｃｏ［ｏｒｉｎｇ ｐｈａｓｅ ｉｓ ｓｉｇｎｉｆｉ— ｃａｎｔｌｙ ｄｅｌａｙｅｄ ａｃｃｏｕｎｔ ｆｏｒ ３８％．Ｔｈｅ ｒｅｌａｔｉｏｎｓｈｉｐ ｂｅｔｗｅｅｎ ｕｌｍｕｓ ｐｕｍｉｌａ ｐｈｅｎｏｌｏｇｙ ｐｈａｓｅｓ ａｎｄ ａｉｒ ｔｅｍｐｅｒａｔｕｒｅ ｉｓ ｓｉｇｎｉｆｉｃａｎｔ．Ｔｈｅ ｌａｅｖｅ ｕｎｆｏｌｄｉｎｇ ｐｈａｓｅ ｈａｓ ｎｅｇａｔｉｖｅ ｒｅｌａｔｉｏｎｓｈｉｐｓ ｗｉｔｈ ｍｅａｎ ａｉｒ ｔｅｍｐｅｒａｔｕｒｅ ｏｆ Ｍａｒｃｈ ａｎｄ Ａｐｒｉｌ ａｎｄ ｉｔｓ ｃｕｍｕｌａｔｉｖｅ ｔｅｍｐｅｒａｔｕｒｅ．Ｔｈｅ ｒｅｌａｔｉｏｎｓｈｉｐ ｂｅｔｗｅｅｎ １ｅ．ａｖｅ ｃｏｌｏｒｉｎｇ ｐｈａｓｅ ｉｎ ｅａｃｈ ｓｔａｔｉｏｎ ａｎｄ ｍｅａｎ ａｉｒ ｔｅｍ— Ｄｅｒａｔｕｒｅ ｏｆ Ａｕｇｕｓｔ ａｎｄ Ｓｅｐｔｅｍｂｅｒ ｉｓ ｎｏｔ ｓｉｇｎｉｆｉａｃｎｔ，ｗｈｉｌｅ ｉｔ ｉｓ ａ！ｘ￣ｓｉｔｉｖｅ ｒｅｌａｔｉｏｎ ｉｎ ｔｈｅ ｗｈｏｌｅ ｓｔｕｄｙ ａｒｅａｓ．Ｔｈｅ ｒｅ— ｓｕｌｔｓ ａｌｓｏ ｉｎｄｉｃａｔｅ ｔｈａｔ ｔｈｅ ｐｈｅｎｏｌｏｇｙ ｐｈａｓｅｓ ｏｆ ｕｌｍ ｕｓ ｐｕｍｉｌａ ａｒｅ ｎｏｔ ｓｉｇｎｉｆｉｃａｎｔｌｙ ｒｅｌａｔｅｄ ｔｏ ｐｒｅｃｉｐｉｔａｔｉｏｎ ａｎｄ ｓｕｎｓｈｉｎｅ ｄｕｒａｔｉｏｎ．Ｉｎ ｔｅｒｍｓ ｏｆ ｔｈｅ ｄａｙ ｏｆ ｍａｘｉｍｕｍ ｄａｉｌｙ ｍｅ２ｌｎ ａｉｒ ｔｅｍｐｅｒａｔｕｒｅ，ｔｈｅ ｃｕｍｕｌａｔｉｖｅ ｔｅｍｐｅｒａｔｕｒｅ ｏｆ ｌａｅｖｅ ｃｏｌｏｒｉｎｇ ｐｈａｓｅ ｃａｎ ｂｅ ｄｉｖｉｄｅｄ ｉｎｔｏ ｔｗｏ ｓｔａｇｅｓ．Ｐｌａｎｔ ｓｐｅｃｉｅｓ ｓｔａｒｔ ｔｏ ｂｅｃｏｍｅ ｙｅｌｌｏｗ ｗｈｅｎ ｔｈｅ ｃｕｍｕｌａｔｉｖｅ ｔｅｍ— Ｄｅｒａｔｕｒｅ ｖａｌｕｅ ｏｆ ｃｏｎｔｉｎｕｏｕｓ １０ ｄａｙｓ ｆｒｏｍ ｔｈｅ ｈｉｇｈｅｓｔ ｄａｉｌｙ ｍｅａｎ ａｉｒ ｔｅｍｐｅｒａｔｕｒｅ ｄａｙ ｉｓ ｌｏｗｅｒ ｔｈａｎ ｔｈｅ ｆｉｘｅｄ ｖａｌ— Ｕｅ． Ｋｅｙ ｗｏｒｄｓ：Ｎｏｒｔｈｅａｓｔ Ｃｈｉｎａ；Ｕｌｍｕｓ ｐｕｍｉｌａ ｐｈｅｎｏｌｏｇｙ；Ｇｒｏｗｉｎｇ ｓｔａｇｅ；Ｐｒｅｃｉｐｉｔａｔｉｏｎ；Ａｉｒ ｔｅｍｐｅｒａｔｕｒｅ