基于氮磷农田利用的黄淮海地区畜禽粪尿土地承载力研究

为评估黄淮海地区畜禽粪尿土地承载力,以《畜禽养殖业源产排污系数手册》为基础,采用2018年国家及地方统计年鉴数据,测算了黄淮海地区畜禽粪尿总量、耕地畜禽粪尿氮磷负荷和畜禽养殖环境风险指数,并基于氮磷农田利用估算了该地区畜禽粪尿土地承载力及发展潜力。结果表明:(1)2017年黄淮海地区畜禽养殖的粪尿、氮和磷产生量分别为31 434.18万、172.90万和31.10万t,单位面积耕地的畜禽粪尿、氮和磷负荷分别为12.00 t·hm~(-2)、66.05 kg·hm~(-2)和11.88 kg·hm~(-2);(2)畜禽污染物耕地负荷总体上处于较低水平,81.82%的城市畜禽粪尿耕地负荷低于15 t·hm~(-2),70.91%的城市畜禽氮负荷低于75 kg·hm~(-2),78.18%的城市畜禽磷负荷低于15 kg·hm~(-2);(3)从畜禽养殖环境风险看,北京属削减养殖区,天津、秦皇岛、烟台属严格控制区,威海、日照、德州、滨州、洛阳、许昌、漯河属适当控制区,邢台、安阳、商丘、枣庄、亳州、淮北、徐州、连云港、盐城、淮安属重点发展区,其余地区属一般发展区;(4)黄淮海地区基于氮磷农田利用的土地承载力为42 262.70万头猪当量,畜禽养殖潜力值为26 469.82万头猪当量,其中北京需削减66.34万头猪当量,天津、安徽、江苏、河北、河南和山东畜禽发展潜力分别为104.56万、3 569.14万、3 691.77万、4 189.90万、8 371.26万和6 609.53万头猪当量。     

北京密云水库上游太师屯镇非点源污染损失估算

选取密云水库上游太师屯镇为研究区,在利用输出系数法计算非点源污染负荷的基础上,运用环境经济学方法,定量估算了非点源污染损失.研究结果表明:本区养殖污染负荷最高,总氮占全镇总氮污染负荷的53.0%.其中规模化养殖占畜禽排泄污染负荷的66.5%;总磷占全镇总磷污染负荷的73.9%,其中规模化养殖占畜禽排泄污染负荷的66.5%.农村生活污染的总氮负荷占全镇总氮污染负荷的31.3%,总磷占全镇总磷污染负荷的15.8%.土地利用污染负荷最低,总磷和总氮分别占全镇总负荷的15.7%和10.3%.区域非点源污染损失折合人民币约407.57万元,其中畜禽养殖污染损失最大,约占总损失的40.7%;土地利用流失次之,占33.5%,其中63.2%的损失是耕地污染造成的;农村生活污染损失最小,但总氮损失达总损失的25%.氮、磷养分流失是主要的非点源污染形式,损失约占79.5%.土地废弃和泥沙淤积造成的污染对水库存在长期影响,也不可忽视.     

珠江三角洲典型地物空间分布特征研究

珠江三角洲改革开放30年以来经历了巨大的发展,但也导致了城市的迅速扩张以及大量宝贵农田的流失,引起了诸多环境问题和社会问题。因此探讨珠江三角洲土地利用的演变规律对于未来土地利用的发展变化具有重要的借鉴和指导意义。空间分析因其能有效揭示隐藏在区域数据之后的规律信息以及提供直观的可视化信息表达,已引起许多学者的关注。引入空间整体统计参数中的平均中心、标准距离和标准偏移椭圆,探讨珠江三角洲耕地和建设用地两种典型地物1990—2005年间空间整体特征变化规律。通过对研究区及研究区各市(区)空间中心转移、空间分布扩散和空间分布偏移的计算,表明珠江三角洲耕地空间特征在1990—2005整个阶段均呈现扩散状态,而建设用地空间特征在1990—2000阶段和2000—2005阶段呈现收缩和扩散两种特征。两类地物空间特征变化强度2000—2005阶段均大于1990—2000阶段。典型地物空间特征的发展变化与社会经济发展相一致。通过对比发现,耕地与建设用地空间中心迁移呈现反向变化,揭示了经济发展中大量的耕地转化为经济价值更高的建设用地。研究区中各市(区)典型地物在不同阶段、不同区域土地利用发展也存在空间差异性。空间统计指标能直观有效地揭示地物空间分布格局变化规律,有利于深入探讨研究区土地利用/覆盖变化。     

干旱区绿洲城市生态系统服务价值空间自相关格局分析与模拟

基于1995—2015年乌鲁木齐土地利用遥感解译数据,运用CA-Markov模型模拟预测2025年土地利用变化并进行生态系统服务价值(ESV)评估,应用格网分析、空间自相关、Getis-Ord热点区分析等空间统计方法,进一步分析生态系统服务价值与土地利用空间格局的分布特征,讨论两者空间自相关关系及年际变化的冷热点空间分布变化原因。结果表明,近30 a研究区各土地利用类型的生态系统服务价值呈波动变化态势,但总体上表现为减少趋势,主要原因在于耕地、林草地提供的生态系统服务价值减少。生态系统服务价值空间分布格局差异明显,高值区域多年来主要分布在乌鲁木齐县林地面积较多的地区,低值区域集中分布于乌鲁木齐中心四区(天山区、沙依巴克区、水磨沟区、头屯河区)的周边区域。研究区生态系统服务价值具有显著的空间正自相关性与空间聚集特征,高-高、低-低聚集分布分别与生态系统服务价值高值、低值区域高度重合。生态系统服务价值热冷点区呈逐年增加趋势,热点区增加与近年来乌鲁木齐大力推行生态保护及修复等措施有极大关联,冷点区增加与城市建成区无序扩张及土地开发等一系列人类活动有密切关系。     

云南不同类型农村沟渠底泥中氮赋存形态分布

由于人为干扰强度较大,农村沟渠底泥氮形态的分布与河流、湖泊等存在较大差异,研究农村沟渠底泥氮形态的分布特征具有重要意义.本研究以云南省4种典型农村(畜禽养殖型、生态休闲型、集镇型和传统型)为对象,研究其沟渠底泥氮形态的分布特征.结果表明,全氮(TN)含量为2285.26—5713.51 mg·kg-1,空间分布呈畜禽养殖型集镇型生态休闲型传统型;各形态可转化态氮(TTN)的含量由大到小依次表现为强氧化剂可提取态氮(SOEF-N)强碱可提取态氮(WAEF-N)≈弱酸可提取态氮(SAEF-N)离子交换态氮(IEF-N);释放程度从易到难依次为IEF-NWAEF-NSAEF-NSOEF-N;除IEF-N外,各形态氮的含量最大值均出现在畜禽养殖型农村.方差分析表明,w(TN)在集合点和分散点间存在显著差异性;传统型农村的w(IEF-N)、w(WAEF-N)明显低于其余3个类型;畜禽养殖型农村的w(SAEF-N)、w(SOEF-N)明显高于其余3个类型.对农村环境氮污染风险贡献最大的是IEF-N和SOEF-N,其次为WAEF-N和SAEF-N.畜禽养殖型农村氮污染风险最高,其次是集镇型和生态休闲型,传统型最低.     

基于日平均养殖量的畜禽养殖排污系数法的改进

传统的畜禽养殖排污系数污染负荷计算方法受养殖规模浮动影响,计算结果精度不高,只能计算整年污染负荷且无法计算某时段污染负荷。在分析传统畜禽产污系数法缺陷的基础上,探讨了日平均养殖量取代年养殖量的计算模式。依据不同类型畜禽养殖的周期特点,提出了各种类型畜禽养殖排污系数污染负荷计算公式和方法。最后,利用2008—2009年天津市宁河县生猪养殖污染负荷计算实例,对比分析了两种计算方法。结果表明,基于日平均养殖量的产污系数计算方法可以有效地减少由于养殖规模波动造成的计算误差,同时还适用于年内某时段污染负荷的计算,克服了传统方法只能以年为计算基本单元的弊端。     

鄱阳湖生态经济区畜禽养殖污染负荷分析

以2011年鄱阳湖生态经济区各类畜禽养殖数量为基础,采用参数估算法估算了该区畜禽养殖产生的粪便及其中的污染物化学需氧量(CODCr)、氨氮(NH3-N)、总氮(TN)和总磷(TP)质量,结合各县(市、区)耕地面积计算得到畜禽粪便耕地负荷量,并以畜禽粪便耕地负荷警报值评价各县(市、区)畜禽粪便负荷量的承受程度。结果显示,2011年该经济区畜禽养殖共产生3 129.42万t粪尿,83.34万t CODCr,6.97万t NH3-N,17.19万t TN,5.79万t TP,该区这些污染物产生量约占全省的40%。对畜禽粪便耕地负荷警报值的分析表明,该经济区有9个县(市、区)畜禽粪便负荷量尚未对环境构成潜在威胁,其余26个县(市、区)畜禽粪便负荷量从理论上来讲都对环境构成了潜在威胁,尤以南昌市青云谱区和九江市浔阳区畜禽粪便负荷量对环境的污染威胁最大。     

兰州—西宁城市群城乡建设用地开发强度空间格局变化及影响因素分析

明确土地开发强度空间格局特征和演变规律,对城市发展边界、生态安全红线和耕地保护红线划定具有重要的指导意义。以1995、2005、2015和2018年兰州—西宁(兰西)城市群39个县级单元为研究区,运用趋势面分析、空间自相关和地理探测器等方法,分析兰西城市群城乡建设用地开发强度空间格局特征及其影响因素。结果表明:(1)兰西城市群城乡建设用地开发强度均值总体呈增长态势,由1995年的2.74%增长到2018年的4.45%,在空间上呈现兰州—西宁2个中心高、外围地区依次降低的"核心—外围"空间分布格局。(2)东西方向和南北方向空间差异较大,且呈现出扁平的倒"U"型曲线分布态势。(3)城乡建设用地开发强度空间上存在显著正相关性且集聚特征明显,冷热点区的空间分布具有相对稳定性。(4)1995—2018年城乡建设用地开发强度受人口集聚水平、经济发展水平、产业结构水平和海拔等要素综合影响。研究结果可为兰西城市群以及同类型地区城乡建设用地管控和高质量发展提供决策参考。     

中国焦化场地近20年时空演变特征及驱动因素

落后的焦化工艺产生的粉尘、烟气和化学废气进入空气,可对场地周边及较远区域造成严重的大气污染;剩余氨水、粗苯分离水等废水排入水体或农田会造成区域水体与土壤污染;废水泄露和固体废物(如焦油渣)管理处置不当,污染物直接进入土壤也会造成场地土壤污染。文章通过国家企业信用信息公示系统、天眼查和文献数据库等共享数据平台收集中国存在较大潜在污染的焦化场地样本信息,采用空间单元统计、标准差椭圆、最邻近指数和核密度估算等空间分析方法研究了2000—2020年中国焦化场地时空格局与演变特征,采用多因素时空关联分析方法探讨了其驱动因素。结果显示,中国焦化场地总体分布极不均衡,主要分布在华北、西南和西北地区,东南地区分布较少;近20年来,在数量特征上总体呈现先快速增加(2000—2005)后缓慢增加(2005—2015)再缓慢减少(2015—2020)的变化趋势;在空间格局上,分布走向有从西南—东北方向往西—东方向变化的趋势,分布重心有向西北方向移动的趋势,集聚程度逐渐增强,呈现显著集聚性,形成了华北、西北新疆以及西南云贵等高集聚区。这种时空演变特征是自然地理因素、社会经济发展、国家行业政策及公众环保意识等多因素共同作用的结果,其中自然地理因素决定了中国焦化场地的基本空间格局,而社会经济发展、国家行业政策和公众环保意识是其空间格局随时间不断演变的主要推动因素。本研究结果可为中国土壤、水和大气环境的空间治理和生态环境建设的宏观战略决策提供科学依据,也可为焦化及相关产业的发展规划提供信息支持。     

东北典型黑土区农村居民点对水土流失影响研究

从人地关系视角探究典型区域农村居民点布局对生态环境的影响,对合理优化农村居民点分布格局、保障粮食安全具有重要战略意义。以东北典型黑土区克山县为例,综合运用分散度统计、核密度估计、空间关联测度模型等空间分析方法,揭示1988—2018年克山县农村居民点时空演变规律;结合修正通用RULSE模型,采用景观指数分析和广义变异系数的灰色关联度量化模型,探讨农村居民点分布特征对水土流失的影响。结果表明:(1)近30 a来,克山县农村居民点数量急剧增加,分散化趋势明显,以高密度多核心集聚、低密度广泛分布为主要空间分布特征,分布重心向东北部迁移。(2)研究区水土流失以微度、轻度和中度侵蚀为主,水土流失强度空间分布格局与农村居民点的空间分布格局呈现一致性。(3)研究区农村居民点分布特征因子与水土流失面积的关系存在差异,农村居民点分布范围与密度对水土流失具有较强作用。     

太湖流域跨界区农业面源污染特征

采用输出系数法,从农村生活污水、农田径流和畜禽养殖3个来源估算太湖流域跨界区农业面源COD、TN和TP污染负荷,并分析其时空变化特征;采用极值标准化法和均方差赋权法,进一步对区域内各地区污染负荷的入水强度进行评价。结果表明,太湖流域跨界区2010—2013年COD、TN和TP入水负荷总体呈下降趋势,4 a间高COD负荷区分布在湖州市和嘉兴市,高TN和TP负荷区分布在苏州市、湖州市和嘉兴市,畜禽养殖是3种污染负荷的主要来源,均占3种污染物总负荷量的80%以上。各地区污染负荷的入水强度评价结果显示,2010—2013年跨界区农业面源污染负荷单位GDP综合评价值分别为0.217 6、0.240 2、0.236 4和0.222 2,无明显变化规律;流域农业面源污染负荷单位面积综合评价值分别为0.407 8、0.289 9、0.289 2和0.281 6,呈下降趋势。     

村级养殖种植园区碳素物质流分析——以北京市平谷区西柏店村为例

将西柏店村畜禽养殖规模折合为1.5万头猪场当量污染负荷,并将整个园区生产工艺分为养殖、废弃物处理和种植3个阶段,不考虑隐藏流的情况下,以1年为系统边界,通过数据调查、已有资料研究和小区种植试验,采用物质流分析方法分析了西柏店村养殖种植园区在整个生产工艺的碳素流动,以期为村级养殖种植园区大力发展低碳经济提供新的方法和视角,为村级区域循环经济及可持续发展提供减少环境压力解决方案的科学依据。通过园区养殖种植过程的C素分析表明,养殖阶段年输入C素总量为112.52×10^4 kg,其中猪身总固碳量为40.04×10^4 kg,粪碳和尿碳总量为49.29×10^4 kg,以CO2形式代谢排出的C为23.19×10^4 kg。废弃物处理阶段输入的碳主要为粪碳和尿碳,其总量为49.29×10^4 kg,其中9.79×10^4 kg尿碳直接进入种植阶段,39.50×10^4 kg粪碳进入沼气站处理,沼气转化出的碳为11.02×10^4 kg,其中CH4为8.43×10^4 kg,CO2为2.59×10^4 kg,养殖污水中通过CH4排放再加上其他途径释放的碳约有23.66×10^4 kg,占粪碳量的59.89%。进入种植阶段的碳素主要为尿碳、沼渣和沼液的碳素,合计为14.61×10^4 kg,假设该村43 hm^2耕地能全部施用沼肥,不计其他作物种植,1季玉米种植土壤可库存有机碳为60.50×10^4 kg,为进入种植阶段碳素14.61×10^4 kg的4倍,还可增加植物有机碳27.31×10^4 kg。由C素流动分析可知,西柏店村具有可容纳该村养殖废弃物的环境容量,有较好实现养殖废弃物循环利用的条件,但需大力加强畜禽废弃物的管理和处理,提高园区养殖废弃物循环利用效率。     

近30年来扎龙湿地丹顶鹤繁殖种群空间格局及动态

种群空间格局是种群生态学和空间生态学的核心问题之一。采集扎龙湿地1981—2008年丹顶鹤种群巢址定位数据,采用平均中心模型和最近邻体分析等方法研究了该研究区丹顶鹤种群空间分布格局及其动态,探讨了变化的原因与机制。结果表明:近30 a间,空间格局总体分布由研究区域核心区中部转移至西南部,再转移至东北区域,丹顶鹤个体间的离散程度经历了逐渐增大→逐渐减小→迅速增大→逐渐减小的过程,各时期空间分布型均为聚集分布,但不同时期聚集强度差异明显。生境条件的剧烈变化,即湿地缺水加剧以及人类活动干扰(如湿地内修建的各种大型工程)是导致丹顶鹤空间格局变化的主要原因,突发性干扰(如火灾)对空间格局也产生了极大的影响。     

景观格局演变及其生态效应研究进展

区域景观格局演变及生态效应研究一直是景观生态学研究的重点,对人类、自然的可持续发展具有重要的意义.文章从景观格局-过程-生态效应的相互驱动机制及当前的研究进展进行了系统的分析,综述了国内外景观格局及其生态效应研究相关进展,集中分析了当前研究中的热点和难点.景观分类是景观格局研究的基础,空间统计分析、转移矩阵分析和景观指数分析是研究景观格局演变的主要手段,新兴的元胞自动机理论和方法在景观格局演变研究中有很大的发展空间.景观格局演变驱动力研究需对自然驱动力、社会驱动力进行综合考虑.不同时空尺度下,景观格局与生态过程的内在机制也不一样.在此基础上,提出了今后研究的一些方向,重点是:(1)景观格局演变的内在驱动机制研究;(2)多尺度下景观格局和生态过程(如土壤侵蚀、养分流失等)的耦合研究;(3)"3S"技术、野外调查监测与格局-过程演变模拟模型结合的应用推广潜力.     

基于县域尺度的稻田土壤碱解氮空间异质性研究

掌握土壤氮素的空间分异特征是实现养分优化管理的重要基础。以江苏省扬中市为研究区,基于田间采样的99个土壤碱解氮样本数据,运用地统计学与缓冲区分析方法,研究县域稻田土壤碱解氮的空间变异特征及其影响因素。结果表明,研究区土壤碱解氮平均含量为114.30 mg·kg~(-1),处于较丰富水平。半方差分析结果显示,土壤碱解氮空间分布具有中等程度的空间自相关性。空间插值结果显示,扬中市土壤碱解氮空间分布总体上呈斑块状,大致表现为中部高、两端低,高值区主要位于研究区中部3镇,在工业区与养殖基地附近均出现局部高值区。影响因素分析表明,城镇化、土壤质地、产业分布、农户经营等因子对研究区土壤碱解氮空间分布均具有较大影响。距离养殖合作社1 km以内区域的土壤碱解氮易于积累,工厂、企业对周边600 m范围内的土壤碱解氮含量影响显著(P0.05)。     

鄱阳湖生态经济区土地利用转型热点区域识别及其动力机制研究

识别土地利用转型的热点区域并探究其动力机制可为保护土地生态环境安全提供科技支撑,促进区域经济社会的可持续发展。利用空间热点分析模型探测鄱阳湖生态经济区2000—2014年期间土地利用转型的热点区域,通过典型相关分析定量识别土地利用转型的主导驱动因子,并从供给与需求的经济学视角阐释了土地利用转型的驱动机制。研究表明,(1)滨湖地区2000—2014年发生了剧烈的土地利用转型,耕地、林地、草地、滩地大量向城乡建设用地转化,土地生态面临严重威胁。在距鄱阳湖水面10 km之内,耕地、林地和草地面积分别减少15 740、15 152 hm~2,建设用地面积增加了47 138 hm~2。(2)城市建成区周围及主要道路对土地利用转型具有显著的辐射作用,区内土地利用变化高频区呈现"井"字形空间格局。在距离道路10 km之内,2000—2014年期间土地利用综合动态度由2.765‰衰减至1.059‰。(3)土地利用转型的热点区域数量呈增加趋势,并具有明显的轮换与跃迁特征,基本形成以南昌、九江、景德镇、鹰潭为中心的4个快速转型区。(4)经济、社会发展导致的土地供需矛盾是引发土地利用形态变化的主要动因。亟待建立滨湖生态环境缓冲区,严格控制滨湖地区的土地资源开发利用活动;对土地利用转型的热点区域进行科学规划,通过土地供需调控促进土地资源的集约、高效利用,增加土地经济供给。     

基于CA-Markov多情景模拟的海南岛土地利用变化及其生态环境效应

以海南岛为研究区域,基于2010年和2020年两期土地利用数据,利用CA-Markov模型预测了3种模拟情景下2030年和2040年土地利用空间分布格局,定量分析了海南岛2010—2040年土地利用变化及其生态环境效应。结果表明,(1)2010—2020年间,海南岛土地利用类型以林地为主,面积占比达62%以上,其次耕地。该时期土地利用类型转移的主要特征是林地、耕地、水域转为建设用地,建设用地面积增加483.55 km2,增长率高达53.27%。(2)利用CA-Markov模型模拟3种情景下海南岛2030年和2040年的建设用地均呈现较大幅度的增长,其中耕地保护情景模拟下的耕地面积最多,生态保护情景模拟下的林地和水域面积有所增加。(3)海南岛生态环境质量呈现出中部高、四周低的阶梯式空间分布格局,2010—2040年间海南岛生态环境质量逐期降低,由2010年的0.606 4下降至2040年的0.587 8。在3种模拟情景中,生态保护情景下的海南岛生态环境质量最高,其次是耕地保护情景,基准情景最差。(4)2010—2040年间,林地转为建设用地是导致海南岛生态环境质量恶化的主要类型和主导因素,耕地转为林地则是致使未来生态环境质量改善的主要土地利用变化类型,生态贡献率超过70%。总体而言,海南岛生态环境质量恶化的趋势大于生态环境质量改善的趋势。     

苏州花桥国际商务城生态用地时空动态特征

以快速城镇化的小尺度区域——苏州花桥国际商务城为研究区,基于高分辨率遥感影像和GIS技术,利用土地利用转移矩阵、重心计算、景观格局指数和平均等效面积法分别对研究区生态用地时空动态、景观格局动态和土地生态质量空间格局动态特征进行研究。结果表明:(1)2005—2015年,研究区建设用地在快速城镇化进程中迅速扩张,而生态用地总面积减少913.11 hm~2,同时生态用地重心逐渐向南迁移。(2)2005—2009年,研究区耕地优势度下降,其他生态用地优势度未出现较大波动;2009—2015年,随着政府生态建设力度加大,研究区斑块破碎化程度降低,景观复杂程度降低,除耕地外其他生态用地优势度迅速提升。(3)2005—2015年研究区平均等效面积由0.28增至0.31,总体生态质量得到改善;全区有57.01%的土地生态质量获得提升,且主要分布于南北两端吴淞江附近,仅有17.22%的土地生态质量下降。这表明近年来在城镇生态用地紧缺的情况下,研究区城镇生态建设对生态环境改善起到重要作用。建议研究区未来重点提升中部区域土地生态质量,进一步优化土地资源配置,以实现土地生态质量全面提升。     

洞庭湖区景观生态风险评价及其时空演化

洞庭湖是长江流域重要的调蓄湖泊,该区生态系统服务极为重要,也是中国重要的农产品提供区域。科学评估洞庭湖区景观生态风险,揭示其时空演化特征,对探究其生态格局形成过程与机理,降低生态风险胁迫具有重要意义。利用1990、2000、2010、2015年4期LandsatTM多光谱遥感影像,分析了1990—2015年洞庭湖区土地利用演变特征;基于景观生态学理论引入景观生态风险评价模型,将研究区划分为1 165个生态风险评价单元,基于GIS和地统计学,构建景观生态风险指数模型,对洞庭湖区土地利用动态变化下的景观生态风险时空分布特征及空间关联格局进行了研究。结果表明:(1)1990—2015年,洞庭湖区各土地利用类型均发生了显著变化,其中建设用地面积增加较大,年均增幅达3.61%,主要由林地和耕地的转入;(2)1990—2015年,洞庭湖区的景观生态风险状况以较低生态风险、中等生态风险为主,其面积占景观总面积的70%以上。但高生态风险区域呈增加趋势,面积占比由1990年的4.95%,增加到2015年的8.70%;与2010年相比,中等风险区面积由46.09%降到2015年的34%,这期间低生态风险区和较低生态风险区的面积明显增加;(3)研究区各时期景观生态风险呈现正的空间相关性,且在空间上趋于集聚。景观生态风险分布与区域地理环境具有较高的耦合性,人类活动的集聚程度与风险的空间集聚程度是对应的。以湖滨地带为核心区和远离湖滨区为两大群体,高值-高值(H-H)风险聚集区主要分布湖滨周边,低值-低值(L-L)风险聚集区主要分布于远离湿地核心区的县市区域。研究结论以期为洞庭湖区生态保护与环境整治,调控生态风险提供参考。     

九龙江流域规模化养殖环境风险评价

在GIS支持下采用养分收支平衡法对九龙江流域规模化生猪养殖的环境风险进行评价。规模化养殖污染负荷在不同类型养殖场的分配、空间分布规律以及规模化养殖的环境风险评价结果表明：规模化养猪业年产粪尿量达107万t，以小型规模化养殖场产生的污染负荷最高；整体来看，规模化生猪养殖的环境风险不大，但由于养殖场的地理分布集中，造成近郊区域规模化养殖的环境风险极高；由于过多施用化肥，中远郊规模化养殖粪尿未得到合理利用而流失。位于水系附近的养殖场对水体影响较大，根据“福建省畜禽养殖污染防治管理办法”，为减少规模化生猪养殖对九龙江水质的直接威胁，需搬迁主要支流沿岸1km范围内的规模化养猪场250个。