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

    內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量空間特征及其驅(qū)動(dòng)力

    2016-11-28 08:55:35歐陽志云王莉雁饒恩明
    生態(tài)學(xué)報(bào) 2016年19期
    關(guān)鍵詞:覆蓋度生物量內(nèi)蒙古

    肖 洋,歐陽志云,王莉雁,饒恩明,江 凌,張 路

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

    ?

    內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量空間特征及其驅(qū)動(dòng)力

    肖 洋1,2,歐陽志云1,*,王莉雁1,2,饒恩明1,2,江 凌1,2,張 路1,2

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

    植被作為生態(tài)系統(tǒng)的重要組成部分,聯(lián)結(jié)著大氣、水分和土壤等自然過程,其變化將直接影響該區(qū)域氣候水文和土壤等狀況,是區(qū)域生態(tài)系統(tǒng)質(zhì)量變化的重要指示器。植被狀況的好壞,主要通過生物量和植被覆蓋度因子來表示。內(nèi)蒙古自治區(qū)是我國北方生態(tài)環(huán)境問題十分嚴(yán)重的省份,弄清當(dāng)前區(qū)域生態(tài)系統(tǒng)質(zhì)量狀況與變化及其近10年來變化的驅(qū)動(dòng)因素,對分析與制定區(qū)域生態(tài)環(huán)境保護(hù)決策具有十分重要的意義?;?000—2010年生物量和植被覆蓋度,并結(jié)合地區(qū)植被區(qū)劃數(shù)據(jù),對內(nèi)蒙古植被生態(tài)系統(tǒng)質(zhì)量狀況進(jìn)行分析,并評估其與氣候(降水、溫度),人類活動(dòng)(交通密度、農(nóng)業(yè)發(fā)展、生態(tài)恢復(fù)工程)的相關(guān)關(guān)系,在此基礎(chǔ)上探討了氣候和人類活動(dòng)對近年來內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量變化的影響。結(jié)果表明:(1) 內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量狀況整體偏低,其中森林生態(tài)系統(tǒng)平均質(zhì)量最高,灌叢、草原生態(tài)系統(tǒng)次之??臻g分布呈明顯的經(jīng)度地帶性,由東向西,質(zhì)量逐漸降低。2000—2010年內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量總體上呈現(xiàn)緩慢增長趨勢,但局部地區(qū)生態(tài)系統(tǒng)質(zhì)量仍存在惡化,其中在107°E以東的草原和森林區(qū)域,生態(tài)系統(tǒng)質(zhì)量變化十分劇烈。(2) 近10年來內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量的變化與氣候和人類活動(dòng)的關(guān)系非常密切,其與降水、GDP1、化肥施用量、天保工程和退耕還草工程呈現(xiàn)明顯的正相關(guān)。而與溫度、道路密度和京津風(fēng)沙治理工程呈現(xiàn)明顯的負(fù)相關(guān)。其中,生態(tài)保護(hù)工程實(shí)施區(qū)域內(nèi)和區(qū)域外的相關(guān)性存在顯著的差異性。隨著內(nèi)蒙古社會(huì)經(jīng)濟(jì)的快速發(fā)展,人類活動(dòng)對生態(tài)系統(tǒng)質(zhì)量的影響逐漸加強(qiáng),但降水仍是該地區(qū)生態(tài)系統(tǒng)質(zhì)量的主要影響因子。(3) 在內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量變化典型區(qū)域內(nèi),質(zhì)量的增長主要是由于降水的增加、溫度的降低、農(nóng)業(yè)的發(fā)展、退耕還草工程的作用和交通發(fā)展的放緩。質(zhì)量的降低則是因?yàn)榻邓臏p少、溫度的增加、農(nóng)業(yè)發(fā)展緩慢和交通發(fā)展的加快所致。

    生物量;植被覆蓋度;生態(tài)系統(tǒng)質(zhì)量;氣候;人類活動(dòng)

    植被作為生態(tài)系統(tǒng)的重要組成部分,聯(lián)結(jié)著大氣、水分和土壤等自然過程,其變化將直接影響該區(qū)域氣候、水文和土壤等狀況,對區(qū)域能量循環(huán)及物質(zhì)的生物化學(xué)循環(huán)具有重要的影響,是區(qū)域生態(tài)系統(tǒng)質(zhì)量變化的重要指示器。目前反映植被生長狀況及生長活力的因子較多,有植被歸一化因子NDVI、葉面積指數(shù)LAI、凈第一生產(chǎn)力NPP、植被覆蓋度VC、生物量Biomass等。由于各因子所具有的生態(tài)學(xué)意義不同,所以其適用范圍也存在差異。其中生物量作為生態(tài)系統(tǒng)中植物有機(jī)物總量,不僅是整個(gè)生態(tài)系統(tǒng)運(yùn)行的能量基礎(chǔ)和營養(yǎng)物質(zhì)來源,還能反映生態(tài)系統(tǒng)的生態(tài)演替和人為干擾過程[1]。其中,森林地上生物量不僅能充分體現(xiàn)森林生態(tài)系統(tǒng)的質(zhì)量,衡量森林生產(chǎn)力,也是評估森林碳收支的重要參數(shù)[2- 4]。生物量在生態(tài)系統(tǒng)結(jié)構(gòu)功能、碳儲(chǔ)量和碳動(dòng)態(tài)方面的有著廣泛應(yīng)用[5- 8]。植被覆蓋度是指植被(包括葉、莖、枝)在地面的垂直投影面積占統(tǒng)計(jì)區(qū)總面積的百分比[9-10]。它是全球及區(qū)域氣候模型、水土流失監(jiān)測、土地沙漠化評價(jià)和分布式水文模型的重要輸入?yún)?shù),是描述生態(tài)系統(tǒng)的重要基礎(chǔ)數(shù)據(jù)[11- 13]。在分析生態(tài)系統(tǒng)環(huán)境狀態(tài),揭示其變化趨勢方面,植被覆蓋度應(yīng)用廣泛[14-15]。植被因子作為表征生態(tài)系統(tǒng)變化的綜合指示器,對評價(jià)陸地生態(tài)系統(tǒng)質(zhì)量、調(diào)節(jié)生態(tài)過程具有重要的理論和實(shí)際意義。

    內(nèi)蒙古位處中國北部邊疆,是我國北部重要的生態(tài)屏障,同時(shí)也是我國重要的農(nóng)牧業(yè)生產(chǎn)基地,人類活動(dòng)強(qiáng)度大,且大部分植被地處干旱、半干旱農(nóng)牧交錯(cuò)帶地區(qū),生態(tài)環(huán)境脆弱,是全球氣候變化最為敏感的區(qū)域之一[16]。特別是,近年來內(nèi)蒙古發(fā)生著巨大的變化,經(jīng)濟(jì)持續(xù)快速增長、城鎮(zhèn)化進(jìn)程加快、資源開發(fā)力度的增大、生態(tài)環(huán)境受到一定的沖擊,同時(shí)生態(tài)保護(hù)與建設(shè)工程 (天然林保護(hù)、退耕還林還草) 又在一定程度上減緩著人類活動(dòng)造成的壓力。這些變化必然會(huì)影響內(nèi)蒙古的生態(tài)系統(tǒng)質(zhì)量狀況,進(jìn)而影響了地區(qū)經(jīng)濟(jì)和生態(tài)環(huán)境的可持續(xù)發(fā)展。因此,全面了解內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量狀況和10年來的變化趨勢,客觀認(rèn)識(shí)生態(tài)系統(tǒng)結(jié)構(gòu)與功能,確定重點(diǎn)保護(hù)區(qū)域以及目前存在的潛在威脅,探討生態(tài)系統(tǒng)質(zhì)量變化的驅(qū)動(dòng)作用,實(shí)現(xiàn)內(nèi)蒙古的可持續(xù)發(fā)展具有重要現(xiàn)實(shí)意義。本研究結(jié)合內(nèi)蒙古的生態(tài)環(huán)境特征,選擇能有效反映森林生態(tài)系統(tǒng)質(zhì)量的生物量因子和表征草地生態(tài)系統(tǒng)質(zhì)量的植被覆蓋度因子來評估地區(qū)生態(tài)系統(tǒng)質(zhì)量狀況的好壞。

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

    1.1 研究區(qū)概況

    圖1 研究區(qū)位置Fig.1 Location of the study area

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

    地上生物量數(shù)據(jù)采用植被指數(shù)-生物量法(對于森林/草地生態(tài)系統(tǒng))和累積NPP法(對于草地/農(nóng)田生態(tài)系統(tǒng))對植被生物量進(jìn)行估算[20]。其中植被指數(shù)-生物量法是采用實(shí)地測量的植被生物量的數(shù)據(jù)和遙感數(shù)據(jù)建立經(jīng)驗(yàn)統(tǒng)計(jì)模型,然后在遙感數(shù)據(jù)的基礎(chǔ)上反演得到區(qū)域范圍內(nèi)植被生物量;而累積NPP 法是通過草地或農(nóng)田的生長期(開始生長時(shí)間與結(jié)束生長時(shí)間)的確定,對生長期內(nèi)的NPP進(jìn)行累加以計(jì)算地上生物量,NPP算法選用CASA模型[21]。植被覆蓋度數(shù)據(jù)是基于像元二分模型[22-24]通過Modis影像反演得到。像元二分模型是一種簡單實(shí)用的遙感估算模型,它假設(shè)一個(gè)像元的地表由有植被覆蓋部分與無植被覆蓋部分地表組成,而遙感傳感器觀測到的光譜信息也由這2個(gè)組分因子線性加權(quán)合成,各因子的權(quán)重是各自的面積在像元中所占的比率,其中植被覆蓋度可以看作是植被的權(quán)重。生態(tài)系統(tǒng)類型圖主要基于LandsatTM數(shù)據(jù)采用的面向?qū)ο蟮姆诸惣夹g(shù),引入非影像光譜信息強(qiáng)化目標(biāo)的識(shí)別能力;實(shí)現(xiàn)快速、高效的分類技術(shù)運(yùn)作。數(shù)字高程模型DEM,空間分辨率為90 m, 來源于國際科學(xué)數(shù)據(jù)服務(wù)平臺(tái)。降水與溫度數(shù)據(jù),基于局部薄盤樣條函數(shù)插值理論,來源于中國國家計(jì)量信息中心/中國氣象局(NMIC/CMA)。詳細(xì)數(shù)據(jù)信息,見表1。

    表1 主要數(shù)據(jù)來源

    2 研究方法

    2.1 生態(tài)系統(tǒng)質(zhì)量

    本文選取地上生物量和植被覆蓋度因子,結(jié)合植被區(qū)劃數(shù)據(jù),以期反映內(nèi)蒙古森林和草地生態(tài)系統(tǒng)的質(zhì)量狀況。基于內(nèi)蒙古2000與2010年生態(tài)系統(tǒng)分類圖,結(jié)合原始森林和草地調(diào)查樣方,選取各植被地帶內(nèi)的多個(gè)原始森林樣方的生物量的均值和草地樣方的覆蓋度均值來代表該范圍內(nèi)頂級群落的生物量和覆蓋度。其中,頂級群落的樣方是經(jīng)過設(shè)定嚴(yán)格的條件,如生物量高、凋落物厚、樹齡高等,篩選得到。以各個(gè)地帶的頂級群落的生物量和覆蓋度為分母,生物量和覆蓋度為分子,得到各個(gè)植被地帶的生態(tài)系統(tǒng)質(zhì)量。公式如下:

    式中,RBDi為森林或草地生態(tài)系統(tǒng)i像元的相對生態(tài)系統(tǒng)質(zhì)量密度;Bi為森林生態(tài)系統(tǒng)i像元的生物量,通過遙感反演獲得;Ci為草地生態(tài)系統(tǒng)i像元的覆蓋度,通過遙感反演獲得;CCB為森林生態(tài)系統(tǒng)頂級群落像元的生物量;CCC為草地生態(tài)系統(tǒng)頂級群落像元的覆蓋度;頂級群落樣本均運(yùn)用生態(tài)系統(tǒng)長期定位觀測數(shù)據(jù),或樣地調(diào)查數(shù)據(jù)。

    2.2 數(shù)據(jù)分析

    2.2.1 變化趨勢分析

    為研究10a來我國生態(tài)系統(tǒng)質(zhì)量變化對區(qū)域氣候、人類活動(dòng)變化的響應(yīng)以及這種變化在時(shí)空上的差異,選用一元線性回歸分析來處理多年氣候和人類活動(dòng)統(tǒng)計(jì)數(shù)據(jù)[25-26],得到其變化趨勢,計(jì)算公式為:

    墓穴被封上了。老巴被扶上去敬了香。然后他點(diǎn)燃了鞭炮。響聲瞬間震撼整個(gè)墓區(qū),碎紙隨風(fēng)揚(yáng)起,陽光仿佛也在抖動(dòng)。

    y=a+b×t+ε

    式中,y為因變量,t為年份,a和b為擬合參數(shù) (b為坡度斜率 代表趨勢特征;a為截距),ε為殘差。斜率為負(fù)值的區(qū)域,呈減小趨勢;斜率為正值地區(qū),呈增加趨勢;斜率絕對值越大,變化的幅度越大,反之則變化的幅度越小。

    2.2.2 變化相關(guān)性分析

    為了探索生態(tài)系統(tǒng)質(zhì)量與氣候和人類活動(dòng)之間的關(guān)系,選用皮爾遜相關(guān)性方法進(jìn)行分析。其中氣候因素主要考慮降水和溫度,人類活動(dòng)主要為放牧壓力(牛羊存欄數(shù)),人口 (農(nóng)村和城鎮(zhèn)人口),國內(nèi)生產(chǎn)總值GDP(農(nóng)業(yè)生產(chǎn)總值GDP1, 工業(yè)生產(chǎn)總值GDP2, 服務(wù)業(yè)生產(chǎn)總值GDP3),交通密度、建設(shè)用地面積、農(nóng)藥化肥施用量、礦區(qū)面積、生態(tài)保護(hù)工程等等。相關(guān)性方分析均以縣域?yàn)榻y(tǒng)計(jì)單元,在SPSS軟件中進(jìn)行。

    3 結(jié)果與分析

    3.1 生態(tài)系統(tǒng)質(zhì)量時(shí)空特性

    內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量整體偏低。按等分方法分類,2010年質(zhì)量為RBD>80的生態(tài)系統(tǒng)占總面積的3.1%,質(zhì)量為60

    森林生態(tài)系統(tǒng)2010年平均質(zhì)量達(dá)到44.83,主要包括位于呼倫貝爾高原的大興安嶺和西遼河平原的山區(qū)。地處歐亞大陸中高緯度地帶的大興安嶺,受人類活動(dòng)影響較少,植被組成以寒溫帶針葉林為主,生態(tài)系統(tǒng)質(zhì)量多為49.11,而處于東北平原與內(nèi)蒙古高原的過渡區(qū)域的山區(qū),是我國北方主要的農(nóng)牧交錯(cuò)區(qū),受人類活動(dòng)影響較大,生態(tài)系統(tǒng)質(zhì)量多為32.26。位于溫帶半干旱氣候帶的草地生態(tài)系統(tǒng),2010年平均質(zhì)量為31.70。其中質(zhì)量高的區(qū)域主要分布在呼倫貝爾草原和錫林郭勒草原。鄂爾多斯高原草地質(zhì)量偏低,為16.83。

    內(nèi)蒙古2000—2010年的生態(tài)系統(tǒng)質(zhì)量總體上呈現(xiàn)增長狀態(tài),質(zhì)量增加的面積約為120900.43 km2。10a間生態(tài)系統(tǒng)質(zhì)量明顯增加的區(qū)域主要分布在賀蘭山、毛烏素沙地西部及渾善達(dá)克沙地北部地區(qū)、錫林郭勒高原西北部、科爾沁沙地東部、呼倫貝爾高原西部地區(qū)和東北部大興安嶺部分地區(qū);其中錫林郭勒附近的生態(tài)系統(tǒng)質(zhì)量增加趨勢最為明顯,增加面積約為31808.01 km2(圖2)。生態(tài)系統(tǒng)質(zhì)量明顯減少的區(qū)域分布相對較集中,除少數(shù)零散分布于科爾沁沙地西部及東北部大興安嶺的部分地區(qū)外,主要沿鄂爾多斯高原-渾善達(dá)克沙地南緣-錫林郭勒高原東南部一線呈條帶狀分布。其中鄂爾多斯和呼和浩特附近的生態(tài)系統(tǒng)質(zhì)量下降趨勢明顯,降低面積約為34922.75 km2(圖2)。

    圖2 內(nèi)蒙古2010年生態(tài)系統(tǒng)質(zhì)量空間分布和內(nèi)蒙古2000—2010年生態(tài)系統(tǒng)質(zhì)量變化Fig.2 Spatial pattern of ecosystem quality in Inner Mongolia in 2010 and variation of ecosystem quality area in Inner Mongolia from 2000—2010

    從內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量的隨經(jīng)度變化規(guī)律可看出 (圖3),內(nèi)蒙古2010年生態(tài)系統(tǒng)質(zhì)量自西向東的呈逐漸增加趨勢,具有明顯的經(jīng)度地帶性(R2= 0.87,P<0.001)。內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量變化方面,100°—106°N在地理空間上基本對應(yīng)荒漠草原生態(tài)系統(tǒng)類型,質(zhì)量變化趨勢較平緩;位于107°—117°N的平原地區(qū),由于氣候變化和人類活動(dòng)的影響,生態(tài)系統(tǒng)質(zhì)量變化劇烈,主要呈降低趨勢;118°N以東地區(qū),生態(tài)系統(tǒng)質(zhì)量主要呈增加趨勢。

    圖3 內(nèi)蒙古2000—2010年生態(tài)系統(tǒng)質(zhì)量經(jīng)向變化規(guī)律Fig.3 Spatial pattern of ecosystem quality and its variation with longitude in Inner Mongolia from 2000—2010

    3.2 生態(tài)系統(tǒng)質(zhì)量變化與氣候和人類活動(dòng)之間的關(guān)系

    由于氣候的作用和人類活動(dòng)的影響,內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量發(fā)生著巨大的變化。本文首先對生態(tài)系統(tǒng)質(zhì)量變化與所有可能的影響因子做相關(guān)性分析,發(fā)現(xiàn)與降水、溫度、農(nóng)業(yè)、道路密度、生態(tài)工程因子相關(guān)性明顯(表2)。生態(tài)系統(tǒng)質(zhì)量與降水、GDP1、化肥施用量、天保工程和退耕還草工程呈現(xiàn)明顯的正相關(guān)(r=0. 543,P<0.01;r=0. 308,P<0.01;r=0. 417,P<0.01;r=0.292,P<0.01;r=0.395,P< 0.01)。然而,與溫度、道路密度和京津風(fēng)沙治理工程呈現(xiàn)明顯的負(fù)相關(guān)(r=-0.381,P<0.01;r=-0.299,P<0.01;r=-0.346,P<0.01)。其中,生態(tài)保護(hù)工程實(shí)施區(qū)域內(nèi)和區(qū)域外的相關(guān)性存在顯著的差異性。例如,在實(shí)施天然林保護(hù)和退耕還草工程區(qū)域內(nèi),生態(tài)系統(tǒng)質(zhì)量變化與降水、溫度、GDP1、道路密度等因子的相關(guān)性不顯著;但在區(qū)域外(未實(shí)施),它們之間的相關(guān)性非常顯著。說明在工程區(qū)域內(nèi),生態(tài)系統(tǒng)質(zhì)量的變化主要與植樹造林、退耕還草、圍封轉(zhuǎn)移和有效管理等密切相關(guān),而受降水、溫度等其他因素的影響較??;而工程區(qū)域外,生態(tài)系統(tǒng)質(zhì)量的變化主要受降水、溫度等因素的影響。然而,京津風(fēng)沙工程區(qū)的狀態(tài)剛好與之相反,區(qū)域內(nèi)主要受到降水、溫度的影響,說明該區(qū)域人類活動(dòng)的影響較小,氣候仍然是主導(dǎo)因素。

    表2 生態(tài)系統(tǒng)質(zhì)量變化與氣候和人類活動(dòng)之間的相關(guān)系

    ** 在 0.01 水平(雙側(cè))上顯著相關(guān);* 在 0.05 水平(雙側(cè))上顯著相關(guān)

    圖4 內(nèi)蒙古2000—2010年生態(tài)系統(tǒng)質(zhì)量變化與影響因子的相關(guān)性 (降水、溫度、GDP1、道路密度、化肥用量)Fig.4 Relation between ecosystem quality and driving factors (preicipitation, temperature, gross domestic product (agricultural product), amount of fertilizer use, road density) from 2000 to 2010 in Inner Mongolia (Each dot represents a county; The line represents the fitted curve)

    圖4(以縣域?yàn)榻y(tǒng)計(jì)單位的散點(diǎn)圖)清楚的說明了這些相關(guān)信息。近10年來內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量的變化與降水的相關(guān)性最大(r=0.543,P<0.01),反映了降水是該地區(qū)生態(tài)系統(tǒng)質(zhì)量變化的主要影響因子。由于內(nèi)蒙古主要的生態(tài)系統(tǒng)類型為草地,同時(shí)草原區(qū)熱量條件相對充足,能夠滿足植被生長的需要,降雨量成為控制植被的生長狀況主要?dú)夂蛞蜃?。張戈麗等[27]研究表明降水是驅(qū)動(dòng)呼倫貝爾草地植被年際變化的主要因素。雖然生態(tài)系統(tǒng)質(zhì)量變化與溫度的相關(guān)關(guān)系不如降水那樣明顯,但總體而言相關(guān)較為顯著。溫度的增加對植被生長有著正負(fù)效應(yīng):正效應(yīng)是延長植被生長季節(jié),提高光合作用效率和水分利用率,從而促進(jìn)植被的生長;而負(fù)效應(yīng)主要在于增加水分消耗從而引發(fā)干旱,不利于植被生長[28]。內(nèi)蒙古地區(qū)屬于干旱、半干旱地區(qū),溫度增加使該地區(qū)缺少水分,從而不利于植被生長。除了氣候變化的影響,人類活動(dòng)對生態(tài)系統(tǒng)存在一定的干擾。在人口稀少的內(nèi)蒙古地區(qū),植被生長狀況主要受自然因素的影響較大,而受人類活動(dòng)的影響較小。但是隨著社會(huì)經(jīng)濟(jì)的快速發(fā)展,人口的增長對資源的需求量的加大,促進(jìn)交通運(yùn)輸?shù)陌l(fā)展,嚴(yán)重制約著該地區(qū)生態(tài)系統(tǒng)質(zhì)量的恢復(fù)與改善。與此同時(shí),近年來國家實(shí)施了天然林保護(hù)、退耕還草、京津風(fēng)沙源治理、禁牧輪牧和圍封轉(zhuǎn)移等一系列生態(tài)恢復(fù)工程,在一定程度上又減少了人類活動(dòng)的負(fù)面影響。此外,農(nóng)業(yè)的發(fā)展、合理的管理和利用草地資源,在一定程度上利于草地質(zhì)量的改善。因此,該區(qū)域的生態(tài)系統(tǒng)質(zhì)量的恢復(fù)與改善,需要重點(diǎn)考慮水分、溫度、農(nóng)業(yè)發(fā)展和生態(tài)工程的影響。

    為了進(jìn)一步了解內(nèi)蒙古近10年來生態(tài)系統(tǒng)質(zhì)量的變化狀態(tài),將生態(tài)系統(tǒng)質(zhì)量變化區(qū)域分為主要增加和降低兩個(gè)典型區(qū)域。增加典型區(qū)域位于毛烏素沙地西部、錫林郭勒高原西北部和科爾沁沙地東部地區(qū)(圖5);相對于增長地區(qū),降低典型區(qū)域分布相對較為集中,主要沿鄂爾多斯高原-渾善達(dá)克沙地南緣-錫林郭勒高原東南部一線呈條帶狀分布(圖5)。統(tǒng)計(jì)分析增加和降低區(qū)域內(nèi)各影響因子狀況(表3),發(fā)現(xiàn)增加區(qū)域內(nèi)生態(tài)系統(tǒng)質(zhì)量改善非常明顯,為5.22 遠(yuǎn)高于總區(qū)域均值0.16。其中降水量的增加非常明顯,遠(yuǎn)高于總區(qū)域均值0.28 mm/a。同時(shí),退耕還草工程主要在該地區(qū)實(shí)施(占69.57%的面積),并且效果明顯。這里的溫度也呈一定的降低趨勢,每年降低0.02 ℃。此外,道路密度的增長速度較慢,低于平均水平。農(nóng)業(yè)的發(fā)展和化肥施用量的增加,在一定程度上利于該區(qū)域植被的恢復(fù)。由于內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量的變化與降水、農(nóng)業(yè)、退耕還草成正相關(guān)性,而與溫度和道路成負(fù)相關(guān)性。所以該區(qū)域生態(tài)系統(tǒng)質(zhì)量的改善主要是由于區(qū)域內(nèi)降水的增加、溫度的降低、農(nóng)業(yè)的發(fā)展、退耕還草工程的作用和交通發(fā)展的放緩。與此同時(shí),降低典型區(qū)域內(nèi)生態(tài)系統(tǒng)質(zhì)量退化非常嚴(yán)重,為-7.43,遠(yuǎn)低于總區(qū)域的均值。其中降水量的減少非常明顯,大約為-0.37 mm/a。這里的溫度也呈現(xiàn)一定的上升趨勢,每年增加0.02 ℃。此外,道路密度的增長速度較快,農(nóng)業(yè)發(fā)展比較緩慢。結(jié)果表明該區(qū)域生態(tài)系統(tǒng)質(zhì)量的下降主要是因?yàn)榻邓臏p少、溫度的增加、農(nóng)業(yè)發(fā)展緩慢和交通發(fā)展的加快所致。雖然京津風(fēng)沙工程主要分布在該地區(qū)(占68.42%的面積),但是治理效果并不明顯,可能是由于降水與溫度的變化導(dǎo)致的植被退化程度要大于生態(tài)恢復(fù)工程的治理效果。

    圖5 內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量主要增加區(qū)域和內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量主要降低區(qū)域Fig.5 The mainly increase region of ecosystem quality in Inner Mongolia and the mainly decrease region of ecosystem quality in Inner Mongolia

    區(qū)域Regions生態(tài)系統(tǒng)質(zhì)量Ecosystemquality降水/(mm/a)Precipitation溫度/(°/a)TemperatureGDP1Grossdomesticproduc/(×108元)道路密度Roaddensity/(km/km2)化肥總量Amountofchemicalfertilizer/(×104t)天保工程N(yùn)aturalforestprotectionproject/%退耕還草工程GreenforGrainProject/%京津風(fēng)沙工程Windandsandstormcontrolproject/%增加區(qū)域ImprovementAreas5.22(+)0.47(+)-0.02(-)10.04(+)0.18(-)2.52(+)30.43(-)69.57(+)26.09(-)降低區(qū)域Degradationareas-7.43(-)-0.37(-)0.02(+)4.92(-)0.28(+)0.29(-)21.05(-)15.79(-)68.42(+)總區(qū)域Wholeregion0.160.280.008.850.201.1841.1344.8232.51

    (+) 相對于總區(qū)域平均值而言,呈增加趨勢; (-) 相對于總區(qū)域平均值而言,呈減少趨勢

    近些年來,國家越來越重視生態(tài)環(huán)境的保護(hù),在全國范圍內(nèi)實(shí)施了一系列的生態(tài)恢復(fù)工程,而內(nèi)蒙古又是這些工程的重點(diǎn)實(shí)施區(qū)域之一。在工程實(shí)施區(qū)域,人類活動(dòng)強(qiáng)度的加大對植被狀況的時(shí)空演變產(chǎn)生了重要影響,同時(shí)導(dǎo)致研究區(qū)植被對氣候因子的變化敏感性降低,這與黃土高原植被變化監(jiān)測結(jié)論類似[29]。大規(guī)模的植被建設(shè)促進(jìn)了區(qū)域植被恢復(fù),區(qū)域內(nèi)生態(tài)系統(tǒng)質(zhì)量顯著增加,同時(shí)對氣候等自然因子的變化敏感度降低。本研究發(fā)現(xiàn)2000—2010年間內(nèi)蒙古一些主要的沙塵源地 (毛烏素沙地西部和科爾沁沙地東部)的植被狀況呈持續(xù)好轉(zhuǎn)趨勢,植樹種草、禁牧輪牧和防沙治沙生態(tài)恢復(fù)措施的廣泛實(shí)施在這一過程中發(fā)揮重要作用,與許多研究結(jié)果一致[30- 31];此外,內(nèi)蒙古中部地區(qū)受降水、溫度等氣候影響仍很敏感,降水的減少、溫度的升高和交通的發(fā)展,使草原區(qū)中部的大部分區(qū)域10年間植被質(zhì)量呈現(xiàn)下降趨勢。

    4 結(jié)論

    本文基于2000—2010年生物量和植被覆蓋度,并結(jié)合植被區(qū)劃數(shù)據(jù),對內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量狀況進(jìn)行分析,并評估其與氣候(降水、溫度),人類活動(dòng)(農(nóng)業(yè)、交通、生態(tài)恢復(fù)工程等)的相關(guān)關(guān)系,在此基礎(chǔ)上探討了氣候和人類活動(dòng)對內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量的影響。得出以下結(jié)論:

    (1) 內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量狀況整體偏低,2010年平均生態(tài)系統(tǒng)質(zhì)量為33.80,其中森林生態(tài)系統(tǒng)平均質(zhì)量最高,灌叢、草原生態(tài)系統(tǒng)次之??臻g分布呈明顯的經(jīng)度地帶性,由東向西,質(zhì)量逐漸降低。2000—2010年內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量總體上呈現(xiàn)緩慢增長趨勢,但局部地區(qū)生態(tài)系統(tǒng)質(zhì)量仍存在惡化,其中質(zhì)量增加的面積約為120900.43 km2。在位于100°—106°E的荒漠地區(qū),生態(tài)系統(tǒng)質(zhì)量變化趨勢較平緩;在位于107°—117°E的草原地區(qū)由于氣候變化和人類活動(dòng)的影響,生態(tài)系統(tǒng)質(zhì)量變化劇烈,主要呈降低趨勢;118°E以東地區(qū),植被生態(tài)系統(tǒng)質(zhì)量主要呈增加趨勢。

    (2) 近10年來內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量的變化與氣候和人類活動(dòng)的關(guān)系非常密切,其與降水、GDP1、化肥施用量、天保工程和退耕還草工程呈現(xiàn)明顯的正相關(guān)。而與溫度、道路密度和京津風(fēng)沙治理工程呈現(xiàn)明顯的負(fù)相關(guān)。其中,生態(tài)保護(hù)工程實(shí)施區(qū)域內(nèi)和區(qū)域外的相關(guān)性存在顯著的差異性。說明在工程區(qū)域內(nèi),生態(tài)系統(tǒng)質(zhì)量的變化主要與植樹造林、退耕還草、圍封轉(zhuǎn)移和有效管理等人類活動(dòng)密切相關(guān),而受降水、溫度等其他因素的影響較??;而工程區(qū)域外,生態(tài)系統(tǒng)質(zhì)量的變化主要受降水、溫度等因素的影響。但是,京津風(fēng)沙工程區(qū)的狀態(tài)剛好與之相反,主要受到降水、溫度的影響,而人類活動(dòng)的影響相對較小。隨著內(nèi)蒙古社會(huì)經(jīng)濟(jì)的快速發(fā)展,人類活動(dòng)對生態(tài)系統(tǒng)質(zhì)量的影響逐漸加強(qiáng),但降水仍是該地區(qū)生態(tài)系統(tǒng)質(zhì)量的主要影響因子。

    (3) 內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量增加典型區(qū)域主要位于毛烏素沙地西部、錫林郭勒高原西北部和科爾沁沙地東部地區(qū);降低典型區(qū)域主要沿鄂爾多斯高原-渾善達(dá)克沙地南緣-錫林郭勒高原東南部一線呈條帶狀分布。典型區(qū)內(nèi)質(zhì)量的增長主要是由于降水的增加、溫度的降低、農(nóng)業(yè)的發(fā)展、退耕還草工程的作用和交通發(fā)展的放緩。質(zhì)量的降低則是因?yàn)榻邓臏p少、溫度的增加、農(nóng)業(yè)發(fā)展緩慢和交通發(fā)展的加快所致。

    [1] 趙敏, 周廣勝. 基于森林資源清查資料的生物量估算模式及其發(fā)展趨勢. 應(yīng)用生態(tài)學(xué)報(bào), 2004, 15(8): 1468- 1472.

    [2] Brown S, Sathaye J, Cannell M, Kauppi P E. Mitigation of carbon emissions to the atmosphere by forest management. Commonwealth Forestry Review, 1996, 75(1): 80- 90, 109- 109, 111- 112.

    [3] Fang J Y, Chen A P, Peng C H, Zhao S Q, Ci L J. Changes in forest biomass carbon storage in China between 1949 and 1998. Science, 2001, 292(5525): 2320- 2322.

    [4] 方精云, 陳安平, 趙淑清, 慈龍駿. 中國森林生物量的估算: 對Fang等Science一文(Science, 2001, 291: 2320—2322)的若干說明. 植物生態(tài)學(xué)報(bào), 2002, 26(2): 243- 249.

    [5] 方精云, 劉國華, 徐嵩齡. 我國森林植被的生物量和凈生產(chǎn)量. 生態(tài)學(xué)報(bào), 1996, 16(5): 497- 508.

    [6] 劉國華, 傅伯杰, 方精云. 中國森林碳動(dòng)態(tài)及其對全球碳平衡的貢獻(xiàn). 生態(tài)學(xué)報(bào), 2000, 20(5): 733- 740.

    [7] 王效科, 馮宗煒, 歐陽志云. 中國森林生態(tài)系統(tǒng)的植物碳儲(chǔ)量和碳密度研究. 應(yīng)用生態(tài)學(xué)報(bào), 2001, 12(1): 13- 16.

    [8] 劉曉梅, 布仁倉, 鄧華衛(wèi), 胡遠(yuǎn)滿, 劉志華, 吳志偉. 基于地統(tǒng)計(jì)學(xué)豐林自然保護(hù)區(qū)森林生物量估測及空間格局分析. 生態(tài)學(xué)報(bào), 2011, 31(16): 4843- 4890.

    [9] Adams J E, Arkin G F. A light interception method for measuring row crop ground cover. Soil Science Society of America Journal, 1977, 41(4): 789- 792.

    [10] Gitelson A A, Kaufman Y J, Stark R, Rundquist D. Novel algorithms for remote estimation of vegetation fraction. Remote Sensing of Environment, 2002, 80(1): 76- 87.

    [11] Purevdorj T S, Tateishi R, Ishiyama T, Honda Y. Relationships between percent vegetation cover and vegetation indices. International Journal of Remote Sensing, 1998, 19(18): 3519- 3535.

    [12] Wang G, Wente S, Gertner G Z, Anderson A. Improvement in mapping vegetation cover factor for the universal soil loss equation by geostatistical methods with Landsat Thematic Mapper images. International Journal of Remote Sensing, 2002, 23(18): 3649- 3667.

    [13] Moreau S, Bosseno R, Gu X F, Baret F. Assessing the biomass dynamics of Andeanbofedalandtotorahigh-protein wetland grasses from NOAA/AVHRR. Remote Sensing of Environment, 2003, 85(4): 516- 529.

    [14] 吳云, 曾源, 趙炎, 吳炳方, 武文波. 基于MODIS數(shù)據(jù)的海河流域植被覆蓋度估算及動(dòng)態(tài)變化分析. 資源科學(xué), 2010, 32(7): 1417- 1424.

    [15] 徐涵秋, 何慧, 黃紹霖. 福建省長汀縣河田水土流失區(qū)植被覆蓋度變化及其熱環(huán)境效應(yīng). 生態(tài)學(xué)報(bào), 2013, 33(10): 2954- 2963.

    [16] 孫根年, 王美紅. 內(nèi)蒙古植被覆蓋與土地退化關(guān)系及空間結(jié)構(gòu)研究. 干旱區(qū)資源與環(huán)境, 2008, 22(2): 140 - 144.

    [17] 孫艷玲, 郭鵬, 延曉冬, 趙天寶. 內(nèi)蒙古植被覆蓋變化及其與氣候、人類活動(dòng)的關(guān)系. 自然資源學(xué)報(bào), 2010, 25(3): 407- 414.

    [18] 時(shí)忠杰, 高吉喜, 徐麗宏, 馮朝陽, 呂世海, 尚建勛. 內(nèi)蒙古地區(qū)近25年植被對氣溫和降水變化的影響. 生態(tài)環(huán)境學(xué)報(bào), 2011, 20(11): 1594- 1601.

    [19] Guo Q, Hu Z M, Li S G, Li X R, Sun X M, Yu G R. Spatial variations in aboveground net primary productivity along a climate gradient in Eurasian temperate grassland: effects of mean annual precipitation and its seasonal distribution. Global Change Biology, 2012, 18(12): 3624- 3631.

    [20] 婁雪婷, 曾源, 吳炳方, 蘇常紅, 徐延達(dá), 趙炎. 延河流域闊葉林地上生物量遙感監(jiān)測及空間分布特征分析. 資源科學(xué), 2010, 32(11): 2229- 2238.

    [21] 樸世龍, 方精云, 郭慶華. 利用CASA模型估算我國植被凈第一性生產(chǎn)力. 植物生態(tài)學(xué)報(bào), 2001, 25(5): 603- 608.

    [22] Leprieur C, Verstraete M M, Pinty B. Evaluation of the performance of various vegetation indices to retrieve vegetation cover from AVHRR data. Remote Sensing Review, 1994, 10(4): 265- 284.

    [23] Jiapaer G, Chen X, Bao A M. A comparison of methods for estimating fractional vegetation cover in arid regions. Agricultural and Forest Meteorology, 2011, 151(12): 1698- 1710.

    [24] Qi J, Marsett R C, Moran M S, Goodrich D C, Heilman P, Kerr Y H, Dedieu G, Chehbouni A, Zhang X X. Spatial and temporal dynamics of vegetation in the San Pedro River basin area. Agricultural and Forest Meteorology, 2000, 105(1/3): 55- 68.

    [25] 于小飛, 孫睿, 陳永俊, 劉紹民, 周會(huì)珍, 張忠山. 烏審旗植被覆蓋度動(dòng)態(tài)變化及其與降水量的關(guān)系. 資源科學(xué), 2006, 28(4): 31- 37.

    [26] Loyarte M M G, Menenti M. Impact of rainfall anomalies on Fourier parameters of NDVI time series of northwestern Argentina. International Journal of Remote Sensing, 2008, 29(4): 1125- 1152.

    [27] 張戈麗, 徐興良, 周才平, 張宏斌, 歐陽華. 近30年來呼倫貝爾地區(qū)草地植被變化對氣候變化的響應(yīng). 地理學(xué)報(bào), 2011, 66(1): 47- 58.

    [28] Mooney H A, Canadel J, Chapin III F S, Ehleringer J R, K?rner C, McMurtrie R E, Parton W J, Pitelka L F, Schulze E D. Ecosystem physiology responses to global change // Walker B, Steffem W, Canadell J, eds. The Terrestrial Biosphere and Global Change. Cambridge: Cambridge University Press, 1999: 190- 228.

    [29] 張寶慶, 吳普特, 趙西寧. 近30a黃土高原植被覆蓋時(shí)空演變監(jiān)測與分析. 農(nóng)業(yè)工程學(xué)報(bào), 2011, 27(4): 287- 293.

    [30] 王娟, 李寶林, 余萬里. 近30年內(nèi)蒙古自治區(qū)植被變化趨勢及影響因素分析. 干旱區(qū)資源與環(huán)境, 2012, 26(2): 132- 138.

    [31] 殷賀, 李正國, 王仰麟, 蔡福. 基于時(shí)間序列植被特征的內(nèi)蒙古荒漠化評價(jià). 地理學(xué)報(bào), 2011, 66(5): 653- 661.

    Spatial patterns of ecosystem quality in Inner Mongolia and its driving forces analysis

    XIAO Yang1,2, OUYANG Zhiyun1,*, WANG Liyan1,2, RAO Enming1,2, JIANG Ling1,2ZHANG Lu1,2

    1StateKeyLaboratoryofUrbanandRegionalEcology,ResearchCenterforEco-EnvironmentalSciences,ChineseAcademyofSciences,Beijing100085,China2GraduateUniversityofChineseAcademyofSciences,Beijing100039,China

    Vegetation is an important component of ecosystems, and it is involved in natural processes such as regional climate circulation, the hydrologic cycle, soil conditions, and other natural processes. It is also an important indicator of regional ecosystem quality change. Aboveground biomass and vegetation cover are important variables for describing vegetation quality and documenting ecosystem changes. Moreover, they are controlling factors in transpiration, photosynthesis, and other terrestrial processes, and are sensitive indicators of land degradation and desertification in arid and semi-arid regions. Aboveground biomass and vegetation cover have been used widely in the analysis of vegetation quality and changes. Inner Mongolia, an important province of northern China that is located mostly in arid and semi-arid areas, is undergoing a process of significant warming and drying, which could lead to serious ecological problems such as land degradation and desertification. Determining the driving cause of these problems has become a focus of vegetation research. Many researchers have attributed the degradation and desertification in northern China to long-term over-grazing, extensive cutting, and widespread conversion of grassland to cropland. Other studies have found that poor climate conditions, such as drought, severe wind erosion, temperatureuctuation, and winter precipitation are the primary cause of degradation and desertification. Nevertheless, some recent studies have shown that human activities can control the degradation in selected study regions of northern China. It is necessary to investigate current ecosystem quality and its response to climatic variation and human activities for a better understanding of the accumulated cons`equences of changes in climate and human activities. In this study, combining remote sensing products (aboveground biomass and vegetation cover) with vegetation regionalization data, we quantified ecosystem quality in Inner Mongolia, and analyzed the effects of climate and human activities on spatial variations in quality change from 2000 to 2010. Our results indicated that overall ecosystem quality in Inner Mongolia is very low. Among the ecosystem types, the highest quality was exhibited in forest ecosystems, followed by shrub and grassland ecosystems, respectively. Ecosystem quality displayed evident spatial heterogeneities, divided according to latitude, and reduced gradually from east to west. There was a slight increase in ecosystem quality over the study period. However, degradation and desertification also existed in some regional areas. The spatial variation characteristics of ecosystem quality are formed by interactions between the inuencing factors. The partial correlations among ecosystem quality and precipitation, GDP1, fertilizer, the natural forest protection project, and the returning farmland to grassland project were positive and statistically significant. By contrast, ecosystem quality was negatively correlated with temperature, road density, and the treatment of sand and wind sources. With the rapid development of the social economy of Inner Mongolia, the influence of human activities on ecosystem quality gradually increased, but precipitation was still the dominant factor that influenced the spatial variation characteristics of ecosystem quality. The main reasons for the improvement of ecosystem quality in a typical area were increasing precipitation, decreasing temperature, development of agriculture, returning farmland to grassland, and a slowing of traffic development. In addition, decreasing precipitation, increasing temperature, a slowdown of agricultural development, and a speeding up of traffic development led to degradation and desertification. These results provide specific information that may serve to strengthen the necessary public awareness about protecting and restoring ecosystems.

    aboveground biomass; vegetation cover; ecosystem quality; climatic interaction; human activities

    國家自然科學(xué)基金重點(diǎn)項(xiàng)目(40901289);國家科技支撐計(jì)劃課題(2011BAC09B08-02)

    2015- 01- 29;

    日期:2016- 01- 15

    10.5846/stxb201501290245

    *通訊作者Corresponding author.E-mail: zyouyang@rcees.ac.cn

    肖洋,歐陽志云,王莉雁,饒恩明,江凌,張路.內(nèi)蒙古生態(tài)系統(tǒng)質(zhì)量空間特征及其驅(qū)動(dòng)力.生態(tài)學(xué)報(bào),2016,36(19):6019- 6030.

    Xiao Y, Ouyang Z Y, Wang L Y, Rao E M, Jiang L, Zhang L.Spatial patterns of ecosystem quality in Inner Mongolia and its driving forces analysis.Acta Ecologica Sinica,2016,36(19):6019- 6030.

    猜你喜歡
    覆蓋度生物量內(nèi)蒙古
    呼和浩特市和林格爾縣植被覆蓋度變化遙感監(jiān)測
    基于NDVI的晉州市植被覆蓋信息提取
    輪牧能有效促進(jìn)高寒草地生物量和穩(wěn)定性
    在內(nèi)蒙古,奶有一百種吃法
    低覆蓋度CO分子在Ni(110)面的吸附研究
    這是美麗的內(nèi)蒙古
    草原歌聲(2019年2期)2020-01-06 03:11:12
    可愛的內(nèi)蒙古
    草原歌聲(2018年2期)2018-12-03 08:14:06
    愛在內(nèi)蒙古
    草原歌聲(2017年3期)2017-04-23 05:13:44
    生物量高的富鋅酵母的開發(fā)應(yīng)用
    基于SPOT-5遙感影像估算玉米成熟期地上生物量及其碳氮累積量
    成人高潮视频无遮挡免费网站| 在线观看免费视频日本深夜| 特级一级黄色大片| 国产成人精品久久二区二区免费| 亚洲av片天天在线观看| 欧美日韩国产亚洲二区| 亚洲专区字幕在线| 欧美绝顶高潮抽搐喷水| 欧美绝顶高潮抽搐喷水| 在线观看66精品国产| 每晚都被弄得嗷嗷叫到高潮| 一级a爱片免费观看的视频| 日韩欧美三级三区| 久久九九热精品免费| 亚洲男人的天堂狠狠| 一本大道久久a久久精品| 神马国产精品三级电影在线观看 | 夜夜夜夜夜久久久久| 国产av麻豆久久久久久久| 亚洲中文字幕一区二区三区有码在线看 | 欧美日本视频| 午夜免费成人在线视频| 色综合站精品国产| 女警被强在线播放| netflix在线观看网站| 大型黄色视频在线免费观看| 黄色视频,在线免费观看| 日本一区二区免费在线视频| 国产激情欧美一区二区| 国产爱豆传媒在线观看 | 9191精品国产免费久久| 午夜福利在线在线| 国产黄色小视频在线观看| 桃色一区二区三区在线观看| 99国产精品一区二区三区| 午夜影院日韩av| 亚洲狠狠婷婷综合久久图片| 精品久久蜜臀av无| 国产精品98久久久久久宅男小说| 丰满人妻一区二区三区视频av | 黄色女人牲交| 一本久久中文字幕| 精品国内亚洲2022精品成人| 欧美又色又爽又黄视频| 男女之事视频高清在线观看| 成年女人毛片免费观看观看9| 搡老熟女国产l中国老女人| 久久久久亚洲av毛片大全| 人人妻人人澡欧美一区二区| 一夜夜www| 久久精品国产99精品国产亚洲性色| 黑人操中国人逼视频| 亚洲av成人不卡在线观看播放网| 俄罗斯特黄特色一大片| 欧美另类亚洲清纯唯美| 久久久久国产精品人妻aⅴ院| 18禁黄网站禁片免费观看直播| 免费看美女性在线毛片视频| 法律面前人人平等表现在哪些方面| 国产亚洲精品av在线| 久久久久亚洲av毛片大全| 久久久精品欧美日韩精品| aaaaa片日本免费| 中文字幕av在线有码专区| 国内少妇人妻偷人精品xxx网站 | 精品无人区乱码1区二区| 波多野结衣巨乳人妻| 大型黄色视频在线免费观看| 亚洲无线在线观看| 午夜老司机福利片| 日韩欧美在线乱码| 国产又色又爽无遮挡免费看| 男人的好看免费观看在线视频 | 久久中文字幕一级| 欧美3d第一页| 免费在线观看视频国产中文字幕亚洲| 色噜噜av男人的天堂激情| 69av精品久久久久久| 亚洲精品国产一区二区精华液| 老熟妇乱子伦视频在线观看| 久久中文看片网| 午夜福利视频1000在线观看| 99国产精品99久久久久| 村上凉子中文字幕在线| 免费人成视频x8x8入口观看| 国产精品一区二区免费欧美| 欧美又色又爽又黄视频| 日韩av在线大香蕉| xxx96com| 久久欧美精品欧美久久欧美| 亚洲av熟女| www.精华液| 婷婷精品国产亚洲av| 精品国产美女av久久久久小说| 日韩高清综合在线| 天堂av国产一区二区熟女人妻 | 给我免费播放毛片高清在线观看| 国产精品电影一区二区三区| 亚洲国产精品999在线| 久久人妻福利社区极品人妻图片| 亚洲国产精品久久男人天堂| 悠悠久久av| 琪琪午夜伦伦电影理论片6080| 免费在线观看黄色视频的| 一级毛片高清免费大全| 国产精品久久久av美女十八| 99riav亚洲国产免费| 精品一区二区三区av网在线观看| 日韩精品青青久久久久久| 嫩草影院精品99| 成人国语在线视频| 午夜两性在线视频| 亚洲熟妇中文字幕五十中出| 国产成人一区二区三区免费视频网站| 国产精品国产高清国产av| 日韩欧美一区二区三区在线观看| 麻豆国产av国片精品| 一级毛片精品| 国产真人三级小视频在线观看| 母亲3免费完整高清在线观看| 黄色丝袜av网址大全| 亚洲精品中文字幕在线视频| 制服丝袜大香蕉在线| 成人三级做爰电影| 精品一区二区三区视频在线观看免费| 精品久久久久久久毛片微露脸| 首页视频小说图片口味搜索| 亚洲人成网站在线播放欧美日韩| 一二三四社区在线视频社区8| 可以在线观看的亚洲视频| 国产熟女xx| 国产又色又爽无遮挡免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 99国产精品一区二区三区| 婷婷六月久久综合丁香| 免费在线观看黄色视频的| 一本综合久久免费| 黄片小视频在线播放| 99精品久久久久人妻精品| 搡老熟女国产l中国老女人| 国产精品香港三级国产av潘金莲| 国产熟女午夜一区二区三区| 欧美 亚洲 国产 日韩一| 国产精品av视频在线免费观看| 高清毛片免费观看视频网站| 亚洲av中文字字幕乱码综合| 亚洲成人精品中文字幕电影| 国产亚洲av嫩草精品影院| 看黄色毛片网站| 成人三级黄色视频| 久久午夜亚洲精品久久| 99国产精品一区二区三区| 国产精品久久久久久精品电影| 久久香蕉激情| 曰老女人黄片| 一二三四社区在线视频社区8| aaaaa片日本免费| 国产精品国产高清国产av| 久久久久久久久久黄片| 久久久国产成人精品二区| 亚洲欧美精品综合一区二区三区| 亚洲无线在线观看| 一本精品99久久精品77| 国内少妇人妻偷人精品xxx网站 | 国产精品,欧美在线| 99国产综合亚洲精品| 国产精品美女特级片免费视频播放器 | 国产精品 国内视频| 欧美色视频一区免费| 99热这里只有精品一区 | 观看免费一级毛片| 激情在线观看视频在线高清| 国产精品久久久av美女十八| 成在线人永久免费视频| 国产成人av教育| 精品国内亚洲2022精品成人| 精品国产亚洲在线| 1024香蕉在线观看| 亚洲国产欧美一区二区综合| 国产激情久久老熟女| 久久久久国产精品人妻aⅴ院| 久久久国产欧美日韩av| 欧美一区二区国产精品久久精品 | 91九色精品人成在线观看| 男人舔女人下体高潮全视频| 精品一区二区三区四区五区乱码| 欧美不卡视频在线免费观看 | 又黄又粗又硬又大视频| 长腿黑丝高跟| 亚洲国产精品久久男人天堂| 国产三级在线视频| 麻豆国产av国片精品| 超碰成人久久| 视频区欧美日本亚洲| 亚洲男人的天堂狠狠| 欧美一区二区精品小视频在线| 一夜夜www| 91老司机精品| av欧美777| 欧美人与性动交α欧美精品济南到| 女人爽到高潮嗷嗷叫在线视频| 他把我摸到了高潮在线观看| 精品欧美一区二区三区在线| 伊人久久大香线蕉亚洲五| 人成视频在线观看免费观看| 日韩欧美在线乱码| 国产高清视频在线观看网站| 成年免费大片在线观看| 久久久久久免费高清国产稀缺| 99国产综合亚洲精品| 成人18禁高潮啪啪吃奶动态图| 亚洲av成人不卡在线观看播放网| 久久精品91无色码中文字幕| 国产成+人综合+亚洲专区| 黄片小视频在线播放| 久久久久久亚洲精品国产蜜桃av| 午夜视频精品福利| 亚洲五月天丁香| 中文字幕av在线有码专区| 免费av毛片视频| 免费在线观看亚洲国产| 老汉色∧v一级毛片| 一夜夜www| 啦啦啦观看免费观看视频高清| 国产成人系列免费观看| 午夜老司机福利片| 精品不卡国产一区二区三区| а√天堂www在线а√下载| 50天的宝宝边吃奶边哭怎么回事| 国产主播在线观看一区二区| 欧美3d第一页| 久久欧美精品欧美久久欧美| 给我免费播放毛片高清在线观看| av免费在线观看网站| 91成年电影在线观看| 黄色丝袜av网址大全| 2021天堂中文幕一二区在线观| 18美女黄网站色大片免费观看| 色综合婷婷激情| 美女黄网站色视频| 法律面前人人平等表现在哪些方面| 大型黄色视频在线免费观看| 真人一进一出gif抽搐免费| 岛国在线观看网站| 免费观看人在逋| 欧美一级a爱片免费观看看 | 日日摸夜夜添夜夜添小说| 亚洲精品久久成人aⅴ小说| 国产成人av教育| 熟女电影av网| 国产精品久久视频播放| 天天一区二区日本电影三级| 听说在线观看完整版免费高清| 黑人欧美特级aaaaaa片| 两个人免费观看高清视频| 很黄的视频免费| 亚洲中文字幕日韩| 亚洲五月天丁香| 国产伦在线观看视频一区| 亚洲18禁久久av| 蜜桃久久精品国产亚洲av| www.自偷自拍.com| 久久香蕉激情| 三级国产精品欧美在线观看 | svipshipincom国产片| 精品人妻1区二区| 毛片女人毛片| 日韩中文字幕欧美一区二区| 色精品久久人妻99蜜桃| 亚洲一区二区三区不卡视频| 一个人免费在线观看的高清视频| 免费观看精品视频网站| 国产亚洲精品第一综合不卡| 亚洲精品中文字幕一二三四区| 97超级碰碰碰精品色视频在线观看| 又黄又粗又硬又大视频| 99久久精品热视频| 精品欧美国产一区二区三| 波多野结衣巨乳人妻| 美女免费视频网站| a级毛片在线看网站| 99热这里只有是精品50| 国产三级在线视频| 黄色 视频免费看| 亚洲av成人一区二区三| 色在线成人网| 在线观看午夜福利视频| 免费一级毛片在线播放高清视频| 亚洲乱码一区二区免费版| 男女之事视频高清在线观看| 久久午夜综合久久蜜桃| 精品电影一区二区在线| 男女视频在线观看网站免费 | 国产精品98久久久久久宅男小说| 19禁男女啪啪无遮挡网站| 国产亚洲欧美98| 亚洲成av人片在线播放无| 亚洲精品一区av在线观看| 最近最新免费中文字幕在线| 国产激情久久老熟女| 成人特级黄色片久久久久久久| 日韩欧美三级三区| 成人手机av| 国产1区2区3区精品| 亚洲成人中文字幕在线播放| 啪啪无遮挡十八禁网站| 色尼玛亚洲综合影院| 不卡一级毛片| 听说在线观看完整版免费高清| 午夜成年电影在线免费观看| 黄色 视频免费看| 18禁美女被吸乳视频| 久久久久精品国产欧美久久久| 妹子高潮喷水视频| 最近视频中文字幕2019在线8| 美女扒开内裤让男人捅视频| 久久精品91蜜桃| 亚洲性夜色夜夜综合| 国产成人啪精品午夜网站| 12—13女人毛片做爰片一| 免费看美女性在线毛片视频| 中文在线观看免费www的网站 | 操出白浆在线播放| 757午夜福利合集在线观看| 久久中文字幕人妻熟女| 国产亚洲精品一区二区www| 色综合站精品国产| 午夜日韩欧美国产| 听说在线观看完整版免费高清| 一级毛片高清免费大全| 一个人免费在线观看的高清视频| 亚洲午夜理论影院| av在线天堂中文字幕| 精品久久久久久久人妻蜜臀av| 老汉色av国产亚洲站长工具| 欧美最黄视频在线播放免费| 国产一区二区在线av高清观看| xxx96com| 免费高清视频大片| 免费在线观看日本一区| 黄色视频不卡| 一级毛片高清免费大全| 嫩草影院精品99| 可以免费在线观看a视频的电影网站| 天天躁夜夜躁狠狠躁躁| 18禁裸乳无遮挡免费网站照片| 又黄又粗又硬又大视频| 免费高清视频大片| 久久久久国产一级毛片高清牌| 少妇熟女aⅴ在线视频| 精华霜和精华液先用哪个| 亚洲av熟女| 九色国产91popny在线| 一进一出抽搐gif免费好疼| 午夜亚洲福利在线播放| 91大片在线观看| 亚洲美女黄片视频| 在线观看免费日韩欧美大片| 国产成+人综合+亚洲专区| 美女黄网站色视频| 成人国语在线视频| 精品一区二区三区视频在线观看免费| 亚洲性夜色夜夜综合| 最近最新中文字幕大全电影3| 亚洲精品一区av在线观看| 久久九九热精品免费| 欧美绝顶高潮抽搐喷水| www.999成人在线观看| 精品午夜福利视频在线观看一区| 我的老师免费观看完整版| 大型av网站在线播放| 色综合婷婷激情| 18禁国产床啪视频网站| 成年免费大片在线观看| 日本黄大片高清| 成人午夜高清在线视频| 欧美成人免费av一区二区三区| 午夜激情av网站| 18禁国产床啪视频网站| 免费在线观看亚洲国产| 三级男女做爰猛烈吃奶摸视频| 舔av片在线| 99热这里只有是精品50| 丁香欧美五月| 久99久视频精品免费| 欧美不卡视频在线免费观看 | 精品国内亚洲2022精品成人| 国内揄拍国产精品人妻在线| 18禁黄网站禁片午夜丰满| 变态另类成人亚洲欧美熟女| 叶爱在线成人免费视频播放| 亚洲在线自拍视频| 一a级毛片在线观看| 热99re8久久精品国产| 亚洲中文字幕日韩| а√天堂www在线а√下载| 男男h啪啪无遮挡| 麻豆av在线久日| 国产精品乱码一区二三区的特点| 99在线人妻在线中文字幕| 成年免费大片在线观看| 国产乱人伦免费视频| 国产午夜精品久久久久久| 国产欧美日韩一区二区精品| 国产亚洲精品av在线| 亚洲国产精品999在线| 久久天躁狠狠躁夜夜2o2o| 日日夜夜操网爽| 国产精品久久久av美女十八| 国产亚洲av嫩草精品影院| 99国产精品一区二区三区| 久久中文看片网| 婷婷丁香在线五月| 成人亚洲精品av一区二区| 中文字幕av在线有码专区| 香蕉av资源在线| 麻豆成人午夜福利视频| 一级毛片高清免费大全| 欧美午夜高清在线| 给我免费播放毛片高清在线观看| 成年免费大片在线观看| 麻豆国产97在线/欧美 | 757午夜福利合集在线观看| 亚洲人成电影免费在线| 成年人黄色毛片网站| 久久午夜亚洲精品久久| 国产真实乱freesex| 毛片女人毛片| tocl精华| 国产av又大| 久久草成人影院| 欧美中文综合在线视频| 久久热在线av| 久久精品亚洲精品国产色婷小说| 国产精品久久久人人做人人爽| 日韩欧美精品v在线| 国产又色又爽无遮挡免费看| 波多野结衣高清作品| 一二三四社区在线视频社区8| 久久精品影院6| 国产黄色小视频在线观看| 欧美色视频一区免费| 免费在线观看亚洲国产| 亚洲国产精品久久男人天堂| 亚洲 欧美一区二区三区| 成人精品一区二区免费| 亚洲av五月六月丁香网| 一级黄色大片毛片| 巨乳人妻的诱惑在线观看| av福利片在线| 老司机深夜福利视频在线观看| 男人舔女人下体高潮全视频| 久久精品国产清高在天天线| 一级a爱片免费观看的视频| 免费av毛片视频| 天堂√8在线中文| 51午夜福利影视在线观看| 又大又爽又粗| bbb黄色大片| 午夜久久久久精精品| 一级毛片高清免费大全| 免费无遮挡裸体视频| 此物有八面人人有两片| 在线观看一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 制服丝袜大香蕉在线| 毛片女人毛片| 狂野欧美激情性xxxx| 亚洲精品av麻豆狂野| 欧美精品啪啪一区二区三区| 成人av一区二区三区在线看| av欧美777| 国产精品日韩av在线免费观看| 亚洲五月天丁香| 两人在一起打扑克的视频| 久久久精品大字幕| 好看av亚洲va欧美ⅴa在| 国产亚洲av高清不卡| 无遮挡黄片免费观看| 99热只有精品国产| av视频在线观看入口| 在线观看日韩欧美| 国产精品久久电影中文字幕| 亚洲一区二区三区不卡视频| 午夜福利欧美成人| 桃红色精品国产亚洲av| 亚洲av片天天在线观看| 亚洲国产看品久久| 18禁裸乳无遮挡免费网站照片| 免费一级毛片在线播放高清视频| 村上凉子中文字幕在线| 人成视频在线观看免费观看| 亚洲欧洲精品一区二区精品久久久| 国产爱豆传媒在线观看 | 精品日产1卡2卡| 午夜激情av网站| www.999成人在线观看| 无人区码免费观看不卡| 天天一区二区日本电影三级| 日本精品一区二区三区蜜桃| 婷婷丁香在线五月| 麻豆久久精品国产亚洲av| 国产黄色小视频在线观看| 人妻丰满熟妇av一区二区三区| 无遮挡黄片免费观看| 欧美成狂野欧美在线观看| 国产成+人综合+亚洲专区| 91麻豆精品激情在线观看国产| 精品熟女少妇八av免费久了| 曰老女人黄片| 亚洲精品国产精品久久久不卡| svipshipincom国产片| 午夜福利在线在线| 亚洲精品久久国产高清桃花| 亚洲自偷自拍图片 自拍| www.999成人在线观看| 亚洲国产精品久久男人天堂| 男女午夜视频在线观看| 国产精华一区二区三区| www日本在线高清视频| 亚洲精品国产精品久久久不卡| 国产片内射在线| 国产野战对白在线观看| 制服诱惑二区| 久久精品亚洲精品国产色婷小说| 国产成年人精品一区二区| 特大巨黑吊av在线直播| 12—13女人毛片做爰片一| 亚洲精华国产精华精| 一本综合久久免费| 男女床上黄色一级片免费看| 午夜福利成人在线免费观看| 日韩免费av在线播放| 桃红色精品国产亚洲av| 日韩精品中文字幕看吧| 琪琪午夜伦伦电影理论片6080| 国产又黄又爽又无遮挡在线| 亚洲欧美日韩无卡精品| 久久久国产成人免费| 国产精品一区二区精品视频观看| 国产伦人伦偷精品视频| 真人做人爱边吃奶动态| 香蕉av资源在线| 国产97色在线日韩免费| 12—13女人毛片做爰片一| 中国美女看黄片| 岛国在线观看网站| 国内毛片毛片毛片毛片毛片| 中文字幕久久专区| 久久天堂一区二区三区四区| 一卡2卡三卡四卡精品乱码亚洲| 国产精品久久久久久亚洲av鲁大| 午夜免费激情av| www日本在线高清视频| 中文字幕熟女人妻在线| 国内精品久久久久精免费| 丰满人妻一区二区三区视频av | 中文字幕高清在线视频| 欧美黑人欧美精品刺激| 伦理电影免费视频| 久久精品夜夜夜夜夜久久蜜豆 | 一区二区三区国产精品乱码| 亚洲,欧美精品.| 丰满人妻熟妇乱又伦精品不卡| 国产成人系列免费观看| bbb黄色大片| 老司机靠b影院| 亚洲人成网站在线播放欧美日韩| 久久久久久大精品| 国产99久久九九免费精品| 欧美午夜高清在线| 午夜成年电影在线免费观看| 在线视频色国产色| av片东京热男人的天堂| 日韩免费av在线播放| 国产av又大| 不卡一级毛片| 黄片大片在线免费观看| 精品福利观看| 岛国在线免费视频观看| 老司机深夜福利视频在线观看| 中文字幕熟女人妻在线| 丁香六月欧美| 天堂√8在线中文| 精品一区二区三区av网在线观看| 久久精品国产亚洲av香蕉五月| 欧美一级毛片孕妇| 国产一区在线观看成人免费| 成熟少妇高潮喷水视频| 久久国产乱子伦精品免费另类| 国产亚洲欧美98| 神马国产精品三级电影在线观看 | 又紧又爽又黄一区二区| 国产精品久久久久久人妻精品电影| 欧美乱色亚洲激情| 香蕉丝袜av| 欧美激情久久久久久爽电影| 国产精品久久久久久人妻精品电影| 在线观看66精品国产| 男女视频在线观看网站免费 | av有码第一页| 真人做人爱边吃奶动态| 90打野战视频偷拍视频| 久久 成人 亚洲| 天堂动漫精品| 99精品在免费线老司机午夜| 久久香蕉国产精品| 亚洲 欧美一区二区三区| 亚洲中文av在线| 成人精品一区二区免费| 90打野战视频偷拍视频| 国产1区2区3区精品| 国产三级黄色录像| 可以在线观看毛片的网站| 国产一区二区在线av高清观看|