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

    基于生態(tài)系統(tǒng)服務(wù)價(jià)值的生態(tài)足跡計(jì)算與分析*
    ——以瑪納斯河流域?yàn)槔?/h1>
    2022-03-25 03:54:28李鵬輝張茹倩徐麗萍
    關(guān)鍵詞:瑪納斯當(dāng)量足跡

    李鵬輝,張茹倩,徐麗萍

    (石河子大學(xué)理學(xué)院,新疆石河子 832000)

    0 引言

    自然資本和生態(tài)服務(wù)是人類賴以生存和現(xiàn)代文明得以延續(xù)的基礎(chǔ),是社會(huì)經(jīng)濟(jì)可持續(xù)發(fā)展的重要支持系統(tǒng)[1]。隨著工業(yè)化、城鎮(zhèn)化、全球化的快速推進(jìn),人類對(duì)于自然資源的過速消耗所帶來的一系列生態(tài)環(huán)境效益縮減和區(qū)域發(fā)展瓶頸問題愈加凸顯。據(jù)IPBES 全球評(píng)估報(bào)告,受氣候變化和人類活動(dòng)共同影響,全球正面臨自然衰退“史無前例”和物種滅絕率“加速”的局面[2]。構(gòu)建科學(xué)的評(píng)價(jià)方法和指標(biāo)體系,充分反映區(qū)域的資源利用和生態(tài)狀況已成為當(dāng)前亟待解決的熱點(diǎn)問題。生態(tài)足跡將人類發(fā)展所消耗的資源及其產(chǎn)生的廢棄物轉(zhuǎn)化為相應(yīng)的生物生產(chǎn)性土地面積,實(shí)現(xiàn)了對(duì)自然資本消耗的統(tǒng)一性核算,為可持續(xù)評(píng)價(jià)提供了有效途徑[3]。憑借其框架簡明、思路清晰、操作簡便的獨(dú)特優(yōu)勢(shì),生態(tài)足跡在不同空間尺度上都展示了良好的普適性,并被廣泛地應(yīng)用在生態(tài)安全[4]、生態(tài)補(bǔ)償[5]、生態(tài)效率[6]、綠色發(fā)展[7]、土地優(yōu)化[8]等相關(guān)領(lǐng)域,但該方法仍存在部分缺憾。理論層面,生態(tài)足跡理論側(cè)重于土地的生物生產(chǎn)性功能,僅考慮土地初級(jí)產(chǎn)品供給能力,忽略了生態(tài)系統(tǒng)的氣候調(diào)節(jié)、土壤保持、美學(xué)景觀等其他服務(wù)功能;方法層面,生態(tài)足跡模型忽視了區(qū)域社會(huì)經(jīng)濟(jì)和地理環(huán)境的異質(zhì)性,易導(dǎo)致間接折算誤差[9],且通用因子在局域小尺度上的適用性不佳[10]。

    生態(tài)系統(tǒng)服務(wù)是可持續(xù)發(fā)展研究的核心內(nèi)容之一,單一的生態(tài)足跡或生態(tài)系統(tǒng)服務(wù)價(jià)值方法均難以充分反映區(qū)域生態(tài)環(huán)境質(zhì)量的實(shí)際狀況[11],如何將生態(tài)系統(tǒng)服務(wù)與生態(tài)足跡模型有機(jī)結(jié)合,仍在進(jìn)一步探索中。閔慶文等[12]在對(duì)傳統(tǒng)生態(tài)足跡局限性分析的基礎(chǔ)上,提出基于生態(tài)系統(tǒng)服務(wù)的生態(tài)足跡概念,并構(gòu)建污染足跡模型進(jìn)一步完善。張義和焦雯珺等[13,14]以水生態(tài)系統(tǒng)為突破點(diǎn),基于水生態(tài)系統(tǒng)的多種生態(tài)服務(wù)功能,構(gòu)建了基于生態(tài)系統(tǒng)服務(wù)的水生態(tài)足跡。曹智等[15]認(rèn)為生態(tài)系統(tǒng)服務(wù)是生態(tài)系統(tǒng)提供給人類生存和發(fā)展的限制性資源,以“生態(tài)系統(tǒng)—生態(tài)系統(tǒng)服務(wù)—人口和經(jīng)濟(jì)(承載力)”為研究主線,構(gòu)建了基于生態(tài)系統(tǒng)服務(wù)的生態(tài)承載力評(píng)估模型。曹瑞芬等[16]和陳芳淼等[17]將生態(tài)服務(wù)價(jià)值引入到土地資源供需平衡研究中,并展開生態(tài)補(bǔ)償和土地資源承載力評(píng)估。王恒博等[18]將生態(tài)足跡與生態(tài)服務(wù)價(jià)值法相結(jié)合,提出生態(tài)足跡-服務(wù)價(jià)值法,并對(duì)山西省縣域生態(tài)承載力進(jìn)行時(shí)空動(dòng)態(tài)分析。朱文博等[19]基于生態(tài)系統(tǒng)服務(wù)足跡分析了中國省域間食物供給足跡、淡水供給足跡、固碳服務(wù)足跡的動(dòng)態(tài)流動(dòng)規(guī)律。郭慧等[20]基于生態(tài)系統(tǒng)服務(wù)價(jià)值測算了北京市門頭溝區(qū)耕地、林地、水域的均衡因子和產(chǎn)量因子。上述研究對(duì)生態(tài)系統(tǒng)服務(wù)與生態(tài)足跡的結(jié)合進(jìn)行了有益探索,極大地豐富了相關(guān)研究成果,但同時(shí)也存在部分缺陷,如放棄均衡因子、產(chǎn)量因子,只能分類比較,無法形成整體判斷,或僅展開單一類型生態(tài)系統(tǒng)研究,無法實(shí)現(xiàn)對(duì)區(qū)域生態(tài)狀況的統(tǒng)一核算。

    瑪納斯河流域?qū)俚湫偷母珊祬^(qū)內(nèi)陸河流域,該流域是“天山北坡經(jīng)濟(jì)帶”和“絲綢之路”經(jīng)濟(jì)帶核心區(qū)的重要組成部分。隨著人類活動(dòng)的不斷增強(qiáng),人工綠洲規(guī)模持續(xù)擴(kuò)張,流域濕地退化嚴(yán)重,生態(tài)服務(wù)價(jià)值整體下降[21],區(qū)域發(fā)展面臨著巨大的現(xiàn)實(shí)困境及潛在的生態(tài)安全風(fēng)險(xiǎn),亟需對(duì)該流域的生態(tài)狀況進(jìn)行系統(tǒng)評(píng)價(jià)。文章在生態(tài)足跡的模型框架下,以生態(tài)系統(tǒng)服務(wù)價(jià)值理論為指導(dǎo)思想,構(gòu)建了基于生態(tài)系統(tǒng)服務(wù)價(jià)值的生態(tài)足跡模型,對(duì)瑪納斯河流域的生態(tài)足跡進(jìn)行統(tǒng)一核算,以期為流域的社會(huì)經(jīng)濟(jì)發(fā)展和綠洲生態(tài)穩(wěn)定提供參考依據(jù),同時(shí)有助于進(jìn)一步完善和補(bǔ)充基于生態(tài)系統(tǒng)服務(wù)的生態(tài)足跡模型研究。

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

    1.1 研究區(qū)概況

    瑪納斯河流域地處天山北麓中段、準(zhǔn)噶爾盆地南緣,地理位置處于43°05′N~45°58′N,85°01′E~86°32′E,行政區(qū)劃上包括石河子市、瑪納斯縣、沙灣縣、兵團(tuán)第六師的新湖總場以及克拉瑪依市的小拐鄉(xiāng)等(圖1)。流域發(fā)育有干旱區(qū)典型的“山地—綠洲—荒漠”生態(tài)系統(tǒng),上游山區(qū)是流域的集水匯流區(qū)和水源涵養(yǎng)地,中游綠洲區(qū)是流域水土資源的精華,撫育了89.4 萬人,創(chuàng)造的經(jīng)濟(jì)總值年均增長12.4%[22],下游荒漠區(qū)約占流域總面積的1/3。該流域灌溉農(nóng)業(yè)發(fā)達(dá),是新疆重要的商品棉和糧食生產(chǎn)基地,但長期以來農(nóng)業(yè)用水嚴(yán)重?cái)D占工業(yè)、生活和生態(tài)用水,進(jìn)一步加劇了流域水資源的供需矛盾,同時(shí)也誘發(fā)了濕地退化、鹽漬化擴(kuò)大、生態(tài)服務(wù)價(jià)值下降等生態(tài)環(huán)境效應(yīng)[23,24]。21世紀(jì)以來,瑪納斯河流域工業(yè)化和城鎮(zhèn)化發(fā)展進(jìn)入快速階段,人類干擾活動(dòng)愈加強(qiáng)烈,綠洲規(guī)模持續(xù)擴(kuò)張,土地利用類型轉(zhuǎn)換復(fù)雜,荒漠化與綠洲化演替頻繁[25],目前該流域已成為我國典型的生態(tài)脆弱區(qū)[26]??傮w而言,瑪納斯河流域是一個(gè)生態(tài)單元完整、經(jīng)濟(jì)發(fā)展快速而潛在生態(tài)安全風(fēng)險(xiǎn)巨大,且相對(duì)獨(dú)立的干旱區(qū)內(nèi)陸河流域。

    圖1 研究區(qū)概況

    1.2 數(shù)據(jù)來源及處理

    該研究采用的數(shù)據(jù)主要包括瑪納斯河流域土地利用數(shù)據(jù)和社會(huì)經(jīng)濟(jì)統(tǒng)計(jì)數(shù)據(jù)兩類,前者主要用來測算流域的生態(tài)承載力,后者主要用來測算模型參數(shù)和流域生態(tài)足跡。流域土地利用數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(數(shù)據(jù)獲取地址:http://www.resdc.cn/),空間分辨率為30m,結(jié)合流域景觀特征和研究需要,依據(jù)“中國土地利用/土地覆蓋遙感監(jiān)測數(shù)據(jù)分類系統(tǒng)”對(duì)2000年、2005年、2010年、2015年、2018 年研究區(qū)的土地利用類型進(jìn)行重分類,分為耕地、林地、草地、水域、建設(shè)用地和未利用地6類。在均衡因子和產(chǎn)量因子計(jì)算過程中,全國不同生態(tài)系統(tǒng)的當(dāng)量因子采用謝高地等的研究成果[27],全國糧食作物的種植和生產(chǎn)數(shù)據(jù)來自《中國統(tǒng)計(jì)年鑒》,各類糧食作物的價(jià)格來自于相應(yīng)年份的《中國農(nóng)產(chǎn)品價(jià)格調(diào)查年鑒》。流域生態(tài)足跡的核算內(nèi)容及社會(huì)統(tǒng)計(jì)基礎(chǔ)數(shù)據(jù)來源見表1。

    2 研究方法

    2.1 基于生態(tài)系統(tǒng)服務(wù)價(jià)值的生態(tài)足跡模型

    2.1.1 生態(tài)足跡模型

    生態(tài)足跡的核心思想是將區(qū)域發(fā)展消耗的資源及排放的廢棄物轉(zhuǎn)化為一定面積的生物生產(chǎn)性土地,進(jìn)而從供需兩個(gè)層面的比較來判斷區(qū)域可持續(xù)發(fā)展?fàn)顩r[28]。雖然有效地測度了人類對(duì)自然環(huán)境的生態(tài)占用狀況,但忽略了土地的多功能性,無法體現(xiàn)資源消耗所帶來的生態(tài)服務(wù)功能變化,難以精準(zhǔn)追蹤生態(tài)系統(tǒng)服務(wù)價(jià)值的流動(dòng)過程[29]。因此,該文基于生態(tài)足跡和生態(tài)系統(tǒng)服務(wù)價(jià)值兩大理論,構(gòu)建了基于生態(tài)系統(tǒng)服務(wù)價(jià)值的生態(tài)足跡模型(圖2)。此時(shí),生態(tài)足跡反映生態(tài)系統(tǒng)的需求,表示為維持一個(gè)區(qū)域人類生存和發(fā)展所需要的具有完整生態(tài)服務(wù)功能的地域面積;生態(tài)承載力反映生態(tài)系統(tǒng)的供給,表示一個(gè)區(qū)域能夠提供給人類社會(huì)發(fā)展的各類生態(tài)系統(tǒng)面積總和。

    圖2 基于生態(tài)系統(tǒng)服務(wù)價(jià)值的生態(tài)足跡模型框架

    式(1)(2)(3)中,EF表示生態(tài)足跡,N為區(qū)域人口數(shù)量,ef為人均生態(tài)足跡,ri為第i類生態(tài)系統(tǒng)的均衡因子,aai為人均占用交易商品折算后的第i類生態(tài)系統(tǒng)面積,EC表示生態(tài)承載力,ec為人均生態(tài)承載力,ai為實(shí)際人均占有的第i類生態(tài)系統(tǒng)面積,yi為第i類生態(tài)系統(tǒng)的產(chǎn)量因子,ER表示生態(tài)盈余,ED表示生態(tài)赤字。當(dāng)EC>EF時(shí),表示生態(tài)系統(tǒng)服務(wù)供給能力>耗損能力,區(qū)域處于生態(tài)盈余狀況,EC-EF=ER;當(dāng)EC<EF時(shí),表示生態(tài)系統(tǒng)服務(wù)供給能力<耗損能力,區(qū)域處于生態(tài)赤字狀況,EC-EF=ED。

    2.1.2 均衡因子

    在生態(tài)足跡測算過程中,由于不同類型生物生產(chǎn)性土地的生產(chǎn)能力不同,無法直接比較,因此通過均衡因子將其轉(zhuǎn)化為同一生產(chǎn)性質(zhì)的土地面積,進(jìn)行統(tǒng)一核算。在基于生態(tài)系統(tǒng)服務(wù)價(jià)值的生態(tài)足跡框架中,均衡因子表示某一生態(tài)系統(tǒng)類型的生態(tài)服務(wù)能力與區(qū)域所有生態(tài)系統(tǒng)的平均生態(tài)服務(wù)能力的差異,可由研究區(qū)某一生態(tài)系統(tǒng)類型單位面積提供的生態(tài)系統(tǒng)服務(wù)價(jià)值與該區(qū)域所有生態(tài)系統(tǒng)類型平均單位面積所能提供的生態(tài)系統(tǒng)服務(wù)價(jià)值的比值表示[20]。計(jì)算公式為:

    式(4)中,Pi為研究區(qū)第i類生態(tài)系統(tǒng)單位面積的生態(tài)系統(tǒng)服務(wù)價(jià)值量,單位:元/(hm2·年)為研究區(qū)所有生態(tài)系統(tǒng)單位面積的平均生態(tài)系統(tǒng)服務(wù)價(jià)值量,單位:元/(hm2·年);Dt為t年度研究區(qū)1個(gè)標(biāo)準(zhǔn)當(dāng)量因子的生態(tài)系統(tǒng)服務(wù)價(jià)值量,單位:元/hm2;Fi為研究區(qū)第i類生態(tài)系統(tǒng)的服務(wù)價(jià)值當(dāng)量因子之和;Si為研究區(qū)第i類生態(tài)系統(tǒng)的面積,單位:hm2。

    2.1.3 產(chǎn)量因子

    在生態(tài)承載力測算過程中,由于同一生物生產(chǎn)性土地在不同區(qū)域之間的生產(chǎn)能力存在差異,因此通過產(chǎn)量因子的橋梁作用,使區(qū)域之間的生物生產(chǎn)性土地面積具有可比性。在基于生態(tài)系統(tǒng)服務(wù)價(jià)值的生態(tài)足跡模型中,產(chǎn)量因子表示區(qū)域某一生態(tài)系統(tǒng)類型的生態(tài)服務(wù)能力與全國同類生態(tài)系統(tǒng)的平均生態(tài)服務(wù)能力的差異,可由研究區(qū)某一生態(tài)系統(tǒng)類型單位面積提供的生態(tài)系統(tǒng)服務(wù)價(jià)值與全國同類的生態(tài)系統(tǒng)平均單位面積所能提供的生態(tài)系統(tǒng)服務(wù)價(jià)值的比值表示[20]。計(jì)算公式為:

    式(5)中,為全國第i類生態(tài)系統(tǒng)單位面積的生態(tài)系統(tǒng)服務(wù)價(jià)值量,單位:元/(hm2·年);為t年度全國1 個(gè)標(biāo)準(zhǔn)當(dāng)量因子的生態(tài)系統(tǒng)服務(wù)價(jià)值量,單位:元/hm2;為全國第i類生態(tài)系統(tǒng)的服務(wù)價(jià)值當(dāng)量因子之和;Pi、Dt、Fi含義同上。

    2.2 單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值測算

    依托謝高地等[27]提出的中國單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子,借鑒已有相關(guān)成果,結(jié)合干旱區(qū)內(nèi)陸河流域不同土地類型的生態(tài)系統(tǒng)特征和綠洲實(shí)際生產(chǎn)活動(dòng)狀況,修正了瑪納斯河流域的單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子表(表2)。由于瑪納斯河流域年均降水量不足200mm,蒸散量高達(dá)1 521mm[30],水文調(diào)節(jié)能力顯著低于全國平均水平[31],同時(shí)流域冰川積雪面積達(dá)6.76×104hm2,占水域面積的33.26%(冰川積雪的水文調(diào)節(jié)當(dāng)量僅7.13[27]),結(jié)合朱增云等[32]的研究成果,最終將研究區(qū)水域的水文調(diào)節(jié)當(dāng)量確定為44.53。謝高地等[33]提出將單位面積農(nóng)田生態(tài)系統(tǒng)糧食生產(chǎn)的凈利潤作為1 個(gè)標(biāo)準(zhǔn)當(dāng)量因子的生態(tài)系統(tǒng)服務(wù)價(jià)值量,但查閱《全國農(nóng)產(chǎn)品成本收益資料匯編》發(fā)現(xiàn)部分年份(如2016年、2017年和2018年)糧食生產(chǎn)凈利潤為負(fù)值,導(dǎo)致該方法失效。且已有研究表明,沒有生產(chǎn)力投入的自然土地利用/覆被類型生態(tài)系統(tǒng)所提供的服務(wù)價(jià)值是現(xiàn)有農(nóng)田單位面積所能提供的食物生產(chǎn)價(jià)值的1/7[34,35]。因此標(biāo)準(zhǔn)當(dāng)量因子的生態(tài)系統(tǒng)服務(wù)價(jià)值量可由下式得到:

    式(6)(7)(8)中,Dt為t年度1個(gè)標(biāo)準(zhǔn)當(dāng)量因子的生態(tài)系統(tǒng)服務(wù)價(jià)值量,j為糧食作物種類,mj為第j類糧食作物的播種面積,單位:hm2;pj為第j類糧食作物的價(jià)格,單位:元/kg;qj為第j類糧食作物的單位面積產(chǎn)量,單位:kg/hm2;Mt為t年度n種糧食作物的總播種面積。Fi為第i類生態(tài)系統(tǒng)的服務(wù)價(jià)值當(dāng)量因子之和,eik表示第i類生態(tài)系統(tǒng)的第k種生態(tài)服務(wù)價(jià)值的當(dāng)量因子(表2)。Pi為第i類生態(tài)系統(tǒng)單位面積的生態(tài)系統(tǒng)服務(wù)價(jià)值量。

    表2 瑪納斯河流域單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子

    2.3 敏感性指數(shù)

    基于生態(tài)系統(tǒng)服務(wù)價(jià)值的生態(tài)足跡模型中均衡因子和產(chǎn)量因子的確定均與研究區(qū)單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值的當(dāng)量因子有關(guān),為反映測算結(jié)果對(duì)當(dāng)量因子的依賴程度,該文引入基于經(jīng)濟(jì)學(xué)中彈性系數(shù)的敏感性指數(shù)(Coefficient of Sensitive,CS)分析方法[36],用兩個(gè)變量變動(dòng)的比率,來驗(yàn)證當(dāng)量因子對(duì)生態(tài)足跡測算結(jié)果的影響程度。

    式(9)中,CS為敏感性指數(shù),Y為生態(tài)足跡的測算結(jié)果,e為當(dāng)量因子,0、1 分別表示初始值和調(diào)整50%之后的數(shù)值。當(dāng)CS=0 時(shí),表明Y對(duì)e完全無彈性;當(dāng)0<CS<1 時(shí),表明Y對(duì)e缺乏彈性;當(dāng)CS=1時(shí),表明Y對(duì)e單位彈性;當(dāng)CS>1時(shí),表明Y對(duì)e富有彈性。CS值越大,說明Y對(duì)e的依賴度越高,當(dāng)量因子的準(zhǔn)確性對(duì)測算結(jié)果越關(guān)鍵。

    3 結(jié)果與分析

    3.1 瑪納斯河流域生態(tài)足跡模型參數(shù)分析

    3.1.1 均衡因子

    瑪納斯河流域耕地、林地、草地、水域和未利用地生態(tài)系統(tǒng)的均衡因子存在年際差異,均是在2015年達(dá)到最高值,在2000 年處于最低值,但總體波動(dòng)幅度不大(表3)。世界自然基金會(huì)(WWF)發(fā)布的《Living Planet Repot 2012》中將能源用地定義為“用于吸收化石能源燃燒排放的溫室氣體的森林”[37],但實(shí)際上草地和耕地也同樣具有碳匯功能,考慮到瑪納斯河流域林地僅占總面積的4%左右,若僅以林地為測算對(duì)象,難以反映流域?qū)嶋H情況,因此能源用地的均衡因子平均值由林地、草地和耕地按比例確定為1.39。隨著城鎮(zhèn)化的快速發(fā)展,建設(shè)用地多由耕地轉(zhuǎn)入,因此建設(shè)用地的均衡因子選取與傳統(tǒng)模型一樣,由耕地?cái)?shù)值代替。各類生態(tài)系統(tǒng)的均衡因子平均值按照大小排序?yàn)椋核颍?6.87)>林地(4.62)>能源用地(1.39)>草地(1.37)>耕地(1.07)=建設(shè)用地(1.07)>未利用地(0.30)。其中水域和林地的均衡因子顯著高于其他土地類型,這與干旱區(qū)內(nèi)陸河流域的資源稟賦特征密切相關(guān),充分體現(xiàn)了水資源和森林資源對(duì)于綠洲生態(tài)穩(wěn)定和可持續(xù)發(fā)展的重要意義。

    表3 2000—2018年瑪納斯河流域不同生態(tài)系統(tǒng)的均衡因子

    3.1.2 產(chǎn)量因子

    瑪納斯河流域耕地、林地、草地、水域和未利用地生態(tài)系統(tǒng)的產(chǎn)量因子存在年際差異,其中耕地的波動(dòng)幅度最大,未利用地的波動(dòng)幅度最?。ū?)。由于自然系統(tǒng)中沒有專門用來生產(chǎn)能源的土地,傳統(tǒng)生態(tài)足跡模型將能源用地的產(chǎn)量因子確定為0,建設(shè)用地的產(chǎn)量因子則由耕地?cái)?shù)值代替[38]。為使核算內(nèi)容完整統(tǒng)一,該研究對(duì)能源用地和建設(shè)用地產(chǎn)量因子的處理方法與傳統(tǒng)生態(tài)足跡模型保持一致。各類生態(tài)系統(tǒng)的產(chǎn)量因子平均值按照大小排序?yàn)椋焊兀?.42)=建設(shè)用地(1.42)>水域(1.31)>林地(1.04)>草地(0.95)>未利用地(0.45),可以發(fā)現(xiàn)耕地、水域和林地生態(tài)系統(tǒng)單位面積的生態(tài)服務(wù)價(jià)值高于全國水平,而草地和未利用地則低于全國平均水平。其主要原因是瑪納斯河流域是新疆重要的商品棉和糧食生產(chǎn)基地,糧食單產(chǎn)較高于全國平均水平,從而引起標(biāo)準(zhǔn)當(dāng)量因子值較高。對(duì)比發(fā)現(xiàn),各類生態(tài)系統(tǒng)產(chǎn)量因子的最大值出現(xiàn)在2005年,最小值出現(xiàn)在2015年。這主要與兩個(gè)年份標(biāo)準(zhǔn)當(dāng)量因子的生態(tài)系統(tǒng)服務(wù)價(jià)值量有關(guān),受作物收成和市場價(jià)格影響,2005 年瑪納斯河流域的標(biāo)準(zhǔn)當(dāng)量因子數(shù)值是全國的1.6倍,而在2015年是1.2倍。

    表4 2000—2018年瑪納斯河流域不同生態(tài)系統(tǒng)的產(chǎn)量因子

    3.2 瑪納斯河流域生態(tài)足跡測算結(jié)果分析

    3.2.1 瑪納斯河流域生態(tài)赤字呈擴(kuò)大之勢(shì)

    由于流域和全國的單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值均存在年度波動(dòng),難以厘清兩種變動(dòng)對(duì)產(chǎn)量因子的具體貢獻(xiàn),為避免可能導(dǎo)致的“生態(tài)足跡幻覺”,該研究用均衡因子和產(chǎn)量因子的多年平均值對(duì)瑪納斯河流域的生態(tài)足跡和生態(tài)承載力進(jìn)行測算。結(jié)果顯示,瑪納斯河流域生態(tài)足跡持續(xù)增加,生態(tài)承載力波動(dòng)變化,可持續(xù)狀態(tài)由生態(tài)盈余轉(zhuǎn)為生態(tài)赤字,且生態(tài)赤字呈擴(kuò)大之勢(shì)(圖3)。流域生態(tài)足跡在2000—2018 年呈增加趨勢(shì),由2000 年的2.58×106hm2增加到2018 年的11.25×106hm2,增長速度達(dá)0.48×106hm2/年;而生態(tài)承載力則呈波動(dòng)變化趨勢(shì),在均值3.68×106hm2附近波動(dòng),變化幅度在-0.41×106~0.49×106hm2,最高值出現(xiàn)在2005 年,為4.18×106hm2,最低值出現(xiàn)在2015 年,為3.27×106hm2;流域2000年處于生態(tài)盈余(0.95×106hm2),之后進(jìn)入生態(tài)赤字狀況,2018 年赤字高達(dá)-7.49×106hm2,年均增速為-0.47×106hm2/年。從人均水平來看,瑪納斯河流域的人均生態(tài)足跡穩(wěn)定增加,由2000 年的3.76hm2/人增加至2018年的12.59hm2/人,18年間增加了2.35倍;而人均生態(tài)承載力呈現(xiàn)波動(dòng)變化,最高值出現(xiàn)在2005年為5.99 hm2/人,最低值出現(xiàn)在2015 年為3.88 hm2/人,總體上呈略有下降的趨勢(shì),由2000 年的5.14 hm2/人下降至2018 年的4.21 hm2/人;人均生態(tài)赤字快速擴(kuò)大,由2000 年的生態(tài)盈余(1.38 hm2/人)增大到2018年的-8.38 hm2/人,年均擴(kuò)張速度達(dá)-0.54 hm2/(人·年)。其中,生態(tài)足跡的變化主要受人口增長和經(jīng)濟(jì)高速發(fā)展影響較大,而生態(tài)承載力的變化則與流域的耕地?cái)U(kuò)張和造林活動(dòng)關(guān)系密切。

    圖3 2000—2018年瑪納斯河流域生態(tài)足跡與生態(tài)承載力

    3.2.2 瑪納斯河流域生態(tài)足跡呈不均衡變化

    瑪納斯河流域生態(tài)足跡和生態(tài)承載力均呈現(xiàn)不均衡變化(表5)。生態(tài)足跡各組分差異懸殊,各地類對(duì)流域生態(tài)足跡總量的平均貢獻(xiàn)率依次為:草地(54.71%)>能源用地(30.34%)>耕地(8.16%)>水域(4.27%)>林地(1.51%)>建設(shè)用地(1.01%);變化趨勢(shì)上,除林地、能源用地在2010 年有所下降隨后上升外,其余類型的生態(tài)足跡均呈現(xiàn)出逐年增加的趨勢(shì),各類型生態(tài)足跡增長速度依次為:草地(33.54×104hm2/年)>能源用地(9.94×104hm2/年)>耕地(2.26×104hm2/年)>水域(1.82×104hm2/年)>林地(0.38×104hm2/年)>建設(shè)用地(0.22×104hm2/年)。生態(tài)承載力各組分變化差異較大,水域和草地的生態(tài)承載力呈下降趨勢(shì),下降速度分別為3.55×104hm2/年和0.25×104hm2/年;耕地持續(xù)增加,增長速度達(dá)3.92×104hm2/年;林地和建設(shè)用地在2010 年略有下降隨后上升,增長速度分別為0.68×104hm2/年和0.33×104hm2/年;未利用地呈波動(dòng)變化,總體增加速度為0.16×104hm2/年。各土地類型對(duì)流域生態(tài)承載力總量的貢獻(xiàn)能力也不相同,以2018 年為例,耕地(37.39%)>草地(29.46%)>水域(13.06%)>林地(11.13%)>未利用地(5.83%)>建設(shè)用地(3.13%)。

    表5 2000—2018年瑪納斯河流域生態(tài)足跡與生態(tài)承載力計(jì)算結(jié)果 萬hm2

    3.3 瑪納斯河流域生態(tài)足跡敏感性指數(shù)分析

    瑪納斯河流域基于生態(tài)系統(tǒng)服務(wù)價(jià)值的生態(tài)足跡測算結(jié)果對(duì)當(dāng)量因子比較敏感,尤其是各土地利用類型的生態(tài)赤字對(duì)當(dāng)量因子的準(zhǔn)確性依賴度較高。計(jì)算結(jié)果顯示:生態(tài)足跡的敏感性指數(shù)為0,生態(tài)承載力的敏感性指數(shù)為1,生態(tài)赤字的敏感性指數(shù)逐漸降低,由2000 年的1.78 下降至2018 年的0.26(圖4),即在當(dāng)量因子變化1%時(shí),生態(tài)足跡未發(fā)生變化,生態(tài)承載力增加或減少100%,生態(tài)赤字的變化量則呈遞減趨勢(shì),由2000 年的變化178%下降到2018 年的變化26%。表明生態(tài)足跡對(duì)當(dāng)量因子完全無彈性,生態(tài)承載力對(duì)當(dāng)量因子為單位彈性,生態(tài)赤字對(duì)當(dāng)量因子的敏感度在下降,在2000 年富有彈性,在2005年、2010年、2015年和2018年缺乏彈性。不難發(fā)現(xiàn),這一結(jié)果主要是由計(jì)算原理決定的,由式(4)(5)可知當(dāng)各土地利用類型的當(dāng)量因子調(diào)整一定比例時(shí),該地類的均衡因子并未改變,而產(chǎn)量因子卻發(fā)生相同比例變化,因而生態(tài)足跡值沒有變化,生態(tài)承載力同比例變化。由圖3可知流域生態(tài)足跡持續(xù)增加,生態(tài)承載力基本不變,二者差距持續(xù)拉大,生態(tài)承載力發(fā)生100%的變化對(duì)生態(tài)赤字的影響逐漸減弱,因此生態(tài)赤字的敏感性指數(shù)呈下降趨勢(shì)。

    圖4 2000—2018年瑪納斯流域生態(tài)足跡、承載力、赤字的敏感性指數(shù)分析

    對(duì)各地類生態(tài)赤字的敏感性指數(shù)進(jìn)行計(jì)算發(fā)現(xiàn),除草地的敏感性指數(shù)在2010年、2015年和2018年低于1 外,其他類型始終≥1,表明各土地利用類型的生態(tài)赤字對(duì)當(dāng)量因子富有彈性(圖5)。從變化趨勢(shì)上看,水域呈顯著上升趨勢(shì),由1.16上升至41.80;耕地呈緩慢增加趨勢(shì),由1.78增至2.02;林地基本不變,在均值1.37 附近波動(dòng),波動(dòng)范圍在-0.14~0.09;建設(shè)用地在均值3.62 附近呈波動(dòng)變化,波動(dòng)幅度為-40.5%~27.0%;草地的敏感性指數(shù)呈下降趨勢(shì),由3.03 下降至0.19;未利用地始終為1。由于能源用地只產(chǎn)生生態(tài)足跡,未利用地只提供生態(tài)承載力,合計(jì)流域總體生態(tài)盈虧時(shí),易將不同地類間生態(tài)盈余與赤字進(jìn)行抵消,因而導(dǎo)致總生態(tài)赤字對(duì)當(dāng)量因子的敏感性低于各地類生態(tài)赤字對(duì)當(dāng)量因子的敏感性。

    圖5 2000—2018年瑪納斯流域各地類生態(tài)赤字的敏感性指數(shù)分析

    4 結(jié)論與討論

    4.1 結(jié)論

    (1)瑪納斯河流域的均衡因子和產(chǎn)量因子均存在年際差異,均衡因子在2015 年達(dá)到最高值,在2000年處于最低值,產(chǎn)量因子的最大值出現(xiàn)在2005年,最小值出現(xiàn)在2015年,從多年平均值看,水域的均衡因子最高,達(dá)16.87,耕地的產(chǎn)量因子最高,為1.42,而未利用地的均衡因子和產(chǎn)量因子均最低,分別為0.30和0.45。

    (2)瑪納斯河流域生態(tài)赤字不斷擴(kuò)大,區(qū)域發(fā)展面臨的生態(tài)壓力持續(xù)增大,18 年間流域生態(tài)足跡增加了3.35 倍,而生態(tài)承載力無顯著變化,生態(tài)赤字年均增速為-0.47hm2/年。其中,生態(tài)足跡呈明顯的不均衡增長,草地和能源用地對(duì)流域生態(tài)足跡貢獻(xiàn)率最高,分別為54.71%和30.34%,二者的生態(tài)足跡增速同樣最快,分別為33.54×104hm2/年和9.94×104hm2/年。

    (3)瑪納斯河流域基于生態(tài)系統(tǒng)服務(wù)價(jià)值的生態(tài)足跡測算結(jié)果對(duì)流域單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子比較敏感,受計(jì)算方法影響,生態(tài)足跡敏感性指數(shù)為0,生態(tài)承載力敏感性指數(shù)為1,生態(tài)赤字的敏感性指數(shù)不斷下降,在各類型的生態(tài)赤字敏感性指數(shù)中,除草地敏感性指數(shù)下降至0.19,未利用地為1外,其余類型敏感性指數(shù)均大于1。因此研究區(qū)域生態(tài)服務(wù)價(jià)值當(dāng)量因子的準(zhǔn)確性對(duì)生態(tài)足跡測算結(jié)果尤其是在揭示地類差異方面起關(guān)鍵作用。

    4.2 討論

    4.2.1 與傳統(tǒng)生態(tài)足跡模型的比較

    將該研究測算的模型參數(shù)平均值與國家參數(shù)和省域參數(shù)相比,發(fā)現(xiàn)存在明顯差異(圖6)。其中國家參數(shù)來自于全球生態(tài)足跡網(wǎng)(Global Footprint Network)2017 年發(fā)布的《Working Guidebook to the National Footprint Accounts》數(shù)據(jù)[39],省域參數(shù)來自于劉某承等基于凈初級(jí)生產(chǎn)力對(duì)新疆均衡因子和產(chǎn)量因子的測算結(jié)果[38,40]。就均衡因子而言,水域的均衡因子高達(dá)16.87,遠(yuǎn)遠(yuǎn)高于國家參數(shù)和省域參數(shù);林地和草地的均衡因子均高于國家參數(shù)和省域參數(shù),分別為4.62和1.37;能源用地的均衡因子高于國家參數(shù)而低于省域參數(shù),為1.39;耕地和建設(shè)用地的均衡因子均低于國家參數(shù)和省域參數(shù),均為1.07;未利用地的均衡因子為0.30。該結(jié)果主要與測算方法有關(guān),基于生態(tài)系統(tǒng)服務(wù)價(jià)值獲得的均衡因子反映的是區(qū)域不同土地類型的生態(tài)服務(wù)能力差異,而傳統(tǒng)生態(tài)足跡模型的均衡因子反映的是區(qū)域不同土地類型的生物生產(chǎn)能力差異。干旱區(qū)內(nèi)陸河流域水域的生態(tài)服務(wù)功能顯著,而生物資源產(chǎn)品提供能力相對(duì)較弱,因而基于生態(tài)系統(tǒng)服務(wù)能力的均衡因子較高,而基于生物生產(chǎn)能力的均衡因子較低。傳統(tǒng)模型中國家參數(shù)的水域均衡因子約為該研究中水域均衡因子的1/48,這一數(shù)值與郭慧等[20]的研究結(jié)果相近。就產(chǎn)量因子而言,耕地、建設(shè)用地和水域的產(chǎn)量因子均高于國家參數(shù)和省域參數(shù),分別為1.42、1.42和1.31;林地的產(chǎn)量因子略高于省域參數(shù)而低于全國參數(shù),為1.04;草地的產(chǎn)量因子明顯高于省域參數(shù)而低于全國參數(shù),為0.95;未利用地的產(chǎn)量因子為0.45。基于生態(tài)系統(tǒng)服務(wù)價(jià)值獲得的產(chǎn)量因子中,耕地的產(chǎn)量因子與國家參數(shù)較接近,林地的產(chǎn)量因子與省域參數(shù)近似,表明耕地的生態(tài)服務(wù)價(jià)值與其生產(chǎn)力相關(guān),而林地的生態(tài)服務(wù)價(jià)值與其凈初級(jí)生產(chǎn)力密切相關(guān)。

    圖6 與傳統(tǒng)生態(tài)足跡模型因子對(duì)比

    基于生態(tài)系統(tǒng)服務(wù)價(jià)值的瑪納斯河流域生態(tài)足跡測算結(jié)果與國家參數(shù)和省域參數(shù)的測算結(jié)果相比,生態(tài)足跡顯著高于省域參數(shù)和國家參數(shù)的測算結(jié)果;生態(tài)承載力高于省域參數(shù),而與國家參數(shù)的測算結(jié)果近似,在2000 年、2005 年和2010 年略高于國家參數(shù)的測算結(jié)果,在2015 年和2018 年略低于國家參數(shù)的測算結(jié)果;生態(tài)赤字在2000年和2005年高于國家參數(shù)的測算結(jié)果而低于省域參數(shù)的測算結(jié)果,在2015年和2018 年略高于省域參數(shù)的測算結(jié)果而低于國家參數(shù)的測算結(jié)果(圖7)。從平均值上看,該研究測算的流域生態(tài)足跡是國家參數(shù)的1.60 倍,是省域參數(shù)的1.24 倍,生態(tài)承載力是國家參數(shù)的1.02 倍,是省域參數(shù)的2.20倍;生態(tài)赤字是國家參數(shù)的3.81倍,是省域參數(shù)的0.85倍。對(duì)比已有研究成果[41],該研究的測算結(jié)果符合流域生態(tài)服務(wù)功能的實(shí)際變化。

    圖7 2000—2018年不同參數(shù)下瑪納斯河流域生態(tài)足跡測算結(jié)果對(duì)比

    4.2.2 方法內(nèi)涵與結(jié)果指引

    相較于傳統(tǒng)的生態(tài)足跡模型,該文構(gòu)建的基于生態(tài)系統(tǒng)服務(wù)價(jià)值的生態(tài)足跡模型具有以下優(yōu)勢(shì)。

    (1)強(qiáng)調(diào)了生態(tài)系統(tǒng)服務(wù)功能的多樣性。傳統(tǒng)生態(tài)足跡模型的均衡因子和產(chǎn)量因子是由各地類生物生產(chǎn)性能力決定,僅考慮了土地的糧食生產(chǎn)和原材料供給的初級(jí)功能;而該研究模型中的均衡因子和產(chǎn)量因子是由生態(tài)系統(tǒng)服務(wù)價(jià)值決定,完整地考慮了不同土地利用類型的多種生態(tài)系統(tǒng)服務(wù)功能。

    (2)生態(tài)承載力核算內(nèi)容更加豐富準(zhǔn)確。傳統(tǒng)的生態(tài)足跡模型依據(jù)世界環(huán)境與發(fā)展委員會(huì)(WCED)的建議[42],生態(tài)承載力需扣除12%的生物多樣性保護(hù)面積,由于不同區(qū)域生境質(zhì)量差異明顯,這種“一刀切”的籠統(tǒng)做法,難以準(zhǔn)確反映區(qū)域?qū)嶋H承載能力,而該研究模型中已對(duì)生物多樣性保護(hù)功能進(jìn)行核算,因此無需再次扣除。此外,傳統(tǒng)生態(tài)足跡模型忽略了荒漠、沼澤等未利用地的生態(tài)供給能力,而對(duì)于西北干旱區(qū)而言,荒漠面積遼闊,其生態(tài)服務(wù)能力不容忽視,該研究模型將未利用地納入生態(tài)承載力的測算對(duì)象,使結(jié)果更加合理客觀。

    (3)模型參數(shù)本地化,增強(qiáng)了普適性。由于區(qū)域氣候、地形、生物生產(chǎn)力等不同,傳統(tǒng)生態(tài)足跡模型中的全球或全國的因子難以真實(shí)反映小尺度(特定區(qū)域)的生態(tài)狀況[43],而該研究模型參數(shù)是根據(jù)研究區(qū)各類生態(tài)系統(tǒng)的服務(wù)價(jià)值確定,有效解決了這一問題,同時(shí)也使得區(qū)域間的測算結(jié)果具有直接可比性。

    基于生態(tài)系統(tǒng)服務(wù)價(jià)值的生態(tài)足跡模型可量化人地系統(tǒng)間生態(tài)資產(chǎn)的供需情況和可持續(xù)程度,有助于揭示生態(tài)系統(tǒng)服務(wù)價(jià)值流動(dòng)規(guī)律及平衡機(jī)理,探索社會(huì)經(jīng)濟(jì)活動(dòng)與生態(tài)服務(wù)功能的耦合關(guān)聯(lián)機(jī)制,同時(shí)測算結(jié)果還可服務(wù)于生態(tài)補(bǔ)償、綠洲適度規(guī)模等相關(guān)研究。該文構(gòu)建的基于生態(tài)系統(tǒng)服務(wù)價(jià)值的生態(tài)足跡模型既反映了人類對(duì)于生物生產(chǎn)性土地的占用面積,也體現(xiàn)了區(qū)域發(fā)展對(duì)于生態(tài)系統(tǒng)服務(wù)功能的消耗狀況,能夠有效地反映區(qū)域的生態(tài)狀況。然而建設(shè)用地的模型參數(shù)沿用傳統(tǒng)生態(tài)足跡模型的處理方法,雖實(shí)現(xiàn)了對(duì)生態(tài)系統(tǒng)服務(wù)供給能力和人類對(duì)之占用程度的全面核算,但同時(shí)一定程度上也增加了測算結(jié)果的不確定性。因此,未來對(duì)于模型中建設(shè)用地的均衡因子與產(chǎn)量因子的確定有待進(jìn)一步深入研究。

    猜你喜歡
    瑪納斯當(dāng)量足跡
    成長足跡
    足跡
    心聲歌刊(2018年6期)2018-01-24 00:56:12
    民族史詩入夢(mèng)來——歌劇《瑪納斯》觀感
    歌劇(2017年11期)2018-01-23 03:41:04
    近60年瑪納斯河徑流變化特征分析
    瑪納斯河流域水文特征值及演變規(guī)律——以瑪納斯河為例
    足跡
    黃河之聲(2016年24期)2016-02-03 09:01:52
    談?wù)劇艾敿{斯”和“瑪納斯奇”這兩個(gè)詞
    語言與翻譯(2015年1期)2015-07-18 11:10:08
    春的足跡
    超壓測試方法對(duì)炸藥TNT當(dāng)量計(jì)算結(jié)果的影響

    夜夜爽天天搞| 男女视频在线观看网站免费 | 国内精品久久久久精免费| 少妇被粗大的猛进出69影院| 日韩精品免费视频一区二区三区| 国产日本99.免费观看| 丁香欧美五月| 五月伊人婷婷丁香| 欧美最黄视频在线播放免费| 国产精品精品国产色婷婷| 无限看片的www在线观看| 久久精品国产亚洲av香蕉五月| 三级男女做爰猛烈吃奶摸视频| 90打野战视频偷拍视频| 老汉色av国产亚洲站长工具| 国产亚洲精品av在线| 一夜夜www| 一进一出抽搐gif免费好疼| 1024手机看黄色片| 国产免费男女视频| 无遮挡黄片免费观看| 无遮挡黄片免费观看| 麻豆成人午夜福利视频| 成人一区二区视频在线观看| a级毛片a级免费在线| 国产黄a三级三级三级人| 在线观看舔阴道视频| 日日干狠狠操夜夜爽| 999久久久精品免费观看国产| 国产片内射在线| 国产亚洲av片在线观看秒播厂 | kizo精华| 亚洲人成网站在线观看播放| 18禁在线无遮挡免费观看视频| 三级男女做爰猛烈吃奶摸视频| 亚洲精品国产av成人精品| 日日摸夜夜添夜夜爱| 黄色日韩在线| 国产精品一区www在线观看| 日本色播在线视频| 真实男女啪啪啪动态图| 久久6这里有精品| 亚洲第一区二区三区不卡| 午夜爱爱视频在线播放| 久久精品夜夜夜夜夜久久蜜豆| 日日啪夜夜撸| 亚洲av成人精品一区久久| 少妇人妻精品综合一区二区 | 一边亲一边摸免费视频| 伊人久久精品亚洲午夜| 听说在线观看完整版免费高清| 欧美日本视频| 国模一区二区三区四区视频| 亚洲欧美日韩无卡精品| 乱码一卡2卡4卡精品| 日本三级黄在线观看| 国产黄a三级三级三级人| 国产免费男女视频| 亚洲av免费高清在线观看| 亚洲精品日韩av片在线观看| 日韩一区二区视频免费看| 久久亚洲国产成人精品v| 亚洲精品久久国产高清桃花| 日日啪夜夜撸| 日韩大尺度精品在线看网址| 一进一出抽搐gif免费好疼| 国内久久婷婷六月综合欲色啪| 国产精品,欧美在线| 成人一区二区视频在线观看| 日韩成人av中文字幕在线观看| 网址你懂的国产日韩在线| 春色校园在线视频观看| 丝袜喷水一区| 超碰av人人做人人爽久久| 啦啦啦啦在线视频资源| 久久久久久久久大av| 国产午夜精品论理片| 欧美成人一区二区免费高清观看| 日产精品乱码卡一卡2卡三| 又爽又黄a免费视频| 听说在线观看完整版免费高清| 男女边吃奶边做爰视频| 欧美激情久久久久久爽电影| 精品人妻熟女av久视频| 日本色播在线视频| 中文字幕久久专区| 亚洲av熟女| 国产乱人视频| 免费人成在线观看视频色| 免费不卡的大黄色大毛片视频在线观看 | av天堂在线播放| 成人毛片60女人毛片免费| 国产精品久久久久久精品电影| 国产精品野战在线观看| 99久国产av精品| 秋霞在线观看毛片| 久久午夜亚洲精品久久| 国产在线精品亚洲第一网站| 在线观看免费视频日本深夜| 少妇高潮的动态图| 久久精品久久久久久噜噜老黄 | 欧美色视频一区免费| 日韩高清综合在线| 最近手机中文字幕大全| 欧美成人一区二区免费高清观看| 99riav亚洲国产免费| 亚洲四区av| 国产精品久久久久久精品电影小说 | 色尼玛亚洲综合影院| 在线观看66精品国产| 国产蜜桃级精品一区二区三区| 岛国在线免费视频观看| 国内精品久久久久精免费| 日本免费一区二区三区高清不卡| 97在线视频观看| 日韩一区二区三区影片| 亚洲国产精品合色在线| 日韩av不卡免费在线播放| 日本黄大片高清| 有码 亚洲区| 热99re8久久精品国产| 色5月婷婷丁香| 男人舔女人下体高潮全视频| 国产 一区 欧美 日韩| 亚洲色图av天堂| 性欧美人与动物交配| 伦精品一区二区三区| 国产中年淑女户外野战色| 亚洲av中文字字幕乱码综合| 国产午夜精品久久久久久一区二区三区| 啦啦啦啦在线视频资源| 婷婷精品国产亚洲av| 蜜桃亚洲精品一区二区三区| av福利片在线观看| 别揉我奶头 嗯啊视频| .国产精品久久| 亚洲精品色激情综合| 久久久久久久亚洲中文字幕| 亚洲精品自拍成人| 99精品在免费线老司机午夜| 22中文网久久字幕| 久久精品国产亚洲av涩爱 | 在线免费观看不下载黄p国产| 搡女人真爽免费视频火全软件| 日本av手机在线免费观看| 国产三级在线视频| 搡老妇女老女人老熟妇| 又黄又爽又刺激的免费视频.| 欧美一区二区亚洲| 别揉我奶头 嗯啊视频| 国模一区二区三区四区视频| 在线播放国产精品三级| а√天堂www在线а√下载| 极品教师在线视频| 亚洲人成网站在线播放欧美日韩| 国产日本99.免费观看| 欧美人与善性xxx| 国产老妇伦熟女老妇高清| 成人综合一区亚洲| 亚洲精品成人久久久久久| 啦啦啦韩国在线观看视频| 久久久午夜欧美精品| 欧美潮喷喷水| 自拍偷自拍亚洲精品老妇| 免费一级毛片在线播放高清视频| 中文字幕av成人在线电影| 国产黄片视频在线免费观看| 黄片wwwwww| 高清毛片免费观看视频网站| 色哟哟·www| 日韩 亚洲 欧美在线| 少妇的逼好多水| 日韩 亚洲 欧美在线| 看免费成人av毛片| 如何舔出高潮| 国产精品久久久久久久电影| 国产黄色小视频在线观看| 2022亚洲国产成人精品| 好男人视频免费观看在线| 中文字幕久久专区| 草草在线视频免费看| 久久久久国产网址| 亚洲自偷自拍三级| 国产91av在线免费观看| 国产精品久久久久久亚洲av鲁大| 十八禁国产超污无遮挡网站| 人人妻人人澡欧美一区二区| 亚洲四区av| 国产午夜精品论理片| 色播亚洲综合网| 欧美一级a爱片免费观看看| 2021天堂中文幕一二区在线观| 国产激情偷乱视频一区二区| 亚洲欧美精品专区久久| 亚洲国产精品成人久久小说 | 免费看美女性在线毛片视频| 亚洲精品亚洲一区二区| 欧美xxxx黑人xx丫x性爽| 天堂av国产一区二区熟女人妻| 亚洲乱码一区二区免费版| 97超碰精品成人国产| 国产精品无大码| 久久国内精品自在自线图片| 99久久九九国产精品国产免费| 毛片女人毛片| 欧美不卡视频在线免费观看| 狂野欧美激情性xxxx在线观看| 极品教师在线视频| 变态另类丝袜制服| 亚洲第一电影网av| 日韩欧美精品v在线| 伦精品一区二区三区| 午夜a级毛片| 人妻少妇偷人精品九色| 三级男女做爰猛烈吃奶摸视频| 国产黄a三级三级三级人| 免费av观看视频| 菩萨蛮人人尽说江南好唐韦庄 | 国产老妇女一区| 长腿黑丝高跟| 99久久成人亚洲精品观看| 2022亚洲国产成人精品| 午夜福利视频1000在线观看| 亚洲成人中文字幕在线播放| 禁无遮挡网站| 伦精品一区二区三区| 国产精品人妻久久久影院| 国产毛片a区久久久久| 亚州av有码| 91精品国产九色| 国产免费一级a男人的天堂| 免费一级毛片在线播放高清视频| av国产免费在线观看| 在线免费十八禁| 一个人看的www免费观看视频| 免费观看的影片在线观看| 人人妻人人看人人澡| 精品一区二区免费观看| 国产精品久久久久久久久免| 亚洲精品久久久久久婷婷小说 | 老师上课跳d突然被开到最大视频| 日韩一本色道免费dvd| 久久久欧美国产精品| 国产真实伦视频高清在线观看| 能在线免费看毛片的网站| 免费观看a级毛片全部| 欧美成人一区二区免费高清观看| 国产成年人精品一区二区| 国产中年淑女户外野战色| 国产熟女欧美一区二区| 国产伦精品一区二区三区四那| 久久99精品国语久久久| 日本一二三区视频观看| av又黄又爽大尺度在线免费看 | 国产av麻豆久久久久久久| 精品久久久久久久久亚洲| 少妇熟女aⅴ在线视频| 99热只有精品国产| 久久人人精品亚洲av| 夜夜爽天天搞| 亚洲国产精品合色在线| 欧美日韩乱码在线| 欧洲精品卡2卡3卡4卡5卡区| av在线蜜桃| 欧美在线一区亚洲| 国产精品免费一区二区三区在线| 国产精品美女特级片免费视频播放器| 看免费成人av毛片| 日本免费a在线| 久久国产乱子免费精品| 乱系列少妇在线播放| 在线观看免费视频日本深夜| 国产色婷婷99| 亚洲美女视频黄频| 中出人妻视频一区二区| 亚洲成人精品中文字幕电影| 夫妻性生交免费视频一级片| 一个人观看的视频www高清免费观看| 日本色播在线视频| 亚洲欧美中文字幕日韩二区| 免费观看在线日韩| 插阴视频在线观看视频| 精品久久国产蜜桃| 欧美一级a爱片免费观看看| 99久久精品热视频| 亚洲精品乱码久久久久久按摩| 最近中文字幕高清免费大全6| 美女国产视频在线观看| 99热全是精品| 日韩人妻高清精品专区| 夫妻性生交免费视频一级片| 青春草亚洲视频在线观看| 亚洲三级黄色毛片| 人妻久久中文字幕网| 久久久久久九九精品二区国产| 亚洲人成网站在线观看播放| 国产精品国产三级国产av玫瑰| 免费在线观看成人毛片| 国产亚洲精品久久久com| 国产探花极品一区二区| 久久久久久大精品| 悠悠久久av| 一区二区三区免费毛片| 亚洲国产欧美人成| 中文在线观看免费www的网站| 人人妻人人澡欧美一区二区| 国产白丝娇喘喷水9色精品| 色综合站精品国产| 精品久久久久久久人妻蜜臀av| 一级av片app| 三级国产精品欧美在线观看| 国产亚洲av嫩草精品影院| 久久精品综合一区二区三区| 国产精品麻豆人妻色哟哟久久 | 久久久久久久午夜电影| 啦啦啦观看免费观看视频高清| 麻豆av噜噜一区二区三区| 国产午夜福利久久久久久| 尾随美女入室| 国产黄a三级三级三级人| 国产在线精品亚洲第一网站| 晚上一个人看的免费电影| 悠悠久久av| 午夜福利视频1000在线观看| 欧美性感艳星| 国产在线精品亚洲第一网站| 亚洲自偷自拍三级| 少妇被粗大猛烈的视频| 久久99热6这里只有精品| 久久久久久九九精品二区国产| 丝袜喷水一区| 久久婷婷人人爽人人干人人爱| 亚洲内射少妇av| 日本三级黄在线观看| 国产真实乱freesex| 久久久久久大精品| 国产成人91sexporn| 麻豆av噜噜一区二区三区| 午夜福利在线观看吧| 综合色av麻豆| 久久久久久伊人网av| 人妻夜夜爽99麻豆av| 亚洲精品影视一区二区三区av| 深夜a级毛片| 午夜福利在线观看免费完整高清在 | 国产成人影院久久av| 乱人视频在线观看| av专区在线播放| 亚洲av成人av| 成人一区二区视频在线观看| 99久久成人亚洲精品观看| 精品人妻熟女av久视频| 日本黄大片高清| 欧美日韩国产亚洲二区| 亚洲在线观看片| 九九爱精品视频在线观看| 熟女人妻精品中文字幕| 欧美激情在线99| 精品久久久久久久久av| 国产成人a区在线观看| 校园春色视频在线观看| 乱系列少妇在线播放| av在线播放精品| 欧美高清成人免费视频www| 成熟少妇高潮喷水视频| 嫩草影院入口| 成人特级av手机在线观看| 热99在线观看视频| 久久热精品热| 欧美色视频一区免费| 18禁在线无遮挡免费观看视频| 亚洲av男天堂| 国产真实伦视频高清在线观看| 久久午夜亚洲精品久久| 91久久精品国产一区二区成人| 国产一区亚洲一区在线观看| 乱人视频在线观看| 91aial.com中文字幕在线观看| 伦理电影大哥的女人| 啦啦啦韩国在线观看视频| 日日撸夜夜添| 天堂中文最新版在线下载 | 亚洲色图av天堂| 超碰av人人做人人爽久久| 美女大奶头视频| 国产成人午夜福利电影在线观看| 久久精品久久久久久久性| 一边摸一边抽搐一进一小说| 少妇高潮的动态图| 亚洲性久久影院| 丝袜喷水一区| 国产精品免费一区二区三区在线| 在线观看av片永久免费下载| 熟女人妻精品中文字幕| 欧美日本亚洲视频在线播放| 中出人妻视频一区二区| 狂野欧美白嫩少妇大欣赏| 身体一侧抽搐| 网址你懂的国产日韩在线| 亚洲精华国产精华液的使用体验 | 亚洲精品久久国产高清桃花| 国产精品国产高清国产av| 黄色欧美视频在线观看| 99久久成人亚洲精品观看| 免费看日本二区| 免费看a级黄色片| 久久精品91蜜桃| 成年版毛片免费区| av视频在线观看入口| 丰满乱子伦码专区| 亚洲av男天堂| 午夜久久久久精精品| 久久精品久久久久久久性| 欧美色视频一区免费| 欧美激情在线99| 国产精品爽爽va在线观看网站| 国产久久久一区二区三区| 91麻豆精品激情在线观看国产| 91av网一区二区| 国产女主播在线喷水免费视频网站 | 亚洲一级一片aⅴ在线观看| 国产三级在线视频| 99热这里只有是精品在线观看| 插阴视频在线观看视频| 免费搜索国产男女视频| or卡值多少钱| 国产成人午夜福利电影在线观看| 国产色婷婷99| 一边亲一边摸免费视频| 亚洲18禁久久av| 国产日韩欧美在线精品| 人妻久久中文字幕网| 如何舔出高潮| 精品久久久久久久久亚洲| 国产精品福利在线免费观看| 欧美一区二区精品小视频在线| 欧美成人一区二区免费高清观看| 少妇裸体淫交视频免费看高清| 菩萨蛮人人尽说江南好唐韦庄 | 熟女人妻精品中文字幕| 一本一本综合久久| 欧美成人一区二区免费高清观看| 超碰av人人做人人爽久久| 久久久久免费精品人妻一区二区| 久久久久久久亚洲中文字幕| 人妻系列 视频| 午夜免费激情av| 熟女电影av网| 婷婷精品国产亚洲av| 久久精品国产清高在天天线| 精品久久久久久久末码| 国产老妇女一区| 麻豆成人午夜福利视频| 国产午夜精品久久久久久一区二区三区| 国产麻豆成人av免费视频| 久久99蜜桃精品久久| 亚洲久久久久久中文字幕| 国产精品国产三级国产av玫瑰| 亚洲av成人精品一区久久| 欧美一区二区国产精品久久精品| 久久久久免费精品人妻一区二区| 日日干狠狠操夜夜爽| 成人亚洲欧美一区二区av| 精品一区二区三区视频在线| 亚洲自拍偷在线| 国内久久婷婷六月综合欲色啪| 免费电影在线观看免费观看| 欧美潮喷喷水| 91在线精品国自产拍蜜月| av福利片在线观看| 国产精品一区www在线观看| 国产午夜精品久久久久久一区二区三区| 国产一级毛片在线| 久久精品国产清高在天天线| 精品久久久久久久人妻蜜臀av| 99九九线精品视频在线观看视频| 伊人久久精品亚洲午夜| 99久久精品国产国产毛片| 波多野结衣高清无吗| 人妻久久中文字幕网| 黑人高潮一二区| 小蜜桃在线观看免费完整版高清| 午夜亚洲福利在线播放| 国产精品久久久久久久电影| 久久精品国产亚洲网站| 久久这里只有精品中国| av专区在线播放| 国产黄片视频在线免费观看| 免费观看在线日韩| 麻豆av噜噜一区二区三区| www.色视频.com| 99热精品在线国产| 国产黄片视频在线免费观看| 永久网站在线| 亚洲成a人片在线一区二区| 精品国产三级普通话版| 欧美潮喷喷水| 69av精品久久久久久| 色哟哟哟哟哟哟| 成人鲁丝片一二三区免费| 草草在线视频免费看| 成人一区二区视频在线观看| 91在线精品国自产拍蜜月| 午夜福利在线观看免费完整高清在 | 久久久久久久久久成人| 男插女下体视频免费在线播放| 1000部很黄的大片| 国产精品久久久久久久久免| 国产高清视频在线观看网站| 搞女人的毛片| 久久九九热精品免费| 人人妻人人澡欧美一区二区| 乱码一卡2卡4卡精品| 中文字幕熟女人妻在线| 小说图片视频综合网站| 少妇熟女欧美另类| 99热网站在线观看| 日韩在线高清观看一区二区三区| 免费搜索国产男女视频| 国产精品综合久久久久久久免费| 免费看日本二区| av又黄又爽大尺度在线免费看 | h日本视频在线播放| 国产人妻一区二区三区在| 久久国产乱子免费精品| 亚洲欧美日韩高清在线视频| 日本免费a在线| 日本-黄色视频高清免费观看| 国产麻豆成人av免费视频| 在线观看一区二区三区| 国产伦精品一区二区三区四那| 一夜夜www| 国产一级毛片七仙女欲春2| 国产单亲对白刺激| 国产男人的电影天堂91| 成人高潮视频无遮挡免费网站| 18禁裸乳无遮挡免费网站照片| 在线观看av片永久免费下载| 亚洲精品乱码久久久v下载方式| 精品午夜福利在线看| 伊人久久精品亚洲午夜| 国产淫片久久久久久久久| 国产不卡一卡二| 青春草视频在线免费观看| 欧美zozozo另类| 两个人的视频大全免费| 精品久久久久久成人av| 欧美性猛交黑人性爽| av天堂中文字幕网| 麻豆精品久久久久久蜜桃| 国产精品,欧美在线| 成人永久免费在线观看视频| 在现免费观看毛片| 精华霜和精华液先用哪个| 啦啦啦韩国在线观看视频| 日韩av在线大香蕉| 日韩av不卡免费在线播放| 高清在线视频一区二区三区 | 91久久精品国产一区二区成人| 日本黄大片高清| 特大巨黑吊av在线直播| 国产精品三级大全| 亚洲色图av天堂| 热99re8久久精品国产| 久久久久久伊人网av| 欧美日韩国产亚洲二区| 日韩亚洲欧美综合| 国产不卡一卡二| 亚洲人成网站在线观看播放| 18禁裸乳无遮挡免费网站照片| 欧美又色又爽又黄视频| 日韩欧美国产在线观看| 少妇裸体淫交视频免费看高清| 你懂的网址亚洲精品在线观看 | 欧美激情国产日韩精品一区| 精品一区二区三区视频在线| 日韩精品有码人妻一区| 18+在线观看网站| 国产精华一区二区三区| www日本黄色视频网| 精品久久久久久久久久免费视频| 欧美精品国产亚洲| 国产成人午夜福利电影在线观看| 久久婷婷人人爽人人干人人爱| 国产一区二区在线av高清观看| 亚洲国产精品久久男人天堂| 丰满的人妻完整版| 国产探花在线观看一区二区| 成年女人永久免费观看视频| 国产在线男女| 国产黄a三级三级三级人| 国产成人午夜福利电影在线观看| 青春草亚洲视频在线观看| 日韩,欧美,国产一区二区三区 | 有码 亚洲区| av在线亚洲专区| 亚洲婷婷狠狠爱综合网| 国产成人精品久久久久久| 热99re8久久精品国产| 日韩中字成人| 国产成人a区在线观看| 国产精品爽爽va在线观看网站| 尾随美女入室| 精品久久久久久久末码| 12—13女人毛片做爰片一| 久久99热6这里只有精品| 热99在线观看视频| 青春草视频在线免费观看| 免费看美女性在线毛片视频| 国语自产精品视频在线第100页| 一区二区三区免费毛片| 一区二区三区高清视频在线| 久久精品久久久久久久性| 国产伦精品一区二区三区视频9| 最近最新中文字幕大全电影3|