2. 中国气象局气象影视中心,北京 100081;
3. 中国科学院大气物理研究所,北京 100029
2. CMA Audio and Video Centre, Beijing 100081;
3. Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029
北京市的旅游资源极为丰富,一年四季都有不少好去处,尤其是秋天的北京最具特色。香山是北京秋季旅游不可错过的景点,其红叶驰名中外,是我国四大赏枫胜地之一,也是“新北京十六景”之一。香山红叶主要有8个科涉及14个树种,总株数达14万株,其中黄栌有10万余株,是香山公园红叶的主体树种。黄栌树叶大概从10月上旬开始慢慢变红,观赏期一直持续到11月上旬甚至中旬,能为公园及周边地区带来可观的经济收益(刘杰等,2010;姚秀萍等,2011)。葛雨萱等(2011)研究指出:黄栌叶片绿色褪去、红色显现的主要原因是叶绿素含量迅速下降以及花青素含量急速上升,发生过程比较迅速。
根据物候学理论,前期气象条件如光、温、水等对植物物候早晚有重要影响(刘建栋等,2002; 刘玲等,2006;竺可桢等,1973),其中气温是影响中国木本植物物候的主要因子,海拔高度以及地形地貌因可影响温度高低有时也被作为影响因子。陈正洪等(2008)的研究指出前期2月、冬季平均气温每升高1℃,武汉的樱花始花期分别提前1.66和2.86 d。黄栌的生长发育及叶片变色与气象条件也存在密切的关系。在春夏季节,由于温度适宜,黄栌的树叶光合作用强烈,叶绿素掩盖着花青素,因此呈现出一片绿色。一到深秋,树叶受到低温霜冻的侵袭,叶绿素被破坏,叶子里的水分减少了,不能及时运输的淀粉变成了葡萄糖,糖分就逐渐转化为花青素,于是绿叶变成了红叶。周肖红等(2009)曾进行科学试验,发现在昼/夜温度为17/9℃循环一周的条件下,黄栌的嫩叶和成熟叶都可以积累高浓度的花青素,从而呈现红色。
以前对香山红叶变色日(CCLD-Day)气象条件的研究多是经验总结(赵颖彦,2011),缺乏利用物候观测资料的科学分析和验证,得到的阈值标准也不统一。本文将以黄栌树叶为代表,分析香山红叶变色日的基本特征,找出影响变色早晚的气象条件阈值和关键气候因子,并建立香山红叶变色日的气象统计预测模型。
1 资料及方法本文采用的物候资料为香山公园红叶物候观测组监测的逐日黄栌树叶变色率,设定全园变色率达到5%~10%的日期为红叶变色日。资料年份为1998、1999、2000、2004、2005、2006、2009、2010、2011和2012年(Phonological Observation Data, POD-10)。方法及机理研究主要使用前9年监测数据,其中前8年数据(POD-8) 用作诊断分析,2011年的监测数据作为独立样本进行验证, 而2012年用作个例预报试验。
根据海拔高度和到香山的距离,本文选用昌平和门头沟气象站的逐日观测平均值作为代表,资料年份为1998—2012年,要素包括最高气温、最低气温、平均气温和降水量等(杨萍等,2011)。本文所采用的方法包括相关分析和回归分析等。
2 黄栌树叶变色的气象条件阈值曾有人根据经验总结出香山红叶变色的阈值条件(赵颖彦,2011),但缺乏科学验证和统一标准。图 1给出了黄栌树叶(POD-8) 变色前10 d最低气温、气温日较差和平均气温的变化曲线,可以发现红叶变色前各温度要素变化非常明显,且离散度非常小,表现出一致的趋势。由此可以得出,黄栌树叶变色的气象条件阈值:(1) 最低气温(Tmin)低于14℃;(2) 气温日较差(DT)大于10℃;(3) 日平均气温(Tmean)在13~20℃;此3项条件保持3 d,第4 d为红叶变色日。
![]() |
图 1 香山黄栌树叶(POD-8) 变色前10 d的最低气温(a)、气温日较差(b)和日平均气温(c) (横坐标0代表树叶变色日,横虚线代表各阈值条件) Fig. 1 The minimum temperature (a), diurnal temperature range (b) and daily mean temperature (c) of the 10 days before CCLD (The "0" on x axis denotes CCLD-Day, and horizontal dashed lines denote threshold) |
香山黄栌树叶变色的物候监测资料非常珍贵,也非常缺乏,POD-9资料并不连续,中间有5年空白,不利于气象统计预测模型的建立。为得到连续的黄栌树叶变色日,采用气象条件阈值对物候资料进行插补。即根据3个气象阈值条件共同约束的指标,得到空白年份的黄栌树叶变色日。从表 1对已有资料年份的回报误差统计效果看,插补方法是可信的,误差在可接受范围内。插补结果如图 2所示,得到1998—2011年共14年连续的黄栌树叶变色日(Phonological Interpolation Data, PID-14)。对PID-14序列进行自相关检验(图略),发现PID-14序列在延迟1年之后,自相关便不再显著,为典型的非自相关序列。
![]() |
图 2 香山黄栌树叶变色日 (空心为插补资料) Fig. 2 The CCLD-days from 1998 to 2011 (Hollow bars are interpolation data) |
黄栌树叶平均变色日的日序数是277.36,对应的日期是10月4—5日。最早的红叶变色日是9月27日(1999年),在9月变色的还有2001年。其他年份的红叶变色日均出现在10月,最晚红叶变色日为10月13日(2006年)。从变化趋势上看,黄栌树叶变色日以2006年为界,存在先延迟后提前的趋势。
4 黄栌变色日与月尺度气象条件的相关性为分析黄栌变色早晚与前期月尺度气象要素之间的关系,表 2给出了变色日与前期气象要素的相关系数。从中可发现树叶变色日与7月的各温度要素均存在显著的负相关,与最低气温、平均气温和最高气温的相关系数均通过0.01的显著性水平检验。由此可以认为前期7月的最高温度和平均气温越高,当年黄栌树叶变色越早。
![]() |
表 2 红叶变色日与前期6—9月气象要素的相关系数 Table 2 Correlation coefficients between CCLD-day (PID-14) and meteorological condition from June to September |
与基本量的相关性不同,黄栌树叶变色日的年际增量与7月降水的年际增量之间存在显著的正相关(表 3),相关系数通过0.01的显著性水平检验,当7月降雨较去年多时,红叶变色期较去年偏晚。充沛的降雨非常有利于树叶的生长和叶绿素的生成,所以导致叶绿素被破坏,叶片绿色褪去所需要的时间相对较长。同时,黄栌树叶变色日的年际增量与平均气温和最高气温之间依然保持显著的负相关性,但是与气温日较差之间有显著的正相关。
![]() |
表 3 红叶变色日年际增量与前期6—9月气象要素年际增量的相关系数 Table 3 Correlation coefficient between annual increment of CCLD-day (PID-14) and that of meteorological condition from June to September |
年际增量在进行短期气候预测时具有放大对象预测信号的优势(范可等,2007;2008)。黄栌树叶变色日基本量及其年际增量的方差分别为19.2和33.2,变换为年际增量后,信号被放大了70%以上。根据以上分析,在进行短期气候预测时,以变色日的年际增量为预测对象。经过多次预报试验,选择7月平均温度和最高温度的年际增量作为预报因子,建立回归预测模型(张德山等,2010)如下:
$Y{\rm{ = 0}}{\rm{.01 + 6}}{\rm{.8320}} \times \Delta {T_{{\rm{mean}}}} - 8.4392 \times \Delta {T_{{\rm{max}}}}$ | (1) |
图 3给出了回报与PID-14两条曲线,相关系数达到0.73, 通过α=0.01的显著性水平检验,1999—2010年回报的平均误差为3 d左右,2011年作为独立样本的预报变色日序为275.75,与实际的变色日(日序:275) 非常接近。由此可以认为香山红叶变色日短期气候预测模型预测结果可信,模型具有一定的实用价值。
![]() |
图 3 1999—2011年回报的变色日(虚线)及PID-14(实线)变色日 Fig. 3 Forecast results used former data of CCLD-day (dashed line) and PID-14 series (solid line) from 1999 to 2011 |
根据短期气候预测的结果,划定进行中期预报的时间段,一般应该比短期气候预测值提前15 d开始根据气象条件阈值进行中期预报。温度条件可来自于BJ-RUC精细化数值预报和北京市气象台每天发布的未来7 d天气预报,也可以是两者的综合结果。预报时应采取迭代预报、逐步订正的方式进行。
5.3 2012年预报试用香山公园红叶物候观测小组的观测资料显示:由于夏季雨水丰沛,人工养护力度加大,2012年红叶长势和观赏效果要好于往年,呈现出叶色红、润、艳的特点。由观测的逐日红叶变色率,确定10月6日为香山红叶变色日。2012年7月昌平和门头沟两站平均的最高气温和平均气温分别为31.2℃和26.7℃。根据香山红叶变色日气象统计模型(式1),采用平均气温和最高气温的年际增量作为预报因子进行短期气候预测,预测的变色日序为278.7,变色日为10月6日左右,与观测结果(日序:279) 相比,误差可以忽略。
根据短期气候预测的变色日,提前15 d(9月21日)进行中期预报,温度曲线如图 4所示。9月29日至10月1日的温度条件达到阈值要求,第一次中期预报的变色日为10月2日,误差为-4 d。按照模型的约定,中期预报采取迭代预报,逐步订正的方式进行。10月4—6日的温度条件再次达到阈值条件要求,第二次预报的变色日为10月7日,误差为1 d,比较接近观测结果。
![]() |
图 4 2012年9月21日至10月14日两站平均的最低气温(a)、日较差(b)和平均气温(c) (虚线为对应的阈值条件) Fig. 4 Minimum temperature (a), diurnal temperature range (b) and daily mean temperature (c) from 21 September to 14 October 2012 (Horizontal dashed lines denote threshold) |
由2011年独立样本检验和2012年预报检验的结果,可以认为香山红叶变色日气象统计模型具有比较高的准确率和可操作性。当然需要注意的是,虽然2012年短期气候预测的结果好于中期预报,但并不能说明两者之间的优劣。在使用过程中,应当根据实际情况,综合使用两种不同时间精度的预报方法。
6 结论与讨论本文采用香山公园逐日黄栌树叶变色率以及相应的气象资料研究了黄栌树叶变色的气象条件,并建立了香山红叶变色日的气象统计预测模型。
(1) 黄栌树叶变色的气象条件阈值:最低气温低于14℃;日较差在10℃以上;日平均气温在13~20℃之间;此3项条件保持3 d,第4 d为红叶变色日。
(2) 黄栌树叶平均变色日的日序数是277.36,对应的日期是10月4—5日。最早的红叶变色日是9月27日(1999年),最晚红叶变色日为10月13日(2006年)。黄栌树叶变色日以2006年为界,存在先延迟后提前的趋势。
(3) 黄栌树叶变色日与7月最低气温、平均气温和最高气温均存在显著的负相关。与基本量的相关性不同,黄栌树叶变色日的年际增量与7月降水的年际增量之间存在显著的正相关。
(4) 选择7月的平均温度和最高温度的年际增量作为预报因子,建立模型,1999—2010年回报的平均误差为3 d左右。同时,根据短期气候预测的结果,可以利用气象条件阈值进行中期预报。
(5) 由2011年独立样本检验和2012年预报检验的结果,可知香山红叶变色日气象统计模型具有比较高的准确率和可操作性。
需要指出的是,用直接观测的香山红叶变色率与气象条件相结合进行研究在国内尚属首次,且物候资料年份跨度较短,还有一些年份的中断,显然其结果是初步的。随着观测样本增多,一些新的特征可能会被揭示出来,对本文的结论也会有很好的补充。同时,香山红叶最佳观赏期与气象条件的关系也是值得进一步研究的问题。
陈正洪, 肖玫, 陈璇, 2008. 樱花花期变化特征及其与冬季气温变化的关系[J]. 生态学报, 28(11): 5209-5217. DOI:10.3321/j.issn:1000-0933.2008.11.002 |
范可, 林美静, 高煜中, 2008. 用年际增量方法预测华北汛期降水[J]. 中国科学地球科学, 38(11): 1452-1459. |
范可, 王会军, ChoiYoung-Jean, 2007. 一个长江中下游夏季降水的物理统计预测模型[J]. 科学通报, 52(24): 2900-2905. DOI:10.3321/j.issn:0023-074x.2007.24.014 |
葛雨萱, 王亮生, 周肖红, 等, 2011. 香山黄栌叶色和色素组成的相互关系及时空变化[J]. 林业科学, 47(4): 38-42. DOI:10.11707/j.1001-7488.20110407 |
刘建栋, 于强, 2002. 温度对夏玉米光合生产力影响的数值模拟研究[J]. 应用气象学报, 13(4): 397-405. |
刘杰, 许小峰, 罗慧, 2010. 气象条件影响我国农业经济产出的计量经济分析[J]. 气象, 36(10): 46-51. DOI:10.7519/j.issn.1000-0526.2010.10.007 |
刘玲, 郭建平, 高素华, 2006. 低温、干旱并发对玉米影响的评估研究[J]. 气象, 32(4): 116-120. DOI:10.7519/j.issn.1000-0526.2006.04.021 |
杨萍, 刘伟东, 仲跻芹, 等, 2011. 北京地区自动气象站气温观测资料的质量评估[J]. 应用气象学报, 22(6): 706-715. DOI:10.11898/1001-7313.20110608 |
姚秀萍, 吕明辉, 范晓青, 等, 2011. 气象服务效益评价研究进展[J]. 气象, 37(6): 749-755. DOI:10.7519/j.issn.1000-0526.2011.06.014 |
赵颖彦, 2011. 今日寒露催红京城[J]. 法制晚报, 2011-10-8(A03). |
张德山, 海玉龙, 冯涛, 等, 2010. 北京地区1~4天花粉浓度预报的应用研究[J]. 气象, 36(5): 128-132. DOI:10.7519/j.issn.1000-0526.2010.05.020 |
周肖红, 葛雨萱, 王亮生, 等, 2009. 黄栌叶片变色期生理变化及植物生长调节剂对叶色的影响[J]. 林业科学, 45(7): 59-62. DOI:10.11707/j.1001-7488.20090710 |
竺可桢, 宛敏渭, 1973. 物候学[M]. 北京: 科学出版社, 51-58.
|