综合考察光、温、水等农业气候资源和土壤资源同作物生长发育和产量形成所必需条件的匹配性，是作物适宜性评价及作物适宜种植区划的主要手段，也是农业气象学的重要研究内容。小麦作为中国三大主要粮食之一，国内外学者对冬小麦和春小麦适宜性评价做了大量研究。如Tuan等^[1]，Mendas等^[2]，吴洪颜等^[3]，千怀遂等^[4]，蒲金涌等^[5]，Wang Dacheng等^[6]，王琦等^[7]对冬小麦做了适宜性的研究，李自珍等^[8]针对春小麦生态位适宜度模型进行了研究。近年来，冬小麦适宜性研究由定性向客观定量化方向发展，特别在指标的选择上，由传统的定性方法转向集成多种地理因子的综合环境指标定量数理统计方法。如钱存鸣等^[9]从江苏不同生态区选择7个地点，根据它们的气候条件、土壤类型和品质特点，将全省分成3个品质区域。王龙俊等^[10]根据生态环境因子、土壤、小麦消费习惯等将江苏省分为4个一级种植区，12个二级种植区，并提出框架性区划方案。以上研究应用少数的气候资料观测站点，并不能精确反映整个江苏的小麦气候生态状况。一些学者将GIS技术应用到农作物精细化评价中，如代立芹等^[11]利用模糊数据方法，建立河北省冬小麦各生育期气候适宜度评价模型，并分析了30 a中各生育期气温、降水、日照和气候适宜度的时空变化特征。马晓群等^[12]对安徽省冬小麦品种做了生态气候适宜性分析和精细化区划。这些研究定量到了空间面上，但主要考虑了气候对冬小麦生长的影响，以气候适宜性评价为主。也有少数学者综合考虑了气候、土壤等多种资源参与评价，如李奇峰等^[13]利用GIS技术构建评价单元，选取气候、土壤、地形等关键的作物生态环境因子，采用空间聚类分析法构建作物生态适宜性评价体系。陆洲等^[14]采用生态距离和二次阈值分级方法解决作物生态适宜性评价中分级缺乏物理意义的问题，并以北京的冬小麦为例做了验证。上述研究为深入开展冬小麦适宜性评价提供了较好的研究思路。

江苏省是中国重要的冬小麦产区之一，而影响冬小麦生长的自然资源要素不仅包括气候因素、还有土壤以及地形等因素，因此正确认识江苏省冬小麦种植的自然资源状况就需要尽可能全面地评价江苏省自然资源各因子与冬小麦生长需求的匹配程度。目前从模糊数学角度全面考虑气候、土壤中多个生态因子进行江苏省冬小麦种植适宜性评价的研究较少，从冬小麦具体品种层次上的研究报道也较少。本文根据模糊数学的方法和生态适宜度理论，从冬小麦品种扬麦16号的各生育期的气候指标出发，综合考虑农业气候资源和土壤资源对江苏省冬小麦品种的综合适宜性量化评价，实现区域空间上冬小麦品种的适宜性区划，以期为充分利用江苏省生态资源、指导及科学制定冬小麦品种区域种植规划提供依据。

1 研究数据及其处理 1.1 研究数据

由国家气候中心提供的江苏省及周边40个气象站1980—2010年日观测资料，包括平均气温、降水量、日照时数数据；江苏省土壤资料如土层厚度、土壤pH值和土壤有机质含量，由南京土壤所提供；从STRM网站下载90 m空间分辨率的数字高程模型DEM(digital elevation model)；江苏省冬小麦扬麦16号品种的物候数据及产量数据，由江苏省农科院提供。

1.2 数据处理

表 1给出了冬小麦品种扬麦16的7个生育期的划分标准。首先，将江苏省及周边40个气象站30 a的气温、降水、日照时数等气象数据按照表 1的生育期整理，并计算各生育期30 a气候平均值。其次，为了保证站点气象资料空间化的精度，对7个生育期分别采用反距离加权平均法(invers distance weight，IDW)、普通克里金(Ordinary Kriging，OK)、协同克里金(Co-Kriging，CK)以及径向基函数(radial basis function，RBF)中的规则样条、张力样条和薄板样条6种空间插值方法实现空间化，并对结果进行交叉检验，使用平均绝对误差(mean absolute error)和均方根误差(root mean square error)作为评价插值精度的依据。通过7个生育期6种插值方法的结果对比可得，气温最佳的插值方法为CK，降水最佳的插值方法为OK，日照时数插值最佳的方法为RBF的规则样条插值。

2 研究方法 2.1 适宜性评价模型

采用多准则综合评价模型计算每个评价单元的综合适宜度，作物适宜性评价模型^[15]为：

式中, IFI为综合适宜度；F_i为第i个评价指标的隶属度；W_i为第i个评价指标的权重。

参考文献^[16]，将江苏省冬小麦的种植适宜性划分为最适宜(＞0.75)、适宜(≥0.5—0.75)、次适宜(≥0.25—0.5)和不适宜(＜0.25)4个等级。

2.2 权重W_i设置

作物的生态适应性指数是由光照、气温、降水、土壤类型、有机质、pH值、氮磷钾含量等因素共同决定的。对于某种生物，各因子的重要性不同。根据最小因素定律，生物的生长发育是由数量最缺的因子决定的^[17]。考虑到不同因子在生态适应评价中决定强度的差异，本文采用考虑限制因子的相对权重法^[18]确定各评价因子的权重。具体权重公式(2)如下：

式中, W_i，含义同上；U_i为评价单元一个因子的隶属度；U_m为评价单元中隶属度最小的因子的隶属度。该公式表明在评价单元格中寻找隶属度最小的因子，若是小于0.2明显的限制因子，则赋予较大的相对权重，否则由其他因子的限制性共同决定的。

2.3 模糊数学隶属函数

作物对某环境因子的适应有很适宜、适宜或不适宜之分，这些都为定性描述，而利用模糊数学隶属函数的方法可以将各因子对作物的适宜程度定量化，此隶属度取值为0—1之间，隶属度值越大表明作物在该区域的生态适应性越高。依据生态位适宜度理论模型，作物对生态资源环境的要求主要有3类型：适宜区间型(如气温，pH值)、越多越好型(如土层厚度)、越低越好型(如土壤污染程度)。

2.3.1 气温隶属函数

气温存在适宜区间，作物光合作用有最低、最适和最高气温。对同一种作物来说，辐射强度的变化，会导致其同化率的变化，从而引起光合速率的气温曲线的变化。同样，不同环境气温下光合速率反应曲线也是不同的。在正常的环境气温下，净光合速率的最适气温不是一点，而是一个较宽的气温区域。作物通过光合作用，完成作物生长发育和产量形成的全过程。气温对光合作用强度的影响有两种效应：一是气温增高时光化学过程加快而使总光合作用强度加强，二是气温增高时呼吸消耗增加。因此净光合产物在初期随气温增加而加强，而在超过最适气温以后，净光合产物则随气温增加而降低。以往学者建立的线形等模型并不能充分表达气温对作物的影响变化特点，根据Larcher和Mirza等气温与作物光合作用的关系曲线，气温与光合速率反应曲线近似呈S型^[19-20]，结合研究地区实际情况，据此建立如式(3)的气温隶属函数S_T：

式中，T 是某一生育期内的实际日平均气温(℃)；T_L、T_H、T₂、T₁分别为作物在该时段内生长发育的下限气温、上限气温和最适气温的上下限(℃)。此处参数依据参考资料^[21-22]与江苏省农业科学院提供品种小麦栽培种植资料获取。

2.3.2 日照隶属函数

作物利用太阳辐射及其自身具有的光合作用功能，将太阳能转化为有机干物质。本文根据Camilo L.Medina^[23]、黄璜^[24]的研究，把日照时数达到可照时数的70%，即日照百分率70%以上的光照条件作为冬小麦的适宜状态。日照时数的隶属函数S_SH为^[25]：

式中，SH为实际日照时数(h)；SH₀为日照百分率为70%的日照时数(h)；b为系数。

2.3.3 降水隶属函数

评价降水对作物生长的影响，按照作物产量与降水量关系曲线图^[26]，建立作物的降水隶属函数S_R，计算公式为：

式中, R为某生育期内的平均降水量(mm)；R_L、R_H、R₂、R₁分别为作物在该时段内生长发育的下限、上限和最适需水量的上下限(mm)。k=1，n=2。此处参数依据参考资料^[21-22]与江苏省农业科学院提供品种小麦栽培种植资料获取。

2.3.4 土壤pH值隶属函数

建立土壤pH值的隶属度函数S_pH，如式(6)：

式中，pH为实际土壤pH值；pH_SL和pH_SH分别为适宜于作物生长的下限值和上限值；pH_OL和pH_OH分别为作物最佳生长的下限值和上限值，此处参数依据参考资料^[21-22]，分别取6.8和7.0。

2.3.5 土壤有机质隶属函数

土壤有机质隶属度函数S_OM^[27]：

式中, OM为实际土壤有机质含量(%)；OM_SL和OM_OL 分别为作物生长的土壤有机质含量适宜下限值及最佳下限值(%)，最佳下限值取1.3%。由于土壤有机质因子在最不利情况下，作物仍有可能生存，因此这里将生态位适宜度模型中的最低适应值由0换为 a，a∈[0，1)。

2.3.6 土层厚度隶属函数

土层厚度隶属度函数S_SD^[28]：

式中, SD 为实际土层厚度(cm)；SD_OL作物生长的土层厚度下限，取20cm。

3 结果与分析 3.1 各要素生态适宜度评价

根据以上气温、降水量和日照时数等各隶属函数及生态指标，可计算得到不同生育期各因子的适宜度。由气候因子各生育期的适宜度与相应生育期的权重乘积累和计算得到气候因子整个生育期的适宜度。冬小麦各生育期内与生育期间气象要素的权重结合文献^[11]和江苏省农科院品种小麦专家意见，得出整个生育期上各生育期的相对权重(表 2)，通过空间分析的栅格计算得到整个生育期的各气象因子适宜度(图 1)。因为土壤各因子随生育期的变化很小，可忽略，即只算整个生育期的土壤各因子的适宜度。

气温适宜度的范围为0.54—0.81(图 1)，由北往南逐渐降低，但在最南端如昆山、苏州等地适宜度较高；同一纬度上，东部地区的适宜度略高于西部地区；最适宜区域在江苏省北部，适宜度较低的区域在江苏省南部。由于播种期昆山、吴江地区适宜度为1，与周边适宜度相差较大；拔节期南京、溧水、溧阳、宜兴和吴江一带气温适宜度为1，与周边相差较大，使得将各生育期的气温取平均值

后，南京与吴江地区显得异常。

降水适宜度的范围为0.25—0.63(图 1)，基本上呈现由南部向北部逐渐降低；高邮、东台区域略低于其南北两侧；最适宜区域在江苏省南部，如无锡、苏州等地；适宜度较低的区域在江苏省北部，如徐州、连云港。

日照时数适宜度的范围为0.53—0.91(图 1)，江苏省大部分地区的日照时数都适宜冬小麦的生长。由北往南逐渐降低；同一纬度上，大部分西部地区的适宜度略低于东部地区；最适宜区域在江苏省北部，如连云港等；适宜度较低的区域在江苏省南部，如苏州等地。

经资料查阅，冬小麦的耕作层深度一般为20cm以上，土壤厚度的最低下限为40cm，计算的土壤厚度适宜度范围较大，见图 1。除去河流、湖泊，全省土壤厚度的适宜度大部分为1，即厚度大于40cm；江苏省的土壤厚度均适宜冬小麦的生长。

江苏省土壤pH值的适宜度的范围见图 1，适宜度大于0.6的区域占全省的61%。北部、苏中的西部和东部、苏南部分地区pH值适宜度较高，如徐州、宿迁、连云港、苏州等地；中部与苏南东部的部分地区适宜度为0，土壤酸性超出冬小麦适宜范围。全省有大约50%的面积pH值适宜度大于0.8。

土壤有机质含量的值达到1.3%为适宜^[29]，根据有机质隶属函数，计算的江苏省有机质适宜度的范围见图 1，适宜度大于0.8的区域占全省64%，大于0.6的区域占全省75%。大部分地区有机质含量较高，中、东南部如苏州、南通、泰州等地有机质适宜度最高；北部地区适宜度低于南部，但基本大于0.2，属于适宜范围。所以，除水域外，江苏省大部分地区的土壤有机质适宜冬小麦生长。

3.2 气候生态适宜度评价和土壤生态适宜度评价 3.2.1 冬小麦气候生态适宜度

单从气候的角度考查江苏省冬小麦品种的适宜性，可用气候适宜度来定量评价。利用式(1)，首先由各生育期的气温、降水、日照时数适宜度和冬小麦生育期内权重计算得到各生育期气候适宜度，再将不同生育期气候适宜度和相应的生育期间权重乘积求和，即得到冬小麦整个生育期的气候适宜度范围为0.68—0.81(图 2)，江苏省中北部大部分区域适宜度大于0.73，全省的气候适宜度均适合种植冬小麦，大体上呈现由北往南逐渐减小的变化。镇江、南通一带适宜度略低于其南部与北部，但都大于0.60。

3.2.2 冬小麦土壤生态适宜度

单从土壤的角度考查江苏省冬小麦品种的适宜性，可用土壤适宜度来定量评价。利用式(1)，综合土层厚度、土壤pH值、土壤有机质单因子的适宜度，其中3因子权重赋值依据文献^[30]，分别为0.29、0.19和0.52计算得到江苏省冬小麦土壤适宜度分布图 3。由图 3看出江苏省土壤适宜度为0—1，中部、南部与西北部等大部分区域土壤适宜度大于0.70，适合种植冬小麦，其中0.7以上区域占全省的82%，适宜度大于0.89的区域占全省的35.7%，主要分布在长江、淮河沿岸、太仓等区域。空间分布图上呈现南部优于北部的趋势。

3.3 冬小麦综合适宜度评价

冬小麦的生长是气候、土壤等多个环境因素的综合作用，因此全面考虑气候、土壤等因素的影响更为合理。根据以上气温、降水、光照、土壤类型、土壤有机质、土壤pH值各因子适宜度，通过权重公式(2)可得各因子的权重图层，利用式(1)计算得到冬小麦综合适宜度空间分布图 4。

从图 4中颜色深浅变化表示种植冬小麦适宜程度的高低。可以看出在各气候因子综合作用下江苏省冬小麦种植适宜度呈现北高南低。依据文献^[16]将适宜性分为最适宜、适宜、次适宜和不适宜4个等级，同时参考研究区土地利用图和冬小麦的产量等调查资料，做出分级调整，由冬小麦产量数据将最适宜区为适宜度0.7以上，适宜度0.65为适宜区和次适宜区的分界点，江河湖泊等水域都低于0.45，将适宜度低于0.45的划为不适宜区。由图 4可看出江苏省大部分区域适宜度为0.50—0.86，适宜种植冬小麦。最适宜的区域主要在江苏北部和中部区域，南部区域主要为适宜区，不适宜的区域主要分布在苏南的太湖、滆湖，苏中的高邮湖以及苏北的洪泽湖、长江等水域。该分级结果和吴洪颜等^[3]、金善宝^[17]学者研究结论相一致。

4 结论与讨论

本文利用多种插值方法对比选取最佳方法实现气象要素的空间化，利用模糊数学建立适宜度模型，分析江苏省冬小麦的气候、土壤适宜度情况，并将各要素的适宜度赋予权重，进行综合分析，得到江苏省冬小麦适宜种植区划。结果表明：

1)江苏省气温适宜度为0.54—0.81，日照适宜度为0.53—0.91，降水适宜度为0.25—0.63，总体上气候条件较适宜冬小麦生长。气温适宜度总体上呈从北往南逐渐降低的趋势；降水适宜度总体上呈现从南往北逐渐降低的趋势；日照时数适宜度从北往南逐渐降低。

2)江苏省大部分区域的土壤厚度能满足冬小麦的生长，适宜度为1；土壤有机质基本满足冬小麦生长，适宜度大于0.8的区域占全省64%，大于0.6的区域占全省75%；土壤pH值适宜度较适宜冬小麦生长，适宜度大于0.6的区域占全省的61%。

3)江苏省冬小麦整个生育期的气候适宜度范围为0.68—0.81，中北部大部分区域适宜度大于0.73，全省的气候适宜度均适合种植冬小麦，大体上呈现由北往南逐渐减小的变化。江苏省土壤适宜度为0—1，中部、南部与西北部等大部分区域土壤适宜度大于0.7，适合种植冬小麦，其中0.7以上区域占全省的82%，适宜度大于0.89的区域占全省的35.7%，主要分布在长江、淮河沿岸、太仓等区域。空间分布图上呈现南部优于北部的趋势。

4)江苏省冬小麦综合适宜度看，整个江苏省基本都适宜种植冬小麦，空间上看北部与中部最适宜冬小麦生长，南部较适宜冬小麦生长，不适宜区域主要分布在江河湖泊等水体区域。

本文研究中以冬小麦优质品种扬麦16号为例，该品种属于春性中熟小麦品种，在各生态要素的隶属度计算时，由于缺少该品种准确的三基点资料，根据专家意见，各气候要素的三基点参考中熟小麦的三基点资料，土壤pH值、土层厚度以及土壤有机质含量等三基点的取值参考冬小麦大类的三基点资料。这可能与该品种的实际三基点取值有偏差。

致谢: 江苏省农业科学院提供冬小麦资料, 范晓梅、孙德勇老师帮助写作，特此致谢。