大气污染对水稻产量的影响研究(3)

杭州   双季晚稻单产   1999—2014   　　杭州统计信息网 http://www.hzstats.gov.cn/

南县   晚稻单产   1981—2014   　　湖南省统计局 http://www.hntj.gov.cn/
成都   单季稻单产   1981—2012   　　　　四川省统计局http://www.sc.stats.gov.cn/

重庆   单季稻单产   1981—2014   重庆省统计信息局http://www.cqtj.gov.cn/

武汉   一季晚稻单产   1981—2014   湖北省统计局http://www.stats-hb.gov.cn/

南昌   一季晚稻单产   1985—2009   　　　江西省统计局 http://www.jxstj.gov.cn/

合肥   单季稻单产   1981—2014   　　　安徽省统计局http://www.ahtjj.gov.cn/

牡丹江   单季稻单产   1981—2014   中国气象数据网 http://data.cma.cn

全部   日均温、月降水量   2000—2012   中国气象数据网 http://data.cma.cn

全部   日API值   2000—2012   中华人民共和国环境保护部http://www.mep.gov.cn/

1.2 数据处理
根据各站点水稻的抽穗期和收获期，划定各地的主要研究时段，即水稻关键生长期，其中南昌、益阳为8月-10月，武汉、牡丹江、合肥为7月-9月、成都为6月-8月，重庆为5月-7月，并通过计算得到各站点每年水稻关键生长期内的平均API、降水总量（mm）以及有效积温（℃•d）。平均API为日API的简单平均，降水量总和（P）为直接相加求值。
有效积温（GDD）计算方法采用“考虑无效高温的平均气温法”[17]，其中，水稻的生长下限温度为10℃，生长上限温度为40℃ [15-16]。
       N2
Tsum = ∑ △T i
     i= N1
         0         当Ta i≤Tl 或 Ta i≥Th
△T i =   Ta i - Tl       当Ta i＞Tl
         Th - Ta i       当Ta i ＜Th
Tsum为生育期积温，单位为℃•d。N1、N2为生育期开始与结束的日期，△TI为每日温度对作物生长发育的贡献值，Ta i为日平均气温，Th为作物的生长上限温度，Tl为作物的生长下限温度，超过Th的高温或低于Th的低温将视为作物生长发育的无效温度。
1.3 分析方法
1.3.1 API变化趋势
    气候倾向率的公式为:设某站点某气象要素时间序列为y1、y2、yi、yn，它可以用一个多项式来表示：
y(t)=a0+a1t+a2t2+amtn (m<n)
    式中t为时间，单位为a。一般来讲，API的气候趋势公式用一次直线方程和二次曲线方程就能满足。这里用一次直线方程来定量描述，即
y(t)=a0+a1t
则API趋势变化率方程为
dy(t)/dt=a1
    即a1×10称做API的气候倾向率，其单位为1/10a，方程中的系数可用Excel或最小二乘法或经验正交多项式来进行确定。大气污染对水稻产量的影响研究(3):http://www.youerw.com/shengwu/lunwen_21103.html