【摘要】【目的】通过对1987—2016年海南岛农牧生产体系磷元素流动时空特征及环境效应进行定量分析,研究其流动过程和规律,探讨农牧生产体系磷素的优化管理途径,为海南岛农牧业可持续发展提供科学依据。【方法】研究基于食物链养分流动模型(NUtrient flows in Food chains, Environment and Resources use,NUFER),通过统计数据、文献检索、实地调研,并结合Origin等软件,定量计算海南岛农牧生产体系的磷素输入、输出、利用率及其环境效应,并通过情景分析探索海南岛农牧生产体系磷素的可持续利用途径。【结果】30年间海南岛农田生产子系统磷素总投入量从21.34 Gg增至81.19 Gg,总输出量由6.20 Gg增至18.20 Gg,化肥作为该系统磷素主要来源,输入量由19.01 Gg增至79.23 Gg,作物产品作为农田磷素主要输出项,30年间由5.25 Gg增至15.48 Gg。动物生产子系统磷素总输入量由11.40 Gg增至15.31 Gg,总输出量由9.63 Gg增至11.90 Gg,其中外源饲料磷素输入量由1987年的10.97 Gg增至2016年的14.77 Gg,动物产品输出量30年间增长了4.95 Gg。秸秆还田量和作物饲用量分别增加了0.37和0.26 kg·hm~(-2),粪尿还田量则减少了0.80 kg·hm~(-2)。空间分布上,澄迈、定安等地30年来磷素输入和输出量较高,五指山、琼中等地较低。就磷素损失情况来看,1987—2016年,海南岛单位耕地面积土壤磷盈余量由35.00 kg·hm~(-2)增至147.40 kg·hm~(-2),2016年土壤磷素盈余量较大的是琼海、澄迈、保亭和临高,分别为372.79、279.82、194.14和181.09 kg·hm~(-2)。磷的其他损失途径为土壤侵蚀、径流和淋洗,损失量在1.21—5.85 kg·hm~(-2)。畜禽粪便单位耕地面积承载量维持在3.83—5.77 kg·hm~(-2)。30年来,磷素利用率增长缓慢,其中农田生产子系统磷素利用率由13.01%增至13.86%,动物生产子系统磷素利用率由4.78%增至7.62%,农牧结合体系磷素利用率由10.78%增至13.09%。情景分析结果显示,保证农牧生产体系各子系统间的协调稳定发展以及通过科学的养分管理方式提高资源的循环利用率对促进海南岛农牧业发展意义重大。【结论】受农牧生产体系规模、区域发展以及管理方式等因素影响,海南岛农牧生产体系环境损失情况严重,磷素利用率较低,体系出现了较严重的分离。因此,在海南岛未来的农牧生产中,应优化技术手段和管理措施,如控制磷素的过量输入,减少粪尿的直接排放,提高秸秆和粪尿循环利用率。同时也应促进农田生产子系统与动物生产子系统间的协调关系,走农牧结合的可持续发展道路。