磷(P)是维持生物体基本机能不可或缺的重要营养元素，其在自然界中主要的存在形式——磷矿石则是重要的、难以再生的非金属矿物资源。人们通过磷矿开采、磷化工生产、农业种植和动物养殖等过程，将矿物磷转化为植物磷和动物磷，以满足人类生存和社会运转的磷需求。除了制造磷肥，磷矿石还可以用来制造黄磷、磷酸、有机磷及磷酸盐等各种产品，广泛地用于医药、食品、火柴、染料、制糖、陶瓷、国防等工业部门。随着人口增长、饮食结构的变化以及经济的高速发展，磷资源的需求量不断增加，磷矿的开采和消费量逐年增长，磷短缺也与气候变化、水短缺、氮管理一起并列成为21世纪重要的全球问题^[1]。根据美国地质勘探局(United States Geological Survey，USGS)的最新数据显示，全球2013 年的磷矿石储量约为670 亿t^[2]，折合80—90 亿t磷，乐观估计可以支持全球200—300a的需求，但是如果考虑开采成本，经济储量则仅能够支持全球约50—100年的需求^{[3, 4, 5, 6]}。中国查明的磷矿资源储量约186.3 亿t(折标矿)，居世界第二位，但是中国的磷矿资源平均品位低(仅为17%，世界平均为30%)，80%以上为中低品位磷矿，经济储量仅能维持国内50—70a的消费需求，形势十分严峻^{[7, 8]}。另一方面，磷矿开采、磷化工、以及各种含磷产品消费过程产生的含磷废水和废物的排放也带来了一系列的环境问题，包括地表水的富营养化和含磷固体废弃物的核辐射污染问题，磷资源的可持续利用及其环境影响已成为国内外学者广泛关注的重要科学问题。

物质流分析(MFA)通过定量分析特定时间和空间范围内物质(或元素)的迁移转化路径，识别其循环流动特征和回收利用的路径，定量分析人类社会经济系统与自然环境之间的物质交换，测度物质使用的环境影响，揭示不同时间和空间尺度内资源的流动特征和转换效率，可以为资源的高效利用和管理提供定量的决策信息，是经济、产业和资源管理等部门可持续发展评估相关研究中的重要分析工具之一^{[9, 10, 11]}。为了应对磷短缺、促进磷资源的可持续利用和管理，国内外广泛开展了针对磷元素的物质流分析，以分析不同尺度磷元素从开采、生产、加工、消费到最终处置、排放环节的环境影响和资源利用效率，促进农业和环境的可持续发展。本研究综述了磷元素物质流分析的最新研究进展，分析了现有磷元素物质流研究的特点和不足，以展望未来磷元素物质流分析及其相关研究的发展趋势，为国内可持续的磷资源管理提供方法学和信息基础。

物质流分析(MFA)指按照质量守恒定律，以实物的物理质量为单位，在一定的时间范围和空间边界内，对预先定义的社会经济系统进行的从物质采掘、生产、转换、消费、循环使用直至最终处置过程中流量和存量的系统性定量评估，它从早期的社会代谢、工业代谢的相关研究发展而来，已成为包括过程控制、资源管理、废物处理、环境管理、产品设计及生命周期评价等各个领域最常用的决策支持工具之一^{[9, 10, 11]}。当物质流分析的对象为特定元素(例如氮、磷、硫、铝等等)或者单一物质时，也称为元素流分析(SFA)，目前元素周期表中已经有约95种元素在不同尺度上进行了元素流分析，其中进行物质流分析的大多数元素为金属元素^{[12, 13]}。

物质流分析模型主要包括只考虑流量变化的静态模型和考虑存量变化的动态模型两大类，其中前者主要用于判别自然环境退化和污染物的来源、预测不同社会经济因素导致的环境影响，后者主要用于分析支撑社会经济发展的资源总需求、发展模式和长期变化趋势^{[11, 13]}。随着相关研究的不断深入，除了传统的包括总物质需求和输出(TMRO)和元素流分析(SFA)分析方法，物质流分析现在已经衍生出或者与多种研究手段和工具结合，包括实物投入产出表(PIOT)、生态足迹分析(EFA)、网络分析、和生命周期分析(LCA)、局部均衡模型等等，其分析问题的深度广度和对决策政策的支持功能也因此在逐渐强化^[14]。

由于其对人类社会的重要支持作用以及利用过程中的重要环境影响，磷元素一直物质流分析相关研究关注的热点元素之一^[12]。目前世界各国在不同尺度、针对不同问题开展了大量磷元素的物质流分析^[15]，根据研究尺度的不同，现有磷元素的物质流分析可以划分全球尺度、国家尺度、区域与城市尺度、企业及产品尺度四个层次，其中国家尺度的研究是目前磷物质流分析研究的主流。

从全球尺度看，现有研究主要集中于核算人为活动导致的磷流动，以及这些人为活动导致的磷流动对磷资源的耗竭速度和环境影响的加速和扰动效应(表1)。Smil^[16]最早分析了全球磷的自然及人为循环，认为化肥的施用、生活和工业废水的排放、持续的土壤腐蚀与流失以及作物秸秆和粪便的排放和利用加速了全球人为磷代谢(P metabolism)强度，并估算到2000 年全球人为活动导致的磷代谢强度将是其自然循环过程的3倍。Villalba等^[17]和Liu等^[18]分别从生产和消费的角度量化了全球磷代谢过程，提出全球开采80%的磷资源用于生产磷酸和生产磷肥，但是全球耕地总体仍然面临着磷亏损问题。Cornell等^{[4, 19, 20]}分析了全球食品生产和消费体系的磷流动，比较了不同饮食偏好(diet preference)的磷需求，断言全球的磷峰值(P peak)将于2030 年左右达到，现有磷资源仅能满足全球50—100 年的需求，世界各国将普遍面临着“磷短缺(P scarcity)”问题。van Vuuren等其它学者^{[21, 22, 23]}却认为现有磷资源至少可满足21世纪的全球磷需求，全球磷生产的峰值可能出现在21世纪末期，到24世纪缓慢下降60%。Scholz^[24]认为现有评估都基于静态指标，如静态使用寿命(Static Lifetime)、哈伯特曲线(Hubbert Curve)及赫芬达尔-赫希曼指数(Herfindahl-Hirschman Index)等等，而这些静态指标实际上无法预测磷资源的物理性稀缺，因此必须构建考虑技术变革、生产链潜力、资源耗散特征、存量变化、效率提升和循环利用的磷流动的动态模型来分析长期磷稀缺问题。Koppelaar和Weikard^[25]利用2009年的全球磷流动模型和局部均衡模型模拟分析了不同磷循环的技术和措施对全球磷资源耗竭速度的影响，结果表明现有的磷矿资源可以维持至少到22世纪磷资源的高消费使用，循环利用和减少磷资源使用的技术措施将大幅提高磷资源储量的寿命，并将磷资源(包括潜在的磷资源储量)耗竭时间延长至23世纪以后。

从国家尺度看，目前已经有10余个国家(包括跨国区域)开展了磷流动研究，包括欧盟^[26]、美国^{[27, 28]}、澳大利亚^[29]、英国^[30]、法国^{[31, 32]}、芬兰^{[33, 34, 35]}、瑞典^{[36, 37]}、荷兰^[38]、奥地利^{[12, 39]}、土耳其^[12]、日本^{[40, 41]}、韩国^[42]、马来西亚^[43]和中国^{[44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56]}(表2)。从研究问题看，国家尺度的研究更多关注农业行业或者食品生产和消费活动对磷流动的扰动及其造成的环境影响，例如Suh和Yee^[27]、Cordell 等^[29]、Senthilkumar 等^{[31, 32]}、Chen 等^[45]等。其他的热点问题还包括社会经济系统的内涵(embodied P demand)、虚拟磷需求(virtual P requirement)和磷足迹(P footprint)^{[28, 37, 41]}；以及其他部门，例如林业部门^[33]、能源部门^[35]和钢铁行业^{[40, 42]}的磷流动特征以及不同社会经济驱动力(经济发展、人口增长、饮食偏好的变化等等)对磷代谢的影响^[44]。从模型的特征看，国家尺度的磷流动模型以静态模型为主，但是也有一些模型开始考虑磷的相关存量变化和动态特征^{[12, 39, 55]}。

区域和城市尺度的研究包括在流域、城市和小区域(不跨国，例如省、县、区等等)水平进行的磷元素物质流动研究。与全球和国家层次的研究相比，该层次的研究比较少见。目前开展了磷流动研究的流域包括中国的滇池流域^{[57, 58]}，巢湖流域^{[59, 60, 61, 62, 63]}和北京密云水库^[64]；城市包括美国凤凰城^[65]，澳大利亚悉尼市^[66]、瑞典Linkoping市^{[67, 68]}，瑞典的Gothenburg市^[69]，越南的Haiphong市^[70]，中国的北京^[71]、天津^[71]和合肥市^[72]。另外Wu 等^[73]还针对安徽省开展了省水平磷利用效率的研究。从模型特征来看，现有区域和城市尺度的磷流动研究相对简单，以静态模型为主，主要关注农业和食品消费对水体的影响，尤其是对水体富营养化的影响。

企业和产品尺度即针对具体企业、产品或者技术系统的磷物质流研究，目前相关研究很少，但是这类能够传递的信息也最为具体，比较成熟的研究多与生命周期分析结合。例如，Bai 等^[74]研究了2010年典型城市排水系统不同技术选择对磷流动的影响，Asmala和Saikku研究了^[75]虹鳟鱼生产和消费体系的磷流动，分析水产养殖业对水体富营养化的贡献。Ridoutt 等^[76]评估了具有磷高效利用特征的水稻系统的环境影响，表明该水稻可以减少68%的磷肥投入并大大减少磷肥向水体的流失。

一般而言，现有研究将磷元素的流动途径区分为自然流动和人为流动，其中自然流动主要考虑大气沉降、侵蚀、径流和淋失等过程，人为流动则按照研究尺度、研究问题、研究目的的不同区别为不同的子系统之间的投入和支出。然而，尽管研究尺度、问题和目的有所不同，现有磷的物质流动模型的结构有着很大的相似性，其人为流动模块包括的子系统不外乎磷矿开采、磷肥生产、农业生产(作物种植、畜禽养殖、渔业、林业等)、其他工业生产(例如洗涤剂、食品添加剂、饲料、饲料添加剂)、居民消费(包括食品和其他相关磷产品)以及废水废物处理处置和排放。由于更多关注人为活动对磷元素流动的扰动及其相应带来的环境影响，现有研究对自然活动的磷流动过程(例如风蚀、水蚀、大气沉降、地表挥发、底泥释放等)的关注和描述较少，更无法考虑其他元素流动对磷元素流动的影响。因此，自然活动和人为活动的耦合研究，尤其是磷元素与其他重要元素的耦合研究和比较将成为今后研究的重点。

从研究方法来看，现有的模型多数都是静态或者比较静态模型，同时考虑流量和存量变化的动态研究非常少，与产业生态学中的其他方法，例如投入产出、生命周期分析和生态足迹等，结合也较少，但是模型的动态化和与其他方法结合的相关研究正在增加。例如Suh和Yee^[27]，Matsubae-Yokoyama 等^[40] 分别将物质流分析与生命周期分析、投入产出结合，分析了美国和日本食品生产和消费过程的磷利用效率和回收潜力。Metson等^[77]和Straat^[78]分别针对全球和瑞典养殖业构建了磷足迹模型。因此，今后磷的物质流分析必须更多与产业经济学及其他领域的各种模型和方法结合。例如经济模型、生命周期分析、投入产出分析、生态足迹分析、网络分析、可计算的一般均衡模型、个体行为模型等等结合，发展考虑存量变化的动态模型工具，以预测全球变化、社会经济发展、技术变革、政策发展等长期变化对磷流动的扰动及其相应的环境影响。

从研究问题看，现有研究主要关注磷资源管理的六大问题：(1)磷资源需求峰值、耗竭时间(depletion time)和对社会经济发展的关键作用；(2)农业、食品生产和消费对磷元素循环的扰动及其环境影响；(3)社会经济系统的磷资源利用效率和总需求；(4)含磷废物废水的排放及其环境影响，尤其是对地表水体富营养化的影响；(5)含磷废物废水的处理处置和回收潜力；以及(6)不同的管理和技术选择对磷资源可持续利用及其向环境排放的影响。其中研究最多的是第二类和第三类问题，研究最多的经济部门为农业部门。然而，由于不同国家和区域具有不同的磷资源禀赋，其关注的部门和磷管理问题也有很大差异。例如日本和韩国作为磷资源的进口国，格外关注从不同途径回收和利用各种磷资源，并开展了钢铁炉渣磷回收的潜力研究，结果表明钢铁炉渣中的磷存量可以达到两个国家磷消费量的10%以上，是两国最具回收潜力的磷资源^{[40, 42]}。随着生物燃料的发展加速，美国^[28]和瑞典^[36] 在磷流动模型中区分生物燃料部门，以特别考虑生物燃料发展对磷流动的影响。但是在全球“磷短缺”和环境剧烈变化的背景下，如何分析磷物质流动的长期变化、不同社会驱动力对磷需求的影响以及全球变化对磷流量和存量的影响变得尤为重要。

综上所述，磷物质流分析及其相关研究未来关注的主要问题包括五大方面：(1)考虑社会驱动力(人口增长、经济发展、技术变革和消费偏好)和存量变化的磷物质流动的长期动态模拟，以分析人类社会磷素需求的峰值和拐点及其相应的环境影响。(2)不同层次的磷足迹研究，包括社会经济发展的磷足迹以及不同产品的磷足迹^[78]，以分析不同社会经济发展模式和不同产品对磷素总需求的影响，从而能够从磷素管理的角度分析不同经济发展模式或者不同产品的优先性。(3)磷元素与其他元素和物质相比(包括氮、碳、金属、水等)，对于人类社会经济发展的重要性^{[79, 80]}，即分析在现有的需求和资源储量格局下，哪一种元素的供给必须优先保障，一般而言元素的不可替代的作用越强，其供给越应该得到优先保障。(4)全球变化(例如气候变化)背景下，食品体系磷依赖的脆弱性^[1]，即分析全球变化是否会影响食品体系磷需求及其环境影响，例如温度升高是否会降低磷肥的肥效从而导致磷肥的用量增加？降雨空间格局和时间格局(极端降雨的增加)的变化是否会增加磷的流失？(5)磷和其他元素的耦合流动研究，包括其他的重要营养元素(氮、碳、钾、硫)和金属元素(铁、铀^[81]、稀土金属等)，以分析不同政策和战略对多种污染物或者多种资源管理的共生效益或者权衡。