应用CNMM-DNDC模拟小流域土壤水力侵蚀和颗粒态碳氮磷迁移

水力侵蚀是我国分布最广、危害最严重的水土流失类型,由此导致的土壤和生源要素流失会引起土壤质量下降、河流湖泊富营养化等生态环境问题。水文生物地球化学模型水力侵蚀过程的优化,对定量评估水力侵蚀强度并识别其空间分布特征具有重要意义。因此,本研究以土壤水力侵蚀物理模型ROSE为基础,在分布式流域水碳氮耦合模型CNMM-DNDC中拓展了土壤水力侵蚀流失和颗粒态碳氮磷元素富集模块。以川中丘陵区盐亭截流紫色土农林复合小流域(简称“截流小流域”)为例,开展了地块尺度验证和小流域尺度的土壤水力侵蚀和颗粒态碳氮磷流失的高时空分辨率模拟应用。结果表明:增加了土壤水力侵蚀和颗粒态碳氮磷富集模块的CNMM-DNDC模型,能够较好地模拟截流小流域玉米-小麦轮作径流小区产沙量和颗粒态氮流失量的季节动态特征,模型验证阶段模拟值和观测值均表现出显著线性拟合关系(R^(2)分别为0.83和0.85),改进的CNMM-DNDC模型能够用于模拟评估流域水土流失和颗粒态碳氮磷流失强度空间分布。截流小流域由地表径流引起的土壤水力侵蚀、颗粒态碳氮磷流失强度与土地利用和地形特征密切相关,其中坡耕地和居民区土壤水力侵蚀最严重,同时也是截流小流域颗粒态碳氮磷流失的主要来源。模型的模拟结果显示,2004—2006年截流小流域栅格尺度土壤侵蚀模数均值分别为400 t∙km^(−2)∙a^(−1)、701 t∙km^(−2)∙a^(−1)和1550 t∙km^(−2)∙a^(−1),属于轻度水力侵蚀区域。2004—2006年截流小流域由水力侵蚀引起的颗粒态碳年总流失量的平均值分别为63.9 kg(C)∙hm^(−2)∙a^(−1)、107.2 kg(C)∙hm^(−2)∙a^(−1)和200.2 kg(C)∙hm^(−2)∙a^(−1),同期颗粒态氮年总流失量的平均值分别为6.3 kg(N)∙hm^(−2)∙a^(−1)、10.5 kg(N)∙hm^(−2)∙a^(−1)和19.5 kg(N)∙hm^(−2)∙a^(−1),颗粒态磷年总流失量平均值分别为0.9 kg(P)∙hm^(−2)∙a^(−1)、1.5 kg(P)∙hm^(−2)∙a^(−1)和2.8 kg(P)∙hm^(−2)∙a^(−1)。本研究拓展了CNMM-DNDC模型对土壤水力侵蚀产沙和随之的颗粒态碳氮磷迁移的模拟功能,为水土流失和面源污染的研究和防控提供先进的模型工具。