中国稻田生态系统固碳效应模拟研究.pdf

引文格式：卢欣晴,张秀英,汪振,等.中国稻田生态系统固碳效应模拟研究J.云南农业大学学报(自然科学),2024,39(1):141152.DOI:10.12101/j.issn.1004-390X(n).202210029中国稻田生态系统固碳效应模拟研究*卢欣晴1,2，张秀英1,3*，汪振1,3，李升峰4，郭文勇5(1.南京大学国际地球系统科学研究所，江苏南京210023；2.中山大学大气科学学院，广东珠海519000；3.江苏地理信息资源开发与应用协同创新中心，江苏南京210046；4.南京大学地理与海洋科学学院，江苏南京210023；5.鲁西化工集团有限公司，山东聊城252000)摘要:【目的】研究稻田生态系统固碳量的动态变化，为制定与推广耕地固碳减排措施提供依据。【方法】基于 3 年稻田受控试验数据，校验生物地球化学脱氮分解作用(denitrification-decomposition，DNDC)模型，并估算 20182020 年间中国稻田生态系统的固碳量；分析稻田生态系统土壤及作物固碳量的时空分布，评估农田生态系统在实现“碳中和”中的贡献。【结果】DNDC 模型模拟的水稻产量及土壤有机碳含量与田间观测结果一致性较高(相关系数均大于等于 0.77，平均绝对误差和归一化均方根误差均小于 13%)；20182020 年中国水稻固碳总量为 523.29Tg，作物对总固碳量的贡献率达 97%以上。江西、湖南、广东和广西的双季稻固碳能力较强，而云南、浙江和海南的相对较弱。单季稻固碳量分布呈现北高南低的空间分布特征，其中黑龙江、江苏和湖北的固碳量较高。湖南、福建、江西和广西的双季稻固碳量大于单季稻，双季稻固碳量占总固碳量的 60%以上。【结论】DNDC 模型可用于模拟稻田生态系统的碳收支，中国稻田生态系统是重要的碳汇，且呈现出较强的空间异质性。关键词:固碳量；脱氮分解(DNDC)模型；稻田生态系统；土壤有机碳中图分类号:Q178.516文献标志码:A文章编号:1004390X(2024)01014112Study of Carbon Sequestration Effect of PaddyEcosystem in ChinaLUXinqing1,2，ZHANGXiuying1,3，WANGZhen1,3，LIShengfeng4，GUOWenyong5(1.InternationalInstituteforEarthSystemScience,NanjingUniversity,Nanjing210023,China;2.SchoolofAtmosphericSciences,SunYat-SenUniversity,Zhuhai519000,China;3.JiangsuCenterforCollaborativeInnovationinGeographicalInformationResourceDevelopmentandApplication,Nanjing210046,China;4.SchoolofGeographyandOceanScience,NanjingUniversity,Nanjing210023,China;5.LuxiChemicalGroupCo.,Ltd.,Liaocheng252000,China)Abstract:Purpose To study the dynamics of carbon sequestration in the paddy ecosystem,providingabasisforformulatingandpromotingmeasuresofcarbonsequestrationandemissionre-ductionincultivatedland.MethodsBasedonthedataofthree-yearcontrolledexperimentsinpaddyfields,thedenitrification-decomposition(DNDC)modelwereverified.Thecarbonsequestra-收稿日期：2022-10-21修回日期：2024-01-29网络首发日期：2024-03-27*基金项目：国家自然科学基金项目(42371378，41471343)；江苏省自然科学基金项目(BK20211156)；自然资源部城市国土资源监测与仿真重点实验室开放课题(KF-2021-06-028)。作者简介：卢欣晴(1997)，女，广东东莞人，在读博士研究生，主要从事生态过程模型与稻田生态系统碳循环研究。E-mail：Xinqing_L*通信作者 Correspondingauthor：张秀英(1977)，女，河北唐山人，博士，教授，主要从事大气环境遥感与农业生态系统碳氮循环过程模拟研究。E-mail：网络首发地址：https:/ and Guangxi,while the regions with lower carbon sequestration fordoublecroppingricewereYunnan,ZhejiangandHainan.Thecarbonsequestrationofsinglecrop-pingricewaslowerinthesouthandhigherinthenorth,withhighercarbonsequestrationinHei-longjiang,JiangsuandHubei.ThecarbonsequestrationofdoublecroppingriceinHunan,Fujian,Ji-angxiandGuangxiwashigherthanthatofsinglecroppingrice,andthecarbonsequestrationofdoublecroppingriceaccountedformorethan60%ofthetotalcarbonsequestration.ConclusionTheDNDCmodelcanbeusedtosimulatecarbonbudgetsinthericepaddyecosystem.TheChinesepaddyfieldsplayanimportantroleincarbonsequestration,withastrongspatialvariation.Keywords:carbonsequestration;denitrification-decomposition(DNDC)model;paddyecosystem;soilorganiccarbon农田是陆地生态系统的重要组成部分1，水稻田在农田生态系统中占据十分重要的地位。中国水稻种植面积超过 2.7107hm2，其碳汇主要包括土壤有机碳(soilorganiccarbon，SOC)及植株生长的碳吸收量。SOC 作为全球碳库中最活跃的部分，是土壤肥力的重要表征，其微小的变化将导致大气 CO2浓度剧烈变化，甚至严重影响全球碳平衡2。水稻 SOC 容易受到人为活动的强烈干扰，同时在较短时间内可对其进行人为调节3。稻田生态系统是重要的碳汇资源，合理利用秸秆资源，增加外源有机肥投入，可提高农田有机质含量，增加农田生态系统碳汇总量3-4。在当前CO2排放增长的环境背景下，利用农田生态系统固定大气 CO2是应对全球气候变化的重要途径，也是中国实现双碳目标的重要手段之一4-6。目前，农田生态系统碳收支核算方法主要包括生命周期7-8、遥感反演9-10、涡度相关通量观测11-13和生态系统过程模型14-15等。作为获取地表信息的有力手段，卫星遥感能准确、及时地获取多时空尺度的净初级生产力时空分布特征，弥补传统地面观测的不足16；但受卫星观测成像频率的影响，数据时空分辨率不高，同时大多数生态参数反演模型缺乏明确的生理生态机制，难以清晰解释区域碳平衡的变化机制17。涡度相关观测虽然具有较高准确性，但受制于有限的观测点，往往局限于特定的生态系统或局部区域18。生命周期法的计算过程详细而准确，包含系统生命周期的各个环节，适用于微观层面的碳足迹核算19-20，但在划定系统边界时可能具有主观性，且需收集大量数据，易产生数据遗漏现象21，在核算区域尺度农田碳源汇过程中，估算作物固碳量主要根据经济系数、碳吸收率和产量，衡量固碳措施通常只采用同一系数来估算碳汇提升的能力。稻田生态系统的碳循环过程比较复杂，除受作物类型影响外，还与作物生长固定的 SOC、作物生长吸收和排放的碳、气温和降水等气象条件以及当地的土壤理化性质相关22-23。因此，利用单一的碳吸收率和产量数据核算农田系统碳吸收量，存在较高的不确定性，核算过程中还需要考虑气象条件和土壤理化性质的时空异质性对农田生态系统碳收支时空分布格局的影响。基于过程的模型已成为国内外学者研究农田生态系统碳循环的重要工具，过程模型可以定量模拟植物在一定气象条件和土壤环境下的生长过程24，并在区域乃至全球尺度上准确模拟农田生态系统碳迁移过程及碳收支情况25-27。生物地球化学脱氮分解作用(denitrification-decompositi-on，DNDC)模型能够准确模拟作物生长情况及产142云南农业大学学报第39卷量、SOC 及其动态变化以及温室气体排放等28-31，已在世界范围内得到验证并开展了广泛的应用32-39。DNDC 模型由利用生态驱动因子模拟土壤环境动态变化的土壤气候、农作物生长和土壤有机质分解模型以及用于计算温室气体和痕量气体排放情况的硝化作用、反硝化作用和发酵作用模型共同构建35,37,40-42。随着模型的发展及参数优化，DN-DC 模型已被广泛应用于湿地、草原和森林等生态系统26,43-44，其碳循环过程考虑了来自大气、土壤母质、施肥以及有机肥中的碳和氮，根据气温计算植物的水分胁迫情况以及受温度、水分和土壤碳氮供给量影响的植物实际生长量，并将植物实际生长量分配至根、茎、叶和籽粒中，最终输出作物生长吸收的碳、土壤中淋失和残留的碳以及排放到大气中的碳14。稻田生态系统的碳汇主要来自水稻生长期间的植株固碳和土壤固碳45。水稻的植株固碳是指水稻通过光合作用吸收 CO2，合成有机物并将碳固定在作物内的过程。该过程中水稻的自身固碳仅限于水稻从播种到收获的生长期，局限在稻田生态系统内，而不包括水稻收获后作为副产品处理或快速分解和燃烧的问题。土壤碳库作为陆地生态系统中最大、最活跃的碳库之一，在全球碳循环过程中扮演着重要的角色46。土壤碳主要包括两类，一类是不断与地球大气进行碳交换且相对不稳定的有机碳，另一类是相对稳定的无机碳47，土壤固碳以及土壤碳库的研究主要集中在前者。土壤碳固存是指土壤从大气中捕获碳并稳定贮存于碳库中45,48。目前稻田生态系统碳汇的研究尺度主要以区域尺度为主49-51，而全国尺度稻田生态系统碳汇估测及固碳量计算的相关报道仍较少52。本研究基于 3 年的水稻田受控试验数据，对 DNDC 模型参数进行校验，利用 DNDC 模型估算 20182020年中国稻田生态系统的固碳量，并分析其时空分布格局，以期为区域稻田生态系统固碳减排提供重要决策依据。1 材料与方法1.1稻田受控试验受控试验于 2015 年 7 月2018 年 7 月在江苏省句容市淮源农场(3148N，11913E)进行。研究区属于典型的亚热带气候，试验期间平均年降雨量 1609.5mm，日平均最低气温 20.7，最高气温 27.9。根据联合国粮农组织土壤分类53，该水稻田土壤为水耕人为土。试验田土壤(020cm)的理化性质为：pH 值 5.8，土壤容重 1.25g/cm3，总碳含量 9.12g/kg，总氮含量 0.96g/kg，黏粒占14.63%，粉粒占 78.62%，砂粒占 6.74%。稻田受控试验期间，水稻秧龄达 30d 后移栽至试验田。水分管理采用中期排水晒田措施54：淹灌晒田复淹干湿交替。在水稻移栽前7d 进行翻地和淹水，并维持淹水状态 1 个月；季中排水的时间通常为 714d；复淹后，采用干湿交替灌溉直至水稻收获前7d。作物管理按照当地水稻冬小麦轮作制度种植。水稻移栽和收获时间见表 1。每年水稻成熟时收割，收割时尽可能保证水稻根系的完整性，并将根系清洗干净。水稻晾干水分后将其分为根、茎、叶和穗4 个部分，测定每个部分的干质量。表 1 水稻田受控试验设计Tab.1Designofcontrolledexperimentsinpaddyfields移栽日期(yyyy-mm-dd)transplantingdate收获日期(yyyy-mm-dd)harvestdate生长时间/dgrowthtime2015-07-012015-11-121352016-06-302016-11-071312017-06-202017-11-111501.2DNDC 模型参数输入及验证1.2.1DNDC 模型模拟的参数输入DNDC 模型模拟包括点位模拟和区域模拟，模拟的前提是输入对应的气象数据、土壤属性、耕作参数等54。区域模拟前，首先对 DNDC 模型进行参数优化和验证，从而确保该模型适用于模拟研究区域的生物地球化学过程，详细的参数优化方案及点位验证结果见文献 55-56。本研究运用 DNDC 模型(9.5 版本)对中国稻田固碳量进行区域模拟，研究设定的模拟单元为 1km1km的规则网格。气象数据来源于中国气象局 20182020 年699 个气象站点的地面气候资料日值数据集(ht-tps:/ 20182020 年间，全第1期卢欣晴，等：中国稻田生态系统固碳效应模拟研究143国气温为 10.110.3，较常年偏高 0.50.8；全国降水量为 645.5673.8mm，比常年偏多 3%7%，其中，东北、西北、华南年降水量偏多，而华北和长江中下游相对偏少，西南略偏少。总体而言，主要粮食产区光、温、水匹配较好，主要粮食作物生产期间的气候条件有利于农业生产57。土壤属性数据来源于协调世界土壤数据库(1.2 版本)(http:/webarchive.iiasa.ac.at/Research/L-UC/External-World-soil-database/)，空间分辨率为1km1km。土地利用类型数据基于 Landsat 卫星获得的 LandsatThematicMapper/EnhancedThem-aticMapperPlus(TM/ETM+)遥感影像(http:/)产生。水稻种植区内各模拟单元的土壤属性由土地利用类型数据和土壤属性数据叠加分析生成。总体而言，中国水稻田 SOC含量相对较低，大多数 SOC 含量低于 15g/kg，变化范围介于 351g/kg，其中黑龙江省和吉林省SOC 含量最高；土壤黏粒含量变异较大，其范围在 4%50%间；土壤容重的变化范围介于 1.21.4g/cm3，广西壮族自治区和贵州省的土壤容重相对较低。水稻种植面积及产量数据来自各省统计年鉴。2018 年中国水稻种植总面积为 3.019107hm2，早稻、晚稻和单季稻分别为 4.79106、5.27106和2.013107hm2，占比分别为16%、17%和67%，平均产量分别为 5.967、5.958 和 7.559t/hm2；2019年水稻种植总面积减少至 2.969107hm2，早稻、晚稻和单季稻的种植面积分别为 4.45106、4.97106和 2.027107hm2，占比分别为 15%、17%和68%，平均产量分别为 5.902、6.049 和 7.561t/hm2。单季稻、早稻和晚稻的平均施氮量分别为150.53、139.91 和 138.98kg/hm258。尿素为中国水稻生产中最常用的氮肥类型，故在 DNDC 模型中以尿素作为输入的氮肥。由于目前缺乏完整的关于中国不同地区水稻田粪肥施用情况的统计资料，故本研究没有考虑粪肥的施用。在水稻收获后，15%的秸秆施入稻田以实现土壤改良。本研究采用单季稻和双季稻 2 种常规的种植制度进行模型模拟。早稻、晚稻和单季稻的平均生育期分别为 87、94 和 109d59-61。1.2.2模型验证本研究依据相关系数(r)、归一化均方根误差(normalizedrootmeansquareerror，nRMSE)和平均绝对误差(meanabsoluteerror，MAE)评价DNDC 模型模拟结果的有效性62-63。为在全国范围内评估校准后 DNDC 模型的性能，本研究收集了各省统计年鉴中的水稻产量数据，并以“有机碳/有机质”和“20102020”为关键词在中国知网上收集了 40 组土壤有机碳观测数据(图 1)，用于验证区域尺度上水稻的产量及 SOC。1.3稻田生态系统固碳量计算1.3.1水稻固碳量水稻固碳量是指水稻光合作用形成的净初级生产量，即收获的籽粒、茎叶和根系生物量中固碳量的总和。水稻单位面积固碳量计算公式为：C=Cg+Cs+Cr。式中：C 为区域水稻单位面积固碳量，kg/hm2；Cg为水稻籽粒固碳量，kg/hm2；Cs为水稻茎叶固碳量，kg/hm2；Cr为水稻根系固碳量，kg/hm2。水稻固碳量的计算公式为：TC=ni=1(SiCi)。SiCi式中：TC 为水稻固碳量，kg；为每个模拟单元的水稻田面积，hm2；为每个模拟单元的水稻单位面积固碳量，kg/hm2；n 为模拟单元数。1.3.2稻田土壤固碳总量稻田土壤固碳总量是指一段时间内稻田生态系统土壤碳储量的变化量，反映了稻田生态系统土壤碳库的变化，其计算公式为：TSOC=ni=1(SiSOCi)。SOCi式中：TSOC 为稻田土壤固碳总量，kg；为每个模拟单元的 SOC 年变化量，kg/hm2。2 结果与分析2.1DNDC 模型模拟结果验证2.1.1点位尺度上的验证DNDC 模型模拟和基于受控试验实测的水稻根(n=15)、茎(n=17)、叶(n=17)、穗(n=25)和总生物量(n=17)的 MAE 均低于 8%，分别为 7.3%、7.7%、7.3%、4.5%和 4.7%。水稻根、茎、叶、穗和总生物量的模拟值与实测值的 r 值均大于0.8(P0.01)，说明 DNDC 模型能够很好地反映水稻各组分生物量的变化，并对其进行模拟。144云南农业大学学报第39卷2.1.2区域尺度上的验证双季稻(图 2a)和单季稻(图 2b)观测和模拟的产量nRMSE 和MAE 较小，且r0.77(P0.01)，表明DNDC 模型在研究区具有较高的模拟精度，可用于模拟水稻产量的动态变化。此外，SOC 的模拟结果与实测数据间的 MAE和 nRMSE 均小于 10%，说明其在水稻 SOC 动态变化的模拟中具有较高的精度和较好的应用前景(图 2c)。50N40N30N20N90E100E110E120E130E0500 km南海诸岛South China Sea islands钓鱼岛Diaoyu Dao赤尾屿Chiwei Yucontrolled experimental site有机碳验证点SOC verification site气象站点meteorological stations水稻种植区rice area受控试验点注：基于自然资源部标准地图服务系统网站下载的 GS(2020)4619 号标准地图制作，底图边界无修改；下同。Note:BasedonthestandardmapofGS(2020)4619downloadedfromthewebsiteofthestandardmapservicesystemoftheMinistryofNaturalResources,theboundaryofthebasemapisnotmodified;thesameasbelow.图 1 气象站点、试验点、水稻种植区和验证点位置图Fig.1Locationmapofmeteorologicalstations,experimentalsites,riceareaandverificationsite7 500a)双季稻 double cropping ricer=0.77 P0.01MAE=5.17%nRMSE=6.77%r=0.84 P0.01MAE=10.78%nRMSE=12.09%r=0.97 P0.01MAE=8.62%nRMSE=9.86%b)单季稻 single cropping ricec)土壤有机碳 soil organic carbon模拟产量/(kghm2)simulated yield模拟产量/(kghm2)simulated yield模拟含量/(kghm2)simulated content实测产量/(kghm2)observed yield实测产量/(kghm2)observed yield实测含量/(kghm2)observed content7 0006 5006 0005 5005 0004 5004 00050 00040 00030 00020 00010 00008 0007 0006 0005 0004 0005 0006 0007 0005 0007 0009 000020 00040 00060 000图 2 区域尺度 DNDC 模型模拟结果验证Fig.2AccuracyassessmentofthesimulationresultsbyDNDCmodelonregionalscale第1期卢欣晴，等：中国稻田生态系统固碳效应模拟研究1452.2中国双季稻固碳量的空间分布作物固碳量和土壤固碳量存在显著的区域差异(图 3a 和 3b)。江西、湖南和广东 3 省的总固碳量排在前列，该 3 省的作物固碳量占双季稻作物总固碳量的 67.65%(图 3c)、土壤固碳量占双季稻土壤总固碳量的 69.27%(图 3d)，主要原因是这 3 省的播种面积占双季稻播种面积的 67.96%。云南、浙江和海南 3 省的作物和土壤固碳能力均较弱，作物平均固碳量仅为 2.49Tg，远低于全国平均值（15.29Tg）；土壤固碳量均低于 80Gg，主要原因是这 3 省的平均播种面积仅占全国双季稻的 5.49%。此外，作物固碳总量还受水稻单产的影响。海南省受地形影响，适宜播种水稻的面积较小，且同时存在产量较低的旱田或涝田，加之水稻管理措施粗放以及病虫害的影响，其水稻单产远低于全国平均水平。50N40N30N20N50N40N30N20N500 km0南海诸岛South ChinaSea islands作物固碳量/Ggcarbon sequestration of crop作物固碳量carbon sequestration of crop土壤固碳量carbon sequestration of soil作物固碳量/Tgcarbon sequestration of crop土壤固碳量/Ggcarbon sequestration of soil云南 Yunnan浙江 Zhejiang海南 Hainan安徽 Anhui福建 Fujian湖北 Hubei广西 Guangxi广东 Guangdong湖南 Hunan江西 Jiangxi云南 Yunnan浙江 Zhejiang海南 Hainan安徽 Anhui福建 Fujian湖北 Hubei广西 Guangxi广东 Guangdong湖南 Hunan江西 Jiangxi云南 Yunnan浙江 Zhejiang海南 Hainan安徽 Anhui福建 Fujian湖北 Hubei广西 Guangxi广东 Guangdong湖南 Hunan江西 Jiangxia)c)e)d)70E80E90E100E110E120E130E140E5.030.1424.2699.68010203040总固碳量/Tgtotal carbon sequestration01020304002004006008001 0001 200500 km0南海诸岛South ChinaSea islandsb)70E80E90E100E110E120E130E140E0.401.851.092.78土壤固碳量/Ggcarbon sequestration of soil钓鱼岛Diaoyu Dao赤尾屿Chiwei Yu钓鱼岛Diaoyu Dao赤尾屿Chiwei Yu图 3 中国双季稻生态系统固碳量空间分布Fig.3SpatialdistributionofcarbonsequestrationindoublecroppingriceinChina146云南农业大学学报第39卷双季稻区的北部省份(湖北和安徽)及东部省份(浙江和福建)的土壤固碳量较低(图 3b)；湖南、江西和广东的土壤固碳量较高；华南地区由于水稻土分布面积广泛，在有机碳固定中的贡献较大，土壤固碳量相对较高。江西、湖南、广东和广西的双季稻固碳能力较强，江西和湖南的固碳量接近 40Tg，广东和广西的固碳量均大于 25Tg(图 3e)。另外，作物固碳量大于土壤固碳量，作物固碳量占双季稻总固碳量的 97.21%。2.3中国单季稻固碳量的空间分布由图 4 可知：单季稻作物固碳量和土壤固碳量总体上呈现北高南低的特征。东北地区(黑龙50N40N30N20N500 km0南海诸岛South ChinaSea islands作物固碳量/Ggcarbon sequestration of cropa)c)e)d)70E80E90E100E110E120E130E140E25.200.0280.60 214.0550N40N30N20N500 km0南海诸岛South ChinaSea islands土壤固碳量/Ggcarbon sequestration of soilb)70E80E90E100E110E120E130E140E0.475.703.7615.85北京 Beijing西藏 Xizang山西 Shanxi甘肃 Gansu天津 Tianjin河北 Hebei宁夏 Ningxia新疆 Xinjiang陕西 Shaanxi上海 Shanghai广西 Guangxi山东 Shandong内蒙古 Inner Mongolia福建 Fujian浙江 Zhejiang辽宁 Liaoning贵州 Guizhou重庆 Chongqing云南 Yunnan河南 Henan吉林 Jilin江西 Jiangxi湖南 Hunan安徽 Anhui四川 Sichuan湖北 Hubei江苏 Jiangsu黑龙江 Heilongjiang北京 Beijing西藏 Xizang山西 Shanxi甘肃 Gansu天津 Tianjin河北 Hebei宁夏 Ningxia新疆 Xinjiang陕西 Shaanxi上海 Shanghai广西 Guangxi山东 Shandong内蒙古 Inner Mongolia福建 Fujian浙江 Zhejiang辽宁 Liaoning贵州 Guizhou重庆 Chongqing云南 Yunnan河南 Henan江西 Jiangxi吉林 Jilin湖南 Hunan安徽 Anhui四川 Sichuan湖北 Hubei江苏 Jiangsu黑龙江 Heilongjiang北京 Beijing西藏 Xizang山西 Shanxi甘肃 Gansu河北 Hebei宁夏 Ningxia天津 Tianjin广西 Guangxi陕西 Shaanxi新疆 Xinjiang上海 Shanghai福建 Fujian山东 Shandong内蒙古 Inner Mongolia浙江 Zhejiang云南 Yunnan重庆 Chongqing贵州 Guizhou江西 Jiangxi辽宁 Liaoning河南 Henan安徽 Anhui湖南 Hunan吉林 Jilin四川 Sichuan湖北 Hubei江苏 Jiangsu黑龙江 Heilongjiang作物固碳量carbon sequestration of crop土壤固碳量carbon sequestration of soil作物固碳量/Ggcarbon sequestration of crop020 00040 00060 000总固碳量/Ggcarbon sequestration020 00040 00060 000土壤固碳量/Ggcarbon sequestration of soil01 0002 0003 0004 000钓鱼岛Diaoyu Dao赤尾屿Chiwei Yu钓鱼岛Diaoyu Dao赤尾屿Chiwei Yu图 4 中国单季稻生态系统固碳量空间分布Fig.4Spatialdistributionofcarbonsequestrationinsinglecroppingrice第1期卢欣晴，等：中国稻田生态系统固碳效应模拟研究147江、吉林)作物和土壤固碳能力较强，主要原因为东北地区气候湿润，降水充沛，适合作物生长，净初级生产力较高，作物固碳量较大，凋落物输入量大；同时温度较低导致凋落物分解速度较慢，利于积累有机物。北京、西藏和山西播种面积仅占单季稻播种面积的 0.01%，故作物和土壤固碳量较低；此外，地处高原的西藏由于海拔高、气温低，气候干燥加之自然灾害较多导致水稻单产较低，水稻作物固碳能力较弱。单季稻总固碳量为 366.04Tg，省份平均固碳量为 13.07Tg。黑龙江、江苏、湖北、四川和安徽的固碳能力较强，平均固碳量为 45.62Tg，远高于全国平均水平。另外，作物固碳量大于土壤固碳量，作物固碳量占单季稻总固碳量的 96.82%。2.4中国稻田生态系统总固碳量的空间分布水稻分布地域广阔，自然条件复杂多样，20182020 年中国水稻作物固碳量的空间分布高度异质，总体呈现东北、华中地区高，东南、西南地区低的空间分布特点(图 5a)。土壤固碳量空间分布大体呈现北高南低的地域分异规律(图 5b)：北方地区(黑龙江)土壤固碳能力强，南方地区(江苏、湖北和广东)土壤固碳能力处于中等水平，而西南、华东及华北地区土壤固碳能力处于较低水平。水稻总固碳量为 523.29Tg，省份平均固碳量为 17.44Tg。从空间分布上看(图 5c)，固碳量较大的区域集中在黑龙江、河南南部、湖南、江西、安徽、江苏等；而固碳量较小的区域位于浙江、宁夏北部、陕西南部等。黑龙江、湖南、江西、湖北和江苏的固碳能力较强，其作物固碳量占全国作物固碳量的一半以上(图 5d)，土壤固碳量占全国土壤总固碳量的 62.96%(图 5e)，固碳总量达到 285.02Tg，固碳能力排在全国前列(图 5f)。湖南、福建、江西和广西的双季稻固碳量大于单季稻，双季稻固碳量约占总固碳量的 60.26%92.19%。除广西双季稻播种面积大于单季稻外，湖南、福建和江西的双季稻播种面积均小于单季稻。3 讨论3.1与其他研究的估算结果比较水稻作物固碳量通常在 0.1150.917Pg，占农田作物固碳的 40%46%64-65，本研究估算的水稻作物固碳量处于该估算范围内，与王雅楠等66计算的结果相近。模型的选取以及输入参数的差异会得到不同的计算结果。本研究计算得到的稻田土壤碳储量为 0.78Pg，略低于 PAN 等67计算得到的 0.85Pg，分析其可能原因包括以下 2 个方面：(1)DNDC 模型计算的是土壤深度 020cm 的碳储量，而 PAN 等67估算的是基于实测土壤剖面资料所记录的土壤深度；(2)PAN 等67估算的有机碳含量是将有机质含量乘以 0.58 转化而得，而本研究是直接以碳(C)为计算单位，故两者略有区别，但总体而言，本研究估算的结果与其他相关研究基本一致。本研究对 DNDC 模型进行了参数优化、校准和验证，模拟值与实测值在点位模拟和区域模拟上均具有较好的一致性，表明 DNDC 模型具有模拟稻田生态系统作物和土壤固碳量的能力。然而，利用 DNDC 模型模拟全国水稻固碳量的时空分布格局也存在一定的不确定性68。首先，在区域尺度上，DNDC 模型模拟的不确定性主要来源于土壤属性的空间异质性，模型默认每个模拟单元内的土壤属性是均一的，这与现实情况不符68，因此，在今后的研究中还需要大量田间试验数据来优化模型内部参数，减少系统误差，从而更全面和准确地提高模拟精度69。其次，长期定位试验数据还比较欠缺，全国水稻分为单、双季稻不够精细等因素也会导致模型模拟的不确定性，因此，试验数据的逐步积累以及不断增加对稻田生态系统固碳能力及其微生物土壤作物交互作用机制的认知才能逐步减少估算的不确定性。3.2农田作物碳储量的碳源汇效应与森林和草地生态系统相比，农田生态系统的碳汇较弱。方精云等70、赵宁等71和杨元合等72学者认为：农作物收获期短，作物生物量的固碳作用不显著，因此把农作物生物量的碳汇设为零。然而，近年研究结果表明农田作物生物量的碳汇为 0.62.0Pg64-65,73，且通过估算作物碳储量动态变化可知：随着作物栽培技术的逐步优化和农业投入的快速增加，作物生产和固碳能力不断增强74。农作物自身的固碳能力在产量不断增加的背景下迅速提升，作物生产系统正逐步成为一个庞大的碳库，对于维持农业土壤固碳、增加土壤有机碳储量、确保粮食安全等至关重要4。稳定的土壤有机碳形成时间较长，且农作物大多为一年生75，因此在农田生态系统中土壤碳148云南农业大学学报第39卷储量(密度)普遍大于作物植被碳储量(密度)74，但在区域尺度上，作物具有显著的碳汇功能，中国农作物年均固碳量高达 6108t76。随着农业生产技术的提升，农作物生产能力和农作物碳蓄积量也在不断提高，同时，作物的生产能力以及秸秆还田模式也会影响土壤碳库水平74。因此，在农田生态系统碳储量及碳收支估算中应充分重视作物碳储量，并具体分析不同时空尺度的碳源汇效应。稻田生态系统在中国陆地生态系统碳循环中占据着重要地位，稻田的耕作制度和田间管理模式不断优化，尤其是土地利用变化和水稻新品种的推广和运用将影响稻田生态系统的固碳过程和固碳能力77。在今后的研究中，阐明各要素对稻田生态系统固碳功能的影响机制，提升核算与预测稻田碳中和的能力78，对于固碳减排和稻田可持续发展具有重要意义。50N40N30N20N500 km0南海诸岛South ChinaSea islands作物固碳量/Ggcarbon sequestration of cropa)c)e)f)70E80E90E100E110E120E130E140E27.720.0278.92 214.0550N40N30N20N500 km0南海诸岛South ChinaSea islands总固碳量/Ggtotal carbon sequestration90E100E110E单季稻single cropping rice双季稻double cropping rice单季稻single cropping rice双季稻double cropping rice单季稻single cropping rice双季稻double cropping rice120E130E28.950.0285.06223.5750N40N30N20N500 km0南海诸岛South ChinaSea islands土壤固碳量/Ggcarbon sequestration of soilb)d)70E80E90E100E110E120E130E140E0.555.703.5915.85北京 Beijing西藏 Xizang山西 Shanxi甘肃 Gansu天津 Tianjin河北 Hebei宁夏 Ningxia新疆 Xinjiang陕西 Shaanxi上海 Shanghai山东 Shandong内蒙古 Inner Mongolia海南 Hainan福建 Fujian辽宁 Liaoning贵州 Guizhou浙江 Zhejiang重庆 Chongqing云南 Yunnan河南 Henan吉林 Jilin广西 Guangxi广东 Guangdong四川 Sichuan安徽 Anhui湖北 Hubei江苏 Jiangsu江西 Jiangxi湖南 Hunan黑龙江 Heilongjiang北京 Beijing西藏 Xizang山西 Sh