• <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é)果的影響

    一区二区三区高清视频在线| 后天国语完整版免费观看| 日韩欧美国产一区二区入口| 久久国产精品人妻蜜桃| 如日韩欧美国产精品一区二区三区| 99国产综合亚洲精品| 久久中文看片网| 97碰自拍视频| 亚洲一区二区三区不卡视频| 免费不卡黄色视频| 一个人观看的视频www高清免费观看 | 丰满的人妻完整版| 亚洲精品中文字幕在线视频| 韩国av一区二区三区四区| 日韩大尺度精品在线看网址 | 国产三级在线视频| 国产亚洲精品久久久久久毛片| 啦啦啦观看免费观看视频高清 | 纯流量卡能插随身wifi吗| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美日韩亚洲国产一区二区在线观看| 成人手机av| 香蕉丝袜av| 久久中文字幕人妻熟女| 亚洲成av人片免费观看| 久久久久国产一级毛片高清牌| 亚洲情色 制服丝袜| 精品久久久久久久人妻蜜臀av | 久久精品aⅴ一区二区三区四区| 一个人免费在线观看的高清视频| 国产欧美日韩一区二区三| www.999成人在线观看| 成人亚洲精品一区在线观看| 久久人妻福利社区极品人妻图片| 午夜福利高清视频| 午夜福利欧美成人| 亚洲精品久久成人aⅴ小说| 国产成人系列免费观看| 亚洲全国av大片| 精品一区二区三区视频在线观看免费| 乱人伦中国视频| 亚洲精品一区av在线观看| 日韩欧美一区视频在线观看| 亚洲激情在线av| 男女做爰动态图高潮gif福利片 | 日韩欧美一区二区三区在线观看| 韩国精品一区二区三区| 黄色成人免费大全| 波多野结衣巨乳人妻| 亚洲av成人一区二区三| 在线国产一区二区在线| 丝袜美足系列| 啦啦啦观看免费观看视频高清 | 国产精品香港三级国产av潘金莲| 国产成人欧美| 变态另类成人亚洲欧美熟女 | 天堂√8在线中文| 亚洲熟妇熟女久久| 美女高潮喷水抽搐中文字幕| 免费少妇av软件| 黄色 视频免费看| 无限看片的www在线观看| 成人18禁高潮啪啪吃奶动态图| 狂野欧美激情性xxxx| 中文字幕人妻熟女乱码| 国产伦人伦偷精品视频| 亚洲av成人av| 国产高清激情床上av| 亚洲少妇的诱惑av| 国产一级毛片七仙女欲春2 | 日韩欧美国产在线观看| 无人区码免费观看不卡| 电影成人av| 国产黄a三级三级三级人| 一级a爱片免费观看的视频| 国产欧美日韩一区二区三区在线| 久久精品亚洲熟妇少妇任你| 人人澡人人妻人| 欧美另类亚洲清纯唯美| 日本一区二区免费在线视频| 久久亚洲真实| 国产三级黄色录像| 国产私拍福利视频在线观看| 久久人人精品亚洲av| 麻豆久久精品国产亚洲av| 欧美性长视频在线观看| 精品国内亚洲2022精品成人| 精品欧美一区二区三区在线| 久久午夜亚洲精品久久| 欧美成人性av电影在线观看| 女生性感内裤真人,穿戴方法视频| 免费少妇av软件| 午夜福利高清视频| 精品不卡国产一区二区三区| 在线观看舔阴道视频| 九色国产91popny在线| 99久久综合精品五月天人人| 在线观看日韩欧美| 18禁裸乳无遮挡免费网站照片 | 一本综合久久免费| 精品国内亚洲2022精品成人| 国产精品 欧美亚洲| 桃色一区二区三区在线观看| av福利片在线| 成年人黄色毛片网站| 国产精品 国内视频| 欧美激情高清一区二区三区| 88av欧美| 少妇被粗大的猛进出69影院| 午夜影院日韩av| 亚洲少妇的诱惑av| 日本免费一区二区三区高清不卡 | 一个人观看的视频www高清免费观看 | 亚洲性夜色夜夜综合| 十八禁网站免费在线| 亚洲久久久国产精品| 日韩 欧美 亚洲 中文字幕| 久久国产乱子伦精品免费另类| 亚洲人成电影观看| 自线自在国产av| 12—13女人毛片做爰片一| 国产视频一区二区在线看| 一级a爱视频在线免费观看| 久久精品亚洲熟妇少妇任你| 好看av亚洲va欧美ⅴa在| 俄罗斯特黄特色一大片| 久久精品国产亚洲av香蕉五月| 亚洲精品久久成人aⅴ小说| 亚洲国产精品999在线| 午夜福利18| 乱人伦中国视频| 乱人伦中国视频| 好男人电影高清在线观看| 成人av一区二区三区在线看| 免费在线观看影片大全网站| 动漫黄色视频在线观看| 少妇 在线观看| 一二三四在线观看免费中文在| 成年版毛片免费区| 女性被躁到高潮视频| 亚洲最大成人中文| 欧美绝顶高潮抽搐喷水| 亚洲国产高清在线一区二区三 | 无人区码免费观看不卡| 一级作爱视频免费观看| 亚洲色图 男人天堂 中文字幕| 国产成人影院久久av| 又大又爽又粗| 黄片大片在线免费观看| 高清毛片免费观看视频网站| 嫁个100分男人电影在线观看| 在线观看午夜福利视频| 日韩高清综合在线| tocl精华| 女警被强在线播放| 一级作爱视频免费观看| 欧美中文综合在线视频| 国产片内射在线| 国产精品一区二区在线不卡| 久久国产精品人妻蜜桃| 高潮久久久久久久久久久不卡| 黄色a级毛片大全视频| 国产99久久九九免费精品| 丝袜美足系列| 国产1区2区3区精品| 精品久久久久久久人妻蜜臀av | 国产精品98久久久久久宅男小说| 欧美日韩精品网址| 欧美精品亚洲一区二区| 亚洲欧美日韩另类电影网站| 不卡一级毛片| 国产欧美日韩一区二区三区在线| 亚洲片人在线观看| 久久精品国产亚洲av高清一级| 女人精品久久久久毛片| 窝窝影院91人妻| 国产乱人伦免费视频| 精品电影一区二区在线| 一区福利在线观看| 99在线人妻在线中文字幕| 18禁裸乳无遮挡免费网站照片 | 久久中文看片网| 90打野战视频偷拍视频| 18禁美女被吸乳视频| 别揉我奶头~嗯~啊~动态视频| 精品无人区乱码1区二区| 一卡2卡三卡四卡精品乱码亚洲| 国产成人免费无遮挡视频| 老司机深夜福利视频在线观看| 久久香蕉精品热| 老熟妇仑乱视频hdxx| 亚洲第一av免费看| 亚洲无线在线观看| 亚洲黑人精品在线| 男人的好看免费观看在线视频 | 精品国产超薄肉色丝袜足j| 国产欧美日韩一区二区三区在线| 国产av在哪里看| 国产xxxxx性猛交| 国产精品香港三级国产av潘金莲| 两个人免费观看高清视频| 90打野战视频偷拍视频| 男女下面插进去视频免费观看| 91国产中文字幕| 国产亚洲欧美精品永久| 色播在线永久视频| 国产精品一区二区精品视频观看| 黑人巨大精品欧美一区二区蜜桃| 久久婷婷人人爽人人干人人爱 | 波多野结衣巨乳人妻| 欧美不卡视频在线免费观看 | 嫁个100分男人电影在线观看| 两人在一起打扑克的视频| 国产亚洲av高清不卡| 波多野结衣高清无吗| 国产在线观看jvid| 波多野结衣高清无吗| 1024香蕉在线观看| 丝袜美腿诱惑在线| 18美女黄网站色大片免费观看| 亚洲人成电影免费在线| 国产精品 欧美亚洲| 国产99久久九九免费精品| 一级a爱视频在线免费观看| 亚洲全国av大片| 国产精品香港三级国产av潘金莲| 亚洲精品粉嫩美女一区| 天天添夜夜摸| 99香蕉大伊视频| 精品久久久久久成人av| 久久青草综合色| 不卡av一区二区三区| 国产精品99久久99久久久不卡| 精品久久久久久,| 欧美日本中文国产一区发布| 久久精品国产亚洲av高清一级| 欧美乱妇无乱码| 美女扒开内裤让男人捅视频| 中文字幕高清在线视频| 亚洲熟女毛片儿| 亚洲伊人色综图| 久久天堂一区二区三区四区| 亚洲av第一区精品v没综合| 国产成人一区二区三区免费视频网站| 一级毛片精品| 好看av亚洲va欧美ⅴa在| www.自偷自拍.com| 岛国视频午夜一区免费看| 狠狠狠狠99中文字幕| 别揉我奶头~嗯~啊~动态视频| 国产精品 欧美亚洲| 色播在线永久视频| 在线观看www视频免费| 美女午夜性视频免费| 十八禁网站免费在线| 国产亚洲精品久久久久5区| 亚洲精品国产精品久久久不卡| 免费看美女性在线毛片视频| 中文字幕人妻熟女乱码| 亚洲免费av在线视频| 18禁黄网站禁片午夜丰满| 性少妇av在线| 在线av久久热| 精品国内亚洲2022精品成人| 亚洲成人免费电影在线观看| 在线视频色国产色| 啦啦啦 在线观看视频| 在线观看午夜福利视频| 久久草成人影院| 久久 成人 亚洲| 国产精品一区二区精品视频观看| 一进一出抽搐gif免费好疼| 国产精品免费一区二区三区在线| 久久人妻熟女aⅴ| 别揉我奶头~嗯~啊~动态视频| 精品日产1卡2卡| 午夜日韩欧美国产| 一二三四在线观看免费中文在| 久久久久亚洲av毛片大全| 国产欧美日韩精品亚洲av| 久久久久久久精品吃奶| 国产国语露脸激情在线看| 亚洲在线自拍视频| 在线十欧美十亚洲十日本专区| 精品电影一区二区在线| 免费在线观看完整版高清| 国产成人精品无人区| 黄色片一级片一级黄色片| 国产在线精品亚洲第一网站| 国产乱人伦免费视频| 亚洲成人精品中文字幕电影| 亚洲熟妇熟女久久| 久久精品91无色码中文字幕| 婷婷丁香在线五月| 欧美乱码精品一区二区三区| 免费在线观看影片大全网站| 久热这里只有精品99| 久久久久久久精品吃奶| 99久久精品国产亚洲精品| 好看av亚洲va欧美ⅴa在| 老司机靠b影院| 国产精品二区激情视频| 夜夜爽天天搞| 国产精品亚洲av一区麻豆| 成人精品一区二区免费| 国产欧美日韩一区二区三区在线| 国产精品,欧美在线| 日本a在线网址| 精品一区二区三区四区五区乱码| 怎么达到女性高潮| 久久国产精品影院| av欧美777| 亚洲av第一区精品v没综合| 首页视频小说图片口味搜索| 日日摸夜夜添夜夜添小说| 亚洲一区二区三区色噜噜| 精品人妻在线不人妻| a级毛片在线看网站| 可以免费在线观看a视频的电影网站| 免费一级毛片在线播放高清视频 | 亚洲黑人精品在线| 免费在线观看亚洲国产| 国产成人免费无遮挡视频| 亚洲国产欧美日韩在线播放| 精品国产美女av久久久久小说| 国产成人精品无人区| 欧美丝袜亚洲另类 | 88av欧美| 别揉我奶头~嗯~啊~动态视频| 精品卡一卡二卡四卡免费| 在线观看免费视频网站a站| 热99re8久久精品国产| 欧美精品啪啪一区二区三区| 少妇被粗大的猛进出69影院| 成年人黄色毛片网站| 两性午夜刺激爽爽歪歪视频在线观看 | 夜夜爽天天搞| 国产99久久九九免费精品| 免费少妇av软件| 99国产精品一区二区蜜桃av| 亚洲精品国产一区二区精华液| 动漫黄色视频在线观看| 午夜福利免费观看在线| 国产97色在线日韩免费| 岛国在线观看网站| 国产蜜桃级精品一区二区三区| 给我免费播放毛片高清在线观看| 亚洲天堂国产精品一区在线| 国产精品综合久久久久久久免费 | 亚洲国产欧美一区二区综合| 成熟少妇高潮喷水视频| 精品福利观看| 日韩 欧美 亚洲 中文字幕| 免费看美女性在线毛片视频| 99久久精品国产亚洲精品| 丝袜美足系列| 天天躁夜夜躁狠狠躁躁| 91av网站免费观看| 国产精品久久电影中文字幕| 美女国产高潮福利片在线看| 天天躁狠狠躁夜夜躁狠狠躁| 日本免费a在线| 国产97色在线日韩免费| 91在线观看av| 久久久久久人人人人人| 免费看美女性在线毛片视频| 欧洲精品卡2卡3卡4卡5卡区| 999精品在线视频| 黄色视频不卡| 人人妻人人澡人人看| netflix在线观看网站| 天天躁狠狠躁夜夜躁狠狠躁| 俄罗斯特黄特色一大片| 欧美另类亚洲清纯唯美| 日本免费一区二区三区高清不卡 | 精品高清国产在线一区| 久久国产精品男人的天堂亚洲| 嫩草影院精品99| 老司机午夜福利在线观看视频| 日韩精品中文字幕看吧| av视频免费观看在线观看| 欧美色视频一区免费| 一区二区三区国产精品乱码| 亚洲国产毛片av蜜桃av| 久久中文字幕人妻熟女| 一卡2卡三卡四卡精品乱码亚洲| 久久精品国产综合久久久| 欧美日本亚洲视频在线播放| 男女之事视频高清在线观看| 国产激情久久老熟女| 欧美乱色亚洲激情| 99国产极品粉嫩在线观看| 午夜亚洲福利在线播放| 法律面前人人平等表现在哪些方面| 又紧又爽又黄一区二区| 男女床上黄色一级片免费看| 91麻豆av在线| 老司机靠b影院| 国内毛片毛片毛片毛片毛片| www国产在线视频色| 国产免费男女视频| 国产91精品成人一区二区三区| www.www免费av| 50天的宝宝边吃奶边哭怎么回事| 国产欧美日韩一区二区三区在线| 日本a在线网址| 精品日产1卡2卡| 午夜两性在线视频| 亚洲精品粉嫩美女一区| 欧美丝袜亚洲另类 | 亚洲片人在线观看| 嫁个100分男人电影在线观看| 91老司机精品| 国产av在哪里看| 国产成人精品久久二区二区91| 欧美另类亚洲清纯唯美| 久热爱精品视频在线9| 欧美黄色片欧美黄色片| 在线观看免费午夜福利视频| 精品国产国语对白av| 久久久国产欧美日韩av| 日本a在线网址| 亚洲精品国产色婷婷电影| 亚洲精品在线美女| 国产精品九九99| 两人在一起打扑克的视频| 老司机在亚洲福利影院| 88av欧美| 国产精品秋霞免费鲁丝片| 熟妇人妻久久中文字幕3abv| avwww免费| 91国产中文字幕| 免费一级毛片在线播放高清视频 | 日韩欧美一区二区三区在线观看| 国产蜜桃级精品一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 男男h啪啪无遮挡| 日韩成人在线观看一区二区三区| 午夜福利在线观看吧| 看免费av毛片| 男女午夜视频在线观看| 亚洲黑人精品在线| 国产精品一区二区在线不卡| 自线自在国产av| 欧美成人性av电影在线观看| 九色国产91popny在线| 韩国精品一区二区三区| 久久久久久久精品吃奶| 91老司机精品| 美女大奶头视频| 日韩高清综合在线| 真人做人爱边吃奶动态| 亚洲欧美日韩无卡精品| 国产成人精品无人区| 黄色a级毛片大全视频| av中文乱码字幕在线| 母亲3免费完整高清在线观看| 啦啦啦观看免费观看视频高清 | 1024视频免费在线观看| 中文字幕av电影在线播放| 亚洲欧美日韩无卡精品| 久久久久九九精品影院| 91在线观看av| 琪琪午夜伦伦电影理论片6080| www.精华液| 日韩免费av在线播放| 亚洲国产精品成人综合色| 免费观看精品视频网站| 成人三级黄色视频| 给我免费播放毛片高清在线观看| 在线观看舔阴道视频| 男女做爰动态图高潮gif福利片 | 99久久综合精品五月天人人| 一区二区三区激情视频| 黄色视频,在线免费观看| 国产精品久久电影中文字幕| 亚洲精品在线美女| 黑丝袜美女国产一区| 校园春色视频在线观看| 欧美亚洲日本最大视频资源| 欧美成人性av电影在线观看| 99精品在免费线老司机午夜| 国产又色又爽无遮挡免费看| 国产成人av教育| 成年版毛片免费区| 国产亚洲欧美在线一区二区| 日本在线视频免费播放| 人人妻,人人澡人人爽秒播| 亚洲黑人精品在线| 一级毛片女人18水好多| 美女免费视频网站| 欧美日韩福利视频一区二区| 亚洲精品一卡2卡三卡4卡5卡| 大陆偷拍与自拍| 亚洲天堂国产精品一区在线| 日韩精品青青久久久久久| 男女下面插进去视频免费观看| 精品国产乱子伦一区二区三区| 搡老岳熟女国产| 18禁黄网站禁片午夜丰满| 97人妻天天添夜夜摸| 91麻豆精品激情在线观看国产| 这个男人来自地球电影免费观看| 欧美另类亚洲清纯唯美| 亚洲午夜精品一区,二区,三区| 美女午夜性视频免费| 欧美成人性av电影在线观看| 成熟少妇高潮喷水视频| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利18| 免费女性裸体啪啪无遮挡网站| 老司机午夜十八禁免费视频| 国产精品野战在线观看| 老熟妇仑乱视频hdxx| 免费高清在线观看日韩| 人人妻人人澡人人看| 亚洲av片天天在线观看| 午夜福利欧美成人| 免费人成视频x8x8入口观看| 校园春色视频在线观看| 97人妻天天添夜夜摸| 国产av一区二区精品久久| 一个人观看的视频www高清免费观看 | 一二三四在线观看免费中文在| 国产成人系列免费观看| 欧美激情高清一区二区三区| 国产99久久九九免费精品| 国产主播在线观看一区二区| 国产亚洲欧美98| 少妇粗大呻吟视频| 久久精品成人免费网站| 欧美绝顶高潮抽搐喷水| 亚洲av日韩精品久久久久久密| 中文字幕最新亚洲高清| 美女免费视频网站| 亚洲人成伊人成综合网2020| 欧美日韩黄片免| 国产精品久久视频播放| 日韩精品中文字幕看吧| 制服丝袜大香蕉在线| 久久国产亚洲av麻豆专区| 深夜精品福利| √禁漫天堂资源中文www| 最好的美女福利视频网| 欧美人与性动交α欧美精品济南到| 欧美日韩精品网址| 国内毛片毛片毛片毛片毛片| 可以在线观看毛片的网站| 国产精品久久久久久精品电影 | 亚洲五月婷婷丁香| 成年版毛片免费区| 一个人免费在线观看的高清视频| 精品一区二区三区av网在线观看| 国产蜜桃级精品一区二区三区| 久久久久国产一级毛片高清牌| 好男人电影高清在线观看| 国产色视频综合| 国产精品亚洲一级av第二区| 久久久久久久久久久久大奶| 欧美绝顶高潮抽搐喷水| 精品久久久久久久毛片微露脸| 久久 成人 亚洲| 99热只有精品国产| 久久国产精品人妻蜜桃| 国产高清激情床上av| 大码成人一级视频| 国产野战对白在线观看| 热99re8久久精品国产| 精品电影一区二区在线| 成人av一区二区三区在线看| 最近最新中文字幕大全免费视频| 国产成+人综合+亚洲专区| 久久久国产精品麻豆| 丁香六月欧美| 天堂动漫精品| 黄色a级毛片大全视频| 欧美激情 高清一区二区三区| 午夜久久久在线观看| 亚洲国产看品久久| 欧美另类亚洲清纯唯美| 美女午夜性视频免费| 成人三级做爰电影| 亚洲欧美日韩无卡精品| 最新在线观看一区二区三区| 国产亚洲精品第一综合不卡| 十分钟在线观看高清视频www| 波多野结衣高清无吗| 黄色片一级片一级黄色片| 国产成人欧美| 午夜精品久久久久久毛片777| 动漫黄色视频在线观看| 欧美成人一区二区免费高清观看 | 女人爽到高潮嗷嗷叫在线视频| 日韩视频一区二区在线观看| 日日摸夜夜添夜夜添小说| 亚洲精品国产精品久久久不卡| 999久久久精品免费观看国产| 中亚洲国语对白在线视频| 少妇熟女aⅴ在线视频| 日本免费a在线| 国产亚洲精品综合一区在线观看 | www.熟女人妻精品国产| 免费av毛片视频| 日本在线视频免费播放| 一进一出好大好爽视频| 亚洲天堂国产精品一区在线| 久久中文看片网| 久久久久久久久免费视频了| АⅤ资源中文在线天堂| 久久精品国产综合久久久| 日本a在线网址| 亚洲人成网站在线播放欧美日韩| 91国产中文字幕| 老司机靠b影院| 国产麻豆成人av免费视频|