基于文献调研的我国部分地区单位产量氮流失量影响因子研究

王 芊

(中国农业科学院 农业环境与可持续发展研究所/农业农村部农业环境重点实验室，北京 100081)

农业面源污染已成为世界范围内普遍关注的问题。在欧美等发达国家，农业面源污染成为河流污染的重要原因，美国59%的受污染河流源于农业源污染，欧洲河流氮污染中农业源占50%～70%[1]；我国太湖流域农业源污染占30%～40%，洱海流域、重庆三峡库区可达到68%以上，全国农业源氮、磷污染占50%以上[2-6]。2015年我国制定了“到2020年化肥使用量零增长行动方案”，大力推广测土配方施肥、秸秆还田、有机肥施用，降低化肥的施用量，以减轻对环境的污染[7]。长期以来单位面积污染负荷一直是衡量面源污染水平的主要指标，但由于该指标设定没有考虑到生产经济效益，不能实现环境效益与经济效益的统筹兼顾。近年来，已有研究提出了将氮流失量与产量相结合的单位产量氮流失量指标，作为指导不同管理措施开展的重要依据：如Zhou等[8]通过收集世界范围内种植玉米、小麦的30余篇氮淋洗文献，分析了施氮量与单位产量氮淋洗量间的关系，发现当施氮量降低至100～200 kg/hm2时，单位产量氮淋洗量平均值为可由5.40～5.41 kg/t降低至3.2～4.5 kg/t；Zhao等[9]基于北美地区30余篇玉米氮淋洗试验结果，分析了单位产量氮淋洗量与氮盈余量、降雨量、土壤质地、耕作措施间的关系，发现氮盈余量、降雨量是影响单位产量氮淋洗量的显著因子，土壤质地具有一定的影响，而耕作措施的影响较小。

现有作物产量尺度的氮流失研究多是基于国外研究数据，针对我国不同区域和施肥管理措施下的开展的研究鲜有报道，并且主要为面积尺度下的氮、磷流失研究[10-17]。因此，为研究我国单位产量氮流失量的影响因素，本研究拟选取长江中下游、黄淮海和东北平原这三大我国粮食主产区，以大田生产中较为常见的“化肥”、“秸秆还田”、“缓控释肥”的3类施肥管理措施等为关键词，搜索CNKI、ScienceDirect等国内外主要数据库，收集并进行数据提取有关氮流失文献工作，并基于典型区域计算不同作物、管理措施间的单位产量氮流失量(以下简称新指标)，通过逐步回归方法建立新指标与影响因子间的线性模型，同时对不同区域、管理措施间的新指标差异进行比较，进而为制订区域尺度下农田管理措施提出建议，也为我国农田面源污染防控政策提供技术支撑和理论依据。

1 材料与方法

1.1 数据来源

基于CNKI、Science Direct、Web of Science等数据库开展面源污染文献收集，以“面源污染”为主题，以“N”、“淋洗”、“径流”、“化肥”、“缓控释肥”、“秸秆还田”等关键词分别进行搜索，论文发表年份为1995—2016年，共收集近200篇文献。调查作物主要为小麦、水稻、玉米。文献中涉及东北平原和黄淮海平原的氮流失以硝态氮为主，长江中下游平原氮流失主要为全氮。对于淋洗测定深度而言，黄淮海平原测定的土层较深，集中在100～200 cm深处。长江中下游平原和东北平原深度较浅，集中在50～100 cm深处。

提取的指标包括水肥管理措施、土壤养分、质地及产量等，主要包含：氮肥施用量(kg/hm2)、磷肥施用量(kg/hm2)、砂粒含量(%)、粉粒含量(%)、粘粒含量(%)、土壤全N(g/kg)、全P(g/kg)、碱解氮(mg/kg)、速效磷(mg/kg)、降雨量(mm)、灌溉量(mm)、施肥方式、作物产量(t/hm2)等。土壤质地按美国分类标准，即粘粒(<0.002 mm)、粉粒(0.002～0.050 mm)和砂粒(0.05～2.00 mm)。磷肥施用量以P2O5计，土壤养分数据取自0～20 cm耕层。对于无法直接获取的文献数据，利用Engauge Digitizer软件提取。

1.2 逐步回归方法

本研究所利用的新指标的计算公式如下：

NI=NL/Y

式中：NI为单位产量氮径流量或淋洗量，kg/t；NL为氮径流量或淋洗量，kg/hm2；Y为作物籽粒产量,t/hm2。

本研究利用SPSS 13.0软件开展逐步回归运算。由于不同文献中提供的指标的完整度不尽相同，为了确保在逐步回归建模型过程中尽可能多地纳入数据齐全的自变量，有相当一部分文献在筛选过程中被排除在建模数据之外，因而模型的建模样本量也相应减少。逐步回归计算结果中显示了各步骤的方差，前后步骤方差的差视为新引入的影响因子的方差，该方差与最终模型总方差之比即为相应影响因子的权重值。对于未采取逐步回归方法的个别后向模型则未计算权重。在构建逐步回归模型时为尽可能考虑更多因素，在前述收集的文献中进一步筛选出23篇开展大田试验的文献，各省份的文献数量如表1所示。文中箱式图用Sigmaplot 12.5(Systat Software Inc.)制作，利用Microsoft Excel 2003软件制作曲线拟合图和计算回归曲线。

表1 参与计算文献省份分布情况Table 1 Adopted literature distribution in three major plains

1.3 水氮管理方式

在长江中下游平原：1)水稻。主要采取了泡田和晒田交替的常规的灌溉方式，安徽巢湖采取了浅湿灌溉方式。氮肥按播前、苗期、分蘖期、拔节期、穗期施肥3～4次，巢湖和江苏宜兴均施用3次氮肥，施用比例各不相同。嘉定地区主要追施碳铵，而浙江杭州和宜兴均为尿素。水稻施用缓控释肥涉及江苏常熟和浙江杭州，施肥方式均为底施，灌溉方式为连续灌溉和节水灌溉2种；2)小麦。氮肥按播前、苗期、拔节期、穗期施用分2～3次施用，生育期一般不灌溉。

在黄淮海平原：1)小麦。在播种、分蘖、起身、拔节、穗期等时期灌溉2～4次，每次灌溉量为75mm。氮肥按底肥、拔节肥施用2次；2)玉米。按出苗期和拔节期灌溉1～2次，氮肥按底肥、追肥施用1～2次；3)秸秆还田涉及河北清苑、栾城与山东南四湖区，水肥管理方式与常规方式相仿。

在东北平原：1)玉米。雨养种植，施肥方式多为底肥，玉米生长季降雨量为280～608 mm；2)水稻。基于黑龙江省齐齐哈尔市的试验，灌溉方式设置了常规灌溉、控制灌溉和浅湿灌溉3 种方式，氮肥按播种、返青期2 个时期各一半的比例施入。

2 结果与分析

2.1 长江中下游平原新指标影响因子

在长江中下游平原，氮流失新指标中以缓控释肥条件下的水稻氮淋洗值最低(0.77 kg/t)。施用化肥条件下，氮流失新指标按从大到小的顺序为小麦径流、水稻径流、小麦淋洗、水稻淋洗，相应的数值分别为4.35、3.07、2.26、0.89 kg/t(图1)。从氮流失量与作物产量关系看，与上述5 类新指标对应的两者的抛物线模型均未达显著水平(P>0.05)(图2)。水稻的抛物线曲线在产量约为7 t/hm2对应的淋洗量和径流量值最高(图2(a)、(c)和(e))，以径流量值最高，约为25 kg/hm2(图2(c))。小麦抛物线曲线在产量处于4～5 t/hm2时对应的淋洗量和径流量值达到最高(图2(b)和(d))，小麦淋洗最高值较水稻值高(图2(a)和(b))，最高值分别约为13和8 kg/hm2。

箱式图从下到上的横实线分别代表10%、25%、50%、75%、90%的分位线，位于分位线之外的点为所对应分布函数<10%或>90%的点。虚线为平均值。RL, RR, WL, WR分别代表施用化肥条件下水稻氮淋洗、水稻氮径流、小麦氮淋洗、小麦氮径流。RSCRFL为施用缓控释肥条件下的水稻氮淋洗。n为样本量。图3和图5同。Solid horizontal lines of the box plot from the bottom to the top represent 10%, 25%, 50%, 75% and 90% percentiles, the dots outside the 10% and 90% percentiles represent data lower than 10% percentile or higher than 90% percentile. Dotted lines represent mean values. RL, RR, WL, WR respectively represent N leaching of rice, N runoff loss of rice, N leaching of wheat, N runoff of wheat after applying chemical fertilizer. RSCRFL represents N leaching of rice by applying slow and controlled-release fertilizer. n represents number of samples. The same in Fig.3 and Fig.5.图1 长江中下游平原单位产量氮流失量(新指标)Fig.1 Yield-scaled N loss new index in Middle-Lower Yangtze Plain

在长江中下游平原，施用化肥条件下水稻氮径流、淋洗新指标模型分别达到了显著、极显著水平，模型自变量分别为施用P2O5量、土壤全N，权重值分别为0.15、0.49。降雨量与2 个因变量间的相关系数均未达到显著水平(P>0.05)(表2)。小麦氮径流新指标模型达到极显著水平(P<0.01)，自变量为土壤全N含量和降雨量，权重值分别为0.52、0.41。小麦氮淋洗新指标模型未达显著水平(P>0.05)，从相关系数看降雨量的影响相对较高，且呈正相关(数据未列出)。缓控释肥条件下，水稻氮淋洗量新指标模型的自变量仅为土壤全N含量，模型达到极显著水平，且与土壤全N含量间呈正相关(表3)。从水肥管理措施看，水稻管理存在着较大差异，小麦管理措施基本相同。施用化肥条件下，水稻氮径流模型的区域间灌溉和施用微肥管理存在着一定差异(安徽巢湖施用了锌肥和采取了浅湿灌溉)，水稻氮淋洗模型主要是施用氮肥种类的差异(嘉定主要追施碳铵，而杭州和宜兴均为尿素)。缓控释肥条件下，主要差异为普通灌溉和节水灌溉(杭州采用了节水灌溉处理)。

表2 长江中下游平原施用化肥的新指标模型系数[18-26]Table 2 Yield-scaled N loss model by applying chemical fertilizer in Middle-Lower Yangtze Plain[18-26]

表3 长江中下游平原缓控释肥条件下水稻氮淋洗新指标模型系数[23,27-28]Table 3 Yield-scaled N leaching loss model by applying slow and controlled-release fertilizer in Middle-Lower Yangtze Plain[23,27-28]

(a)、(b)、(c)和(d)分别为施用化肥条件水稻氮淋洗量、小麦氮淋洗量、水稻氮径流量及小麦氮径流量，(e)为缓控释肥条件下水稻氮淋洗量。R2代表决定系数，ns代表模型未达到显著水平，*、**代表模型分别达到P=0.05、P=0.01显著性水平。下同。(a), (b), (c) and (d) represent RL, WL, RR, WR and respectively. (e) represents RSCRFL.R2 represents determination coefficient. ns represents the model doesn’t reach 0.05 significant level. * and ** represent the model reaches P=0.05 and P=0.01 significance level, respectively. The same below.图2 长江中下游平原氮流失量与作物产量关系Fig.2 Relationships between N loss and crop yield in Middle-Lower Yangtze Plain

2.2 黄淮海平原新指标影响因子

在黄淮海平原，小麦玉米秸秆还田周年轮作的氮淋洗新指标值比仅施用化肥的小麦或玉米单价低，平均值为 2.35 kg/t，小麦、玉米值分别为5.87 和4.66 kg/t[28-36](图3)。从产量和氮淋洗量之间的变化趋势看，秸秆还田措施下两者呈指数曲线变化(图4(c))，仅施用化肥的玉米氮淋洗量与产量呈线性变化(图4(a))，这2 个模型的因变量与自变量均呈正相关，模型分别达到极显著水平(P<0.01)和显著水平(P<0.05)。对于仅施用化肥的小麦而言，两者间抛物线模型未达显著水平，抛物线在产量约7 t/hm2时氮淋洗量最高，约为60 kg/hm2(图4(b))。

图3 黄淮海平原氮淋洗新指标Fig.3 Yield-scaled N leaching loss in North China Plain

(a)、(b)分别为施用化肥条件玉米氮淋洗量、小麦氮淋洗量。(c)为秸秆还田条件下小麦玉米周年轮作氮淋洗量。(a) and (b) represent maize N leaching and wheat N leaching by applying chemical fertilizer, respectively. (c) N leaching in the case of WMRYRSR.图4 黄淮海平原氮淋洗量与作物产量间的曲线图Fig.4 Relationship curves between N leaching loss and crop yield in North China Plain

在黄淮海平原，仅施用化肥条件下小麦、玉米氮淋洗新指标模型均达到极显著水平。小麦氮淋洗新指标与施氮量呈正相关，与土壤粘粒含量呈负相关。玉米则与施氮量、全N含量呈正相关(表4)。玉米新指标模型中权重值较高的自变量为全N含量，权重值为0.24。小麦新指标模型采用后向回归方法建模，施氮量、粘粒含量与新指标的相关性均未达显著水平，但二元回归模型达到极显著水平，施氮量表现出正向作用，粘粒含量为负向作用。虽然灌溉与降雨量因子未入选模型，但其与玉米氮淋洗新指标呈显著正相关，表明了小麦季与玉米季的不同降雨分布特点。据文献数据显示，小麦季的灌溉降雨量之和变幅较小，均值和标准差分别为407、56 mm；而在玉米季，降雨期较集中且降雨量变幅较大，水分供应以降雨为主，均值和标准差分别为370、138 mm。在秸秆还田条件下，由于文献数据较少，故本研究利用一个轮作周期内小麦和玉米两季数据进行计算，氮淋洗新指标模型仅含施氮量因子，且达极显著相关(表4)。

表4 黄淮海平原仅施用化肥和秸秆还田条件下氮淋洗新指标模型系数[29-37]Table 4 Yield-scaled N leaching model by only applying chemical fertilizer and the model of wheat-maize round year rotation by implementing straw returning tillage in North China Plain[29-37]

2.3 东北平原新指标影响因子

在东北平原，施用化肥条件下玉米氮淋洗新指标值比水稻高，玉米平均值为2.98 kg/t，水稻值为1.59 kg/t[38-40](图5)。从玉米的氮淋洗量与产量间的曲线关系看，两者间存在开口向下的抛物线变化，模型达到显著水平(P<0.05)，而水稻的则无明显的变化趋势，模型未达到显著水平(图6)。

图5 东北平原氮淋洗新指标Fig.5 Yield-scaled N leaching loss in Northeast Plain

(a)、(b)分别为施用化肥条件水稻氮淋洗量、玉米氮淋洗量。(a) and (b) represent rice N leaching and maize N leaching by applying chemical fertilizer, respectively. 图6 东北平原氮淋洗量与作物产量间的曲线图Fig.6 Relationship curves between N leaching loss and crop yield in Northeast Plain

由于东北平原的文献相对较少，未在文中专门列出玉米、水稻的单位产量氮淋洗量模型。模型结果显示，玉米单位产量氮淋洗量与降雨量呈负相关，与化肥氮施用量呈正相关，降雨量的权重值高于施氮量。而水稻单位产量氮淋洗量与灌溉降雨之和、化肥氮施用量均呈正相关，灌溉降雨之和是主要的影响因子。

3 讨 论

3.1 不同区域对新指标的影响

我国制订化肥、农药零增长政策的出发点是为了更科学有效地施用化肥、农药，减少不合理的施用量，通过开展有针对性的管理技术措施，进而减少面源污染流失量和提高农业生产效率[41-42]。这与本研究中新指标的背景相契合，即在提高效率的基础上不断提高农产品产量，将效率优势转换为经济优势。

由于不同区域的气象、土壤等自然环境条件不同，必然会导致区域间新指标影响因素的差异。江苏环太湖地区是长江中下游平原的重要农业产区，该区为北亚热带季风性湿润气候区，气候温和湿润，热量充足，降雨丰沛，雨热同季，年平均气温15.4～16 ℃，>10 ℃的积温为4 800～5 000 ℃，年总降水量1 010～1 200 mm。主要土壤类型为水稻土，因长期稻、麦、棉水旱轮作，剖面中沙黏间层明显，故多漏水漏肥，形成漏水水稻田[43]。在长江中下游平原，影响新指标的主要因素不是施氮量，而是土壤全N含量、降雨、施磷量等因素，应加强土壤养分管理和田间水利基础设施建设，稳步减少施氮量。同样，Zhang 等[44]在江苏常熟开展的2 个周期稻麦轮作试验得出类似的结果，在减少施用氮肥的条件下对水稻的产量影响并不明显，水稻季的后作小麦的产量则与其前作的土壤氮盈余呈显著相关。

黄淮海平原属于暖温带半湿润、湿润气候区，年平均气温8～16 ℃，>0 ℃的积温在4 200～5 500 ℃，无霜期较长。光照充足，年日照为 2 800～2 100 h。年降水量为700～900 mm。潮土是该区域的主要土壤类型，有机质含量普遍偏低。潮土土体深厚，结构较松，宜于耕作管理。该区域影响新指标的主要因素是施氮量，同时也存在土壤质地和养分含量等因素，表明该地区存在满足减少氮流失和维持现有产量条件的最优施氮量。应重点围绕优化施氮量来开展田间管理，并考虑土壤养分管理和质地改良措施，减少氮流失量。这与一项基于黄淮海平原农田的多年的淋洗试验的结论相似，当施氮量介于180～240 kg/hm2可兼顾作物产量和淋溶水水质[45]。

东北平原区属温带季风气候类型，年降水量300～900 mm，年均气温0.7～9.0 ℃，无霜期为90～185 d，热量适于一年一熟。黑土是该区域的主要土壤类型，黑土区的地形主要是山麓平原台地及山前洪积阶地，地势平缓辽阔，呈波状漫岗起伏，坡度一般在1°～5°。影响该区域新指标的影响因子主要是水分，另外还有施氮量。研究结果表明水肥管理对于东北黑土区玉米、水稻氮利用率的影响[46-47]。玉米产量与氮肥利用率呈正相关，随着降雨量的增加促进了玉米的生长，提高了氮肥利用率，从而降低了氮的淋洗风险。对于黑土区水稻籽粒氮含量与氮肥的施用方法和施用量具有统计相关性，且氮肥利用率还受到水氮运筹方法的影响，不当的水肥管理方式会造成氮损失的增加。

本研究结果显示，施化肥条件下黄淮海平原的玉米、小麦氮淋洗新指标达到4.66～5.87 kg/t，而长江中下游平原的水稻、小麦值为0.89～2.26 kg/t，这与两地的环境有关：黄淮海平原的高施氮量造成较高的土壤氮残留，极易随灌溉和降雨而淋失[30]；而长江中下游平原发生的氮反硝化、氨挥发和径流损失较高，而淋洗量相对较少，施肥不会带来过多的氮残留[25]。东北平原水稻、玉米的氮淋洗新指标值为1.59～2.98 kg/t，介于长江中下游平原与黄淮海平原之间，这与当地一年一熟的种植制度及氮肥的激发效应存在关系[47]。相对于一年两熟制，一年一熟制，肥料氮投入减少，因土壤中的残留氮升高引起的氮淋洗风险也有所降低，但是东北平原地区增加氮肥时存在着正激发效应，又促进了土壤氮残留的增加。

3.2 不同管理措施间新指标影响因子差异

1)施用化肥

灌溉和降雨是引起农田氮流失的2个主要的水分因子。长江中下游平原的小麦氮径流新指标模型的影响因子包括了降雨量因子，但土壤全N含量的权重更高。而水稻氮径流新指标模型未包括降雨因子，这与水分管理措施间差异有关。小麦仅靠降雨补给水分，而水稻采用了泡田、晒田等灌溉措施，减弱了降雨对水稻氮径流的影响。已有测定水稻氮径流量的方法既有仅考虑降雨引的损失量(NI=16.3 kg/t，n=16)[19,20,22,25]的类型，也有同时考虑排水和降雨二者之和的损失量(NI=26.3 kg/t，n=14)[21]的类型。东北平原模型数据主要由单一地点试验数据构建，涉及到玉米和水稻2个氮淋洗新指标模型的主要因子均为水分因子。由于东北玉米主要是雨养种植方式，在降雨量未能造成显著氮淋洗的条件下，随着降雨量的升高玉米产量的增加幅度超过了氮淋洗量，就会导致新指标与降雨量呈负相关的现象。这也得到了玉米氮淋洗量与产量间关系图的进一步印证，当产量增加到一定程度后玉米氮淋洗量开始下降，表明了降雨增加导致了玉米强烈的氮吸收作用，减少了氮的淋洗(图6)。

施氮量及土壤养分对新指标具有重要影响。黄淮海平原和东北平原2个北方地区新指标模型均包括施氮量因子，但南方长江中下游平原新指标模型未包括施氮量。长江中下游平原氮肥反硝化损失严重，产量因基础肥力较高受施氮量的影响较小，施氮量对氮流失和产量的影响有限[21]。长江中下游平原的小麦氮径流、水稻氮淋洗和黄淮海平原玉米氮淋洗新指标模型均包括了全N含量因子，这表明除肥料外土壤氮残留或矿化作用对新指标具有重要影响，通过适当减少施氮量水平仍能保证较高的环境经济效益。黄淮海平原小麦氮淋洗新指标模型包括了土壤粘粒含量、施氮量2 个因子，表明存在施肥与土壤质地间的互作作用，应根据土壤条件合理施肥，减少氮淋洗损失。

此外，长江中下游平原水稻氮径流新指标与施磷量呈显著正相关。施磷量与氮径流量和籽粒产量分别呈正相关和负相关，却均未达到显著水平(数据未列出)，这表明增加施磷量可能通过促进植物根系生长和土壤氮的矿化增加了氮径流量，施磷量过高也会影响土壤铁和锌等微量元素的供应进而降低作物产量[48]，从而使施磷量成为显著影响因子。

2)秸秆还田

开展秸秆还田措施对于氮淋洗的影响研究结果不尽相同：有研究认为秸秆覆盖减少地面蒸发，在增加水分下渗的同时加剧了氮淋洗[49]；也有研究发现秸秆覆盖条件下存在水肥互作效应，能促进小麦对氮肥的吸收利用，减少了土壤中氮残留[50-51]。本研究数据显示，施氮量是小麦玉米周年轮作体系新指标的主要影响因子，相对于仅施用化肥秸秆还田减轻了氮淋洗。但是随着小麦玉米的周年产量达到较高水平后，氮淋洗量急剧增加(图4(c))。这表明在较高的水肥条件下，伴随着秸秆还田措施的疏松土壤和增加秸秆氮的矿化的双重作用，促进了氮淋洗的发生。在华北平原黑龙港流域开展的秸秆还田试验结果就表明，降雨量的增加会使0～90 cm土壤的硝态氮残留量显著增加[52]。

3)缓控释肥

施用缓控肥条件下，长江中下游平原的水稻新指标模型仅含土壤全N含量因子，而施氮量的影响不显著，表明了缓控释肥符合作物的吸肥规律，可减轻氮淋洗发生。与施用化肥条件下的模型相比，缓控释肥的土壤全N的权重更高，二者分别为0.49、0.80。这表明缓控释肥条件下由于化肥氮逐步释放，水稻生长需氮更依赖于土壤肥力，从而加速了土壤矿化进程，使得土壤全N含量与新指标相关性更为显著。与在湖北江汉平原开展的水稻-油菜轮作的缓释肥试验结果相似，缓释肥能够稳定水稻产量，并减少氮的淋洗量[53]。在太湖直湖港流域开展水稻氮淋洗试验结果发现，缓控释肥相对于常规施肥增加了氮肥恢复效率，显著增加了氮肥农学效益，进而减少了氮淋洗损失量[54]。

本研究采用了3 个平原部分地区的文献数据进行分析，文献量相对较少，不足以全面反映区域总体情况，分析结果仅限于典型区域的分析。另外，由于区域内水肥管理措施间的差异，也会对结果造成一定影响。今后有待于进一步增加研究样本量，减少管理措施差异影响，提高研究结果的针对性。

4 结 论

本研究围绕单位产量氮流失量这一新指标开展了文献收集和整合，比较了不同区域、施肥管理措施间的新指标差异，重点围绕长江中下游平原和黄淮海平原开展了新指标影响因子分析，并提出了区域尺度下的水肥管理优化措施：1)在长江中下游平原，由于其影响因子不包括施氮量，需通过生态沟渠、缓冲带来拦截、回收利用田间径流氮应扩大缓控释肥施用面积；2)在黄淮海平原，影响新指标的影响因子中均包括施氮量，因此应该优化减少施氮量，控制砂质土壤区的水肥施用量，利用在适宜施氮量条件下秸秆还田措施兼具减少新指标值和维持周年产量的优势；3)在东北平原，应加强水氮运筹以提高氮利用率和减少氮淋洗损失。

致谢

感谢生态环境部环境规划院郑利杰、中国环境科学研究院熊向艳和青岛农业大学资源与环境学院的周奕廷、卢家森在数据处理中的细致工作。