中国农作物气候适宜性研究进展

吕子鹤, 丁松爽, 卢瑞琳

中国农学通报. 2020, 36(24): 78-84

PDF(1116 KB)
PDF(1116 KB)
中国农学通报 ›› 2020, Vol. 36 ›› Issue (24) : 78-84. DOI: 10.11924/j.issn.1000-6850.casb20190700470
资源·环境·生态·土壤·气象

中国农作物气候适宜性研究进展

作者信息 +

Climate Suitability of Crops in China: Research Progress

Author information +
History +

摘要

探讨中国农作物气候适宜性研究现状,为更好地利用气候资源,以及为农作物种植方案和保护措施的制定提供理论依据。本研究总结了目前农作物气候适宜性研究的常用方法,讨论了各方法的原理、优缺点和适用范围,并在此基础上综述了中国水稻、小麦、玉米、棉花、烟草等5种农作物气候适宜性的研究进展。结果发现,中国农作物气候适宜性研究已逐渐深入,在充分利用气候资源、合理安排生产措施等方面发挥重要作用。MAXENT模型、GIS技术、积分回归法、层次分析法等技术各有优缺点,现有方法的综合运用及新技术的不断完善将为农作物气候适宜性研究提供新思路。

Abstract

To discuss the status of climate suitability of crops in China, make better use of climate resources, and provide a theoretical basis for planting plans and protection measures for crops, we summarized the common methods used in the study of climate suitability of crops, and discussed the principles, advantages and disadvantages of each method and the scope of their application; on the basis of this, we reviewed the research progress on the climate suitability of 5 crops, including rice, wheat, corn, cotton and tobacco in China. The results show that: the research on the climate suitability of crops in China has been gradually deepened, and played an important role in making full use of climate resources and rational production measures. MAXENT model, geographic information system, integral regression method, analytic hierarchy process and other technologies have their own advantages and disadvantages. The comprehensive application of these methods and the continuous improvement of new technologies will provide new ideas for the study of climate suitability of crops.

关键词

气候变化 / 气候适宜性 / 研究方法 / MAXENT模型 / GIS技术 / 积分回归法 / 层次分析法

Key words

climate change / climate suitability / research method / MAXENT model / geographic information system / integral regression method / analytic hierarchy process

引用本文

导出引用
吕子鹤 , 丁松爽 , 卢瑞琳. 中国农作物气候适宜性研究进展. 中国农学通报. 2020, 36(24): 78-84 https://doi.org/10.11924/j.issn.1000-6850.casb20190700470
Lv Zihe , Ding Songshuang , Lu Ruilin. Climate Suitability of Crops in China: Research Progress. Chinese Agricultural Science Bulletin. 2020, 36(24): 78-84 https://doi.org/10.11924/j.issn.1000-6850.casb20190700470

参考文献

[1]
秦大河, 罗勇, 陈振林, 等. 气候变化科学的最新进展:IPCC第四次评估综合报告解析[J]. 气候变化研究进展, 2007(06):311-314.
[2]
居煇, 许吟隆, 熊伟. 气候变化对我国农业的影响[J]. 环境保护, 2007(11):71-73.
农业是气候变化最敏感的领域之一.气候变化对我国农业会产生或利或弊的影响,但基本上是以不利影响为主,并由于作物种类、区域、环境条件等因素的不同而不同.中国是一个农业大国,也是一个人口大国,农业生产特别是粮食生产直接关系到社会的稳定和可持续发展.
[3]
Davidson Thomas. On the Climate of South California and Its Suitability as a Health Resort[J]. Glasgow medical journal, 1893,39(6).
[4]
郭柏林. 长兴岛柑桔气候适宜性的探讨[J]. 上海农业科技, 1983(06):20-21.
[5]
王雪娥. 玉米气候适宜度动态模型的建立和应用[J]. 南京气象学院学报, 1992(02):63-72.
[6]
赵峰, 千怀遂, 焦士兴. 农作物气候适宜度模型研究——以河南省冬小麦为例[J]. 资源科学, 2003(06):77-82.
[7]
聂志刚, 任新庄, 李广, 等. 基于APSIM的黄土丘陵区旱地小麦气候适宜性评价[J]. 中国农业气象, 2017,38(06):369-377.
[8]
陈炳东. 气候变化对农作物适宜度影响研究进展[J]. 北京农业, 2012(18):153.
[9]
冶明珠, 郭建平, 蒋跃林, 等. 气候变化对农作物气候适宜度影响研究进展[J]. 安徽农业科学, 2011,39(15):9104-9105,9134.
[10]
罗怀良, 陈国阶, 朱波. 农业生态气候适宜度研究进展[J]. 中国农业资源与区划, 2004(01):29-33.
[11]
Phillips S J, Anderson R P, E, Schapire R. Maximum entropy modeling of species geographic distributions[J]. Ecological Modelling, 2006,190:231-259.
[12]
王运生, 谢丙炎, 万方浩, 等. 相似穿孔线虫在中国的适生区预测[J]. 中国农业科学, 2007(11):2502-2506.
【目的】相似穿孔线虫(Radopholus similis)是中国一类检疫对象,随着国际贸易日益加剧,该线虫传入中国的风险也日益加大,急需对其适生性进行分析。【方法】利用已广泛用于物种潜在分布预测的GARP与MAXENT两种生态位模型预测相似穿孔线虫在中国的适生区,并对预测结果进行了阈值依赖和非阈值依赖比较。【结果】两个模型均能较好地预测相似穿孔线虫的潜在地理分布;预测在中国适生区主要分布在南方各省,包括海南、广东、福建、云南、台湾等省区;经检验MAXENT预测结果较GARP好。【结论】相似穿孔线虫可随花卉苗木的国际运输传入中国,并可在中国南方沿海地区及花卉进出口大省云南省等地区发生;上述地区的进出口口岸应加强对相似穿孔线虫的检疫工作。
[13]
Phillips S J, Miroslav D. Modeling of species distributions with Maxent: new extensions and a comprehensive evaluation[J]. Ecography, 2008,31(2):161-175.
[14]
高蓓, 卫海燕, 郭彦龙, 等. 应用GIS和最大熵模型分析秦岭冷杉潜在地理分布[J]. 生态学杂志, 2015,34(03):843-852.
[15]
吴庆明, 王磊, 朱瑞萍, 等. 基于MAXENT模型的丹顶鹤营巢生境适宜性分析——以扎龙保护区为例[J]. 生态学报, 2016,36(12):3758-3764.
明确物种繁殖栖息地的潜在分布对于受胁物种的保护与科学管理具有十分重要的意义。扎龙保护区是国际濒危物种丹顶鹤西部迁徙种群的主要繁殖地之一,近年来繁殖种群数量波动较大。通过GIS空间技术平台,利用MAXENT模型,结合2012-2013年丹顶鹤(Grus japonensis)营巢分布点和环境特征变量,对扎龙保护区丹顶鹤营巢生境的适宜性进行了分析。结果表明:模型的评价效果达到优秀水平(训练集AUC=0.901);绿度(34.9%)、距道路距离(17.1%)、土壤湿度(16.3%)、海拔(15.7%)是丹顶鹤营巢生境的主要环境特征变量;核心区中部和南部是丹顶鹤主要的适宜营巢分布区,缓冲区和实验区有少量的适宜营巢区分布,丹顶鹤适宜营巢生境面积为35857.35hm2,占保护区总面积的17.07%。建议在重点管理核心区的同时,对缓冲区和实验区的小面积芦苇沼泽给予更多关注。
[16]
王茹琳, 王明田, 罗家栋, 等. 基于MaxEnt模型的美味猕猴桃在中国气候适宜性分析[J]. 云南农业大学学报:自然科学, 2019,34(03):522-531.
[17]
朱耿平, 乔慧捷. Maxent模型复杂度对物种潜在分布区预测的影响[J]. 生物多样性, 2016,24(10):1189-1196.
[18]
贾翔, 王超, 金慧, 等. 基于优化的MaxEnt模型评价红松适宜分布区[J/OL]. 生态学杂志,1-9. Https://doi.org/10.13292/j.1000-4890.201908.017.
[19]
常晋义. GIS技术与应用[J]. 大同高等专科学校学报, 1994(04):95-97.
[20]
梁轶, 柏秦凤, 李星敏, 等. 基于GIS的陕南茶树气候生态适宜性区划[J]. 中国农学通报, 2011,27(13):79-85.
[21]
汤思遥. 基于RS与GIS技术的景观格局分析及生态适宜性评价——以黄石市阳新县为例[J]. 价值工程, 2019,38(12):139-141.
[22]
李东颖, 周苏刚, 王莉霞, 等. 基于MR-GIS的油橄榄种植气候适宜性分析[J]. 地理空间信息, 2018,16(12):68-70,10.
[23]
宛公展, 杨光斌. 水稻产量与气候条件的积分回归分析[J]. 天津农业科学, 1982(03):22-25.
[24]
裴顺祥, 郭泉水, 贾渝彬, 等. 保定市8种乔灌木开花始期对气候变化响应的积分回归分析[J]. 北京林业大学学报, 2015,37(07):11-18.
[25]
肖云云, 欧阳金琼. 塔里木河流域农业气候资源综合评价与区域比较——基于模糊层次分析法[J]. 塔里木大学学报, 2018,30(03):87-94.
[26]
赵辉, 米鸿涛, 杜子璇. 基于AHP的模糊综合评判在茶叶气候种植区划中的应用[J]. 气象科技, 2013,41(06):1134-1137.
[27]
田辉, 金洪涛, 孙岐发, 等. 基于层次分析法的盘锦湿地生态评价[J]. 地质与资源, 2018,27(03):268-271,287.
[28]
Tang H J, Van Ranst E. Testing of fuzzy set theory in land suitability assessment for rainfed grain maize production[J]. Pedologie, 1992,42:129-147.
[29]
Burrough P A. Fuzzy mathematical methods for soil survey and land evaluation[J]. Journal of Soil Sciences, 1989,40:477-492.
[30]
李凡修, 陈武. 模糊集对分析在城市绿地景观生态综合评价中的应用[J]. 新疆环境保护, 2002(02):32-34.
[31]
姚宁, 周元刚, 宋利兵, 等. 不同水分胁迫条件下DSSAT-CERES-Wheat模型的调参与验证[J]. 农业工程学报, 2015,31(12):138-150.
作物模型为人们认识旱区农业生境过程并对其进行调控提供了一种有效的工具。为了探讨小麦生长模拟模型DSSAT-CERES-Wheat能否准确模拟水分胁迫条件下旱区冬小麦的生长发育和产量形成过程,同时确定参数估计和模型验证的最优方案,该研究进行了连续两季(2012.10-2013.06和2013.10-2014.06)的冬小麦分段受旱田间试验。试验将冬小麦整个生育期划分为越冬、返青、拔节、抽穗和灌浆5个主要生长阶段,每相邻两个生长阶段连续受旱,形成4个不同的受旱时段水平(D1-D4),根据小麦生育期的需水量,设置灌水定额分别为40和80 mm 2个水平(I1和I2),共形成8个处理,每处理3次重复,在遮雨棚内采用裂区试验布置,此外在旁边设置1个各生育期全灌水的对照处理。文中设置了5套不同的参数估计和验证方案,利用DSSAT-GLUE参数估计模块得到不同的参数估计结果。通过对比分析冬小麦物候期、单粒质量、生物量、产量、以及土壤水分含量的模拟值和实测值之间的差异,以确定利用DSSAT-CERES-Wheat模型模拟旱区冬小麦生境过程的精度。结果表明,参数P1V(最适温度条件下通过春化阶段所需天数)和G3(成熟期非水分胁迫下单株茎穂标准干质量)具有较强的变异性,变异系数分别为19.07%和16.34%,受基因型-环境互作的影响较大,而其他参数的变异性则较弱,变异系数均小于10%;DSSAT-GLUE参数估计工具具有较好的收敛性,不同参数估计方案所得的参数值具有一定的一致性;不同的参数估计方案所得的模型输出结果有较大差异,其中参数估计方案1(利用两季试验中的充分灌溉处理CK数据进行参数估计,其他不同阶段受旱处理数据进行验证)的模型校正和验证精度最高,其中模型校正的绝对相对误差(absolute relative error,ARE)和相对均方根误差(relative root mean squared error, RRMSE)分别为4.89%和5.18%。在冬小麦抽穗期和灌浆期受旱时,DSSAT-CERES-Wheat模型可以较好地模拟小麦的生长发育过程以及土壤水分的动态变化,但是在越冬期和返青期受旱时,模拟结果相对较差,并且随着受旱时段提前和受旱程度的加重,模拟精度将变得更低。此外,该模型无法模拟由不同水分胁迫造成的冬小麦物候期差异,需要对模型进行相应的改进。交叉验证表明DSSAT-CERES-Wheat模型模拟该研究中不同水分胁迫条件下冬小麦生长和产量的总体性误差在15%~18%左右。总之,DSSAT-CERES-Wheat模型在模拟旱区冬小麦生境过程时存在着一定的局限性,若要更广泛地将该模型应用在中国干旱半干旱地区的冬小麦生产管理和研究,有必要对冬小麦营养生长阶段前期的水分胁迫响应机制和模拟方法进行进一步的深入研究。
[32]
聂志刚, 李广, 王钧, 等. 基于APSIM模型旱地小麦叶面积指数相关参数的优化[J]. 中国农业科学, 2019,52(12):2056-2068.
【目的】 模型的有效应用依赖于参数的快速、准确估算。本研究拟解决作物生长模型参数本土化率定过程中运算量大、耗时长、精度低、效率低的问题。【方法】 依据甘肃省定西市安定区2个试验点(李家堡镇麻子川村和凤翔镇安家沟村),多年(2002—2005年和2015—2017年)大田试验数据以及定西市安定区1971—2017年气象资料,利用混合蛙跳算法智能的迭代搜索原理,对APSIM模型旱地小麦叶面积指数相关参数进行了优化,并采用相关性分析方法对模型校正结果进行检验。【结果】 利用青蛙群体即相对独立又合作协调的子群内局部深度搜索与子群间全局信息交流生物进化学习策略,有效提高了运算的速度,实现了对APSIM模型中与旱地小麦叶面积指数相关参数的快速、准确估算。相关参数主要包括:主茎上节出现所需的热时间间隔、小麦出苗后初始化的节数、小麦出苗后初始化的叶片数、小麦出苗后初始化的叶面积指数、某日正在生长的节数和最大比叶面积。分别使用穷举试错法所得参数值和混合蛙跳算法所得参数值模拟叶面积指数,参数优化后,叶面积指数模拟值和实测值之间的RMSE(root mean square error)平均值由0.069降低到0.027,NRMSE(normalized root mean square error)平均值由8.09%降低到4.56%,ME(model effective index)平均值由0.979提高到0.993。【结论】 相对于参数率定常用穷举试错法,混合蛙跳算法具有自发学习特征的智能迭代行为,实现了参数的自动率定,提高了效率。基于该算法进行APSIM模型旱地小麦叶面积指数相关参数的优化,使得模型对叶面积指数的模拟精度显著提高,证明该算法的使用对作物生理生态系统复杂模型的校正效果良好,为改善模型参数率定过程存在的运算量大、耗时长、精度低、效率低的缺点提供了一种行之有效的方法。
[33]
成林, 刘荣花, 王信理. 气候变化对河南省灌溉小麦的影响及对策初探[J]. 应用气象学报, 2012,23(05):571-577.
[34]
卢为国, 李卫东, 梁慧珍, 等. DTOPSIS法综合评价大豆新品种的初步探索[J]. 中国油料作物学报, 1998(03):23-27.
[35]
郝君明, 赵军, 刘红兵, 等. 基于FIS与GIS方法的山西省玉米气候适宜性评价[J]. 安徽农业科学, 2008(23):10129-10130,10142.
在ArcGIS和Matlab软件的支持下,采用网格分析与模糊推理集成的方法实现了山西省玉米气候适宜性评价.结果表明,模糊推理与GIS格网分析相结合,在农业气候适宜性评价中具有良好的应用效果.
[36]
徐中强, 陈俊华, 邵俊明, 等. 基于AHP和GIS的重庆市避暑自然适宜性评价[J]. 贵州师范大学学报:自然科学版, 2016,34(01):14-21.
[37]
章秀福, 王丹英, 方福平, 等. 中国粮食安全和水稻生产[J]. 农业现代化研究, 2005(02):85-88.
[38]
段居琦, 周广胜. 我国单季稻种植区的气候适宜性[J]. 应用生态学报, 2012,23(02):426-432
[39]
段居琦, 周广胜. 中国双季稻种植区的气候适宜性研究[J]. 中国农业科学, 2012,45(02):218-227.
[40]
凌霄霞, 张作林, 翟景秋, 等. 气候变化对中国水稻生产的影响研究进展[J]. 作物学报, 2019,45(03):323-334.
[41]
郭夏宇, 艾治勇, 龙继锐. 气候变化对中国水稻生产的影响[J]. 湖南农业科学, 2018(07):51-54.
[42]
徐敏, 徐乐, 巫丽君, 等. 江苏省水稻产量结构变化特征及与综合气候适宜度的关系[J/OL]. 江苏农业科学,1-6. Https://doi.org/10.15889/j.issn.1002-1302.2019.15.024.
[43]
韩芳玉, 张俊飚, 程琳琳, 等. 气候变化对中国水稻产量及其区域差异性的影响[J]. 生态与农村环境学报, 2019,35(03):283-289.
[44]
千怀遂, 焦士兴, 赵峰. 河南省冬小麦气候适宜性变化研究[J]. 生态学杂志, 2005(05):503-507.
[45]
宋洪秋, 千怀遂, 俞芬 等. 全球气候变化下淮河流域冬小麦气候适宜性评价[J]. 自然资源学报, 2009,24d(5):890-897.
[46]
成林, 李彤霄, 刘荣花. 主要生育期气候变化对河南省冬小麦生长及产量的影响[J]. 中国生态农业学报, 2017,25(06):931-940.
[47]
许彦平, 万信, 姚晓红, 等. 甘肃河东地区小麦条锈病越冬气候适宜性研究[J]. 中国农学通报, 2018,34(12):142-145.
[48]
魏瑞江, 宋迎波, 王鑫. 基于气候适宜度的玉米产量动态预报方法[J]. 应用气象学报, 2009,20(05):622-627.
[49]
何奇瑾, 周广胜. 我国春玉米潜在种植分布区的气候适宜性[J]. 生态学报, 2012,32(12):3931-3934,3937-3939.
[50]
赵锦, 杨晓光, 刘志娟, 等. 气候变化对东北三省春玉米气候适宜性的影响[J]. 中国农业科学, 2014,47(16):3143-3156.
【目的】研究气候变化背景下中国东北三省春玉米气候适宜性的变化特征,为东北地区春玉米种植的合理布局提供科学依据。【方法】本文以1981年为时间节点,把1961—2010年分为两个时间段,基于东北三省74个气象站点1961—2010年的观测资料,根据农业气象学指标,在分析1981—2010年较1961—1980年东北三省春玉米可能种植北界变化的基础上,利用APSIM-Maize模型模拟春玉米可能种植区域内各气象站点逐年的雨养产量,结合统计学方法,分析春玉米雨养产量高产性和稳产性的变化特征,并综合得到气候变化背景下中国东北三省春玉米的气候适宜区分布的变化特征。【结果】(1)与1961—1980年相比,1981—2010年春玉米的可能种植北界向北移动了158.3—285.8 km,春玉米可能种植面积增加了3.87×104 km2,占东北三省土地面积的4.91%。(2)1981—2010年,东北三省春玉米雨养产量的最高产区、高产区和次高产区面积占可能种植面积比例由81.14%增加为86.66%,其中最高产区和高产区面积比例由36.61%减少为34.82%,而次高产区则由44.53%增加到51.85%,低产区面积占研究区域面积的比例由18.86%减少为13.34%。研究区域内雨养产量的单产减少40 kg•hm-2,但由于春玉米可能种植区域总面积的扩大,特别是高产区和次高产区面积的增加,研究区域内雨养产量的总产增加了7.0%。(3)东北三省春玉米雨养产量的最稳产区、稳产区和次稳产区面积占可能种植面积的比例由80.20%增加为89.28%,且最稳产区和稳产区面积比例由40.97%增加为49.97%,而低稳产区面积占研究区域面积的比例由19.80%减少为10.72%。(4)东北三省春玉米气候适宜区和次适宜区面积占可能种植面积的比例由61.09%增加为83.00%,但最适宜区面积比例由18.83%减少为6.67%,可种植区面积占研究区域面积的比例由20.08%减少为10.33%。研究区域内春玉米可稳定获得的雨养产量的单产总体下降了171 kg•hm-2,但由于春玉米可能种植区域总面积的增加,特别是适宜区和次适宜区面积的增加,研究区域内可稳定获得的雨养产量的总产增加了2.6%。【结论】全球气候变暖背景下,东北三省春玉米种植北界明显向北向西移动,春玉米的可能种植面积增加;在春玉米可能种植区域内,如果不考虑品种和栽培管理措施的适应,雨养产量的最高产区面积所占比例缩小,春玉米雨养产量的单产下降,但由于可能种植面积的增加,东北三省春玉米雨养产量的总产增加;春玉米雨养产量的稳定性增加,其中最稳产区和稳产区面积占研究区域面积的比例增加;在春玉米可能种植的区域内,春玉米的气候最适宜区面积减少明显,研究区域内春玉米可稳定获得的雨养产量的单产下降,但由于春玉米可能种植面积的增加,特别是适宜区和次适宜区面积的增加,可稳定获得的雨养产量的总产增加。
[51]
初征, 郭建平. 未来气候变化对东北玉米品种布局的影响[J]. 应用气象学报, 2018,29(02):165-176.
[52]
孙良斌. 中国棉花主产区的优势动态比较[J]. 湖北农业科学, 2014,53(13):3195-3198.
[53]
程维明, 柴慧霞, 周成虎, 等. 新疆地貌空间分布格局分析[J]. 地理研究, 2009,28(05):1157-1169.
[54]
韩慧君. 气候生态因素对棉花产量与纤维品质的影响[J]. 中国农业科学, 1991(05):23-29.
[55]
陶满德. 试论中国棉花主产区形成的历史原因[J]. 江西师范大学学报, 1994(04):55-60.
[56]
徐德源. 新疆农业气候资源特点及其利用和开发(一)[J]. 新疆农业科技, 1986(04):8-12.
[57]
徐德源. 新疆农业气候资源特点及其利用和开发(二)[J]. 新疆农业科技, 1986(05):30-32.
[58]
宋艳玲, 张强, 董文杰. 气候变化对新疆地区棉花生产的影响[J]. 中国农业气象, 2004(03):15-20.
[59]
赵黎, 李文博, 吾米提·居马泰, 等. 气候变化对新疆棉花种植布局与生长发育的影响[J]. 新疆农垦科技, 2018,41(07):7-10.
[60]
李迎春, 谢国辉, 王润元, 等. 北疆棉区棉花生长期气候变化特征及其对棉花发育的影响[J]. 干旱地区农业研究, 2011,29(2):253-258.
[61]
新疆农学会. 北疆棉区划分及相应对策[J]. 新疆农业科学, 1995,32(1):5-8.
[62]
姚源松. 新疆棉花区划新论[J]. 中国棉花, 2001,28(2):2-5.
[63]
黄滋康, 崔读昌. 中国棉花生态区划[J]. 棉花学报, 2002,14(3):185-190.
[64]
严志丹, 王兰, 李美华, 等. 基于DTOPSIS法的塔里木河流域棉花种植适宜性区划研究[J]. 中国农学通报, 2015,31(11):210-216.
[65]
陈海生, 刘国顺, 刘大双, 等. GIS支持下的河南省烟草生态适宜性综合评价[J]. 中国农业科学, 2009,42(07):2425-2433.
[66]
刘国顺. 烟草栽培学[M]. 北京: 中国农业出版社,年: 2003:58-67.
[67]
胡雪琼, 黄中艳, 朱勇, 等. 云南烤烟气候类型及其适宜性研究[J]. 南京气象学院学报, 2006(04):563-568.
[68]
邵岩. 基于GIS的云南烤烟种植生态适宜性区划[D]. 长沙:湖南农业大学, 2008:72-75.
[69]
莫静静, 刘国顺, 叶协锋, 等. 基于GIS和AHP的南阳市烟草种植生态适宜性评价[J]. 河南农业大学学报, 2009,43(03):331-334,348.
[70]
陈海生, 刘国顺, 刘大双, 等. GIS支持下的河南省烟草生态适宜性综合评价[J]. 中国农业科学, 2009,42(07):2425-2433.
[71]
张久权, 梁洪波, 董建新, 等. 基于GIS和模糊集理论的四川烤烟生态适应性评价[J]. 中国烟草科学, 2016,37(03):8-14.

基金

湖北省烟草公司科技项目“鄂西南优质雪茄烟叶生产关键技术研究”(027Y2018-003)

版权

版权所有,未经授权,不得转载、摘编本刊文章,不得使用本刊的版式设计。
PDF(1116 KB)

文章所在专题

水稻

Accesses

Citation

Detail

段落导航
相关文章

/