1. 引言
在大豆产量方面,当下我国总产量占据世界各国的前列。但若从国内供需关系方面看,存在着远远不能满足我国内部自主需求的情况。而每年我国进口量持续增加,已成为全球第一大豆进口国[1]。目前如何提高大豆平均亩产值已经成为了一项影响国家农业进出口贸易的问题,利用现代科学技术便是解决这一问题的一个很好举措。
黑龙江省作为我国的大豆种植大省,其产量更是占据了我国总产量的百分之四十,世界大豆总产量的百分之六。一直以来产量都远远高于其他种植大省,但近年来存在减产情况。近10年来黑龙江省大豆平均单产减产的主要原因是大豆重迎茬所致[2]。当下先一步的实现提高黑龙江省大豆亩产值的目标,对实现提高全国大豆亩产值有示范性的作用,对提高全国大豆产值及品质有重要意义[3]。
本文为进一步提高黑龙江省大豆产量,基于GIS技术对黑龙江省进行第一积温带大豆种植区域规划。使用MAPGIS软件系统工具,结合相关气象数据,对黑龙江省第一积温带大豆种植区域进行划分,希望这一举措能对黑龙江省大豆产量及品质的提高有促进作用[4]。
作为大豆适宜生长的黑龙江省第一积温带区域,在全球变暖的趋势下已出现不规则移动,若仍按照过往的积温带划分来对大豆种植区域进行规划已不适宜当下情况,不利于大豆产值及品质的提高,甚至导致产值及品质下降的情况发生。为避免这一情况的出现,对黑龙江第一积温带大豆种植进行区域规划便势在必行。地理信息系统(GIS)软件较之其他软件在气候区域规划方面有极大优势,使用地理信息系统(GIS)能有效收集相关气候资料及地理环境信息能对大豆种植进行科学指导,利于对黑龙江省第一积温带大豆种植区域的重新规划,能有效剔除黑龙江省当下不适宜的大豆种植的区域,能有效将气候影响减到能控制的最小范围,在减少无效种植投入及相关人工维护成本的同时,还能对未来的气候变迁情况进行预测,做到提早发现、提早防备、杜绝隐患,进一步提高大豆产量及品质,本文以小范围内对黑龙江省第一积温带大豆种植区域规划为例展开研究,以期实现提高我省提高大豆平均产值的目的。
2. 资料与方法
2.1. 研究区
研究区主要位于哈尔滨、齐齐哈尔、大庆、东宁、肇东等五个主要区域进行自然地理概况的介绍,并针对大豆种植的影响因素如日照强度、降水、活动积温等因素[5],对其历年的变化趋势进行汇总分析。区域内自然条件相似,耕作习惯、管理方式和水平接近[6]。黑龙江省第一积温带所选的五个区域的地理信息(表1)及气候因素(表2)。
Table 1. Geographical information of five regions in the first accumulated temperature zone of Heilongjiang Province
表1. 黑龙江省第一积温带五大区域地理信息
地区 |
地理位置 |
总面积 |
海拔 |
哈尔滨 |
黑龙江省中南部 |
53,100 km2 |
127 m |
大庆 |
黑龙江省西南部 |
21,219 km2 |
146 m |
齐齐哈尔 |
黑龙江省西南部松嫩平原 |
42,286 km2 |
149 m |
东宁 |
牡丹江市东南部 |
7139 km2 |
111 m |
肇东 |
松嫩平原中部 |
4332 km2 |
198 m |
Table 2. Climatic factors of five regions in the first accumulated temperature zone of Heilongjiang Province
表2. 黑龙江省第一积温带五大区域气候因素
地区 |
气候类型 |
春最高/低温 |
夏最高/低温 |
秋最高/低温 |
冬最高/低温 |
年日照时数 |
年平均降水 |
哈尔滨 |
温带季风 |
12℃/1℃ |
27℃/18℃ |
11℃/0℃ |
−10℃/−20℃ |
1490 h |
554 mm |
大庆 |
温带季风 |
12℃/1℃ |
27℃/18℃ |
10℃/1℃ |
−10℃/−19℃ |
1470 h |
522 mm |
齐齐哈尔 |
中温带季风 |
12℃/0℃ |
27℃/18℃ |
10℃/0℃ |
−10℃/−21℃ |
1419 h |
475 mm |
东宁 |
大陆性季风 |
13℃/0℃ |
26℃/16℃ |
13℃/1℃ |
−5℃/−16℃ |
1408 h |
611 mm |
肇东 |
寒温带 |
12℃/0℃ |
27℃/17℃ |
11℃/0℃ |
−10℃/−22℃ |
1467 h |
524 mm |
2.2. 数据来源
本研究所采用数据主要来自地理空间数据云,在数据库中找到所需的黑龙江省矢量地图数据,下载至个人电脑,完成黑龙江省矢量地图数据的准备,之后再利用MapGIS K9软件结合黑龙江省矢量地图数据,逐一地划分出哈尔滨市、齐齐哈尔市、大庆市、东宁市、肇东市这5个城市的矢量边界。具体工具通过查询–按位置选择实现,最后再将数据导出,进行分析即可[7]。
再将在黑龙江省气象台网站收集整理的这五个城市2013年至2017年这五年的降水总量,平均日照强度以及≥10℃的总活动积温的文本数据导入MapGIS K9中,即完成数据准备阶段的工作。
2.3. 分析方法
将划分出的5个城市属性中建立“积温”“日照强度”“降雨”共3个字段,在系统将5年的文本数据录入完成后转换为矢量数据。之后需要使用到系统工具箱中的“要素转栅格”工具,以达到将“积温”、“日照强度”、“降雨”三者分别由矢量数据转为栅格数据的目的。利用软件的“重分类”工具将区域规划分为:不适宜区、一般适宜区、适宜区。
使用“重分类”工具对“积温”进行了重分类处理,因为大豆这一作物是喜温作物,在温暖的自然环境下才能生长良好,所以在区域划分过程中,将1 d∙℃至2400 d∙℃范围划分为不适宜区(赋值1),2200 d∙℃至2600 d∙℃范围为一般适宜区(赋值2),大于2600 d∙℃范围为适宜区(赋值3);
由于大豆是喜光作物,故在此需要以光照强度作为依据,使用“重分类”工具对“日照强度”进行重分类,规定1 w/m2至195 w/m2范围为不适宜区(赋值1),195 w/m2至200 w/m2为一般适宜区(赋值2),200 w/m2以上为适宜区(赋值3);
大豆生长发育过程中对水的需求量较高,每形成一克营养物质就需消耗水600克到1000克,即需要使用“重分类”工具对“降雨量”进行重分类,将1 mm至35 mm为不适宜区(赋值1),35 mm至40 mm为一般适宜区(赋值2),40 mm至53 mm为适宜区(赋值3)。
最后利用“栅格计算器”工具,将以上三个重分类完成的要素进行相加计算,处理后结果共有:4 (不适宜区)、5 (适宜区)、6 (最适宜区)共3类。
三要素的重分类过程如(图1)和(图2)所示。
Figure 1. Schematic diagram of reclassification process 1
图1. 重分类过程示意图1
Figure 2. Schematic diagram of reclassification process 2
图2. 重分类过程示意图2
3. 研究结果
3.1. 影响因素变化趋势分析
3.1.1. 历年降水量变化
哈尔滨市、齐齐哈尔市、大庆市、东宁市、肇东市历年降水量如(表3)所示;哈尔滨市、齐齐哈尔市、大庆市、东宁市、肇东市历年降水量变化趋势如(图3)所示。
Table 3. 2013~2017 annual precipitation
表3. 2013~2017历年降水量
区域/时间 |
2013年 |
2014年 |
2015年 |
2016年 |
2017年 |
哈尔滨 |
618.13 mm |
432.36 mm |
440.01 mm |
569.24 mm |
491.18 mm |
大庆 |
608.06 mm |
432.43 mm |
440.24 mm |
352.27 mm |
300.08 mm |
齐齐哈尔 |
529.42 mm |
441.27 mm |
462.41 mm |
567.18 mm |
369.05 mm |
东宁 |
685.33 mm |
508.39 mm |
663.37 mm |
721.16 mm |
640.34 mm |
肇东 |
552.25 mm |
448.28 mm |
468.61 mm |
610.33 mm |
398.38 mm |
Figure 3. The trend map of precipitation change in five regions from 2013 to 2017
图3. 2013~2017五大区域降水量变化趋势图
由表3历年降水量及图3五大区域降水量变化趋势图,可看出2013年为五大区域降水量最多的一年,其中东宁市降水量为五大区域中降水量最大值,数值为685.33 mm;2017年为五大区域降水量最少的一年,其中齐齐哈尔市降水量为五大区域中降水量最小值,数值为300.08 mm选择2013年至2017年五年来历年降水量资料数据,对五大区域降水量变化进行趋势预测分析(图中仅标注了东宁市及肇东市的预测趋势线),可看出在五年中五大区域中仅东宁市降水量总体呈上升趋势,且历年降水量都高于其他四个区域,除东宁区域外,剩余四大区域降水量虽然在局部有轻微的上升趋势,但总体上仍呈现下降趋势,且其中齐齐哈尔市降水量下降趋势最为明显,2016年及2017年降水量为五大区域中最少的。
3.1.2. 历年日照强度变化
哈尔滨市、齐齐哈尔市、大庆市、东宁市、肇东市历年月平均日照强度数值如(表4)所示;哈尔滨市、齐齐哈尔市、大庆市、东宁市、肇东市历年月平均变化趋势如(图4)所示。
Table 4. 2013~2017 average monthly sunshine intensity in five regions over the years
表4. 2013~2017五大区域历年月平均日照强度
区域/时间 |
2013年 |
2014年 |
2015年 |
2016年 |
2017年 |
哈尔滨 |
182.97 w/m2 |
195.96 w/m2 |
189.10 w/m2 |
191.09 w/m2 |
193.79 w/m2 |
大庆 |
191.27 w/m2 |
204.80 w/m2 |
200.29 w/m2 |
199.93 w/m2 |
206.43 w/m2 |
齐齐哈尔 |
191.07 w/m2 |
206.24 w/m2 |
199.93 w/m2 |
198.67 w/m2 |
207.51 w/m2 |
东宁 |
186.41 w/m2 |
193.07 w/m2 |
192.89 w/m2 |
188.77 w/m2 |
192.71 w/m2 |
肇东 |
187.84 w/m2 |
202.09 w/m2 |
195.60 w/m2 |
194.52 w/m2 |
202.27 w/m2 |
Figure 4. 2013~2017 trend chart of monthly average sunshine intensity in five major regions over the years
图4. 2013~2017五大区域历年月平均日照强度变化趋势图
由表4历年月平均日照强度及图4五大区域历年月平均日照强度变化趋势图可看出五大区域历年来月平均日照强度一直呈上升趋势,2013年为五大区域历年中月平均日照强度最低年,其中月均日照强度最小值为哈尔滨市日照强度值,数值为182.97 w/m2;2017年为五大区域历年中月平均日照强度最高的年,其中月均日照强度最大值为大庆市日照强度值,数值为207.51 w/m2,选择2013年至2017年五年来历年月平均日照强度资料数据,对五大区域历年月平均日照强度进行趋势预测分析(图中仅标注了斜率最小的东宁市及斜率最大的肇东市的预测趋势线)得出历年中东宁市月平均日照强度趋势最为缓慢,哈尔滨市次之,大庆市与齐齐哈尔市月平均日照强度增加趋势都及其明显,其中齐齐哈尔市月平均日照强度增加趋势最大。
3.1.3. 历年日照强度变化
哈尔滨市、齐齐哈尔市、大庆市、东宁市、肇东市历年活动积温数值如(表5)所示;哈尔滨市、齐齐哈尔市、大庆市、东宁市、肇东市历年降水活动积温变化趋势如(图5)所示。
由表5历年≥10℃活动积温及图5五大区域历年≥10℃活动积温变化趋势图,可看出五大区域中仅东宁市活动积温处于3000 d∙℃水平以下,其最低的一年为2017年,活动积温数值为2424.12 d∙℃;其余四区域活动积温都处于3000 d∙℃水平以上,且哈尔滨≥10℃活动积温最高,明显高于其他四个区域,其最高的一年为2015年,活动积温数值为3309.40 d∙℃;选择2013年至2017年五年来历年活动积温资料数据,对五大区域历年活动积温进行趋势预测分析,其中东宁市区域积温水平有明显的下滑趋势,剩余四大区域中大庆市区域的积温水平增加趋势最高,哈尔滨积温有轻微的减少趋势,肇东及齐齐哈尔两大区域积温无明显的变化趋势。
Table 5. Accumulated temperature ≥10℃ over the years from 2013 to 2017
表5. 2013~2017历年≥10℃积温
区域\时间 |
2013年 |
2014年 |
2015年 |
2016年 |
2017年 |
哈尔滨 |
3235.40 d∙℃ |
3235.50 d∙℃ |
3309.40 d∙℃ |
3246.30 d∙℃ |
3202.53 d∙℃ |
大庆 |
3104.50 d∙℃ |
3197.90 d∙℃ |
3107.70 d∙℃ |
3315.65 d∙℃ |
3192.13 d∙℃ |
齐齐哈尔 |
3086.73 d∙℃ |
3163.30 d∙℃ |
3119.80 d∙℃ |
3310.60 d∙℃ |
3192.24 d∙℃ |
东宁 |
2471.85 d∙℃ |
2484.40 d∙℃ |
2554.40 d∙℃ |
2479.70 d∙℃ |
2424.12 d∙℃ |
肇东 |
3197.61 d∙℃ |
3260.33 d∙℃ |
3184.82 d∙℃ |
3208.01 d∙℃ |
3221.63 d∙℃ |
Figure 5. Trend diagram of accumulated active temperature ≥10℃ in five major regions over the years in 2013~2017
图5. 2013~2017五大区域历年≥10℃活动积温变化趋势图
3.2. 黑龙江省第一积温带大豆适宜性种植区域划分成果
使用MapGIS K9软件,最终完成基于GIS的黑龙江省第一积温带大豆种植气候区域规划,黑龙江省第一积温带大豆种植适宜性分布图如(图6)所示。
Figure 6. The distribution map of soybean planting suitability in the first accumulated temperature zone of Heilongjiang Province
图6. 黑龙江省第一积温带大豆种植适宜性分布图
如图6黑龙江省第一积温带大豆种植适宜性分布图所示结果,在这五大区域划分中,东宁市因其气候条件:≥10℃活动积温在五区域中最低,日照强度偏低,2013年至2017年五年平均日照强度为190.77 w/m2被系统划分为不适宜区;系统结合肇东市气候条件:≥10℃活动积温较高,2400 < T ≤ 2600将其划分为一般适宜区;系统结合气候条件:≥10℃活动积温高,T > 2600,且齐齐哈尔市、大庆市、哈尔滨市三个区域日照强度高,降水多,系统将齐齐哈尔市、大庆市、哈尔滨市划分为最适宜区。
4. 结论
在当前的研究中,我们完成了基于GIS的黑龙江省第一积温带大豆种植气候区域规划,系统地收集和分析了相关气候气象资料,包括降水量、日照强度和活动积温的历年变化趋势。通过使用MapGIS软件,我们不仅完成了对五个区域的划分,还明确了大豆种植的适应性区域。这一过程体现了我们对已有研究成果的继承与发展,但也存在可进一步优化的空间。
首先,与前人的研究相比,我们在数据收集和处理上进行了更为细致的分析,尤其是在气候因素的多维度评估上[8]。尽管之前的工作在应用基本工具上取得了一定成果,但往往缺乏对复杂气候条件的深入解析[9]。本研究通过对积温、降雨和日照强度进行综合分析,力求提供更加全面的区域适应性评估,为未来的种植决策提供更可靠的依据。
其次,前期的研究多依赖于MapGIS的基本功能,而本研究则计划在未来的工作中引入更高级的GIS工具和方法。例如,我们希望结合多元统计分析和机器学习技术,利用大数据和实时气象数据进行动态评估。这不仅可以提高区域划分的精确度,还能使规划结果更具适应性和时效性,满足不同气候条件下的种植需求。
此外,我们还意识到,与前人的工作相比,系统工具的运用尚显单一[10]。因此,未来的研究将致力于开发更为复杂的分析模型,整合更多数据源,以实现对气候变化的快速响应和调整。这一目标不仅将推动区域规划的科学化、精细化,也为其他地区的类似研究提供了可借鉴的经验。
通过这些改进,我们的研究不仅是对现有工作的延续,更是对其不足之处的积极回应。我们期待在未来的研究中,能够深入探索气候变化对农业的影响,实现更为准确和动态的气候区域规划,从而推动黑龙江省大豆种植的可持续发展。这不仅是对我们研究价值的提升,更是为地方农业发展提供科学支持的重要举措。