林业科学  2009, Vol. 45 Issue (8): 76-81   PDF    
0

文章信息

吴家森, 姜培坤, 盛卫星, 谢秉楼.
Wu Jiasen, Jiang Peikun, Sheng Weixing, Xie Binglou
雷竹集约栽培对周边河流水质的影响
Effect of Bamboo (Phyllostachys praecox f. prevernalis) Intensive Cultivation on Water Quality in Water Systems
林业科学, 2009, 45(8): 76-81.
Scientia Silvae Sinicae, 2009, 45(8): 76-81.

文章历史

收稿日期:2007-12-13

作者相关文章

吴家森
姜培坤
盛卫星
谢秉楼

雷竹集约栽培对周边河流水质的影响
吴家森, 姜培坤, 盛卫星, 谢秉楼     
浙江林学院环境科技学院 临安 311300
摘要: 为探讨雷竹集约化栽培对周边河流水质的动态变化的影响,在雷竹栽培区与天然林区河流各设置3个水质监测断面,于2006年5月—2007年5月每2个月进行采样,对主要水化指标进行分析。结果表明:雷竹栽培区河流水质状况明显差于天然林区,水体中所有水化指标均表现为雷竹栽培区高于天然林区。总氮、硝氮、氨氮和有机氮含量丰水期高于枯水期,而总磷、可溶性碳、COD、BOD5含量则表现为丰水期小于枯水期。雷竹栽培区水体中总氮含量年均达4.24 mg·L-1,是天然林区水体的3.15倍,硝氮含量为2.01 mg·L-1,是天然林区的5.07倍,氨氮含量为0.07 mg·L-1,比栽培区水体高0.66倍;水体中总磷含量雷竹栽培区平均达0.05 mg·L-1,是天然林区的2.39倍;可溶性碳含量平均达17.78 mg·L-1,比天然林区高出2.87倍。雷竹集约化栽培引起了周边水体氮、磷和可溶性碳的污染。
关键词:水质    动态变化    雷竹    
Effect of Bamboo (Phyllostachys praecox f. prevernalis) Intensive Cultivation on Water Quality in Water Systems
Wu Jiasen, Jiang Peikun , Sheng Weixing, Xie Binglou    
School of Environment Science and Technology, Zhejiang Forestry College Lin'an 311300
Abstract: Water samples from two near rivers respectively around bamboo stands and natural forests were collected in three sampling points in two months interval from May 2006 to May 2007. It was found that the quality of water systems under the bamboo stands was much worse than that under the natural forest. The concentration of total N, nitrate N, and organic N in two rivers appeared higher in wet season (during March and May), than the other months that had less rainfall. However, the concentration of total P and dissolved C and the values of COD and BOD were much higher in dry season (during September and January), and the peak values of water from river around bamboo stands appeared in November. The average of concentrations of total N in water from river around bamboo stands was 4.24 mg·L-1 and 2.15 times greater than that around the natural forest. The corresponding values for NO3-N was 2.01 mg·L-1 and 4.07 times, for NH4+-N was 0.07 mg·L-1 and 0.66 times, for total P was 0.05 mg·L-1 and 1.39 times, for dissolved organic C was 17.78 mg·L-1 and 2.87 times. The results showed that the river around bamboo stands was polluted by N, P and DOC.
Key words: water quality    dynamics    Phyllostachys praecox f. prevernalis    

雷竹(Phyllostachys praecox f. prevernalis)是优良的笋用竹种。近20年来,通过集约化经营,采用冬季覆盖和重施肥栽培技术,使雷竹笋实现了反季节生产,并使产量上升,从而使竹农获得可观的经济效益(方伟等,1994)。随着效益林业的实施,雷竹在浙江的栽培面积不断扩大,目前已成为许多地区农民的经济支柱产业。雷竹栽培中,肥料施用量很大,每年施入化肥超过3 t·hm-2,有机肥(厩肥)达100 t·hm-2(汪祖潭等,1995)。氮肥施用量普遍超过水资源免受污染国际限量标准的4倍(崔玉亭,1999),并且超量施肥的势头有增无减,有的农民几乎已达到滥用肥料的地步。超量施用化学肥料已使大部分雷竹林土壤化学养分严重积累(杨芳等,2003),磷素积累尤其突出(姜培坤等,2000a),化学肥料超量施用也必然会造成氮、磷、钾等化学元素在环境中富集,最主要是污染周边水体。前人对集约经营雷竹林栽培技术、林地土壤特征及对笋体质量的影响均作了较多的研究(姜培坤等,19992000b2003a2004a2004b2005a2006),但该技术的应用对栽培区周边地表水的性质有何影响,至今尚未见报道。

本文通过采用标准水质监测方法,探讨了雷竹栽培区与天然林区流域地表水水质的动态变化及差异,为雷竹栽培区水体污染特征、水污染控制提供理论与试验依据。

1 材料与方法 1.1 研究区概况

研究区位于浙江省临安市太湖源镇,属中纬度北亚热带季风气候,四季分明,气候温和,雨量充沛,全年降雨量平均1 628 mm。多年平均气温为15.8 ℃,7月为最热月,历年平均为28.1 ℃,1月为最冷月,历年平均为3.4 ℃,极端高温41.9 ℃,极端低温-13.3 ℃。历年平均日照时数1 939 h,无霜期234 d。

雷竹集约经营栽培模式的主要技术是:每年11月下旬至12月上旬给雷竹林地表覆盖有机物以增温避冷保湿,日增地温达4~5 ℃。通常是在雷竹林地表先覆10~15 cm稻草,上面再覆盖10~15 cm的砻糠,稻草用量40 t·hm-2,砻糠用量55 t·hm-2,次年3,4月揭去未腐烂的砻糠;下层的稻草经过一个冬春的发酵和雨雪水的淋泡已基本腐烂入土。每年施3次肥,时间为5月中旬、9月中旬和覆盖前。施肥量为无机复合肥(N:P:K=16:16:16)2.25 t·hm-2和尿素(含N 46%)1.125 t·hm-2,或施等量养分的有机肥。管理较好的林地可在当年的12月下旬出笋,比常规栽培的出笋时间提早2个月,而且产量高质量好,可在春节前上市,大大提高了价格和收益。

1.2 采样与分析方法

根据野外调查及当地政府的建议,在地理、人口等条件基本一致的情况下,选择雷竹栽培区东天目溪和天然林区东坑设置水质监测定位试验。栽培区河流周边的农田、山坡自1995年起均种植了雷竹,并实行集约化栽培,面积达300 hm2,而天然林区河流周边的农田没有进行雷竹栽培,山坡上的植被以天然常绿阔叶幼林为主。在东天目溪设立水质监测段面3个(A1:30°22.056′N,119°35.399′E;A2:30°21.053′N,119°35.296′E;A3:30°21.042′N,119°35.268′E),同时在东坑也设立水质监测段面3个(B1:30°15.242′ N,119°33.164′ E;B2:30°16.232′ N,119°33.512′ E;B3:30°15.582′ N,119°34.157′E)。并于2006年5月至2007年5月,每隔2个月进行1次采样,采样方法按照国家标准进行(GBT 5750.2-2006),监测项目主要包括:总氮、氨氮、硝氮、总磷、可溶性碳、COD、BOD5等指标。

水样采回后分成2份,1份用0.45 μm滤膜抽滤,用于测定总氮(TN)、氨氮、硝氮、可溶性碳(DOC)含量;另1份样品不抽滤,用于测定总磷(TP)、COD和BOD5含量。分析方法如下:TN直接在岛津TOC分析仪中的TNM分析仪测定;氨氮、硝氮直接在戴安公司生产的ICS-1500离子色谱上测定;总磷,硫酸和高氯酸消化后用钼蓝比色法测定(鲁如坤,2002GB/T 5750.5-2006);DOC,直接在岛津TOC分析仪上测定,COD测定采用高锰酸钾法,BOD5测定采用稀释与接种法(GB/T 5750.7-2006)。有机氮通过计算求得,其计算方法:有机氮为总氮减去氨氮和硝氮。所有数据均用DPS软件进行统计分析(唐启义等,1997)。

2 结果与分析 2.1 雷竹栽培区与天然林区水体中氮的年动态变化

图 1可知,雷竹栽培区水体中总氮、硝氮、氨氮、有机氮含量在1年中的各月都明显高于天然林区,但总的变化规律基本一致,总氮、硝氮、有机氮含量均表现为5月为最高,而氨氮则以7月为最高,5月,雷竹栽培区河流总氮含量为8.65 mg·L-1,天然林区河流为2.90 mg·L-1,而后逐渐降低,到第2年3月后又逐渐升高。雷竹栽培区河流中各种氮形态含量在各月份波动较大,而天然林区则相对较平稳。

图 1 雷竹栽培区与天然林区水体氮的年动态变化 Figure 1 Dynamics of N in water system A:总氮Total-N;B:硝氮Nitrate-N;C:氨氮Ammonia-N;D:有机氮Organic-N.
2.2 雷竹栽培区与天然林区水体中磷的年动态变化

图 2结果显示,雷竹栽培区与天然林区河流中总磷含量的年变化规律基本一致,均表现为11月最高,水体中总磷含量雷竹栽培区河流为0.09 mg·L-1,而天然林区河流仅为0.04 mg·L-1,而后水体中磷含量均下降。

图 2 雷竹栽培区与天然林区水体中总磷的年动态变化 Figure 2 Dynamics of P in water system
2.3 雷竹栽培区与天然林区水体中可溶性碳的年动态变化

雷竹栽培区与天然林区水体中可溶性碳的年动态变化如图 3所示。雷竹栽培区水体中可溶性碳含量在不同月份差异较大,以11月的枯水期为最大,达25.31 mg·L-1,而后其含量降低并趋于稳定,至7月达最低,而后又逐渐上升。而天然林区水体中可溶性碳含量的年动态变化规律相对稳定,其含量也远远低于雷竹栽培区水体。

图 3 雷竹栽培区与天然林区水体中DOC的年动态变化 Figure 3 Dynamics of DOC in water system
2.4 雷竹栽培区与天然林区水体中COD和BOD5的年动态变化

化学需氧量(chemical oxygen demand,简称COD),它是表示水中还原性物质多少的一个指标。因此,COD又往往作为衡量水中有机物质含量多少的指标。五日生化需氧量(BOD5),是指在规定条件下水中有机物和无机物在生物氧化作用下所消耗的溶解氧。

已有研究表明,地表水中污染物质量浓度季节变化较大。一般来说,较干旱的季节因地表水流量少,水中污染物质量浓度相对较高;而在雨季,由于稀释作用,地表水中的污染物质量浓度相对较低(Dassenakis et al., 1998; Brown et al., 1998)。雷竹栽培区COD和BOD5的年动态变化也符合以上规律,总体表现为,从2006年5月开始BOD5含量逐渐升高,至11月达最大值,而后缓慢下降。而天然林区水体中COD含量以2006年5月最高,而后逐渐降低,趋于稳定,而BOD5的变化规律则与栽培区水体中的相似。从图 45可知,雷竹栽培区水体中COD和BOD5含量均高于天然林区,但均低于Ⅰ类地表水环境质量标准。说明雷竹集约化栽培尚未造成水体COD、BOD5的污染。

图 4 雷竹栽培区与天然林区水体中COD的年动态变化 Figure 4 Dynamics of COD in water system
图 5 雷竹栽培区与天然林区水体中BOD5的年动态变化 Figure 5 Dynamics of BOD5 in water system
2.5 雷竹栽培区与天然林区水体中各指标年均值比较

将雷竹栽培区与天然林区水体中各水化指标在各月份的含量进行平均,得到年平均值,结果如图 6所示。从图 6可知,雷竹栽培区水体中各水化指标均显著高于天然林区水体。雷竹栽培区水体中总氮(TN)含量平均达4.24 mg·L-1,是天然林区水体的3.15倍;雷竹栽培区水体中硝氮含量为2.01 mg·L-1,是天然林区的5.07倍;雷竹栽培区水体中可溶性碳(DOC))含量达17.78 mg·L-1,比天然林区水体高出2.87倍;雷竹栽培区水体中COD和BOD5含量分别3.57,2.01 mg·L-1,是天然林区的2.34和1.82倍。水体中氨氮和总磷含量较低,难以在图上显示,水体中总磷含量雷竹栽培区平均达0.05 mg·L-1,是天然林区水体的2.39倍;水体中氨氮含量雷竹栽培区平均达0.07 mg·L-1,比天然林区水体高0.66倍。

图 6 雷竹栽培区与天然林区水体中各指标年均值比较 Figure 6 Average of different indexes in water system
3 结论与讨论

雷竹栽培区河流水质状况明显差于天然林区,水体中所有水化指标均表现为雷竹栽培区高于天然林区。氮和磷是引起水体富营养化的2种主要营养盐,它们在水体中的过量富集可使藻类和其他水生植物数量增加,降低水体水质(Smith et al., 1999)。碳是微生物生命活动的能源及生物体结构材料,水体丰富的有机碳可导致水体微生物的大量繁殖,使水体富营养化比氮磷污染更加严重。

3.1 雷竹集约化栽培对水体中氮的影响

据研究(黄芳等,2007)雷竹集约化栽培中纯N的投入(包括施肥、覆盖及氮素沉降等)为1 062 kg·hm-2a-1,通过挥发、反硝化、渗漏淋失和地表径流等损失N素289 kg·hm-2a-1,占氮素投入的27%。氮素通过淋失和地表径流进入地表水和地下水环境系统,使水质受到污染。地表径流和渗漏淋失占氮素投入的11%~16%(段亮等,2007王家玉等,1996高效江等,2001陈国军等,2003),如以14%计,雷竹集约化栽培进入水体的氮素将达到40.46 kg·hm-2a-1。定位试验所在的东天目溪小流域一年中就有12.138 t纯氮进入水体,浙江省雷竹栽培面积达20 000 hm2,则每年进入水体的氮素将达到809.200 t,从而加剧了水体的污染,对生态环境产生危害。本研究表明,雷竹栽培区水体中总氮(TN)含量平均达4.24 mg·L-1,是天然林区水体的3.15倍,水体中总氮、硝氮、氨氮和有机氮等污染物含量则表现为丰水期大于枯水期。水体中氮的年变化规律与该区域的降雨特征、雷竹的经营措施是相吻合。5月,雷竹林地上层覆盖物砻糠已揭去,林地表面留着上一年的半腐烂的稻草,施用复合肥和有机肥,然后翻耕。从4月下旬开始,本地区开始大量降水,到9月下旬雨季基本结束,开始进入干旱季节,降雨较少。由于降雨使氮素随地表径流和渗漏淋失进入河流,从而表现为5月最高。这与张金镇等(2004)对流溪河的研究结果相似。雷竹栽培区施用了大量的氮肥,一遇下雨大量的肥料随水流入河流,导致栽培区水体中各种氮含量明显高于天然林区。与国家地表水环境质量标准(GB3838-2002)相比,雷竹栽培区水体中总氮含量除11月为Ⅳ类外,其余月份均为劣Ⅴ类,这与吴家森(2008)研究结果相似。

3.2 雷竹集约化栽培对水体中磷的影响

磷肥施入土壤后,非常容易与土壤结合,因而其迁移距离和速度均非常低,但雷竹栽培中大量磷肥的使用,一方面在土壤中累积,而另一方面也对周边水体的质量产生了影响。Heckrath等(1995)研究了洛桑实验站土壤有效磷含量为7~110 mg·kg-1田块渗漏水的磷含量,结果表明当地土壤有效磷含量低于57 mg·kg-1时,渗漏水中未检测到全磷和可溶性磷,当土壤有效磷含量大于57 mg·kg-1时,渗漏水中磷质量浓度随土壤有效磷含量的增加而呈线性增加。王彩绒等(2005)将太湖典型地区蔬菜地土壤磷素淋失风险划分为3个等级,最低等级土壤磷几乎不发生淋失,即土壤速效磷上限60 mg·kg-1、中等强度磷淋失指标,土壤速效磷70 mg·kg-1、高强度淋失指标,土壤速效磷大于70 mg·kg-1。对于速效磷含量大于60 mg·kg-1的区域不但要控制土壤径流向水体迁移磷,对磷以淋失、侧渗方式向水体迁移也要引起重视。而雷竹栽培区土壤中有效磷含量平均超过300 mg·kg-1,将会对周边水体产生极大的风险。本研究表明,水体中磷含量在1年的监测中均表现为雷竹栽培区远远大于天然林区,总磷含量在雷竹栽培区河流中年均值达0.05 mg·L-1,是天然林区河流中的2.39倍。对照国家地表水环境质量标准,雷竹栽培区水体中磷含量,11月为Ⅳ类水,其他月份均为Ⅲ类水。

3.3 雷竹集约化栽培对水体中可溶性碳的影响

可溶性有机碳(DOC)代表了水中所含有机物的含量,它反映了水体被有机物质污染的程度,影响着水体的质量(叶麟等,2006)。碳是微生物生命活动的能源及生物体结构材料,水体丰富的有机碳可导致水体微生物的大量繁殖,使水体富营养化比氮磷污染更加严重。雷竹集约经营的重要措施之一是冬季覆盖,每年11月大量的稻草及砻糠覆盖在竹林地表上,在长达4~5个月的冬季地表覆盖过程中,稻草、砻糠等有机物料部分腐烂入土,增加了土壤有机碳总量。有机碳总量的增加一方面使可溶性碳增加,另一方面也会刺激土壤微生物使生物代谢加剧。微生物代谢加剧的结果使土壤可溶性有机物数量上升(Haynes et al., 1993),结果是土壤可溶性有机碳含量会进一步增加。姜培坤等(20022003b2005b)分析了27个雷竹栽培区土壤样品,结果表明,高效雷竹林土壤水溶性有机碳含量平均达155 mg·kg-1,随着覆盖年份增加土壤水溶性有机碳含量明显上升,覆盖1年比未覆盖林地土壤水溶性有机碳含量增加了24.18%,连续覆盖5年的土壤水溶性有机碳含量达231.00 mg·kg-1,比覆盖1年的增加了约1倍。杨芳等(2006)研究结果也表明,雷竹林覆盖处理土壤水溶性有机碳含量是未覆盖处理的1.56倍。土壤水溶性有机碳易溶于水,由于降雨作用,将随着地表径流和下渗水进入周边水体或地下水系统。本研究表明,雷竹栽培区河流中可溶性碳含量明显比天然林区高,月平均含量达17.78 mg·L-1,比天然林区水体高出2.87倍。

参考文献(References)
陈国军, 曹林奎, 陆贻通, 等. 2003. 稻田氮素流失规律测坑研究. 上海交通大学学报:农业科学版, 21(4): 320-324.
崔玉亭. 1999. 化肥与生态环境保护. 北京: 化学工业出版社.
段亮, 段增强, 常江. 2007. 地表管理与施肥方式对太湖流域旱地氮素流失的影响. 农业环境科学学报, 26(3): 813-818. DOI:10.3321/j.issn:1672-2043.2007.03.002
方伟, 何钧潮, 卢学可, 等. 1994. 雷竹早产高效栽培技术. 浙江林学院学报, 11(2): 121-128.
高效江, 胡雪峰, 王少平, 等. 2001. 淹水稻田中氮素损失及其对水环境影响的试验研究. 农业环境保护, 20(4): 196-198, 205.
黄芳, 蔡荣荣, 孙达, 等. 2007. 集约经营雷竹林土壤氮素状况及氮平衡的估算. 植物营养与肥料学报, 13(6): 1193-1196. DOI:10.3321/j.issn:1008-505x.2007.06.034
姜培坤, 徐秋芳, 储家淼, 等. 2006. 雷竹早产高效栽培过程中土壤养分质量分数的变化. 浙江林学院学报, 23(3): 242-247. DOI:10.3969/j.issn.2095-0756.2006.03.002
姜培坤, 徐秋芳, 钱新标, 等. 1999. 雷竹林地覆盖增温过程中土壤化学性质的动态变化. 浙江林学院学报, 16(2): 123-130. DOI:10.3969/j.issn.2095-0756.1999.02.004
姜培坤, 徐秋芳, 罗煦钦, 等. 2004a. 雷竹笋重金属含量及其与施肥的关系. 浙江林学院学报, 21(4): 424-427.
姜培坤, 徐秋芳. 2004b. 雷竹笋硝酸盐含量及其与施肥的关系. 浙江林学院学报, 21(1): 10-14.
姜培坤, 徐秋芳. 2005a. 不同施肥雷竹林土壤重金属含量的动态分析. 水土保持学报, 19(2): 168-170, 180.
姜培坤, 徐秋芳. 2005b. 施肥对雷竹林土壤活性有机碳的影响. 应用生态学报, 16(2): 253-256.
姜培坤, 叶正钱, 徐秋芳. 2003a. 高效栽培雷竹林土壤重金属含量的分析研究. 水土保持学报, 17(4): 61-63, 74.
姜培坤, 徐秋芳, 杨芳. 2003b. 雷竹土壤水溶性有机碳及其与重金属的关系. 浙江林学院学报, 20(1): 8-11.
姜培坤, 俞益武, 金爱武, 等. 2000a. 丰产雷竹林地土壤养分分析. 竹子研究汇刊, 19(4): 50-53.
姜培坤, 俞益武. 2000b. 雷竹叶营养元素含量与土壤养分的关系. 浙江林学院学报, 17(4): 360-363.
姜培坤, 周国模, 徐秋芳. 2002. 雷竹高效栽培措施对土壤碳库的影响. 林业科学, 38(6): 6-11. DOI:10.3321/j.issn:1001-7488.2002.06.002
鲁如坤. 2000. 土壤农业化学分析方法. 北京: 中国农业科技出版社.
唐启义, 冯明光. 1997. 实用统计分析及其计算机处理平台. 北京: 中国农业出版社, 15-70.
汪祖潭, 方伟, 何钧潮, 等. 1995. 雷竹笋用林高产高效栽培技术. 北京: 中国林业出版社.
王彩绒, 胡正义, 杨林章, 等. 2005. 太湖典型地区蔬菜地土壤磷素淋失风险. 环境科学学报, 25(1): 76-80. DOI:10.3321/j.issn:0253-2468.2005.01.013
王家玉, 王胜佳, 陈义, 等. 1996. 稻田土壤中氮素淋失的研究. 土壤学报, 33(1): 28-35. DOI:10.3321/j.issn:0564-3929.1996.01.004
吴家森.2008.雷竹产区水体污染特征及其与施肥的关系.浙江林学院硕士学位论文.
杨芳, 吴家森, 姜培坤, 等. 2006. 不同施肥雷竹林土壤水溶性有机碳动态变化. 浙江林业科技, 26(3): 34-37. DOI:10.3969/j.issn.1001-3776.2006.03.008
杨芳, 徐秋芳. 2003. 不同栽培历史雷竹林土壤养分与重金属含量的变化. 浙江林学院学报, 20(2): 111-114. DOI:10.3969/j.issn.2095-0756.2003.02.001
叶麟, 韩新芹, 蔡庆华. 2006. 香溪河库湾春季水华期间可溶性碳动力学研究. 水生生物学报, 30(1): 80-83. DOI:10.3321/j.issn:1000-3207.2006.01.015
张金镇, 邓熙, 顾继光. 2004. 流溪河水质的动态特征. 生态科学, 23(3): 231-235. DOI:10.3969/j.issn.1008-8873.2004.03.011
GB/T5750.5-2006.2007.生活饮用水标准检验方法: 无机非金属指标.北京: 中国标准出版社. http://www.csres.com/detail/181192.html
GB/T5750.7-2006.2007.生活饮用水标准检验方法: 有机物综合指标.北京: 中国标准出版社. http://www.csres.com/detail/181194.html
GB3838-2002.2002.地表水环境质量标准.北京: 中国标准出版社.
GBT 5750.2-2006.2007.生活饮用水标准检验方法: 水样的采集和保存.北京: 中国标准出版社. http://www.csres.com/detail/181189.html
Brown P L, Markich S J. 1998. Relative importance of natural and anthropogenic influences on the fresh surface water chemistry of the Hawkesbury-Nepenan river, south-eastern Australia. Sci Total Environ, 217(3): 201-230. DOI:10.1016/S0048-9697(98)00188-0
Dassenakis M, Scoullos M, Foufa E, et al. 1998. Effects of multiple source pollution on a small Mediterranean river. Appli Geochem, 13(2): 197-211. DOI:10.1016/S0883-2927(97)00065-6
Haynes R J, Francis G S. 1993. Changes in microbial biomass C, soil carbohyd rate composition and aggregate stability induced by growth of selected crop and forage species under field conditions. J Soil Sci, 44: 665-675. DOI:10.1111/ejs.1993.44.issue-4
Heckrath G, Brookes P C, Poulton P R, et al. 1995. Phosphorus leaching from soils containing different phosphorus concentrations in the broad balk experiment. Journal of Environment Quality, 24: 904-910.
Smith V H, Tilman G D, Nekola J C. 1999. Eutrophication: impacts of excess nutrient inputs on freshwater, marine, and terrestrial ecosystems. Environ Pollute, 100: 179-196. DOI:10.1016/S0269-7491(99)00091-3