大气污染对水稻产量的影响研究(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或最小二乘法或经验正交多项式来进行确定。