PDF(2144 KB)
Quality Evaluation of Cultivated Land and Soil Nutrient Status in Ningjinpo-Daluze Low Plain Area
SUNYuqiang, FANGZhu, XIAOHui, ZHANGJianfa, JIAYing, LIDanlin, ZHAOYanna
Journal of Agriculture ›› 2025, Vol. 15 ›› Issue (9) : 21-29.
PDF(2144 KB)
PDF(2144 KB)
Quality Evaluation of Cultivated Land and Soil Nutrient Status in Ningjinpo-Daluze Low Plain Area
Through the 2022 Hebei Province Farmland Quality Evaluation Project, data from farmland quality monitoring points in eight counties of the Ningjinbo-Daluze low plain area were summarized and analyzed to evaluate the quality of farmland and soil nutrient status in the region. Using the grid method to evenly arrange sample evaluation units, the five-point method was adopted to collect cultivated soil samples from the surface layer (0-20 cm), and laboratory tests were conducted according to the soil testing series standards. Referring to the National Farmland Quality Grade Evaluation Index System of the Ministry of Agriculture and Rural Affairs, the weight of the indicators was determined. For numerical indicators of organic matter, total nitrogen, available phosphorus, and available potassium, the membership degree was determined by establishing functional relationships. For other conceptual indicators, the Delphi method expert scoring was used to give the membership degree, and the cumulative method was used to calculate the comprehensive index of arable land quality to obtain the arable land quality grade. The results showed that the average quality grade of cultivated land in the area was 3.51, with the main cultivated land grades being third and fourth grade. The third grade cultivated land was the main body, with an area of 140191 hm2, accounting for 49.44% of the total cultivated land area in the area. The overall level of nutrients in the cultivated land in this area was relatively high, all of which were above level three. Nutrients were evenly distributed within the region, and the maximum values of organic matter, total nitrogen, available phosphorus, and available potassium appeared around the ancient Daluze marsh land. By analyzing the factors affecting the quality evaluation of cultivated land, such as irrigation guarantee, the texture and configuration of the plow layer, and the reserve of organic matter content, measures and suggestions for maintaining the quality level of cultivated land in the region were proposed, providing reference for soil nutrient monitoring, improving cultivated land quality, protecting ecological balance, and maintaining sustainable agricultural development in other low plain areas in Hebei Province.
Ningjinbo-Daluze / low plain area / quality of cultivated land / soil nutrient / comprehensive index
| [1] |
邹逸麟. 黄淮海平原历史地理[M]. 合肥: 安徽教育出版社,1993.
|
| [2] |
王会昌. 河北平原的古代湖泊.地理集刊(第18号)[M]. 北京: 科学出版社,1987.
|
| [3] |
郭盛乔, 王玉海, 杨丽娟, 等. 宁晋泊地区冰消期以来的气候变化[J]. 第四纪研究, 2000(5):490.
|
| [4] |
石超艺. 明代以来大陆泽与宁晋泊的演变过程[J]. 地理科学, 2007, 27(3):414-419.
采用历史文献及地名学的研究方法,对明代以后海河南系平原地区的湖沼演变过程进行分析。结果表明,明代至16世纪后期,该区主要有南起任县、北至宁晋的广阔湖泊——大陆泽。17世纪初,大陆泽分解成南泊大陆泽与北泊宁晋泊。分解时,大陆泽比宁晋泊更为深广。此后大陆泽不断缩小,宁晋泊却持续扩大。至19世纪初,宁晋泊已3倍于大陆泽。19世纪二三十年代之后,两湖泊同时发生显著萎缩。19世纪末,大陆泽就已走向消涸,宁晋泊此时也已"不显泊形",并在20世纪初同样走向终结。
|
| [5] |
常全旺. 明清时期大陆泽、宁晋泊湖区垦殖情况及影响[J]. 邢台学院学报, 2011, 26(1):8-10.
|
| [6] |
李俊杰, 李建平, 梅冬. 新形势下高标准农田建设管理政策存在的问题及建议[J]. 中国农业资源与区划, 2022, 43(5): 84-92.
|
| [7] |
农业农村部. 2019年全国耕地质量等级情况公报[EB/OL].(2020-02-06)[2024-3-29]. http://www.ntjss.moa.gov.cn/zcfb/202006/P020200622573390595236.pdf.
|
| [8] |
农业农村部. 耕地质量调查监测与评价办法[EB/OL].(2016-06-21)[2024-04-29]. http://www.moa.gov.cn/govpublic/ZZYGLS/201607/t20160722_5215391.htm.
|
| [9] |
段友春. 基于GIS的鲁东南典型丘陵平原区农田土壤养分空间变异特征及影响因素研究[J]. 中国农学通报, 2021, 37(23):61-68.
以鲁东南丘陵平原区的某镇为例,利用采集的200个农田表层土壤样点数据,运用地统计学方法研究了土壤中有机质(SOM)、全氮(TN)、碱解氮(AN)、速效磷(AP)和速效钾(AK)含量的空间变异特征及其影响因素。结果表明:SOM、TN、AN、AP和AK的平均含量分别为13.52 g/kg、0.75 g/kg、87.27 mg/kg、50.55 mg/kg和103.20 mg/kg,土壤养分含量总体属于稍缺水平,变异系数依次按大小排序为AP>AN>AK>TN>SOM,全部为中等变异程度。地统计分析表明,5种土壤养分的变异函数符合指数模型,通过普通克里格法对土壤养分进行空间插值来绘制了养分含量等级空间分布图,总体来说研究区土壤的有机质和全氮含量水平都较低,其他有效组分在含量等级及分布范围相对于有机质和全氮的状况要稍好。地形地貌、土壤类型和土地利用方式等因素对研究区土壤养分的空间变异程度有显著影响,其中地形地貌、土壤类型等结构性因素的影响程度较强。
|
| [10] |
李娜, 陈署晃, 赖宁, 等. 基于GIS和地统计的伊犁州农田土壤养分空间变异特征分析[J]. 新疆农业科学, 2019, 56(11):2096-2103.
【目的】 研究土壤养分的空间变异程度及分布规律,为干旱区农业施肥区划及精准农业的推广提供理论依据。【方法】 以新疆伊犁州为研究区,运用GIS技术与地统计学相结合的方法,研究伊犁州农田土壤养分的空间分布特征及变异规律。【结果】 (1)AP、SOM均属于强变异强度,变异系数分别为0.82、0.79;其它土壤肥力指标变异系数处于0.60~0.67,均属于中等变异强度。(2)伊犁州土壤TN、SOM含量总体较高,较丰富以上的样点比例分别占66.66%、61.07%;其次为AK、AP,较丰富以上的样点比例分别占50.85%、49.37%,AN含量相对较低,缺乏和较缺乏的样点比例超过65.72%。(3)TN、SOM、AN表现为极强的空间相关性,其块基比分别为0.107、0.002、0.000 07;AP的块基比为0.403,在区域内表现为中等强度的空间自相关性;速效钾的块基比为0.755,在区域内表现为较弱的空间自相关性。(4)该地区土壤养分含量在空间分布上多呈斑块状格局。5个养分指标中,TN、AN含量的空间分布大体相似,含量较高的区域主要分布伊犁州的东部;SOM、AP含量的空间分布大体相似,含量较高的区域主要分布伊犁州的南部;AK含量的空间分布相对均衡。【结论】 研究区AK、SOM、AP平均含量均处于较高水平;TN、AN平均含量均处于较低水平。土壤养分含量空间分布总体呈现南高北低、东高西低的格局。
|
| [11] |
赵志尚, 胡伟艳, 魏安奇. 湖北省耕地多功能变化及障碍因素诊断[J]. 江苏农业科学, 2018, 46(5):268-272.
|
| [12] |
李开丽, 檀满枝, 密术晓, 等. 基于土壤剖面质地构型的土壤质量评价-以河南省封丘县为例[J]. 土壤, 2016, 48(6):1253-1260.
|
| [13] |
谢英凯, 杨颢, 胡月明, 等. 县域耕地质量等别监测分区布点研究[J]. 农业资源与环境学报, 2020, 37(6):845-855.
|
| [14] |
桂苗, 汪甜甜, 马友华, 等. 合肥市农田土壤质地空间特征分析与耕地质量提升[J]. 中国农学通报, 2022(38):114-119.
|
| [15] |
马倩倩, 董博, 许旺旺, 等. 干旱区耕地质量等级评价及土壤养分与盐渍化的分析研究-以民勤绿洲为例[J]. 干旱区地理, 2021, 44(2):514-524.
为整体掌握民勤绿洲耕地质量、土壤养分以及盐渍化状况等,充分利用县域耕地资源,分析研究民勤绿洲耕地质量等级划分及其耕地土壤养分和盐渍化现状,为及时准确掌握绿洲耕地地力状况提供科学依据,同时为绿洲耕地土壤养分、盐渍化障碍提供改良措施。基于农业农村部《耕地质量调查监测与评价办法》（农业部2016年第2号）和《耕地质量等级》国家标准（GB/T33469-2016）的相关要求,构建民勤绿洲县域评价体系,通过综合运用空间分析、模糊数学、层次分析和综合指数等方法,对民勤绿洲耕地质量等级进行综合评价。结果表明：（1） 民勤绿洲耕地质量平均等级为3.26等;分布有一至七等级地,八、九、十等级地没有分布;以三、四等级地为民勤绿洲耕地的主要等级分布,面积占比分别为31.63%、28.27%。泉山区中部和坝区灌溉农业区多分布一、二等级地。湖区灌溉区多分布三、四等地。在绿洲边缘带较多分布五、六、七等地。（2） 民勤绿洲耕地土壤有机质平均含量在10.00~15.00 g·kg<sup>-1</sup>之间,全氮平均含量在0.50~1.00 g·kg<sup>-1</sup>之间,总体处于较低水平。有效磷和速效钾总体处于中等水平。昌宁区、泉山区与坝区交界东北部以及湖区镇（昌宁镇、东湖镇、南湖镇、泉山镇、西渠镇）耕地土壤养分平均含量较均低。（3） 民勤绿洲盐渍化耕地的面积为32240.36 hm<sup>2</sup>,占民勤绿洲耕地总面积的30.09%;主要以轻度盐渍化为主。轻度盐渍化多分布在湖区红沙岗镇、红沙梁镇、西渠镇等镇;中度盐渍化分布在昌宁镇以及环河区沿风沙线一带。其评价结果对民勤绿洲耕地资源的科学管理和可持续利用以及绿洲农业的发展具有重要意义。
|
| [16] |
曹婧, 陈怡平, 毋俊华, 等. 近40年陕西省农田土壤有机质时空变化及其影响因素[J]. 地球环境学报, 2022, 13(3):331-343.
|
| [17] |
范铭丰. 三明市耕地土壤交换性镁含量空间分布特征及影响因素分析[J]. 中国农学通报, 2021, 37(29):71-77.
为了解三明市耕地土壤交换性镁含量空间分布特征及影响因素,利用三明市4833个土壤样品化验结果和相关数据,应用统计学和地理统计学方法对三明市耕地土壤交换性镁含量空间分布特征及其影响因素进行了分析。结果表明：三明市耕地土壤样品交换性镁平均含量处于中等水平,平均值为62.7 mg/kg。各县（市、区）土壤样品交换性镁含量差异达到显著水平。土壤交换性镁含量缺乏样品所占比例和缺乏的耕地面积所占比例分别达到44.15%和40.72%。除土壤类型分组的交换性镁含量差异不显著外,其他影响因素分组的土壤交换性镁含量差异均显著。因此施用有机肥、调节土壤pH、客土改良砂土类质地土壤、补充镁肥等可以作为调控三明市土壤交换性镁含量的主要措施。
|
| [18] |
陈强, 杨劲松, 姚荣江, 等. 河套平原典型县域耕地土壤养分空间变异特征研究[J]. 中国农学通报, 2020, 36(10):102-108.
为研究河套平原地区土壤养分状况,探究县域尺度土壤养分的空间变异性与分布特征,以杭锦后旗县域耕地为对象,采用实地调查和地统计学相结合的方法,研究该区耕层土壤有机质、全氮、有效磷、速效钾养分含量状况,分析土壤养分空间分布及其变异特征。结果表明,研究区耕层土壤有机质、全氮、速效磷和速效钾平均含量依次为13.90 g/kg、0.93 g/kg、16.25 mg/kg和170.79 mg/kg,土壤有机质和全氮处于缺乏状态,速效磷含量中等,而速效钾的含量较为丰富。从空间变异及相关性看,各土壤养分的变异程度都为中等变异,有机质和速效钾为中等相关性,全氮和有效磷均为弱相关性。从空间分布格局看,土壤有机质与全氮分布规律相近,呈现西北部高、东南部和东北部低的特点;土壤速效钾整体分布较为均衡;速效磷空间差异明显,呈东高西低格局。针对该地区土壤养分含量及分布规律,在进行养分管理过程中,应以控氮、稳磷、稳钾、全区规划与重点应对相结合的原则,坚持以有机肥为主、有机与无机相结合的施肥策略。对县域尺度耕地土壤养分状况的调查研究和统计分析,明晰了区域土壤养分空间分布特征及其与土地利用方式之间的关系,为河套平原耕地土壤养分管理、地力培育以及合理施肥运筹等提供科学依据。
|
| [19] |
王淇韬, 孔祥斌, 辛芸娜, 等. 基于利用效率限制的县域耕地质量定级方法研究:以河北省平山县为例[J]. 北京师范大学学报(自然科学版), 2018, 54(3):292-299.
|
| [20] |
何东娟, 郭燕华, 张青璞, 等. 基于灌溉条件的作物生产潜力内插方法研究-以河北省为例[J]. 河北省科学院学报. 2019, 36(2):71-78.
|
| [21] |
李旭光, 邢凤丽, 袁雪, 等. 基于GIS的河北省太行山山前平原区耕地质量等级评价[J]. 山西农业科学, 2023, 51(7):793-801.
|
| [22] |
殷守强, 王鑫, 贺文龙. 耕地质量监测中不同抽样方法的精度比较-以河北省黄骅市为例[J]. 资源科学, 2016, 38(11):2049-2057.
耕地质量监测是及时掌握耕地质量和产能变化的必要手段,科学合理的抽样方法不仅能保证监测精度,还可以减少监测成本。本文以河北省黄骅市为研究区域,结合区域耕地质量的实际状况,运用6种抽样方法抽取耕地单元作为监测样点,并计算和比较不同抽样均值的最大相对误差和相对抽样精度。结果表明,抽样精度随着样本容量的增加而提升,速度由快变慢。与简单随机抽样相比,在抽样精度方面,系统分层抽样和系统等距抽样的优势极大,并随着样本数量的增加优势扩大;等别分层抽样的优势较大并且平稳;乡镇分层抽样的优势不大。网格分层抽样不但没有优势,反而可能造成精度的损失。不同的抽样方法在满足相同的监测精度时所需的样本量差异显著,合理的抽样方法在满足监测精度的同时,可以明显地减少样本量。本研究有助于选择出最合适的抽样方法,以期为耕地质量监测中不同抽样方法的精度比较提供方法指导。
|
| [23] |
吴克宁, 焦雪瑾, 梁思源, 等. 基于标准样地国家级汇总的耕地质量动态监测点构架研究[J]. 农业工程学报, 2008, 24(10):74-79.
|
| [24] |
余述琼, 张蚌蚌, 相慧, 等. 基于因素组合的耕地质量等级监测样点布控方法[J]. 农业工程学报2014, 30(24):288-297.
|
| [25] |
张玉臻, 孔祥斌, 刘炎, 等. 基于标准样地的省级耕地质量监测样地布设方法-以内蒙古自治区为例[J]. 资源科学, 2016, 38(11):2037-2048.
耕地质量是国家粮食安全的基础,科学布设耕地质量监测样地是支撑实时掌握国家耕地质量的关键技术环节。当前对耕地质量监测的研究主要集中在国家级和县级尺度,缺乏省级尺度。本文以内蒙古自治区为例,提出在省级标准样地基础上进行耕地监测样地数量和质量代表性两方面的优化,进而形成省级耕地质量监测布点的技术方法。本文运用最邻近点指数和面积代表性指数对标准样地的空间分布和数量进行分析和验证,结果表明,原144个标准样地在空间上呈均匀分布,但在数量上不能满足监测布样的要求;运用地统计学方法确定理论监测样地数量为194个;以“耕作制度二级区-土壤类型-地貌-土地利用现状”为基础,结合《中国土地资源图集》中的中国土地潜力分区图,将内蒙古划分为14个监测控制分区;依据“等别-监测分区”的控制顺序和优先选择标准样地的原则修正监测样地数量和确定其空间位置,最终确定省级监测样地数量为219个,其中标准样地124个,新增监测样地95个。依据本文提出的方法进行省级监测样地布设,一方面可以建立起省级尺度符合耕地质量空间变异特征的耕地质量监测布设样地体系,另一方面可以促进国家和县级耕地质量监测工作之间的衔接。
|
| [26] |
卫新东, 王宁, 王筛妮, 等. 省域尺度耕地质量监测样点布控-以陕西省为例[J]. 地球科学与环境学报, 2021, 43(3):549-558.
|
| [27] |
李玮, 杨建波, 温锦盼, 等. 基于GIS的耕地质量评价单元确定方法对比研究[J]. 中国农业信息, 2022, 34(5):11-23.
|
| [28] |
杨帆, 徐阳, 崔勇, 等. 近30年中国农田耕层土壤有机质含量变化[J]. 土壤学报, 2017, 54(5):1047-1056.
|
| [29] |
沈仁芳, 王超, 孙波. “藏粮于地、藏粮于技”战略实施中的土壤科学与技术问题[J]. 中国科学院院刊, 2018, 33(2):135-144.
|
| [30] |
徐明岗, 卢昌艾, 张文菊, 等. 我国耕地质量状况与提升对策[J]. 中国农业资源与区划, 2016, 37(7):8-14.
|
| [31] |
陈文广, 孔祥斌. 中国耕地保护科技创新战略探索[J]. 土壤通报, 2023, 54(4):947-954.
|
| [32] |
高君亮, 罗凤敏, 高永, 等. 农牧交错带不同土地利用类型土壤碳氮磷生态化学计量特征[J]. 生态学报, 2019, 39(15):5594-5602.
|
| [33] |
姚喜军, 徐进才, 刘静, 等. 伊金霍洛旗不同土地利用方式下的土壤养分特征[J]. 中国水土保持科学, 2017, 15(5):111-116.
|
| [34] |
张珍明, 林绍霞, 张清海, 等. 不同土地利用方式下草海高原湿地土壤碳、氮、磷分布特征[J]. 水土保持学报, 2013, 27(6):199-204.
|
| [35] |
杨静涵, 刘梦云, 张杰, 等. 黄土高原沟壑区小流域土壤养分空间变异特征及其影响因素[J]. 自然资源学报, 2020, 35(3):743-754.
了解表层土壤有机质与全量养分的空间变异规律及其影响因素,能够为黄土高原生态脆弱带的土壤养分管理提供参考。基于王东沟的93个野外采样点,综合多种地统计学方法,分析了土壤养分的空间变异特征及其影响因素。结果表明:研究区表层土壤养分含量处于中等水平,空间变异大小依次为全磷>全氮>有机质>全钾,均由结构性因素主导;有机质、全氮较全磷、全钾变异尺度小、空间自相关性弱、空间复杂程度高。有机质、全氮的分布格局为南高北低,主要影响因素包括海拔、坡度、曲率和土地利用;全磷呈现相反的北高南低分布,海拔、坡度、土地利用、土壤机械组成和人类活动对其影响较大;全钾则为内部高四周低,分布较破碎,坡向和土壤机械组成作用较强。据此建立了9个环境因子与土壤养分之间的回归预测模型,以期为研究区土壤资源的可持续利用提供数据支持。
|
| [36] |
姚冬辉, 吴建富, 卢志红, 等. 江西省耕地土壤pH值和养分状况调查[J]. 湖南农业大学学报(自然科学版), 2020, 46(5):574-579.
|
| [37] |
木合塔尔·艾买提, 肖鹏南, 周勇, 等. 山地丘陵交错地带耕层土壤pH和养分空间变异特征及影响因素分析[J]. 南方农业学报, 2019, 50(7):1432-1441.
|
| [38] |
张京华. 燕赵文化[M]. 沈阳: 辽宁教育出版社,1998.
|
| [39] |
李振旭, 郑涛. 邢台古地地理环境与商周文明[J]. 邢台学院学报, 2022, 37(2):5-9.
|
| [40] |
陈敬安, 殷超, 曾艳, 等. 湖泊有机质降解过程生源要素活化行为研究进展与展望[J]. 地球与环境, 2023, 51(06):655-666.
|
| [41] |
石建省, 李国敏, 梁杏, 等. 华北平原地下水演变机制与调控[J]. 地球学报, 2014, 35(5):527-534.
|
| [42] |
谢花林, 程玲娟. 地下水漏斗区农户冬小麦休耕意愿的影响因素及其生态补偿标准研究-以河北衡水为例[J]. 自然资源学报, 2017, 32(12):2012-2022.
近年来,国家大力提倡耕地休养生息,推行土地休耕制度,以提高地力,保障农业可持续发展。由于水资源短缺已经成为制约我国华北地区经济社会发展的重要因素。华北地下水漏斗区作为土地休耕的三个试点之一,主要目标是退耕高耗水的冬小麦,实现“一季休耕、一季雨养”,以减少地下水用量,修复地下水环境。农户作为农村最主要的生产与经营主体,如何给予农户合理的补偿,使其有效地实施冬小麦休耕政策成为缓解华北地区水资源短缺的关键因素。论文运用问卷调查法、机会成本法与Logistic模型对意愿与补偿标准进行分析,探讨了河北衡水地下水漏斗区农户冬小麦休耕意愿的主要影响因素及其生态补偿标准的合理性,以便为我国科学合理地制定地下水漏斗区耕地轮作休耕政策提供理论依据。研究结果表明：1）农户的文化程度、家庭农业劳动力总数、抚养比、农户对地下水漏斗区冬小麦休耕态度、农户对政府实行政策的信任度与农户休耕意愿呈显著的正向影响;2）务农工日、耕地质量、人均耕地面积与农户休耕意愿呈显著的负向影响;3）只考虑冬小麦对地下水影响,不考虑农产品市场价格变化的影响,建议河北衡水地区休耕农户补偿标准为518元/亩（7 770元/hm<sup>2</sup>）。最后论文提出了因地制宜地制定地下水漏斗区冬小麦休耕的补偿政策,发展绿色产业,转移剩余劳动力,并做好休耕政策宣传力度等政策建议。
|
| [43] |
王学, 李秀彬, 辛良杰, 等. 华北地下水超采区冬小麦退耕的生态补偿问题探讨[J]. 地理学报, 2016, 71(5):829-839.
压减冬小麦种植面积是有效减少华北平原地下水用量最为有效的方法。近期,国家和地方政府均出台了土地休耕政策,拟在华北地下水超采区退耕冬小麦,实现“一季休耕、一季雨养”。然而,华北平原已经出现农户主动退耕冬小麦的现象。在此背景下,明确冬小麦退耕的瞄准目标、制订合理适度的补偿标准是土地休耕政策有效实施所面对的核心问题。本文利用农户问卷数据,构建多层次logit模型分析农户种植制度决策的影响因素,并计算冬小麦退耕的机会成本,旨在确定华北平原地下水超采区土地休耕政策的瞄准目标,并为农户补偿标准的制订提供建议。结果表明：① 地块层次的土地质量和灌溉条件是解释农户种植制度差异的关键因素,近70%的四等耕地和90%以上的旱地已退出冬小麦耕作,土地休耕政策应瞄准一、二、三等土地质量的水浇地;② 相似耕作条件的地块种植春、夏玉米的净收益大致相同,土地休耕的机会成本等于种植冬小麦的净收益;③ 仅考虑冬小麦对地下水资源的影响,土地休耕政策初期以地下水回升和地下水环境恢复为主要目标时,建议以350元/亩作为补偿标准参考值,后期目标转为维持地下水资源采补平衡时,可调整至280元/亩左右。
|
| [44] |
赵霞, 黄瑞冬, 唐保军, 等. 潮土区不同土体构型对夏玉米生长与产量的影响[J]. 土壤通报, 2013, 44(3):538-542.
|
| [45] |
檀满枝, 李开丽, 史学正, 等. 华北平原土壤剖面质地构型对小麦产量的影响研究[J]. 土壤, 2014, 46(5):913-919.
|
| [46] |
赵执, 吴克宁, 王海. 我国农田整治中典型土壤障碍因素的机械改土技术[J]. 广东农业科学, 2012, 39(21):182-185+199.
|
| [47] |
韩迪, 李蓬勃, 徐东瑞. 河北省土壤有机质时空变化分析[J]. 土壤通报, 2016, 47(1):65-72.
|
| [48] |
李可, 杨勇, 刘亚军, 等. 基于回归克里格的河北省土壤有机质时空变化特征研究[J]. 中国土壤与肥料, 2020(3):1-7.
|
/
| 〈 |
|
〉 |
