南京站近60年农业界限温度变化特征分析

南京站近60年农业界限温度变化特征分析（南京信息工程大学应用气象学院，南京210044）摘要：用南京站1951—2011年三0°C、三10°C界限温度和初终日的观测资料，采用5d滑动平均法和三日连续偏低法对比分析日平均气温稳定通过各界限温度（三0C、三10C）的初、终日和初终间持续日数的差异，以及积温的年际和年代的变化。结果发现：各级界限温度初日均呈提前趋势，终日均呈推迟趋势，初终间持续日数呈增加趋势，活动积温和有效积温均表现为增大趋势。20世纪80年代，三0°C、三10C界限温度的积温偏少，90年代后逐年代增加的趋势，近几年热量资源最充沛，比80年代分别增加了501.1C,534.8C。关键词：界限温度初终日；积温；变化引言农业界限温度是标志某些重要物候现象或农事活动的开始、终止或转折点的日平均温度,对确定地区布局、耕作制度、结构调整等都具有十分重要的意义［1-2］。南京位于31.4—32.6N,118.22—119.4E,气候属于亚热带湿润气候，四季分明，光能充足，年平均气温16C，春季多晴天，夏季高温闷热，秋季干燥凉爽，冬季寒冷、干燥。南京春秋短、冬夏长，冬夏温差显著,在全球气候变暖的背景下［3-6］,相关界限温度变化及其对农业生产影响方面研究论文较多"-10］。本文采用南京站1951—2011年20°C、三10C界限温度和初终日的观测资料，讨论界限温度的变化特征，以及对农业气候资源、农作物种植结构和生产布局的影响。活动积温，是指某一时段内逐日平均气温累积之和，它是研究作物生长发育对热量的要求，同时也是评价热量资源的一种指标。10C是春季喜温作物开始播种与生长，喜凉作物开始迅速生长的界限温度。随着全球气候变化，各地作物关键生长季的界限温度和相应的热量资源已发生显著变化，李元华U1通过研究河北省近50a0C界限温度积温变化，得出界限温度初日提前，终日推迟，且积温增加的结论。1.资料及处理方法资料本文采用的基础资料是南京站（58238）1951.1.1—2005.1231逐日平均温度观测记录，为了保证数据的完整性，对于1951.1.1缺测的数据采用多年求平均的方法进行估计填充。资料处理方法界限温度初终日的确定采取5d滑动平均法和三日连续偏低法计算稳定通过各界限温度（20°C、三10°C）的初终日［12］。（1）5d滑动平均法：由春（秋）季开始出现高于（低于）某界限温度之日起，在依次算出的每个连续5日平均气温中，挑选出第1个其后不再出现低于（高于）该界限温度的五天，并在这五天中挑选出第1个（最后1个）三该界限温度的日期，作为稳定通过该界限温度的初（终）日；（2）三日连续偏低法：在春（秋）季，以最后（先）一次出现的连续3天或以上各天日平均气温都低于某界限温度的后（前）1天，作为稳定通过该界限温度的初（终）日。多年平均初终日及界限温度持续日数的确定（1）多年平均初终日［12］:①定基本月（n天），计算各年累积天数，从历年初日中选择最早出现的月份作为基本月，以基本月的第一天作为起始日，计算各年距起始日的天数；②求多年平均累积天数N；③确定多年平均初日，当多年平均累积天数N小于基本月的天数（n）时，基本月的第N天就是多年平均初日；当N大于n时，基本月之后第1个月的（N-n）日就是多年平均初日；④采用同样方法求多年平均终日。（2）界限温度持续日数：初日日期与终日日期间隔的天数。积温的求算方法活动积温：（；当时，以0计） （1）有效积温：（；当时，以0计）（2）式中：一一活动积温和有效积温（°c）n——植物某发育期的天数；一一发育期中第i天的日平均气温；B 界限温度（针对植物生长发育而言，B为生物学零度）。数据分析在Excel中求算界限温度的初、终日、初终间持续日数、活动积温和有效积温，在Matlab中进行图形和回归方程模拟,并使用5d滑动平均法和三日连续偏低法进行对比比较初日、终日、初终间日数、积温的年际以及年代变化特征。2.两种方法对比南京站界限温度变化特征2.1界限温度及其生物学意义适合南京气候特点和实际生产的界限温度及其生物学意义如下：0C表明春季地面冰雪开始融化，土壤开始解冻，冬小麦开始返青生长；秋季稳定通过0C,表明土地开始冻结，植物枯萎或休眠越冬，大田农事活动基本结束；20C期间为农耕期，期间的活动积温是衡量地区热量资源状况的重要指标。春季日平均气温稳定通过10C某些多年作物开始迅速积累有机物，杏花盛开，桃花怒放，小麦拔节，油菜花开,喜温作物可以开始播种［12］秋季日平均气温稳定通过10C喜温作物开始停止生长，水稻为喜温作物，日平均气温三10C的积温可以反映水稻生长的热量状况。$0°C、$10°C初日、终日、初终间持续日数、积温的年际变化特征2.2.1界限温度初日变化南京站日平均气温稳定通过各界限温度初日的变化趋势如图1,2所示。采用5d滑动平均法和三日连续偏低法模拟的三0C初日的回归方程分别为y=-0.2781x+47.6094,y=-0.4150x+37.9936采用5d滑动平均法和三日连续偏低法模拟的三10C初日的回归方程分别为y=-0.1874x+92.6178,y=-0.248x+92.8955。对比分析可得，南京站三0C、三10C初日均呈现出提前趋势，且表现出界限温度越高，初日提前更加显著的规律性（图1,2）。其中，三10C初日是水稻播种的时期，而三10C终日则是喜温作物停止生长的日期，水稻为喜温作物，故初日大幅度提前有利于农作物产量的增加。：wc）图1南京站界限温度（三0C）初日年际变化图2南京站界限温度（三10°C）初日年际变化2.2.2界限温度终日变化南京站日平均气温稳定通过各界限温度终日的变化趋势如图3,4所示。采用5d滑动平均法和三日连续偏低法模拟的三0C终日的回归方程分别为y=0.0492x+360.481,y=0.1444x+354.3081；采用5d滑动平均法和三日连续偏低法模拟的三10C终日的回归方程分别为y=0.0728x+323.9686,y=0.0775x+314.2691。对比分析可得，南京站三0°C、三10°C终日均呈现出推迟趋势，比较而言，三10°C终日推迟趋势更为明显。350M340360-a32019601970200020102020I9601970加20102020370AA330360AAAAA31030020002010202019601970200020102020AAA«310AA%乱 AA师 1990轴19® 19901990 1990年份1980 1990轴卫350M340360-a32019601970200020102020I9601970加20102020370AA330360AAAAA31030020002010202019601970200020102020AAA«310AA%乱 AA师 1990轴19® 19901990 1990年份1980 1990轴卫A『沁M2A.A A.AA2A尸0«92"®Q1子册糊聪A•三日蹦魅r&-力捕时聪］y=0Mx+3142691SAAAAAAAAAMAAAAAAAAAAAAAAAAAAMAAAAAAAA图3南京站界限温度（三O°C）终日年际变化图4南京站界限温度（三1O°C）终日年际变化2.2.3界限温度初终间持续日数变化日平均气温大于某界限温度的初、终日之间的持续日数如图5,6所示。界限温度初终间持续日数能够较好地表示相应农事季节或可供某类作物生长发育时间的长短，是确定某种作物类型和调整种植结构的重要依据。采用5d滑动平均法和三日连续偏低法模拟的三0C初终间持续日数的回归方程分别为 y=0.3273x+312.6005,y=0.5594x+316.3082；采用5d滑动平均法和三日连续偏低法模拟的三1O°C初终间持续日数的回归方程分别为y=0.2602x+231.350,y=0.3255x+221.3736对比分析可得，各界限温度期间的持续日数均呈增加趋势，表明伴随温度升高，界限温度初终日之间的天数也随之增加。比较而言，2O°C的持续日数更为明显。肌）ISO M I960 1990 2000 2010 2020的持续日数更为明显。肌）ISO M I960 1990 2000 2010 2020图5南京站界限温度（三OC）初终间持续日数年际变化图6南京站界限温度（三1OC）初终间持续日数年际变化界限温度积温变化南京站活动积温、有效积温的年际变化如图7,8所示。三O°C、三1O°C的活动积温回归方程分别为y=80x+55670,y=80x+47130（注：y为活动积温，但是包含了一位小数）；2O°C、三1O°C的有效积温回归方程分别为y=80x+55670,y=70x+2900（注:y为活动积温，但是包含了一位小数）。分析可得：（1）南京站界限温度积温随年际均表现为增加趋势，且有效积温和活动积温基本一致。（2）温度升高对积温的影响，十分明显。这表明，积温增多的变化趋势有利于种植业的稳定，也为种植业的结构调整提供了热量资源。

1®2H01S60KJ20102020训h19501960ffl200020102020瓯m2010$阴 m俶 1990ISO1®2H01S60KJ20102020训h19501960ffl200020102020瓯m2010$阴 m俶 1990ISOinm冋”26300图7南京站界限温度（三O°C、三1O°C）活动积温年际变化图8南京站界限温度（三O°C、三1O°C）有效积温年际变化$O°C、$1O°C初日、终日、初终间持续日数、积温的年代际变化特征（一）南京站2O°C界限温度的年代际变化来看（表1）。具有以下特点：1） 2O°C的初日出现明显的年代际变化，随着年代变化有所提前（表1）；三OC平均初始日期为2月8日（5d滑动平均法）或1月25日（三日连续偏低法）（表2）；2） 2O°C的终日，随着年代变化有所推迟，但变化幅度不大（表1）；最晚终止日期为12月31日（5d滑动平均法和三日连续偏低法）（表2）；3） 2O°C的初终间日数，随着年代变化增加，近几年比50年代平均提前了22d，近几年比50年代增加了15—29d（表1）；4） 各级界限温度期间活动积温和有效积温的年代际变化具有一致的规律性：即，NTC活动积温和有效积温，在80年代后呈现逐年代增加的趋势，近几年热量资源最充沛，比80年代增加501.1C（表1）表1日平均气温稳定通过OC的初日、终日、初终间日数、积温变化年代初日（日/月）终日（日/月）持续日数（天）5日滑动平均3日连续偏低活动积温（C）5日滑动平均法3日连续偏低法5日滑动平均法3日连续偏低法法法50年代13-210-228-1221-123183145596.060年代10-231-127-1223-123203265692.370年代1-227-129-1224-123313315626.280年代14-223-125-1222-123143335562.190年代3-216-129-1227-123293455816.7近几年1-220-131-1229-123333436063.2

表2南京站稳定通过界限温度0°C日期初日终日持续日数

日/月） （日/月） （天）项目5日滑动平均法3日连续偏低法5日滑动平均法3日连续偏低法5日滑动平均法3日连续偏低法平均8-225-127-1224-12322333最早1-11-127-1127-11330330最晚22-315-331-1231-12284277（二）南京站三io°c界限温度的年代际变化来看（表3）。表现为初日提前、终日推迟、持续日数延长、积温增加的趋势。其中，①三10C初始日期变化不大，三10C初始日期在60年代最晚，为4月2日；90年代和2000年以后比50年代提前了2d和9d；但60年代三10C初始日期比50年代延迟了Id；三10C平均初始日期为3月28日（5d滑动平均法）或3月26日（三日连续偏低法）（表4）；②与210C初始日期提前程度相比较，终止日期延迟程度小于初始日期提前程度，且终止日期延迟程度趋势不明显（表3）。最晚终止日期为12月14日（5d滑动平均法）或12月1日（三日连续偏低法）（表4）；210C平均终止日期在11月23日（5d滑动平均法）或11月14日（三日连续偏低法）（表4）。③210C平均持续日数为240d（5d滑动平均法）或233d（三日连续偏低法）（表4）,并且，在此期间大多数年份持续日数以增加为主（表3）；④210C积温，50年代到80年代略有下降，80年代后随着年代的变化积温均在增加，近几年增加最为明显，比80年代增加了534.8C（表3）。表3日平均气温稳定通过10C的初日、终日、初终间日数、积温变化年代初日（日/月）终日（日/月）持续日数（天）5日滑动平均3日连续偏低活动积温（c）5日滑动平均法3日连续偏低法5日滑动平均法3日连续偏低法法法50年代1-431-326-1114-112392285280.060年代2-41-422-1111-112342245170.070年代3-328-317-119-112322265050.580年代26-327-323-1112-112422305002.690年代30-328-324-1117112392345245.7近几年21-317-328-1117-112522455537.4表4南京站稳定通过界限温度10C日期初日 终日 持续日数日/月） （日/月） （天）项目5日滑动平均法3日连续偏低法5日滑动平均法3日连续偏低法5日滑动平均法3日连续偏低法平均28-326-323-1114-11240233最早24-224-27-1125-10256243最晚20-415-414-121-122382303结论1）南京站稳定通过三0C初日在2月25日，终日在12月25日；南京站稳定通过三10C初日在3月13日,终日在11月18日。20°C初日呈现出提前趋势，平均每10a提前2.4d；终日呈现出推迟趋势，平均每10a推迟0.6d；初终间