廖文婷,鄧紅兵,李若男,鄭 華,*

1 中國生態(tài)環(huán)境研究中心城市與區(qū)域生態(tài)國家重點實驗室,北京 100085 2 中國科學(xué)院大學(xué)資源與環(huán)境學(xué)院,北京 100049

生態(tài)系統(tǒng)的水文調(diào)節(jié)服務(wù)是生態(tài)系統(tǒng)對自然界中水的各種運動變化所發(fā)揮的作用,表現(xiàn)為通過生態(tài)系統(tǒng)對水的利用、濾過等影響和作用以后,水在時間、空間、數(shù)量等方面發(fā)生變化的現(xiàn)象和過程[1]。生態(tài)系統(tǒng)通過蒸散和利用消耗水資源、截留等方式改變水的時間和空間分配特征,導(dǎo)致水文情勢發(fā)生變化,從而產(chǎn)生有利(如：調(diào)節(jié)洪峰、增加枯水季節(jié)徑流)或不利(如：加劇洪澇和干旱災(zāi)害)的影響。研究生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù)可為生態(tài)水文格局的合理構(gòu)建和水資源的持續(xù)利用提供理論支持[2]。

在生態(tài)系統(tǒng)和景觀尺度,許多學(xué)者圍繞不同植被類型和景觀格局對水文調(diào)節(jié)服務(wù)的影響開展了大量研究[3- 5],森林植被能夠影響水文過程、促進降雨再分配、影響土壤水分運動以及改變產(chǎn)流匯流條件等,從而緩和地表徑流,增加土壤徑流和地下徑流,在一定程度上起到了削峰補枯、控制土壤侵蝕、改善河流水質(zhì)等作用[6]。植被變化、農(nóng)業(yè)開發(fā)活動、道路建設(shè)以及城鎮(zhèn)化等土地利用變化也對水文調(diào)節(jié)服務(wù)產(chǎn)生顯著影響[7]。部分研究還認為：水利工程建設(shè)對削減洪峰、增加枯水季徑流發(fā)揮一定的作用[8- 10]。這些生態(tài)系統(tǒng)或景觀尺度上的研究對于認識生態(tài)系統(tǒng)水文調(diào)節(jié)過程或機理發(fā)揮了重要作用。在流域尺度上,生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù)受到生態(tài)系統(tǒng)特征、水利工程設(shè)施和人類用水等多個因素的綜合影響,并且這些影響存在尺度效應(yīng)[1,6-7],現(xiàn)階段,由于缺乏系統(tǒng)的流域尺度生態(tài)水循環(huán)過程的綜合觀測,從宏觀角度評估生態(tài)系統(tǒng)對區(qū)域水文過程的影響仍存在較大難度[11]。如何表征和認識流域尺度上生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù)特征、辨識不同影響因素對流域尺度生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù)的作用,仍然研究很少。

本研究以我國面積最大、降雨時空分布不均、洪澇災(zāi)害頻發(fā)的長江流域為對象,以平水年2013年為基礎(chǔ),基于降雨-地表徑流同步性特征研究子流域尺度上長江流域生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù)的空間格局以及影響因素,其目標(biāo)是：揭示長江流域生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù)的空間特征,闡明導(dǎo)致長江流域生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù)空間差異的主要因素及其相對貢獻,為長江流域生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù)管理提供基礎(chǔ)信息。

1 研究區(qū)域與方法

1.1 研究區(qū)概況

長江流域地處中國南部(24°—35°N,90°—122°E),全長6280 km,流域面積180萬km2(圖1)。流域受季風(fēng)性氣候影響降雨豐富但相對集中且時空分布不均,流域水源總量約為9960億m3,占全國35%以上[12],居全國首位,但人均占有量還處于較低水平。雨季容易出現(xiàn)洪澇災(zāi)害,平均約10年發(fā)生一次大洪水。如1998年長江流域發(fā)生特大洪水災(zāi)害；2007年,整個四川盆地在7 月經(jīng)歷半個月的大暴雨,100多萬人受災(zāi);2010年7月,長江流域境內(nèi)江蘇、江西、湖北、湖南、重慶、四川等多省市遭受暴雨洪澇災(zāi)害。在枯水季來水較少時很容易出現(xiàn)季節(jié)性缺水,如2006年川渝大旱部分縣持續(xù)干旱超過80 d,三分之二的溪河斷流,471座水庫干涸;2011年4月下旬至5月中旬,長江中下游地區(qū)出現(xiàn)重度以上氣象干旱,部分縣市一度達到特旱等級[13]。

圖1 長江流域水系及水文站分布圖Fig.1 Water system and hydrologic distribution in Yangtze River Basin

1.2 研究方法

1.2.1 水文調(diào)節(jié)服務(wù)表征方法

生態(tài)系統(tǒng)對降雨-徑流過程的所有影響都可以稱為水文調(diào)節(jié),例如降水-徑流關(guān)系、徑流過程的延遲、流域的蓄滲能力。本研究采用與江河徑流過程密切相關(guān)的變異系數(shù)法表征流域尺度上生態(tài)系統(tǒng)對水文過程的調(diào)節(jié)作用。其基本假設(shè)是：在無任何自然和人為干擾情況下,降雨過程與地表徑流過程具有高度同步性,降雨與河流徑流的年內(nèi)變異系數(shù)差異較小。但自然與人為影響會導(dǎo)致降雨與河流徑流的變異系數(shù)產(chǎn)生差異,二者差值(CVp-CVR)的正負號表明徑流年內(nèi)變化相對降雨更為平穩(wěn)或離散,二者差值(CVp-CVR)大小表明徑流相對降雨離散性變化的幅度,即水文調(diào)節(jié)的幅度。

1.2.2 水文調(diào)節(jié)功能影響因素選擇

人類活動對生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù)的影響主要通過改變下墊面類型(即：土地利用與覆蓋類型)、水利工程設(shè)施和社會經(jīng)濟發(fā)展的用水[14- 16]。本研究據(jù)此選擇自然植被(草地+灌叢+森林)覆蓋比例、農(nóng)田面積比例代表主要的土地覆蓋因子；選擇水庫庫容代表水利工程設(shè)施因子；選擇單位面積人口代表社會經(jīng)濟因子,探討不同因子對長江流域生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù)空間差異的相對影響。

1.2.3 數(shù)據(jù)收集與統(tǒng)計分析

子流域劃分：基于長江水利委員會對流域的劃分,選取每個子流域出水口水文站所控制的流域為研究區(qū),形成34個評價單元(圖1)。

降雨與徑流數(shù)據(jù)：為了使研究不受極端氣候的影響,參照中國水資源公報[17]選取平水年2013年作為研究時段。收集了2013年長江流域148個氣象臺站降雨數(shù)據(jù),采用普通克里格方法對流域降水進行空間插值,形成12個月面降雨量數(shù)據(jù)。

自然與社會數(shù)據(jù)：2010年生態(tài)系統(tǒng)數(shù)據(jù)來源于中國科學(xué)院遙感與數(shù)字地球研究所；水庫庫容數(shù)據(jù)來源于現(xiàn)有水利工程統(tǒng)計數(shù)據(jù)、人口數(shù)據(jù)來源于中國農(nóng)科院農(nóng)業(yè)信息研究所。

統(tǒng)計分析：采用主成分分析和線性回歸相結(jié)合的方法,建立自然地理因素(自然植被比例(X1))、人類活動(庫容(X2)、農(nóng)田比例(X3)、單位面積人口(X4))與水文調(diào)節(jié)服務(wù)回歸模型,分析自然地理因素和人類活動對水文調(diào)節(jié)服務(wù)的影響；對影響水文調(diào)節(jié)服務(wù)的4個因素進行主成分分析,提取主要成分,主成分需要反映原始指標(biāo)80%以上的變異信息。使用提取的主成分與因變量進行回歸,建立主成分與因變量回歸模型,然后將主成分還原成最初自變量得到初始自變量與因變量回歸模型(CVp-CVR=0.00100X1+0.19410X2-0.14506X3-0.00027X4-0.06439,R2=0.327,F=7.531)。

運用相對貢獻率的算法計算各影響因素的重要性大小,其主要思想是計算回歸模型中每添加一個因子引起R2變化的近似值,具體計算方法見參考文獻[18-19]。

2 結(jié)果與分析

2.1 降雨與徑流的空間分布

2013年是平水年,長江流域各地區(qū)年降雨量約為300—1700 mm。受季風(fēng)氣候和地形的影響,整個流域年降水量空間分布并不均勻,年降雨量由東南向西北遞減。降雨集中在東南部沿海地區(qū)及成都平原,金沙江上游、漢江上游降雨相對較少(圖2)。

長江東南區(qū)的流域降雨變異系數(shù)較西北區(qū)低,變異系數(shù)在0.83以下,即東南區(qū)的流域降雨年內(nèi)分配要比西北區(qū)的更加均勻。青衣江、岷江干流、涪江和嘉陵江干流降雨年內(nèi)分配極不均勻,變異系數(shù)大于1(圖2)。

長江屬于雨源性河流,徑流十分豐富。徑流深度大多在300—800 mm之間,分布在流域的東南部、雅礱江、大渡河和嘉陵江流域。徑流深低于300 mm的流域分布在金沙江、漢江和赤水河,屬于水資源相對匱乏的區(qū)域(圖3)。

徑流的年內(nèi)分布離散性空間變異較大,徑流變異系數(shù)變動范圍在0.26—1.65之間。下游的漢江、洞庭湖、鄱陽湖徑流年內(nèi)分布比較均勻,沱江和嘉陵江水系徑流變異系數(shù)較大,變異系數(shù)大于1,徑流年內(nèi)分布極不均勻(圖3)。

圖2 長江流域2013年降雨量分布與降雨年內(nèi)分配變異系數(shù)空間分布Fig.2 Rainfall distribution and rainfall coefficient of variation distribution in Yangtze River Basin

圖3 長江流域2013年徑流深度空間分布與徑流年內(nèi)分配變異系數(shù)空間分布Fig.3 Runoff distribution and runoff coefficient of variation distribution in Yangtze River Basin

2.2 水文調(diào)節(jié)服務(wù)空間特征

長江流域自西向東各子流域水文調(diào)節(jié)功能呈現(xiàn)出明顯的空間異質(zhì)性。長江流域源區(qū)和上游地區(qū)降雨和徑流變異系數(shù)差異較大,下游地區(qū)徑流相對降雨變幅較小,即水文調(diào)節(jié)服務(wù)主要作用地區(qū)在中上游地區(qū)(圖4)。通過對水文服務(wù)表征指標(biāo)聚類分析可以看出,水文調(diào)節(jié)服務(wù)可分為3類,一類是嘉陵江水系；第二類是金沙江以上的長江源區(qū)和漢江水系,剩下的流域形成第三類,且第三類與第一類更為相似(圖5)。

圖4 長江流域水文調(diào)節(jié)服務(wù)分布Fig.4 Hydrological regulating service distribution in Yangtze River Basin

2.3 水文調(diào)節(jié)服務(wù)的影響因素

自然地理要素和人類活動均對自流域尺度水文調(diào)節(jié)服務(wù)產(chǎn)生影響。水庫庫容和自然植被面積比例對水文調(diào)節(jié)服務(wù)具有正向促進作用,而農(nóng)田面積比例和單位面積人口對水文調(diào)節(jié)服務(wù)的作用正好相反。研究區(qū)域內(nèi)水庫庫容、自然植被面積比例、農(nóng)田面積比例和單位面積人口均對長江流域子流域水文調(diào)節(jié)服務(wù)的空間異質(zhì)性產(chǎn)生影響：水庫庫容對水文調(diào)節(jié)服務(wù)空間變異的貢獻率最大(58.55%)；其次是農(nóng)田面積比例(18.79%),自然植被面積比例和單位面積人口對水文調(diào)節(jié)服務(wù)空間變異的貢獻率均在10%左右(圖6)。

3 討論

研究區(qū)域內(nèi)表征流域尺度生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù)的CVp-CVR均值為0.085,但長江流域生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù)呈現(xiàn)出明顯的空間異質(zhì)性：水文調(diào)節(jié)服務(wù)最強的子流域位于長江源區(qū)和漢江水系,在水文調(diào)節(jié)作用下該子流域徑流年內(nèi)分配趨于平坦(CVp-CVR為0.477)；水文調(diào)節(jié)服務(wù)最弱的子流域位于嘉陵江水系(涪江、渠江和嘉陵江干流),降雨徑流的變異系數(shù)遠遠大于降水(CVp-CVR均值為-0.474)；長江中下游地區(qū)各子流域水文調(diào)節(jié)服務(wù)不明顯(CVp-CVR均值為0.033)(圖4,圖5)。

長江流域各子流域生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù)影響因素主要包括水庫庫容、農(nóng)田面積比例、自然植被面積比例和單位面積人口,但在不同子流域影響生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù)的主導(dǎo)因素各不相同。(1)第一類嘉陵江流域生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù)主要影響因素是農(nóng)田面積比例和單位面積人口,農(nóng)田面積比例高達39.80%,人口分布密集(280人/km2),導(dǎo)致水文調(diào)節(jié)服務(wù)弱。水源涵養(yǎng)功能是森林生態(tài)系統(tǒng)服務(wù)功能的重要組成部分,通過測評秦巴山區(qū)植被涵養(yǎng)水源總物質(zhì)量為33.82×108t/a,其中林冠截流量10.99×108t/a,枯枝落葉截流量3.20×108t/a,土壤截留量19.63×108t/a[20]。因而在嘉陵江流域,可以通過恢復(fù)植被提高水文調(diào)節(jié)服務(wù)作用。(2)第二類長江源區(qū)和漢江水系生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù)主要影響因素是水庫庫容和自然植被面積比例。水庫的調(diào)蓄功能能直接改變壩下河道的徑流量季節(jié)分配方式,丹江口水庫從1968年蓄水以后,汛期徑流量大幅減少,枯水季節(jié)流量增加,即使最枯的月份也能有20×108m3左右的水量[15]。長江源區(qū)和漢江流域水庫庫容(162.497×1012m3)遠遠超過了全流域平均水平,加上最高的自然植被覆蓋率(71.14%),因而水文調(diào)節(jié)服務(wù)作用最強。(3)第三類長江中下游地區(qū)各子流域較高的人口密度(230 人/km2)和農(nóng)田面積比例(24.10%)導(dǎo)致生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù)較弱。

本研究也還存在一些不足：變異系數(shù)只能代表降雨和徑流的離散程度,不能明確降雨和徑流極值出現(xiàn)的時間和大小,還需要其他輔助指標(biāo)才能進一步描述水文調(diào)節(jié)過程；此外,將庫容作為一個因子探討水電開發(fā)與流域水文調(diào)節(jié)服務(wù)的關(guān)系,但未考慮水庫不同調(diào)節(jié)方式對水文調(diào)節(jié)服務(wù)的影響,進一步結(jié)合水庫調(diào)節(jié)方式才能充分闡明水庫的水文調(diào)節(jié)服務(wù),例如：在嘉陵江流域,盡管水庫庫容量較大(67.2×1012m3),但水庫削峰補枯的作用發(fā)揮很小,原因一方面是渠江、沱江流域存在通過水庫進行流域間引水如引大濟岷,加大了徑流的變異,另一方面水庫的調(diào)節(jié)方式?jīng)Q定了水利工程的調(diào)節(jié)作用,如徑流式水庫對徑流的年內(nèi)分配影響較小。渠江、涪江7月份徑流量的比例遠遠大于降雨量的比例(圖7),表明水庫在此流域枯水期幾乎沒有調(diào)節(jié)作用可見,因此水庫的水文調(diào)節(jié)服務(wù)效用與水庫調(diào)節(jié)方式有密切關(guān)系；另外,將自然植被覆蓋比例作為一個因子沒有考慮不同植被類型的水源涵養(yǎng)能力差異性,例如本研究中自然植被包含森林、灌從、和草地三種土地利用模式,而相關(guān)研究表明這三種植被類型持水量有差異。張學(xué)龍等在祁連山寺大隆林區(qū)監(jiān)測數(shù)據(jù)分析結(jié)果顯示,各類林地的年平均持水量以灌叢林最高,為391.92 mm；其次是云杉林、圓柏林、分別為324.82 mm和298.25 mm；牧坡草地最少,為180.29 mm[21]。

圖7 嘉陵江水系月降雨量比例和月徑流比例圖Fig.7 The proportion of rainfall and runoff in Jialing River Basin

本研究在揭示長江流域子流域尺度上生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù)空間特征基礎(chǔ)上,研究結(jié)果也可為流域生態(tài)保護與恢復(fù)、生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù)提升提供科學(xué)依據(jù)。(1)在長江源區(qū)和漢江水系,水庫庫容和自然植被面積比例最高,導(dǎo)致生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù)最強。盡管研究表明水庫和自然植被對于水文調(diào)節(jié)服務(wù)都具有正向促進作用,但其生態(tài)學(xué)機制和效用有明顯差異。森林調(diào)節(jié)是自然調(diào)節(jié)過程,不會對河流連通性,魚類洄游和濕地候鳥棲息地產(chǎn)生不利影響,與水利工程調(diào)節(jié)相比,自然植被調(diào)節(jié)更加突出,需要更加重視。表明：在長江源區(qū)和漢江流域需要加強自然生態(tài)系統(tǒng)保護與恢復(fù),同時需要關(guān)注大壩建設(shè)的負面效應(yīng)。盡管水利工程建設(shè)增強了生態(tài)系統(tǒng)水文調(diào)節(jié)服務(wù),但一方面,由于水庫攔截,河流流速減慢,河水自凈能力被削弱,水質(zhì)遭到破壞。研究表明,攀枝花市在建庫前取水口的飲用水水源一級保護區(qū)內(nèi)COD和氨氮水質(zhì)均能達標(biāo),觀音巖電站和金沙、銀江電站建庫后,位于金沙庫區(qū)的精煤公司水廠和陶家渡水廠的飲用水水源COD有超標(biāo)現(xiàn)象[22]；另一方面,水壩建成后,破壞了河流生物原有棲息環(huán)境。上游河段由激流、緩流與回水塘等多樣化的流水型生境,演變成河面寬闊、深度變淺、水流緩慢的湖泊型水域,適合湍蛙等棲息的激流生境消失了,而適于緩流生活的大齒蟾和華西蟾蜍在水壩上游大量繁殖[23]。此外,水庫修建還改變了下游河道的流量過程。水庫不僅存蓄了汛期洪水,而且還截流了非汛期基流,往往會使下游河道水位大幅度下降甚至斷流[24]。據(jù)湖南省水文局2016年12月11日分析,長江入洞庭湖水量銳減9成,長江三口河系基本斷流,其中藕池河康家嘴站斷流達279 d[25]。(2)嘉陵江流域由于最高的農(nóng)田面積比例和最高人口密度,導(dǎo)致生態(tài)系統(tǒng)水文調(diào)節(jié)功能最弱。在這些流域內(nèi)需要加強自然植被的恢復(fù)、調(diào)控人類活動對生態(tài)系統(tǒng)水文調(diào)節(jié)功能的影響,提升生態(tài)系統(tǒng)水文調(diào)節(jié)功能。長江中下游流域也需通過上述生態(tài)保護與恢復(fù)措施增強生態(tài)系統(tǒng)水文調(diào)節(jié)功能。

[1] 呂一河, 胡健, 孫飛翔, 張立偉. 水源涵養(yǎng)與水文調(diào)節(jié): 和而不同的陸地生態(tài)系統(tǒng)水文服務(wù). 生態(tài)學(xué)報, 2015, 35(15): 5191- 5196.

[2] 黃奕龍, 傅伯杰, 陳利頂. 生態(tài)水文過程研究進展. 生態(tài)學(xué)報, 2003, 23(3): 580- 587.

[3] 楊國靖, 丁永建, 肖篤寧. 干旱區(qū)山地典型流域不同景觀結(jié)構(gòu)的水文調(diào)節(jié)能力分析. 冰川凍土, 2005, 27(5): 701- 708.

[4] 肖篤寧, 裴鐵凡, 趙羿. 遼河三角洲濕地景觀的水文調(diào)節(jié)與防洪功能. 濕地科學(xué), 2003, 1(1): 21- 25.

[5] 曹云. 紅壤丘陵區(qū)不同森林恢復(fù)類型水文調(diào)節(jié)功能的機制研究[D]. 北京: 中國科學(xué)院生態(tài)環(huán)境研究中心, 2007.

[6] 張彪, 李文華, 謝高地, 肖玉. 森林生態(tài)系統(tǒng)的水源涵養(yǎng)功能及其計量方法. 生態(tài)學(xué)雜志, 2009, 28(3): 529- 534.

[7] Bronstert A, Niehoff D, Bürger G. Effects of climate and land-use change on storm runoff generation: present knowledge and modelling capabilities. Hydrological Processes, 2002, 16(2): 509- 529.

[8] Yang Z S, Wang H J, Saito Y, Milliman J D, Xu K, Qiao S, Shi G. Dam impacts on the Changjiang (Yangtze) River sediment discharge to the sea: the past 55 years and after the Three Gorges Dam. Water Resources Research, 2006, 42(4): W04407.

[9] Yin H, Li C G. Human impact on floods and flood disasters on the Yangtze River. Geomorphology, 2001, 41(2/3): 105- 109.

[10] Zeilhofer P, de Moura R M. Hydrological changes in the northern Pantanal caused by the Manso dam: impact analysis and suggestions for mitigation. Ecological Engineering, 2009, 35(1): 105- 117.

[11] 王根緒, 李娜, 胡宏昌. 氣候變化對長江黃河源區(qū)生態(tài)系統(tǒng)的影響及其水文效應(yīng). 氣候變化研究進展, 2009, 5(4): 202- 208.

[12] 常福宣. 長江流域水資源配置的幾個關(guān)鍵問題思考. 長江科學(xué)院院報, 2011, 28(10): 54- 58.

[13] 陳鮮艷, 周兵, 鐘海玲, 柳艷菊. 2011年長江中下游春旱的氣候特征分析. 長江流域資源與環(huán)境, 2014, 23(1): 139- 145.

[14] 徐光來, 許有鵬, 徐宏亮. 城市化水文效應(yīng)研究進展. 自然資源學(xué)報, 2010, 25(12): 2171- 2178.

[15] 張俊宏. 丹江口水庫下游河道對漢江調(diào)水響應(yīng)機制及航道整治對策研究[D]. 武漢: 武漢大學(xué), 2014.

[16] 郭旭寧, 胡鐵松. 農(nóng)田流域水文響應(yīng)特征分析及模型考慮. 灌溉排水學(xué)報, 2010, 29(4): 26- 29.

[17] 中華人民共和國水利部. 2013年中國水資源公報. 北京: 中華人民共和國水利部, 2013.

[18] 郭呈全, 陳希鎮(zhèn). 主成分回歸的SPSS實現(xiàn). 統(tǒng)計與決策, 2011, (5): 157- 159.

[19] Kabacoff R I. R語言實戰(zhàn). 高濤, 肖楠, 陳鋼, 譯. 北京: 人民郵電出版社, 2013: 194- 196.

[20] 任志遠, 李晶. 陜南秦巴山區(qū)植被生態(tài)功能的價值測評. 地理學(xué)報, 2003, 58(4): 503- 511.

[21] 張學(xué)龍, 車克鈞, 王金葉, 閆文德, 張虎. 祁連山寺大隆林區(qū)土壤水分動態(tài)研究. 西北林學(xué)院學(xué)報, 1998, 13(1): 1- 9.

[22] 阮婭, 肖彩. 金沙江攀枝花河段水電開發(fā)對水質(zhì)影響的分析. 人民長江, 2009, 40(20): 81- 84.

[23] 李成, 顧海軍, 戴強, 劉志君, 王躍招. 草坡河流域小水電開發(fā)對無尾兩棲動物的影響. 長江流域資源與環(huán)境, 2008, 17(S1): 117- 121.

[24] 曹永強, 倪廣恒, 胡和平. 水利水電工程建設(shè)對生態(tài)環(huán)境的影響分析. 人民黃河, 2005, 27(1): 56- 58.

[25] 《長江流域資源與環(huán)境》編輯部. 長江入洞庭湖水量銳減9成. (2016- 12- 15) [2017-01-04]. http://changjiang.whlib.ac.cn/changjlyhp/dongth/201312/t20131226_7016.htm.