• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于PSR模型的烏梁素海生態(tài)系統(tǒng)健康分區(qū)評(píng)價(jià)*

    2013-09-25 03:05:24毛旭鋒崔麗娟張曼胤
    湖泊科學(xué) 2013年6期
    關(guān)鍵詞:烏梁素海湖區(qū)富營(yíng)養(yǎng)化

    毛旭鋒,崔麗娟,張曼胤

    (1:中國(guó)林業(yè)科學(xué)研究院濕地研究所,北京100091)

    (2:青海師范大學(xué)生命與地理科學(xué)學(xué)院,西寧810006)

    中國(guó)湖泊眾多,面積在1 km2以上的天然湖泊就有2800多個(gè).湖泊不僅提供了涵養(yǎng)水源、抵御洪水、調(diào)節(jié)氣候、凈化污染、野生動(dòng)植物棲息地等生態(tài)功能,而且提供了豐富的產(chǎn)品資源、歷史跡地和休閑觀光場(chǎng)所等社會(huì)服務(wù)功能,對(duì)于地區(qū)生態(tài)環(huán)境維持具有不可替代的作用[1].

    在自然和人類擾動(dòng)日益加劇、湖泊生態(tài)環(huán)境日益退化的背景下,定量分析湖泊生態(tài)健康特征,評(píng)價(jià)湖泊生態(tài)系統(tǒng)健康狀態(tài),對(duì)合理利用有限的湖泊資源、實(shí)現(xiàn)湖泊生態(tài)系統(tǒng)的科學(xué)管理和有效恢復(fù)具有重要的理論和現(xiàn)實(shí)意義,生態(tài)系統(tǒng)健康已經(jīng)成為環(huán)境管理的新目標(biāo)[2].早在1980s,Rapport等就提出用生態(tài)系統(tǒng)壓力指標(biāo)評(píng)估湖泊生態(tài)系統(tǒng)健康程度[3],而不同壓力狀態(tài)下湖泊的生產(chǎn)力和生物網(wǎng)絡(luò)的完整性也可以用來(lái)判斷湖泊的健康狀態(tài)[4-5].Constanza等提出的生態(tài)系統(tǒng)健康概念得到廣泛認(rèn)可[2],共涵蓋6個(gè)方面:自我平衡、沒有病征、多樣性、有恢復(fù)力、有活力和能夠保持系統(tǒng)組分間的平衡,不同學(xué)者提出的指標(biāo)多涵蓋于其中.Ul-anowicz從食物網(wǎng)絡(luò)結(jié)構(gòu)中營(yíng)養(yǎng)物質(zhì)流角度評(píng)價(jià)了湖泊生態(tài)系統(tǒng)健康[2].熱力學(xué)指標(biāo)也因簡(jiǎn)單快速的特點(diǎn)被廣泛和深入地應(yīng)用到湖泊生態(tài)健康評(píng)價(jià)[6],研究案例區(qū)包括巢湖等國(guó)內(nèi)外不同類型湖泊[7-10].隨著研究的深入,除生態(tài)系統(tǒng)以外的包括社會(huì)、經(jīng)濟(jì)和人口健康等方面的指標(biāo)也逐漸被關(guān)注[11].生態(tài)系統(tǒng)不同功能組的狀態(tài)也被用來(lái)指示湖泊生態(tài)系統(tǒng)的健康[12].在借鑒國(guó)外研究的基礎(chǔ)上,劉永等[13]從生態(tài)系統(tǒng)完整性和湖泊內(nèi)物質(zhì)循環(huán)等角度研究了湖泊生態(tài)系統(tǒng)的健康.胡志新等[14]用系統(tǒng)能、系統(tǒng)能結(jié)構(gòu)、生態(tài)緩沖容量和湖泊營(yíng)養(yǎng)狀態(tài)指數(shù)表征湖泊生態(tài)系統(tǒng)健康.由于指標(biāo)眾多,對(duì)現(xiàn)有的評(píng)價(jià)指標(biāo)的總結(jié)和分析研究也開展起來(lái)[15].

    在評(píng)價(jià)方法上,J?rgensen提出了一套初步評(píng)價(jià)方法[6];而生態(tài)系統(tǒng)健康指數(shù)法也被提出并被廣泛應(yīng)用[7-10,16].此外,人工神經(jīng)網(wǎng)絡(luò)方法[17]、層次分析法[18]等也不斷應(yīng)用到湖泊生態(tài)系統(tǒng)健康評(píng)估中.隨著研究的深入,人們發(fā)現(xiàn)生態(tài)系統(tǒng)健康并非反映在系統(tǒng)本身,湖泊系統(tǒng)的外在壓力下系統(tǒng)呈現(xiàn)的狀態(tài)以及人類對(duì)系統(tǒng)壓力的反應(yīng),都是衡量湖泊生態(tài)系統(tǒng)健康的重要因子.因此,一些學(xué)者將PSR模型應(yīng)用到生態(tài)系統(tǒng)健康評(píng)估中[19-20].然而,湖泊由于植被分布、地形變化、水文特征的差異,不同湖區(qū)健康狀態(tài)也會(huì)不同,尤其是面積較大的湖泊,區(qū)別更是明顯[21].用一個(gè)評(píng)價(jià)指標(biāo)反映整個(gè)湖區(qū)健康狀態(tài),并不能客觀、準(zhǔn)確地反映不同湖區(qū)的健康狀態(tài).因此,本研究在PSR模型的基礎(chǔ)上,以內(nèi)蒙古烏梁素海為例,開展湖泊生態(tài)系統(tǒng)健康分區(qū)評(píng)價(jià)研究,以期推動(dòng)湖泊生態(tài)系統(tǒng)健康評(píng)價(jià)體系和方法的完善,為湖泊生態(tài)系統(tǒng)科學(xué)有效的管理提供有效信息.

    1 研究區(qū)概況

    烏梁素海(40°36'~41°03'N,108°43'~108°57'E)位于內(nèi)蒙古自治區(qū)巴彥淖爾市烏拉特前旗境內(nèi)(圖1).烏梁素海現(xiàn)有水域面積292 km2,湖面高程海拔1018 m,湖水深度多數(shù)區(qū)域在0.5~2.5 m之間,多年平均水深為0.7 m.湖泊呈南北長(zhǎng)、東西窄的狹長(zhǎng)形態(tài),其中南北長(zhǎng)35~40 km,東西寬5~10 km,蓄水2.5×108~3.0×108m3.它是我國(guó)第八大淡水湖,也是黃河流域同緯度最大的淡水湖,對(duì)改善周邊地區(qū)氣候、氣象條件、增加大氣降水、維持當(dāng)?shù)氐纳鷳B(tài)平衡有著極為重要的作用.

    近幾十年來(lái),隨著上游工、農(nóng)業(yè)生產(chǎn)的發(fā)展和城鎮(zhèn)人口的增加,每年都有相當(dāng)數(shù)量的農(nóng)田退水、工業(yè)污水和城鎮(zhèn)生活污水通過(guò)總排干進(jìn)入湖中,使烏梁素海的水質(zhì)受到了嚴(yán)重污染、水體富營(yíng)養(yǎng)化嚴(yán)重、生態(tài)系統(tǒng)退化迅速,嚴(yán)重影響到當(dāng)?shù)氐纳鷳B(tài)安全.由于湖泊主要的水源為排干排入的上游農(nóng)田退稅、工業(yè)廢水和生活污水,由湖區(qū)南部出口排出.所以湖區(qū)由北向南的水質(zhì)污染程度和富營(yíng)養(yǎng)化程度是逐漸降低的,不同湖區(qū)生態(tài)系統(tǒng)健康呈現(xiàn)不同的水平.根據(jù)烏梁素海環(huán)境污染的特征和消減規(guī)律[22],將湖區(qū)劃分為Ⅰ、Ⅱ和Ⅲ3個(gè)區(qū).Ⅰ區(qū)為湖區(qū)西北部入水口附近,Ⅱ區(qū)位于湖區(qū)中部,Ⅲ區(qū)為湖區(qū)西南部出水口附近區(qū)域.開展烏梁素海濕地生態(tài)健康分區(qū)評(píng)價(jià)研究,以期找出其湖泊健康退化的主要原因,并提出有效保護(hù)、恢復(fù)及合理利用烏梁素海濕地的措施.

    2 評(píng)價(jià)指標(biāo)與方法

    2.1 評(píng)價(jià)指標(biāo)體系

    綜合考慮烏梁素海社會(huì)、經(jīng)濟(jì)和環(huán)境現(xiàn)狀,從湖泊的壓力、狀態(tài)和響應(yīng)3個(gè)角度選取湖泊健康評(píng)價(jià)指標(biāo),確定生態(tài)環(huán)境特征、生態(tài)功能和社會(huì)環(huán)境3類不同但又相互聯(lián)系的25項(xiàng)指標(biāo)及其量化標(biāo)準(zhǔn),構(gòu)建出適合該地區(qū)湖泊生態(tài)系統(tǒng)健康評(píng)價(jià)的指標(biāo)體系.壓力指標(biāo)反映人類活動(dòng)給濕地生態(tài)系統(tǒng)所帶來(lái)的影響與脅迫,反映資源的利用強(qiáng)度和變化趨勢(shì);狀態(tài)指標(biāo)反映濕地生態(tài)環(huán)境要素的變化,體現(xiàn)環(huán)境保護(hù)政策的最終目標(biāo);響應(yīng)指標(biāo)反映了社會(huì)或個(gè)人為停止、減輕、預(yù)防或恢復(fù)不利于濕地生態(tài)系統(tǒng)變化采取的措施.該指標(biāo)體系由目標(biāo)層、子系統(tǒng)層和評(píng)價(jià)指標(biāo)層構(gòu)成.詳細(xì)的指標(biāo)如表1所示.

    2.2 評(píng)價(jià)方法

    2.2.1 指標(biāo)權(quán)重確定方法 采用層次分析法(AHP)定性、定量結(jié)合,借助專家經(jīng)驗(yàn)判斷定量化各指標(biāo)權(quán)重.層次分析法權(quán)重確定過(guò)程如下,首先請(qǐng)專家們?cè)u(píng)價(jià)指標(biāo)的相對(duì)重要性.將初步結(jié)果經(jīng)數(shù)學(xué)驗(yàn)證和處理后的結(jié)果反饋給數(shù)位專家,請(qǐng)專家確認(rèn),經(jīng)多次征求意見形成判斷矩陣.例如,子系統(tǒng)B1、B2和B3的最終相對(duì)重要性判斷矩陣為A.如B2相對(duì)B1的重要性均值為3,則B1相對(duì)B2的重要性均值為1/3;B2相對(duì)B3的重要性均值為5,則B3相對(duì)B2的重要性均值為1/5.

    表1 烏梁素海湖泊健康診斷指標(biāo)Tab.1 Health diagnosis indices of the Lake Wuliangsuhai

    將對(duì)矩陣A進(jìn)行歸一化(矩陣元素除以每一列元素之和),得到歸一化的矩陣B.

    將矩陣B各行相加結(jié)果分別為3/13+5/23+3/9=0.774、9/13+15/23+5/9=1.911、1/13+3/23+1/9=0.315.矩陣 B 中,λmax=3;CI=0;RI=0.58;CR=0 <0.10.其中,CR 為判斷矩陣的隨機(jī)一致性比率,CI為判斷矩陣一致性指標(biāo),RI為判斷矩陣的平均隨機(jī)一致性指標(biāo).子系統(tǒng)的權(quán)重為B中行元素之和與λmax的比值,即0.774/3=0.258、1.911/3=0.637 和0.315/3=0.105,即權(quán)重矩陣為 W=(0.258 0.637 0.105)T.其他指標(biāo)的權(quán)重同樣采用上述方法,剔除無(wú)效結(jié)果后經(jīng)多次驗(yàn)證確定.

    2.2.2 湖泊綜合評(píng)價(jià)方法 本研究采用基于PSR模型的綜合評(píng)價(jià)方法[23].該模型最初由加拿大統(tǒng)計(jì)學(xué)家Rapport和Friend提出,后由經(jīng)濟(jì)合作與發(fā)展組織(OECD)和聯(lián)合國(guó)環(huán)境規(guī)劃署(UNEP)于20世紀(jì)八九十年代共同發(fā)展起來(lái)的用于研究環(huán)境問(wèn)題的框架體系.該方法從社會(huì)經(jīng)濟(jì)與環(huán)境有機(jī)統(tǒng)一的觀點(diǎn)出發(fā),精確地反映了生態(tài)系統(tǒng)健康的自然、經(jīng)濟(jì)、社會(huì)因素間的關(guān)系,綜合考慮了環(huán)境受到的壓力和環(huán)境退化之間的因果關(guān)系,壓力、狀態(tài)、響應(yīng)3個(gè)環(huán)節(jié)相互制約、相互影響,以及環(huán)境決策和制定對(duì)策措施的全過(guò)程.湖泊生態(tài)系統(tǒng)的健康狀態(tài),采用綜合健康指數(shù)(Comprehensive Health Index,CHI)定量評(píng)價(jià):

    式中,CHI為湖泊生態(tài)系統(tǒng)綜合健康指數(shù);Wi為各指標(biāo)的綜合權(quán)重;Xi為各指標(biāo)得分.按照綜合指數(shù)由高到低將湖泊生態(tài)系統(tǒng)健康分為5 個(gè)等級(jí):①[0.8,1.0],良好狀態(tài);②[0.6,0.8],較好狀態(tài);③[0.4,0.6],警戒狀態(tài);④[0.2,0.4],較差狀態(tài);⑤[0,0.2],極差狀態(tài).

    其他評(píng)價(jià)方法包括湖泊富營(yíng)養(yǎng)化程度評(píng)價(jià),主要采用綜合營(yíng)養(yǎng)狀態(tài)指數(shù)(Trophic Level Index,TLI)表示[25]:

    式中,TLI(∑)為綜合營(yíng)養(yǎng)狀態(tài)指數(shù);Wj為第j種參數(shù)的營(yíng)養(yǎng)狀態(tài)指數(shù)的相關(guān)權(quán)重;TLI(j)代表第j種參數(shù)的營(yíng)養(yǎng)狀態(tài)指數(shù).依照TLI(∑)指數(shù)大小,可將營(yíng)養(yǎng)狀態(tài)劃分為貧營(yíng)養(yǎng)(TLI(∑)<30)、中營(yíng)養(yǎng)(30≤TLI(∑)<50)、輕度富營(yíng)養(yǎng)(50≤TLI(∑)<60)、中度富營(yíng)養(yǎng)(60≤TLI(∑)≤70)和重度富營(yíng)養(yǎng)(TLI(∑)>70).參考相關(guān)資料,Chl.a、TP、TN、SD 和 CODMn等指標(biāo)權(quán)重分別為 0.2663、0.2237、0.2183、0.2210 和 0.2210[24].

    2.3 數(shù)據(jù)來(lái)源

    水質(zhì)類別、富營(yíng)養(yǎng)化程度、湖泊淤積程度、土壤有機(jī)質(zhì)、濕地植物、初級(jí)生產(chǎn)力、水質(zhì)凈化功能等數(shù)據(jù)來(lái)自于課題組2011-2012年監(jiān)測(cè)和調(diào)查數(shù)據(jù),監(jiān)測(cè)點(diǎn)遍布全湖20個(gè)地區(qū).水質(zhì)類別參照《中國(guó)人民共和國(guó)地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB/T 3838-2002)要求,對(duì)不同水體指標(biāo)進(jìn)行監(jiān)測(cè)評(píng)估,水質(zhì)分類數(shù)據(jù)主要來(lái)自烏拉特后旗環(huán)境保護(hù)監(jiān)測(cè)站的定期監(jiān)測(cè)結(jié)果.烏拉特后旗濕地管理水平和公共環(huán)境意識(shí)等數(shù)據(jù)來(lái)自2011年對(duì)游客、當(dāng)?shù)鼐用竦某闃诱{(diào)查以及網(wǎng)絡(luò)調(diào)查數(shù)據(jù),共計(jì)456份.魚類豐富度、優(yōu)勢(shì)種覆蓋度、補(bǔ)水保證率、珍稀鳥類種類、棲息地維持功能、物質(zhì)生產(chǎn)功能、水文調(diào)節(jié)功能數(shù)據(jù)主要來(lái)自烏梁素海漁場(chǎng)和烏梁素海鳥類管理站,其他數(shù)據(jù)主要來(lái)源于內(nèi)蒙古河套灌區(qū)管理總局[25]和《巴彥淖爾統(tǒng)計(jì)年鑒》等.

    3 結(jié)果與討論

    3.1 各指標(biāo)權(quán)重結(jié)果

    基于專家打分法構(gòu)建判斷矩陣,在層次分析法的基礎(chǔ)上,分析不同子系統(tǒng)中各指標(biāo)的權(quán)重.其中,壓力指標(biāo)、狀態(tài)指標(biāo)和響應(yīng)指標(biāo)的權(quán)重分別為0.258、0.637和0.105,其中狀態(tài)指標(biāo)中生態(tài)狀態(tài)、環(huán)境狀態(tài)和功能狀態(tài)具有同等重要性,權(quán)重均為0.333,其他指標(biāo)層指標(biāo)的權(quán)重和綜合權(quán)重如表2所示.

    表2 各評(píng)價(jià)指標(biāo)權(quán)重Tab.2 Weights of various assessment indices

    根據(jù)專家打分的結(jié)果,狀態(tài)指標(biāo)是反映湖泊健康的最重要的指標(biāo),壓力指標(biāo)次之,響應(yīng)指標(biāo)最低.壓力指標(biāo)中農(nóng)藥使用強(qiáng)度(C3)和化肥使用強(qiáng)度(C4)具有最高的權(quán)重,反映了湖泊高污染和高營(yíng)養(yǎng)狀態(tài).相應(yīng)指標(biāo)中污水處理指數(shù)(C21)的權(quán)重最大.狀態(tài)層又進(jìn)步細(xì)化劃分為生態(tài)狀態(tài)、環(huán)境狀態(tài)和功能狀態(tài),根據(jù)專家打分結(jié)果,三者具有相同的權(quán)重.而凈初級(jí)生產(chǎn)力(C8)、富營(yíng)養(yǎng)化程度(C12)和水質(zhì)凈化功能(C18)分別是3類狀態(tài)指標(biāo)中權(quán)重最大的評(píng)價(jià)指標(biāo).

    3.2 各評(píng)價(jià)指標(biāo)綜合得分

    由于指標(biāo)體系中的各項(xiàng)評(píng)價(jià)指標(biāo)的類型復(fù)雜,各系數(shù)之間的量綱不統(tǒng)一,各指標(biāo)之間缺乏可比性.為了便于比較和計(jì)算,依據(jù)當(dāng)?shù)亟?jīng)濟(jì)、社會(huì)和自然環(huán)境特點(diǎn),并參考其他研究結(jié)果[26-27],按指標(biāo)實(shí)際數(shù)值等級(jí)劃分為5級(jí),標(biāo)準(zhǔn)化分值設(shè)定在0~1.0之間(0~0.2代表極差,0.2 ~0.4代表較差,0.4~0.6 代表警戒,0.6 ~0.8代表較好,0.8~1.0代表良好,主要參考其他數(shù)據(jù)和歷史數(shù)據(jù)確定),每項(xiàng)指標(biāo)得分見表3.

    表3 各評(píng)價(jià)指標(biāo)得分及綜合得分Tab.3 Scores of different assessment indices

    從子系統(tǒng)層來(lái)看,壓力指標(biāo)層中得分最高和最低的分別為C5(萬(wàn)元GDP COD排放強(qiáng)度)和C2(水資源開發(fā)利用率)指標(biāo),但綜合得分最高和最低的分別為C4(化肥使用強(qiáng)度)和C2,說(shuō)明指標(biāo)得分與綜合得分趨勢(shì)并不一致.狀態(tài)指標(biāo)中綜合得分較高的指標(biāo)包括C8(初級(jí)生產(chǎn)力)、C12(富營(yíng)養(yǎng)化程度)、C18(水質(zhì)凈化功能)和C19(棲息地維持功能),而綜合得分較低的指標(biāo)包括C9(魚類豐富度)、C13(土壤有機(jī)質(zhì))、C14(補(bǔ)水保證率)、C17(水文調(diào)節(jié)功能)和C20(科考旅游功能).值得一提的是,雖然烏梁素海湖區(qū)為富營(yíng)養(yǎng)化湖泊,但富營(yíng)養(yǎng)化指數(shù)獲得較高的得分,說(shuō)明不能僅通過(guò)分析富營(yíng)養(yǎng)化狀態(tài)判斷湖泊的健康狀況.與壓力指標(biāo)層得分類似,指標(biāo)得分和綜合得分情況并不完全一致.響應(yīng)指標(biāo)中綜合得分最高和最低的指標(biāo)分別為C21(廢水處理指數(shù))和C22(公共環(huán)境意識(shí)).對(duì)第Ⅲ區(qū)不同指標(biāo)進(jìn)行聚類分析,劃分為2~4類群集結(jié)果(表4).

    對(duì)指標(biāo)的聚類分析結(jié)果顯示,按照指標(biāo)對(duì)生態(tài)系統(tǒng)健康的貢獻(xiàn)值不同,可將指標(biāo)劃分為4類.從分析結(jié)果來(lái)看,C8、C12和C18是反映該湖泊Ⅲ區(qū)健康的第4類指標(biāo),而在Ⅰ區(qū)中的第4類指標(biāo)為C12,Ⅲ區(qū)中則為C12和C18.從不同湖區(qū)得分來(lái)看,3個(gè)湖區(qū)生態(tài)系統(tǒng)健康指標(biāo)在壓力層和響應(yīng)層的得分上沒有區(qū)分,差異主要體現(xiàn)在狀態(tài)層的指標(biāo)上,包括優(yōu)勢(shì)種覆蓋度(C10)、水質(zhì)類別(C11)、富營(yíng)養(yǎng)化程度(C12)、湖泊淤積程度(C15)、物質(zhì)生產(chǎn)功能(C16)和水質(zhì)凈化功能(C18).造成這些指標(biāo)得分差異的主要原因是湖區(qū)地形、水文特征等差異引起水質(zhì)、水生植被種類和分布、水生動(dòng)物分布的差異,進(jìn)而使得不同湖區(qū)在物質(zhì)生產(chǎn)功能、泥沙淤積程度、富營(yíng)養(yǎng)化程度呈現(xiàn)不同的狀態(tài).同一指標(biāo)Ⅰ區(qū)和Ⅲ區(qū)的差異最大達(dá)到5倍.指標(biāo)呈現(xiàn)的基本規(guī)律是,隨著水體由北向南流動(dòng),污染物濃度逐漸減小,湖泊泥沙的淤積程度和富營(yíng)養(yǎng)化程度均逐漸降低,物質(zhì)生產(chǎn)功能逐漸升高,優(yōu)勢(shì)種的覆蓋度也逐漸增大,湖泊的水質(zhì)由入口處的Ⅴ類甚至劣Ⅴ類水升至出口處Ⅲ類水,水質(zhì)凈化功能逐漸提升.因此,3個(gè)湖區(qū)的健康綜合指數(shù)基本上呈現(xiàn)遞增的趨勢(shì).

    表4 聚類分析結(jié)果Tab.4 Results of cluster analysis

    3.3 烏梁素海綜合健康指數(shù)

    3個(gè)湖區(qū)的綜合健康指數(shù)分別為0.3395、0.3866和0.4494(圖1a),均處于警戒狀態(tài)或者較差狀態(tài).這說(shuō)明烏梁素海生態(tài)系統(tǒng)健康已經(jīng)出現(xiàn)嚴(yán)重的問(wèn)題,亟需采取相關(guān)的措施改善湖區(qū)的生態(tài)環(huán)境,恢復(fù)烏梁素海的生態(tài)功能.從分區(qū)來(lái)看,健康指數(shù)分別是Ⅲ區(qū)>Ⅱ區(qū)>Ⅰ區(qū)(圖1a),說(shuō)明隨著湖泊對(duì)污染物的消減和凈化過(guò)程,湖泊生態(tài)系統(tǒng)健康指數(shù)逐漸增加,同時(shí),這也說(shuō)明不同湖區(qū)的生態(tài)系統(tǒng)健康狀況不盡相同,對(duì)不同地區(qū)應(yīng)實(shí)施不同的湖泊恢復(fù)措施.

    3個(gè)區(qū)域各指標(biāo)對(duì)綜合健康指數(shù)的貢獻(xiàn)率都是狀態(tài)指標(biāo)>壓力指標(biāo)>響應(yīng)指標(biāo)(圖1b、c、d).從不同湖區(qū)來(lái)看,生態(tài)系統(tǒng)的狀態(tài)指標(biāo)的得分逐漸升高,對(duì)健康指數(shù)的貢獻(xiàn)率由67%上升到70%、最后升至74%;而壓力指標(biāo)和響應(yīng)指標(biāo)對(duì)健康指數(shù)的貢獻(xiàn)率是逐漸降的,二者在Ⅰ、Ⅱ和Ⅲ區(qū)的貢獻(xiàn)率分別表現(xiàn)為30%和12%、26%和10%、22%和9%.由上文可知,3類指標(biāo)的權(quán)重分別為0.2578、0.6359和0.3081,狀態(tài)指標(biāo)的權(quán)重明顯大于其他兩類指標(biāo),這也是造成上述結(jié)果的主要原因.然而,從3類指標(biāo)得分占各自權(quán)重的比例(得分率)來(lái)看,Ⅰ、Ⅱ和Ⅲ區(qū)的壓力指標(biāo)和響應(yīng)指標(biāo)得分率均為39%,而狀態(tài)指標(biāo)在3個(gè)區(qū)域的得分率分別為31%、37%和47%.

    總體而言,3類指標(biāo)的得分率均低于50%,說(shuō)明在湖泊健康的恢復(fù)上,還需要加大源頭治理的力度,進(jìn)一步減少農(nóng)藥化肥的使用強(qiáng)度;同時(shí),應(yīng)提升入湖污水處理率,增加環(huán)保的投資.分區(qū)對(duì)比可以發(fā)現(xiàn),對(duì)Ⅰ區(qū)和Ⅱ區(qū)而言,壓力指標(biāo)和響應(yīng)指標(biāo)的得分率略大于狀態(tài)指標(biāo)得分率;而在Ⅲ區(qū),壓力指標(biāo)和響應(yīng)指標(biāo)的得分率要小于狀態(tài)指標(biāo)的得分率.因此,針對(duì)不同湖區(qū)的不同情況,還應(yīng)采取不同的措施恢復(fù)湖泊的健康.例如,對(duì)于污染和富營(yíng)養(yǎng)化較重的Ⅰ區(qū),應(yīng)注重提高生態(tài)狀態(tài)和環(huán)境狀態(tài)指標(biāo)的得分,如采用效果較快的物理和化學(xué)開展水質(zhì)恢復(fù);而對(duì)于Ⅱ區(qū),應(yīng)注重提高其功能狀態(tài)指標(biāo)的得分,如科研旅游功能、物質(zhì)生產(chǎn)功能和棲息地維持功能等.

    3.4 敏感性分析

    PSR模型具有很強(qiáng)的因果關(guān)系,可作為濕地生態(tài)系統(tǒng)健康評(píng)價(jià)的指標(biāo)分類體系,但它未能對(duì)所有指標(biāo)進(jìn)行嚴(yán)格的分類,也不能判斷不同指標(biāo)之間的相互關(guān)系和影響,只是對(duì)相關(guān)狀態(tài)的綜合評(píng)價(jià)結(jié)果.因此,對(duì)模型的不確定性進(jìn)行分析,測(cè)試模型結(jié)果對(duì)指標(biāo)變化的反映,分析不同指標(biāo)的敏感性,是分析檢驗(yàn)?zāi)P涂茖W(xué)性的重要途徑之一.通過(guò)調(diào)整系數(shù)K(-10%,+10%,-20%,+20%),采用公式ΔK)核算敏感性,相關(guān)結(jié)果如圖2所示.

    圖1 烏梁素海健康評(píng)價(jià)結(jié)果Fig.1 The CHI of the Lake Wuliangsuhai

    壓力指標(biāo)變化 ±10%,綜合健康指數(shù)變化 ±2.7%;狀態(tài)指標(biāo)變化 ±10%,綜合健康指數(shù)變化約±6.2%;而響應(yīng)指標(biāo)變化±10%,綜合健康指數(shù)變化±1.0%左右.當(dāng)壓力指標(biāo)變化±20%,綜合健康指數(shù)變化±5.5%;狀態(tài)指標(biāo)變化±20%,綜合健康指數(shù)變化±12.4%左右;響應(yīng)指標(biāo)變化±20%,綜合健康指數(shù)變化±2.1%左右,模型敏感性排序?yàn)闋顟B(tài)指標(biāo)>壓力指標(biāo)>響應(yīng)指標(biāo)(圖2).總體而言,模型對(duì)指標(biāo)變化的反映較為敏感,說(shuō)明本模型具有較好的適用性.

    從敏感性分析結(jié)果來(lái)看,烏梁素海綜合健康指數(shù)受狀態(tài)指標(biāo)的影響最大,這也與烏梁素海當(dāng)前富營(yíng)養(yǎng)化嚴(yán)重、湖泊生態(tài)系統(tǒng)的生態(tài)環(huán)境受到破壞、生態(tài)系統(tǒng)功能退化有關(guān)[22,28].狀態(tài)指標(biāo)與生態(tài)系統(tǒng)初級(jí)生產(chǎn)、水質(zhì)類別、富營(yíng)養(yǎng)化程度、水文調(diào)節(jié)、水質(zhì)凈化等狀態(tài)密切相關(guān),這與其他研究成果一致,即強(qiáng)調(diào)對(duì)生態(tài)系統(tǒng)水質(zhì)、生態(tài)和水文的修復(fù),恢復(fù)湖泊調(diào)節(jié)、供給和支撐功能.因此,首先需改善水質(zhì)、保證水量充足并提高湖泊整體聯(lián)通性和流動(dòng)性[28].此外,水華暴發(fā)對(duì)湖泊狀態(tài)指標(biāo)的影響巨大,其與營(yíng)養(yǎng)鹽的空間分布具有很強(qiáng)的相關(guān)性[29].由于營(yíng)養(yǎng)鹽呈現(xiàn)北高南低、西高東低的規(guī)律,加上該區(qū)域湖體較為封閉,水動(dòng)力條件較為穩(wěn)定,導(dǎo)致烏梁素海黃苔集中在湖區(qū)北部和西岸排干口附近區(qū)域暴發(fā)[30],這也導(dǎo)致第Ⅰ湖區(qū)的綜合健康指數(shù)最低.隨著水流由北向南的流動(dòng)過(guò)程,營(yíng)養(yǎng)鹽濃度逐漸降低,富營(yíng)養(yǎng)化導(dǎo)致的水華現(xiàn)象減少,湖泊健康指數(shù)在第Ⅱ區(qū)有所升高,但仍應(yīng)加強(qiáng)該區(qū)域底泥疏浚和水體疏通,保持湖泊連通性是提高該區(qū)域湖泊健康指數(shù)的重要措施[31].在水質(zhì)相對(duì)改善的第Ⅲ區(qū)域,應(yīng)注重其生態(tài)狀態(tài)和功能狀態(tài)指標(biāo)的恢復(fù).提高水生態(tài)系統(tǒng)中的生物多樣性,如采用投放魚苗等方式恢復(fù)魚類多樣性.同時(shí),應(yīng)盡量減少人為影響和破壞,尤其應(yīng)控制游船數(shù)量和旅游范圍,減少對(duì)動(dòng)物棲息地干擾,提高棲息地維持功能的得分.

    4 結(jié)論

    烏梁素海作為干旱區(qū)典型的湖泊濕地,其生態(tài)系統(tǒng)具有敏感性和脆弱性的特點(diǎn),基于PSR模型的生態(tài)系統(tǒng)健康評(píng)價(jià),能客觀描述湖泊生態(tài)生態(tài)系統(tǒng)退化的原因和程度,為科學(xué)有效地恢復(fù)湖泊提供依據(jù).尤其是對(duì)湖泊健康的分區(qū)評(píng)價(jià),能根據(jù)湖泊不同區(qū)域理化特征和生態(tài)差異評(píng)價(jià)湖泊健康,更具針對(duì)性和準(zhǔn)確性.通過(guò)本次研究,初步得到以下幾點(diǎn)結(jié)論:

    1)烏梁素海的綜合健康指數(shù)得分較低,濕地整體正處于警戒狀態(tài)或者較差狀態(tài),說(shuō)明該湖泊正受到較大的外界壓力,湖泊生態(tài)系統(tǒng)結(jié)構(gòu)和功能正處于較差的狀態(tài).

    圖2 Ⅰ、Ⅱ、Ⅲ區(qū)敏感性分析結(jié)果Fig.2 Results of sensitive analysis of three regions in the Lake Wuliangsuhai

    2)烏梁素海不同湖區(qū)的綜合健康指數(shù)存在較為明顯的區(qū)別,總體表現(xiàn)為Ⅲ區(qū)>Ⅱ區(qū)>Ⅰ區(qū),而各區(qū)域中不同指標(biāo)的得分情況依次為狀態(tài)指標(biāo)>壓力指標(biāo)>響應(yīng)指標(biāo).

    3)對(duì)Ⅰ區(qū)和Ⅱ區(qū)而言,壓力指標(biāo)和響應(yīng)指標(biāo)的得分率略大于狀態(tài)指標(biāo)得分率;而在Ⅲ區(qū),壓力指標(biāo)和響應(yīng)指標(biāo)的得分率要小于狀態(tài)指標(biāo)的得分率.針對(duì)不同湖區(qū)的不同情況,應(yīng)關(guān)注不同的指標(biāo)得分,采取不同的措施恢復(fù)湖泊的健康.

    致謝:感謝內(nèi)蒙古烏拉特梭梭林—蒙古野驢國(guó)家級(jí)自然保護(hù)區(qū)管理局、烏梁素海保護(hù)區(qū)管理站在野外監(jiān)測(cè)、資料共享等方面提供的幫助.

    [1]Mao XF,Yang ZF.Functional assessment of associated aquatic ecosystems in the Baiyangdian Basin:an ecological network analysis based approach.Ecological Modelling,2011,222(23/24):3811-3820.

    [2]Costanza R,Norton BG,Hasbell BD.Ecosystem health:New goals for environmental management.Washington DC:Island Press,1992:190-206,239-256.

    [3]Rapport DJ,Regier HA,Hutchinson TC.Ecosystem behavior under stress.American Naturalist,1985,116:326-346.

    [4]Hannon B.Ecosystem flow analysis.Canadian Bulletin of Fisheries and Aquatic Sciences,1985,213:97-118.

    [5]Karr JR,F(xiàn)ausch KD,Angermeier PL et al.Assessing biological integrity in running waters:A method and its rationale.Champaign:Illinois Natural History Survey,Special Publication 5,1986.

    [6]J?rgensen SE.Exergy and ecological buffer capacitiesas measures of ecosystem health.Ecosystem Health,1995,1(3):150-160.

    [7]Xu FL.Ecosystem health assessment for Lake Chao,a shallow eutrophic Chinese lake.Lakes& Reservoirs:Research&Management,1996,2:101-109.

    [8]Xu FL,J?rgensen SE,Tao S.Ecological indicators for assessing fresh water ecosystem health.Ecological Modeling,1999,116(1):77-106.

    [9]Xu FL,Tao S,Dawson PW.Lake ecosystem health assessment:indicators and methods.Water Research,2001,35(13):3157-3167.

    [10]Xu FL,Dawson RW,Tao S.A method for lake ecosystem health assessment:an ecological modelling method and its application.Hydrobiologica,2001,443(1/2/3):159-175.

    [11]Shear H.The development and use of indicators to assess ecosystem health state in the Great Lake.Ecosystem Health,1996,2(2):241-258.

    [12]Connor RJ,Walls TE,Hughes RM.Using multiple taxonomic groups to index the ecological condition of lakes.Environmental Monitoring and Assessment,2000,61:207-228.

    [13]劉 永,郭懷成,戴永立.湖泊生態(tài)系統(tǒng)健康評(píng)價(jià)方法研究.環(huán)境科學(xué)學(xué)報(bào),2004,24(4):723-729.

    [14]胡志新,胡維平,谷孝鴻等.太湖湖泊生態(tài)系統(tǒng)健康評(píng)價(jià).湖泊科學(xué),2005,17(3):256-262.

    [15]崔保山,楊志峰.濕地生態(tài)系統(tǒng)健康評(píng)價(jià)指標(biāo)體Ⅱ:方法與案例.生態(tài)學(xué)報(bào),2002,22(8):1231-1239.

    [16]趙臻彥,徐福留,詹 巍等.湖泊生態(tài)系統(tǒng)健康定量評(píng)價(jià)方法.生態(tài)學(xué)報(bào),2005,25(6):1466-1474.

    [17]Mo MH,Wang XL,Wu HJ.Ecosystem health assessment of Honghu Lake wetland of China using Artificial Neural Network Approach.Chinese Geographical Science,2009,19(4):125-129.

    [18]Ma LG,Cao YR,Li XT.Assessment of ecosystem health in Lashihai Lake using AHP method.Journal of Geo-Information Science,2011,13(2):234-239.

    [19]許文杰,曹升樂(lè).基于PSR熵權(quán)綜合健康指數(shù)法的城市湖泊生態(tài)系統(tǒng)健康評(píng)價(jià).水文,2010,30(5):64-68.

    [20]Chen WG.Study on urban lake ecosystem health assessment.Applied Mechanics and Material,2012,178-181:304-308.

    [21]胡志新,胡維平,陳永根等.太湖不同湖區(qū)生態(tài)系統(tǒng)健康評(píng)價(jià)方法研究.農(nóng)村生態(tài)環(huán)境,2005,21(4):28-32.

    [22]毛旭鋒,崔麗娟,李勝男等.烏梁素海水體營(yíng)養(yǎng)鹽削減的空間變化規(guī)律研究.干旱區(qū)資源與環(huán)境,2012,11:116-121.

    [23]Rapport DJ,F(xiàn)riend AM.Towards a comprehensive framework for environmental statistics:a stress-response approach.Statistics Canada,Ottawa,1979:11-510.

    [24]金相燦.中國(guó)湖泊環(huán)境.北京:海洋出版社,1995.

    [25]內(nèi)蒙古河套灌區(qū)管理總局.淺析巴彥淖爾市水資源開發(fā)與利用.黃河流域部分城市水環(huán)境論壇論文集,2008.

    [26]顏 利,王金坑,黃 浩.基于PSR框架模型的東溪流域生態(tài)系統(tǒng)健康評(píng)價(jià).資源科學(xué),2008,30(1):107-112.

    [27]許文杰,許士國(guó).湖泊生態(tài)系統(tǒng)健康評(píng)價(jià)的熵權(quán)綜合健康指數(shù)法.水土保持研究,2008,15(1):125-127.

    [28]姜忠峰,李暢游,張 生等.改進(jìn)AHP法在烏梁素海生態(tài)系統(tǒng)服務(wù)功能評(píng)價(jià)中的應(yīng)用.干旱區(qū)資源與環(huán)境,2011,25(1):135-139.

    [29]王 艷,黃永梅,于長(zhǎng)水等.基于MODIS數(shù)據(jù)的2000-2010年烏梁素?!包S苔”時(shí)空變化.湖泊科學(xué),2012,24(4):519-527.

    [30]楊志巖,李暢游,張 生等.烏梁素海氮磷濃度與葉綠素a時(shí)空分布關(guān)系研究.干旱環(huán)境監(jiān)測(cè),2009,(1):16-19.

    [31]厚福祥,鄧 芳.內(nèi)蒙古烏梁素海生態(tài)恢復(fù)工程試驗(yàn)研究.內(nèi)蒙古大學(xué)學(xué)報(bào):自然科學(xué)版,2006,37(1):105-110.

    猜你喜歡
    烏梁素海湖區(qū)富營(yíng)養(yǎng)化
    基于臨界點(diǎn)的杭州灣水體富營(yíng)養(yǎng)化多年變化研究
    大通湖區(qū)河蟹產(chǎn)業(yè)發(fā)展綜述
    洪口水庫(kù)近年富營(yíng)養(yǎng)化程度時(shí)間分布的研究
    生活在湖區(qū)
    海峽旅游(2018年4期)2018-06-01 11:20:00
    湖區(qū)航道風(fēng)速預(yù)警監(jiān)測(cè)點(diǎn)布設(shè)研究
    江西建材(2018年4期)2018-04-10 12:37:24
    十字路口的烏梁素海
    洞庭湖典型垸內(nèi)溝渠水體富營(yíng)養(yǎng)化評(píng)價(jià)
    烏梁素海地區(qū)近50年氣候變化特征分析
    富營(yíng)養(yǎng)化藻的特性與水熱液化成油的研究
    風(fēng)景、地方與浪漫主義——華茲華斯的湖區(qū)書寫
    天天躁夜夜躁狠狠躁躁| 电影成人av| 日韩一卡2卡3卡4卡2021年| 我的亚洲天堂| 亚洲精品在线美女| 亚洲精品一二三| 免费女性裸体啪啪无遮挡网站| 午夜老司机福利片| 亚洲色图综合在线观看| 少妇 在线观看| 亚洲国产看品久久| 一夜夜www| 一边摸一边抽搐一进一小说 | 露出奶头的视频| 久久久水蜜桃国产精品网| 久久精品人人爽人人爽视色| 日本av免费视频播放| 99国产精品一区二区蜜桃av | 免费一级毛片在线播放高清视频 | 国产97色在线日韩免费| 欧美 亚洲 国产 日韩一| 欧美 日韩 精品 国产| 18禁国产床啪视频网站| 后天国语完整版免费观看| 亚洲精品国产色婷婷电影| 美女国产高潮福利片在线看| 精品视频人人做人人爽| 国产免费福利视频在线观看| 国产成人免费观看mmmm| 精品久久久精品久久久| 在线观看一区二区三区激情| 欧美在线一区亚洲| 欧美大码av| 五月天丁香电影| 亚洲精品国产区一区二| 老司机靠b影院| 变态另类成人亚洲欧美熟女 | 久久久欧美国产精品| 美女午夜性视频免费| 日韩欧美一区二区三区在线观看 | 国产视频一区二区在线看| 欧美+亚洲+日韩+国产| 亚洲精品久久成人aⅴ小说| 91精品三级在线观看| 极品教师在线免费播放| 午夜福利在线观看吧| 日韩欧美国产一区二区入口| 色综合欧美亚洲国产小说| 亚洲avbb在线观看| 国产三级黄色录像| 国产一卡二卡三卡精品| 成年人免费黄色播放视频| 国产伦人伦偷精品视频| 超色免费av| 两人在一起打扑克的视频| 人人妻人人添人人爽欧美一区卜| 亚洲精品在线美女| 国产午夜精品久久久久久| 麻豆国产av国片精品| 丰满迷人的少妇在线观看| 欧美成人午夜精品| 国产伦理片在线播放av一区| 香蕉国产在线看| 国产一区二区 视频在线| 亚洲欧美激情在线| 一二三四社区在线视频社区8| 人成视频在线观看免费观看| 欧美大码av| 丰满饥渴人妻一区二区三| 久久av网站| 超色免费av| 天堂8中文在线网| 婷婷成人精品国产| 国产成人精品在线电影| 亚洲综合色网址| 成年人午夜在线观看视频| 日韩精品免费视频一区二区三区| 成年人黄色毛片网站| 日本av手机在线免费观看| 12—13女人毛片做爰片一| 国产在线视频一区二区| 国产极品粉嫩免费观看在线| 色综合婷婷激情| 涩涩av久久男人的天堂| 一级片免费观看大全| 黑人猛操日本美女一级片| 亚洲欧美日韩另类电影网站| 国产欧美亚洲国产| av片东京热男人的天堂| 黄色成人免费大全| 国产欧美亚洲国产| 777米奇影视久久| 九色亚洲精品在线播放| 99国产精品免费福利视频| 麻豆国产av国片精品| 午夜福利视频精品| av天堂在线播放| 国产精品av久久久久免费| 国产麻豆69| 午夜福利影视在线免费观看| 精品一区二区三卡| 久久午夜亚洲精品久久| 一区二区三区乱码不卡18| 亚洲少妇的诱惑av| 亚洲精华国产精华精| av天堂久久9| 超碰97精品在线观看| 19禁男女啪啪无遮挡网站| 国产麻豆69| av不卡在线播放| 精品视频人人做人人爽| 免费观看a级毛片全部| 丝袜美足系列| 国产xxxxx性猛交| 久久免费观看电影| 黄片播放在线免费| 精品亚洲成国产av| 99riav亚洲国产免费| 精品国内亚洲2022精品成人 | 免费在线观看影片大全网站| 肉色欧美久久久久久久蜜桃| 久久热在线av| 免费看a级黄色片| 18禁美女被吸乳视频| √禁漫天堂资源中文www| 黄网站色视频无遮挡免费观看| 亚洲熟妇熟女久久| 亚洲综合色网址| 国产精品亚洲一级av第二区| 岛国毛片在线播放| 亚洲精品国产一区二区精华液| 欧美成人午夜精品| 国内毛片毛片毛片毛片毛片| 久久午夜综合久久蜜桃| 大型黄色视频在线免费观看| 中亚洲国语对白在线视频| 黄色a级毛片大全视频| 伊人久久大香线蕉亚洲五| 国产成人影院久久av| 久久久精品94久久精品| 老熟妇仑乱视频hdxx| 叶爱在线成人免费视频播放| 黑人欧美特级aaaaaa片| 亚洲精品在线观看二区| 视频在线观看一区二区三区| netflix在线观看网站| 亚洲国产看品久久| 国产精品美女特级片免费视频播放器 | 久久久久久免费高清国产稀缺| av有码第一页| 丰满人妻熟妇乱又伦精品不卡| 久久婷婷成人综合色麻豆| 黑丝袜美女国产一区| 亚洲第一欧美日韩一区二区三区 | videosex国产| 久久久精品94久久精品| 午夜免费成人在线视频| 亚洲精华国产精华精| 精品国内亚洲2022精品成人 | a级片在线免费高清观看视频| 满18在线观看网站| 亚洲第一欧美日韩一区二区三区 | 一本一本久久a久久精品综合妖精| 亚洲五月婷婷丁香| 中文字幕精品免费在线观看视频| 亚洲国产看品久久| 国产免费福利视频在线观看| 国产精品影院久久| 国产无遮挡羞羞视频在线观看| 国产日韩欧美在线精品| 肉色欧美久久久久久久蜜桃| 热99re8久久精品国产| 人人妻人人添人人爽欧美一区卜| 欧美日韩黄片免| 新久久久久国产一级毛片| 国产一卡二卡三卡精品| 51午夜福利影视在线观看| 十分钟在线观看高清视频www| 亚洲色图av天堂| 国产欧美日韩精品亚洲av| 国产精品秋霞免费鲁丝片| 中文字幕精品免费在线观看视频| 国产一区二区三区综合在线观看| 欧美另类亚洲清纯唯美| 亚洲伊人色综图| 老鸭窝网址在线观看| 亚洲精品av麻豆狂野| 日韩成人在线观看一区二区三区| 国产黄频视频在线观看| 香蕉久久夜色| 亚洲精品美女久久久久99蜜臀| 18禁裸乳无遮挡动漫免费视频| 黑人巨大精品欧美一区二区mp4| 成人免费观看视频高清| 成人18禁高潮啪啪吃奶动态图| 韩国精品一区二区三区| 成人手机av| 一区二区三区乱码不卡18| 亚洲精品一二三| 免费不卡黄色视频| 蜜桃国产av成人99| 久久亚洲真实| 国产成+人综合+亚洲专区| 国产主播在线观看一区二区| av福利片在线| 日韩中文字幕欧美一区二区| 69av精品久久久久久 | 久久久欧美国产精品| 久久这里只有精品19| 成年版毛片免费区| 叶爱在线成人免费视频播放| 啦啦啦在线免费观看视频4| 搡老岳熟女国产| 如日韩欧美国产精品一区二区三区| 咕卡用的链子| 亚洲熟妇熟女久久| 久久午夜综合久久蜜桃| 日韩精品免费视频一区二区三区| 久久久久久亚洲精品国产蜜桃av| 女人被躁到高潮嗷嗷叫费观| 汤姆久久久久久久影院中文字幕| 亚洲午夜精品一区,二区,三区| 黄色视频不卡| 亚洲欧洲日产国产| 91成人精品电影| 一区二区三区精品91| 中文字幕精品免费在线观看视频| 久久久国产精品麻豆| 久久久精品94久久精品| 国产真人三级小视频在线观看| 波多野结衣一区麻豆| 亚洲一区中文字幕在线| 极品教师在线免费播放| 91精品国产国语对白视频| 国产主播在线观看一区二区| 狂野欧美激情性xxxx| 久久久久久免费高清国产稀缺| 亚洲av第一区精品v没综合| 人人妻人人添人人爽欧美一区卜| 一进一出好大好爽视频| 亚洲avbb在线观看| 国产成人系列免费观看| 亚洲人成伊人成综合网2020| 悠悠久久av| 纵有疾风起免费观看全集完整版| 咕卡用的链子| 欧美激情极品国产一区二区三区| 久久热在线av| 久久久欧美国产精品| 18禁裸乳无遮挡动漫免费视频| 婷婷成人精品国产| 桃红色精品国产亚洲av| 在线观看66精品国产| 丝袜美腿诱惑在线| 日本av手机在线免费观看| 亚洲精品久久成人aⅴ小说| 菩萨蛮人人尽说江南好唐韦庄| 天堂俺去俺来也www色官网| 狠狠狠狠99中文字幕| 国精品久久久久久国模美| 在线观看免费视频网站a站| 日本vs欧美在线观看视频| 成人特级黄色片久久久久久久 | 亚洲av片天天在线观看| 欧美黄色片欧美黄色片| 欧美人与性动交α欧美精品济南到| 精品亚洲成国产av| 精品人妻熟女毛片av久久网站| 久久免费观看电影| 国产一区二区 视频在线| 一级a爱视频在线免费观看| videos熟女内射| 交换朋友夫妻互换小说| www.熟女人妻精品国产| 亚洲精品中文字幕在线视频| 一二三四社区在线视频社区8| 午夜精品久久久久久毛片777| 亚洲av片天天在线观看| 国产精品影院久久| 啦啦啦 在线观看视频| 首页视频小说图片口味搜索| 欧美另类亚洲清纯唯美| 欧美精品亚洲一区二区| 国产一卡二卡三卡精品| 国产熟女午夜一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 日韩 欧美 亚洲 中文字幕| 久久九九热精品免费| 国产99久久九九免费精品| 极品少妇高潮喷水抽搐| 欧美精品人与动牲交sv欧美| 国产在线精品亚洲第一网站| 亚洲专区国产一区二区| 亚洲第一av免费看| 久久精品国产亚洲av高清一级| 免费av中文字幕在线| 丰满饥渴人妻一区二区三| 国产av精品麻豆| 脱女人内裤的视频| 捣出白浆h1v1| 久久午夜综合久久蜜桃| 宅男免费午夜| 久久久久精品国产欧美久久久| tube8黄色片| 精品一区二区三区四区五区乱码| 自线自在国产av| 亚洲熟女精品中文字幕| 久热爱精品视频在线9| 久久人妻福利社区极品人妻图片| 咕卡用的链子| 日韩熟女老妇一区二区性免费视频| 久久人人爽av亚洲精品天堂| 日韩欧美三级三区| 欧美精品av麻豆av| 久久毛片免费看一区二区三区| 大码成人一级视频| 国产精品麻豆人妻色哟哟久久| aaaaa片日本免费| 少妇粗大呻吟视频| 国产精品成人在线| 亚洲精品国产区一区二| 麻豆国产av国片精品| 久久久精品区二区三区| 欧美乱妇无乱码| 女同久久另类99精品国产91| a在线观看视频网站| 亚洲精品一卡2卡三卡4卡5卡| 国产亚洲精品久久久久5区| 色播在线永久视频| 日韩视频一区二区在线观看| 精品久久久精品久久久| 免费观看av网站的网址| 不卡av一区二区三区| 香蕉丝袜av| 咕卡用的链子| 久久人妻熟女aⅴ| 咕卡用的链子| 日韩欧美三级三区| 一级片免费观看大全| videosex国产| 一级片'在线观看视频| 女同久久另类99精品国产91| 美女午夜性视频免费| 久久精品国产亚洲av高清一级| 亚洲色图综合在线观看| 国产精品久久久久久精品古装| 久久中文字幕人妻熟女| 午夜福利影视在线免费观看| 51午夜福利影视在线观看| 夜夜骑夜夜射夜夜干| 精品国产一区二区三区久久久樱花| 精品少妇内射三级| 天天添夜夜摸| 一本综合久久免费| 黄频高清免费视频| 成年人午夜在线观看视频| 在线看a的网站| 色婷婷久久久亚洲欧美| 久久久久久久久久久久大奶| 黄色成人免费大全| 黄片大片在线免费观看| 欧美日韩黄片免| 两性夫妻黄色片| 妹子高潮喷水视频| 99香蕉大伊视频| 久久久久久久大尺度免费视频| 亚洲熟女精品中文字幕| 亚洲全国av大片| 国产极品粉嫩免费观看在线| 99国产极品粉嫩在线观看| 老司机午夜十八禁免费视频| 成人国产一区最新在线观看| 99热网站在线观看| 亚洲欧美色中文字幕在线| 精品久久蜜臀av无| 男女无遮挡免费网站观看| 国产精品久久电影中文字幕 | 午夜福利视频在线观看免费| 大码成人一级视频| 一区二区日韩欧美中文字幕| 国产aⅴ精品一区二区三区波| 久久毛片免费看一区二区三区| 伦理电影免费视频| 亚洲第一青青草原| 国产精品电影一区二区三区 | 三级毛片av免费| 亚洲av成人不卡在线观看播放网| 亚洲免费av在线视频| 久久久国产欧美日韩av| 人妻久久中文字幕网| 国产精品香港三级国产av潘金莲| 一级毛片电影观看| 99精品欧美一区二区三区四区| 天堂8中文在线网| 日本黄色日本黄色录像| 色播在线永久视频| 亚洲欧美一区二区三区黑人| 国产精品免费视频内射| tocl精华| 成人国产一区最新在线观看| 国产成人欧美在线观看 | av超薄肉色丝袜交足视频| 久热爱精品视频在线9| 超色免费av| 国产无遮挡羞羞视频在线观看| 69精品国产乱码久久久| 嫩草影视91久久| e午夜精品久久久久久久| 午夜福利视频在线观看免费| 一级片'在线观看视频| 久久精品国产亚洲av香蕉五月 | 丰满饥渴人妻一区二区三| 亚洲熟女精品中文字幕| 人人妻,人人澡人人爽秒播| 十分钟在线观看高清视频www| 亚洲欧美色中文字幕在线| 国产精品免费大片| 18禁国产床啪视频网站| 国产成人免费无遮挡视频| 90打野战视频偷拍视频| 欧美精品啪啪一区二区三区| 最近最新免费中文字幕在线| 999精品在线视频| 99久久国产精品久久久| 免费少妇av软件| 丰满人妻熟妇乱又伦精品不卡| 日韩欧美三级三区| 中文字幕高清在线视频| 国产午夜精品久久久久久| 欧美精品av麻豆av| 日本精品一区二区三区蜜桃| 热re99久久国产66热| 一二三四社区在线视频社区8| √禁漫天堂资源中文www| 亚洲视频免费观看视频| 亚洲色图综合在线观看| 精品国产乱码久久久久久男人| 国产精品国产高清国产av | 亚洲情色 制服丝袜| 捣出白浆h1v1| 一本久久精品| 中亚洲国语对白在线视频| 国产无遮挡羞羞视频在线观看| 99久久国产精品久久久| 水蜜桃什么品种好| 久久性视频一级片| 人人妻,人人澡人人爽秒播| 久久久久国产一级毛片高清牌| 欧美变态另类bdsm刘玥| 久久99热这里只频精品6学生| 欧美日韩中文字幕国产精品一区二区三区 | av欧美777| 丰满迷人的少妇在线观看| 国产三级黄色录像| 两性午夜刺激爽爽歪歪视频在线观看 | 搡老岳熟女国产| 女人精品久久久久毛片| videos熟女内射| 久久这里只有精品19| 男人舔女人的私密视频| 18禁黄网站禁片午夜丰满| 亚洲专区中文字幕在线| 国产欧美日韩综合在线一区二区| 丝袜喷水一区| 国产精品一区二区在线不卡| 国产一区二区激情短视频| 一本大道久久a久久精品| 夜夜爽天天搞| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品乱久久久久久| 国产不卡av网站在线观看| 99精品在免费线老司机午夜| 丝袜美腿诱惑在线| 黑人巨大精品欧美一区二区mp4| 美女主播在线视频| 免费少妇av软件| 制服人妻中文乱码| 免费高清在线观看日韩| 国产精品一区二区精品视频观看| 美女主播在线视频| 久久国产精品男人的天堂亚洲| 欧美大码av| 久久精品人人爽人人爽视色| 欧美日韩av久久| 欧美日韩视频精品一区| tube8黄色片| 欧美大码av| 色94色欧美一区二区| 我要看黄色一级片免费的| 这个男人来自地球电影免费观看| 91字幕亚洲| 久久久国产一区二区| 天堂动漫精品| 女人爽到高潮嗷嗷叫在线视频| 午夜福利在线观看吧| 日韩免费av在线播放| 婷婷丁香在线五月| 高清视频免费观看一区二区| 亚洲欧美一区二区三区久久| 视频在线观看一区二区三区| 日韩欧美国产一区二区入口| 黑人猛操日本美女一级片| 国产午夜精品久久久久久| 欧美性长视频在线观看| 成年人黄色毛片网站| 亚洲美女黄片视频| 欧美久久黑人一区二区| 99精国产麻豆久久婷婷| 精品少妇一区二区三区视频日本电影| 一区二区av电影网| 国产高清激情床上av| 亚洲精品av麻豆狂野| 国产成人精品久久二区二区91| 亚洲精品中文字幕在线视频| 亚洲第一欧美日韩一区二区三区 | 黑人猛操日本美女一级片| 美女高潮到喷水免费观看| a级毛片黄视频| 在线观看免费日韩欧美大片| 亚洲精品国产区一区二| 精品国产乱码久久久久久小说| 亚洲国产成人一精品久久久| 一区在线观看完整版| 成人特级黄色片久久久久久久 | 精品欧美一区二区三区在线| 国产精品影院久久| 女人精品久久久久毛片| 波多野结衣av一区二区av| 久久ye,这里只有精品| 国产男女超爽视频在线观看| 欧美成人午夜精品| 国产97色在线日韩免费| 深夜精品福利| 在线十欧美十亚洲十日本专区| 久久九九热精品免费| 日韩成人在线观看一区二区三区| 黄色a级毛片大全视频| 一夜夜www| 老司机靠b影院| av网站在线播放免费| 欧美日本中文国产一区发布| 曰老女人黄片| 18禁黄网站禁片午夜丰满| 久久精品国产亚洲av香蕉五月 | 极品少妇高潮喷水抽搐| 成人18禁高潮啪啪吃奶动态图| 久久中文字幕一级| 亚洲精品国产区一区二| 欧美日韩亚洲国产一区二区在线观看 | 日韩精品免费视频一区二区三区| 国产免费视频播放在线视频| 一级片免费观看大全| 99久久人妻综合| 制服人妻中文乱码| 久久人妻熟女aⅴ| 亚洲国产欧美日韩在线播放| 老司机福利观看| 国产亚洲欧美在线一区二区| 亚洲欧美精品综合一区二区三区| 精品国产一区二区久久| 久久久欧美国产精品| 国产在视频线精品| 最新的欧美精品一区二区| 色在线成人网| 精品国内亚洲2022精品成人 | 日韩有码中文字幕| 在线永久观看黄色视频| 黑人操中国人逼视频| 成人手机av| 亚洲国产看品久久| 色综合婷婷激情| 黄片大片在线免费观看| 巨乳人妻的诱惑在线观看| 国产精品电影一区二区三区 | 制服人妻中文乱码| 精品免费久久久久久久清纯 | 一个人免费在线观看的高清视频| 欧美成人午夜精品| 狂野欧美激情性xxxx| 岛国毛片在线播放| 精品少妇一区二区三区视频日本电影| 欧美黄色淫秽网站| 国产一区二区三区综合在线观看| 国产精品99久久99久久久不卡| 性少妇av在线| 91精品三级在线观看| 精品福利观看| e午夜精品久久久久久久| 十八禁网站网址无遮挡| 亚洲精品一卡2卡三卡4卡5卡| e午夜精品久久久久久久| 一二三四社区在线视频社区8| 欧美日韩av久久| 两人在一起打扑克的视频| 国产欧美日韩精品亚洲av| 国产精品1区2区在线观看. | 91大片在线观看| 亚洲av第一区精品v没综合| 91字幕亚洲| 国产精品一区二区免费欧美| 人成视频在线观看免费观看| 日本五十路高清| 精品国产乱子伦一区二区三区| 日韩人妻精品一区2区三区| 午夜福利视频在线观看免费| 日日爽夜夜爽网站| 超色免费av| 日韩精品免费视频一区二区三区| 丰满少妇做爰视频| 久久香蕉激情| 91麻豆精品激情在线观看国产 | 国产亚洲精品一区二区www | 久久国产亚洲av麻豆专区| 欧美日韩成人在线一区二区| 亚洲国产看品久久|