气候变化是当前热点问题,众多学者对区域气候变化进行过深入研究,如刘裕禄等[1]应用EOF进行降水时空分布研究,效果很好。但现有相关文献大多是针对单一气候要素气温或降水进行的[1-6],而关于光、温、水三者相结合的农业气候变化研究报道甚少。地处低纬高原的云南是南亚和东亚季风影响的过渡区,它立体气候突出,干湿季分明,冬夏半年日照差异大[7]。在全球气候变化背景下,云南气候变化特点鲜明。植物生态学研究和生产实践认为,气候对植物生长、产量形成和农业的影响包括光、温、水总量的匹配,气候要素的时段变化和气象灾害的综合作用[8-9],而非仅仅单一要素的影响。云南大部农业生产一年两熟(与干、雨季对应),秋收农业产量和产值占全年总产的80%以上[10],雨季5—10月还是森林等植物的主要生长期。分析研究云南农业主要生长季气候变化规律,对生产实践和气象服务意义重大。作者基于3个时段光、温、水匹配和农业气候主成分逐年变化的深入分析,研究云南农业主要生长季气候变化特征,揭示其主要农业气候变化规律和影响。
1 应用资料和农业气候因子的选取 1.1 资料选择数据完整性和连续性最好,分布于滇中、滇东北、滇东南、滇西北和滇西南的8个站点(昆明、昭通、沾益、蒙自、丽江、大理、思茅、腾冲),使用1954—2007年逐年月、旬和候平均气温、降水量和日照时数进行分析研究。
1.2 因子的选择云南主要农业气候影响时段为3—10月(秋收作物生长季)。其中3—5月为播种移栽和幼苗期,6—8月为作物营养生长中后期和生殖生长中期,7—8月为产量形成关键期,9—10月为经济产量形成期(和晚秋作物重要生长期)。因作物各生育阶段对气候要求不同[8-9]以及气候的季节性变化,不同阶段气候对作物的利弊影响差异大。根据相关研究[7],云南农业气候主要问题可归结为:3—5月春旱、短时低温冷害和部分地区3—4月热量较差;6—8月低纬高原气温偏低、日照偏少、雨量偏多同时部分区域受旱涝影响;9—10月常受阴雨寡照、热量强度不足和短时低温困扰。本研究以反映此3个时段光、温、水总量和上述农业气候主要问题为依据,在分析各时段光温水常规要素与阶段性农业气候问题相关性基础上,筛选17个因子来共同表征云南作物主要生长期的农业气候。选定的气候因子3—5月时段为:3—4月平均气温(T3-4)、3—4月最小候平均气温(Th3-4min )、4月至5月中旬雨量(R4-5)、4—5月日照时数(S4-5)、4—5月2旬滑动平均日照时数最小值(S4-5s2min);6—8月时段为:总雨量(R6-8)、3旬滑动平均雨量最小值(R6-8s3min)和最大值(R6-8s3max)、总日照时数(S6-8)、2旬滑动平均日照时数最小值(S6-8s2min)、7—8月平均气温(T7-8)、7—8月最小候平均气温(Th7-8min);9—10月时段为:平均气温(T9-10)、最小候平均气温(Th9-10min)、总雨量(R9-10)、日照时数(S9-10)、2旬滑动日照时数最小值(S9-10s2min)。
在初选气候因子中,最小候平均气温反映短时低温,连续2旬滑动平均日照时数表示寡照阴雨,连续3旬滑动平均雨量最小值和最大值分别指示干旱和洪涝。3月日照充足、5—6月气温较高,三者对农作物产量影响很小,不予考虑。云南雨季一般5月中下旬开始,5月下旬降雨普遍陡增量大,选用4月至5月中旬雨量更能反映春旱影响程度。云南夏秋季降水充沛、日照较少,2旬滑动日照和3旬滑动雨量最小值与同期干旱相关性极小;其余月段日照时数、平均气温和降水量均代表 3个特定时段或作物敏感时段的光、温、水总量。限于篇幅和研究主题,各初选气候因子的物理意义不予详述。
数据预处理先统计8站点逐年各因子值,再分别对各因子求8站点逐年平均值,最后应用17个变量、54年的样本序列,分析研究云南主要农业气候变化。
2 研究方法 2.1 应用EOF分析研究3个时段农业气候变化特征主成分分析(EOF)是从多原始变量的样本序列中归纳提练出潜在的、相互独立的、少数不同类别“公共因子”的一种统计学方法[11-12],每个“公共因子”都由少数原始变量组成并具有一定物理意义。EOF分析是研究变量场主要和次要变化规律的有效方法。分析流程[11-12]可简述为:
使用变量数为m、样本数为n的原始变量序列,构建Xn×m阶矩阵;求其标准变量场的方差协差阵Lxx=XTX及其所有特征根和特征向量矩阵;提取特征根>1的前k个主成分因子[11](即“公共因子”),特征根值越大,其主成分因子方差对原始变量场总方差的解释度越大;用前k个主成分的特征根和其特征向量构造新矩阵;根据Fn×k=Xn×mAm×k得到前k个主成分因子表达式,即使用原始变量可计算得到逐年各主成分因子值。根据变量场的主成分因子构成、各主成分特征根大小和主成分因子值的逐年变化,可分析提练出变量场的时间或空间变化特征。
2.2 应用系统聚类法研究农业气候匹配年型及其年际变化气候对农业的影响不是多个气候因子影响的简单叠加,而是相同时段多气候因子匹配的综合作用和多个时段不同气候匹配的共同影响[8-9],即不同的农业气候年型具有不同的农业效应。系统聚类法适用于农业气候的年型划分,是对多变量(或样本)按照相似程度(以距离表示)进行客观分类或归类的一种成熟统计方法[11-13]。在此以EOF分析得到的3时段所有主成分因子作为新变量组,经对历年新变量样本标准化处理后,针对多个主成分因子历年值的匹配进行系统聚类和年份归类。依据划分的年型种类和各分类的主成分因子多年平均值可确定主要农业气候匹配年型的标准参照系,而某个农业气候年型的物理意义可由其主成分因子物理意义和因子值来解释。按照农业气候匹配年型划分和逐年归类结果,可分析关注年份和各年代农业气候特点或评判其对农业的影响。
3 结果和分析 3.1 6—8月时段EOF分析经分析,EOF对6—8月时段5个变量可提取84%~95%的原始信息,对另外2变量(3旬滑动平均雨量最大值和最小值)可分别提取77.8%和68.2%的原始信息;提取的3个主成分因子(见表 1)其特征根对原始变量场总方差的贡献率分别为46.37%、22.33%和15.89%,三者累计贡献率为84.58%。
![]() |
表 1 6—8月各成分因子特征根及对原始变量场总方差的解释度 |
表 2为6—8月时段EOF分析得到的两个矩阵,主成分因子负荷矩阵表示3个主成分因子对各原始变量(标准化值)变异的解释度,而主成分因子得分系数矩阵反映各原始变量(标准化值)对3主成分的贡献大小。
![]() |
表 2 6—8月主成分因子特征向量负荷矩阵和因子得分系数矩阵 |
由于变量经标准化处理后消除了量纲的影响,故表 2中主成分因子得分系数绝对值越大表示相应变量对其主成分因子的影响越大。根据EOF分析和表 2,6—8月第1主成分因子可简化表达为:
$ \begin{array}{l} Z{F_{1 - 6 - 8}} \approx 0.44Z{R_{6 - 8}} + 0.381Z{R_{6 - 8s3{\rm{max}}}} + \\ 0.356Z{R_{6 - 8s3{\rm{min}}}} + 0.178\mathit{Z}{\mathit{T}_{h7 - 8{\rm{min}}}} \end{array} $ | (1) |
式中各变量代号前加Z表示为变量标准化处理值(以下相同),等式(1)右边项系数值大的前3项均为雨量项。可见,第1主成分因子主要反映3个雨量项的贡献,即雨量越大,第1主成分因子得分越高,洪涝越重;反之,夏旱影响越明显。同样,第2主成分因子表达式可列出(略),F2主要表示2个日照项的变化,日照项越大,F2因子得分越高,说明夏季光照条件越好。同理,第3主成分因子表示热量条件,6—8月平均气温和最低候平均温度越高,F3得分越高,代表夏季热量条件越好。
3.2 3—5月和9—10月两时段EOF分析EOF分析提取了3—5月5个变量原始信息的78.6%~90.2%,和9—10月5个变量原始信息的78.5%~86.6%。对3—5月提取了2个主成分因子,其特征根对总方差贡献率分别为49.01%、3 1.952%,两者累计贡献率80.96%。对9—10月提取的2个主成分因子,其特征根贡献率分别为53.51%和27.811%,累计贡献率为81.32%。
3—5月时段第1主成分因子表达式为:
$ \begin{array}{l} Z{F_{1 - 3 - 5}} \approx 0.385Z{S_{4 - 5s2{\rm{min}}}} + 0.379\mathit{Z}{\mathit{S}_{4 - 5}} - \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;0.341Z{R_{4 - 5}} \end{array} $ | (2) |
即3—5月F1综合反映日照和雨量的变化,其中日照项和雨量项分别为正、负效应;F1值大,表示3—5月光照好、阴雨寡照影响小;因春季属云南气候干旱最重季节,F1值大也意味着春季旱情偏重。第2主成分表达式为:
$ \begin{array}{l} Z{F_{2 - 3 - 5}} \approx 0.556Z{T_{3 - 4}} + \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;0.554Z{T_{h3 - 4{\rm{min}}}} \end{array} $ | (3) |
根据式(3),3—5月F2表征春季热量条件和低温影响,F2值大表示春季气温较高、低温影响小。
9—10月时段第1主成分因子表达式为:
$ \begin{array}{l} Z{F_{1 - 9 - 10}} \approx 0.435Z{S_{9 - 10}} + \\ 0.384Z{S_{9 - 10s2{\rm{min}}}} - 0.344Z{R_{9 - 10}} \end{array} $ | (4) |
显然,9—10月F1主要反映秋季阴雨寡照影响,两个日照项为正效应,雨量项为负效应,F1值大,说明9—10月光照条件好、阴雨寡照影响小。9—10月第2主成分表达式为:
$ \begin{array}{l} Z{F_{2 - 9 - 10}} \approx 0.551T{h_{9 - 10{\rm{min}}}} + \\ 0.549Z{T_{9 - 10}} - 0.141Z{S_{9 - 10}} \end{array} $ | (5) |
即9—10月F2主要反映秋季热量条件,虽然有日照项,但其负效应小;F2值大,表示9— 10月气温较高、低温影响不明显。
3.3 三个时段各主成分因子的逐年变化根据各主成分因子逐年值,绘制了其逐年变化图(图 1)。图中还给出了由正交多项式模拟的变化趋势线。
![]() |
图 1 云南3—10月农业气候主成分因子得分逐年变化 (a)3—5月第1主成分(阴雨寡照和春旱),(b)3—5月第2主成分(春季低温);(c)6—8月第1主成分(夏季雨量和涝旱),(d)6—8月第2主成分(光照);(e)6—8月第3主成分(气温),(f)9—10月第1主成分(阴雨寡照);(g)9—10月第2主成分(气温) |
根据EOF理论,云南农业气候7主成分因子是相互独立的。同时图 1说明:
(1) 图 1(b)、(e)、(g)反映3时段云南气温的逐年变化。三者都有增温变暖的趋势,夏、秋季变暖较显著;进入1990年代中期后各时段气温均进入快速上升期,这与全球气候变暖趋势相呼应;温度的低值区春季出现在1990年代中期、夏季为1970年代中期,而9—10月多出现于1950年代末期。
(2) 根据图 1(a),1960年代以来春旱强度有缓慢减弱的趋势,1960年代春旱最重,1990年代春旱次重;1999—2007年春旱强度最弱。图 1(c)显示夏季雨量有减少之势,1990年代后期以来递减趋势尤为显著。
(3) 根据图 1(d),1960年代以来夏季日照时数逐渐递减趋势明显,1960年代和1970年代前期为峰值期,1990年代为低谷期,近几年有所回升。图 1(c)和(d)比较,夏季日照时数年代际振幅大于雨量的。
(4) 图 1(f)说明,进入1970年代以来9—10月阴雨寡照逐渐加重,1960年代前期和1970年代影响最小,1980年代中期至1990年代初期影响最重,近年仍处于比常年偏重期。
3.4 系统聚类分析对7个主成分因子历年值进行农业气候匹配年型的系统聚类分析,结果见表 3。
![]() |
表 3 3—10月云南农业气候匹配年型划分及其主成分因子平均值 |
根据表 3,云南农作物主要生长季有6类农业气候匹配年型,各年型气候特点鲜明。比如第1类(Type1),春季雨量正常稍多旱情略轻、气温正常稍高,夏季雨量稍多旱涝正常、日照接近常年但热量条件特差(有夏季低温),9—10月光热条件好(气温偏高突出)。第5类,春季少雨多光照旱情重、低温影响小,夏季雨量稍少、日照特少、气温略偏高,9—10月阴雨寡照典型偏重、热量条件较好。第6类,春季雨量偏多日照特少干旱不明显、低温影响小,夏季雨量典型偏少、日照略偏差、热量条件好,9—10月阴雨寡照典型突出、热量条件优越。其余各类年型略述。
基于表 3,最近54年云南有10个年份主要生长季农业气候年型与其上年相似,但连续3年年型相似的只出现过1次(1999—2001)。1960年代以第4、3、2类为主,1970年代以第2、3类为主,1980年代以第4、5、6类为主,1995年后以第5、6类为主(概率为11/13);2001—2007年代主要为第6类。
4 讨论(1) 基于气候变化的农业和生物学意义,针对区域农业气候特点来划分研究时段和选定分析要素项,所选17个因子考虑了3时段气候三因素的总量,主要不利短期气候和中等时段要素值3个方面。
(2) 3个时段都分析了气候三因素总量,代表了整个农业生长季的光温水总量;又根据文献[14],云南各气候要素全生长季的相对变异比各月段的明显要小,故在因子选择上不必单独考虑全生长季气候要素总量。另一方面,不但分3个时段提取相互独立主成分因子,而且从主成分因子的不同匹配上进行了研究。因子选择和研究方法具有合理性和科学性。
(3) 进入1990年代中期后,第5、6类云南农业气候年型出现概率很高,根据植物生态学原理,此二类气候年型的农业效应差异大[15];进入21世纪后,第6类农业气候年型出现的概率达3/4,其气候特点不仅是温度高,而且春旱特轻,夏季雨量典型偏少、秋季阴雨寡照重。说明云南农业主要生长季气候变化加剧了,或许正进入一个与1950—1980年代时期偏差大的阶段。
(4) 3个时段提取的第1主成分反映的都不是温度,而是日照或雨量;春、秋季日照对第1主成分的影响大于雨量,即光照和降水的变化比气温变化更显著,很值得关注。
(5) 关于农业气候匹配年型的划分及其变化特点,具有明确的农业和生物学意义,对做好云南气候影响评价、产量预测和农业年景展望,具有重要指导意义。
5 结论基于前面的分析研究,可得到如下主要结论:
(1) 云南农业气候变化问题复杂,其主要生长季农业气候变化特征可用7个相互独立的主成分因子描述。夏季降水变化最显著,气温相对变化较小;春、秋季日照年际变化最明显。
(2) 云南农业主要生长季各时段都有气候变暖趋势,夏、秋季变暖特征比春季明显;进入1990年代中期后气候变暖加剧,同时日照变幅大、降水量夏季雨量递减显著。
(3) 1960年代以来云南春旱强度有缓慢减弱的趋势,1999—2007年为最近54年春旱最弱期。1960年代以来夏季日照递减趋势明显;1970年代以来9—10月阴雨寡照影响逐渐加重。
(4) 夏季雨量年际波动比日照显著,但日照时数年代际变幅大于雨量的;春、秋季日照年际变化都比雨量显著。
(5) 作物主要生长季云南农业气候有6种年型。只有少数年份的农业气候年型与其上一年相似;进入1990年代后云南农业气候变化明显加剧,主要出现两种农业气候年型;1999—2007年云南农业气候年型独特:各时段气温高、低温影响小;夏季雨量典型偏少;春季干旱影响小;秋季阴雨寡照突出。
刘裕禄, 许柏林, 许剑勇, 2007. 黄山地区汛期(5-9月)降水时空特征[J]. 气象, 33(8): 47-51. DOI:10.7519/j.issn.1000-0526.2007.08.008 |
刘燕, 程正泉, 叶萌, 2008. 广州市气温的气候变化特征及其成因分析[J]. 气象, 34(2): 52-60. DOI:10.7519/j.issn.1000-0526.2008.02.008 |
罗伯良, 张超, 林浩, 2008. 近40年湖南省极端强降水气候变化趋势与突变特征[J]. 气象, 34(1): 80-85. DOI:10.7519/j.issn.1000-0526.2008.01.012 |
姚愚, 李晓鹏, 闫丽萍, 等, 2006. 近44年云南年平均气温的时空变化特征[J]. 气象, 32(10): 81-87. |
尤卫红, 傅抱璞, 林振山, 1997. 云南近百年气温变化与8月低温冷害天气[J]. 高原气象, 16(1): 63-72. |
段玮, 肖子牛. 近33年云南降水及其构成的气候变化[C]. 中国气象学会2006年年会"首届研究生年会"分会场论文集: 32-36. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGQX200610011036.htm
|
王宇, 1990. 云南省农业气候资源及区划[M]. 北京: 气象出版社, 10-12.
|
韩湘玲, 曲曼丽, 1991. 作物生态学[M]. 北京: 气象出版社, 27-39, 42-46, 51-63, 67-70.
|
崔读昌, 1994. 世界农业气候与作物气候[M]. 杭州: 浙江科学技术出版社, 3-5, 8-9.
|
云南省统计局, 2007. 云南统计年鉴[M]. 北京: 中国统计出版社, 239.
|
张文彤, 2002. SPSS11统计分析教程[M]. 北京: 北京希望电子出版社, 171-177, 191-202.
|
严华生, 王学仁, 1991. 多因变量及要素场统计预报[M]. 北京: 气象出版社, 100-104.
|
高惠璇, 2005. 应用多元统计分析[M]. 北京: 北京大学出版社, 102-118.
|
黄中艳, 朱勇, 王树会, 等, 2007. 云南烤烟内在品质与气候的关系[J]. 资源科学, 29(2): 83-90. |
冯秀藻, 陶炳炎, 1991. 农业气象学原理[M]. 北京: 气象出版社, 86-92, 189-191, 193-199.
|