中国稻田土壤有机质时空变化及其驱动因素

奥鹏网院作业

中国稻田土壤有机质时空变化及其驱动因素
李冬初1,2，黄晶2，马常宝3，薛彦东3，高菊生2，王伯仁2，张杨珠1，柳开楼2,4，韩天富2，张会民2

（1湖南农业大学资源环境学院，长沙 410128；2中国农业科学院农业资源与农业区划研究所/耕地培育技术国家工程实验室，北京 100081；3农业农村部耕地质量监测保护中心，北京 100125；4江西省红壤研究所/国家红壤改良工程技术研究中心，江西进贤 331717）

摘要：【目的】评价中国稻田土壤有机质时空变化特征，为提高耕地质量，应对种植结构调整及气候变化等提供支撑。【方法】基于1988—2017年开展的338个国家级定位监测点，分析稻田耕层土壤有机质含量变化特征、驱动因素以及对土壤容重影响。【结果】近30年全国稻田耕层土壤有机质平均提高了3.49 g·kg-1，年均增速0.09—0.12 g·kg-1。稻田耕层土壤有机质含量年均增速高低依次为东北、长江下游、长江中游、华南和西南。目前，全国稻田耕层土壤有机质平均含量32.4 g·kg-1，从高到低依次为长江中游、华南、东北、西南和长江下游。气候、土壤类型、氮肥投入以及种植制度等对土壤有机质产生影响。西南稻区和高纬度的东北稻区，稻田土壤有机质含量与年均气温显著负相关（P＜0.05），东部地区和低纬度的地区稻田土壤有机质含量与年均气温显著正相关（P＜0.05）。潜育型水稻土有机质平均含量显著高于其他类型水稻土。合适的氮肥投入量（200—300 kgN·hm-2·a-1）有利于土壤有机质累积。土壤容重及耕层深度与土壤有机质存在显著响应关系（P＜0.01）。【结论】我国稻田耕层土壤有机质含量整体呈上升趋势，年均增速呈现从南到北依次增加的趋势。年均气温和降水量，土壤类型，氮肥用量和种植制度等管理措施是区域土壤有机质含量分异的主要驱动因子。

关键词：稻田土壤；土壤有机质；国家级耕地监测点；驱动因子；氮肥投入量；中国

0 引言
【研究意义】土壤有机质与土壤质量密切相关，有机质数量的耗竭和质量的恶化可直接导致土壤生态功能的衰退[1]。土壤有机质的变化受多种因素的影响，具有较强的时空变异性[2-4]。水稻土是中国面积最大、分布最广的耕地土壤，占全国耕地面积的1/5，维持稻田土壤质量和生态功能对我国粮食生产和安全至关重要[5-6]。研究稻田土壤有机质变化对全面了解我国稻田合理施肥以及保障作物高产稳产具有重要指导意义。【前人研究进展】土壤有机质变化方面已开展了较多研究，包括国家尺度、省域尺度，以及典型地区和气候带，土壤类型区等方面。黄耀等基于1980—2000年文献数据的研究表明，全国近60%的农田表层土壤有机质含量上升，有机碳储量总体增加310—400 Tg，其中水稻土和潮土有机质含量增加明显[3]。杨帆等比较测土施肥和全国第二次土壤普查数据，结果表明近30年全国农田耕层土壤有机质提高4.85 g·kg-1，并指出我国农田有机质上升原因可能包括秸秆还田、农作物产量的快速提升、免耕少耕技术推广以及堆肥和绿肥的施用等[4]。TANG等通过14 371个野外监测点，估算出全国陆地生态系统不同土地利用方式碳库特征，认为其与年均气温和降水关系显著，并指出82.9%的碳库存在于0—1 m的土层中[7]。PAN等[8]通过比较国家耕地土壤监测数据和全国第二次土壤普查数据，研究表明全国稻田表层土壤有机碳密度显著高于旱地，其固碳潜力巨大层中。顾成军等分别在不同省域尺度上研究土壤有机质变化特征，认为不同土地利用方式，肥料的大量施用，秸秆还田及其土壤类型等显著影响土壤有机质含量[9-16]。杨玉盛等分别分析了中亚热带山区、黄土丘陵沟壑区、东北黑土区等地带性土壤有机质变化特征，表明不同土地利用方式、土壤类型、土壤覆盖、土壤坡度等影响土壤有机质含量[17-19]。任丽等报道了苹果、柑橘等果树种植区土壤有机质空间分布及其影响因素，并建立了相关模型[20-21]。此外还有较多关于不同点位上的土壤有机质变化研究，多数报道施用有机肥显著提升土壤有机质含量，常规施肥土壤有机质含量呈上升趋势[22-24]。在有机质变化的研究手段上有基于多点位的统计学比较分析[3-4, 6-10]，也有基于地统计学方法及GIS技术的时空变异特征分析[12,25]，以及模型分析[11, 26-27]。目前，研究较多集中在对土壤有机质驱动因子的探讨，以及土壤有机质结构、形态及功能等研究[28-31]。【本研究切入点】目前国家尺度的农田土壤有机质含量变化研究一方面基于布置在全国各地的长期定位试验点。比如中国农业科学院的肥力与肥料效应监测网长期定位试验[32]，中国科学院CERN台站长期定位试验，以及各地农业科学院系统布置的大量长期定位试验，其研究多基于不同施肥、轮作、耕作等对土壤有机质的影响[22]。另一方面，基于区域尺度上的采样数据或者文献调研数据，通过与第二次土壤普查数据比较，或者通过模型模拟获取全国农田土壤有机质含量变化[3, 8]。农业农村部自1988年开始建立国家级耕地质量监测点，其中在中国主要水稻种植区稻田土壤上布置了大量的长期定位监测点。李建军等依托国家耕地质量1988—2012年数据研究了中国稻田土壤基础地力变化特征和长江中下游稻田土壤养分变化特征，结果表明中国稻田土壤基础地力总体呈上升趋势，基础地力与水稻产量同步上升。长江中下游稻田土壤养分含量基本呈上升趋势，土壤肥力总体改善[33-34]。韩天富等依托国家耕地质量1988—2017年数据研究认为施肥尤其是有机肥与无机肥配合施用是维持和提高水稻高产的重要措施[35]。近年来，我国农业结构调整，尤其是种植结构优化和有机替代，化肥减施增效等政策密集影响，土壤面临更加复杂的作物、肥料投入和气候变化。另一方面，基于国家级耕地质量监测点的全国稻田土壤有机质长期、动态、定位等的系统研究和分析相对较少，其变化特征、速率和趋势不明确。【拟解决的关键问题】本研究拟依托338个国家级稻田土壤长期定位监测点，系统分析中国主要水稻种植区稻田土壤有机质变化特征，旨在揭示我国稻田土壤有机质演变规律及区域分异特征，分析稻田土壤有机质变化驱动因素。评价土壤有机质对土壤容重和耕层深度等影响，为稻田土壤合理施肥和作物稳产高效提供科学指导，为耕地地力提升和建设高质量稻田提供科学依据。

1 材料与方法
1.1 监测点分布概况
基于农业农村部在全国不同稻作区布置的土壤监测点及其数据库[35]，选取全国水稻种植区域稻田土壤监测点共338个，根据水稻监测点分省（区）分布情况，结合自然地理区域，本研究监测点位按照数量高低，分布于长江中游（111个），长江下游（82个），华南（80个），西南（47个），东北（18个）。在全国稻田总的338个土壤监测点中，有140个为长期定位监测点，其中长江中游、长江下游、华南、西南和东北分别有52个、29个、31个、22个和6个。

全国稻田土壤长期监测点自1988年开始监测，分别设置无肥对照区和常规施肥区（农民习惯施肥）两个处理。种植制度根据不同区域气候特点分别包括①东北，一年一熟制，种植水稻为主；②长江中游，一年两熟制，种植制度主要为双季稻，部分为水稻-油菜（小麦、烟叶、蔬菜、绿肥等）轮作；③长江下游，一年两熟制，以水稻-小麦（油菜、绿肥等）水旱轮作为主；④西南，一年两熟制，以水稻-小麦（油菜、玉米、冬闲等）水旱轮作为主；⑤华南，一年两熟制（部分三熟制），种植制度主要为双季稻，部分为水稻-蔬菜（烟叶、绿肥等）轮作。各点分别监测肥料养分氮磷钾投入、主要作物产量以及土壤肥力状况。

1.2 样品采集与分析
各监测点按照统一的规范，分别于每年度最后一季作物收获后采集耕层土壤样品。耕层指农业耕作、施肥、灌溉影响区域和作物根系分布的集中层段，其颜色、结构、紧实度都有明显的特征和界线。每个处理采用“S”型采集5—10钻土壤样品，3次重复。土壤样品自然风干，人工除去植物残体和石块后过2 mm 筛，混匀后备用。土壤各指标测定方法依据《土壤分析技术规范》[36]。各监测点每年记录施肥种类、数量，作物籽粒和秸秆产量以及管理措施等，现场根据作物根系分布和土壤颜色、紧实度等指标调查耕层厚度。

1.3 数据处理与统计分析
各稻作区监测点位数量及其代表面积不同，本文利用各省（区）稻田土壤有机质含量的平均值与该省稻田土壤面积比例加权平均后得到各稻作区土壤有机质平均含量[37-38]。

pagenumber_ebook=98,pagenumber_book=2412
式中，SOM为稻作区土壤有机质平均含量（g·kg-1）；Mi为该区域i省份土壤有机质平均含量（g·kg-1），Ai为i省份稻作土壤面积（104 hm2）。

各点数据采用Excel 2016整理计算，运用SigmaPlot 12.0分析并作图，IBM SPSS 22.0进行差异显著性检验（LSD，Duncan (D)，P＜0.05）。

2 结果
2.1 土壤有机质含量空间变化
全国稻田耕层土壤有机质平均含量32.4 g·kg-1，变幅在11.3—65.0 g·kg-1（表1）。其中湖南、云南、贵州、广西、广东、江西、黑龙江、吉林、浙江和福建等10省（区）稻田耕层土壤有机质平均含量在30.0 g·kg-1以上。安徽、海南和上海低于25.0 g·kg-1。各稻作区域耕层土壤有机质平均含量高低顺序为长江中游35.1 g·kg-1、华南34.2 g·kg-1、东北33.6 g·kg-1、西南31.0 g·kg-1和长江下游26.3 g·kg-1。由图1可知，根据箱线图统计的各稻作区耕层土壤有机质含量置信区间分别为长江中游20.6—51.2 g·kg-1，华南15.2—48.6 g·kg-1，东北11.9—51.6 g·kg-1，西南18.1—45.4 g·kg-1和长江下游17.2—43.9 g·kg-1。各稻作区稻田耕层土壤有机质含量差异显著。长江中游地区耕层土壤有机质含量显著高于西南地区和长江下游地区（P＜0.05）。从经度分布来看，各区域稻田耕层土壤有机质含量无显著差异。从纬度分布看，N25°—30°和N20°—25°区域稻田耕层土壤有机质含量显著高于＜N20°和N30°—40°区域（P＜0.05）。

表1 不同区域稻田耕层土壤有机质平均含量
Table 1 The mean SOM content of paddy soil in plough layer in different regions

pagenumber_ebook=99,pagenumber_book=2413
NE: Northeast of China; MYR: Middle of the Yangtze River; LYR: Lower of Yangtze River; SW: Southwest of China; SC: South of China

2.2 土壤有机质含量的时间变化
2.2.1 全国稻田土壤有机质含量的时间变化 全国140个长期监测点结果表明，常规施肥条件下近30年全国稻田耕层土壤有机质含量总体呈现上升趋势（图2-a）。通过对所有监测点338个点位数据土壤有机质含量的统计，与1988年相比，近30年全国稻田耕层土壤有机质平均含量上升3.49 g·kg-1（图2-b）。回归方程表明稻田耕层土壤有机质年增加速率0.09—0.12 g·kg-1（P＜0.01）。

2.2.2 不同区域稻田土壤有机质含量的时间变化 不同稻作区域长期监测点稻田耕层土壤有机质含量变化见图3。各稻区稻田耕层土壤有机质含量均随耕作时间呈现上升趋势。土壤有机质年增长速率从高到低依次为东北（0.31 g·kg-1·a-1，P＜0.01）、长江下游（0.22 g·kg-1·a-1，P＜0.01）、长江中游（0.19 g·kg-1·a-1，P＜0.05）、华南（0.16 g·kg-1·a-1，P＜0.05）和西南（0.11 g·kg-1·a-1，P＜0.05）。土壤有机质含量年均增速总体表现为从北到南依次降低的趋势。

pagenumber_ebook=100,pagenumber_book=2414
图1 不同区域稻田耕层土壤有机质含量及随经度和纬度变化
Fig.1 The SOM content of paddy soil in different rice cultivation regions and different regions with longitude and latitude

箱式图内中间实线、虚线以及箱体的下边缘线、上边缘线、下误差线、上误差线分别表示中值、平均值、第25、75、5、95 百分位数。不同大写字母表示不同区域差异显著（P＜0.05）。下同
Solid line and dotted line, lower edge and upper edge, bars represent median and mean values, 25th and 75th, 5th and 95th percentiles of all data, respectively.Capital letter indicates the difference significance among different regions at the 0.05 level.The same as below

pagenumber_ebook=100,pagenumber_book=2414
图2 常规施肥下稻田土壤有机质变化(a，长期监测点；b，所有监测点1988、1998、2008、2017年统计)
Fig.2 Change of SOM content in paddy soil under conventional fertilization within period of 1988 to 2017(a, data comes from long term monitoring points; b, data comes from all monitoring points in 1988, 1998, 2008, and 2017)

pagenumber_ebook=101,pagenumber_book=2415
图3 不同稻作区域稻田耕层土壤有机质含量变化
Fig.3 Change of SOM content in paddy soil in different region within period of 1988 to 2017

拟合曲线上下虚线为95%置信区间
The dotted lines at fitted curve up and down is the 95% confidence interval curve fitting

表2 不同区域长期监测点土壤有机质含量变化统计
Table 2 Variation of SOM in monitoring points in different region

pagenumber_ebook=101,pagenumber_book=2415
2.2.3 不同区域有机质含量变化的点位特征 通过对不同区域长期监测点有机质含量的定位跟踪，1988—2017年间全国140个长期监测点中，58个监测点稻田土壤有机质含量显著上升，占总量的41.4%（表2）。44个监测点稻田土壤有机质含量显著下降，占总量的31.4%，另外有27.1%的长期监测点稻田土壤有机质含量无显著上升或下降趋势。受到长期监测点地理位置、气候、种植制度、肥料投入等影响，1988—2017年间长江中游地区52个长期定位监测点中，48.1%的监测点土壤有机质含量上升，42.3%的监测点下降，9.1%的监测点无显著变化。华南地区31个长期定位监测点中32.3%监测点的有机质含量上升，25.8%下降，41.9%无显著变化。长江下游地区29个长期定位监测点中41.4%上升，24.1%下降，34.5%无显著变化。西南地区22个长期定位监测点中36.4%上升，22.7%下降，40.9%无显著变化。东北地区6个长期定位监测点中50.0%上升，33.3%下降，16.7%无显著变化。

2.3 土壤有机质含量变化影响因素
土壤有机质含量受土壤类型、气候特征、种植制度以及施肥等的影响。不同区域土壤有机质含量对气温和降水量响应关系不同（图4）。从东到西随经度变化以及从南到北随纬度变化，土壤有机质含量与年均温度相关系数呈下降趋势，两者由正相关逐渐变化为显著负相关（P＜0.05）。从东到西随经度变化，土壤有机质含量与年均降水量在＜E105°、E110°—115°、＞125°区间表现为负相关，在E105°—110°、E115°—120°、E120°—125°区域为正相关（P＜0.05），其相关系数整体表现为两头低，中间高。从南到北随纬度变化，土壤有机质含量与年均降水量由正相关关系逐渐变化为负相关（P＜0.05），相关系数随纬度有下降趋势。

pagenumber_ebook=102,pagenumber_book=2416
图4 稻田土壤有机质含量与年均温度和降水量相关关系随经度和纬度变化
Fig.4 The relationships between the SOM content and annual average air temperature, precipitation following change with longitude and latitude

从现有监测点位水稻土类型来看（图5-a），潜育型水稻土的稻田耕层土壤有机质平均含量为38.9 g·kg-1，显著高于其他类型水稻土（P＜0.05）。耕层土壤有机质含量受种植制度影响相对较少，一年三熟制稻田耕层土壤有机质平均含量要略高于一年两熟和一年一熟制稻田（图5-b）。土壤有机质与氮肥投入响应关系表明，在投入量0—200 kg·hm-2·a-1下，土壤有机质与氮肥投入呈负相关（图5-c，P＜0.05）。年投入量200—300 kg·hm-2·a-1下，土壤有机质与氮肥投入呈极显著正相关（P＜0.01）。年投入量大于300 kg·hm-2·a-1下，土壤有机质与氮肥年投入量极显著负相关（P＜0.01）。

2.4 土壤有机质变化与土壤容重和耕层深度关系
土壤有机质含量与土壤容重及耕层深度存在响应关系。通常情况下有机质含量高，其土壤疏松，结构好，容重低，耕层深厚。通过对我国稻田耕层土壤602组有机质含量与容重数据做线性回归获得容重与土壤有机质含量的经验方程：y（容重g·cm-3）=-0.0049×SOM(g·kg-1) + 1.3739 （R2 = 0.1198, n = 602）（图6-a）。同样，对670组耕层深度数据与土壤有机质含量获得经验方程：y（SOM g·kg-1）=0.2879×T（耕层厚度 cm）+ 25.3651（R2 = 0.0120, n = 670）（图6-b）。

3 讨论
3.1 土壤有机质时空变化特征
20世纪80年代以来，随着我国农业投入的增加和农业科技快速发展，水稻总产量持续上升，农作物秸秆资源数量和还田比例持续增加，我国稻田耕层土壤有机质含量整体呈现上升趋势。本研究得出近30年全国稻田耕层有机质含量平均增加约3.49 g·kg-1，年均增速0.09—0.12 g·kg-1。与前人观点基本相似[3-4]。从区域来看，受到各地水热条件、种植模式和肥料投入等影响，稻田土壤有机质含量上升速率存在差异，年均递增速率呈现从南到北依次增加趋势，结果与李建军等的观点相似[39]。我国稻田土壤也存在提升和发展不平衡状况，在全国耕地质量长期监测点位中，约2/3以上的点位稻田土壤有机质含量保持稳定或呈增加趋势，其结果略低于黄耀等文献调研结果[3]。受时间跨度、点位和数据来源等影响，以往研究大多基于文献调研或模型模拟研究国家尺度农田土壤有机质动态变化[2-4]。本研究依托国家级耕地质量长期定位监测点，具有长期稳定，定位精准，能真实反应农业生产方式和投入变化情况，数据可靠性高等特点。

pagenumber_ebook=103,pagenumber_book=2417
图5 不同水稻土类型（a）、种植制度（b）以及氮肥投入量（c）稻田的耕层土壤有机质含量变化
Fig.5 The change of SOM content with different paddy soil types (a), different cropping systems (b) and N fertilizer applications (c)

pagenumber_ebook=103,pagenumber_book=2417
图6 稻田土壤有机质与土壤容重（a）及耕层深度（b）的关系
Fig.6 Correlations between bulk density and SOM content (a), SOM and plough depth (b) content for the plough layers of the paddy soil

我国稻田土壤有机质提升技术措施主要包括农作物秸秆还田、施用有机肥或化肥有机肥配施以及种植绿肥[40-42]。稻田耕层土壤有机质含量明显高于全国耕地土壤以及旱地耕作土壤[4]，其中10个省份超过30.0 g·kg-1。在不同区域，稻田土壤有机质含量呈现明显的时空变异性[2-4]。长江中游和华南为我国典型双季稻产区，水稻生产水平普遍较高，具有较好的水热条件，有利于秸秆腐解，其稻田土壤有机质含量整体较高。长江中游、华南和东北稻田耕层土壤有机质含量较西南和长江下游高，原因可能是与后两地为我国典型水旱轮作区域，稻田水旱交替加速土壤有机质分解有关[43-44]。

3.2 土壤有机质变化因素分析
稻田土壤有机质主要来源作物根茬、秸秆、绿肥还田，有机物料投入等。土壤有机质含量在固持和分解中维持动态平衡，影响稻田土壤有机质含量的因素有气候、土壤类型、地形、施肥和耕作措施等[2, 45-47]。温度和水分是决定土壤有机质输入和分解的气候因子，其一方面影响作物产量，制约土壤有机质输入量；另一方面对土壤水热状况和微生物活动产生深远影响[48]。本研究表明，高纬度地区，稻田土壤有机质含量与年均温度显著负相关，说明温度越高，有机质分解越快。低纬度地区，稻田土壤有机质含量与年均温度呈正相关，温度越高，越有利于作物生长，累积更多的碳投入[48-49]。从东到西，稻田土壤有机质含量与年均温度由正相关，转变为负相关，长江下游、华南地区东部年均气温对土壤有机质积累起促进作用，而到西南地区，年均气温则对土壤有机质分解起主导作用。从稻田土壤有机质与年降水量关系来看，在＜E105°、E110°—115°、＞E125°区间，降水加剧了有机质分解，两者表现为负相关，在E105°—110°、E115°—120°、E120°—125°区域，随着降水量增加，稻田淹水时间和淹水量增加，产生还原环境，有利于缓解土壤有机质矿化分解，增加土壤有机质积累[46]。随着纬度增加，稻田土壤有机质与年降水量之间相关系数有下降趋势。总体而言，温度和水分二者的综合配置影响稻田土壤有机质含量的地带性分布。

施肥是影响土壤有机质至关重要的因素。长期施肥显著提高土壤有机质含量[50-52]。全国耕地质量监测数据揭示稻田土壤有机质与肥料投入存在响应关系。合适的氮肥年投入量（200—300 kgN·hm-2·a-1）能提升土壤有机质含量，较低和过高氮肥年投入量均不利于土壤有机质含量提升。目前南方的长江下游、长江中游和华南等区域稻田氮肥年投入量普遍较高。因此，减肥是必然选择，土壤有机质含量提升要与化肥减施有机结合，协调发展。就我国现有肥料投入水平来看，生产上应该以“减氮”为原则，采用秸秆还田和有机肥替代部分化肥，达到减施、提质和增效目的。水稻土类型和种植制度等主要通过土壤水热状况、微生物种群和外源投入等对稻田土壤有机质含量产生影响。

3.3 有机质与土壤容重及耕层深度的关系
土壤有机质对土壤质量及功能的调节起关键作用，较低的土壤有机质会降低土壤结构的稳定性。研究表明，土壤有机质提升有利于增加土壤孔隙，降低土壤容重，增厚土壤耕层 [53-55]。WU和CALLESEN等[55-56]分别利用中国矿质土壤和加拿大森林土壤拟合土壤容重与有机质含量。本研究得出稻田土壤有机质含量与土壤容重响应关系，与PAN等[8]结果基本一致。通过拟合的经验方程，可用来补充缺失土壤容重值，为计算稻田土壤碳库储量提供支撑。

土壤有机质有利于构建肥沃土壤耕层。相反，合适土壤耕层深度有利于土壤微生物活动，影响土壤有机质累积和分解。本研究结果也表明，稻田土壤有机质含量与耕层深度呈正相关关系（P＜0.01），但两者相关性不如土壤容重，耕层深度可能与稻田管理措施，比如土壤翻耕深度、次数、以及机械化耕作管理等有关[56]。

4 结论
近30年全国稻田耕层土壤有机质含量平均上升3.49 g·kg-1，从南到北土壤有机质含量年均增速依次增加。不同区域稻田土壤有机质对年均气温和降水量响应关系不一样。从东部到西部以及纬度从低到高，土壤有机质含量与年均气温相关性由正相关转变为负相关。氮肥年投入量、水稻土类型以及种植制度等对稻田土壤有机质产生影响。稻田土壤有机质含量与土壤容重及耕层深度有响应关系，提升土壤有机质有利于降低土壤容重，土壤有机质含量随土壤耕层增加而呈上升趋势。

References

[1] LARSON W E, PIERCE F J.The dynamics of soil quality as a measure of sustainable management// DORAN J W, COLEMAN D C,BEZEDICK D F, STEWART D A.Defining Soil Qualit for a Sustainable Environment.Madison, Wisconsin: SSSA Special Publication, 1994: 37-51.

[2] 周莉, 李保国, 周广胜.土壤有机碳的主导影响因子及其研究进展.地球科学进展, 2005, 20(1): 99-105.ZHOU L, LI B G, ZHOU G S.Advances in controling factors of soil organic carbon.Advances in Earth Science, 2005, 20(1): 99-105.(in Chinese)

[3] 黄耀, 孙文娟.近20 年来中国大陆农田表土有机碳含量的变化趋势.科学通报, 2006, 51(7): 750-763.HUANG Y, SUN W J.The change trend of organic carbon content in farmland topsoil in mainland China in recent 20 years.Chinese Science Bulletin, 2006, 51(7): 750-763.(in Chinese)

[4] 杨帆, 徐洋, 崔勇, 孟远夺, 董燕, 李荣, 马义兵.近30 年中国农田耕层土壤有机质含量变化.土壤学报, 2017, 54(5): 1047-1056.YANG F, XU Y, CUI Y, MENG Y D, DONG Y, LI R, MA Y B.Variation of soil organic matter content in croplands of China over the last three decades.Acta Pedologica Sinica, 2017, 54(5): 1047-1056.(in Chinese)

[5] 龚子同.中国土壤分类.北京: 科学出版社, 2003.GONG Z T.Chinese Soil Classification.Beijing: Science Press, 2003.(in Chinese)

[6] 武红亮, 王士超, 闫志浩, 槐圣昌, 马常宝, 薛彦东, 徐明岗, 卢昌艾.近30年我国典型水稻土肥力演变特征.植物营养与肥料学报,2018, 24(6): 1416-1424.WU H L, WANG S C, YAN Z H, HUAI S C, MA C B, XUE Y D,XU M G, LU C A.Evolution characteristics of fertility of typical paddy soil in China in recent 30 years.Journal of Plant Nutrition and Fertilizers, 2018, 24(6): 1416-1424.(in Chinese)

[7] TANG X L, ZHAO X, BAI Y F, TANG Z Y, WANG W T, ZHAO Y C, WAN H W, XIE Z Q, SHI X Z, WU B F, WANG G X, YAN J H,MA K P, DU S, LI S G, HAN S J, MA Y X, HU H F, HE N P, YANG Y H, HAN W X, HE H L, YU G R, FANG J Y, ZHOU G Y.Carbon pools in China’s terrestrial ecosystems: New estimates based on an intensive field survey.Proceedings of the National Academy of Sciences, 2018, 115: 4021-4026.

[8] PAN G X, LI L Q, WU L S, ZHANG X H.Storage and sequestration potential of topsoil organic carbon in China’s paddy soils.Global Change Biology, 2004, 10(1): 79-92.

[9] 顾成军, 史学正, 于东升, 徐胜祥, 孙维侠, 赵永存.省域土壤有机碳空间分布的主控因子——土壤类型与土地利用比较.土壤学报, 2013, 50(3): 425-432.GU C J, SHI X Z, YU D S, XU S X, SUN W X, ZHAO Y C.Main factor controlling SOC spatial distribution at the province scale as affected by soil type and land use.Acta Pedologica Sinica, 2013,50(3): 425-432.(in Chinese)

[10] 赵明松, 张甘霖, 吴运金, 李德成, 赵玉国.江苏省土壤有机质含量时空变异特征及驱动力研究.土壤学报, 2014, 51(3): 448-458.ZHAO M S, ZHANG G L, WU Y J, LI D C, ZHAO Y G.Temporal and spatial variability of soil organic matter and its driving force in Jiangsu Province, China.Acta Pedologica Sinica, 2014, 51(3):448-458.(in Chinese)

[11] 姜赛平, 张怀志, 张认连, 李兆君, 谢良商, 徐爱国.基于三种空间预测模型的海南岛土壤有机质空间分布研究.土壤学报, 2018,55(4): 1007-1017.JIANG S P, ZHANG H Z, ZHANG R L, LI Z J, XIE L S, XU A G.Research on spatial distribution of soil organic matter in Hainan island based on three spatial prediction models.Acta Pedologica Sinica,2018, 55(4): 1007-1017.(in Chinese)

[12] 姜赛平, 张认连, 张维理, 徐爱国, 张怀志, 谢良商, 冀宏杰.近30年海南岛土壤有机质时空变异特征及成因分析.中国农业科学,2019, 52(6): 1032-1044.JIANG S P, ZHANG R L, ZHANG W L, XU A G, ZHANG H Z, XIE L S, JI H J.Spatial and temporal variation of soil organic matter and cause analysis in Hainan Island in resent 30 years.Scientia Acricultura Sinica, 2019, 52(6): 1032-1044.(in Chinese)

[13] 胡克林, 余艳, 张凤荣, 王茹.北京郊区土壤有机质含量的时空变异及其影响因素.中国农业科学, 2006, 39(4): 764-771.HU K L, YU Y, ZHANG F R, WANG R.The spatial-temporal variability of soil organic matter and its influencing factors in suburban area of Beijing.Scientia Agricultura Sinica, 2006, 39(4):764-771.(in Chinese)

[14] 史利江, 郑丽波, 梅雪英, 俞立中, 贾正长.上海市不同土地利用方式下的土壤碳氮特征.应用生态学报, 2010, 21(9): 2279-2287.SHI L J, ZHENG L B, MEI X Y, YU L Z, JIA Z C.Characteristics of soil organic carbon and total nitrogen under different land use types in Shanghai.Chinese Journal of Applied Ecology, 2010, 21(9): 2279-2287.(in Chinese)

[15] 曾招兵, 汤建东, 刘一峰, 张满红, 林碧姗.广东耕地土壤有机质的变化趋势及其驱动力分析.土壤, 2013, 45(1): 84-90.ZENG Z B, TANG J D, LIU Y F, ZHANG M H, LIN B S.Changes and driving forces of farmland organic matter in Guangdong Province,China.Soils, 2013, 45(1): 84-90.(in Chinese)

[16] 周睿, 潘贤章, 王昌坤, 刘娅, 李燕丽, 石荣杰, 解宪丽.上海市城郊土壤有机质的时空变异特征及其影响因素.土壤, 2014, 46(3):433-438.ZHOU R, PAN X Z, WANG C K, LIU Y, LI Y L, SHI R J, XIE X L.Spatial-temporal variation characteristics of soil organic matter and its impact factors in suburban Shanghai.Soils, 2014, 46(3): 433-438.(in Chinese)

[17] 杨玉盛, 谢锦升, 盛浩, 陈光水, 李旭.中亚热带山区土地利用变化对土壤有机碳储量和质量的影响.地理学报, 2007, 62(11):1123-1131.YANG Y S, XIE J S, SHENG H, CHEN G S, LI X.The impact of land use/cover change on soil organic carbon stocks and quality in mid-subtropical mountainous area of southern China.Acta Geographica Sinica, 2007, 62(11): 1123-1131.(in Chinese)

[18] 连纲, 郭旭东, 傅伯杰, 虎陈霞.黄土丘陵沟壑区县域土壤有机质空间分布特征及预测.地理科学进展, 2006, 25(2): 112-122.LIAN G, GUO X D, FU B J, HU C X.Spatial variability and prediction of soil organic matter at county scale on the Loess Plateau.Progress in Geography, 2006, 25(2): 112-122.(in Chinese)

[19] YANG O, ALAIN N R, WANG L X, YAN B X.Spatio-temporal patterns of soil organic carbon and pH in relation to environmental factors-A case study of the Black Soil Region of Northeastern China.Agriculture, Ecosystems and Environment, 2017, 245: 22-31.

[20] 任丽, 杨联安, 王辉, 杨粉莉, 陈卫军, 张林森, 徐瑾昊.基于随机森林的苹果区土壤有机质空间预测.干旱区资源与环境, 2018,32(8): 141-146.REN L, YANG L A, WANG H, YANG F L, CHEN W J, ZHANG L S, XU J H.Spatial prediction of soil organic matter in apple region based on random forest.Journal of Arid Land Resources and Environment, 2018, 32(8): 141-146.(in Chinese)

[21] 段丽君, 张海涛, 郭龙, 杜佩颖, 陈可, 琚清兰.典型柑橘种植区土壤有机质空间分布与含量预测.华中农业大学学报, 2019, 38(1):73-81.DUAN L J, ZHANG H T, GUO L, DU P Y, CHEN K, JU Q L.Spatial distribution and content prediction of soil organic matter in typical citrus growing area.Journal of Huazhong Agricultural University, 2019, 38(1): 73-81.(in Chinese)

[22] 徐明岗, 张文菊, 黄绍敏.中国土壤肥力演变.北京: 中国农业科学技术出版社, 2015.XU M G, ZHANG W J, HUANG S M. Soil Fertility Evolution in China.Beijing: China Agricultural Science and Technology Press,2015.(in Chinese)

[23] 张淑香, 张文菊, 沈仁芳, 徐明岗.我国典型农田长期施肥土壤肥力变化与研究展望.植物营养与肥料学报, 2015, 21(6): 1389-1393.ZHANG S X, ZHANG W J, SHEN R F, XU M G.Variation of soil quality in typical farmlands in China under long-term fertilization and research expedition.Journal of Plant Nutrition and Fertilizer, 2015,21(6): 1389-1393.(in Chinese)

[24] ZHANG W J, XU M G, WANG X J, HUANG Q H, NIE J, LI Z Z, LI S L., HWANG S W, LEE K B.Effects of organic amendments on soil carbon sequestration in paddy fields of subtropical China.Journal of Soils and Sediments, 2012, 12: 457-470.

[25] CHEN D, CHANG N J, XIAO J F, ZHOU Q B, WU W B.Mapping dynamics of soil organic matter in croplands with MODIS data and machine learning algorithms.Science of the Total Environment, 2019,669: 844-855.

[26] BEGUM K, KUHNERT M, YELURIPATI J, OGLE S, PARTON W,KADER M A, SMITH P.Soil organic carbon sequestration and mitigation potential in a rice cropland in Bangladesh – a modelling approach.Field Crops Research, 2018, 226: 16-27.

[27] 段丽君, 郭龙, 张海涛, 琚清兰.基于改进OK模型的土壤有机质空间分布预测——以宜都市红花套镇为例.中国生态农业学报,2019, 27(1): 131-141.DUAN L J, GUO L, ZHANG H T, JU Q L.Prediction of spatial distribution of soil organic matter based on improved OKmodels: A case study of Honghuatao Town in Yidu City.Chinese Journal of Eco-Agriculture, 2019, 27(1): 131-141.(in Chinese)

[28] PAN G X, SMITH P, PAN W N.The role of soil organic matter in maintaining the productivity and yield stability of cereals in China.Agriculture, Ecosystems and Environment, 2009, 129: 344-348.

[29] LI M F, WANG J, GUO D, YANG R R, FU H.Effect of land management practices on the concentration of dissolved organic matter in soil: A meta-analysis.Geoderma, 2019, 344: 74-81.

[30] NATH A J, BRAHMA B, SILESHI G W, DAS A K.Impact of land use changes on the storage of soil organic carbon in active and recalcitrant pools in a humid tropical region of India.Science of The Total Environment, 2018, 624: 908-917.

[31] WIESMEIER M, URBANSKI L, HOBLEY E, LANG B, VON LÜTZOW M, MARIN-SPIOTTA E, VAN WESEMAEL B, RABOT E, LIEß M, GARCIA-FRANCO N, WOLLSCHLÄGER U, VOGEL H J, KÖGEL-KNABNER I.Soil organic carbon storage as a key function of soils - A review of drivers and indicators at various scales.Geoderma, 2019, 333: 149-162.

[32] 赵秉强, 张夫道.我国的长期肥料定位试验研究.植物营养与肥料学报, 2002, 8(增刊): 3-8.ZHAO B Q, ZHANG F D.Long-term fertilizer experiments in China.Plant Nutrition and Fertilizer Science, 2002, 8(Suppl.): 3-8.(in Chinese)

[33] 李建军, 徐明岗, 辛景树, 段建军, 任意, 李冬初, 黄晶, 申华平,张会民.中国稻田土壤基础地力的时空演变特征.中国农业科学,2016, 49(8): 1510-1519.LI J J, XU M G, XIN J S, DUAN J J, REN Y, LI D C, HUANG J,SHEN H P, ZHANG H M.Spatial and Temporal Characteristics of Basic Soil Productivity in China.Scientia Agricultura Sinica, 2016,49(8): 1510-1519.(in Chinese)

[34] 李建军, 辛景树, 张会民, 段建军, 任意, 孙楠, 徐明岗.长江中下游粮食主产区25年来稻田土壤养分演变特征.植物营养与肥料学报, 2015, 21(1): 92-103.LI J J，XIN J S，ZHANG H M，DUAN J J，REN Y，SUN N，XU M G.Evolution characteristics of soil nutrients in the main rice production regions，the middle-lower reach of Yangtze River of China.Journal of Plant Nutrition and Fertilizer, 2015, 21(1): 92-103.(in Chinese)

[35] 韩天富, 马常宝, 黄晶, 柳开楼, 薛彦东, 李冬初, 刘立生, 张璐,刘淑军, 张会民.基于Meta分析中国水稻产量对施肥的响应特征.中国农业科学, 2019, 52(11): 1918-1929.HAN T F, MA C B, HUANG J, LIU K L, XUE Y D, LI D C, LIU L S,ZHANG L, LIU S J, ZHANG H M.Variation in rice yield response to fertilization in China: Meta-analysis.Scientia Agricultura Sinica,2019, 52(11): 1918-1929.(in Chinese)

[36] 全国农业技术推广服务中心.土壤分析技术规范.北京: 中国农业科学技术出版社, 2006.The Center of Extending and Service of Agricultural Technique in China.Soil Analysis Technical Specifications.Beijing: China Agricultural Science and Technology Press, 2006.(in Chinese)

[37] 沈善敏.中国土壤肥力.北京: 中国农业出版社, 1998.SHEN S M. Soil Fertility in China. Beijing: China Agriculture Press,1998.(in Chinese)

[38] 中华人民共和国国家统计局.2018年中国统计年鉴.北京: 中国统计出版社, 2019.China National Bureau of Statistics.China Statistical Yearbook, 2018.Beijing: China Statistics Press, 2019.(in Chinese)

[39] 李建军.我国粮食主产区稻田土壤肥力及基础地力的时空演变特征[D].贵阳: 贵州大学, 2015.LI J J.Temporal and spatial evolution characteristics of paddy soil fertility and basic soil fertility in main grain producing areas of China[D].Guiyang: Guizhou University, 2015.(in Chinese)

[40] GAO S J, GAO J S, CAO W D, ZOU C Q, HUANG J, BAI J S, DOU F G.Effects of long-term green manure application on the content and structure of dissolved organic matter in red paddy soil.Journal of Integrative Agriculture, 2018, 17(5): 60345-60347.

[41] LI Z Q, ZHAO B Z, ZHANG J B.Effects of maize residue quality and soil water content on soil labile organic carbon fractions and microbial properties.Pedosphere, 2016, 26(6): 829-838.

[42] GUO Z C, ZHANG J B, FAN J, YANG X Y, YI Y L, HAN X R,WANG D Z, ZHU P, PENG X H.Does animal manure application improve soil aggregation? Insights from nine long-term fertilization experiments.Science of the Total Environment, 2019, 660: 1029-1037.

[43] 黄杰.水旱轮作体系下水—旱转换过程中土壤养分变化规律研究[D].成都: 四川农业大学, 2015.HUANG J.Study on the Change of Soil Nutrients in the Process of Water and Dry Rotation System[D].Chengdu: Sichuan Agricultural University, 2015.(in Chinese)

[44] 曾希柏, 孙楠, 高菊生, 王伯仁, 李莲芳.双季稻田改制对作物生长及土壤养分的影响.中国农业科学, 2007, 40(6): 1198-1205.ZENG X B, SUN N, GAO J S, WANG B R, LI L F.Effects of cropping system change for paddy field with double harvest rice on the crops growth and soil nutrient.Scientia Agricultura Sinica, 2007,40(6): 1198-1205.(in Chinese)

[45] 李小涵, 郝明德, 王朝辉, 李利利.农田土壤有机碳的影响因素及其研究.干旱地区农业研究, 2008(3): 176-181.LI X H, HAO M D, WANG Z H, LI L L.Factors affecting soil organic carbon in cropland and their regulation.Agricultural Research in the Arid Areas, 2008(3): 176-181.(in Chinese)

[46] 周萍.南方典型稻田土壤有机碳固定机制研究[D].南京: 南京农业大学, 2009.ZHOU P.A study on soil carbon sequestration fate in typical paddy soils from south China[D].Nanjing: Nanjing Agricultural University,2009.(in Chinese)

[47] XU S X, SHI X Z, ZHAO Y C, YU D S, LI, C S, WANG S H, TAN M Z, SUN W X.Carbon sequestration potential of recommended management practices for paddy soils of China, 1980-2050.Geoderma, 2011, 166: 206-213.

[48] 周晓宇, 张称意, 郭广芬.气候变化对森林土壤有机碳贮藏影响的研究进展.应用生态学报, 2010(7): 1867-1874.ZHOU X Y, ZHANG C Y, GUO G F.Effects of climate change on forest soil organic carbon storage: A review.Chinese Journal of Applied Ecology, 2010(7): 1867-1874.(in Chinese)

[49] 郭广芬, 张称意, 徐影.气候变化对陆地生态系统土壤有机碳储量变化的影响.生态学杂志, 2006(4): 435-442.GUO G F, ZHANG C Y, XU Y.Effects of climate change on soil organic carbon storage in terrestrial ecosystem.Chinese Journal of Ecology, 2006(4): 435-442.(in Chinese)

[50] LIN Y X, YE G P, KUZYAKOV Y, LIU D Y, FAN J B, DING W X,Long-term manure application increases soil organic matter and aggregation, and alters microbial community structure and keystone taxa.Soil Biology and Biochemistry, 2019, doi: https://doi.org/10.1016/j.soilbio.2019.03.030.

[51] YANG F, TIAN J, MEERSMANS J, FANG H J, YANG H, LOU Y L,LI Z F, LIU K L, ZHOU Y, BLAGODATSKAYA E, KUZYAKOV Y.Functional soil organic matter fractions in response to long-term fertilization in upland and paddy systems in South China. Catena,2018, 162: 270-277.

[52] 逯非, 王效科, 韩冰, 欧阳志云, 段晓男, 郑华.中国农田施用化学氮肥的固碳潜力及其有效性评价.应用生态学报, 2008, 19(10):2239-2250.LU F, WANG X K, HAN B, OUYANG Z Y, DUAN X N, ZHENG H.Assessment on the availability of nitrogen fertilization in improving carbon sequestration potential of China's cropland soil. Chinese Journal of Applied Ecology, 2008, 19(10): 2239-2250.(in Chinese)

[53] 马麟英, 梁月兰, 韦国钧, 梁运.东兰县林地土壤有机质含量与土壤容重的相关性分析.湖北农业科学, 2014, 53(1): 59-63.MA L Y, LIANG Y L, WEI G J, LIANG Y.Studies of relations between soil organic matter content and soil bulk density in different soil level in Donglan County.Hubei Agricultural Sciences, 2014,53(1): 59-63.(in Chinese)

[54] 田耀武, 黄志霖, 肖文发, 王宁, 刘晶.三峡库区兰陵溪流域森林土壤有机碳、有机质与容重间的回归模型.华南农业大学学报,2006, 37(1): 89-95.TIAN Y W, HUANG Z L, XIAO W F, WANG N, LIU J.Organic carbon, organic matter and bulk density regression models for forest soils in Lanlingxi watershed, Three Gorges Reservoir area.Journal of South China Agricultural University, 2006, 37(1): 89-95.(in Chinese)

[55] CALLESEN I, LISKI J, RAULUND-RASMUSSEN K, OLSSON M T, TAU ‐STRAND L, VESTERDAL L, WESTMAN C J.Soil carbon stores in Nordic well-drained forest soils – relationships with climate and texture class.Global Change Biology, 2002, 9: 358-370.

[56] WU H B, GUO Z T, PENG C H.Land use induced changes of organic carbon storage in soils of China.Global Change Biology, 2003, 9:305-315.

Spatio-Temporal Variations of Soil Organic Matter in Paddy Soil and Its Driving Factors in China

LI DongChu1, 4, HUANG Jing2, MA ChangBao3, XUE YanDong3, GAO JuSheng2, WANG BoRen2,ZHANG YangZhu1, LIU KaiLou2, 4, HAN TianFu2, ZHANG HuiMin2
(1 College of Resources and Environment, Hunan Agricultural University, Changsha 410128; 2Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences/National Engineering Laboratory for Improving Quality of Arable Land, Beijing 100081; 3Center of Arable Land Quality Monitoring and Protection, Ministry of Agriculture and Rural Affairs,Beijing 100125; 4Jiangxi Institute of Red Soil/National Engineering and Technology Research Center for Red Soil Improvement,Jinxian 331717, Jiangxi)

Abstract:【Objective】Spatio-temporal variation characteristics of soil organic matter (SOM) in paddy fields were evaluated in China, which would provide support for the improving the quality of cultivated land, adjustment of planting structure and climate change.【Method】The characteristics of SOM change and its driving factors were analyzed based on 338 national long-term monitoring sites, which were carried out in different rice planting regions from 1988 to 2017.The effects of SOM on soil bulk density were also studied.【Result】The SOM content of the paddy field in China increased by 3.49 g·kg-1 in averaged in the past 30 years, which was increased at the rate of 0.09-0.12 g·kg-1 annually.The average increased rate of SOM content ranked from high to low were Middle of Yangtze River region (MYR), South of China (SC), Northeast of China (NE), Southwest of China (SW) and Lower of Yangtze River (LYR).The average SOM content of the paddy field was 32.4 g·kg-1 in China at present, which was highest in MYR, followed by SC, NE, SW and LYR.The SOM content were responded by climate, soil types, nitrogen (N) fertilizer input,and cropping system.There was a significantly negative correlation relationship between SOM content and annual mean temperature in SW region and high latitude regions (P＜0.05), and it was a positive correlation relationship in East region and low latitudes region (P＜0.05).The average SOM content in gleyic paddy soil was significantly higher than that in other types of paddy soil.N fertilizer input at 200-300 kg N·hm-2·a-1 was conducive to SOM accumulation.There was strong response between soil bulk density and SOM, plough layer depth and SOM significantly (P＜0.01).【Conclusion】The SOM content of the paddy field in China showed an upward trend in the past 30 years.The average increased rate of SOM content annually showed an increasing trend from South to North of China.The average annual temperature, annual precipitation, paddy soil types and nitrogen fertilizer application were the main driving factors to the SOM content of paddy field in China.

Key words: paddy soil; SOM; national cultivated land monitoring sites; driving factor; nitrogen fertilizer applied rate; China