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

    基于RS與GIS的長(zhǎng)江三角洲生態(tài)環(huán)境脆弱性綜合評(píng)價(jià)

    2011-09-23 02:13:02徐慶勇陸佩玲李仁強(qiáng)
    環(huán)境科學(xué)研究 2011年1期
    關(guān)鍵詞:長(zhǎng)江三角洲旱澇脆弱性

    徐慶勇,黃 玫,陸佩玲,李仁強(qiáng)

    1.北京林業(yè)大學(xué)林學(xué)院,北京 100083

    2.中國(guó)科學(xué)院地理科學(xué)與資源研究所生態(tài)系統(tǒng)網(wǎng)絡(luò)觀測(cè)與模擬重點(diǎn)實(shí)驗(yàn)室,北京 100101

    基于RS與GIS的長(zhǎng)江三角洲生態(tài)環(huán)境脆弱性綜合評(píng)價(jià)

    徐慶勇1,2,黃 玫2*,陸佩玲1,李仁強(qiáng)2

    1.北京林業(yè)大學(xué)林學(xué)院,北京 100083

    2.中國(guó)科學(xué)院地理科學(xué)與資源研究所生態(tài)系統(tǒng)網(wǎng)絡(luò)觀測(cè)與模擬重點(diǎn)實(shí)驗(yàn)室,北京 100101

    使用空間主成分分析法構(gòu)建評(píng)價(jià)指標(biāo)體系,采用層次分析法確定指標(biāo)權(quán)重,結(jié)合遙感數(shù)據(jù)與地理信息系統(tǒng)軟件,對(duì)長(zhǎng)江三角洲生態(tài)環(huán)境脆弱性進(jìn)行了綜合評(píng)價(jià),并對(duì)脆弱性成因進(jìn)行了分析.結(jié)果表明,長(zhǎng)江三角洲生態(tài)環(huán)境極度和重度脆弱區(qū)主要分布在其中部的太湖流域和浙江中西部,占整個(gè)研究區(qū)的20.10%;輕度和中度脆弱區(qū)遍布于整個(gè)研究區(qū),占55.25%;微度脆弱區(qū)主要分布在江蘇北部和浙江東部,占24.65%.總體來(lái)看,長(zhǎng)江三角洲大部分區(qū)域的生態(tài)環(huán)境屬輕度和中度脆弱.影響長(zhǎng)江三角洲生態(tài)環(huán)境脆弱性的自然因素有≥35℃日數(shù)、旱澇分布、海拔高度、歸一化植被指數(shù)(NDVI)年累加值、景觀多樣性指數(shù)和土壤侵蝕強(qiáng)度;人為因素有人均耕地面積、人均水資源、人均廢水排放量、人均廢氣排放量、化肥施用強(qiáng)度、土地利用變化、“三廢”綜合利用產(chǎn)品產(chǎn)值、人口密度和GDP.極度和重度脆弱區(qū)生態(tài)環(huán)境的主要特征是自然災(zāi)害發(fā)生頻率大、資源匱乏、污染強(qiáng)度大、土壤侵蝕嚴(yán)重和生物多樣性低.

    生態(tài)環(huán)境;脆弱性;遙感;地理信息系統(tǒng);長(zhǎng)江三角洲

    Abstract:The vulnerability of eco-environments in the Yangtze River Delta was assessed in an integrated manner with remote sensing data and geographic information system software combined with spatial principal component analysis.Evaluation index systems were generated,and an analytic hierarchy process(AHP)was used to determine evaluation factor weightings.The causes of vulnerability were analyzed.The results show that the Taihu Lake drainage basin and the western and the middle parts of Zhejiang Province,which are located in themiddle part of the Yangtze River Delta,are themost extremely and seriously vulnerable regions.These regions account for 20.10% of the whole study area.The lightly and moderately vulnerable areas were scattered throughout the entire study area,and accounted for 55.25%of the study area.The slightly vulnerable areas weremainly located in the north of Jiangsu Province and in the east of Zhejiang Province,and accounted for 24.65%of the study area.On the whole,the eco-environment in most areas of the Yangtze River Delta is lightly and moderately vulnerable.The natural factors which influence the eco-environmental vulnerability of the Yangtze River Delta are:days of temperature not less than 35℃,drought and flood frequency,elevation,yearly integrated NDVI,landscape diversity index,and soil erosion intensity.The human factors which influence the eco-environmental vulnerability are:per capita arable land area,per capita water resources,discharge of waste water per capita,per capita exhaust emissions,amount of chemical fertilizer, land use change, wasteutilization product value,population density,GDP,etc.The extremely and seriously vulnerable regions show the characteristics of high frequency of natural disasters,lack of resources,high pollution level,serious soil erosion and low biodiversity.

    Keywords:eco-environmental;vulnerability;remote sensing;geographic information systems;Yangtze River Delta

    隨著全球變暖引起的生態(tài)環(huán)境退化的加劇,生態(tài)環(huán)境脆弱性研究日益引起學(xué)術(shù)界關(guān)注.有學(xué)者就生態(tài)環(huán)境脆弱性的內(nèi)涵、類型、成因、表現(xiàn)、程度和分布等問(wèn)題進(jìn)行了深入研究,獲得了大量研究成果[1-4].在脆弱性評(píng)價(jià)方面,形成了模糊判定分析方法[5]、主成分分析法[6]和綜合評(píng)價(jià)法[7]等系列評(píng)價(jià)方法.使用這些方法獲得了巖溶山區(qū)[8]、喀斯特區(qū)域[9]、岷江上游流域[10]和吉林西部[11]等地區(qū)的生態(tài)環(huán)境脆弱性研究成果.劉燕華等[12-13]對(duì)生態(tài)環(huán)境脆弱性定量評(píng)價(jià)研究進(jìn)行了比較全面的綜述.早期的生態(tài)環(huán)境脆弱性評(píng)價(jià)多集中于指標(biāo)體系的構(gòu)建,因子權(quán)重的確定和評(píng)價(jià)模型的建立;研究單元多局限在縣、市或其他行政單元.隨著遙感(RS)數(shù)據(jù)與地理信息系統(tǒng)(GIS)軟件的廣泛應(yīng)用,生態(tài)環(huán)境脆弱性研究可以在更大范圍內(nèi)更精細(xì)地進(jìn)行,如 WANG等[6,14]把 RS與GIS相結(jié)合,分別以分辨率500 M和100 M的像元作為評(píng)價(jià)單元,評(píng)價(jià)了黃河流域和青藏高原的生態(tài)環(huán)境脆弱性,精細(xì)地反映出區(qū)域內(nèi)的差異.從研究對(duì)象來(lái)看,過(guò)去的研究多集中于我國(guó)典型生態(tài)環(huán)境脆弱區(qū),如西南巖溶山區(qū)、喀斯特區(qū)域、西北干旱半干旱和北方農(nóng)牧交錯(cuò)帶等,這些區(qū)域經(jīng)濟(jì)相對(duì)落后,以東部經(jīng)濟(jì)發(fā)達(dá)地區(qū)作為評(píng)價(jià)對(duì)象的研究目前還不多見(jiàn).

    長(zhǎng)江三角洲作為舉世矚目的經(jīng)濟(jì)高速發(fā)展地區(qū)[15],在我國(guó)經(jīng)濟(jì)生活中占有舉足輕重的地位,已成為拉動(dòng)我國(guó)經(jīng)濟(jì)增長(zhǎng)的重要貢獻(xiàn)區(qū).長(zhǎng)江三角洲對(duì)我國(guó)經(jīng)濟(jì)未來(lái)的走向乃至全球經(jīng)濟(jì)格局都將產(chǎn)生深遠(yuǎn)影響[16].作為海、陸、氣界面交接的敏感地區(qū),長(zhǎng)江三角洲是我國(guó)旱澇成災(zāi)較為嚴(yán)重的地區(qū)之一[17];作為中國(guó)乃至世界最重要的經(jīng)濟(jì)核心區(qū)之一,在全球氣候變暖的大背景下,隨著工業(yè)化和城市化進(jìn)程的加快,在高強(qiáng)度的人類活動(dòng)影響和干預(yù)下,長(zhǎng)江三角洲生態(tài)環(huán)境呈現(xiàn)出了脆弱性的特征,如洪澇威脅加重、地質(zhì)災(zāi)害增加、酸雨頻率高[18]、水土流失嚴(yán)重和生物多樣性受到威脅[19]等.長(zhǎng)江三角洲已成為我國(guó)新的生態(tài)環(huán)境脆弱帶[20],其生態(tài)環(huán)境脆弱性程度及發(fā)展趨勢(shì)將影響區(qū)域的可持續(xù)發(fā)展.而生態(tài)環(huán)境脆弱性評(píng)價(jià)的目的是為了維持系統(tǒng)的持續(xù)發(fā)展,減輕外部脅迫對(duì)系統(tǒng)的不利影響和為退化系統(tǒng)的綜合整治提供策略依據(jù)[12],故在該區(qū)域進(jìn)行生態(tài)環(huán)境脆弱性評(píng)價(jià)是非常必要的.然而,以往關(guān)于長(zhǎng)江三角洲的研究,大都偏于經(jīng)濟(jì)結(jié)構(gòu)調(diào)整、經(jīng)濟(jì)協(xié)作、區(qū)域規(guī)劃及一體化、資源開(kāi)發(fā)和環(huán)境污染治理等,鮮見(jiàn)對(duì)長(zhǎng)江三角洲生態(tài)環(huán)境脆弱性的研究報(bào)道.

    筆者嘗試結(jié)合RS與GIS,以分辨率1 km的像元作為評(píng)價(jià)單元,用定量化的方法對(duì)長(zhǎng)江三角洲生態(tài)環(huán)境脆弱性進(jìn)行綜合評(píng)價(jià),以期為區(qū)域生態(tài)環(huán)境的治理和建設(shè)提供科學(xué)依據(jù).

    1 研究方法

    1.1 研究區(qū)概況

    研究區(qū)以2008年國(guó)務(wù)院常務(wù)會(huì)議審議并原則通過(guò)的《關(guān)于進(jìn)一步推進(jìn)長(zhǎng)江三角洲地區(qū)改革開(kāi)放和經(jīng)濟(jì)社會(huì)發(fā)展的指導(dǎo)意見(jiàn)》中確定的“兩省一市”(即江蘇省、浙江省和上海市)的區(qū)域范圍為對(duì) 象,介 于 27°57′N ~35°20′N,116°18′E ~122°45′E之間;面積21.07×104km2,約占全國(guó)陸地面積的2%(見(jiàn)圖1);人口1.38×108人(2008年),約占全國(guó)人口總數(shù)的10%,是我國(guó)人口最稠密的地區(qū)之一.長(zhǎng)江三角洲地區(qū)水網(wǎng)密布,河道縱橫,湖泊眾多.地形以平原和丘陵山地為主,平原約占陸域面積的47.7%,主要分布在江蘇南部、浙江東北部和上海;丘陵山地約占40.8%,主要分布在江蘇西南部和北部、浙江中部和西南部[21-25].氣候?qū)俦眮啛釒Ъ撅L(fēng)氣候,寒暑變化顯著、四季分明.水、熱、光條件優(yōu)越,土壤類型多樣,質(zhì)地優(yōu)良,通透性好,養(yǎng)分充足,是我國(guó)著名的高產(chǎn)穩(wěn)產(chǎn)區(qū)[26].植被類型為常綠闊葉與落葉闊葉混交林.

    1.2 數(shù)據(jù)收集

    研究所用基礎(chǔ)數(shù)據(jù)包括空間數(shù)據(jù)和屬性數(shù)據(jù).空間數(shù)據(jù)包括:2002年4月—2003年4月研究區(qū)歸一化植被指數(shù)(NDVI)影像圖及1985年土地利用柵格圖,均由中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心提供;2005年土地利用柵格圖、旱澇分布矢量圖,2004年地區(qū)生產(chǎn)總值(GDP)、人口密度和海拔高度柵格圖,均由國(guó)家科技基礎(chǔ)條件平臺(tái)建設(shè)項(xiàng)目——地球系統(tǒng)科學(xué)數(shù)據(jù)共享網(wǎng)(www.geodata. cn)提供.屬性數(shù)據(jù)包括:1981—2000年氣象臺(tái)站的日均氣溫、日均降水量數(shù)據(jù),均由國(guó)家氣象局提供;景觀多樣性指數(shù)數(shù)據(jù)來(lái)源于1:1 000 000中國(guó)植被圖[27];土壤侵蝕強(qiáng)度數(shù)據(jù)來(lái)源于1∶15 000 000中國(guó)土壤侵蝕圖[28];人均耕地面積、人均水資源、人均廢水排放量、人均廢氣排放量、人均固體廢物產(chǎn)生量、化肥施用強(qiáng)度、人均“三廢”綜合利用產(chǎn)品產(chǎn)值和農(nóng)民人均純收入數(shù)據(jù),均來(lái)自2003—2007年江蘇、浙江和上海三省市的統(tǒng)計(jì)年鑒、環(huán)境狀況公報(bào)、國(guó)民經(jīng)濟(jì)和社會(huì)發(fā)展統(tǒng)計(jì)公報(bào).

    圖1 研究區(qū)位置Fig.1 The study area map

    1.3 數(shù)據(jù)整理

    首先需要將屬性數(shù)據(jù)轉(zhuǎn)化為空間數(shù)據(jù).對(duì)于無(wú)X,Y坐標(biāo)的屬性數(shù)據(jù),在ArcGIS中將屬性數(shù)據(jù)信息輸入研究區(qū)行政區(qū)劃矢量圖的屬性表,實(shí)現(xiàn)屬性數(shù)據(jù)的空間化;對(duì)于有X,Y坐標(biāo)的屬性數(shù)據(jù),使用AddXYData方法實(shí)現(xiàn)空間化.再將矢量圖轉(zhuǎn)換成柵格圖,柵格大小為1 000 m×1 000 m,使用統(tǒng)一的Krasovsky橢球體坐標(biāo)和Albers投影.

    土壤侵蝕強(qiáng)度數(shù)據(jù)通過(guò)對(duì)1∶15 000 000中國(guó)土壤侵蝕圖的數(shù)字化得到;土地利用變化數(shù)據(jù)通過(guò)比較1985年與2005年的土地利用數(shù)據(jù)得到.

    NDVI年累加值的大小表明植被冠層吸收的光合輻射的多少,能反映植被初級(jí)生產(chǎn)力的大小.參照文獻(xiàn)[29]NDVI年累加值計(jì)算方法,計(jì)算了NDVI年累加值(AI),計(jì)算公式:

    式中,n為每年合成的NDVI影像圖數(shù)量,幅;NDVIi為一年中第i幅NDVI影像圖;Ti為一個(gè)波段天數(shù)(10 d).

    景觀多樣性是生物多樣性的重要組成部分.在較大的時(shí)空尺度上,景觀多樣性構(gòu)成了其他層次生物多樣性的背景,并制約著這些層次生物多樣性的時(shí)空格局及其變化過(guò)程[30].研究區(qū)的生物多樣性狀況用景觀多樣性指數(shù)(H)表示,計(jì)算公式[31]:

    式中,Pi為景觀類型i所占面積的比例,%;m為景觀類型的數(shù)目,個(gè).

    1.4 指標(biāo)體系構(gòu)建方法

    影響生態(tài)環(huán)境脆弱性的指標(biāo)有很多,如果把所有指標(biāo)都加以考慮,將會(huì)得到龐大的指標(biāo)體系,且各指標(biāo)間有很強(qiáng)的相關(guān)性,致使評(píng)價(jià)精度降低.主成分分析法是運(yùn)用統(tǒng)計(jì)分析原理將原來(lái)眾多具有一定相關(guān)性指標(biāo)重新組合成一組新的相互無(wú)關(guān)的綜合指標(biāo),而原指標(biāo)所包含的大量信息仍能得到保持的一種統(tǒng)計(jì)方法[32].空間主成分分析是在GIS軟件的支持下,在ARC/INFO的GRID模塊中使用PRINCOMP函數(shù),通過(guò)將原始空間軸旋轉(zhuǎn),把相關(guān)的多變量空間數(shù)據(jù)轉(zhuǎn)化為少數(shù)不相關(guān)的綜合指標(biāo)[33].使用該方法對(duì)評(píng)價(jià)指標(biāo)進(jìn)行篩選.

    1.5 評(píng)價(jià)指標(biāo)的量化分級(jí)方法

    根據(jù)各評(píng)價(jià)指標(biāo)特征,采用3種分級(jí)方法:①根據(jù)數(shù)據(jù)的統(tǒng)計(jì)特征,采用聚類分析方法分級(jí);②在綜合考慮研究區(qū)的自然地理狀況、資源環(huán)境現(xiàn)狀和社會(huì)經(jīng)濟(jì)發(fā)展水平的基礎(chǔ)上,參照國(guó)家或國(guó)際標(biāo)準(zhǔn)并結(jié)合聚類分析方法分級(jí);③參照有關(guān)文獻(xiàn)分級(jí).

    聚類分析(Cluster Analysis)是根據(jù)事物本身的特性研究分類問(wèn)題的一種多元統(tǒng)計(jì)方法,其基本思想是根據(jù)樣品或指標(biāo)之間的距離或相似系數(shù)進(jìn)行分類[32].研究采用R-型聚類分析的最短距離法將指標(biāo)分類,使用歐氏距離公式計(jì)算各指標(biāo)兩兩數(shù)值間的距離,將距離較短的數(shù)值歸為一類.這一計(jì)算過(guò)程使用社會(huì)科學(xué)統(tǒng)計(jì)軟件包(SPSS)完成.

    采用聚類分析方法分級(jí)的指標(biāo)有≥35℃日數(shù),海拔高度,NDVI年累加值,景觀多樣性指數(shù),人均耕地面積,人均廢水排放量,人均廢氣排放量,“三廢”綜合利用產(chǎn)品產(chǎn)值,人口密度和GDP.旱澇分布指標(biāo)根據(jù)研究區(qū)旱澇發(fā)生頻率分級(jí):無(wú)旱澇區(qū)為第1級(jí);少旱澇區(qū)為第2級(jí);次多旱澇區(qū)為第3級(jí);多旱澇區(qū)為第4級(jí);嚴(yán)重旱澇區(qū)為第5級(jí).人均水資源和化肥施用強(qiáng)度指標(biāo)參照國(guó)際標(biāo)準(zhǔn)并結(jié)合聚類分析方法分級(jí),國(guó)際公認(rèn)缺水線為1 000 m3/人,化肥施用強(qiáng)度國(guó)際標(biāo)準(zhǔn)為<25 g/m2.據(jù)此,人均水資源小于1 000的數(shù)值歸為一類,大于1 000的數(shù)值用聚類分析方法分級(jí),化肥施用強(qiáng)度小于25的數(shù)值歸為一類,大于25的數(shù)值用聚類分析方法分級(jí).土地利用變化指標(biāo)根據(jù)土地利用變化的大小和方向分級(jí):首先參照黃方等[33]的土地利用分級(jí)標(biāo)準(zhǔn),把內(nèi)陸水體、近海濕地、河湖灘地和沼澤地分為第1級(jí);有林地分為第2級(jí);中覆蓋度草地、灌木林和疏林地分為第3級(jí);低覆蓋度草地、旱地、水澆地和水田分為第4級(jí);城鎮(zhèn)建設(shè)用地、農(nóng)村聚落、裸巖和裸地分為第5級(jí).然后,將兩期土地利用數(shù)據(jù)進(jìn)行比較,再分級(jí):2005年較1985年土地利用級(jí)別低3級(jí)以上的像元定為土地利用變化1級(jí)(如土地利用級(jí)別由1985年的第5級(jí)變?yōu)?005年的第1或第2級(jí),或者由1985年的第4級(jí)變?yōu)?005年的第1級(jí));低1級(jí)或2級(jí)的像元定為土地利用變化2級(jí);土地利用級(jí)別沒(méi)有變化的定為第3級(jí);2005年高1級(jí)或2級(jí)的像元定為土地利用變化4級(jí);高3級(jí)以上的定為土地利用變化5級(jí).土壤侵蝕強(qiáng)度指標(biāo)依照文獻(xiàn)[28]的分級(jí)方法并考慮研究區(qū)的實(shí)際狀況得到:無(wú)明顯侵蝕為第1級(jí);微度侵蝕為第2級(jí);輕度侵蝕為第3級(jí);中度侵蝕為第4級(jí);強(qiáng)度侵蝕為第5級(jí).生態(tài)環(huán)境脆弱度隨著劃分等級(jí)的增加而增加.

    1.6 脆弱性指數(shù)計(jì)算及脆弱度分類方法

    參照文獻(xiàn)[14]的脆弱生態(tài)環(huán)境定量評(píng)價(jià)方法,對(duì)研究區(qū)的每個(gè)像元計(jì)算脆弱性指數(shù)(EVI),計(jì)算公式:式中,wi為評(píng)價(jià)指標(biāo)的權(quán)重;fi為評(píng)價(jià)指標(biāo)的等級(jí).

    評(píng)價(jià)指標(biāo)權(quán)重(wi)通過(guò)層次分析法確定.層次分析法是一種定性與定量相結(jié)合的多目標(biāo)決策分析方法,基本原理是:將一個(gè)復(fù)雜問(wèn)題看成一個(gè)系統(tǒng),根據(jù)系統(tǒng)內(nèi)部因素之間的隸屬關(guān)系,將一個(gè)復(fù)雜問(wèn)題的各種要素轉(zhuǎn)化為有條理的有序?qū)哟?,并以同一層次的各種要素按照上一層要素為準(zhǔn)則,構(gòu)造判斷矩陣,進(jìn)行兩兩判斷比較,計(jì)算出各要素的權(quán)重[34],并進(jìn)行一致性檢驗(yàn),當(dāng)CR(隨機(jī)一致性比率)<0.1時(shí),判斷矩陣具有滿意的一致性,否則就需要調(diào)整,直到達(dá)到滿意的一致性為止.使用Excel軟件作層次分析,計(jì)算了各評(píng)價(jià)指標(biāo)的權(quán)重.

    在GIS軟件支持下,使用式(3)計(jì)算每個(gè)像元的EVI,得到EVI的空間分布圖.采用自然斷點(diǎn)法(Nature Breaks Classification)將EVI分成5類,分別為微度脆弱、輕度脆弱、中度脆弱、重度脆弱和極度脆弱5個(gè)脆弱性等級(jí).自然斷點(diǎn)法是一種根據(jù)數(shù)值統(tǒng)計(jì)分布規(guī)律分級(jí)和分類的統(tǒng)計(jì)方法,它能使類與類之間的不同最大化[6].

    2 結(jié)果與分析

    2.1 長(zhǎng)江三角洲生態(tài)環(huán)境脆弱性評(píng)價(jià)指標(biāo)

    在詳細(xì)分析長(zhǎng)江三角洲自然、經(jīng)濟(jì)、社會(huì)和環(huán)境狀況的基礎(chǔ)上,借鑒國(guó)內(nèi)外脆弱性評(píng)價(jià)指標(biāo)的內(nèi)容,選擇兩大類(即自然成因和人為成因指標(biāo))共20個(gè)指標(biāo)作為初選指標(biāo).這些指標(biāo)是:≥10℃積溫,≥35℃日數(shù),日降水量,大暴雨日數(shù),旱澇分布,海拔高度,NDVI年累加值,景觀多樣性指數(shù),土壤侵蝕強(qiáng)度,人均耕地面積,化肥施用強(qiáng)度,土地利用變化,“三廢”綜合利用產(chǎn)品產(chǎn)值,人口密度,GDP和農(nóng)民人均純收入.評(píng)價(jià)指標(biāo)的選取主要基于長(zhǎng)江三角洲的生態(tài)環(huán)境問(wèn)題和脆弱性特征,由于該區(qū)域旱澇災(zāi)害頻發(fā),故選取≥35℃日數(shù)、≥10℃積溫、日降水量、大暴雨日數(shù)、旱澇分布和NDVI年累加值等指標(biāo)反映這一生態(tài)環(huán)境問(wèn)題.長(zhǎng)江三角洲主要的地質(zhì)災(zāi)害——地面沉降和地面裂縫,主要是由污染和過(guò)度開(kāi)采地下水導(dǎo)致,酸雨頻率高則主要是由污染導(dǎo)致,選取人均水資源、人均廢水排放量、人均廢氣排放量、人均固體廢物產(chǎn)生量、化肥施用強(qiáng)度、人口密度(人口密度與資源環(huán)境的壓力正相關(guān))和人均耕地面積等指標(biāo)反映這2種生態(tài)環(huán)境問(wèn)題.水土流失、海岸帶蝕退和生物多樣性受到威脅等生態(tài)環(huán)境問(wèn)題,通過(guò)土壤侵蝕強(qiáng)度、海拔高度、景觀多樣性指數(shù)和土地利用變化指標(biāo)來(lái)反映.“三廢”綜合利用產(chǎn)品產(chǎn)值,GDP和農(nóng)民人均純收入3項(xiàng)指標(biāo)表明該區(qū)域的經(jīng)濟(jì)發(fā)展?fàn)顩r和對(duì)環(huán)保的投入.由于評(píng)價(jià)指標(biāo)間具有相關(guān)性,通過(guò)主成分分析,把相關(guān)系數(shù)達(dá)0.9(含)以上的指標(biāo)去掉一個(gè),由此得到15個(gè)相互獨(dú)立的指標(biāo),其中自然因素指標(biāo)有≥35℃日數(shù),旱澇分布,海拔高度,NDVI年累加值,景觀多樣性指數(shù)和土壤侵蝕強(qiáng)度;人為因素指標(biāo)有人均耕地面積、人均水資源、人均廢水排放量、人均廢氣排放量、化肥施用強(qiáng)度、土地利用變化、“三廢”綜合利用產(chǎn)品產(chǎn)值、人口密度和GDP.

    2.2 評(píng)價(jià)指標(biāo)的量化分級(jí)

    根據(jù)上述量化分級(jí)方法,依參評(píng)因子對(duì)環(huán)境脆弱性的正向影響程度,將脆弱度從低到高分為5級(jí)并賦值:1~5分別為微度、輕度、中度、重度和極度脆弱(見(jiàn)表1).

    表1 長(zhǎng)江三角洲生態(tài)環(huán)境脆弱性評(píng)價(jià)指標(biāo)量化分級(jí)Table 1 The classification of eco-environmental vulnerability factors in Yangtze River Delta

    2.3 評(píng)價(jià)因子權(quán)重

    使用層次分析法構(gòu)造判斷矩陣,并經(jīng)過(guò)一致性檢驗(yàn),得到各評(píng)價(jià)因子權(quán)重(計(jì)算方法見(jiàn)文獻(xiàn)[35]),如表2所示.其中,旱澇分布指標(biāo)權(quán)重最大,為0.196;其次為≥35℃日數(shù)和海拔高度指標(biāo),權(quán)重為0.138;GDP和景觀多樣性指數(shù)指標(biāo)權(quán)重最小,為0.014.

    2.4 生態(tài)環(huán)境脆弱性指數(shù)及脆弱度分類

    由式(3)計(jì)算得到長(zhǎng)江三角洲生態(tài)環(huán)境脆弱性指數(shù)的空間分布,整個(gè)區(qū)域脆弱性指數(shù)值在1.7~4.3之間,平均值為 2.828.根據(jù)自然斷點(diǎn)法,將脆弱性指數(shù)分成 5類:1.7~<2.5,2.5~<2.9,2.9~<3.3,3.3~<3.7和3.7~4.3,分別代表微度、輕度、中度、重度和極度脆弱5個(gè)脆弱度等級(jí).圖2為長(zhǎng)江三角洲生態(tài)環(huán)境脆弱度空間分布.由圖2可以看出,長(zhǎng)江三角洲生態(tài)環(huán)境脆弱性大體上呈現(xiàn)出中部重、南北部輕的特征.極度和重度脆弱區(qū)主要分布在太湖流域和浙江中西部;微度脆弱區(qū)主要分布在江蘇北部和浙江東部;輕度和中度脆弱區(qū)遍布于整個(gè)研究區(qū).

    表2 評(píng)價(jià)指標(biāo)的權(quán)重Table 2 The weights of evaluation factors

    圖2 長(zhǎng)江三角洲生態(tài)環(huán)境脆弱度空間分布Fig.2 The spatial distribution of the eco-environmental vulnerability in Yangtze River Delta

    研究區(qū)各脆弱度等級(jí)分布面積和所占比例如表3所示.輕度脆弱區(qū)(28.92%)和中度脆弱區(qū)(26.33%)占整個(gè)長(zhǎng)江三角洲面積的一半多(55.25%),微度、重度和極度脆弱區(qū)分別占24.65%,11.82%和8.28%,因此,總體來(lái)看,長(zhǎng)江三角洲大部分區(qū)域的生態(tài)環(huán)境屬輕度和中度脆弱.

    表3 長(zhǎng)江三角洲各等級(jí)脆弱度分布面積及所占比例Table 3 The total area and the proportion of different eco-environmental vulnerability level in Yangtze River Delta

    2.5 脆弱性成因分析

    自然因素是脆弱生態(tài)環(huán)境形成的內(nèi)因,人為因素是外因,是觸發(fā)性因子.屬于極度或重度脆弱的太湖流域從自然因素看,脆弱性表現(xiàn)在:洪澇災(zāi)害發(fā)生頻率大;NDVI年累加值低,小于90;景觀多樣性指數(shù)低,大部分地區(qū)小于 0.5.從人為因素看,太湖流域的脆弱性主要源于:①污染較嚴(yán)重.太湖流域廢水、廢氣排放量均較多,5年平均廢水排放量大于53.5 t/人,廢氣排放量大于3.6×104m3/人.由于污染水生動(dòng)植物大量死亡,生物多樣性下降,并且因水污染造成的水質(zhì)型缺水問(wèn)題已十分突出.②土地利用變化大.1985—2005年太湖流域大部分土地由濕地和林地變?yōu)榈透采w度草地和城鎮(zhèn)建設(shè)用地,城鎮(zhèn)建設(shè)用地增加,導(dǎo)致產(chǎn)生大量不透水層,致使產(chǎn)水量增加,大大降低流域的蓄水能力.③對(duì)資源的過(guò)度開(kāi)發(fā)造成對(duì)環(huán)境的破壞.由于過(guò)度開(kāi)采地下水,太湖流域人均水資源為整個(gè)研究區(qū)最少,5年平均少于1 000 m3/人,低于國(guó)際公認(rèn)缺水線,對(duì)地下水的過(guò)度開(kāi)采引發(fā)了地面沉降、地面裂縫等生態(tài)環(huán)境問(wèn)題.

    屬于極度脆弱的浙江中西部,從自然因素看其脆弱性表現(xiàn)在:浙江中部為嚴(yán)重易旱區(qū),高溫日數(shù)多,20年平均≥35℃日數(shù)多于16.6 d.浙江西部為嚴(yán)重易澇區(qū),洪澇災(zāi)害發(fā)生頻率大;土壤侵蝕較嚴(yán)重,大部分地區(qū)為中度侵蝕.從人為因素看,隨著經(jīng)濟(jì)的發(fā)展,浙江中西部工業(yè)污染較嚴(yán)重,浙江西部廢水排放量最多,5年平均大于207 t/人;廢氣排放量也較多,5年平均 6.6×104~15×104m3/人;土地利用變化大,1985—2005年浙江中部大部分土地由林地、中覆蓋度草地變?yōu)榈透采w度草地或旱地、水澆地,植被嚴(yán)重破壞,土壤侵蝕加重.

    屬于微度脆弱的江蘇北部和浙江東部,從自然因素看,該地區(qū)極少發(fā)生旱澇災(zāi)害;極端高溫日數(shù)少,20年平均≥35℃日數(shù)少于16.6 d;NDVI年累加值較高,大于120.從人為因素看,江蘇北部和浙江東部廢水廢氣排放均較少,5年平均廢水排放少于53.5 t/人,廢氣少于3.6×104m3/人.浙江東部土地利用類型為有林地,水土流失少,土壤肥力高;人均水資源較多,5年平均大于1 200 m3/人.江蘇北部大部分地區(qū)人均耕地面積較多,5年平均大于750 m2/人.

    由以上分析可知,在自然因素和人為因素的疊加影響下,太湖流域和浙江中西部為極度或重度脆弱區(qū).而江蘇北部和浙江東部由于自然脆弱性較低,人為干擾相對(duì)較弱,故該區(qū)為微度脆弱區(qū).

    3 結(jié)論

    a.在自然因素和人為因素的疊加影響下,長(zhǎng)江三角洲生態(tài)環(huán)境脆弱性大體呈現(xiàn)出中部重、南北部輕的特征.

    b.極度和重度脆弱區(qū)主要分布在中部的太湖流域以及浙江中西部地區(qū),分別占整個(gè)研究區(qū)的8.28%和11.82%;輕度和中度脆弱區(qū)遍布于整個(gè)研究區(qū),占55.25%;微度脆弱區(qū)主要分布在江蘇北部和浙江東部,占24.65%.總體來(lái)看,長(zhǎng)江三角洲大部分區(qū)域的生態(tài)環(huán)境屬輕度和中度脆弱.

    c.極度和重度脆弱區(qū)生態(tài)環(huán)境的主要特征是旱澇災(zāi)害發(fā)生頻率大、資源匱乏、污染強(qiáng)度大、土壤侵蝕嚴(yán)重和生物多樣性低.

    4 討論

    根據(jù)長(zhǎng)江三角洲生態(tài)環(huán)境脆弱性特征及不同脆弱度分布,對(duì)該區(qū)域生態(tài)環(huán)境保護(hù)的建議是:在極度和重度脆弱區(qū),重點(diǎn)治理,大力建設(shè).依靠對(duì)災(zāi)害防治、環(huán)境基礎(chǔ)設(shè)施建設(shè)以及環(huán)保資金的高投入,降低旱澇災(zāi)害發(fā)生頻率,減少人類活動(dòng)對(duì)濕地的圍墾和對(duì)林地的破壞,控制和治理污染,提高對(duì)工業(yè)“三廢”和生活污水、生活垃圾的處理率和利用率,同時(shí)工程修復(fù)和生態(tài)修復(fù)相結(jié)合,減緩目前該區(qū)生態(tài)環(huán)境十分脆弱的態(tài)勢(shì).在中度和輕度脆弱區(qū),生態(tài)環(huán)境形勢(shì)不是十分嚴(yán)峻,但仍要高度重視生態(tài)環(huán)境保護(hù)和建設(shè),減少人為干擾和破壞,防止生態(tài)環(huán)境向高一級(jí)脆弱類型轉(zhuǎn)化.在微度脆弱區(qū),生態(tài)環(huán)境較好,但仍不能忽視保護(hù),堅(jiān)持資源開(kāi)發(fā)前的生態(tài)環(huán)境評(píng)價(jià)并做好保護(hù)規(guī)劃,防止對(duì)生態(tài)環(huán)境造成破壞.

    然而研究結(jié)果還具有一定的不確定性.這種不確定性一方面來(lái)源于屬性數(shù)據(jù)的空間化,如氣象數(shù)據(jù)、社會(huì)經(jīng)濟(jì)數(shù)據(jù)等是按行政區(qū)劃分的,將這些屬性數(shù)據(jù)轉(zhuǎn)換為空間數(shù)據(jù)會(huì)對(duì)評(píng)價(jià)結(jié)果的精確性產(chǎn)生影響;另一方面,在評(píng)價(jià)指標(biāo)選取、指標(biāo)權(quán)重的確定等方面帶有一定的主觀性,會(huì)對(duì)評(píng)價(jià)結(jié)果的準(zhǔn)確度產(chǎn)生影響.探索一種更加科學(xué)和完善的定量評(píng)價(jià)經(jīng)濟(jì)高度發(fā)達(dá)地區(qū)生態(tài)環(huán)境脆弱性的方法將是未來(lái)的研究?jī)?nèi)容.

    [1]劉燕華.中國(guó)脆弱環(huán)境類型劃分與指標(biāo)[C]//生態(tài)環(huán)境綜合整治與恢復(fù)技術(shù)研究(第二集).北京:北京科學(xué)技術(shù)出版社,1995:8-17.

    [2]趙躍龍.中國(guó)脆弱生態(tài)環(huán)境類型分布及其綜合整治[M].北京:中國(guó)環(huán)境科學(xué)出版社,1999.

    [3]申元村.中國(guó)脆弱環(huán)境區(qū)劃的初步研究[C]//生態(tài)環(huán)境綜合整治與恢復(fù)技術(shù)研究(第二集).北京:北京科學(xué)技術(shù)出版社,1995:69-76.

    [4]張建平.云南元謀干熱河谷脆弱環(huán)境及其退化[C]//生態(tài)環(huán)境綜合整治與恢復(fù)技術(shù)研究(第二集).北京:北京科學(xué)技術(shù)出版社,1995:140-145.

    [5]TRAN LIEM T, KNIGHT C GREGORY, O’NEILL ROBERTV,et al. Fuzzy decision analysis for integrated environmental vulnerability assessment of the mid-atlantic region[J].Environ Manage,2002,29(6):845-859.

    [6]WANG S Y,LIU J S,YANG C J.Eco-environmental vulnerability evaluation in the Yellow River Basin,China[J]. Pedosphere,2008,18(2):171-182.

    [7]LI X L,MIN M,TAN C F.The functional assessment of agriculture ecosystems in Hubei Province, China[J]. Ecological Modelling,2005,187:352-360.

    [8]官冬杰,蘇維詞.基于 GIS重慶巖溶地區(qū)生態(tài)環(huán)境脆弱度評(píng)價(jià)[J].中國(guó)巖溶,2006,25(3):211-218.

    [9]張笑楠,王克林,張偉,等.桂西北喀斯特區(qū)域生態(tài)環(huán)境脆弱性[J].生態(tài)學(xué)報(bào),2009,29(2):749-757.

    [10]姚建,丁晶,艾南山.岷江上游生態(tài)脆弱性評(píng)價(jià)[J].長(zhǎng)江流域資源與環(huán)境,2004,13(4):380-383.

    [11]楊育武,湯潔,麻素挺.脆弱生態(tài)環(huán)境指標(biāo)庫(kù)的建立及其定量評(píng)價(jià)[J].環(huán)境科學(xué)研究,2002,15(4):46-49.

    [12]劉燕華,李秀彬.脆弱生態(tài)環(huán)境與可持續(xù)發(fā)展[M].北京:商務(wù)印書(shū)館,2001:98-121.

    [13]趙躍龍,張玲娟.脆弱生態(tài)環(huán)境定量評(píng)價(jià)方法的研究[J].地理科學(xué)進(jìn)展,1998,17(1):67-72.

    [14]WANG X D,ZHONG X H,LIU S Z,et al.Regional assessment of environmental vulnerability in the Tibetan Plateau:development and app lication of a new method[J]. Journal of Arid Environments,2008(72):1929-1939.

    [15]朱帥,馬建中,王堰,等.長(zhǎng)江三角洲地區(qū)春季臭氧異常高值的數(shù)值模擬研究[J].環(huán)境科學(xué)研究,2006,19(6):1-8.

    [16]謝天成,朱曉華.1978年以來(lái)長(zhǎng)江三角洲地區(qū)經(jīng)濟(jì)格局動(dòng)態(tài)變化研究[J].中國(guó)人口·資源與環(huán)境,2008,18(5):7-11.

    [17]蔣薇,宋連春,王式功,等.長(zhǎng)江三角洲夏季降水異常及氣候成因[J].氣象科學(xué),2009,29(3):355-361.

    [18]王志憲,虞孝感,徐科峰,等.長(zhǎng)江三角洲地區(qū)可持續(xù)發(fā)展的態(tài)勢(shì)與對(duì)策[J].地理學(xué)報(bào),2005,60(3):381-390.

    [19]鄭茂坤,駱永明,趙其國(guó),等.應(yīng)用DPSIR體系解決長(zhǎng)江、珠江三角洲地區(qū)環(huán)境問(wèn)題的初步思考[J].土壤,2006,38 (5):662-666.

    [20]王艷,柴發(fā)合,劉厚鳳,等.長(zhǎng)江三角洲地區(qū)大氣污染物水平輸送場(chǎng)特征分析[J].環(huán)境科學(xué)研究,2008,21(1):22-29.

    [21]浙江省統(tǒng)計(jì)局,國(guó)家統(tǒng)計(jì)局浙江調(diào)查總隊(duì).浙江統(tǒng)計(jì)年鑒: 2009[Z].北京:中國(guó)統(tǒng)計(jì)出版社,2010.

    [22]江蘇省統(tǒng)計(jì)局.江蘇統(tǒng)計(jì)年鑒:2009[Z].北京:中國(guó)統(tǒng)計(jì)出版社,2010.

    [23]上海市統(tǒng)計(jì)局.上海統(tǒng)計(jì)年鑒:2009[Z].北京:中國(guó)統(tǒng)計(jì)出版社,2010.

    [24]李松平,嚴(yán)力蛟,婁偉平,等.浙江省農(nóng)業(yè)氣象災(zāi)害特點(diǎn)及防災(zāi)措施[J].中國(guó)農(nóng)學(xué)通報(bào),2006,22(9):483-486.

    [25]嚴(yán)新明,童星.江蘇省自然災(zāi)害風(fēng)險(xiǎn)管理研究[J].江南大學(xué)學(xué)報(bào):人文社會(huì)科學(xué)版,2006(5):26-31.

    [26]彭補(bǔ)拙,高中貴.長(zhǎng)江三角洲地區(qū)土地利用變化及對(duì)策研究[J].第四紀(jì)研究,2004,24(2):506-511.

    [27]中國(guó)科學(xué)院中國(guó)植被圖編輯委員會(huì).中國(guó)植被圖集[M].北京:科學(xué)出版社,2001.

    [28]廖克.中華人民共和國(guó)國(guó)家自然地圖集[M].2版.北京:中國(guó)地圖出版社,1999.

    [29]JUAN PG,JOSEM P,INGRID C B.Land use impacts on the normalized difference vegetation index in temperate Argentina[J].Ecological Applications,2003,13(3):616-628.

    [30]李曉文,胡遠(yuǎn)滿,肖篤寧.景觀生態(tài)學(xué)與生物多樣性保護(hù)[J].生態(tài)學(xué)報(bào),1999,19(3):399-407.

    [31]傅伯杰.景觀多樣性分析及其制圖研究[J].生態(tài)學(xué)報(bào),1995,15(4):345-349.

    [32]于秀林,任雪松.多元統(tǒng)計(jì)分析[M].北京:中國(guó)統(tǒng)計(jì)出版社,1999.

    [33]黃方,劉湘南,張養(yǎng)貞.GIS支持下的吉林省西部生態(tài)環(huán)境脆弱態(tài)勢(shì)評(píng)價(jià)研究[J].地理科學(xué),2003,23(1):95-100.

    [34]程建權(quán).城市系統(tǒng)工程[M].武漢:武漢大學(xué)出版社,2001: 122-135.

    [35]李柏年.模糊數(shù)學(xué)及其應(yīng)用[M].合肥:合肥工業(yè)大學(xué)出版社,2007.

    Integrated Assessm ent o f Eco-Environm enta l Vu lnerability in Yangtze River De lta Based on RS and G IS

    XU Qing-yong1,2,HUANG Mei2,LU Pei-ling1,LIRen-qiang2

    1.Beijing Forestry University,Beijing 100083,China

    2.Key Laboratory of Ecosystem Network Observation and Modeling,Institute of Geographic Sciences and Natural Resources Research,Beijing 100101,China

    X826

    A

    1001-6929(2011)01-0058-08

    2010-05-28

    2010-09-10

    國(guó)家“十一五”科技支撐計(jì)劃項(xiàng)目(2008BAK50B01,2007BAC03A11);國(guó)家自然科學(xué)基金項(xiàng)目(30590384)

    徐 慶 勇 (1975 -),女,黑 龍 江 伊 春 人,xuqingyong1234@yahoo.cn.

    *責(zé)任作者,黃玫(1968-),女,貴州興仁人,副研究員,博士,主要從事氣候變化與生態(tài)系統(tǒng)研究,huangm@igsnrr.ac.cn

    猜你喜歡
    長(zhǎng)江三角洲旱澇脆弱性
    基于PSR-TOPSIS模型的長(zhǎng)三角地區(qū)水資源承載力評(píng)價(jià)研究
    四川省1960-2017年旱澇時(shí)空變化特征分析
    煤礦電網(wǎng)脆弱性評(píng)估
    電子制作(2017年10期)2017-04-18 07:23:09
    長(zhǎng)江三角洲各城市主要經(jīng)濟(jì)指標(biāo)(2016年1-9月)
    殺毒軟件中指令虛擬機(jī)的脆弱性分析
    “一帶一路”戰(zhàn)略下的長(zhǎng)江三角洲地區(qū)對(duì)外貿(mào)易轉(zhuǎn)型及發(fā)展思路
    氣候變化背景下西北干旱區(qū)旱澇的變化規(guī)律
    基于攻擊圖的工控系統(tǒng)脆弱性量化方法
    北方農(nóng)牧交錯(cuò)帶氣候變化與旱澇響應(yīng)特征
    寧夏旱澇災(zāi)害對(duì)氣候變化的響應(yīng)
    18美女黄网站色大片免费观看| 日韩中文字幕欧美一区二区| 色综合婷婷激情| 韩国av一区二区三区四区| 欧美日韩av久久| 久久精品亚洲熟妇少妇任你| 黄网站色视频无遮挡免费观看| 变态另类成人亚洲欧美熟女 | 国产1区2区3区精品| 91麻豆av在线| 热99re8久久精品国产| 日本wwww免费看| e午夜精品久久久久久久| 热99re8久久精品国产| 欧美丝袜亚洲另类 | 欧美在线黄色| 亚洲,欧美精品.| av有码第一页| 免费高清在线观看日韩| bbb黄色大片| 成人18禁高潮啪啪吃奶动态图| 久久人人爽av亚洲精品天堂| 高清毛片免费观看视频网站 | 国产视频一区二区在线看| 无限看片的www在线观看| 50天的宝宝边吃奶边哭怎么回事| 人妻丰满熟妇av一区二区三区| 国产成人精品久久二区二区免费| 欧美 亚洲 国产 日韩一| 9热在线视频观看99| 国产精品99久久99久久久不卡| 动漫黄色视频在线观看| 国产欧美日韩一区二区三区在线| 日韩欧美国产一区二区入口| 亚洲国产精品999在线| 日本 av在线| 亚洲精品美女久久久久99蜜臀| 国产成人精品在线电影| 三级毛片av免费| 国产1区2区3区精品| ponron亚洲| 欧美日本中文国产一区发布| 三级毛片av免费| 亚洲中文日韩欧美视频| 亚洲成人免费av在线播放| 亚洲人成77777在线视频| 99热国产这里只有精品6| 国产99久久九九免费精品| 国产无遮挡羞羞视频在线观看| 他把我摸到了高潮在线观看| 在线观看免费日韩欧美大片| 久久影院123| 成人特级黄色片久久久久久久| 国产主播在线观看一区二区| 99久久久亚洲精品蜜臀av| 欧美另类亚洲清纯唯美| 欧美另类亚洲清纯唯美| 可以在线观看毛片的网站| 欧美人与性动交α欧美软件| 欧美一区二区精品小视频在线| 天天躁夜夜躁狠狠躁躁| 免费少妇av软件| 亚洲激情在线av| 美女高潮到喷水免费观看| 亚洲av第一区精品v没综合| 成人亚洲精品一区在线观看| 久久久久久久精品吃奶| 日本免费a在线| 高清毛片免费观看视频网站 | 日韩一卡2卡3卡4卡2021年| 国产高清videossex| 一级毛片精品| 真人做人爱边吃奶动态| 精品国产乱子伦一区二区三区| 高清欧美精品videossex| 天堂俺去俺来也www色官网| 涩涩av久久男人的天堂| 1024香蕉在线观看| 少妇裸体淫交视频免费看高清 | 久久人妻熟女aⅴ| 最新在线观看一区二区三区| 9色porny在线观看| 国产免费男女视频| 老司机靠b影院| 天天影视国产精品| www.999成人在线观看| 午夜福利在线免费观看网站| 视频在线观看一区二区三区| 99在线视频只有这里精品首页| 国产麻豆69| 黄色丝袜av网址大全| 精品高清国产在线一区| 黄色毛片三级朝国网站| 亚洲自偷自拍图片 自拍| 在线观看66精品国产| www.自偷自拍.com| 精品人妻在线不人妻| 亚洲精品美女久久av网站| 男女之事视频高清在线观看| 波多野结衣一区麻豆| 女人被狂操c到高潮| 淫妇啪啪啪对白视频| 欧美日韩瑟瑟在线播放| 不卡一级毛片| 男人舔女人下体高潮全视频| 午夜福利免费观看在线| a在线观看视频网站| 一区在线观看完整版| 丰满迷人的少妇在线观看| 悠悠久久av| 成人免费观看视频高清| 欧美一级毛片孕妇| 老司机深夜福利视频在线观看| 我的亚洲天堂| 亚洲熟妇熟女久久| 99热只有精品国产| 亚洲人成77777在线视频| 99国产精品一区二区蜜桃av| 欧美乱妇无乱码| 狂野欧美激情性xxxx| 9色porny在线观看| 久久精品亚洲熟妇少妇任你| 夜夜看夜夜爽夜夜摸 | 男人舔女人下体高潮全视频| 久久久国产精品麻豆| 这个男人来自地球电影免费观看| 一级作爱视频免费观看| 女人爽到高潮嗷嗷叫在线视频| 久久国产精品男人的天堂亚洲| 国产精品久久久久久人妻精品电影| 亚洲欧洲精品一区二区精品久久久| 国产无遮挡羞羞视频在线观看| 久久久久久大精品| 久久久久精品国产欧美久久久| 一二三四在线观看免费中文在| 岛国在线观看网站| 久久久久国产一级毛片高清牌| 国产片内射在线| 不卡av一区二区三区| 人人妻,人人澡人人爽秒播| 久久中文看片网| 长腿黑丝高跟| 日韩视频一区二区在线观看| 国产精品爽爽va在线观看网站 | 午夜免费鲁丝| 久久国产精品影院| 巨乳人妻的诱惑在线观看| 淫秽高清视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产aⅴ精品一区二区三区波| 女性被躁到高潮视频| 99国产极品粉嫩在线观看| 国产三级黄色录像| 午夜精品久久久久久毛片777| xxx96com| 国产精品1区2区在线观看.| 一区福利在线观看| 久久精品亚洲av国产电影网| 9191精品国产免费久久| 美女高潮喷水抽搐中文字幕| 婷婷丁香在线五月| 午夜视频精品福利| 午夜亚洲福利在线播放| 夜夜躁狠狠躁天天躁| 国产深夜福利视频在线观看| www.www免费av| 久久 成人 亚洲| 国产真人三级小视频在线观看| 嫁个100分男人电影在线观看| 可以在线观看毛片的网站| 国产成人系列免费观看| 69精品国产乱码久久久| 一边摸一边抽搐一进一出视频| 深夜精品福利| 精品一区二区三区四区五区乱码| 欧美日韩精品网址| 五月开心婷婷网| 色综合欧美亚洲国产小说| 亚洲精品在线美女| 免费搜索国产男女视频| 夫妻午夜视频| 久久久久精品国产欧美久久久| 一个人观看的视频www高清免费观看 | 亚洲熟女毛片儿| 欧美日本亚洲视频在线播放| 一级毛片高清免费大全| 一进一出抽搐动态| 可以在线观看毛片的网站| 97碰自拍视频| 亚洲国产精品一区二区三区在线| 美女扒开内裤让男人捅视频| 黑人巨大精品欧美一区二区mp4| 久久精品国产清高在天天线| 日韩大码丰满熟妇| 精品久久蜜臀av无| 好看av亚洲va欧美ⅴa在| 亚洲全国av大片| 国产一区二区在线av高清观看| 青草久久国产| 久久人人爽av亚洲精品天堂| 在线观看免费视频日本深夜| 婷婷六月久久综合丁香| 黄色视频,在线免费观看| 国产麻豆69| 亚洲午夜理论影院| 国产深夜福利视频在线观看| 亚洲中文av在线| 日本一区二区免费在线视频| 窝窝影院91人妻| 一进一出抽搐gif免费好疼 | 成人国语在线视频| 精品高清国产在线一区| 韩国精品一区二区三区| 露出奶头的视频| av片东京热男人的天堂| 国产免费男女视频| 亚洲国产精品sss在线观看 | xxx96com| 99久久久亚洲精品蜜臀av| 欧美黑人欧美精品刺激| 一级片免费观看大全| 亚洲第一欧美日韩一区二区三区| 麻豆久久精品国产亚洲av | 在线观看66精品国产| videosex国产| 亚洲精品久久午夜乱码| 精品一区二区三卡| 免费不卡黄色视频| 成人av一区二区三区在线看| 久久久久国内视频| a级毛片黄视频| 欧美日韩中文字幕国产精品一区二区三区 | 国产极品粉嫩免费观看在线| 在线观看www视频免费| 天堂中文最新版在线下载| 午夜免费鲁丝| 1024香蕉在线观看| 日韩欧美在线二视频| 99久久人妻综合| 男人的好看免费观看在线视频 | 香蕉国产在线看| 不卡av一区二区三区| 啦啦啦 在线观看视频| svipshipincom国产片| 亚洲av成人不卡在线观看播放网| 午夜福利在线观看吧| 欧美乱码精品一区二区三区| 国产一区二区激情短视频| 丰满人妻熟妇乱又伦精品不卡| av天堂在线播放| 久久精品亚洲av国产电影网| 亚洲欧美日韩另类电影网站| 国产成人精品久久二区二区91| 伦理电影免费视频| 亚洲精品成人av观看孕妇| 久久精品国产亚洲av香蕉五月| 久久久久久久精品吃奶| 欧美中文日本在线观看视频| 欧美不卡视频在线免费观看 | 亚洲黑人精品在线| 男女床上黄色一级片免费看| 国产精品偷伦视频观看了| 国产亚洲欧美精品永久| 神马国产精品三级电影在线观看 | 麻豆成人av在线观看| 大陆偷拍与自拍| 日韩大尺度精品在线看网址 | 丝袜人妻中文字幕| 亚洲av五月六月丁香网| а√天堂www在线а√下载| 黄色毛片三级朝国网站| 亚洲av美国av| 搡老乐熟女国产| 精品一区二区三卡| 精品欧美一区二区三区在线| 亚洲美女黄片视频| 97超级碰碰碰精品色视频在线观看| 日韩中文字幕欧美一区二区| 母亲3免费完整高清在线观看| 搡老熟女国产l中国老女人| 91精品三级在线观看| svipshipincom国产片| 天天躁狠狠躁夜夜躁狠狠躁| 日韩一卡2卡3卡4卡2021年| 丝袜在线中文字幕| 成人黄色视频免费在线看| 在线观看66精品国产| 18禁美女被吸乳视频| 精品高清国产在线一区| 女性生殖器流出的白浆| 欧美久久黑人一区二区| 精品乱码久久久久久99久播| 中文字幕最新亚洲高清| 久久国产精品男人的天堂亚洲| 老司机午夜十八禁免费视频| 好男人电影高清在线观看| 成人手机av| 国产乱人伦免费视频| 久久人妻熟女aⅴ| 亚洲av成人一区二区三| 亚洲av电影在线进入| 亚洲 欧美一区二区三区| 欧美激情 高清一区二区三区| 欧美中文综合在线视频| 精品一品国产午夜福利视频| 国产成人啪精品午夜网站| 无遮挡黄片免费观看| 视频区图区小说| 国产99白浆流出| 久久精品成人免费网站| 又紧又爽又黄一区二区| 亚洲精华国产精华精| 亚洲va日本ⅴa欧美va伊人久久| 露出奶头的视频| 中文字幕av电影在线播放| 日韩欧美三级三区| 久久人人爽av亚洲精品天堂| 欧美黑人欧美精品刺激| 91大片在线观看| 黄色女人牲交| 亚洲精品国产色婷婷电影| 久久狼人影院| 中出人妻视频一区二区| 欧美av亚洲av综合av国产av| 亚洲第一欧美日韩一区二区三区| 亚洲第一av免费看| 欧美成人免费av一区二区三区| av天堂在线播放| 国产成人av教育| 黄色 视频免费看| 久久九九热精品免费| 久久中文字幕一级| 窝窝影院91人妻| 一二三四社区在线视频社区8| xxx96com| 一区二区日韩欧美中文字幕| 国内久久婷婷六月综合欲色啪| 在线观看一区二区三区| www.精华液| 午夜a级毛片| 精品日产1卡2卡| 巨乳人妻的诱惑在线观看| 欧美黄色淫秽网站| 国产成人系列免费观看| 高清av免费在线| 美女大奶头视频| 欧美中文日本在线观看视频| 欧美在线黄色| 国产精品 欧美亚洲| 亚洲片人在线观看| 亚洲欧洲精品一区二区精品久久久| 亚洲欧美日韩无卡精品| 欧美成人免费av一区二区三区| 国产精品久久久av美女十八| 91麻豆av在线| 日韩 欧美 亚洲 中文字幕| www国产在线视频色| 国产亚洲精品第一综合不卡| 久久久久国产一级毛片高清牌| 亚洲国产精品999在线| 久久精品国产99精品国产亚洲性色 | 久久香蕉精品热| www.www免费av| 国产亚洲av高清不卡| 18美女黄网站色大片免费观看| 亚洲专区国产一区二区| 在线观看免费日韩欧美大片| 亚洲成人久久性| 中文字幕人妻熟女乱码| 97碰自拍视频| 久久久国产成人精品二区 | 9色porny在线观看| 久久久水蜜桃国产精品网| 高潮久久久久久久久久久不卡| 亚洲av成人av| 热re99久久国产66热| 国产精品日韩av在线免费观看 | 亚洲精品在线观看二区| 欧美中文日本在线观看视频| avwww免费| 久久伊人香网站| 国产欧美日韩一区二区三| 久久午夜综合久久蜜桃| 国产成年人精品一区二区 | 久久人妻熟女aⅴ| 色婷婷久久久亚洲欧美| 一个人观看的视频www高清免费观看 | 亚洲一区二区三区不卡视频| 亚洲精品中文字幕一二三四区| 女性被躁到高潮视频| 国产精品久久久久久人妻精品电影| 老司机在亚洲福利影院| 在线观看免费日韩欧美大片| 日韩免费av在线播放| 日本wwww免费看| 新久久久久国产一级毛片| 性欧美人与动物交配| 国产成人系列免费观看| 黑人猛操日本美女一级片| 免费av毛片视频| xxx96com| 欧美日韩中文字幕国产精品一区二区三区 | 精品第一国产精品| 亚洲精品美女久久av网站| 一级毛片高清免费大全| 免费女性裸体啪啪无遮挡网站| 老司机亚洲免费影院| 夜夜爽天天搞| 精品国产亚洲在线| 午夜免费激情av| 国产一区二区三区视频了| 在线观看免费午夜福利视频| 免费看a级黄色片| 欧美一级毛片孕妇| 精品人妻在线不人妻| www.熟女人妻精品国产| 亚洲片人在线观看| 国产在线精品亚洲第一网站| 国产精品一区二区三区四区久久 | 极品人妻少妇av视频| 最好的美女福利视频网| 国产激情久久老熟女| 日韩欧美一区二区三区在线观看| 亚洲情色 制服丝袜| 一二三四社区在线视频社区8| 久久久久精品国产欧美久久久| 在线视频色国产色| 男女下面插进去视频免费观看| 免费搜索国产男女视频| 精品高清国产在线一区| 不卡av一区二区三区| 黄色视频,在线免费观看| 久久人妻av系列| 美女大奶头视频| 热99re8久久精品国产| 久久香蕉精品热| 狠狠狠狠99中文字幕| 国产人伦9x9x在线观看| av在线播放免费不卡| 国产成人av教育| 亚洲色图av天堂| 国产精品爽爽va在线观看网站 | 美女国产高潮福利片在线看| 成人亚洲精品一区在线观看| 欧洲精品卡2卡3卡4卡5卡区| 新久久久久国产一级毛片| 在线观看免费日韩欧美大片| 国产激情欧美一区二区| 婷婷六月久久综合丁香| bbb黄色大片| 精品国产超薄肉色丝袜足j| 母亲3免费完整高清在线观看| 精品卡一卡二卡四卡免费| 久久精品91蜜桃| 女人被躁到高潮嗷嗷叫费观| 在线观看午夜福利视频| av电影中文网址| 欧美激情极品国产一区二区三区| 国产av一区在线观看免费| 久久国产乱子伦精品免费另类| 操出白浆在线播放| 超碰97精品在线观看| 日韩大码丰满熟妇| 亚洲av五月六月丁香网| 亚洲国产精品一区二区三区在线| 欧美不卡视频在线免费观看 | 日韩精品青青久久久久久| 老司机午夜福利在线观看视频| 少妇被粗大的猛进出69影院| 日日摸夜夜添夜夜添小说| 免费观看人在逋| 麻豆成人av在线观看| 一级作爱视频免费观看| 在线观看一区二区三区激情| 最新美女视频免费是黄的| 亚洲av日韩精品久久久久久密| 不卡av一区二区三区| 99riav亚洲国产免费| 大型黄色视频在线免费观看| 国产亚洲欧美精品永久| 国产欧美日韩精品亚洲av| 一本综合久久免费| 视频在线观看一区二区三区| 伊人久久大香线蕉亚洲五| 日韩一卡2卡3卡4卡2021年| 亚洲精品在线观看二区| 精品久久久久久电影网| 国产精品一区二区在线不卡| 久久精品91蜜桃| 免费在线观看影片大全网站| 嫁个100分男人电影在线观看| 变态另类成人亚洲欧美熟女 | 最近最新中文字幕大全免费视频| 亚洲色图av天堂| 亚洲aⅴ乱码一区二区在线播放 | 嫩草影视91久久| 视频区图区小说| 亚洲美女黄片视频| 老司机午夜福利在线观看视频| 成人国产一区最新在线观看| 国产欧美日韩综合在线一区二区| 欧美亚洲日本最大视频资源| 制服诱惑二区| 色哟哟哟哟哟哟| tocl精华| 亚洲精品国产区一区二| 久久午夜亚洲精品久久| 精品一区二区三区视频在线观看免费 | 超色免费av| 国产亚洲精品久久久久久毛片| 日韩大码丰满熟妇| 黄网站色视频无遮挡免费观看| 91国产中文字幕| 国产成人欧美| 成人黄色视频免费在线看| 色哟哟哟哟哟哟| av福利片在线| 97碰自拍视频| 99久久精品国产亚洲精品| 久久久久亚洲av毛片大全| av片东京热男人的天堂| 亚洲欧美日韩无卡精品| 婷婷丁香在线五月| 精品国产一区二区三区四区第35| 成人亚洲精品av一区二区 | 两性午夜刺激爽爽歪歪视频在线观看 | 欧美激情久久久久久爽电影 | 成人精品一区二区免费| 亚洲七黄色美女视频| 色哟哟哟哟哟哟| 精品欧美一区二区三区在线| 亚洲精品美女久久av网站| 欧美日韩瑟瑟在线播放| 99精品久久久久人妻精品| 亚洲国产看品久久| 成人国产一区最新在线观看| 欧美中文日本在线观看视频| 99国产精品免费福利视频| 欧美性长视频在线观看| 一区在线观看完整版| 国产真人三级小视频在线观看| 久久久久久久久免费视频了| а√天堂www在线а√下载| 丝袜在线中文字幕| 高清av免费在线| 国产三级黄色录像| 黄色成人免费大全| 99久久精品国产亚洲精品| 久久中文看片网| 无限看片的www在线观看| 如日韩欧美国产精品一区二区三区| 久久国产乱子伦精品免费另类| 久久人妻av系列| 国内毛片毛片毛片毛片毛片| 亚洲欧美激情在线| 两人在一起打扑克的视频| a级片在线免费高清观看视频| 在线观看舔阴道视频| 高清毛片免费观看视频网站 | 国产精品美女特级片免费视频播放器 | 99在线人妻在线中文字幕| 国产精品乱码一区二三区的特点 | 看黄色毛片网站| 色精品久久人妻99蜜桃| 男女午夜视频在线观看| 夜夜躁狠狠躁天天躁| 国产成人精品无人区| 国产欧美日韩一区二区三区在线| 日本三级黄在线观看| 黑丝袜美女国产一区| 热re99久久精品国产66热6| 精品欧美一区二区三区在线| 欧美激情极品国产一区二区三区| 黄色成人免费大全| 国产97色在线日韩免费| 老汉色av国产亚洲站长工具| 9热在线视频观看99| 久久这里只有精品19| 在线免费观看的www视频| 在线十欧美十亚洲十日本专区| 久99久视频精品免费| 一进一出好大好爽视频| 国产免费男女视频| 久99久视频精品免费| 一进一出好大好爽视频| 国产精品久久久久成人av| 亚洲人成电影免费在线| 精品国产美女av久久久久小说| 中文字幕精品免费在线观看视频| 久久亚洲真实| 一区二区三区激情视频| a级毛片在线看网站| e午夜精品久久久久久久| 国产一区在线观看成人免费| 日韩欧美在线二视频| 宅男免费午夜| 久久 成人 亚洲| 久久中文字幕一级| 精品福利观看| 午夜精品在线福利| 91字幕亚洲| 欧美一级毛片孕妇| 在线免费观看的www视频| 男男h啪啪无遮挡| 欧美激情高清一区二区三区| 欧美精品亚洲一区二区| 麻豆国产av国片精品| 国产三级在线视频| 成人18禁在线播放| 99在线人妻在线中文字幕| 一本大道久久a久久精品| 亚洲一区二区三区欧美精品| 午夜精品国产一区二区电影| 老司机午夜福利在线观看视频| 宅男免费午夜|