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

    基于PSR的黃河河口區(qū)生態(tài)系統(tǒng)健康評價

    2017-04-10 12:28:09牛明香徐賓鐸
    生態(tài)學(xué)報 2017年3期
    關(guān)鍵詞:河口區(qū)河口健康狀況

    牛明香, 王 俊, 徐賓鐸

    1 中國水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所, 青島 266071 2 青島海洋科學(xué)與技術(shù)國家實(shí)驗(yàn)室海洋生態(tài)與環(huán)境科學(xué)功能實(shí)驗(yàn)室,青島 266200 3 農(nóng)業(yè)部海洋漁業(yè)可持續(xù)發(fā)展重點(diǎn)實(shí)驗(yàn)室,青島 266071 4 山東省漁業(yè)資源與生態(tài)環(huán)境重點(diǎn)實(shí)驗(yàn)室, 青島 266071 5 中國海洋大學(xué)水產(chǎn)學(xué)院,青島 266003

    基于PSR的黃河河口區(qū)生態(tài)系統(tǒng)健康評價

    牛明香1,3,4, 王 俊1,2,3,4,*, 徐賓鐸5

    1 中國水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所, 青島 266071 2 青島海洋科學(xué)與技術(shù)國家實(shí)驗(yàn)室海洋生態(tài)與環(huán)境科學(xué)功能實(shí)驗(yàn)室,青島 266200 3 農(nóng)業(yè)部海洋漁業(yè)可持續(xù)發(fā)展重點(diǎn)實(shí)驗(yàn)室,青島 266071 4 山東省漁業(yè)資源與生態(tài)環(huán)境重點(diǎn)實(shí)驗(yàn)室, 青島 266071 5 中國海洋大學(xué)水產(chǎn)學(xué)院,青島 266003

    根據(jù)壓力-狀態(tài)-響應(yīng)(PSR)框架模型,從廣義上定義河口區(qū)生態(tài)系統(tǒng),將河口及毗鄰的陸域、海域生態(tài)系統(tǒng)作為一個整體,從壓力指標(biāo)、狀態(tài)指標(biāo)、響應(yīng)指標(biāo)3個方面構(gòu)建了黃河河口區(qū)生態(tài)系統(tǒng)健康評價的指標(biāo)體系,以研究區(qū)1991年數(shù)據(jù)和相關(guān)國家標(biāo)準(zhǔn)為基準(zhǔn),2013年代表現(xiàn)況,利用綜合指數(shù)法(CEI)評價了黃河河口區(qū)的生態(tài)系統(tǒng)健康狀況。結(jié)果顯示:黃河河口區(qū)生態(tài)系統(tǒng)健康評價的響應(yīng)指數(shù)最高(0.9055),壓力指數(shù)居中(0.8288),狀態(tài)指數(shù)最低(0.6458),綜合指數(shù)為0.7427。總體來看,與1991年相比,目前黃河河口區(qū)生態(tài)系統(tǒng)仍處于“健康”狀態(tài),但健康狀況明顯下降,其中狀態(tài)指數(shù)下降最為嚴(yán)重。從區(qū)域輕度開發(fā)到人類活動強(qiáng)烈干擾階段,黃河河口區(qū)存在過度捕撈、濕地不合理開發(fā)、淺海養(yǎng)殖過度及污染物排放等一系列影響生態(tài)系統(tǒng)健康的問題,應(yīng)進(jìn)行區(qū)域的生態(tài)恢復(fù)和科學(xué)管理。

    生態(tài)系統(tǒng);健康評價;PSR模型;綜合評價指數(shù);黃河口

    河口是流域和海洋物質(zhì)交換的主要通道,由于陸、海交互作用及強(qiáng)烈的人類干擾活動,河口區(qū)成為環(huán)境變化的敏感帶和生態(tài)系統(tǒng)的脆弱帶。河口生態(tài)系統(tǒng)健康評價作為區(qū)域生態(tài)環(huán)境評價的一種新方法,能夠客觀反映河口區(qū)生態(tài)系統(tǒng)的總體特征,并探索生態(tài)系統(tǒng)各個層次存在的問題和隱患[1],從而采取相應(yīng)的策略,實(shí)現(xiàn)生態(tài)系統(tǒng)的健康和可持續(xù)發(fā)展。該方法已廣泛應(yīng)用于河口生態(tài)系統(tǒng)環(huán)境研究[2-4]。

    由于河口處在流域的最下游,來自流域、河口海岸帶和海洋的自然壓力和人類活動壓力共同影響著河口生態(tài)系統(tǒng)的健康狀況[5],僅利用河口收集的數(shù)據(jù)進(jìn)行評價對于理解河口生態(tài)系統(tǒng)健康是受限制的,應(yīng)在流域水平上考慮各因素之間的交互影響[6]。黃河口是我國著名的大河河口,但已開展的生態(tài)系統(tǒng)健康評價分別集中于黃河三角洲濕地[7-9]、萊州灣及黃河口水域[10-11],而大空間尺度的綜合性健康評價尚未開展。為了全面理解河口生態(tài)系統(tǒng)健康狀況,本研究選取反映河口及其毗鄰陸域、海域生態(tài)特征的一系列指標(biāo),基于壓力-狀態(tài)-響應(yīng)(PSR)模型建立生態(tài)系統(tǒng)健康評價指標(biāo)體系,利用綜合指數(shù)法(comprehensive evaluation index, CEI)對黃河河口區(qū)生態(tài)系統(tǒng)健康狀況進(jìn)行評價,以期為區(qū)域科學(xué)管理和可持續(xù)發(fā)展提供依據(jù)。

    1 研究區(qū)域與數(shù)據(jù)來源

    1.1 研究區(qū)概況

    圖1 研究區(qū)地理位置示意圖Fig.1 Sketch map of the study area

    河口生態(tài)系統(tǒng)健康受流域、海岸帶及海洋的綜合影響,而流域下游的各種生態(tài)指標(biāo)是對上游生態(tài)環(huán)境的一定反映。因此,本研究中的黃河河口區(qū)是指與黃河口相鄰的、對其生態(tài)環(huán)境最具影響的黃河三角洲濕地、黃河河口及一定范圍內(nèi)的相鄰海域作為一個整體進(jìn)行研究?;趨^(qū)劃的完整性、統(tǒng)計(jì)資料的一致性及海上調(diào)查范圍,以東營市河口區(qū)、東營區(qū)和墾利縣的行政區(qū)域、黃河河口及相鄰海域?yàn)榈湫脱芯繀^(qū)域(圖1)。

    黃河口區(qū)域處于山東半島和遼東半島環(huán)抱的中心地帶,處于京津唐經(jīng)濟(jì)區(qū)與沿黃經(jīng)濟(jì)帶的交匯點(diǎn)。研究區(qū)所屬的黃河三角洲是國際重要濕地,也是世界上暖溫帶保存最廣闊、最完整、最年輕的濕地生態(tài)系統(tǒng);黃河三角洲國家級自然保護(hù)區(qū),是以保護(hù)新生濕地生態(tài)系統(tǒng)和珍稀瀕危鳥類為主的濕地類型自然保護(hù)區(qū)。毗鄰的渤海是多種經(jīng)濟(jì)魚蝦類的產(chǎn)卵場和索餌場,在黃渤海漁業(yè)生產(chǎn)上曾占有極其重要的地位。

    20世紀(jì)80年代以來,黃河斷流、入海河流水質(zhì)污染、三角洲開發(fā)、捕撈強(qiáng)度增大及海水養(yǎng)殖面積迅速擴(kuò)增等多種因素的影響,黃河口及鄰近海域生態(tài)環(huán)境面臨嚴(yán)重威脅。

    1.2 數(shù)據(jù)來源

    研究區(qū)內(nèi),1990年設(shè)立黃河三角洲自然保護(hù)區(qū),1993年黃河三角洲開發(fā)進(jìn)入山東省委決策。本研究選擇1991年和2013年為研究時段,分別代表黃河口區(qū)域輕度開發(fā)及人類活動強(qiáng)烈干擾階段。根據(jù)研究區(qū)的植被特征、物候期特點(diǎn)以及海上調(diào)查時間,以10月份為基準(zhǔn)進(jìn)行數(shù)據(jù)選擇。

    遙感數(shù)據(jù)有美國Landsat 5 TM和Landsat 8 OLI 影像數(shù)據(jù)各一景,成像時間分別為1991年9月23日和2013年10月5日,軌道號為121/34。

    實(shí)地調(diào)查數(shù)據(jù)包括2013年10月河口及鄰近海域的海洋生物資源數(shù)據(jù)、物理化學(xué)等環(huán)境數(shù)據(jù);研究區(qū)1991年相同時段內(nèi)無海上調(diào)查航次,因此海洋生物資源歷史數(shù)據(jù)選用最接近的1992年8月份。

    統(tǒng)計(jì)數(shù)據(jù):1991年、2013年東營市統(tǒng)計(jì)年鑒;1991年、2013年《黃河水沙公報》;2013年《黃河水資源公報》、2013年《中國海洋災(zāi)害公報》。

    其它資料:《海水水質(zhì)標(biāo)準(zhǔn)GB 3097—1997》、中華人民共和國環(huán)境保護(hù)部2014年1月頒發(fā)的《國家生態(tài)文明建設(shè)示范村鎮(zhèn)指標(biāo)(試行)》及公開發(fā)表的文獻(xiàn)資料。

    2 研究方法

    2.1 評價指標(biāo)體系構(gòu)建

    聯(lián)合國OECD和UNEP提出的“P-S-R”概念框架模型,反映人類活動對生態(tài)環(huán)境造成的壓力、生態(tài)系統(tǒng)狀態(tài)變化以及人類對于生態(tài)系統(tǒng)狀態(tài)變化的響應(yīng)三者之間的關(guān)系,目前被應(yīng)用于各個領(lǐng)域及區(qū)域的生態(tài)健康評價中[12-13]。

    在遵循科學(xué)性、整體性、代表性與敏感性、規(guī)范化、簡明性與可操作性、動態(tài)與穩(wěn)定性以及定性與定量相結(jié)合[14]等指標(biāo)體系構(gòu)建原則的基礎(chǔ)上,綜合對研究區(qū)生態(tài)環(huán)境特征和生態(tài)健康影響因素的分析,基于PSR模型,構(gòu)建黃河河口區(qū)生態(tài)系統(tǒng)健康評價指標(biāo)體系(表1)。

    表1 黃河河口區(qū)生態(tài)系統(tǒng)健康評價指標(biāo)體系

    +:該指標(biāo)為正向指標(biāo),值越大,生態(tài)系統(tǒng)越健康;-:該指標(biāo)為負(fù)向指標(biāo),值越小,生態(tài)系統(tǒng)越健康

    2.2 評價指標(biāo)計(jì)算

    利用ENVI軟件,分別選取OLI 5 6 4和TM 4 5 3原始波段進(jìn)行遙感影像合成,在幾何精校正的基礎(chǔ)上,進(jìn)行增強(qiáng)處理、噪音消除等。根據(jù)計(jì)算評價指標(biāo)所需的數(shù)據(jù),采用非監(jiān)督分類和人工目視解譯相結(jié)合,將研究區(qū)陸地景觀分成8種類型,即鹽田、養(yǎng)殖水面、林草地、河流、人工水域(水庫、坑塘、溝渠)、灘涂、人類干擾用地(耕地、建設(shè)用地、道路)和未利用地,生成景觀分類圖,并計(jì)算各類景觀面積。利用抽樣統(tǒng)計(jì),以現(xiàn)狀圖、野外調(diào)查采樣點(diǎn)對比等方法對解譯結(jié)果進(jìn)行檢驗(yàn),分類精度均達(dá)到90%以上。

    利用遙感影像紅色波段和近紅外波段構(gòu)成的歸一化植被指數(shù)進(jìn)行反演,獲取陸地區(qū)域NDVI指數(shù),并計(jì)算平均值。平均彈性度指數(shù)根據(jù)不同景觀面積和景觀彈性度分值計(jì)算獲得,參考相關(guān)文獻(xiàn)[15-17]并結(jié)合實(shí)際情況,研究區(qū)各景觀的彈性度分值見表2。景觀破碎度指數(shù)和景觀多樣性指數(shù)利用Fragstats 3.2景觀指數(shù)軟件計(jì)算獲得。

    利用遙感反演數(shù)據(jù)、實(shí)測數(shù)據(jù)和統(tǒng)計(jì)數(shù)據(jù),根據(jù)數(shù)學(xué)公式計(jì)算,獲取各評價指標(biāo),詳細(xì)計(jì)算方法見文獻(xiàn)[14]。

    表2 不同景觀類型的生態(tài)彈性分值

    2.3 數(shù)據(jù)標(biāo)準(zhǔn)化

    評價指標(biāo)類型復(fù)雜,數(shù)據(jù)性質(zhì)、量綱不同,各指標(biāo)之間沒有可比性,為了消除這種影響,通過標(biāo)準(zhǔn)化處理成無量綱數(shù)據(jù),統(tǒng)一在0—1之間[18],合理確定基準(zhǔn)值是數(shù)據(jù)標(biāo)準(zhǔn)化的前提。健康是一種相對概念,絕對健康的生態(tài)系統(tǒng)是不存在的,絕對的健康標(biāo)準(zhǔn)也是不存在的,區(qū)域生態(tài)系統(tǒng)健康評價,應(yīng)著重于探討區(qū)域生態(tài)系統(tǒng)健康的時間動態(tài),而非人為判定某時某地生態(tài)系統(tǒng)的健康與否,從而保障研究的客觀性[19-20]。因此,本研究中基準(zhǔn)值的選擇按照以下原則:1)若有國家標(biāo)準(zhǔn),則以國家標(biāo)準(zhǔn)為基準(zhǔn)值;2)若沒有國家標(biāo)準(zhǔn),則以1991年的數(shù)據(jù)作為基準(zhǔn);3)既無國家標(biāo)準(zhǔn),又無1991年數(shù)據(jù),則借鑒相關(guān)科研成果。

    評價指標(biāo)分為正向指標(biāo)和負(fù)向指標(biāo),兩者的標(biāo)準(zhǔn)化形式不同。

    對于正向指標(biāo):

    對于負(fù)向指標(biāo):

    2.4 評價等級的確定

    生態(tài)系統(tǒng)健康狀況的劃分沒有統(tǒng)一或公認(rèn)的標(biāo)準(zhǔn),本研究參考相關(guān)文獻(xiàn)資料[8,21],將黃河河口區(qū)生態(tài)系統(tǒng)健康狀況劃分為5個等級(表3)。

    表3 黃河河口區(qū)生態(tài)系統(tǒng)健康評價等級

    2.5 評價方法

    本研究通過專家打分法[22]并參考相關(guān)研究成果[7-8]構(gòu)造判斷矩陣,采用層次分析法確定指標(biāo)因子權(quán)重[23-24],通過加權(quán)求和[2,17]構(gòu)建黃河河口區(qū)生態(tài)健康評價綜合指數(shù)(CEI),獲得其生態(tài)系統(tǒng)健康評價結(jié)果。

    式中,CEI為生態(tài)系統(tǒng)健康評價綜合指數(shù),取值范圍為0—1,n為評價指標(biāo)個數(shù),Wi為第i指標(biāo)的權(quán)重,Ei為第i指標(biāo)的標(biāo)準(zhǔn)化值。各評價因子的權(quán)重與各評價因子的無量綱化標(biāo)準(zhǔn)值進(jìn)行加權(quán)計(jì)算得到上一級的健康指數(shù),依次計(jì)算,分別獲得黃河河口區(qū)生態(tài)系統(tǒng)壓力層、狀態(tài)層和響應(yīng)層健康指數(shù),最終獲得整個研究區(qū)的生態(tài)系統(tǒng)健康評價綜合指數(shù)。

    “有抗節(jié)不屈、臨危致命”作“前明諸臣”的定語,這句可譯為有為保全志節(jié)而寧死不屈、臨危獻(xiàn)身的以前明朝的臣子。

    3 結(jié)果與分析

    3.1 壓力分析

    壓力指數(shù)描述人類活動和自然過程對生態(tài)系統(tǒng)形成的干擾和脅迫,一定程度上反映生態(tài)系統(tǒng)所承受的壓力以及人類向生態(tài)系統(tǒng)索取資源的程度。根據(jù)專家打分和層次分析結(jié)果,各指標(biāo)的權(quán)重值見表4。在眾多壓力指標(biāo)中,黃河入海年徑流量對生態(tài)系統(tǒng)壓力狀況影響最大,其次為黃河斷流天數(shù)。黃河入海徑流每年向渤海輸入巨量淡水,泥沙和各種營養(yǎng)鹽類,并在河口和鄰近海域形成了適宜于海洋生物生長、發(fā)育的高生產(chǎn)力海洋生態(tài)環(huán)境,因此黃河入海徑流量和黃河斷流天數(shù)對研究區(qū)生態(tài)系統(tǒng)健康意義重大。另外,人類干擾度指數(shù)、海洋捕撈產(chǎn)量等人為影響因素對河口區(qū)域的生態(tài)系統(tǒng)壓力也影響顯著。

    從現(xiàn)狀值和基準(zhǔn)值的比較來看,人工鹽田指數(shù)、淺海養(yǎng)殖產(chǎn)量、養(yǎng)殖用地指數(shù)以及海洋捕撈產(chǎn)量變化較大,基準(zhǔn)值均為1991年數(shù)據(jù),這表明從1991年到2013年,隨著海岸帶開發(fā)及海洋過度捕撈,人工鹽田面積、淺海養(yǎng)殖產(chǎn)量、養(yǎng)殖用地面積以及海洋捕撈產(chǎn)量都大幅上升。研究區(qū)的農(nóng)藥使用強(qiáng)度遠(yuǎn)遠(yuǎn)超過規(guī)定標(biāo)準(zhǔn),而生活污水處理率、黃河入海年徑流量和黃河斷流天數(shù)都有所改善。

    根據(jù)各評價因子標(biāo)準(zhǔn)值及其權(quán)重,獲得黃河河口區(qū)壓力指數(shù)為0.8288,處于健康狀態(tài)。

    表4 黃河河口區(qū)壓力評價的各指標(biāo)標(biāo)準(zhǔn)化值及權(quán)重

    3.2 狀態(tài)分析

    狀態(tài)指數(shù)是生態(tài)系統(tǒng)健康評價中最重要的指標(biāo),不僅反映生態(tài)系統(tǒng)在自然和人為因素作用下的結(jié)果,也反映生態(tài)系統(tǒng)結(jié)構(gòu)和功能在外力干擾和脅迫下所處的狀態(tài)。黃河河口區(qū)狀態(tài)評價各指標(biāo)的權(quán)重及標(biāo)準(zhǔn)化值見表5。從狀態(tài)指標(biāo)的組合權(quán)重來看,魚類生物多樣性指數(shù)和魚類生物量的重要性最大,其次為濕地面積指數(shù)和浮游動物種類數(shù),而河口及鄰近海域的石油類、無機(jī)氮相對較小。魚類作為海洋生態(tài)系統(tǒng)食物產(chǎn)出的主要類群,其生物多樣性指數(shù)和生物量代表了海域生態(tài)系統(tǒng)的變化。濕地面積指數(shù)代表了陸地區(qū)域土地利用結(jié)構(gòu)和功能,濕地作為與海洋、森林并稱的全球三大生態(tài)系統(tǒng)之一,在涵養(yǎng)水源、調(diào)節(jié)氣候以及美化景觀等多個方面具有不可替代的作用,對維護(hù)區(qū)域生態(tài)系統(tǒng)健康至關(guān)重要。

    對各指標(biāo)現(xiàn)狀值和基準(zhǔn)值對比發(fā)現(xiàn),正向指標(biāo)中的陸地區(qū)域NDVI指數(shù)、魚類生物多樣性指數(shù)和底棲生物種類數(shù)較1991年均大幅提高。NDVI指數(shù)表征陸地區(qū)域的生產(chǎn)力水平;生物多樣性則描述了生態(tài)系統(tǒng)中物種的豐富程度,體現(xiàn)生態(tài)系統(tǒng)的復(fù)雜性,對維護(hù)生態(tài)系統(tǒng)的穩(wěn)定具有重要作用,這表明該類指標(biāo)對維護(hù)生態(tài)系統(tǒng)的健康起到積極作用。濕地面積指數(shù)、天然濕地面積指數(shù)、林草覆蓋率、魚類生物量較1991年均顯著降低,根據(jù)1991年和2013年遙感數(shù)據(jù)解譯,1991年研究區(qū)濕地面積為2679 km2,而2013年降為2155 km2;天然濕地面積由1991年的2496 km2下降為2013年的1222 km2;林草地由1991年的1240 km2下降為2013年的499 km2;2013年的魚類生物量僅為1991的3.84%。物理化學(xué)指標(biāo)中的重金屬指數(shù)遠(yuǎn)遠(yuǎn)超過《海水水質(zhì)標(biāo)準(zhǔn)GB 3097- 1997》規(guī)定的海洋漁業(yè)水域標(biāo)準(zhǔn)(≤7)。

    經(jīng)計(jì)算,黃河口區(qū)域狀態(tài)指數(shù)為0.6458,處在健康和亞健康的臨界狀態(tài)。

    表5 黃河河口區(qū)狀態(tài)評價的各指標(biāo)標(biāo)準(zhǔn)化值及權(quán)重

    *自然濕地退化指數(shù)、海岸侵蝕,將基準(zhǔn)值看作1,則現(xiàn)狀值等于變化數(shù)量與1991年數(shù)據(jù)的比值加1

    3.3 響應(yīng)分析

    響應(yīng)指數(shù)主要表征人類為防止生態(tài)環(huán)境惡化和恢復(fù)生態(tài)系統(tǒng)健康所做出的反應(yīng)。人是生態(tài)系統(tǒng)的核心,所以人類健康指標(biāo)所占的權(quán)重較大,區(qū)域的人口健康狀況對其生態(tài)系統(tǒng)響應(yīng)評價有重要影響。政策法規(guī)貫徹力度對區(qū)域生態(tài)系統(tǒng)響應(yīng)評價影響次之,這是由于政策法規(guī)的傾斜程度和執(zhí)行力度對區(qū)域社會、經(jīng)濟(jì)和生態(tài)發(fā)展都具有深遠(yuǎn)影響。

    從指標(biāo)現(xiàn)狀值和基準(zhǔn)值的比較可以看出,與1991年相比,2013年人均GDP為1991年的133倍,GDP增長率為1991年的3.5倍,區(qū)域的海洋經(jīng)濟(jì)比重也大幅增長。但人口死亡率上升,人口健康狀況較1991年有所下降。

    經(jīng)計(jì)算,黃河口區(qū)域響應(yīng)指數(shù)為0.9055,為健康狀態(tài)。

    表6 黃河河口區(qū)響應(yīng)評價的各指標(biāo)標(biāo)準(zhǔn)化值及權(quán)重

    海洋經(jīng)濟(jì)比重的基準(zhǔn)值為2006年數(shù)據(jù);國民受教育程度的基準(zhǔn)值為理想值;生態(tài)恢復(fù)工程,基準(zhǔn)值看作1,現(xiàn)狀值在比值的基礎(chǔ)上加上1

    3.4 黃河河口區(qū)健康狀況評價

    根據(jù)區(qū)域特征分析和相關(guān)參考文獻(xiàn),壓力、狀態(tài)、響應(yīng)的權(quán)重分別賦予0.297、0.539、0.164。基于上述分析和權(quán)重值,利用健康評價綜合指數(shù)模型,獲得黃河河口區(qū)綜合健康指數(shù)為0.7427。表明就整體而言,與1991年相比,2013年黃河河口區(qū)生態(tài)系統(tǒng)仍處于健康狀態(tài)(0.6—0.8),但健康狀況卻明顯下降。

    4 討論與結(jié)論

    4.1 關(guān)于研究區(qū)域的選擇

    河口生態(tài)系統(tǒng)健康不僅受自身因素的影響,來自流域、海岸和海洋的多重壓力共同影響著其健康狀況[5- 6],本著評價指標(biāo)全面、評價結(jié)果客觀的原則,該研究著眼于大空間尺度、多生態(tài)類型、綜合性評價河口生態(tài)系統(tǒng)健康狀況。理論上講,整個流域的自然壓力和人為壓力都會影響河口生態(tài)系統(tǒng)健康狀況,但由于數(shù)據(jù)獲得等方面的原因,將整個流域的指標(biāo)納入評價體系是困難的。而與河口緊鄰的三角洲濕地處于咸淡水交匯處,生態(tài)環(huán)境異常脆弱,受外界因素影響極易發(fā)生變化,且其變化對黃河河口的生態(tài)環(huán)境產(chǎn)生重要影響;此外,下游的各種生態(tài)指標(biāo)也是對流域生態(tài)環(huán)境的一定反映。因此本論文將研究范圍定義為對河口生態(tài)環(huán)境影響最大的黃河三角洲濕地、黃河河口及其一定范圍內(nèi)的相鄰海域作為一個整體。同時,評價指標(biāo)中的統(tǒng)計(jì)指標(biāo)都以一定的區(qū)劃單位獲取數(shù)據(jù),因此,黃河三角洲濕地具體指臨海的東營市河口區(qū)、東營區(qū)和墾利縣的行政區(qū)域;考慮海洋相關(guān)數(shù)據(jù)的獲得性,相鄰海域根據(jù)海洋生物資源、海洋環(huán)境海上調(diào)查范圍及其與河口的距離大致確定。

    4.2 評價指標(biāo)體系

    在河口生態(tài)系統(tǒng)健康評價中,選擇合適的評價指標(biāo)是關(guān)鍵。研究區(qū)包含了多個子系統(tǒng),并且各子系統(tǒng)間存在著交互作用,為了能夠客觀評價生態(tài)系統(tǒng)健康狀況,必須按照一套嚴(yán)格的指標(biāo)篩選原則,構(gòu)建一個能覆蓋和衡量河口區(qū)生態(tài)系統(tǒng)各個方面的指標(biāo)體系,全面反映該特定生態(tài)系統(tǒng)實(shí)際狀態(tài)及生態(tài)系統(tǒng)內(nèi)部的聯(lián)系[25]。本研究借助于PSR因果關(guān)系模型,全面分析黃河河口區(qū)生態(tài)系統(tǒng)的“壓力”、“狀態(tài)”和“響應(yīng)”,充分理解影響河口區(qū)生態(tài)系統(tǒng)中各因素的作用過程以及彼此之間的因果關(guān)系,從多個方面遴選了43個指標(biāo)系統(tǒng)地構(gòu)建了黃河河口區(qū)生態(tài)系統(tǒng)健康評價指標(biāo)體系。

    壓力指標(biāo)包含了人口壓力、陸源污染、海岸帶開發(fā)和自然壓力4個方面,為人類和自然對生態(tài)系統(tǒng)的外源性輸入和干擾,反映了一定時期內(nèi)的資源利用強(qiáng)度及其變化趨勢。狀態(tài)指標(biāo)是生態(tài)系統(tǒng)對人類壓力和自然壓力響應(yīng)后的現(xiàn)狀,狀態(tài)指標(biāo)可作為生態(tài)系統(tǒng)異常的早期預(yù)警提示,狀態(tài)指標(biāo)包含了河口區(qū)生態(tài)系統(tǒng)的生產(chǎn)力、生物多樣性、結(jié)構(gòu)和功能、自然災(zāi)害、物理化學(xué)等指標(biāo)。生產(chǎn)力表明了生態(tài)系統(tǒng)的生產(chǎn)能力,是一切生態(tài)系統(tǒng)存在的物質(zhì)基礎(chǔ);生物多樣性描述了生態(tài)系統(tǒng)中物種的豐富程度,表征了生態(tài)系統(tǒng)的復(fù)雜性,對維護(hù)生態(tài)系統(tǒng)的穩(wěn)定具有重要作用;結(jié)構(gòu)和功能指標(biāo)表明生態(tài)系統(tǒng)中不同類型的組成及其服務(wù)和產(chǎn)出功能;物理化學(xué)指標(biāo)則用于衡量河口及鄰近海域生態(tài)系統(tǒng)的非生物環(huán)境狀況。響應(yīng)指標(biāo)反映了社會或個人為了改變或預(yù)防不利于人類生存和發(fā)展的生態(tài)環(huán)境而做出的改變,主要考慮社會經(jīng)濟(jì)、人類健康、文化、政策措施等因素。所有指標(biāo)涵蓋了黃河三角洲濕地、河口、相鄰海域各區(qū)域的主要生態(tài)特征和生態(tài)問題。

    該指標(biāo)體系主要基于宏觀指標(biāo)和景觀格局指標(biāo),將遙感數(shù)據(jù)、監(jiān)測數(shù)據(jù)、統(tǒng)計(jì)數(shù)據(jù)及他人研究成果相結(jié)合進(jìn)行評價指標(biāo)的計(jì)算,與其他類似研究[7- 8]相比較,本研究提出了更全面、更綜合的評價指標(biāo)體系,并將遙感技術(shù)應(yīng)用其中。

    4.3 黃河河口區(qū)生態(tài)系統(tǒng)健康現(xiàn)狀及原因分析

    評價結(jié)果表明,以1991年及相關(guān)國家標(biāo)準(zhǔn)為基準(zhǔn),2013年黃河河口區(qū)生態(tài)系統(tǒng)健康綜合指數(shù)為0.7427,處于健康狀態(tài),但與1991年相比較,健康狀況明顯下降。于洪良[9]從濕地功能整合性、濕地生態(tài)特征及社會環(huán)境3個方面分析了墾利縣濕地生態(tài)健康狀況,總體得分為0.6865;王薇等[8]的研究表明墾利縣濕地生態(tài)系統(tǒng)處于脆弱狀態(tài);安樂生等[7]以現(xiàn)代黃河三角洲為研究區(qū)域進(jìn)行濕地生態(tài)系統(tǒng)健康狀況空間分析,結(jié)果顯示研究區(qū)濕地健康條件處于一般病態(tài)和健康之間。由于評價目的、感興趣區(qū)域、研究關(guān)注點(diǎn)甚至個人學(xué)術(shù)興趣的差異,評價指標(biāo)的選擇會有所不同,評價結(jié)果往往不能進(jìn)行絕對比較[26-27],但本研究結(jié)果反映的黃河河口區(qū)生態(tài)系統(tǒng)健康狀況與上述研究基本一致。

    狀態(tài)指數(shù)和壓力指數(shù)下降是黃河河口區(qū)生態(tài)系統(tǒng)健康綜合指數(shù)下降的主要原因。就各評價指標(biāo)而言,無論是正向指標(biāo)還是負(fù)向指標(biāo),標(biāo)準(zhǔn)化值越大,表明該指標(biāo)的現(xiàn)狀值越接近基準(zhǔn)值,對生態(tài)系統(tǒng)健康狀況貢獻(xiàn)越大;反之,則會引起生態(tài)系統(tǒng)健康狀況退化。對表4和表5各指標(biāo)權(quán)重值和標(biāo)準(zhǔn)化值綜合分析,從1991到2013年,也即黃河口區(qū)域輕度開發(fā)至人類活動強(qiáng)烈干擾階段,過度捕撈、濕地不合理開發(fā)、淺海養(yǎng)殖過度及污染排放是黃河河口區(qū)生態(tài)系統(tǒng)健康下降的主要影響因素。研究區(qū)內(nèi),海洋捕撈產(chǎn)量從1991年的587 774 t增加到2013年的1043 391 t(表4),捕撈強(qiáng)度過大使得渤海魚類小型化和低齡化,魚類資源密度大幅下降[28-29],黃河口鄰近海域魚類生物量大大降低(表5)。遙感數(shù)據(jù)解譯獲得,1991—2013年間,養(yǎng)殖蝦蟹池增加了318.66 km2,鹽田由1.56 km2增加為432.42 km2,人類干擾用地亦大幅增加,由2641.25 km2增加到3104.65 km2,致使研究區(qū)內(nèi)的濕地景觀和結(jié)構(gòu)發(fā)生了巨大變化:天然濕地面積大大下降,林草覆蓋率降低,以養(yǎng)殖水域、人工鹽田為主的人工濕地急劇增加[30-31]。濕地面積減少,導(dǎo)致濕地的平均彈性度指數(shù)和水文調(diào)節(jié)指數(shù)下降(表5)。淺海養(yǎng)殖大力發(fā)展,1991—2013年,淺海養(yǎng)殖產(chǎn)量增加了232%(表4),由于淺海養(yǎng)殖主要位于水交換能力較差的淺海灘涂和內(nèi)灣水域,使得局部水域環(huán)境惡化。陸源污染物的排放及海水養(yǎng)殖過度等原因,使得海水中的重金屬、氮、磷含量嚴(yán)重超標(biāo);海洋污染、海水富營養(yǎng)化嚴(yán)重,海域赤潮頻發(fā),面積不斷擴(kuò)大(表5)。黃河入海徑流量是河口生態(tài)系統(tǒng)的主要影響因子[32],自1999年黃河水量統(tǒng)一調(diào)度以來,黃河保持不斷流,因此黃河斷流天數(shù)和徑流量指標(biāo)都優(yōu)于基準(zhǔn)值;但由于黃河流域工農(nóng)業(yè)引沙量的急劇上升且凈流量仍然偏小,致使河口來沙量呈下降趨勢(表4)。

    黃河河口區(qū)生態(tài)系統(tǒng)健康狀況已明顯下降,應(yīng)進(jìn)行科學(xué)管理和生態(tài)恢復(fù)??刂坪Q蟛稉?養(yǎng)護(hù)海洋生物資源;嚴(yán)禁天然濕地開發(fā),逐步恢復(fù)林草覆蓋,保護(hù)三角洲濕地的生態(tài)環(huán)境;合理規(guī)劃淺海養(yǎng)殖,開展生態(tài)養(yǎng)殖;控制陸源污染等一系列措施亟待實(shí)施,但河口區(qū)生態(tài)系統(tǒng)健康恢復(fù)是復(fù)雜而長期的系統(tǒng)工程,需要建立長期的觀測和評價機(jī)制,更需要社會各界的積極參與與支持。

    需要說明的是,生態(tài)系統(tǒng)健康狀況具有時效性,同一系統(tǒng)在不同的時間內(nèi)所面臨的壓力、呈現(xiàn)的狀態(tài)及其做出的響應(yīng)都有差別。本研究基于PSR框架模型構(gòu)建的評價指標(biāo)體系是以研究區(qū)現(xiàn)有的生態(tài)特征和發(fā)展水平為基礎(chǔ)的,隨著黃河河口區(qū)生態(tài)系統(tǒng)的自然演變、人類活動方式的改變、人類對河口價值期望的改變以及人類對河口科學(xué)認(rèn)識水平的提高,反映河口區(qū)生態(tài)系統(tǒng)健康水平的指標(biāo)也會發(fā)生變化,研究區(qū)的生態(tài)系統(tǒng)健康狀況也會發(fā)生相應(yīng)改變。

    [1] 李純厚, 林琳, 徐姍楠, 戴明, 黃洪輝, 杜飛雁, 劉永, 齊占會. 海灣生態(tài)系統(tǒng)健康評價方法構(gòu)建及在大亞灣的應(yīng)用. 生態(tài)學(xué)報, 2013, 33(6): 1798- 1810.

    [2] 葉屬峰, 劉星, 丁德文. 長江河口海域生態(tài)系統(tǒng)健康評價指標(biāo)體系及其初步評價. 海洋學(xué)報, 2007, 29(4): 128- 136.

    [3] Hallett C S, Valesini F J, Clarke K R, Hesp S A, Hoeksema S D. Development and validation of fish-based, multimetric indices for assessing the ecological health of Western Australian estuaries. Estuarine, Coastal and Shelf Science, 2012, 104- 105: 102- 113.

    [4] Van Niekerk L, Adams J B, Bate G C, Forbes A T, Forbes N T, Huizinga P, Lamberth S J, Mackay C F, Petersen C, Taljaard S, Weerts S P, Whitfield A K, Wooldridge T H. Country-wide assessment of estuary health: An approach for integrating pressures and ecosystem response in a data limited environment. Estuarine, Coastal and Shelf Science, 2013, 130: 239- 251.

    [5] 楊志, 趙冬至, 林元燒. 河口生態(tài)安全評價方法研究綜述. 海洋環(huán)境科學(xué), 2011, 30(2): 296- 300.

    [6] Meng W, Liu L S. On approaches of estuarine ecosystems health studies. Estuarine, Coastal and Shelf Science, 2010, 86(3): 313- 316.

    [7] 安樂生, 劉貫群, 葉思源, 趙全升, 丁喜桂, 張建偉. 黃河三角洲濱海濕地健康條件評價. 吉林大學(xué)學(xué)報: 地球科學(xué)版, 2011, 41(4): 1157- 1165.

    [8] 王薇, 陳為峰, 李其光, 王昕, 馬海燕. 黃河三角洲濕地生態(tài)系統(tǒng)健康評價指標(biāo)體系. 水資源保護(hù), 2012, 28(1): 13- 16.

    [9] 于洪良. 基于模糊AHP的濕地生態(tài)系統(tǒng)健康評價研究——以黃河三角洲高效生態(tài)經(jīng)濟(jì)區(qū)為例. 山東財經(jīng)大學(xué)學(xué)報, 2014, (3): 33- 40.

    [10] 楊建強(qiáng), 崔文林, 張洪亮, 徐子鈞. 萊州灣西部海域海洋生態(tài)系統(tǒng)健康評價的結(jié)構(gòu)功能指標(biāo)法. 海洋通報, 2003, 22(5): 58- 63.

    [11] 孫鵬飛. 萊州灣及黃河口水域漁業(yè)資源結(jié)構(gòu)特征與漁業(yè)生態(tài)系統(tǒng)健康評價的初步分析[D]. 上海: 上海海洋大學(xué), 2014.

    [12] 劉佳. 九龍江河口生態(tài)系統(tǒng)健康評價研究[D]. 廈門: 廈門大學(xué), 2008.

    [13] 周曉蔚, 王麗萍, 鄭丙輝. 長江口及毗鄰海域生態(tài)系統(tǒng)健康評價研究. 水利學(xué)報, 2011, 42(10): 1201- 1217.

    [14] 牛明香, 王俊. 基于PSR模型的黃河河口區(qū)生態(tài)系統(tǒng)健康評價指標(biāo)體系探討. 水資源保護(hù), 2016, 32(1): 57- 63.

    [15] 劉曉丹. 基于遙感圖像的濕地生態(tài)系統(tǒng)健康評價——以大沽河河口濕地為例[D]. 青島: 中國海洋大學(xué), 2006.

    [16] 高吉喜. 可持續(xù)發(fā)展基礎(chǔ)理論探索—生態(tài)承載力理論、方法與應(yīng)用. 北京: 中國環(huán)境科學(xué)出版社, 2001.

    [17] 徐明德, 李靜, 彭靜, 鈕鍵, 曹露. 基于RS和GIS的生態(tài)系統(tǒng)健康評價. 生態(tài)環(huán)境學(xué)報, 2010, 19(8): 1809- 1814.

    [18] 彭濤, 陳曉宏. 海河流域典型河口生態(tài)系統(tǒng)健康評價. 武漢大學(xué)學(xué)報: 工學(xué)版, 2009, 42(5): 631- 634, 639- 639.

    [19] 彭建, 王仰麟, 吳健生, 張玉清. 區(qū)域生態(tài)系統(tǒng)健康評價——研究方法與進(jìn)展. 生態(tài)學(xué)報, 2007, 27(11): 4877- 4885.

    [20] 王一涵, 周德民, 孫永華. RS和GIS支持的洪河地區(qū)濕地生態(tài)健康評價. 生態(tài)學(xué)報, 2011, 31(13): 3590- 3600.

    [21] Department of Water Affairs and Forestry. Water Resource Protection Policy Implementation. Resource Directed Measures for Protection of Water Resources: Estuarine Ecosystems. Version 1.0. Department of Water Affairs and Forestry, Pretoria, 1999, http://www.africanwater.org/dwaf_rdm_page.htm.

    [22] 楊光華, 包安明, 陳曦, 劉海隆, 黃瑩, 史建康, 沈利強(qiáng). 新疆博斯騰湖濕地生態(tài)質(zhì)量的定量評價. 干旱區(qū)資源與環(huán)境, 2009, 23(2): 119- 124.

    [23] 戴新, 丁希樓, 陳英杰, 王貴田, 王東偉. 基于AHP法的黃河三角洲濕地生態(tài)環(huán)境質(zhì)量評價. 資源環(huán)境與工程, 2007, 21(2): 135- 139.

    [24] 喻立, 王建力, 李昌曉, 朱嵬, 李志剛, 李健. 基于DPSIR與AHP的寧夏沙湖濕地健康評價. 西南大學(xué)學(xué)報: 自然科學(xué)版, 2014, 36(2): 124- 130.

    [25] 孔紅梅, 趙景柱, 姬蘭柱, 陸兆華, 鄧紅兵, 馬克明, 張萍. 生態(tài)系統(tǒng)健康評價方法初探. 應(yīng)用生態(tài)學(xué)報, 2002, 13(4): 486- 490.

    [26] Meyer JL. Stream health: incorporating the human dimension to advance stream ecology. Journal of the North American Benthological Society, 1997, 16(2): 439- 447.

    [27] Kingsford R T. Aerial survey of waterbirds on wetlands as a measure of river and floodplain health. Freshwater Biology, 1999, 41(2): 425- 438.

    [28] 金顯仕. 渤海主要漁業(yè)生物資源變動的研究. 中國水產(chǎn)科學(xué), 2000, 7(4): 22- 26.

    [29] 單秀娟, 金顯仕, 李忠義, 陳云龍, 戴芳群. 渤海魚類群落結(jié)構(gòu)及其主要增殖放流魚類的資源量變化. 漁業(yè)科學(xué)進(jìn)展, 2012, 33(6): 1- 9.

    [30] 陳建, 王世巖, 毛戰(zhàn)坡. 1976- 2008年黃河三角洲濕地變化的遙感監(jiān)測. 地理科學(xué)進(jìn)展, 2011, 30(5): 585- 592.

    [31] 孫曉宇, 蘇奮振, 呂婷婷, 仉天宇, 吳迪, 付敏. 黃河三角洲濕地資源時空變化分析. 資源科學(xué), 2011, 33(12): 2277- 2284.

    [32] 黃葵, 石達(dá)扎實(shí), 鄧春蕾, 杜凱. 黃河口近海生態(tài)系統(tǒng)與入海徑流的響應(yīng)關(guān)系. 人民黃河, 2012, 34(11): 66- 69.

    Assessment of the ecosystem health of the Yellow River Estuary based on the pressure-state-response model

    NIU Mingxiang1,3,4, WANG Jun1,2,3,4,*, XU Binduo5

    1YellowSeaFisheriesResearchInstitute,ChineseAcademyofFisherySciences,Qingdao266071,China2LaboratoryforMarineEcologyandEnvironmentalScience,QingdaoNationalLaboratoryforMarineScienceandTechnology,Qingdao266200,China3KeyLaboratoryofSustainableDevelopmentofMarineFisheries,MinistryofAgriculture,Qingdao266071,China4ShandongProvincialKeyLaboratoryofFisheryResourcesandEcologicalEnvironment(SFREE),Qingdao266071,China5CollegeofFisheries,OceanUniversityofChina,Qingdao266003,China

    As the main channel of matter exchange between riverine and marine environments, estuarine ecosystems have been subjected to growing environmental interference and disturbance. Anthropogenic activities and the interactions between land and sea considerably affect estuarine ecosystems and health, and they have resulted in challenges for ecosystem management. Ecosystem health assessment, a new method of ecosystem assessment, determines the current ecosystem structure and function, and provides substantial information necessary for effective ecosystem management. Therefore, for the management of the Yellow River Estuary, a pressure-state-response (PSR) method was used to assess the health status of the ecosystem. The estuary and the adjacent marine and terrestrial ecosystems were considered an entity under the broad definition of estuarine ecosystems. An ecosystem health assessment indicator system was developed, taking external pressure, ecosystem state, and response parameters into consideration to calculate indicators of regional ecosystem health. The health status of the Yellow River estuarine ecosystem was then assessed using a comprehensive evaluation index (CEI), including historical data from 1991. The results showed that the overall index of the Yellow River Estuary was 0.7427 in October 2013. In detail, the response index was 0.9055, which indicates that the Yellow River estuarine ecosystem was “very healthy,” suggesting the policy and behavior of local government, relevant departments, and individuals reduced the pressure on the regional environment to a certain extent. The pressure index was 0.8288, which indicates that the ecosystem was at a “good status,” but has been under growing pressure and its health has deteriorated. The state index was 0.6458, suggesting that the ecosystem is “unhealthy” and is under considerable pressure, requiring improvement. Overall, the current status of the Yellow River estuarine ecosystem was reasonably “healthy” as compared to that in 1991, but it has deteriorated. The state index showed the most serious decline among pressure, state, and response indexes. The main reasons for the health degradation of the Yellow River estuarine ecosystem were investigated by analyzing the index variation and its weight. The factors resulting in the deterioration included overfishing, unreasonable wetland development, over-farming in shallow seawater, and pollution. Therefore, actions should be taken to prevent further deterioration through ecological restoration and appropriate management of the Yellow River estuarine ecosystem. Any future development should be based on scientific planning for the sustainability of the ecosystem.

    ecosystem; health assessment; PSR model; CEI (comprehensive evaluation index); Yellow River Estuary

    公益性行業(yè)(農(nóng)業(yè))科研專項(xiàng)(201303050);黃海水產(chǎn)研究所級基本科研業(yè)務(wù)費(fèi)(20603022013001);國家自然基金委員會-山東省人民政府聯(lián)合資助海洋科學(xué)研究中心項(xiàng)目(U1406403)

    2015- 08- 13;

    日期:2016- 06- 14

    10.5846/stxb201508131702

    *通訊作者Corresponding author.E-mail: wangjun@ysfri.ac.cn

    牛明香, 王俊, 徐賓鐸.基于PSR的黃河河口區(qū)生態(tài)系統(tǒng)健康評價.生態(tài)學(xué)報,2017,37(3):943- 952.

    Niu M X, Wang J, Xu B D.Assessment of the ecosystem health of the Yellow River Estuary based on the pressure-state-response model.Acta Ecologica Sinica,2017,37(3):943- 952.

    猜你喜歡
    河口區(qū)河口健康狀況
    漢江河口區(qū)四明灘段四明上護(hù)岸坍岸修復(fù)技術(shù)
    “傳承—信仰—未來”教育主題的設(shè)計(jì)者
    ——記山東省東營市河口區(qū)河安小學(xué)校長劉波
    教書育人(2017年23期)2017-09-26 06:47:50
    他們?yōu)槭裁催x擇河口
    河口,我們的家
    特殊的河口水
    河口
    閩江下游及河口區(qū)表層沉積物多種微量元素特征
    閩江下游及河口區(qū)表層沉積物粘土礦物含量分布特征及其環(huán)境意義
    高校教師健康狀況調(diào)查
    “中小學(xué)生身體健康狀況下降”問題雜談
    體育師友(2013年6期)2013-03-11 18:52:21
    十八禁高潮呻吟视频| 夜夜骑夜夜射夜夜干| 一区在线观看完整版| 亚洲熟女精品中文字幕| 亚洲欧美成人综合另类久久久| 欧美日韩一级在线毛片| 后天国语完整版免费观看| 日韩人妻精品一区2区三区| 亚洲五月婷婷丁香| 精品国产一区二区久久| 爱豆传媒免费全集在线观看| √禁漫天堂资源中文www| 婷婷丁香在线五月| 色综合欧美亚洲国产小说| 国产成人精品无人区| 99re6热这里在线精品视频| 亚洲国产日韩一区二区| videosex国产| 成人免费观看视频高清| 人体艺术视频欧美日本| 亚洲七黄色美女视频| a级毛片在线看网站| 后天国语完整版免费观看| 99久久精品国产亚洲精品| 国产精品欧美亚洲77777| 一区二区av电影网| 久久人妻熟女aⅴ| 人人妻人人爽人人添夜夜欢视频| 久久av网站| 成年人免费黄色播放视频| 久久女婷五月综合色啪小说| 天天躁夜夜躁狠狠久久av| 国产精品 国内视频| 99国产精品免费福利视频| 1024香蕉在线观看| 国产又爽黄色视频| 国产野战对白在线观看| 亚洲欧美激情在线| 色婷婷久久久亚洲欧美| 久久久久久久久免费视频了| 99久久综合免费| 丝袜脚勾引网站| 亚洲精品第二区| 秋霞在线观看毛片| 91精品三级在线观看| 免费看十八禁软件| 欧美日韩一级在线毛片| 欧美成狂野欧美在线观看| cao死你这个sao货| 亚洲国产欧美一区二区综合| 观看av在线不卡| 日韩中文字幕视频在线看片| 在线观看国产h片| 国产精品熟女久久久久浪| 精品人妻1区二区| 十八禁网站网址无遮挡| 在线观看一区二区三区激情| 18禁观看日本| 午夜两性在线视频| 丝袜在线中文字幕| 少妇精品久久久久久久| 精品人妻在线不人妻| 在线天堂中文资源库| 只有这里有精品99| 午夜免费鲁丝| 久久青草综合色| 狂野欧美激情性xxxx| 久久这里只有精品19| av天堂久久9| 国产黄色视频一区二区在线观看| 另类亚洲欧美激情| 亚洲色图 男人天堂 中文字幕| 久久影院123| 侵犯人妻中文字幕一二三四区| 日韩欧美一区视频在线观看| 人人妻,人人澡人人爽秒播 | 亚洲伊人色综图| 色94色欧美一区二区| 亚洲欧洲国产日韩| 欧美 亚洲 国产 日韩一| 在线观看人妻少妇| 视频在线观看一区二区三区| 中文字幕最新亚洲高清| 国产成人a∨麻豆精品| www.熟女人妻精品国产| 久久ye,这里只有精品| 亚洲国产欧美网| 国产成人欧美| 波野结衣二区三区在线| 久久久久久亚洲精品国产蜜桃av| 亚洲欧美精品综合一区二区三区| 精品人妻熟女毛片av久久网站| 后天国语完整版免费观看| 青春草亚洲视频在线观看| 久久久国产一区二区| 99香蕉大伊视频| 婷婷色麻豆天堂久久| 午夜久久久在线观看| 国产精品一区二区精品视频观看| 99国产精品免费福利视频| 国产熟女午夜一区二区三区| 91九色精品人成在线观看| 国产一区有黄有色的免费视频| 久久99一区二区三区| 国产麻豆69| 久久99精品国语久久久| 成人手机av| 99精国产麻豆久久婷婷| 在线观看免费午夜福利视频| 啦啦啦视频在线资源免费观看| 好男人电影高清在线观看| 成在线人永久免费视频| 欧美日韩亚洲国产一区二区在线观看 | 亚洲精品国产av蜜桃| 一区二区三区激情视频| 亚洲欧美激情在线| 久久毛片免费看一区二区三区| 国产黄频视频在线观看| 日日夜夜操网爽| 韩国精品一区二区三区| 老鸭窝网址在线观看| 国产日韩一区二区三区精品不卡| 在线亚洲精品国产二区图片欧美| 成年美女黄网站色视频大全免费| 日韩av免费高清视频| 日韩一区二区三区影片| 女人久久www免费人成看片| 免费不卡黄色视频| 亚洲专区国产一区二区| 黄色视频不卡| xxx大片免费视频| 国产精品免费视频内射| 午夜福利免费观看在线| 欧美日韩亚洲国产一区二区在线观看 | 国产日韩欧美亚洲二区| 亚洲精品自拍成人| 精品少妇内射三级| 欧美变态另类bdsm刘玥| 国产一区二区三区av在线| 少妇猛男粗大的猛烈进出视频| 涩涩av久久男人的天堂| 亚洲精品久久久久久婷婷小说| 久久精品久久精品一区二区三区| 亚洲欧美色中文字幕在线| 亚洲七黄色美女视频| 国产人伦9x9x在线观看| av线在线观看网站| 中国美女看黄片| 脱女人内裤的视频| 成年动漫av网址| 只有这里有精品99| 91精品伊人久久大香线蕉| 啦啦啦啦在线视频资源| 赤兔流量卡办理| 日韩人妻精品一区2区三区| 精品熟女少妇八av免费久了| 亚洲国产中文字幕在线视频| 欧美亚洲日本最大视频资源| 精品国产一区二区三区久久久樱花| 欧美日韩亚洲国产一区二区在线观看 | 99久久人妻综合| 亚洲欧美一区二区三区黑人| 久久热在线av| 搡老岳熟女国产| 亚洲一区二区三区欧美精品| 午夜日韩欧美国产| 久久天堂一区二区三区四区| 最新在线观看一区二区三区 | 成在线人永久免费视频| 久久影院123| 亚洲图色成人| 欧美亚洲 丝袜 人妻 在线| 午夜免费鲁丝| 一区福利在线观看| 国产麻豆69| 精品福利观看| 九草在线视频观看| 成人三级做爰电影| 亚洲专区中文字幕在线| 免费看十八禁软件| 精品少妇内射三级| 免费在线观看完整版高清| 女警被强在线播放| 手机成人av网站| 在线精品无人区一区二区三| 一级毛片 在线播放| 日韩 亚洲 欧美在线| 中文字幕人妻丝袜一区二区| 国产激情久久老熟女| 丝袜在线中文字幕| a级毛片在线看网站| 日日摸夜夜添夜夜爱| 久久国产精品男人的天堂亚洲| 汤姆久久久久久久影院中文字幕| 一区二区三区乱码不卡18| 国产欧美日韩一区二区三 | 黑人欧美特级aaaaaa片| 亚洲精品乱久久久久久| 国产精品免费视频内射| 精品国产一区二区三区久久久樱花| 一级毛片电影观看| 极品少妇高潮喷水抽搐| 后天国语完整版免费观看| 我要看黄色一级片免费的| a级片在线免费高清观看视频| 国产成人一区二区三区免费视频网站 | 亚洲av成人不卡在线观看播放网 | 美女扒开内裤让男人捅视频| 午夜免费成人在线视频| 国产一区二区三区av在线| 最近手机中文字幕大全| 捣出白浆h1v1| 欧美精品av麻豆av| 日韩视频在线欧美| 精品国产一区二区久久| 亚洲av综合色区一区| 天堂中文最新版在线下载| 我的亚洲天堂| 国产精品麻豆人妻色哟哟久久| 成年人黄色毛片网站| 中文字幕精品免费在线观看视频| 欧美日本中文国产一区发布| 欧美大码av| 波多野结衣av一区二区av| 香蕉国产在线看| 不卡av一区二区三区| 亚洲久久久国产精品| 精品久久蜜臀av无| 国产在线一区二区三区精| 老鸭窝网址在线观看| 亚洲成人手机| 另类亚洲欧美激情| 操美女的视频在线观看| 中文字幕制服av| 欧美精品人与动牲交sv欧美| 色播在线永久视频| 国产一区二区三区综合在线观看| 国产片内射在线| 国产野战对白在线观看| 国产国语露脸激情在线看| 老司机亚洲免费影院| 久久鲁丝午夜福利片| 午夜久久久在线观看| 99久久人妻综合| 只有这里有精品99| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲欧美一区二区三区黑人| 久久精品国产综合久久久| 国产成人一区二区在线| 欧美国产精品一级二级三级| 精品一区二区三卡| 七月丁香在线播放| 久久久久久人人人人人| 欧美成狂野欧美在线观看| 美女国产高潮福利片在线看| 国产有黄有色有爽视频| 日本av手机在线免费观看| 精品国产一区二区三区四区第35| 一区二区av电影网| 中文字幕亚洲精品专区| 久久99精品国语久久久| 久久久久国产精品人妻一区二区| 狠狠婷婷综合久久久久久88av| 色视频在线一区二区三区| 国产在线观看jvid| 日韩av在线免费看完整版不卡| 久久久国产精品麻豆| 国产精品九九99| 七月丁香在线播放| av福利片在线| 欧美精品啪啪一区二区三区 | 久久国产亚洲av麻豆专区| 精品人妻一区二区三区麻豆| 大型av网站在线播放| 18禁裸乳无遮挡动漫免费视频| 久久人人爽人人片av| 国产三级黄色录像| 波多野结衣一区麻豆| 久久国产精品人妻蜜桃| 亚洲图色成人| 人妻一区二区av| 电影成人av| 亚洲,一卡二卡三卡| 国产又爽黄色视频| 超色免费av| 亚洲中文日韩欧美视频| 精品亚洲成a人片在线观看| 久久久国产精品麻豆| 久久久久久亚洲精品国产蜜桃av| 精品福利观看| 中文字幕人妻丝袜一区二区| 国产精品偷伦视频观看了| 丝袜美足系列| 天天躁狠狠躁夜夜躁狠狠躁| 母亲3免费完整高清在线观看| 国产一区二区 视频在线| 色播在线永久视频| 无遮挡黄片免费观看| 精品福利永久在线观看| av一本久久久久| 欧美日韩黄片免| 亚洲人成77777在线视频| 一本色道久久久久久精品综合| 不卡av一区二区三区| 久久天堂一区二区三区四区| 欧美人与性动交α欧美精品济南到| 丝袜脚勾引网站| 免费观看人在逋| 91精品国产国语对白视频| 亚洲,欧美精品.| 一边亲一边摸免费视频| 丝袜喷水一区| 狂野欧美激情性xxxx| 午夜免费成人在线视频| 少妇人妻久久综合中文| 十八禁高潮呻吟视频| 日本一区二区免费在线视频| 飞空精品影院首页| 日本午夜av视频| 久久人人爽人人片av| 交换朋友夫妻互换小说| 老司机在亚洲福利影院| 天天操日日干夜夜撸| 国产在视频线精品| 免费观看av网站的网址| 高清欧美精品videossex| 久久国产精品大桥未久av| 国产免费视频播放在线视频| 日本黄色日本黄色录像| 一级毛片女人18水好多 | 久久 成人 亚洲| 久久精品国产亚洲av涩爱| 亚洲欧洲精品一区二区精品久久久| 欧美97在线视频| 18禁裸乳无遮挡动漫免费视频| 午夜免费成人在线视频| 一级毛片 在线播放| 亚洲av在线观看美女高潮| 女警被强在线播放| 久久久久久久大尺度免费视频| 免费高清在线观看视频在线观看| 黄网站色视频无遮挡免费观看| 1024视频免费在线观看| 亚洲五月色婷婷综合| 永久免费av网站大全| 免费日韩欧美在线观看| 欧美精品一区二区大全| 亚洲精品一卡2卡三卡4卡5卡 | 如日韩欧美国产精品一区二区三区| 九色亚洲精品在线播放| 老鸭窝网址在线观看| 亚洲一区中文字幕在线| 99国产精品一区二区蜜桃av | 亚洲熟女精品中文字幕| 亚洲精品在线美女| 亚洲色图 男人天堂 中文字幕| 热re99久久精品国产66热6| 亚洲精品日韩在线中文字幕| 建设人人有责人人尽责人人享有的| av网站免费在线观看视频| a级毛片黄视频| 高清黄色对白视频在线免费看| av不卡在线播放| 天堂中文最新版在线下载| 亚洲欧美清纯卡通| 精品免费久久久久久久清纯 | 电影成人av| 欧美日韩精品网址| 最近最新中文字幕大全免费视频 | 男女边摸边吃奶| 老司机在亚洲福利影院| 又大又爽又粗| 亚洲成人免费电影在线观看 | 日韩制服丝袜自拍偷拍| 两人在一起打扑克的视频| 日韩制服丝袜自拍偷拍| 亚洲免费av在线视频| av在线app专区| 99久久精品国产亚洲精品| 十八禁高潮呻吟视频| 日本vs欧美在线观看视频| 国产日韩欧美亚洲二区| 亚洲国产欧美在线一区| 精品一品国产午夜福利视频| 高清av免费在线| 久久久久久久国产电影| 亚洲av欧美aⅴ国产| 午夜福利视频精品| 人妻 亚洲 视频| 日韩av不卡免费在线播放| 啦啦啦在线免费观看视频4| 国产淫语在线视频| 亚洲欧洲国产日韩| 国产精品免费视频内射| 香蕉丝袜av| 黄色怎么调成土黄色| 午夜福利乱码中文字幕| 国产一区二区在线观看av| 韩国高清视频一区二区三区| 男女免费视频国产| 日日摸夜夜添夜夜爱| 午夜免费男女啪啪视频观看| 欧美黄色淫秽网站| 欧美亚洲日本最大视频资源| 交换朋友夫妻互换小说| 王馨瑶露胸无遮挡在线观看| 熟女av电影| 啦啦啦在线观看免费高清www| 成年av动漫网址| 99香蕉大伊视频| 成年av动漫网址| 色婷婷av一区二区三区视频| 天天影视国产精品| 国产精品三级大全| 欧美黄色片欧美黄色片| 一区二区三区四区激情视频| 久久女婷五月综合色啪小说| av不卡在线播放| 男的添女的下面高潮视频| 久久久久国产精品人妻一区二区| 午夜av观看不卡| 午夜精品国产一区二区电影| 成人国产av品久久久| 亚洲欧美一区二区三区久久| 国产免费福利视频在线观看| 亚洲专区国产一区二区| 男人操女人黄网站| 久久精品国产a三级三级三级| 亚洲av男天堂| 亚洲精品久久成人aⅴ小说| 亚洲精品国产色婷婷电影| 亚洲男人天堂网一区| 亚洲欧美日韩高清在线视频 | 免费看不卡的av| 50天的宝宝边吃奶边哭怎么回事| 国产亚洲欧美在线一区二区| 国产精品久久久久成人av| 国产免费福利视频在线观看| 国产视频一区二区在线看| 欧美精品一区二区免费开放| av线在线观看网站| 欧美日韩亚洲国产一区二区在线观看 | 一级毛片黄色毛片免费观看视频| 精品福利观看| 后天国语完整版免费观看| 91字幕亚洲| 一本色道久久久久久精品综合| 久久狼人影院| 久久人人爽av亚洲精品天堂| 日韩熟女老妇一区二区性免费视频| 国产在视频线精品| 男女边摸边吃奶| 蜜桃国产av成人99| 国产精品秋霞免费鲁丝片| 亚洲欧美色中文字幕在线| 丰满人妻熟妇乱又伦精品不卡| 亚洲久久久国产精品| 大片免费播放器 马上看| 国产男女超爽视频在线观看| 成年人午夜在线观看视频| 亚洲激情五月婷婷啪啪| 婷婷成人精品国产| 菩萨蛮人人尽说江南好唐韦庄| 国产又爽黄色视频| 黑人欧美特级aaaaaa片| 国产高清不卡午夜福利| 精品福利观看| 国产熟女午夜一区二区三区| 国产片特级美女逼逼视频| 亚洲一码二码三码区别大吗| 精品国产超薄肉色丝袜足j| 又大又黄又爽视频免费| 国产成人系列免费观看| 成人18禁高潮啪啪吃奶动态图| 捣出白浆h1v1| 我要看黄色一级片免费的| 色精品久久人妻99蜜桃| 纯流量卡能插随身wifi吗| 最新在线观看一区二区三区 | av片东京热男人的天堂| 80岁老熟妇乱子伦牲交| 日日爽夜夜爽网站| 亚洲av日韩在线播放| 国产亚洲av高清不卡| 欧美日韩福利视频一区二区| 欧美av亚洲av综合av国产av| 成人18禁高潮啪啪吃奶动态图| 女性生殖器流出的白浆| 亚洲自偷自拍图片 自拍| 婷婷色麻豆天堂久久| 18禁观看日本| 一边亲一边摸免费视频| 捣出白浆h1v1| 欧美国产精品va在线观看不卡| 丝袜人妻中文字幕| 成年av动漫网址| 又粗又硬又长又爽又黄的视频| 欧美黄色片欧美黄色片| 美国免费a级毛片| 亚洲国产精品一区二区三区在线| 免费看av在线观看网站| 99热网站在线观看| 久久天堂一区二区三区四区| 亚洲国产欧美网| 国产女主播在线喷水免费视频网站| 久久久精品区二区三区| 七月丁香在线播放| 亚洲成人国产一区在线观看 | 美女大奶头黄色视频| 蜜桃国产av成人99| 国产一区亚洲一区在线观看| 亚洲男人天堂网一区| 一级毛片我不卡| 亚洲精品日韩在线中文字幕| 老熟女久久久| 欧美另类一区| 人人妻人人澡人人爽人人夜夜| a 毛片基地| 精品熟女少妇八av免费久了| 美女中出高潮动态图| 校园人妻丝袜中文字幕| 男女之事视频高清在线观看 | 男的添女的下面高潮视频| 美女视频免费永久观看网站| 精品欧美一区二区三区在线| 一二三四社区在线视频社区8| 一区二区三区乱码不卡18| 日韩大码丰满熟妇| 亚洲av在线观看美女高潮| 免费av中文字幕在线| 国产精品.久久久| 曰老女人黄片| 欧美日韩福利视频一区二区| 日日夜夜操网爽| 亚洲av综合色区一区| 亚洲成国产人片在线观看| 国产不卡av网站在线观看| 99国产精品99久久久久| 久久久亚洲精品成人影院| 好男人电影高清在线观看| 国产精品欧美亚洲77777| 亚洲av日韩精品久久久久久密 | 亚洲精品美女久久av网站| 欧美性长视频在线观看| 男人操女人黄网站| 日韩 欧美 亚洲 中文字幕| 国产99久久九九免费精品| 天堂中文最新版在线下载| 婷婷成人精品国产| 大片电影免费在线观看免费| 国产欧美亚洲国产| 下体分泌物呈黄色| 精品国产一区二区久久| 91麻豆精品激情在线观看国产 | 国产99久久九九免费精品| 看免费成人av毛片| 日日爽夜夜爽网站| videosex国产| 精品少妇内射三级| 丝袜美腿诱惑在线| 亚洲自偷自拍图片 自拍| 国产亚洲精品久久久久5区| 夜夜骑夜夜射夜夜干| 欧美精品一区二区大全| 水蜜桃什么品种好| 免费高清在线观看视频在线观看| 三上悠亚av全集在线观看| 老司机午夜十八禁免费视频| 新久久久久国产一级毛片| 啦啦啦在线观看免费高清www| av视频免费观看在线观看| 精品第一国产精品| 国产在线视频一区二区| 欧美xxⅹ黑人| 一边摸一边抽搐一进一出视频| 国产成人精品无人区| 亚洲中文字幕日韩| 国产一区二区在线观看av| 美女中出高潮动态图| 欧美日韩视频高清一区二区三区二| 亚洲av日韩精品久久久久久密 | 国产高清视频在线播放一区 | 多毛熟女@视频| 亚洲,欧美,日韩| av国产久精品久网站免费入址| 免费女性裸体啪啪无遮挡网站| 一边摸一边抽搐一进一出视频| 亚洲av国产av综合av卡| 成人三级做爰电影| 国产又爽黄色视频| 中文字幕最新亚洲高清| 亚洲精品美女久久久久99蜜臀 | 久久久久久人人人人人| 国产精品秋霞免费鲁丝片| 日本91视频免费播放| 午夜久久久在线观看| a级片在线免费高清观看视频| 伦理电影免费视频| 国产在线视频一区二区| 精品少妇内射三级| 亚洲av日韩在线播放| 国产精品免费视频内射| 久久九九热精品免费| www.自偷自拍.com| 大片电影免费在线观看免费| 欧美激情 高清一区二区三区| 啦啦啦在线免费观看视频4| 亚洲九九香蕉| a级毛片黄视频| 欧美成狂野欧美在线观看| 国产在线一区二区三区精| 欧美日韩亚洲国产一区二区在线观看 | 欧美国产精品va在线观看不卡| 亚洲激情五月婷婷啪啪| 老汉色∧v一级毛片|