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

    露天礦區(qū)景觀生態(tài)風(fēng)險(xiǎn)空間分異

    2013-12-16 08:20:26吳健生黃秀蘭劉建政潘雅婧
    生態(tài)學(xué)報(bào) 2013年12期
    關(guān)鍵詞:礦區(qū)景觀區(qū)域

    吳健生,喬 娜,彭 建,黃秀蘭,劉建政,潘雅婧

    (1.城市人居環(huán)境科學(xué)與技術(shù)重點(diǎn)實(shí)驗(yàn)室,北京大學(xué)深圳研究生院城市規(guī)劃與設(shè)計(jì)學(xué)院,深圳 518055;2.地表過程分析與模擬教育部重點(diǎn)實(shí)驗(yàn)室,北京大學(xué)城市與環(huán)境學(xué)院,北京 100871)

    我國(guó)疆域遼闊,成礦地質(zhì)條件優(yōu)良、礦種齊全、總量豐富,95%以上的能源、80%以上的工業(yè)原料、70%以上的農(nóng)業(yè)生產(chǎn)資料都來自礦業(yè)[1]。目前,我國(guó)因礦產(chǎn)資源開發(fā)造成的挖損、塌陷、壓占等已損毀土地約21億畝[2],由此引發(fā)的環(huán)境污染、生態(tài)資源破壞以及地質(zhì)災(zāi)害等生態(tài)環(huán)境問題日益嚴(yán)峻。

    生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)作為管理生態(tài)環(huán)境的重要工具,其主要內(nèi)容是評(píng)價(jià)負(fù)面影響發(fā)生的可能性及危害程度[3-4]。當(dāng)前對(duì)生態(tài)風(fēng)險(xiǎn)壓力的描述已從單一化學(xué)因子,擴(kuò)展到多來源、多層次的風(fēng)險(xiǎn)因子及生態(tài)事件[5]。隨著研究尺度的擴(kuò)大,基于流域、城市等尺度進(jìn)行描述和評(píng)估生態(tài)風(fēng)險(xiǎn)的工作日益增多[6-9]。礦區(qū)是以礦山生產(chǎn)作業(yè)區(qū)為核心的一個(gè)獨(dú)立的人工、半人工生態(tài)系統(tǒng),其輻射范圍包括礦山職工及礦區(qū)農(nóng)民所在地,甚至包括依托礦業(yè)演替形成的鄉(xiāng)鎮(zhèn)、縣市及工業(yè)小區(qū)[2,10]。礦區(qū)獨(dú)特的生態(tài)環(huán)境、經(jīng)濟(jì)和社會(huì)發(fā)展規(guī)律促使其風(fēng)險(xiǎn)形成機(jī)理與發(fā)展方向有別于流域、城市、農(nóng)村等綜合性區(qū)域。由于存在諸多生態(tài)環(huán)境問題,礦區(qū)日益成為生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)的熱點(diǎn),但多以定性或半定量化研究為主[11-15]。

    景觀生態(tài)學(xué)強(qiáng)調(diào)空間格局及生態(tài)學(xué)過程與尺度之間的相互作用,并將人類活動(dòng)與生態(tài)系統(tǒng)的結(jié)構(gòu)和功能相整合[16]。隨著景觀生態(tài)學(xué)理論與生態(tài)風(fēng)險(xiǎn)的結(jié)合,景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)成為區(qū)域生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)的熱點(diǎn)方向,但是針對(duì)礦區(qū)景觀生態(tài)風(fēng)險(xiǎn)的研究仍較為少見。而從景觀斑塊的層面考慮,更能體現(xiàn)礦區(qū)生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)的綜合性和系統(tǒng)性[10,17]。因此,本研究以平朔礦區(qū)為例,構(gòu)建綜合指數(shù)定量評(píng)估露天礦區(qū)景觀生態(tài)風(fēng)險(xiǎn),并嘗試采用ESDA方法探析礦區(qū)景觀生態(tài)風(fēng)險(xiǎn)的空間分異特征。

    1 研究區(qū)概況

    平朔礦區(qū)位于東經(jīng) 112°17'—112°26',北緯39°24'—39°32',地處黃土高原晉陜蒙接壤的黑三角地帶,山西省北部的朔州市境內(nèi),屬于我國(guó)煤炭大規(guī)模集中開發(fā)地區(qū),以露天開采為主的作業(yè)方式對(duì)該區(qū)造成嚴(yán)重的土地破壞。研究區(qū)總面積為1764 km2,其中煤礦核心作業(yè)區(qū)面積約380 km2,約占研究區(qū)面積的20%。研究區(qū)屬溫帶半干旱大陸性季風(fēng)氣候,年均降雨量428.2—449.0 mm,而蒸發(fā)量1786.6—2598.0 mm,超過降水量的4倍,導(dǎo)致干旱與雨水侵蝕并發(fā)。研究區(qū)自然環(huán)境呈現(xiàn)風(fēng)蝕、水土流失嚴(yán)重、植被覆蓋度低等特征。

    本研究采用30 m分辨率的2010年7月的TM影像為數(shù)據(jù)源(處于植被生長(zhǎng)茂盛期,利于采用NDVI進(jìn)行植被與非植被和植被種類的區(qū)分),經(jīng)幾何精校準(zhǔn)與影像預(yù)處理,應(yīng)用決策樹算法進(jìn)行遙感解譯。參考相關(guān)專家關(guān)于露天煤礦生態(tài)風(fēng)險(xiǎn)的研究成果,將挖損與壓占等土地破壞視為煤礦區(qū)主要的人類干擾活動(dòng)[1,18],綜合全國(guó)土地利用分類標(biāo)準(zhǔn),將研究區(qū)景觀類型為林地、荒草地、坡耕地、平原耕地、裸地、采掘地、壓占地、水域、居民點(diǎn)九類?;谝延形墨I(xiàn)的研究尺度[19],及研究區(qū)的空間異質(zhì)性、斑塊大小并保證柵格信息的完整,將研究區(qū)劃分為1.5 km×1.5 km的采樣單元,共計(jì)784個(gè)(圖1)。

    同時(shí),根據(jù)人類活動(dòng)對(duì)景觀的干擾程度將礦區(qū)分為人工區(qū)、半自然區(qū)和自然區(qū)(圖1)。其中,人工區(qū)包括礦業(yè)核心區(qū)、城市生活區(qū)與礦業(yè)生活區(qū),依據(jù)明顯的道路及礦區(qū)作業(yè)邊緣等人工線性地物劃分;半自然區(qū)依據(jù)山脊線、道路、河流等自然線性與人工線性地物綜合劃分;自然區(qū)則以山脊線(自然線性地物)為劃分依據(jù)。

    圖1 研究區(qū)景觀分類及功能分區(qū)Fig.1 Landscape classification and different functional zones of opencast mining area

    2 研究方法

    2.1 景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)

    景觀生態(tài)風(fēng)險(xiǎn)指基于景觀尺度的環(huán)境污染、人為活動(dòng)及自然災(zāi)害等干擾源對(duì)生態(tài)系統(tǒng)的結(jié)構(gòu)和功能造成不利影響的可能性和危害程度[20]。評(píng)價(jià)的核心內(nèi)容分為3部分:確定干擾源對(duì)區(qū)域生態(tài)環(huán)境的作用效果,構(gòu)建景觀干擾度指數(shù);評(píng)估景觀要素抵抗外界干擾的能力,構(gòu)建景觀脆弱度指數(shù);應(yīng)用恰當(dāng)?shù)脑u(píng)估模型綜合兩者得到景觀生態(tài)風(fēng)險(xiǎn)指數(shù),評(píng)價(jià)研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)水平。

    2.1.1 景觀干擾度指數(shù)

    不同景觀類型在維護(hù)生態(tài)系統(tǒng)結(jié)構(gòu)和功能完整性的過程中會(huì)面臨外界環(huán)境的干擾[21],而景觀所受外界干擾的程度可由其景觀結(jié)構(gòu)的變化程度予以度量。參考相關(guān)文獻(xiàn)以景觀結(jié)構(gòu)指數(shù)為基礎(chǔ),通過表征景觀結(jié)構(gòu)的各指數(shù)疊加構(gòu)建景觀干擾度指數(shù)Ei[17],表達(dá)式如下:

    式中,Ci、Di、DQi分別為景觀破碎度、景觀分離度、景觀優(yōu)勢(shì)度,公式內(nèi)容及表征意義詳見文獻(xiàn)[16,22]。其中,破碎度值越大,表明景觀單元內(nèi)部穩(wěn)定性越低;分離度指數(shù)越大,表明景觀在地域上分布越分散,受到的干擾程度越大;景觀優(yōu)勢(shì)度表征景觀結(jié)構(gòu)中某一類型支配景觀的程度,反映了該景觀類型對(duì)景觀格局形成和變化影響的大小。

    對(duì)上述指標(biāo)進(jìn)行歸一化,以W1、W2、W3為對(duì)應(yīng)權(quán)重,三者相加為1。權(quán)重值不同,反映各景觀指數(shù)解釋景觀所受干擾的能力不同[17],參考相關(guān)文獻(xiàn)綜合確定破碎度指數(shù)、分離度指數(shù)和優(yōu)勢(shì)度指數(shù)的重要性差異,將三者分別賦值為 0.6、0.3、0.1 的權(quán)值[23,24]。

    2.1.2 景觀脆弱度指數(shù)

    不同的景觀類型在維護(hù)生物多樣性、保護(hù)物種與促進(jìn)景觀結(jié)構(gòu)自然演替等方面的作用是有差別的,同時(shí)抵抗外界干擾的能力也不同[17,21],由此構(gòu)建景觀脆弱度指數(shù)表征各景觀類型內(nèi)部結(jié)構(gòu)的易損性。易損性越大,抗干擾的能力越小,表明生態(tài)風(fēng)險(xiǎn)越大。綜合各類景觀易損性的特點(diǎn)[18-21],將九類景觀按脆弱程度由高到低賦值:采掘地9、壓占地8、裸地7、水域6、坡耕地5、荒草地4、平原耕地3、居民點(diǎn)2、林地1,歸一化得到脆弱度指數(shù)Fi。

    2.1.3 景觀生態(tài)風(fēng)險(xiǎn)指數(shù)

    基于景觀結(jié)構(gòu),引入景觀面積比重,建立景觀結(jié)構(gòu)指數(shù)與礦區(qū)生態(tài)風(fēng)險(xiǎn)之間的聯(lián)系,用于描述一個(gè)樣地內(nèi)整體生態(tài)風(fēng)險(xiǎn)的相對(duì)大小,通過采樣的方法將景觀空間格局轉(zhuǎn)化為空間化的生態(tài)風(fēng)險(xiǎn)變量[17]。其計(jì)算公式如下:

    式中,ERI為景觀生態(tài)風(fēng)險(xiǎn)指數(shù);n為景觀類型數(shù)量;Ei為景觀類型i的干擾度指數(shù);Fi為景觀類型i的脆弱度指數(shù);Aki為第k個(gè)采樣單元i類景觀組分的面積,Ak為第k個(gè)采樣單元的總面積。

    2.2 空間數(shù)據(jù)探索性分析(ESDA)

    ESDA是一系列空間數(shù)據(jù)分析方法和技術(shù)的集合[25],它以空間關(guān)聯(lián)測(cè)度為核心,注重?cái)?shù)據(jù)的空間關(guān)聯(lián)性、集聚性與異質(zhì)性,通過對(duì)事物空間分布格局的描述,揭示空間關(guān)聯(lián)特征與模式[25-26]。國(guó)內(nèi)外研究多集中于區(qū)域經(jīng)濟(jì)[27-29]、社會(huì)和政治研究、城市社會(huì)結(jié)構(gòu)[26,30]、空間結(jié)構(gòu)演變[31]等方面。景觀生態(tài)風(fēng)險(xiǎn)作為空間變量,探析其空間特征有助于理解風(fēng)險(xiǎn)的發(fā)生機(jī)理及變化趨勢(shì)等深層次問題。

    2.2.1 全局空間自相關(guān)

    空間自相關(guān)反映了某一變量在空間上是否相關(guān)及其相關(guān)程度,常用測(cè)度指標(biāo)為Moran's I指數(shù),計(jì)算公式詳見有關(guān)文獻(xiàn)[25-26]。Moran's I系數(shù)取值在-1至+1之間,其絕對(duì)值越趨近于1,表示研究單元的空間自相關(guān)程度越強(qiáng)。一般使用Moran's I的標(biāo)準(zhǔn)化統(tǒng)計(jì)量Z來檢驗(yàn)?;谝欢@著性水平,當(dāng)Moran's I>0,表明存在正的空間自相關(guān),研究單元屬性值呈趨同集聚;Moran's I<0,表明存在負(fù)的空間自相關(guān),研究單元屬性值呈離散分布;Moran's I=0,表明不存在空間自相關(guān),空間單元觀測(cè)值呈隨機(jī)分布[25,27]。

    2.2.2 局部空間自相關(guān)

    全局Moran's I指數(shù)值能夠測(cè)度事物在整體空間上的相關(guān)程度,但對(duì)于局部異?,F(xiàn)象,需引入局部空間自相關(guān)方法進(jìn)行探析[25]。Moran's I散點(diǎn)圖不能反映空間關(guān)聯(lián)類型的顯著性水平,而LISA(Local Indicators of Spatial Association,局域空間關(guān)聯(lián)指數(shù))分析方法可解決此類問題。其計(jì)算公式如下:

    式中,Ii為L(zhǎng)ISA指數(shù)空間單元值;Zi和Zj分別為空間單元i和j上觀測(cè)值的標(biāo)準(zhǔn)化值;Wij為空間權(quán)重矩陣。LISA的實(shí)質(zhì)在于將Moran's I分解并呈現(xiàn)到各個(gè)區(qū)域單元,并形成LISA聚類圖,由此識(shí)別局部空間高高集聚的“熱點(diǎn)”和低低集聚的“冷點(diǎn)”,并探析局部空間異常特征。

    3 結(jié)果分析

    3.1 礦區(qū)景觀生態(tài)風(fēng)險(xiǎn)

    根據(jù)自然斷裂點(diǎn)法分級(jí)原理[32-33]將景觀生態(tài)風(fēng)險(xiǎn)值分為3級(jí)(圖2):風(fēng)險(xiǎn)等級(jí)由高到低呈環(huán)形包圍趨勢(shì),核心區(qū)域?yàn)楦唢L(fēng)險(xiǎn)區(qū),景觀組分以采掘地、壓占地、居民點(diǎn)、水域?yàn)橹?中等風(fēng)險(xiǎn)區(qū)域部分集中于高風(fēng)險(xiǎn)區(qū)域外圍的平原耕地、居民點(diǎn)及部分壓占地,部分集中于人類擾動(dòng)相對(duì)較少的裸地、坡耕地區(qū)域;低風(fēng)險(xiǎn)區(qū)域基本介于兩部分中等風(fēng)險(xiǎn)區(qū)域之間,景觀組分以人類擾動(dòng)最少的林地、荒草地為主。

    由此可見,高風(fēng)險(xiǎn)集中于礦業(yè)生產(chǎn)與城市生活區(qū)域,并且由中心向外呈現(xiàn)水平遞減的趨勢(shì)。與此相應(yīng),風(fēng)險(xiǎn)等級(jí)逐漸過渡到中、低水平,擾動(dòng)過渡為城市生活、農(nóng)業(yè)生產(chǎn)為主,波及范圍較大,景觀逐漸轉(zhuǎn)為半自然、自然類型。如圖2所示,城市生活區(qū)風(fēng)險(xiǎn)最高,次之為礦業(yè)核心區(qū),半自然區(qū)域中等,最低為自然區(qū)。礦業(yè)核心區(qū)除了大部分采礦用地類型還分布有風(fēng)險(xiǎn)值較低的礦區(qū)復(fù)墾林、荒草地及坡耕地,這些景觀對(duì)整體生態(tài)環(huán)境有緩沖、調(diào)節(jié)作用,而城市生活區(qū)由于包含高風(fēng)險(xiǎn)的水域、居民點(diǎn)、裸地使得整體的風(fēng)險(xiǎn)值最高。半自然區(qū)依河流劃分為兩塊區(qū)域,右側(cè)緊鄰城市生活區(qū)的區(qū)域風(fēng)險(xiǎn)值高于左側(cè)區(qū)域。礦業(yè)生活區(qū)的風(fēng)險(xiǎn)水平位于兩者之間,但半自然區(qū)域平均生態(tài)風(fēng)險(xiǎn)值小于礦業(yè)生活區(qū)。自然區(qū)生態(tài)風(fēng)險(xiǎn)最小,主要原因是區(qū)域內(nèi)存在大面積干擾度較低、脆弱度較低的林地、荒草地等景觀。

    綜合研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)特征(表1),可以發(fā)現(xiàn),干擾度、脆弱度、景觀生態(tài)風(fēng)險(xiǎn)各等級(jí)采樣單元的比例差異明顯。其中,干擾度整體以中、低水平為主,說明研究區(qū)承受外界擾動(dòng)的劇烈程度在空間上差異較小;而采樣單元中74.36%為中等脆弱度水平,呈現(xiàn)“兩頭少,中間多”的正態(tài)分布模式,說明研究區(qū)景觀脆弱度總體水平較為均衡;由干擾度與脆弱度得到的景觀生態(tài)風(fēng)險(xiǎn)綜合了兩者的比例特征,高風(fēng)險(xiǎn)區(qū)比例最小,中等風(fēng)險(xiǎn)區(qū)比例最大,各等級(jí)風(fēng)險(xiǎn)的采樣單元比例近似等于同等級(jí)干擾度與脆弱度比例加和的平均值。

    圖2 礦區(qū)采樣單元與功能分區(qū)景觀生態(tài)風(fēng)險(xiǎn)空間分布Fig.2 landscape eco-risk of sample units and the different functional zones of opencast mining area

    表1 礦區(qū)景觀生態(tài)風(fēng)險(xiǎn)統(tǒng)計(jì)Table 1 landscape eco-risk statistics of opencast mining area

    3.2 礦區(qū)景觀生態(tài)風(fēng)險(xiǎn)空間自相關(guān)性

    礦區(qū)景觀生態(tài)風(fēng)險(xiǎn)空間自相關(guān)指數(shù)的計(jì)算結(jié)果表明,Moran's I為0.7521且隨機(jī)分布檢驗(yàn)Z值高達(dá)39.7341,說明礦區(qū)景觀生態(tài)風(fēng)險(xiǎn)值存在顯著的空間正相關(guān),即風(fēng)險(xiǎn)值高的區(qū)域,周邊區(qū)域的風(fēng)險(xiǎn)值亦高;風(fēng)險(xiǎn)值低的區(qū)域,周邊區(qū)域的風(fēng)險(xiǎn)值亦低;空間趨同集聚現(xiàn)象明顯。同時(shí),標(biāo)準(zhǔn)差僅為0.0252,相對(duì)較小,說明風(fēng)險(xiǎn)值大小波動(dòng)不明顯,總體水平較為平穩(wěn)。

    對(duì)比景觀生態(tài)風(fēng)險(xiǎn)指數(shù)值的空間鄰接特征(圖3),可以發(fā)現(xiàn),HH象限內(nèi)的研究單元風(fēng)險(xiǎn)值差距較大,散點(diǎn)較為分散,即高風(fēng)險(xiǎn)的區(qū)域內(nèi)部差異較大;LL象限內(nèi)的研究單元風(fēng)險(xiǎn)值差距較小,散點(diǎn)較密集,即低風(fēng)險(xiǎn)的區(qū)域內(nèi)部差異不明顯;LH象限的研究單元數(shù)目較少,低高離群程度較為顯著。即存在有一部分風(fēng)險(xiǎn)“低洼”區(qū)域(低風(fēng)險(xiǎn)區(qū)被高風(fēng)險(xiǎn)區(qū)包圍);HL象限的研究單元數(shù)目最少且高低離群程度更低。即存在個(gè)別風(fēng)險(xiǎn)“高聳”區(qū)域(高風(fēng)險(xiǎn)區(qū)被低風(fēng)險(xiǎn)區(qū)包圍)。同時(shí),散點(diǎn)接近于回歸線,表明在局部空間上,它們呈現(xiàn)更顯著的“同質(zhì)集聚、異質(zhì)隔離”特征。

    基于LISA指數(shù)分析礦區(qū)景觀生態(tài)風(fēng)險(xiǎn)的集聚特征(圖3),并通過P≤0.05的顯著性檢驗(yàn)。圖中“熱點(diǎn)”、“冷點(diǎn)”十分突出,同時(shí)存在一定數(shù)量的空間“奇異值”。其中,“熱點(diǎn)”主要集聚于礦業(yè)核心區(qū)與城市生活區(qū);“冷點(diǎn)”主要集聚在半自然、自然區(qū)域?!盁狳c(diǎn)”、“冷點(diǎn)”的位置與功能分區(qū)的對(duì)應(yīng)關(guān)系證明人類活動(dòng)的集聚程度與風(fēng)險(xiǎn)的空間集聚程度是對(duì)應(yīng)的。低高離群“奇異點(diǎn)”位于礦業(yè)核心區(qū),原因是該區(qū)域的復(fù)墾林地周邊存在較高風(fēng)險(xiǎn)的裸地、采掘地;而高低離群“奇異點(diǎn)”位于半自然、自然區(qū),原因是裸地區(qū)域被林地與坡耕地包圍??傮w而言,“奇異點(diǎn)”一般為局部風(fēng)險(xiǎn)發(fā)生急劇變化的區(qū)域,也是風(fēng)險(xiǎn)最可能蔓延的位置。

    3.3 景觀生態(tài)風(fēng)險(xiǎn)驅(qū)動(dòng)因子分析

    景觀生態(tài)風(fēng)險(xiǎn)的空間差異的決定因素包括:景觀的干擾度與脆弱度。景觀干擾度可視為外因(人為干擾),景觀脆弱度可視為內(nèi)因(景觀自身維持生態(tài)穩(wěn)定性的能力),而內(nèi)、外因在影響風(fēng)險(xiǎn)的大小及空間分異的作用是有差別的。

    圖3 礦區(qū)景觀生態(tài)風(fēng)險(xiǎn)Moran散點(diǎn)分布與生態(tài)風(fēng)險(xiǎn)值LISA指數(shù)空間特征Fig.3 Moran scatter plot and Cluster map of LISA of landscape eco-risk of opencast mining area

    本研究分析各風(fēng)險(xiǎn)等級(jí)與干擾度及脆弱度水平的對(duì)應(yīng)關(guān)系,理論上共有27種組合模式,但實(shí)際疊加只有14種對(duì)應(yīng)模式(圖4)。其中,低風(fēng)險(xiǎn)采樣單元比例較高,以LLM、LML模式為主,各模式比例差異較大;中等風(fēng)險(xiǎn)采樣單元比例最高,各模式比例差距較小,其中以MMM、MLM與MML模式為主,其余模式均為零星分布;高風(fēng)險(xiǎn)采樣單元以HHM、HHH兩種模式為主??偨Y(jié)而言,低風(fēng)險(xiǎn)且低干擾度的采樣單元比例最大(31.5%),說明決定低風(fēng)險(xiǎn)狀態(tài)的主導(dǎo)因素是景觀干擾度;同理,高風(fēng)險(xiǎn)狀態(tài)下,高干擾度的采樣單元比例最大(7.65%),即景觀干擾度起主導(dǎo)作用;而中等風(fēng)險(xiǎn)狀態(tài)下,中等脆弱度采樣單元比例最大(39.8%),即景觀脆弱度起主導(dǎo)作用。由此可知,低、高風(fēng)險(xiǎn)主要由“外因”驅(qū)動(dòng)為主,而中等風(fēng)險(xiǎn)以“內(nèi)因”驅(qū)動(dòng)為主。

    圖4 礦區(qū)景觀生態(tài)風(fēng)險(xiǎn)結(jié)構(gòu)統(tǒng)計(jì)Fig.4 The statistics of landscape eco-risk structure

    同時(shí),對(duì)比各風(fēng)險(xiǎn)組合模式的空間分布格局(圖5)與景觀功能分區(qū),可以發(fā)現(xiàn),由低風(fēng)險(xiǎn)到高風(fēng)險(xiǎn),各種結(jié)構(gòu)的集聚形式趨于簡(jiǎn)單。其中,HHH模式主要集聚于礦業(yè)核心區(qū),HHM主要集聚于礦業(yè)生活區(qū)外圍與城市生活區(qū)的核心;MML、MMM與MLM3種主要模式大面積分布于自然、半自然區(qū)域、城市生活區(qū)及礦業(yè)核心區(qū)邊緣,整體面積較大、分布較散;LLM、LML模式局部集聚特征明顯,且于自然區(qū)、半自然區(qū)、礦業(yè)核心區(qū)分布較多。綜合景觀生態(tài)風(fēng)險(xiǎn)的數(shù)量結(jié)構(gòu)與空間格局,總體來看,區(qū)域景觀生態(tài)風(fēng)險(xiǎn)越大,景觀類型越單一,風(fēng)險(xiǎn)的空間集聚程度越大。由此可見,人類干擾是景觀多樣性喪失、景觀結(jié)構(gòu)單一化的重要影響因素。

    4 結(jié)論與討論

    本研究基于風(fēng)險(xiǎn)評(píng)價(jià)原理,結(jié)合景觀生態(tài)學(xué)與生態(tài)風(fēng)險(xiǎn)理論,構(gòu)建景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)模型定量表征露天礦區(qū)景觀生態(tài)風(fēng)險(xiǎn)狀況。研究結(jié)果表明,礦區(qū)大部分區(qū)域的景觀生態(tài)風(fēng)險(xiǎn)處于中、低水平,礦業(yè)生產(chǎn)地區(qū)處于高水平;而大規(guī)模的土地破壞、壓占等礦業(yè)生產(chǎn)與生活擾動(dòng)是高風(fēng)險(xiǎn)的主要來源。雖然這種干擾波及范圍小,但呈現(xiàn)出由礦業(yè)核心生產(chǎn)區(qū)向外蔓延的趨勢(shì);中等風(fēng)險(xiǎn)區(qū)域以裸地、坡耕地、平原耕地等景觀類型為主,風(fēng)險(xiǎn)來源主要是城市生活與農(nóng)業(yè)生產(chǎn)。此外,惡劣的自然條件造成水土流失嚴(yán)重、植被退化,形成大面積生態(tài)穩(wěn)定性低的裸地、荒草地,使得平朔礦區(qū)的自然生態(tài)環(huán)境整體較為脆弱。因此,在生態(tài)環(huán)境治理、生態(tài)風(fēng)險(xiǎn)防范的過程中,需要從自然條件出發(fā),因地制宜,采用適宜的植被覆蓋類型、多層次的治理、復(fù)墾手段對(duì)礦區(qū)環(huán)境進(jìn)行分類型、分重點(diǎn)治理。

    圖5 礦區(qū)景觀生態(tài)風(fēng)險(xiǎn)結(jié)構(gòu)空間分布Fig.5 Landscape eco-risk structure and distribution in opencast mine area

    對(duì)比不同功能區(qū)的景觀生態(tài)風(fēng)險(xiǎn),表明礦區(qū)整體呈現(xiàn)人工區(qū)風(fēng)險(xiǎn)>半自然區(qū)>自然區(qū)的空間分異特征。但是,礦業(yè)核心區(qū)不是風(fēng)險(xiǎn)最高的區(qū)域,城市生活區(qū)風(fēng)險(xiǎn)水平反而更高,同樣礦業(yè)生活區(qū)比部分半自然區(qū)域風(fēng)險(xiǎn)值小。這些異??赡苁秋L(fēng)險(xiǎn)值在空間的累積速度小于面積的擴(kuò)張速度(面積越大,其均值就越小),且這種面積與均值的關(guān)系在決定功能分區(qū)的風(fēng)險(xiǎn)平均水平方面呈主導(dǎo)作用。此外,不同區(qū)域風(fēng)險(xiǎn)值的差異一方面表明干擾的空間分異特征較為明顯,另一方面說明人類的生產(chǎn)、生活擾動(dòng)是主要的風(fēng)險(xiǎn)源。因此,從基礎(chǔ)環(huán)節(jié)層面控制風(fēng)險(xiǎn)要從減少人類擾動(dòng)、減輕破壞程度、科學(xué)處理已有破壞痕跡等方面出發(fā)。

    進(jìn)一步的景觀生態(tài)風(fēng)險(xiǎn)空間自相關(guān)分析表明,礦區(qū)景觀生態(tài)風(fēng)險(xiǎn)值呈現(xiàn)明顯的空間正相關(guān)關(guān)系,說明風(fēng)險(xiǎn)具有顯著的空間集聚特征。局部空間的景觀生態(tài)風(fēng)險(xiǎn)值具有明顯的“冷點(diǎn)”、“熱點(diǎn)”區(qū)域,這為針對(duì)性開展生態(tài)環(huán)境治理工作提供了空間參考范圍。同時(shí),明顯的“奇異點(diǎn)”為預(yù)測(cè)區(qū)域的風(fēng)險(xiǎn)水平將要發(fā)生變化提供了依據(jù),對(duì)風(fēng)險(xiǎn)防控起到預(yù)警作用;尤其是“低洼”區(qū)域極易成為下階段風(fēng)險(xiǎn)擴(kuò)散的突破口,應(yīng)在礦區(qū)環(huán)境管理工作中予以重點(diǎn)關(guān)注。而對(duì)比風(fēng)險(xiǎn)驅(qū)動(dòng)因子的組合形式發(fā)現(xiàn),高風(fēng)險(xiǎn)與低風(fēng)險(xiǎn)區(qū)域,干擾度為主要驅(qū)動(dòng)因子,即造成風(fēng)險(xiǎn)的外因起主導(dǎo)作用;而中等風(fēng)險(xiǎn)區(qū)域,脆弱度為主要驅(qū)動(dòng)因子,即景觀內(nèi)在的特質(zhì)為決定性因素。因此,對(duì)于低風(fēng)險(xiǎn)和高風(fēng)險(xiǎn)區(qū)域,宜重點(diǎn)控制外界擾動(dòng)強(qiáng)度,從風(fēng)險(xiǎn)源頭控制;而針對(duì)中等風(fēng)險(xiǎn)區(qū)域,則需加強(qiáng)脆弱度低的景觀類型的覆蓋度。

    但是,需要指出的是,上述研究結(jié)果均為基于一期遙感影像對(duì)未來風(fēng)險(xiǎn)的估測(cè),對(duì)風(fēng)險(xiǎn)驅(qū)動(dòng)的規(guī)律有待于在多期影像的長(zhǎng)時(shí)間序列評(píng)估中進(jìn)一步驗(yàn)證。同時(shí),囿于TM遙感影像空間分辨率,本研究未能進(jìn)一步細(xì)分礦業(yè)核心生產(chǎn)區(qū)域,有待深入探析礦山開采作業(yè)區(qū)內(nèi)部不同礦業(yè)生產(chǎn)單元之間的景觀生態(tài)風(fēng)險(xiǎn)空間分異特征。此外,采用專家打分法對(duì)干擾度權(quán)重及脆弱度進(jìn)行賦值,這種方法盡管具備專家經(jīng)驗(yàn)的優(yōu)點(diǎn),仍存在一定主觀性,對(duì)研究結(jié)果帶來了一定的不確定性。

    [1] Chang Q,Qiu Y,Xie M M,Peng J.Theory and method of ecological risk assessment for mining areas based on the land destruction.Acta Ecologica Sinica,2012,32(16):5164-5174.

    [2] Li J C,Bai Z K.Land Reclamation and Ecological Rehabilitation of Land in Surface Mines.Beijing:Science Press,2000.

    [3] Suter G W II,Norton S B,Barnthouse L W.The evolution of frameworks for ecological risk assessment from the Red Book Ancestor.Human and Ecological Risk Assessment,2003,9(5):1349-1360.

    [4] Chen H,Liu J S,Cao Y,Li S C,Ouyang H.Progresses of ecological risk assessment.Acta Ecologica Sinica,2006,26(5):1558-1566.

    [5] Landis W G.Twenty years before and hence:ecological risk assessment at multiple scales with multiple stressors and multiple endpoints.Human and Ecological Risk Assessment,2003,9(5):1317-1326.

    [6] Fu Z Y,Xu X G.Regional ecological risk assessment.Advance in Earth Sciences,2001,16(2):267-271.

    [7] Fu Z Y,Xu X G,Lin H P,Wang X L.Regional ecological risk assessment of in the Liaohe river delta wetlands.Acta Ecologica Sinica,2001,21(3):365-373.

    [8] Sun H B,Yang G S,Su W Z,Zhu T M,Wang R R.Ecological risk assessment of land use in the area along Changjiang River:a case study of Nanjing,China.Acta Ecologica Sinica,2010,30(20):5616-5625.

    [9] Zeng Y.The regional ecological risk assessment of Hohhot City.Acta Ecologica Sinica,2010,30(3):668-673.

    [10] Pan Y J,Wang Y L,Peng J,Han Y N.Research progress in ecological risk assessment of mining area.Acta Ecologica Sinica,2012,32(20):6566-6574.

    [11] Bai Z K,Zhao J K,Li J C,Wang W Y,Lu C E,Ding X Q,Cai S J,Chen J J.Ecosystem damage in a large opencast coal mine——A case study on PINGSHUO Surface Coal Mine,China.Acta Ecologica Sinica,1999,19(6):870-875.

    [12] Hattemer-Frey H A,Quinlan R E,Krieger G R.Ecological risk assessment case study:impacts to aquatic receptors at a former metals mining superfund site.Risk Analysis,1995,15(2):253-265.

    [13] Michalik B.NORM impacts on the environment:an approach to complete environmental risk assessment using the example of areas contaminated due to mining activity.Applied Radiation and Isotopes,2008,66(11):1661-1665.

    [14] Jia Y,Cao L X.Review on ecological risk assessment method for coal-mining area.Environmental Science and Management,2011,36(4):177-182.

    [15] Ma X.Ecological Risk Assessment of Vulnerable Mine Area[D]:Beijing:China University of Geoscience,2011.

    [16] Wu J G.Landscape Ecology-Pattern,Process,Scale and Grade.2nd ed.Beijing:Beijing Higher Education Press,2000.

    [17] Xie H L.Regional eco risk analysis of based on landscape structure and spatial statistics.Acta Ecologica Sinica,2008,28(10):5020-5026.

    [18] Cheng J L,Lu Z H,F(xiàn)an Y H.Method of ecological risk assessment for opencast mine area.Acta Ecologica Sinica,2004,24(12):2945-2950.

    [19] Hu H B,Liu H Y,Hao J F,An J.The urbanization effects on watershed landscape structure and their ecological risk assessment.Acta Ecologica Sinica,2011,31(12):3432-3440.

    [20] Li X H,Li J Y.Analysis on regional landscape ecological risk based on GIS——A case study along the lower reaches of the Weihe river.Arid Zone Research,2008,25(6):899-903.

    [21] Cheng P,Pan X L.Ecological risk analysis of regional landscape in inland river watershed of arid area——a case study of Sangong river basin in Fukang.Chinese Journal of Ecology,2003,22(4):116-120.

    [22] Wang G X,Chen G D.Study on the landscape pattern of a desert-Oasis ecological system:a spatial grid method and its application.Arid Zone Research,1999,16(3):6-11.

    [23] Gao B,Li X Y,Li Z G,Chen W,He X Y,Qi S Z.Assessment of ecological risk of coastal economic developing zone in Jinzhou Bay based on Landscape pattern.Acta Ecologica Sinica,2011,31(12):3441-3450.

    [24] Lu Y,Su W J,Hua C,Tang C Y.Landscape ecological risk assessment for upper Zuojiang River basin.Tropical Geography,2010,30(5):496-503.

    [25] Anselin L.Local indicators of spatial association-LISA.Geographical Analysis,1995,27(2):93-115.

    [26] Xuan G F,Xu J G,Zhao J.An analysis of urban social space based on ESDA——A case study of the central urban district in Shanghai.Scientia Geographica Sinica,2010,30(1):22-29.

    [27] Lee S II.Spatial Association Measures for an ESDA-GIS Framework:Developments,Significance Tests,and Applications to Spatio-Temporal Income Dynamics of United States Labor Market Areas,1969—1999[D].Ohio:Ohio State University,2001.

    [28] Meng B,Wang J F,Zhang W Z,Liu X H.Evaluation of regional disparity in China based on spatial analysis.Scientia Geographica Snica,2005,25(4):393-400.

    [29] Zhang H F,Bai Y P,Chen Q,Wang B H,Niu D W.Regional economic disparities in Qinghai Province based on ESDA and GIS.Arid Land Geography,2009,32(3):454-461.

    [30] Lü A M,Li C M,Lin Z J,Shi W Z.Population growth rate and its spatial association by providence in China.Acta Geographica Sinica,2002,57(2):143-150.

    [31] Huang H C,Li M Y.Spatial urban structure evolution analysis of the forerunner development zone of the Changchun-Jilin-Tumen River development region based on the ESDA-Take Yanji as the example.Journal of Hebei Normal University:Natural Science Edition,2011,35(6):633-639.

    [32] Jiang W G,Zhu X H,Wu J J,Gu L,Ma G B,Liu X C.Retrieval and analysis of coal fire temperature in Wuda Coalfield,Inner Mongolia,China.Chinese Geographical Science,2011,21(2):159-166.

    [33] Cao X L,Zhang L M,Xue D S,Wang D P.The changes in disparity of urban transportation development level rank in China.Acta Geographica Sinica,2007,62(10):1034-1040.

    參考文獻(xiàn):

    [1] 常青,邱瑤,謝苗苗,彭建.基于土地破壞的礦區(qū)生態(tài)風(fēng)險(xiǎn)評(píng)價(jià):理論與方法.生態(tài)學(xué)報(bào),2012,32(16):5164-5174.

    [2] 李晉川,白中科.露天煤礦土地復(fù)墾與生態(tài)重建.北京:科學(xué)出版社,2000.

    [4] 陳輝,劉勁松,曹宇,李雙成,歐陽(yáng)華.生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)研究進(jìn)展.生態(tài)學(xué)報(bào),2006,26(5):1558-1566.

    [6] 付在毅,許學(xué)工.區(qū)域生態(tài)風(fēng)險(xiǎn)評(píng)價(jià).地理科學(xué)進(jìn)展,2001,16(2):267-271.

    [7] 付在毅,許學(xué)工,林輝平,王憲禮.遼河三角洲濕地區(qū)域生態(tài)風(fēng)險(xiǎn)評(píng)價(jià).生態(tài)學(xué)報(bào),2001,21(3):365-373.

    [9] 曾勇.區(qū)域生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)——以呼和浩特市區(qū)為例.生態(tài)學(xué)報(bào),2010,30(3):668-673.

    [10] 潘雅靜,王仰麟,彭建,韓憶楠.礦區(qū)生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)研究述評(píng).生態(tài)學(xué)報(bào),2012,32(20):6566-6574.

    [11] 白中科,趙景逵,李晉川,王文英,盧崇恩,丁新啟,柴書杰,陳建軍.大型露天煤礦生態(tài)系統(tǒng)受損研究——以平碩露天煤礦為例.生態(tài)學(xué)報(bào),1999,19(6):870-875.

    [14] 賈媛,曹玲嫻.煤炭礦區(qū)生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)方法研究.環(huán)境科學(xué)與管理,2011,36(4):177-182.

    [15] 馬蕭.脆弱性礦區(qū)生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)——以勝利東二號(hào)露天礦為例[D].北京:中國(guó)地質(zhì)大學(xué),2011.

    [16] 鄔建國(guó).景觀生態(tài)學(xué),格局、過程、尺度與等級(jí) (第二版).北京:高等教育出版社,2007.

    [17] 謝花林.基于景觀結(jié)構(gòu)和空間統(tǒng)計(jì)學(xué)的區(qū)域生態(tài)風(fēng)險(xiǎn)分析.生態(tài)學(xué)報(bào),2008,28(10):5020-5026.

    [18] 程建龍,陸兆華,范英宏.露天煤礦區(qū)生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)方法.生態(tài)學(xué)報(bào),2004,24(12):2945-2950.

    [19] 胡和兵,劉紅玉,郝敬鋒,安靜.流域景觀結(jié)構(gòu)的城市化影響與生態(tài)風(fēng)險(xiǎn)評(píng)價(jià).生態(tài)學(xué)報(bào),2011,31(12):3432-3440.

    [20] 李謝輝,李景宜.基于GIS的區(qū)域景觀生態(tài)風(fēng)險(xiǎn)分析——以渭河下游河流沿線區(qū)域?yàn)槔?干旱區(qū)研究,2008,25(6):899-903.

    [21] 陳鵬,潘曉玲.干旱區(qū)內(nèi)陸流域區(qū)域景觀生態(tài)風(fēng)險(xiǎn)分析——以阜康三工河流域?yàn)槔?生態(tài)學(xué)雜志,2003,22(4):116-120.

    [22] 王根緒,程國(guó)棟.荒漠綠洲生態(tài)系統(tǒng)的景觀格局分析——景觀空間方法與應(yīng)用.干旱區(qū)研究,1999,16(3):6-11.

    [23] 高賓,李小玉,李玉剛,陳瑋,何興元,齊善忠.基于景觀格局的錦州灣沿海經(jīng)濟(jì)開發(fā)區(qū)生態(tài)風(fēng)險(xiǎn)分析.生態(tài)學(xué)報(bào),2011,31(12):3441-3450.

    [24] 盧遠(yuǎn),蘇文靜,華璀,湯傳勇.左江上游流域景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià).熱帶地理,2010,30(5):496-503.

    [26] 宣國(guó)富,徐建剛,趙靜.基于ESDA的城市社會(huì)空間研究——以上海市中心城區(qū)為例.地理科學(xué),2010,30(1):22-29.

    [28] 孟斌,王勁峰,張文忠,劉旭華.基于空間分析方法的中國(guó)區(qū)域差異研究.地理科學(xué),2005,25(4):393-400.

    [29] 張海峰,白永平,陳瓊,王保宏,牛定煒.基于ESDA-GIS的青海省區(qū)域經(jīng)濟(jì)差異研究.干旱區(qū)地理,2009,32(3):454-461.

    [30] 呂安民,李成名,林宗堅(jiān),史文中.中國(guó)省級(jí)人口增長(zhǎng)率及其空間關(guān)聯(lián)分析.地理學(xué)報(bào),2002,57(2):143-150.

    [31] 黃煥春,李明玉.基于ESDA的長(zhǎng)吉圖開發(fā)先導(dǎo)區(qū)城市內(nèi)部空間結(jié)構(gòu)演化分析——以延吉市為例.河北師范大學(xué)學(xué)報(bào):自然科學(xué)版,2011,35(6):633-639.

    [33] 曹小曙,張利敏,薛德升,王大鵬.中國(guó)城市交通運(yùn)輸發(fā)展水平等級(jí)差異變動(dòng)特征.地理學(xué)報(bào),2007,62(10):1034-1040.

    猜你喜歡
    礦區(qū)景觀區(qū)域
    景觀別墅
    火山塑造景觀
    加納Amanforom礦區(qū)Ⅲ號(hào)隱伏金礦帶的發(fā)現(xiàn)與評(píng)價(jià)
    加納Amanforom礦區(qū)Ⅲ號(hào)隱伏金礦帶的發(fā)現(xiàn)與評(píng)價(jià)
    湖北省??悼h堰邊上礦區(qū)發(fā)現(xiàn)超大型磷礦
    廣東省蕉嶺縣作壁坑礦區(qū)探明超大型銣礦
    包羅萬象的室內(nèi)景觀
    關(guān)于四色猜想
    分區(qū)域
    基于嚴(yán)重區(qū)域的多PCC點(diǎn)暫降頻次估計(jì)
    亚洲美女视频黄频| 欧美日韩在线观看h| 亚洲国产高清在线一区二区三| 欧美人与善性xxx| 国产精品国产av在线观看| 免费看不卡的av| 91狼人影院| 又爽又黄a免费视频| 狠狠精品人妻久久久久久综合| 亚洲不卡免费看| 美女视频免费永久观看网站| 国产精品久久久久成人av| 嫩草影院入口| 黑人高潮一二区| 色吧在线观看| 少妇人妻一区二区三区视频| 精品熟女少妇av免费看| 亚洲久久久国产精品| 一个人看的www免费观看视频| 看十八女毛片水多多多| 国产男女超爽视频在线观看| 黄色日韩在线| 亚洲精品亚洲一区二区| 国产精品成人在线| 亚洲欧美日韩无卡精品| 国产精品一及| 18禁裸乳无遮挡动漫免费视频| 777米奇影视久久| 国国产精品蜜臀av免费| 自拍欧美九色日韩亚洲蝌蚪91 | 99re6热这里在线精品视频| 赤兔流量卡办理| 亚洲无线观看免费| 午夜视频国产福利| 久久精品国产亚洲av涩爱| 免费看光身美女| 亚洲欧美日韩东京热| 成人美女网站在线观看视频| 老女人水多毛片| 久久久久久久亚洲中文字幕| 涩涩av久久男人的天堂| 亚洲精品色激情综合| 视频区图区小说| 亚洲综合色惰| 国产高清国产精品国产三级 | 99久久精品一区二区三区| 成年女人在线观看亚洲视频| 亚洲av中文字字幕乱码综合| 午夜视频国产福利| 欧美高清成人免费视频www| 亚洲三级黄色毛片| 亚洲国产精品国产精品| 国产大屁股一区二区在线视频| 乱码一卡2卡4卡精品| 午夜视频国产福利| 九草在线视频观看| 久久久久久久久久人人人人人人| 亚洲一区二区三区欧美精品| 欧美日韩精品成人综合77777| 婷婷色麻豆天堂久久| 免费播放大片免费观看视频在线观看| av国产免费在线观看| 亚洲欧美成人综合另类久久久| 国产精品成人在线| 亚洲一区二区三区欧美精品| 午夜日本视频在线| 午夜福利在线观看免费完整高清在| 高清午夜精品一区二区三区| 五月开心婷婷网| 亚洲成人av在线免费| 伊人久久精品亚洲午夜| 亚洲精华国产精华液的使用体验| 又黄又爽又刺激的免费视频.| 久久久精品免费免费高清| 国产精品99久久久久久久久| 国产精品一区二区性色av| 男女啪啪激烈高潮av片| 免费播放大片免费观看视频在线观看| 少妇人妻久久综合中文| 久久99热这里只频精品6学生| 精品少妇久久久久久888优播| 国产高清有码在线观看视频| 卡戴珊不雅视频在线播放| 久久午夜福利片| 美女cb高潮喷水在线观看| 日本欧美国产在线视频| 国产精品99久久99久久久不卡 | 熟妇人妻不卡中文字幕| 国产深夜福利视频在线观看| 精品久久久噜噜| 99久久精品一区二区三区| 国产精品不卡视频一区二区| 亚洲色图av天堂| 欧美区成人在线视频| 国产黄片美女视频| 亚洲,一卡二卡三卡| 2021少妇久久久久久久久久久| 久久精品久久精品一区二区三区| 成人无遮挡网站| 51国产日韩欧美| 热re99久久精品国产66热6| 妹子高潮喷水视频| 韩国高清视频一区二区三区| 香蕉精品网在线| 久久久久国产精品人妻一区二区| 高清欧美精品videossex| 亚洲精品aⅴ在线观看| 久久精品熟女亚洲av麻豆精品| 在线观看美女被高潮喷水网站| 中国三级夫妇交换| 国产亚洲最大av| 免费观看的影片在线观看| 免费在线观看成人毛片| 久久久久久人妻| 老女人水多毛片| 尤物成人国产欧美一区二区三区| 香蕉精品网在线| 亚洲欧美一区二区三区国产| 亚洲内射少妇av| 亚洲精品自拍成人| 国产高潮美女av| 18+在线观看网站| 国产在线男女| 九草在线视频观看| 国产女主播在线喷水免费视频网站| 国产真实伦视频高清在线观看| 搡女人真爽免费视频火全软件| 亚洲综合色惰| 国产成人免费观看mmmm| 大码成人一级视频| 成人国产av品久久久| 中国三级夫妇交换| 久久99热这里只有精品18| 亚洲伊人久久精品综合| 爱豆传媒免费全集在线观看| 成人黄色视频免费在线看| 精品人妻视频免费看| 欧美精品人与动牲交sv欧美| 全区人妻精品视频| 亚洲色图av天堂| 国产精品一区二区三区四区免费观看| 国产日韩欧美亚洲二区| 干丝袜人妻中文字幕| 国产淫片久久久久久久久| 免费在线观看成人毛片| av在线老鸭窝| 国产在线一区二区三区精| 亚洲美女搞黄在线观看| 亚洲精品乱码久久久v下载方式| 国产黄色视频一区二区在线观看| 日韩 亚洲 欧美在线| 国产乱来视频区| 日日摸夜夜添夜夜爱| 国产精品精品国产色婷婷| 在线观看一区二区三区| 22中文网久久字幕| 人妻夜夜爽99麻豆av| 一个人看视频在线观看www免费| 亚洲精品一区蜜桃| 亚洲第一av免费看| a级一级毛片免费在线观看| 亚洲精品乱码久久久v下载方式| 国产淫语在线视频| 亚洲欧美日韩东京热| 午夜福利网站1000一区二区三区| 男人舔奶头视频| 99热国产这里只有精品6| 国产黄频视频在线观看| 99久久中文字幕三级久久日本| 夫妻午夜视频| 91在线精品国自产拍蜜月| 女性被躁到高潮视频| 日本黄色片子视频| 国产精品伦人一区二区| 国产精品麻豆人妻色哟哟久久| 国产成人午夜福利电影在线观看| 女性生殖器流出的白浆| 午夜福利在线在线| 黄片wwwwww| 国产精品免费大片| 久久国内精品自在自线图片| 亚洲国产成人一精品久久久| 国产亚洲av片在线观看秒播厂| kizo精华| 成年美女黄网站色视频大全免费 | 天美传媒精品一区二区| 国产精品爽爽va在线观看网站| 国产成人免费无遮挡视频| 午夜激情福利司机影院| 成人二区视频| 久久精品国产a三级三级三级| 尤物成人国产欧美一区二区三区| 伦精品一区二区三区| 纵有疾风起免费观看全集完整版| 一区在线观看完整版| 久久影院123| 精品人妻熟女av久视频| 亚洲精品一区蜜桃| 女人久久www免费人成看片| 久久热精品热| 久久久久久久久大av| 国产成人a∨麻豆精品| 国产中年淑女户外野战色| 在线亚洲精品国产二区图片欧美 | 纵有疾风起免费观看全集完整版| 在线观看免费日韩欧美大片 | 最黄视频免费看| 3wmmmm亚洲av在线观看| 一区二区av电影网| 成人特级av手机在线观看| 亚洲av电影在线观看一区二区三区| a级毛片免费高清观看在线播放| 91精品国产国语对白视频| 久久午夜福利片| 国产成人a∨麻豆精品| 狂野欧美白嫩少妇大欣赏| 我要看黄色一级片免费的| 18禁在线无遮挡免费观看视频| 精品少妇久久久久久888优播| a 毛片基地| 国产淫语在线视频| 欧美三级亚洲精品| 免费播放大片免费观看视频在线观看| 伦理电影免费视频| 在线 av 中文字幕| 极品教师在线视频| 一级毛片aaaaaa免费看小| h日本视频在线播放| 日本黄大片高清| 菩萨蛮人人尽说江南好唐韦庄| 视频中文字幕在线观看| 国产美女午夜福利| 熟女电影av网| 欧美日韩一区二区视频在线观看视频在线| 国产精品无大码| 中国国产av一级| av女优亚洲男人天堂| 亚洲国产毛片av蜜桃av| 大陆偷拍与自拍| 久久久午夜欧美精品| 亚洲av成人精品一二三区| 观看美女的网站| 亚洲色图综合在线观看| 色婷婷久久久亚洲欧美| 香蕉精品网在线| 精品一品国产午夜福利视频| 亚洲精品aⅴ在线观看| 欧美高清成人免费视频www| 女人久久www免费人成看片| 国产午夜精品久久久久久一区二区三区| 成人毛片60女人毛片免费| 亚洲综合精品二区| 免费黄色在线免费观看| 久久婷婷青草| 看十八女毛片水多多多| 青春草视频在线免费观看| 亚洲av综合色区一区| 国产高潮美女av| 国产黄频视频在线观看| 国国产精品蜜臀av免费| 欧美xxxx性猛交bbbb| 久久久欧美国产精品| 亚洲成人中文字幕在线播放| 深爱激情五月婷婷| 久久精品久久久久久久性| 最近中文字幕2019免费版| 精品午夜福利在线看| 美女视频免费永久观看网站| 国产免费福利视频在线观看| 国产免费视频播放在线视频| 久久精品国产鲁丝片午夜精品| 夜夜看夜夜爽夜夜摸| 亚洲国产精品一区三区| 99久久精品一区二区三区| 亚洲欧洲国产日韩| 亚洲精品成人av观看孕妇| 国产男女超爽视频在线观看| 青青草视频在线视频观看| 国产大屁股一区二区在线视频| 日韩av在线免费看完整版不卡| 免费久久久久久久精品成人欧美视频 | 大片免费播放器 马上看| 夜夜爽夜夜爽视频| 高清日韩中文字幕在线| 极品教师在线视频| 亚洲av欧美aⅴ国产| 六月丁香七月| 伦精品一区二区三区| 一级片'在线观看视频| 干丝袜人妻中文字幕| 成年美女黄网站色视频大全免费 | 狂野欧美激情性xxxx在线观看| 嘟嘟电影网在线观看| 欧美精品人与动牲交sv欧美| 天天躁日日操中文字幕| 亚洲欧美日韩东京热| 多毛熟女@视频| 如何舔出高潮| 免费观看av网站的网址| 亚洲成人手机| 免费大片黄手机在线观看| 黑人猛操日本美女一级片| 新久久久久国产一级毛片| 精华霜和精华液先用哪个| 国产黄色视频一区二区在线观看| 欧美高清成人免费视频www| 国产视频内射| 欧美激情极品国产一区二区三区 | 亚洲欧美一区二区三区国产| 少妇的逼水好多| 大香蕉久久网| 一个人免费看片子| 最近中文字幕2019免费版| 亚洲国产精品专区欧美| 少妇猛男粗大的猛烈进出视频| 久久久久国产精品人妻一区二区| 老熟女久久久| 熟女电影av网| 高清日韩中文字幕在线| 亚洲高清免费不卡视频| 熟妇人妻不卡中文字幕| 亚洲最大成人中文| 黄片wwwwww| 噜噜噜噜噜久久久久久91| 国产精品99久久99久久久不卡 | 最近最新中文字幕大全电影3| 街头女战士在线观看网站| 麻豆成人午夜福利视频| 亚洲欧美一区二区三区黑人 | 一级毛片我不卡| 国产片特级美女逼逼视频| av线在线观看网站| 男女无遮挡免费网站观看| 十八禁网站网址无遮挡 | 特大巨黑吊av在线直播| 亚洲美女搞黄在线观看| 在线精品无人区一区二区三 | 国产成人免费无遮挡视频| 欧美精品亚洲一区二区| 尤物成人国产欧美一区二区三区| 精品少妇黑人巨大在线播放| 边亲边吃奶的免费视频| 色婷婷av一区二区三区视频| 免费观看性生交大片5| 亚洲综合精品二区| 久久99热6这里只有精品| 国产精品99久久久久久久久| 成人漫画全彩无遮挡| 99久久精品一区二区三区| 欧美精品国产亚洲| 亚洲第一av免费看| 亚州av有码| 国产大屁股一区二区在线视频| 亚洲国产av新网站| 午夜精品国产一区二区电影| 黄片wwwwww| 极品教师在线视频| 久久久久精品久久久久真实原创| 纯流量卡能插随身wifi吗| 国产精品成人在线| 美女高潮的动态| 纵有疾风起免费观看全集完整版| 日韩亚洲欧美综合| 亚洲av日韩在线播放| 插阴视频在线观看视频| 高清午夜精品一区二区三区| 亚洲国产最新在线播放| 高清午夜精品一区二区三区| 亚洲四区av| 全区人妻精品视频| 美女cb高潮喷水在线观看| 国产男女超爽视频在线观看| 国产在视频线精品| 肉色欧美久久久久久久蜜桃| 国产精品久久久久成人av| 亚洲成人av在线免费| 久久韩国三级中文字幕| 免费大片黄手机在线观看| 最近中文字幕高清免费大全6| 内射极品少妇av片p| 久久ye,这里只有精品| 人人妻人人看人人澡| 激情五月婷婷亚洲| 欧美高清成人免费视频www| 国产视频首页在线观看| 欧美日本视频| 超碰av人人做人人爽久久| 欧美老熟妇乱子伦牲交| 久久综合国产亚洲精品| 看免费成人av毛片| 91精品国产国语对白视频| 日韩,欧美,国产一区二区三区| 天堂8中文在线网| 久久亚洲国产成人精品v| 欧美性感艳星| 日韩成人伦理影院| 免费大片18禁| 国产在线视频一区二区| 欧美成人a在线观看| 夫妻性生交免费视频一级片| 18禁在线无遮挡免费观看视频| 国产白丝娇喘喷水9色精品| 免费观看无遮挡的男女| 一个人看的www免费观看视频| 99久久综合免费| 亚洲av二区三区四区| 成人无遮挡网站| 欧美日韩综合久久久久久| videossex国产| 欧美日韩综合久久久久久| 午夜福利影视在线免费观看| 大话2 男鬼变身卡| 亚洲欧美一区二区三区国产| 亚洲欧美精品专区久久| 如何舔出高潮| 亚洲综合色惰| 免费人成在线观看视频色| 一级a做视频免费观看| 国产欧美日韩精品一区二区| av一本久久久久| 国产成人午夜福利电影在线观看| 人人妻人人添人人爽欧美一区卜 | 激情五月婷婷亚洲| 在线精品无人区一区二区三 | av在线蜜桃| 亚洲av免费高清在线观看| 在线看a的网站| 欧美精品亚洲一区二区| 亚洲欧美精品专区久久| 99热这里只有精品一区| 久久久久久伊人网av| 99热这里只有是精品50| 激情五月婷婷亚洲| 国产在线视频一区二区| 六月丁香七月| 身体一侧抽搐| 国产乱来视频区| 免费观看a级毛片全部| 熟妇人妻不卡中文字幕| 日韩中文字幕视频在线看片 | 啦啦啦中文免费视频观看日本| 黄色配什么色好看| 精品久久久久久久末码| 成人毛片60女人毛片免费| 水蜜桃什么品种好| 国产精品偷伦视频观看了| 亚洲精品456在线播放app| 国产日韩欧美亚洲二区| 日日撸夜夜添| 日韩欧美一区视频在线观看 | 内射极品少妇av片p| 日日摸夜夜添夜夜爱| 亚洲成色77777| 久久国内精品自在自线图片| 天天躁夜夜躁狠狠久久av| 少妇熟女欧美另类| 26uuu在线亚洲综合色| 直男gayav资源| 日本vs欧美在线观看视频 | 国产爱豆传媒在线观看| 国产成人aa在线观看| 亚洲美女黄色视频免费看| 国产精品蜜桃在线观看| 亚洲综合色惰| 国产av一区二区精品久久 | 中文字幕久久专区| 成人影院久久| 国产探花极品一区二区| videossex国产| 两个人的视频大全免费| 99热全是精品| 天天躁夜夜躁狠狠久久av| 日韩亚洲欧美综合| 免费观看在线日韩| av女优亚洲男人天堂| 熟女电影av网| 亚洲欧美清纯卡通| 久久久久久九九精品二区国产| 久久精品久久久久久久性| 一区二区av电影网| 成人一区二区视频在线观看| 日本av手机在线免费观看| 三级经典国产精品| 日本欧美国产在线视频| 国产精品福利在线免费观看| 九九爱精品视频在线观看| 国产精品国产三级国产av玫瑰| 久久婷婷青草| 国产黄片视频在线免费观看| 美女脱内裤让男人舔精品视频| 黄色视频在线播放观看不卡| 日本av免费视频播放| 免费黄网站久久成人精品| 亚洲国产精品专区欧美| 久久毛片免费看一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 午夜福利在线观看免费完整高清在| 亚洲av免费高清在线观看| 国产精品偷伦视频观看了| 伊人久久国产一区二区| 精品一区在线观看国产| 亚洲经典国产精华液单| 熟女人妻精品中文字幕| 亚洲欧美精品专区久久| 一个人免费看片子| 深爱激情五月婷婷| 亚洲国产高清在线一区二区三| 久久精品国产亚洲网站| 少妇人妻 视频| 精品一品国产午夜福利视频| 日本爱情动作片www.在线观看| 人体艺术视频欧美日本| 99热这里只有是精品在线观看| 又大又黄又爽视频免费| 精品一品国产午夜福利视频| 一二三四中文在线观看免费高清| 纵有疾风起免费观看全集完整版| 特大巨黑吊av在线直播| 日韩av免费高清视频| 十八禁网站网址无遮挡 | 美女中出高潮动态图| 18禁在线播放成人免费| 国产有黄有色有爽视频| 韩国av在线不卡| 少妇人妻 视频| a级一级毛片免费在线观看| 久久6这里有精品| 男女免费视频国产| 男人狂女人下面高潮的视频| 自拍欧美九色日韩亚洲蝌蚪91 | 日韩伦理黄色片| 国产色婷婷99| 天美传媒精品一区二区| 天堂中文最新版在线下载| 最近2019中文字幕mv第一页| 亚洲欧美日韩无卡精品| av黄色大香蕉| 久久久久人妻精品一区果冻| 国产男女超爽视频在线观看| 国产乱来视频区| 我的女老师完整版在线观看| 日韩欧美一区视频在线观看 | av在线app专区| 男女国产视频网站| 亚洲av二区三区四区| 亚洲精品第二区| 亚洲国产av新网站| 熟女av电影| 王馨瑶露胸无遮挡在线观看| 亚洲精品乱码久久久久久按摩| 免费黄网站久久成人精品| av一本久久久久| 国产av精品麻豆| 国产69精品久久久久777片| 日本午夜av视频| 极品少妇高潮喷水抽搐| 久久久久久伊人网av| 久久久久精品性色| 成人一区二区视频在线观看| 一级毛片黄色毛片免费观看视频| 亚洲精品第二区| 中文欧美无线码| 国产精品久久久久久久久免| 91aial.com中文字幕在线观看| 免费看日本二区| 九九久久精品国产亚洲av麻豆| 国产淫语在线视频| 激情 狠狠 欧美| 五月开心婷婷网| 日韩欧美精品免费久久| 午夜福利在线观看免费完整高清在| 下体分泌物呈黄色| 久久热精品热| 久久久久性生活片| 国产一区二区三区av在线| 久热久热在线精品观看| 日韩精品有码人妻一区| 高清在线视频一区二区三区| 一本色道久久久久久精品综合| 国产色爽女视频免费观看| av又黄又爽大尺度在线免费看| 精品亚洲乱码少妇综合久久| 美女国产视频在线观看| 亚洲人成网站高清观看| 亚洲av福利一区| 午夜视频国产福利| 亚洲精品乱久久久久久| 日韩欧美一区视频在线观看 | 精品一区二区免费观看| 夜夜爽夜夜爽视频| 激情五月婷婷亚洲| 欧美人与善性xxx| 久久99热这里只频精品6学生| 成人免费观看视频高清| av国产久精品久网站免费入址| 国产成人精品久久久久久| 男人爽女人下面视频在线观看| 国产黄色免费在线视频| 欧美日韩亚洲高清精品| 97超视频在线观看视频| 国产成人freesex在线| 51国产日韩欧美| 26uuu在线亚洲综合色| 在线播放无遮挡| 国产片特级美女逼逼视频| 精品一区二区免费观看| 一个人看视频在线观看www免费| 欧美精品亚洲一区二区| 七月丁香在线播放| www.av在线官网国产| 亚洲国产精品成人久久小说| 视频区图区小说| 熟女av电影| 简卡轻食公司| 一本—道久久a久久精品蜜桃钙片|