城鎮(zhèn)化對湟水河上游水質(zhì)和底棲動物群落結(jié)構(gòu)的影響

李 寧,陳阿蘭,楊春江,孫瑜旸,馬國良,馬 琪

青海大學,西寧 810016

?

城鎮(zhèn)化對湟水河上游水質(zhì)和底棲動物群落結(jié)構(gòu)的影響

李 寧*,陳阿蘭,楊春江,孫瑜旸,馬國良,馬 琪

青海大學,西寧 810016

湟水河；城鎮(zhèn)化；大型底棲無脊椎動物；不透水地表面積比(PIA)

青藏高原是亞洲多條大江大河的發(fā)源地,境內(nèi)河流縱橫,湖泊星羅棋布,被稱為“中華水塔”。同時青藏高原植被覆蓋度低,土層薄,水土流失嚴重,生態(tài)環(huán)境脆弱、敏感,在城鎮(zhèn)化的過程中尤其要關(guān)注對于高原水生態(tài)環(huán)境的監(jiān)測和保護。

青藏高原城鎮(zhèn)化自20世紀90年代開始發(fā)展迅速,牧區(qū)興建居民點,農(nóng)區(qū)撤村并鄉(xiāng)、撤鄉(xiāng)并鎮(zhèn),城市規(guī)模不斷擴大[1- 2]。然而城市化作為衡量一個國家發(fā)展水平的重要標志,在一定程度上增大了人類社會與生態(tài)環(huán)境之間的相互作用,從而引發(fā)了一系列的社會-環(huán)境-生態(tài)問題。城鎮(zhèn)化過程中造成包括建筑物、道路、廣場、居民區(qū)等不透水地表面積比(PIA)迅速增加,農(nóng)田、林地和草地面積迅速下降,自然河道的渠道化改造。土壤接納降水量急劇下降,地表徑流系數(shù)和徑流量增加,極端降水時間增多。平時基流量大幅減少,暴雨后河水暴漲,沖刷堤岸,城市段河流淤積物增加[3]。降雨后富含氮磷等成分的雨水和城市生活污水通過地表徑流和城市管網(wǎng)排入溪流。加之城市附近的工礦企業(yè)排放的工業(yè)污水,造成河流氮磷等營養(yǎng)鹽及有毒化學品進入河流。隨著城市溪流水文變化的加劇以及溪流生境和水質(zhì)的退化,溪流水生生物豐富度急劇下降,敏感物種大量消失,耐污種類密度急劇增加[4]。

衡量流域內(nèi)城鎮(zhèn)化強度大小的指標可用流域內(nèi)不透水地表面積比例(PIA)表示,即流域內(nèi)城市住宅、工業(yè)、商業(yè)、道路等占用土地面積占流域總面積的比例,用以研究城鎮(zhèn)化對溪流水生生物的影響。研究表明,PIA與底棲動物群落之間存在著線性或非線性的關(guān)系,隨著PIA值的升高,底棲動物多樣性開始持續(xù)下降,當PIA增大到一定閾值后,底棲動物多樣性會保持基本不變或持續(xù)下降[5- 8]。

本文以青藏高原黃河重要支流湟水河為例,探討高原地區(qū)城鎮(zhèn)化對溪流水質(zhì)和生物的影響,并以PIA作為城鎮(zhèn)化強度的量化指標,分析其與底棲動物群落間的關(guān)系,為我國西北高原地區(qū)城鎮(zhèn)化建設(shè)過程中的規(guī)劃管理、溪流生態(tài)保護與修復(fù)提供科學數(shù)據(jù)。

1 研究地區(qū)、內(nèi)容與方法

1.1 研究地區(qū)概況

湟水河為黃河上游最大的支流,源頭為青藏高原祁連山脈的包呼圖山,海拔4400 m,于甘肅永靖縣匯入黃河,全長349 km。其上游為源頭至西寧市小峽口段,河流長184 km,小峽口海拔2178 m,湟水河上游跨越青藏高原和黃土高原(36.55°—37.29°N,100.80°—101.92°E),流域總面積16100 km2,流域內(nèi)地勢西北高、東南低,高山深谷。氣候?qū)儆诘湫痛箨懶詺夂?垂直變化明顯,且地域差異大。愈向上游氣溫愈低,降水量增大,蒸發(fā)量減小,多潮濕沼澤地。流域年平均氣溫0.6℃—7.9℃,年均降水量500 mm,局部地區(qū)可達600 mm。湟水干流谷地6—9月降水占全年降水量的70%左右,且多暴雨。無霜期西北部山區(qū)為31 d,東南部丘陵區(qū)為130—180 d。

圖1 湟水河上游27個采樣點分布 Fig.1 locations of 27 sampling sites in the upstream of Huangshui River

湟水河為羽狀水系,發(fā)源于祁連山脈及其支脈,高山大川,源頭人跡罕至,地理隔絕明顯。為保證多樣性數(shù)據(jù)采樣的完整性,設(shè)湟水河流域源頭參照樣點7個,其中干流樣點1個,一級支流6個。湟水河源頭至西寧市段城鎮(zhèn)化程度逐漸增強,從牧區(qū)村鎮(zhèn)、牧業(yè)區(qū)縣城、農(nóng)區(qū)村鎮(zhèn)、農(nóng)區(qū)縣城、西寧市郊至城鎮(zhèn)化程度最高的西寧市,干擾樣點設(shè)置原則為：在每一城鎮(zhèn)及其上游3—5 km和下游1 km各設(shè)1個采樣點,各支流匯入點設(shè)1個采樣點,共20個干擾樣點(圖1)。

1.2 流域內(nèi)土地利用類型計算

從青海省國土資源廳獲得研究區(qū)域2010年第二次土地利用狀況調(diào)查土地利用類型數(shù)據(jù),結(jié)合30 m分辨率的TM影像監(jiān)督分類結(jié)果,獲得研究區(qū)域內(nèi)不同的土地利用類型。然后將研究區(qū)域DEM格式轉(zhuǎn)換為SWAT格式,提取河網(wǎng)與不同低級流域的界限,再將河網(wǎng)、流域及不同土地利用類型的圖層進行疊加,利用ArcGIS 10.2劃分每個采樣點上游流域的邊界,并利用1∶50000的電子地圖對這些邊界進行校正。將流域內(nèi)土地利用類型按照草地、林地、耕地和城市用地4個類群進行分類匯總,并進行量化,最后計算出草地、林地、耕地和城市用地面積。

1.3 樣品采集與測定

1.3.1 水體理化指標測定

1.3.2 大型底棲無脊椎動物

1.4 數(shù)據(jù)處理

多樣性指數(shù)和豐富度指數(shù)計算參照Shannon[16]及Ganfin[17]的方法,用Duncan檢驗對數(shù)據(jù)進行參照和干擾樣點間環(huán)境和生物參數(shù)的方差分析和差異顯著性檢驗(α=0.05)。環(huán)境變量和生物指標與土地利用類型間的Spearman 秩相關(guān)分析在SPSS 22.0軟件中進行。不透水地表面積比(PIA)與生物指標進行曲線擬合在Origin 8.0軟件中完成。

2 結(jié)果分析

2.1 理化指標與土地利用的關(guān)系

2.2 底棲動物群落與土地利用的關(guān)系

2.2.1 底棲動物群落組成及數(shù)量

調(diào)查共采集到大型底棲動物31科58屬60個分類單元,共5738頭。其中昆蟲綱7目46個分類單元,占所有底棲動物種類的76.67%,昆蟲個體數(shù)共4301頭,占所有底棲動物數(shù)量的74.96%。軟體動物門2個分類單元(3.3%)；環(huán)節(jié)動物門、甲殼綱各1個分類單元(1.67%)。

2.2.2 生物指標與土地利用關(guān)系

總分類單元數(shù)、Shannon多樣性指數(shù)、豐富度指數(shù)、EPT(%)及捕食者(%)、撕食者(%)、濾食者(%)、刮食者(%)、敏感類群(%)與PIA呈顯著負相關(guān),與林地(%)和草地(%)呈顯著正相關(guān)。但BI指數(shù)、集食者(%)、耐污類群(%)、寡毛綱(%)與PIA呈顯著正相關(guān),與草地(%)和林地(%)呈顯著負相關(guān)(表2)。

表1 采樣點環(huán)境變量值及其與土地利用的關(guān)系

*P<0.05; **P<0.01

表2 采樣點底棲動物群落參數(shù)與土地利用的關(guān)系

EPT: 水生昆蟲中蜉蝣目Ephemeroptera、襀翅目Plechopetra和毛翅目Trichopetra;*P<0.05; **P<0.01

總分類單元數(shù)、豐富度指數(shù)、Shannon多樣性指數(shù)和EPT(%)與PIA進行線性或指數(shù)擬合結(jié)果見圖2。隨著PIA值的升高,4項指標呈現(xiàn)明顯的下降趨勢,Shannon多樣性指數(shù)、豐富度指數(shù)與PIA呈顯著的線性關(guān)系,兩項值持續(xù)下降并趨于零；總分類單元數(shù)、EPT(%)則呈非線性關(guān)系。前者下降后穩(wěn)定在一定范圍內(nèi),而EPT(%)在PIA達到約40%是則完全消失。

圖2 總分類單元數(shù)、Shannon多樣性指數(shù)、豐富度及EPT與PIA的關(guān)系Fig.2 Relationships between total number of taxa, Shannon diversity index, richness index, EPT individuals and PIA

除集食者外,底棲動物取食功能類群與PIA呈現(xiàn)顯著的非線性關(guān)系,隨著PIA值升高,集食者數(shù)量持續(xù)增加,捕食者、撕食者、刮食者和濾食者的數(shù)量卻隨PIA值升高快速下降直至全部消失(圖3)。

敏感類群(%)、寡毛綱(%)與PIA同樣呈現(xiàn)顯著的非線性關(guān)系,而BI指數(shù)和耐污類群(%)與PIA 呈現(xiàn)顯著的線性關(guān)系。隨著PIA值的升高,BI指數(shù)、寡毛綱(%)和耐污類群(%)持續(xù)升高。敏感類群(%)在PIA值升至6%時即降低到最低水平并保持穩(wěn)定(圖4)。

3 討論

3.1 理化指標與PIA的關(guān)系

3.2 生物指標與PIA的關(guān)系

本研究表明,首先,總分類單元數(shù)、EPT(%)和捕食者(%)、撕食者(%)、濾食者(%)、刮食者(%)及敏感

圖3 底棲動物取食功能團與PIA的關(guān)系Fig.3 Relationships between functional groups of macrobenthos and PIA

圖 4 敏感類群、BI指數(shù)、寡毛綱和耐污類群與PIA關(guān)系Fig.4 Relationships between sensitive taxa, BI index, oligochaeta, tolerance taxa and PIA

類群(%)與PIA呈顯著負相關(guān),與林地(%)、草地(%)呈顯著正相關(guān)。而BI指數(shù)、集食者(%)、耐污類群(%)和寡毛綱(%)與PIA呈顯著正相關(guān),與草地(%)、林地(%)呈顯著負相關(guān)。這與Walsh等[3]提出的觀點一致。即隨著城鎮(zhèn)化強度的增加,水體中的底棲動物敏感類群減少,耐污類群增加,寡毛綱逐漸成為優(yōu)勢類群。其次,從取食功能類群方面來看,隨著PIA值的升高,溪流中的藻類等周叢生物數(shù)量減少,食物缺乏導致刮食者數(shù)量減少；由于流速降低及水文條件的改變,撕食者、濾食者和捕食者的數(shù)量也不斷減少,但水中有機碎屑沉積和增加導致了集食者增多[20-21]。第三,從生物指標與PIA之間曲線擬合來看,劉東曉[18]的研究顯示均為非線性關(guān)系,而本研究顯示底棲動物生物指標中Shannon多樣性指數(shù)、豐富度指數(shù)、集食者(%)、耐污類群(%)和BI指數(shù)與PIA呈現(xiàn)線性相關(guān),與King[7]等和Walsh[3]等的結(jié)論基本一致。分析原因,青藏高原溪流的特點是海拔高、水溫低、污染少,水體中營養(yǎng)物質(zhì)少,地理隔絕現(xiàn)象明顯,底棲生物種類較少,但數(shù)量相對較高,對于干擾反應(yīng)表現(xiàn)出線性的關(guān)系。而對于非線性的關(guān)系而言,表明PIA對生物指標的影響存在著一個閾值。當PIA值升高到6%時,底棲生物中的敏感類群(%)迅速降低到很低的水平直至消失(圖4)。這與Beach[22]和Stepenuck[23]等提出的PIA大于10%或在8%—12%時,底棲動物群落就有可能發(fā)生突變的結(jié)論基本一致,但本研究顯示其閾值更低,上述結(jié)論都還有待于今后更多的數(shù)據(jù)予以驗證。

4 結(jié)論

[1] 章輝. 青藏高原牧區(qū)城鎮(zhèn)化實證研究[D]. 蘭州: 西北民族大學, 2006.

[2] 馬維勝. 青藏高原生態(tài)城市化模式研究. 青海民族研究: 社會科學版, 2002, 13(4): 6- 10.

[3] 張建云, 宋曉猛, 王國慶, 賀瑞敏, 王小軍. 變化環(huán)境下城市水文學的發(fā)展與挑戰(zhàn)——Ⅰ.城市水文效應(yīng). 水科學進展, 2014, 25(4): 594- 605.

[4] Neal C, Robson A J. A summary of river water quality data collected within the land-ocean interaction study: core data for eastern UK rivers draining to the North Sea. Science of the Total Environment, 2000, 251- 252: 585- 665.

[5] Thomson J D, Weiblen G, Thomson B A, Alfaro S, Legendre P. Untangling multiple factors in spatial distributions: lilies, gophers, and rocks. Ecology, 1996, 77(6): 1698- 1715.

[6] Booth D B, Karr J R, Schauman S, Konrad C P, Morley S A, Larson M G, Burges S J. Reviving urban streams: land use, hydrology, biology, and human behavior. Journal of the American Water Resources Association, 2004, 40(5): 1351- 1364.

[7] King R S, Baker M E, Whigham D F, Weller D E, Jordan T E, Kazyak P F, Hurd M K. Spatial considerations for linking watershed land cover to ecological indicators in streams. Ecological Applications, 2005, 15(1): 137- 153

[8] Walsh C J, Fletcher T D, Ladson A R. Stream restoration in urban catchments through redesigning stormwater systems: looking to the catchment to save the stream. Journal of the North American Benthological Society, 2005, 24(3): 690- 705

[9] 金相燦, 屠清瑛. 湖泊富營養(yǎng)化調(diào)查規(guī)范(第二版). 北京: 中國環(huán)境科學出版社, 1990: 159- 188.

[10] 國家環(huán)保局. 水和廢水監(jiān)測分析方法(第三版). 北京: 中國環(huán)境科學出版社, 1989: 252- 286.

[11] Barbour M T, Gerritsen J, Griffith G E, Frydenborg R, McCarron E, White J S, Bastian M L. A framework for biological criteria for Florida streams using benthic macroinvertebrates. Journal of the North American Benthological Society, 1996, 15(2): 185- 211.

[12] Hauer F R, Lamberti G A. Methods in Stream Ecology. 2nd ed. San Diego: Academic Press, 2007: 51- 78.

[13] Merritt R W, Cummins K W, Berg M B. An Introduction to the Aquatic Insects of North America. Dubuque: Kendall/Hunt Publishing Company, 1996.

[14] 王備新. 大型底棲無脊椎動物水質(zhì)生物評價研究[D]. 南京: 南京農(nóng)業(yè)大學, 2003: 63- 72.

[15] Qin C Y, Zhou J, Cao Y, Zhang Y, Hughes R M, Wang B X. Quantitative tolerance values for common stream benthic macroinvertebrates in the Yangtze River Delta, Eastern China. Environmental Monitoring and Assessment, 2014, 186(9): 5883- 5895.

[16] Shannon C E, Weaver W. The Mathematical Theory of Communication. Urbana: The University of Illinois Press, 1949.

[17] Gaufin A R. The use and value of aquatic insects as indicators of organic enrichment // Tarzwell C M, ed. Biological Problems in Water Pollution. Cincinnati: US Public Health Service, 1957: 136- 143.

[18] 劉東曉, 于海燕, 劉朔孺, 胡尊英, 俞健, 王備新. 城鎮(zhèn)化對錢塘江中游支流水質(zhì)和底棲動物群落結(jié)構(gòu)的影響. 應(yīng)用生態(tài)學報, 2012, 23(5): 1370- 1376.

[19] 謝有仁. 湟水流域水系組成及分布特征. 水利科技與經(jīng)濟, 2011, 17(1): 72- 73.

[20] Miller W, Boulton A J. Managing and rehabilitating ecosystem processes in regional urban streams in Australia. Hydrobiologia, 2005, 552(1): 121- 133.

[21] Paul MJ. Stream ecosystem function along a land-use gradient[D]. Georgia: University of Georgia, 1999.

[22] Beach D. Coastal Sprawl: The Effects of Urban Design on Aquatic Ecosystems in the United States. Arlington VA: Pew Oceans Commission, 2002.

[23] Stepenuck K F, Crunkilton R L, Wang L Z. Impacts of urban landuse on macroinvertebrate communities in southeastern wisconsin streams. Journal of the American Water Resources Association, 2002, 38(4): 1041- 1051.

Impacts of urbanization on water quality and macrobenthos community structure upstream in the Huangshui River

LI Ning*, CHEN Alan, YANG Chunjiang, SUN Yuyang, MA Guoliang, MA Qi

QinghaiUniversity,Xining810016,China

Huangshui River; urbanization; macro-benthic invertebrate; percentage of impervious area (PIA)

國家自然科學地方基金資助項目(31560601)；青海省自然基金資助項目(2014-ZJ- 905)

2016- 03- 12; 網(wǎng)絡(luò)出版日期:2017- 02- 17

10.5846/stxb201603120438

*通訊作者Corresponding author.E-mail: newlining@sina.com

李寧,陳阿蘭,楊春江,孫瑜旸,馬國良,馬琪.城鎮(zhèn)化對湟水河上游水質(zhì)和底棲動物群落結(jié)構(gòu)的影響.生態(tài)學報,2017,37(10):3570- 3576.

Li N, Chen A L, Yang C J, Sun Y Y, Ma G L, Ma Q.Impacts of urbanization on water quality and macrobenthos community structure upstream in the Huangshui River.Acta Ecologica Sinica,2017,37(10):3570- 3576.