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

    基于力學(xué)平衡模型的生態(tài)系統(tǒng)服務(wù)協(xié)同與權(quán)衡關(guān)系研究——以長江中游城市群為例

    2023-11-28 12:56:36梁玉琦危小建謝名睿
    中國環(huán)境科學(xué) 2023年11期
    關(guān)鍵詞:生態(tài)服務(wù)研究

    梁玉琦,危小建,3*,江 平,蔡 進(jìn),謝名睿

    基于力學(xué)平衡模型的生態(tài)系統(tǒng)服務(wù)協(xié)同與權(quán)衡關(guān)系研究——以長江中游城市群為例

    梁玉琦1,2,危小建1,2,3*,江 平4,蔡 進(jìn)1,2,謝名睿1,2

    (1.東華理工大學(xué),自然資源部環(huán)鄱陽湖區(qū)域礦山環(huán)境監(jiān)測與治理重點(diǎn)實(shí)驗室,江西 南昌 330013;2.東華理工大學(xué)測繪與空間信息工程學(xué)院,江西 南昌 330013;3.江西省測繪地理信息工程技術(shù)研究中心,江西 南昌 330025;4.武漢大學(xué)資源與環(huán)境科學(xué)學(xué)院,湖北 武漢 430079)

    在運(yùn)用InVEST模型估算2000~2020年長江中游城市群6種生態(tài)系統(tǒng)服務(wù)的基礎(chǔ)上,通過主成分分析法選出3個作用力構(gòu)建力學(xué)平衡模型,根據(jù)合力的大小和偏離程度來評判研究區(qū)生態(tài)系統(tǒng)服務(wù)時空演變及協(xié)同與權(quán)衡關(guān)系特征.結(jié)果表明:研究期間6種生態(tài)系統(tǒng)服務(wù)時空分異特征明顯,呈現(xiàn)出產(chǎn)水量、土壤保持量和凈初級生產(chǎn)力增加,而糧食生產(chǎn)、碳儲量和生境質(zhì)量減少的趨勢;根據(jù)主成分分析的結(jié)果,將3個主成分分別重命名為調(diào)節(jié)支持服務(wù)、糧食供給服務(wù)、產(chǎn)水供給服務(wù).研究期間生態(tài)系統(tǒng)服務(wù)協(xié)調(diào)度合的均值為5.21,總體上處于基本協(xié)調(diào)狀態(tài),協(xié)同性有所提升且結(jié)構(gòu)趨于合理,同時協(xié)調(diào)度分布的區(qū)域差異明顯;協(xié)調(diào)度的偏離角分布在I、II、VI象限,產(chǎn)水供給服務(wù)主導(dǎo)區(qū)(I、VI)多于調(diào)節(jié)支持服務(wù)主導(dǎo)區(qū)(II).從協(xié)調(diào)性的流動性來看,雖存在跨越式的發(fā)展,但其水平仍有待提高.各縣區(qū)應(yīng)結(jié)合區(qū)域自然條件、資源稟賦和社會經(jīng)濟(jì)狀態(tài),優(yōu)化生態(tài)系統(tǒng)服務(wù)之間的關(guān)系.研究結(jié)果可為研究區(qū)生態(tài)建設(shè)和土地合理利用提供參考依據(jù)和決策輔助.

    生態(tài)系統(tǒng)服務(wù);協(xié)同與權(quán)衡;主成分分析法;力學(xué)平衡模型;長江中游城市群

    生態(tài)系統(tǒng)服務(wù)是指生態(tài)系統(tǒng)與生態(tài)過程所形成及所維持的人類賴以生存與發(fā)展的自然環(huán)境條件與效用[1],是人類從生態(tài)系統(tǒng)中所獲得的利益[2].具體而言,生態(tài)系統(tǒng)能夠為人類提供包括供給、調(diào)節(jié)、文化和支持在內(nèi)的生態(tài)服務(wù),從而滿足人類的各種需求[3].而如今,城鎮(zhèn)化和工業(yè)化的持續(xù)推進(jìn)加劇了生態(tài)環(huán)境的破壞,全球的生態(tài)系統(tǒng)服務(wù)迅速衰退[4].在此背景下,厘清生態(tài)系統(tǒng)服務(wù)之間的關(guān)系,探討生態(tài)系統(tǒng)的管理和保護(hù)顯得尤為緊迫和必要.

    由于人類活動和自然條件的相互作用,生態(tài)系統(tǒng)服務(wù)之間呈現(xiàn)出空間異質(zhì)性,難以達(dá)到平衡,表現(xiàn)為協(xié)同與權(quán)衡的關(guān)系[5].協(xié)同意味著生態(tài)系統(tǒng)服務(wù)之間相互協(xié)作共同增益,而權(quán)衡則是一種生態(tài)系統(tǒng)服務(wù)增加,其他生態(tài)系統(tǒng)服務(wù)隨之減弱的情況.目前,國內(nèi)外的許多學(xué)者對生態(tài)系統(tǒng)服務(wù)協(xié)同權(quán)衡關(guān)系進(jìn)行了大量研究.學(xué)者們通過運(yùn)用InVSET[6]、CASA[7]、RUSLE[8]、ESTD[9]等模型,使用空間疊置[10]、逐像元偏相關(guān)分析[11]、空間自相關(guān)分析[12]、情景模擬[13]、生態(tài)服務(wù)簇[14]等多種方法,從時間尺度[15]、空間尺度[16]、驅(qū)動因素[17]等方面展現(xiàn)生態(tài)系統(tǒng)服務(wù)協(xié)同權(quán)衡關(guān)系的特征,研究的對象及空間范圍囊括了耕地、草地、林地、流域、生態(tài)脆弱區(qū)及城市群.如劉海等[18]基于當(dāng)量因子法,以5年為時間間隔分析了丹江口水源區(qū)10種生態(tài)系統(tǒng)服務(wù)間的協(xié)同權(quán)衡關(guān)系;張靜靜等[19]從區(qū)域、南北坡、垂直帶等多個空間尺度探討了伏牛山地區(qū)森林生態(tài)系統(tǒng)服務(wù)的權(quán)衡協(xié)同效應(yīng);Jiang等[20]通過繪制英國鄉(xiāng)村70年來的生態(tài)系統(tǒng)服務(wù)和生物多樣性變化圖,從國家層面上來研究生態(tài)系統(tǒng)協(xié)同權(quán)衡關(guān)系.這些研究雖取得了豐碩的成果,但大多局限于兩兩生態(tài)系統(tǒng)服務(wù)間的關(guān)系探討,缺少從生態(tài)系統(tǒng)內(nèi)部子系統(tǒng)角度來分析生態(tài)服務(wù)間的動態(tài)協(xié)調(diào)關(guān)系,較少關(guān)注生態(tài)系統(tǒng)服務(wù)的發(fā)展短板.鑒于此,本文引入力學(xué)平衡模型,測度生態(tài)系統(tǒng)內(nèi)部子系統(tǒng)之間的協(xié)調(diào)關(guān)系,從力作用的角度為生態(tài)系統(tǒng)服務(wù)協(xié)同與權(quán)衡關(guān)系分析提供一種較為便捷的方法.

    長江中游城市群作為中國快速城鎮(zhèn)化的典型地區(qū),長期以來區(qū)域經(jīng)濟(jì)迅猛發(fā)展、城市空間格局和土地利用結(jié)構(gòu)劇變,致使自然景觀破壞、生態(tài)系統(tǒng)服務(wù)功能衰退,亟需提升社會經(jīng)濟(jì)與生態(tài)環(huán)境質(zhì)量的適配性[21-22].基于此,本文將以長江中游城市群為研究對象,通過InVEST模型測算2000~2020年6種典型的生態(tài)系統(tǒng)服務(wù)價值量,并借助主成分分析法選出3個作用力建立變化的力學(xué)平衡模型,根據(jù)合力的大小和偏離程度來識別各生態(tài)系統(tǒng)服務(wù)的相互匹配關(guān)系與發(fā)展短板[23-25],客觀地反映出生態(tài)服務(wù)在多個不同方向力量的作用下所產(chǎn)生的協(xié)同權(quán)衡關(guān)系特征,對研究期間長江中游城市群的生態(tài)系統(tǒng)服務(wù)進(jìn)行動態(tài)考量,以期為生態(tài)系統(tǒng)協(xié)調(diào)發(fā)展和科學(xué)制定生態(tài)保護(hù)策略提供參考和依據(jù).

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

    1.1 研究區(qū)概況

    圖1 研究區(qū)位置及2020年土地利用現(xiàn)狀

    長江中游城市群是由武漢城市圈、環(huán)長株潭城市群和環(huán)鄱陽湖城市群組成的以武漢為核心的特大城市群,其跨越了湖北、湖南和江西三省,囊括31個城市,國土面積約31.7萬km2.該區(qū)域?qū)賮啛釒Ъ撅L(fēng)氣候,雨熱同期,水資源豐富,并且擁有我國最大的淡水湖群.區(qū)域內(nèi)稀土、鈾礦等戰(zhàn)略性礦產(chǎn)豐富,頁巖氣資源儲量較大.土地類型多樣,以丘陵、平原和濕地為主.江漢平原、洞庭湖平原和鄱陽湖平原互聯(lián)互通,平原面積廣袤,是我國重要的糧食產(chǎn)區(qū).周邊有大別山、幕阜山和雪峰山,森林覆蓋率高,生態(tài)價值凸顯.長江中游城市群位于浙贛、京廣、武廣及京九鐵路的交匯處,工業(yè)門類較齊全,戰(zhàn)略性新興產(chǎn)業(yè)和服務(wù)業(yè)發(fā)展迅速.該區(qū)域是長江經(jīng)濟(jì)帶的重要組成部分,是中國經(jīng)濟(jì)的新增長極,在全國經(jīng)濟(jì)地理格局中處于承南啟北,連接?xùn)|西的重要地位.

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

    以長江中游城市群范圍內(nèi)的縣級行政單位作為研究單元,使用的數(shù)據(jù)包括:(1)2000年、2010年以及2020年三期土地利用數(shù)據(jù)來源于Zenodo數(shù)據(jù)庫(https://www.zenodo.org),空間分辨率為30m′30m; (2)DEM數(shù)據(jù)來源于地理空間數(shù)據(jù)云(https://www. gscloud.cn/),空間分辨率為30m; (3)2000年、2010年以及2020年的氣象數(shù)據(jù),包括氣溫、降水量數(shù)據(jù),來自國家地球系統(tǒng)科學(xué)數(shù)據(jù)中心(https://www. geodata.cn/),空間分辨率為1km;(4)土壤數(shù)據(jù),包括土壤類型和土壤質(zhì)地數(shù)據(jù),提取自北京大學(xué)城市與環(huán)境學(xué)院地理數(shù)據(jù)平臺(http://geodata.pku.edu.cn/)的世界土壤數(shù)據(jù)庫(HWSD)中國土壤特征數(shù)據(jù)集,空間分辨率為1km;(5)相關(guān)社會經(jīng)濟(jì)數(shù)據(jù)來源于各地的統(tǒng)計年鑒.

    2 研究方法

    本文的研究重點(diǎn)在于生態(tài)系統(tǒng)服務(wù)協(xié)同與權(quán)衡關(guān)系的探討,通過運(yùn)用InVEST模型、主成分分析法、力學(xué)平衡模型和馬爾可夫鏈方法[26],梳理研究區(qū)生態(tài)系統(tǒng)服務(wù)的現(xiàn)狀,探究生態(tài)系統(tǒng)服務(wù)協(xié)調(diào)性的時空分異特征和偏度特征,識別問題區(qū)域,提出改善生態(tài)系統(tǒng)服務(wù)的對策與建議,為長江中游城市群生態(tài)文明建設(shè)發(fā)展提供理論依據(jù).技術(shù)路線如圖2所示.

    圖2 技術(shù)路線

    2.1 生態(tài)系統(tǒng)服務(wù)測算

    根據(jù)研究區(qū)特征和生態(tài)系統(tǒng)服務(wù)研究的現(xiàn)狀與熱點(diǎn),選擇3個類別6種生態(tài)系統(tǒng)服務(wù)進(jìn)行研究,包括糧食生產(chǎn)(grain production,GP)和產(chǎn)水量(water yield,WP)2項供給服務(wù),碳儲量(carbon sequestration, CS)和土壤保持量(soil retention,SR)2項調(diào)節(jié)服務(wù),以及生境質(zhì)量(habitat quality,HQ)和凈初級生產(chǎn)力(net primary production, NPP)2項支持服務(wù).生態(tài)系統(tǒng)服務(wù)測算方法詳見表1.

    表1 生態(tài)系統(tǒng)服務(wù)測算方法

    2.2 主成分分析法

    主成分分析法是基于多指標(biāo)分析,通過合適的數(shù)學(xué)變換,使新變量——主成分成為原來多個變量的線性組合,并選取少數(shù)幾個在變差總信息量中比例較大的主成分來分析事物的一種方法.本文通過對研究區(qū)研究時段內(nèi)6種生態(tài)系統(tǒng)服務(wù)進(jìn)行主成分分析,根據(jù)主成分的因子荷載系數(shù)進(jìn)行總結(jié)與重命名,然后計算主成分得分,公式如下:

    式中:Score為主成分得分;a為各成分所包含的因子得分系數(shù),由SPSS軟件計算得出;Factor為因子值.

    2.3 生態(tài)系統(tǒng)服務(wù)協(xié)調(diào)性判別

    在城鎮(zhèn)化和社會經(jīng)濟(jì)活動活躍發(fā)展的同時,為了實(shí)現(xiàn)生態(tài)系統(tǒng)服務(wù)的協(xié)調(diào)發(fā)展,達(dá)到生態(tài)總體效益最大化,需要深入研究服務(wù)間各項評價指標(biāo)之間的均衡狀態(tài),通過技術(shù)和結(jié)構(gòu)調(diào)整緩和沖突,實(shí)現(xiàn)各類型服務(wù)間的協(xié)同與共生.因此,本文引入力作用的概念,構(gòu)建基于力學(xué)平衡模型的生態(tài)系統(tǒng)服務(wù)協(xié)調(diào)狀態(tài)判別方法,全面分析生態(tài)系統(tǒng)服務(wù)的協(xié)同與權(quán)衡關(guān)系.

    (1)協(xié)調(diào)性判別作用力

    現(xiàn)有的關(guān)于力學(xué)平衡模型的研究對于均衡狀態(tài)的定義,大致可分為兩種:一是在周華等[33]的研究中,將研究時段要素的起始狀態(tài)定義為均衡狀態(tài),同時用與起始狀態(tài)相比的變化量(正向或負(fù)向)表示作用力;二是在項曉敏等[34]的研究中,將要素的標(biāo)準(zhǔn)值定義為均衡狀態(tài),用與標(biāo)準(zhǔn)值相比的值來表示作用力.而本文基于前文的研究,采用主成分1、主成分2和主成分3的重命名結(jié)果作為生態(tài)系統(tǒng)服務(wù)協(xié)調(diào)性判別作用力的名稱.并結(jié)合了上述兩種方法,將變化量(三年主成分得分分別做差值)的標(biāo)準(zhǔn)值定義為均衡狀態(tài),然后與標(biāo)準(zhǔn)值進(jìn)行比較,將結(jié)果用來表示3個作用力的大小.

    (2)協(xié)調(diào)性判別模型

    假設(shè)其他條件均質(zhì),將主成分1、主成分2和主成分3抽象為笛卡爾坐標(biāo)系中三個不同方向的矢量力,以矢量力的合力及合力所在象限的結(jié)果來動態(tài)地檢測指標(biāo)間的協(xié)調(diào)程度及匹配關(guān)系[35-36].如果3個作用力均達(dá)到預(yù)期目標(biāo),則合力為0,表明3種生態(tài)系統(tǒng)服務(wù)均衡發(fā)展,此時合位于均衡點(diǎn),即圖3中點(diǎn).因此根據(jù)作用合力合的大小和偏離角度,可以研究生態(tài)系統(tǒng)服務(wù)的協(xié)調(diào)關(guān)系特征(概念模型見圖3).圖中:、、分別表示研究期間內(nèi)主成分1、主成分2和主成分3所代表的作用力的大小.計算中可用極坐標(biāo)(合,)來表示生態(tài)系統(tǒng)服務(wù)的協(xié)調(diào)性,合為極徑,其值越大說明生態(tài)系統(tǒng)服務(wù)協(xié)調(diào)性越差;極角表示合力與三個力之間的關(guān)系,反映服務(wù)協(xié)調(diào)性的具體特征.將、、的方向角度分別定義為π/2、7π/6、11π/6.

    根據(jù)向量運(yùn)算規(guī)則,極坐標(biāo)(合,)的計算公式如下:

    反向延長、、為、和,以此作為主成分1、主成分2和主成分3所代表力的逆向發(fā)展矢量方向,并將合力空間劃為6個象限(圖3),根據(jù)3種作用力間矢量正負(fù)關(guān)系相互組合情況,形成各象限矢量特征表(表2).

    圖3 基于力學(xué)平衡的生態(tài)系統(tǒng)服務(wù)協(xié)同與權(quán)衡模型

    表2 生態(tài)系統(tǒng)服務(wù)協(xié)調(diào)性偏離動力狀態(tài)及其特征

    3 結(jié)果與分析

    3.1 長江中游城市群生態(tài)系統(tǒng)服務(wù)空間分布格局

    長江中游城市群6種生態(tài)系統(tǒng)服務(wù)的空間差異性(以2020年為例)見圖4.糧食生產(chǎn)量主要成片集中在江漢平原、洞庭湖平原及鄱陽湖平原,這些地區(qū)地形平坦、土壤肥沃、耕地集中,是重要的農(nóng)作物產(chǎn)區(qū);產(chǎn)水量受降水、蒸發(fā)、土壤及植被等多種因素的影響,多集中在城市建設(shè)用地密度較高的地方,如贛江沿線及鄱陽湖周邊縣(市)區(qū)和以長沙為核心的周邊縣(市)區(qū);湘贛之間的羅霄山和贛東的武夷山森林覆蓋率高,碳密度高,碳儲量較高;土壤保持量主要分布在山區(qū)林草地密集區(qū)域,因植被保持水土的功能強(qiáng),故土壤不易被侵蝕.而城市周邊和平原地帶,植被覆蓋率較低,以建設(shè)用地和耕地為主,土壤易被侵蝕;生境質(zhì)量在遠(yuǎn)離城區(qū)的山地和丘陵區(qū)域較高,而中心城區(qū)及交通主干道附近社會經(jīng)濟(jì)發(fā)展迅速,受人類干擾較大,生物多樣性減少;凈初級生產(chǎn)力的高值主要分布在植被茂密的林地附近,而植被稀疏的城區(qū)凈初級生產(chǎn)力較低.

    圖4 2020年長江中游城市群生態(tài)系統(tǒng)服務(wù)空間分布格局特征

    從研究的整體時段來看,產(chǎn)水量、土壤保持量和凈初級生產(chǎn)力增加,而糧食生產(chǎn)、碳儲量和生境質(zhì)量減少(圖5).2000~2010年,江西和湖南交界的羅霄山脈和周邊山區(qū)糧食生產(chǎn)量下降.產(chǎn)水量呈增長趨勢,增長顯著區(qū)域集中在武漢的漢陽區(qū)、岳陽的岳陽縣及南昌、上饒和九江交接的部分縣區(qū),這主要是由于城市蔓延、不透水面面積增大所致.同時,大片的森林和農(nóng)田轉(zhuǎn)變?yōu)榻ㄔO(shè)用地,碳儲量和生境質(zhì)量整體呈下降趨勢,凈初級生產(chǎn)力在江西省內(nèi)下降嚴(yán)重. 2010~2020年,隨著經(jīng)濟(jì)發(fā)展,耕地面積減少,但資本和技術(shù)的投入提高了糧食生產(chǎn)率.產(chǎn)水量持續(xù)增長,但增長幅度有所放緩.碳儲量和土壤保持量增長區(qū)域主要集中在研究區(qū)西北部,凈初級生產(chǎn)力呈增長趨勢.由此說明人們實(shí)施的退耕還林還草、天然林保護(hù)工程、退田還湖等生態(tài)保護(hù)與恢復(fù)政策取得了一定成效,促進(jìn)了森林、草地和濕地的恢復(fù)與增加.

    3.2 長江中游城市群生態(tài)系統(tǒng)服務(wù)協(xié)同與權(quán)衡關(guān)系

    3.2.1 主成分因子分析 運(yùn)用SPSS對前文測算的6種生態(tài)系統(tǒng)服務(wù)值進(jìn)行主成分分析,根據(jù)檢驗結(jié)果,KMO的值為0.76,Barttett球形梯度檢驗對應(yīng)值為0.00,水平上呈現(xiàn)出顯著性,表明各生態(tài)系統(tǒng)服務(wù)間具有相關(guān)性,可以進(jìn)行主成分分析.根據(jù)因子荷載系數(shù)對主成分進(jìn)行重命名,以便提煉結(jié)果的實(shí)際意義.因子荷載系數(shù)結(jié)果見表3,前3個成分的累計方差貢獻(xiàn)率達(dá)91.6%,可解釋大部分變量信息,故確定出3個主成分.

    從因子荷載的絕對值大于0.4的結(jié)果來看,主成分1與碳儲量、土壤保持量、生境質(zhì)量和凈初級生產(chǎn)力的相關(guān)程度較大,故概括為“調(diào)節(jié)支持服務(wù)”;主成分2只與糧食生產(chǎn)這一個變量相關(guān)程度大,故概括為“糧食供給服務(wù)”;產(chǎn)水量在主成分3上荷載系數(shù)最高,故將主成分3概括為“產(chǎn)水供給服務(wù)”.

    表3 生態(tài)系統(tǒng)服務(wù)因子荷載系數(shù)表

    3.2.2 協(xié)調(diào)性時空分異特征 根據(jù)式(9)~(14)計算各縣區(qū)合力,并利用自然斷點(diǎn)法將生態(tài)系統(tǒng)服務(wù)協(xié)調(diào)性由小到大分為4類,分別表示高度協(xié)調(diào)(£3.36)、基本協(xié)調(diào)(3.36~6.45)、過渡類型(6.45~10.31)及失調(diào)類型(310.31).其中協(xié)同包括高度協(xié)調(diào)和基本協(xié)調(diào),權(quán)衡包括過渡類型和失調(diào)類型.據(jù)計算,研究期內(nèi)合的平均值為5.21,可見長江中游城市群的生態(tài)系統(tǒng)服務(wù)總體上處于基本協(xié)調(diào)狀態(tài),生態(tài)系統(tǒng)服務(wù)的改善有一定成效.

    在空間格局上(圖6),協(xié)同區(qū)域主要位于中部平原地帶,該區(qū)域地形平坦,水土肥沃,自然條件較好.高度協(xié)調(diào)的縣(市)區(qū)包括湖北的荊州區(qū)、孝南區(qū)和鄂城區(qū),江西的共青城市、永修縣和彭澤縣,湖南的云溪區(qū)、資陽區(qū)和津市市.位于該區(qū)的縣級單元隨著工業(yè)化和城市化的發(fā)展,經(jīng)濟(jì)發(fā)展較快,但其生態(tài)平衡并未破壞,調(diào)節(jié)支持服務(wù)、糧食供給服務(wù)和產(chǎn)水供給服務(wù)三者之間協(xié)同發(fā)展、彼此促進(jìn).同時,處于基本協(xié)調(diào)的縣(市)區(qū)在空間上多集中于高度協(xié)調(diào)區(qū)域的周圍,分布廣泛且聚集.該類地區(qū)的協(xié)調(diào)度雖稍低于高度協(xié)調(diào)型,但在發(fā)展中若采取積極的生態(tài)保護(hù)措施,便能收獲明顯的生態(tài)成效;權(quán)衡區(qū)域主要位于山脈綿延區(qū),包括研究區(qū)的西北角,江西懷玉山與武夷山附近的玉山縣、廣豐區(qū)和廣信區(qū).受制于自然資源條件和社會經(jīng)濟(jì)條件,該類區(qū)域生態(tài)系統(tǒng)服務(wù)協(xié)調(diào)性偏差,出現(xiàn)了調(diào)節(jié)支持、糧食供給和產(chǎn)水供給某一方發(fā)展過度或欠缺的權(quán)衡狀態(tài),各服務(wù)間相互阻擾,成效低下,是當(dāng)前長江中游城市群生態(tài)文明建設(shè)過程中的問題區(qū)域.總體上看,2000~2020年協(xié)同與權(quán)衡狀況有所改善,表現(xiàn)為協(xié)同區(qū)域的擴(kuò)張和權(quán)衡區(qū)域的縮小,協(xié)調(diào)性類型向高一級狀態(tài)轉(zhuǎn)變.

    從各類型協(xié)調(diào)度的數(shù)量變化來看(表4),總體上高度協(xié)調(diào)區(qū)和基本協(xié)調(diào)區(qū)占比大,協(xié)調(diào)性發(fā)展穩(wěn)中向好.在前兩個研究時段期間協(xié)同區(qū)域數(shù)量增加,權(quán)衡區(qū)域數(shù)量逐步減少.協(xié)調(diào)度均值和中位值下降,協(xié)調(diào)性提升.同時全距和標(biāo)準(zhǔn)差下降,表明各區(qū)域生態(tài)系統(tǒng)服務(wù)的協(xié)調(diào)度差異縮小,結(jié)構(gòu)趨于合理.從2000~2020年長時段來看,協(xié)調(diào)度均值為6.64,總體上處于過渡狀態(tài),協(xié)同區(qū)域多于權(quán)衡區(qū)域.綜上,長江中游城市群生態(tài)系統(tǒng)服務(wù)的協(xié)調(diào)狀況有所改善,生態(tài)狀況良好.

    3.2.3 協(xié)調(diào)性偏度特征 通過力學(xué)平衡模型計算得到偏離角,以偏離角所在象限來分析偏離特征,如圖7.結(jié)果顯示偏離角分布在I、II、VI象限,總體上產(chǎn)水供給服務(wù)主導(dǎo)區(qū)(I、VI)多于調(diào)節(jié)支持服務(wù)主導(dǎo)區(qū)(II).

    產(chǎn)水供給服務(wù)主導(dǎo)區(qū)集中反映為產(chǎn)水供給服務(wù)偏高的同時,相應(yīng)的調(diào)節(jié)支持或糧食供給服務(wù)卻低于既定目標(biāo),制約著生態(tài)系統(tǒng)服務(wù)的協(xié)調(diào)發(fā)展.2000~2010年長江中游城市群21.3%的縣區(qū)處于第Ⅰ象限,61.6%的縣區(qū)處于第VI象限.2010~2020年,第Ⅰ象限占比由21.3%下降到15.6%,第VI象限的縣區(qū)數(shù)量保持不變.隨著生態(tài)文明理念的傳播與貫徹實(shí)施,使得部分縣區(qū)從第Ⅰ象限轉(zhuǎn)向第Ⅱ象限.Ⅰ象限主要分布在傳統(tǒng)農(nóng)區(qū),如江漢平原的武漢、仙桃和孝感,鄱陽湖平原和洞庭湖平原周圍的縣域.這些區(qū)域農(nóng)耕歷史悠久,但隨著人為活動的增多,耕地的生態(tài)功能遭到破壞,糧食供給服務(wù)的增長受到限制.該地區(qū)今后應(yīng)以現(xiàn)代農(nóng)業(yè)為依托,充分發(fā)揮其旅游休閑功能與農(nóng)耕文化傳播功能,從而促進(jìn)經(jīng)濟(jì)與生態(tài)的協(xié)同發(fā)展;第VI象限的縣區(qū)占比最大,這些區(qū)域由于建設(shè)用地的擴(kuò)張,產(chǎn)水供給服務(wù)逐步增長.但其對于生態(tài)保護(hù)政策的要求未落實(shí)到位,相應(yīng)的保護(hù)規(guī)模較小且成效不甚顯著,未能充分挖掘生態(tài)系統(tǒng)的調(diào)節(jié)支持服務(wù)價值.因此這些縣區(qū)后續(xù)應(yīng)注重經(jīng)濟(jì)規(guī)劃和生態(tài)規(guī)劃的一致性,優(yōu)化經(jīng)濟(jì)發(fā)展格局,提升生態(tài)保護(hù)的效果.

    調(diào)節(jié)支持服務(wù)主導(dǎo)區(qū)偏向于區(qū)域內(nèi)調(diào)節(jié)支持服務(wù)和糧食供給服務(wù)的不協(xié)調(diào),即第II象限的縣區(qū)生態(tài)保護(hù)措施實(shí)施良好,而相匹配的糧食供給偏低.2000~2010年第II象限的縣區(qū)數(shù)量占比為17.1%, 2010~2020年占比為22.8%,其中包括湖北的荊州和宜昌、江西的九江和撫州,及湖南的長沙、湘潭和常德等地方.2000~2020年,受經(jīng)濟(jì)轉(zhuǎn)型和耕地面積減少的影響,該區(qū)域糧食生產(chǎn)量呈下降趨勢.今后這些地區(qū)應(yīng)結(jié)合國家和區(qū)域規(guī)劃的要求,明確生態(tài)保護(hù)的目標(biāo)和重點(diǎn),同時加強(qiáng)用途管制,加大保護(hù)耕地.

    圖7 2000~2020年各縣區(qū)生態(tài)系統(tǒng)服務(wù)協(xié)調(diào)度象限分布

    從長江中游城市群生態(tài)系統(tǒng)服務(wù)協(xié)調(diào)性的馬爾可夫鏈結(jié)果來看(表5),協(xié)調(diào)性的流動性存在差異.從前兩個研究時段的對比情況來看,象限流動性排序從高到低依次為第Ⅵ象限、第Ⅱ象限、第Ⅰ象限,其他象限向第Ⅵ象限轉(zhuǎn)移的可能性更高.由此可知協(xié)調(diào)性的轉(zhuǎn)移不僅發(fā)生在相鄰狀態(tài)中(第Ⅰ、Ⅵ象限相互轉(zhuǎn)移),同時也存在跨越式發(fā)展的現(xiàn)象(第Ⅰ、Ⅵ象限轉(zhuǎn)為第Ⅱ象限).相鄰象限轉(zhuǎn)移的最高概率為55%,而跨越象限轉(zhuǎn)移的最高概率為29%,可見產(chǎn)水供給服務(wù)主導(dǎo)區(qū)和調(diào)節(jié)支持服務(wù)主導(dǎo)區(qū)存在著時間上的躍遷,多數(shù)地方的生態(tài)系統(tǒng)服務(wù)協(xié)調(diào)性有所提升.同時,將10年間隔的結(jié)果與20年間隔的結(jié)果對比來看,保持在原有象限不變的概率遠(yuǎn)高于象限流動的概率.因此,從長期來看,研究區(qū)的生態(tài)系統(tǒng)服務(wù)協(xié)調(diào)性較為穩(wěn)定.綜上,長江中游城市群生態(tài)系統(tǒng)服務(wù)的協(xié)調(diào)性雖存在跨越式的發(fā)展,但其水平仍有待提高,同時需持續(xù)采取科學(xué)的措施,突破現(xiàn)有狀態(tài),促使各地的生態(tài)系統(tǒng)服務(wù)協(xié)調(diào)性向高水平類型轉(zhuǎn)移.

    表5 2000~2020年生態(tài)系統(tǒng)服務(wù)協(xié)調(diào)性的馬爾可夫轉(zhuǎn)移矩陣表

    4 討論

    長江中游城市群生態(tài)系統(tǒng)服務(wù)整體上協(xié)同大于權(quán)衡,并且在2000~2010年、2010~2020年和2000~ 2020年三個研究時段下協(xié)調(diào)性(合,)在協(xié)調(diào)度、協(xié)同特征和空間分布上發(fā)生明顯變化,這也在一定程度上說明研究區(qū)生態(tài)系統(tǒng)服務(wù)的協(xié)調(diào)性改善,并持續(xù)向好發(fā)展.生態(tài)系統(tǒng)服務(wù)不僅取決于所處區(qū)域自身的資源環(huán)境狀態(tài),而且會受到鄰近地域的影響,研究區(qū)協(xié)調(diào)性的流動水平仍需繼續(xù)提升,向高水平發(fā)展.因此,建議未來生態(tài)系統(tǒng)服務(wù)管理突破行政區(qū)劃限制,實(shí)施長江中游城市群地區(qū)生態(tài)系統(tǒng)服務(wù)聯(lián)合管理與保護(hù)[37].限制城鎮(zhèn)建設(shè)用地對耕地的占用,推進(jìn)土地整治工程,發(fā)展綠色農(nóng)業(yè),強(qiáng)化區(qū)域糧食安全與生態(tài)安全責(zé)任,以驅(qū)動城鎮(zhèn)化與生態(tài)系統(tǒng)服務(wù)耦合協(xié)調(diào).

    根據(jù)本文的表3可知,前三個主成分的累計方差貢獻(xiàn)率達(dá)到了91.6%,具有顯著的代表性,可以滿足分析問題的需要[38].因此本文運(yùn)用主成分分析法降到三維,使用三個主成分來構(gòu)建力學(xué)平衡模型,這針對本研究區(qū)長江中游城市群是適用的.但針對于其他面積更大的研究區(qū)來說,降到三維的處理方法適不適用,這個問題有待今后做進(jìn)一步的研究.同時,需要指出的是,當(dāng)涉及到更多要素時,將諸多要素降到二維或-1維也可以在模型中實(shí)現(xiàn)平衡狀態(tài).力學(xué)平衡模型可以發(fā)展成四個(或更多)的作用力,但同時各個作用力矢量正負(fù)關(guān)系相互組合成的象限數(shù)量會越來越復(fù)雜,不利于分析,而三個力,數(shù)量適中,也是文章中常用的形式.本文在前人研究的基礎(chǔ)上,通過變化量的標(biāo)準(zhǔn)值來構(gòu)建變化的力學(xué)平衡模型,很好地揭示了生態(tài)系統(tǒng)服務(wù)動態(tài)的演變規(guī)律,清晰地表達(dá)出多種類型服務(wù)之間的空間關(guān)聯(lián)及發(fā)展短板.但該方法難以全面認(rèn)識不同維度(要素)之間的相互作用機(jī)制,并要求參與評價的要素具有較強(qiáng)的獨(dú)立性,這是今后應(yīng)用力學(xué)平衡模型時需要注意的問題.

    研究生態(tài)系統(tǒng)服務(wù)及其相互關(guān)系的動態(tài)變化,度量服務(wù)間協(xié)調(diào)狀況的好壞,可以加強(qiáng)對于生態(tài)系統(tǒng)服務(wù)的理解,實(shí)現(xiàn)從生態(tài)科學(xué)到生態(tài)決策的轉(zhuǎn)型,以此來追求生態(tài)系統(tǒng)總體最優(yōu).在生態(tài)系統(tǒng)服務(wù)評估研究中,InVEST模型以簡單、快速等優(yōu)點(diǎn)被廣泛地應(yīng)用,然而其也存在一定的局限性.由于模型結(jié)果與方法的簡化,加之參數(shù)不統(tǒng)一、缺乏實(shí)測數(shù)據(jù)等問題,因此導(dǎo)致模型模擬結(jié)果具有不確定性[39].今后的研究應(yīng)更關(guān)注參數(shù)本地化,并通過實(shí)地調(diào)研對結(jié)果的合理性進(jìn)行驗證,進(jìn)而提高InVEST模型評估結(jié)果的準(zhǔn)確性.同時,城市群發(fā)展具有復(fù)雜性,梳理生態(tài)系統(tǒng)服務(wù)的權(quán)衡與協(xié)同關(guān)系還應(yīng)考慮響應(yīng)時間和空間尺度,而本文僅對空間權(quán)衡的時空特點(diǎn)做了初步探討.因此,深入量化生態(tài)系統(tǒng)服務(wù)產(chǎn)生協(xié)同與權(quán)衡關(guān)系的原因,并通過設(shè)置不同發(fā)展模式和政策情景,預(yù)測生態(tài)保護(hù)措施的效果[40]也是后續(xù)研究的重要方向.

    5 結(jié)論

    5.1 受不同環(huán)境因素影響,長江中游城市群6種典型的生態(tài)系統(tǒng)服務(wù)時空分異特征明顯,其中糧食生產(chǎn)、碳儲量、土壤保持量、生境質(zhì)量和凈初級生產(chǎn)力均呈現(xiàn)出郊區(qū)高、城區(qū)低的分布格局,而產(chǎn)水量與之相反.生態(tài)系統(tǒng)服務(wù)格局變化主要體現(xiàn)在人類活動干擾強(qiáng)烈的地區(qū),在2000~2020年間呈現(xiàn)出產(chǎn)水量、土壤保持量和凈初級生產(chǎn)力增加,而糧食生產(chǎn)、碳儲量和生境質(zhì)量減少的趨勢.

    5.2 將6種生態(tài)系統(tǒng)服務(wù)的估算值進(jìn)行主成分分析的結(jié)果符合相關(guān)條件,表明各生態(tài)系統(tǒng)服務(wù)間具有相關(guān)性.并且根據(jù)主成分因子荷載系數(shù)的大小,分別將主成分1、2、3重命名為調(diào)節(jié)支持服務(wù)、糧食供給服務(wù)、產(chǎn)水供給服務(wù).

    5.3 研究期間長江中游城市群生態(tài)系統(tǒng)服務(wù)協(xié)調(diào)度合的均值為5.21,總體上處于基本協(xié)調(diào)狀態(tài),協(xié)同性有所提升,結(jié)構(gòu)趨于合理.同時協(xié)調(diào)度分布的區(qū)域差異明顯,協(xié)同區(qū)域主要位于中部平原地帶,權(quán)衡區(qū)域主要位于山脈綿延區(qū).

    5.4 研究期間長江中游城市群生態(tài)系統(tǒng)服務(wù)協(xié)調(diào)度的偏離角分布在I、II、VI象限,總體上產(chǎn)水供給服務(wù)主導(dǎo)區(qū)(I、VI)多于調(diào)節(jié)支持服務(wù)主導(dǎo)區(qū)(II).從協(xié)調(diào)性的流動性來看,雖存在跨越式的發(fā)展,但其水平仍有待提高,同時需持續(xù)采取科學(xué)的措施,突破現(xiàn)有狀態(tài),促使各地的生態(tài)系統(tǒng)服務(wù)協(xié)調(diào)性向高水平類型轉(zhuǎn)移.大多數(shù)縣區(qū)屬于產(chǎn)水供給服務(wù)主導(dǎo)縣區(qū),但各類型區(qū)域生態(tài)系統(tǒng)服務(wù)協(xié)調(diào)發(fā)展所面臨的制約因素有所差異,因此應(yīng)結(jié)合區(qū)域自然條件、資源稟賦和社會經(jīng)濟(jì)狀態(tài),制定差異化的提升方案.

    [1] Daily G C. Nature's services: Societal dependence on natural ecosystems [M]. Washington DC: Island Press, 1997.

    [2] Costanza R, D Arge R, De Groot R, et al. The value of the world's ecosystem services and natural capital [J]. Nature, 1997,387(6630): 253-260.

    [3] 李文華,張 彪,謝高地.中國生態(tài)系統(tǒng)服務(wù)研究的回顧與展望 [J]. 自然資源學(xué)報, 2009,24(1):1-10. Li W H, Zhang B, Xie G D. Research on ecosystem services in China:Progress and perspectives [J]. Journal of Natural Resources, 2009,24(1):1-10.

    [4] Tolessa T, Senbeta F, Kidane M. The impact of land use/land cover change on ecosystem services in the central highlands of Ethiopia [J]. Ecosystem Services, 2017,23:47-54.

    [5] Liu C, Xu Y, Huang A, et al. Spatial identification of land use multifunctionality at grid scale in farming-pastoral area: A case study of Zhangjiakou City, China [J]. Habitat International, 2018,76:48-61.

    [6] 王啟名,楊 昆,李立曉,等.滇池流域水文生態(tài)系統(tǒng)服務(wù)權(quán)衡與協(xié)同時空異質(zhì)性及其歸因分析 [J]. 生態(tài)學(xué)報, 2023,43(12):4876-4891. Wang Q M, Yang K, Li L X, et al. Spatio-temporal heterogeneity and attribution analysis of hydrological ecosystem services tradeoffs and synergies in Dianchi Lake Basin [J]. Acta Ecologica Sinica, 2023, 43(12):4876-4891.

    [7] 孫藝杰,任志遠(yuǎn),趙勝男,等.陜西河谷盆地生態(tài)系統(tǒng)服務(wù)協(xié)同與權(quán)衡時空差異分析 [J]. 地理學(xué)報, 2017,72(3):521-532. Sun Y J, Ren Z Y, Zhao S N, et al. Spatial and temporal changing analysis of synergy and trade-off between ecosystem services in valley basins of Shaanxi Province [J]. Acta Geographica Sinica, 2017, 72(3):521-532.

    [8] 仲俊濤,王 蓓,米文寶,等.農(nóng)牧交錯帶禁牧草地生態(tài)系統(tǒng)服務(wù)空間權(quán)衡與協(xié)同關(guān)系 [J]. 農(nóng)業(yè)工程學(xué)報, 2020,36(12):268-275. Zhong J T, Wang B, Mi W B, et al. Spatial trade-offs and synergies of ecosystem services for grazing-prohibited grassland in agro-pasture ecotone [J]. Transactions of the Chinese Society of Agricultural Engineering, 2020,36(12):268-275.

    [9] 李 理,朱連奇,朱文博,等.生態(tài)系統(tǒng)服務(wù)價值與人類活動強(qiáng)度關(guān)聯(lián)性分析及權(quán)衡——以淇河流域為例 [J]. 中國環(huán)境科學(xué), 2020,40(1): 365-374. Li L, Zhu L Q, Zhu W B, et al. The correlation between ecosystem service value and human activity intensity and its trade-offs—take Qihe River basin for example [J]. China Environmental Science, 2020, 40(1):365-374.

    [10] 劉 洋,畢 軍,呂建樹.生態(tài)系統(tǒng)服務(wù)權(quán)衡與協(xié)同關(guān)系及驅(qū)動力——以江蘇省太湖流域為例 [J]. 生態(tài)學(xué)報, 2019,39(19):7067-7078. Liu Y, Bi J, Lü J S. Trade-off and synergy relationships of ecosystem services and the driving forces: A case study of the Taihu Basin, Jiangsu Province [J]. Acta Ecologica Sinica, 2019,39(19):7067-7078.

    [11] 陳心盟,王曉峰,馮曉明,等.青藏高原生態(tài)系統(tǒng)服務(wù)權(quán)衡與協(xié)同關(guān)系 [J]. 地理研究, 2021,40(1):18-34. Chen X M, Wang X F, Feng X M, et al. Ecosystem service trade-off and synergy on Qinghai-Tibet Plateau [J]. Geographical Research, 2021,40(1):18-34.

    [12] 齊 麟,張 月,許 東,等.東北森林屏障帶生態(tài)系統(tǒng)服務(wù)權(quán)衡與協(xié)同關(guān)系 [J]. 生態(tài)學(xué)雜志, 2021,40(11):3401-3411. Qi L, Zhang Y, Xu D, et al. Trade-offs and synergies of ecosystem services in forest barrier belt of Northeast China [J]. Chinese Journal of Ecology, 2021,40(11):3401-3411.

    [13] Bohensky E L, Reyers B, Van Jaarsveld A S. Conservation in practice: Future ecosystem services in a Southern African River Basin: A scenario planning approach to uncertainty [J]. Conservation Biology, 2006,20(4):1051-1061.

    [14] 祁 寧,趙 君,楊延征,等.基于服務(wù)簇的東北地區(qū)生態(tài)系統(tǒng)服務(wù)權(quán)衡與協(xié)同 [J]. 生態(tài)學(xué)報, 2020,40(9):2827-2837. Qi N, Zhao J, Yang Y Z, et al. Quantifying ecosystem service trade- offs and synergies in Northeast China based on ecosystem service bundles [J]. Acta Ecologica Sinica, 2020,40(9):2827-2837.

    [15] Schroder S A, Tóth S F, Deal R L, et al. Multi-objective optimization to evaluate tradeoffs among forest ecosystem services following fire hazard reduction in the Deschutes National Forest, USA [J]. Ecosystem Services, 2016,22:328-347.

    [16] Vallet A, Locatelli B, Levrel H, et al. Relationships between ecosystem services: Comparing methods for assessing tradeoffs and synergies [J]. Ecological Economics, 2018,150:96-106.

    [17] Palacios-Agundez I, Onaindia M, Barraqueta P, et al. Provisioning ecosystem services supply and demand: The role of landscape management to reinforce supply and promote synergies with other ecosystem services [J]. Land Use Policy, 2015,47:145-155.

    [18] 劉 海,武 靖,陳曉玲.丹江口水源區(qū)生態(tài)系統(tǒng)服務(wù)時空變化及權(quán)衡協(xié)同關(guān)系 [J]. 生態(tài)學(xué)報, 2018,38(13):4609-4624. Liu H, Wu J, Chen X L. Study on spatial-temporal change and trade- off/synergy relationships of ecosystem services in the Danjiangkou water source area [J]. Acta Ecologica Sinica, 2018,38(13):4609-4624.

    [19] 張靜靜,朱文博,朱連奇,等.伏牛山地區(qū)森林生態(tài)系統(tǒng)服務(wù)權(quán)衡/協(xié)同效應(yīng)多尺度分析 [J]. 地理學(xué)報, 2020,75(5):975-988. Zhang J J, Zhu W B, Zhu L Q, et al. Multi-scale analysis of trade- off/synergy effects of forest ecosystem services in the Funiu Mountain Region [J]. Acta Geographica Sinica, 2020,75(5):975-988.

    [20] Jiang M, Bullock J M, Hooftman D A P. Mapping ecosystem service and biodiversity changes over 70years in a rural English county [J]. Journal of Applied Ecology, 2013,50(4):841-850.

    [21] Jiang W, Chen Z, Lei X, et al. Simulation of urban agglomeration ecosystem spatial distributions under different scenarios: A case study of the Changsha-Zhuzhou–Xiangtan urban agglomeration [J]. Ecological Engineering, 2016,88:112-121.

    [22] Gu Q, Wang H, Zheng Y, et al. Ecological footprint analysis for urban agglomeration sustainability in the middle stream of the Yangtze River [J]. Ecological Modelling, 2015,318:86-99.

    [23] 班茂盛,方創(chuàng)琳,劉曉麗,等.北京高新技術(shù)產(chǎn)業(yè)區(qū)土地利用績效綜合評價 [J]. 地理學(xué)報, 2008,63(2):175-184. Ban M S, Fang C L, Liu X L, et al. Comprehensive evaluation of land use performance of Beijing High-tech Industrial Zone [J]. Acta Geographica Sinica, 2008,63(2):175-184.

    [24] 余富祥,胡月明,劉振杰,等.基于力學(xué)平衡模型的耕地多功能演變及其協(xié)調(diào)性分析——以珠江三角洲為例 [J]. 農(nóng)業(yè)資源與環(huán)境學(xué)報, 2019,36(6):728-737. Yu F X, Hu Y M, Liu Z J, et al. The coordination and evolution of farmland multifunctionality based on a mechanical model:A case study of the Pearl River Delta [J]. Journal of Agricultural Resources and Environment, 2019,36(6):728-737.

    [25] 葉士琳,曹有揮,蔣自然,等.基于力學(xué)平衡模型的長三角港口物流發(fā)展協(xié)調(diào)性研究 [J]. 地理科學(xué), 2017,37(11):1624-1631. Ye S L, Cao Y H, Jiang Z R, et al. Coordination of port logistics development based on mechanical model:A case of the Yangtze River Delta port system [J]. Scientia Geographica Sinica, 2017,37(11): 1624-1631.

    [26] 蒲英霞,馬榮華,葛 瑩,等.基于空間馬爾可夫鏈的江蘇區(qū)域趨同時空演變 [J]. 地理學(xué)報, 2005,60(5):817-826. Pu Y X, Ma R H, Ge Y, et al. Spatial-temporal dynamics of Jiangsu regional convergence with spatial markov chains approach [J]. Acta Geographica Sinica, 2005,60(5):817-826.

    [27] Cui F, Tang H, Zhang Q, et al. Integrating ecosystem services supply and demand into optimized management at different scales: A case study in Hulunbuir, China [J]. Ecosystem Services, 2019,39:100984.

    [28] 潘 韜,吳紹洪,戴爾阜,等.基于InVEST模型的三江源區(qū)生態(tài)系統(tǒng)水源供給服務(wù)時空變化 [J]. 應(yīng)用生態(tài)學(xué)報, 2013,24(1):183-189. Pan T, Wu S H, Dai E F, et al. Spatiotemporal variation of water source supply service in Three Rivers Source Area of China based on InVEST model [J]. Chinese Journal of Applied Ecology, 2013,24(1): 183-189.

    [29] 趙武生,石培基.基于InVEST模型的復(fù)合生態(tài)系統(tǒng)耦合協(xié)調(diào)關(guān)系研究——以蘭西城市群為例 [J]. 中國環(huán)境科學(xué), 2023,43(4):1883- 1894. Zhao W S, Shi P J. The coupling and coordination of complex ecosystem based on the InVEST model- A case study in the Lanzhou- Xining urban agglomeration [J]. China Environmental Science, 2023, 43(4):1883-1894.

    [30] Zhang L, Fu B, Lü Y, et al. Balancing multiple ecosystem services in conservation priority setting [J]. Landscape Ecology, 2015,30(3):535- 546.

    [31] 鞏 杰,馬學(xué)成,張玲玲,等.基于InVEST模型的甘肅白龍江流域生境質(zhì)量時空分異 [J]. 水土保持研究, 2018,25(3):191-196. Gong J, Ma X C, Zhang L L, et al. Spatiotemporal variation of habitat quality in Bailongjiang Watershed in Gansu based on InVEST model [J]. Research of Soil and Water Conservation, 2018,25(3):191-196.

    [32] 朱文泉,潘耀忠,張錦水.中國陸地植被凈初級生產(chǎn)力遙感估算 [J]. 植物生態(tài)學(xué)報, 2007,31(3):413-424. Zhu W Q, Pan Y Z, Zhang J S. Estimation of net primary productivity of Chinese terrestrial vegetation based on remote sensing [J]. Chinese Journal of Plant Ecology, 2007,31(3):413-424.

    [33] 周 華,馬春寧,周生路,等.基于力學(xué)平衡模型的鄉(xiāng)村轉(zhuǎn)型均衡發(fā)展判別方法研究 [J]. 長江流域資源與環(huán)境, 2014,23(3):303-310. Zhou H, Ma C N, Zhou S L, et al. The study of discrimination method of rural coordinated transformation development based on mechanical model [J]. Resources and Environment in the Yangtze Basin, 2014, 23(3):303-310.

    [34] 項曉敏,金曉斌,杜心棟,等.基于“強(qiáng)度—潛力—難度”綜合測度的中國農(nóng)用地整治實(shí)施協(xié)調(diào)性分析 [J]. 地理研究, 2016,35(2):285- 298. Xiang X M, Jin X B, Du X D, et al. The coordination of farmland consolidation implementation in China by comprehensive measure of "intensity, potential and difficulty" [J]. Geographical Research, 2016, 35(2):285-298.

    [35] 康 慶,郭青霞,丁 一,等.山西省“三生”功能協(xié)同/權(quán)衡關(guān)系分析 [J]. 自然資源學(xué)報, 2021,36(5):1195-1207. Kang Q, Guo Q X, Ding Y, et al. Tradeoffs/synergies analysis of "Production-Living-Ecological" functions in Shanxi province [J]. Journal of Natural Resources, 2021,36(5):1195-1207.

    [36] 張富剛,劉彥隨.中國區(qū)域農(nóng)村發(fā)展動力機(jī)制及其發(fā)展模式 [J]. 地理學(xué)報, 2008,63(2):115-122. Zhang F G, Liu Y S. Dynamic mechanism and models of regional rural development in China [J]. Acta Geographica Sinica, 2008,63(2):115- 122.

    [37] 賀清云,李慧平,歐陽曉.長江中游城市群城市化與生態(tài)環(huán)境耦合協(xié)調(diào)分析及模擬預(yù)測 [J]. 生態(tài)科學(xué), 2020,39(2):182-190. He Q Y, Li H P, Ouyang X. Analyzing and predicting the urbanization and eco-environment coupling and coordination of the middle reaches of Yangtze River urban agglomerations [J]. Ecological Science, 2020, 39(2):182-190.

    [38] 任俊霖,李 浩,伍新木,等.基于主成分分析法的長江經(jīng)濟(jì)帶省會城市水生態(tài)文明評價 [J]. 長江流域資源與環(huán)境, 2016,25(10):1537- 1544. Ren J L, Li H, Wu X M, et al. Assessment of 11provincial capitals' water ecological civilization of the Yangtze River Economic Belt on the principal component [J]. Resources and Environment in the Yangtze Basin, 2016,25(10):1537-1544.

    [39] 牛麗楠,邵全琴,寧 佳,等.西部地區(qū)生態(tài)狀況變化及生態(tài)系統(tǒng)服務(wù)權(quán)衡與協(xié)同 [J]. 地理學(xué)報, 2022,77(1):182-195. Niu L N, Shao Q Q, Ning J, et al. Ecological changes and the tradeoff and synergy of ecosystem services in western China [J]. Acta Geographica Sinica, 2022,77(1):182-195.

    [40] 祝 萍,劉 鑫,鄭瑜晗,等.北方重點(diǎn)生態(tài)功能區(qū)生態(tài)系統(tǒng)服務(wù)權(quán)衡與協(xié)同 [J]. 生態(tài)學(xué)報, 2020,40(23):8694-8706. Zhu P, Liu X, Zheng Y H, et al. Tradeoffs and synergies of ecosystem services in key ecological function zones in northern China [J]. Acta Ecologica Sinica, 2020,40(23):8694-8706.

    Study on trade-offs and synergies of ecosystem services based on mechanical equilibrium model:A case of the Middle Reaches of the Yangtze River Urban Agglomerations.

    LIANG Yu-qi1,2, WEI Xiao-jian1,2,3*,JIANG Ping4,CAI Jin1,2,XIE Ming-rui1,2

    (1.Key Laboratory of Mine Environmental Monitoring and Improving around Poyang Lake of Ministry of Natural Resources, East China University of Technology, Nanchang 330013, China;2.School of Surveying and Geoinformation Engineering, East China University of Technology, Nanchang 330013, China;3.Jiangxi Province Engineering Research Center of Surveying, Mapping and Geographic Information, Nanchang 330025, China;4.School of Resource and Environmental Sciences ,Wuhan University, Wuhan 430079, China)., 2023,43(11):5974~5986

    Based on the use of the InVEST model to estimate six ecosystem services in the middle reaches of the Yangtze River urban agglomeration from 2000 to 2020, a mechanical equilibrium model was constructed by selecting three forces through Principal Component Analysis (PCA). The spatiotemporal evolution of ecosystem services and the characteristics of trade-offs and synergies within the study area were assessed based on the magnitude and deviation of the resultant forces. The conclusions can be summarized as follows: the spatial and temporal characteristics of the six ecosystem services during the study period were obvious, demonstrating the trend of increasing water yield, soil retention and net primary production, while decreasing grain production, carbon sequestration and habitat quality. On the basis of the findings of PCA, the three principal components were renamed as regulation support service, food supply service and water production supply service. The average of coordination variabletotalwas 5.21, indicating a generally basic coordination state. The synergies among ecosystem services improved and the overall structure tended to be more reasonable. Meanwhile, significant regional differences in coordination degree distribution were observed. Coordination skewnesswas mainly distributed in quadrants I, II and VI, with a greater number of dominant areas for water production supply service (I and VI) compared to regulation support service (II). Regarding coordinated mobility, although there was some leapfrog development, its level still needs to be enhanced. All counties and districts should optimize the relationship between ecosystem services depending on regional natural conditions, resource endowments and socio-economic status. The research results can provide references and decision-making assistance for ecological construction and rational land utilization in the research area.

    ecosystem services;trade-offs and synergies;principal component analysis;mechanical equilibrium model;the middle reaches of the Yangtze River urban agglomerations

    X171

    A

    1000-6923(2023)11-5974-13

    梁玉琦(2000-),女,湖南長沙人,東華理工大學(xué)碩士研究生,主要研究方向為土地利用規(guī)劃與評價.liang_yuqi22@163.com.

    梁玉琦,危小建,江 平,等.基于力學(xué)平衡模型的生態(tài)系統(tǒng)服務(wù)協(xié)同與權(quán)衡關(guān)系研究——以長江中游城市群為例 [J]. 中國環(huán)境科學(xué), 2023,43(11):5974-5986.

    Liang Y Q, Wei X J,Jiang P, et al. Study on trade-offs and synergies of ecosystem services based on mechanical equilibrium model:A case of the Middle Reaches of the Yangtze River Urban Agglomerations [J]. China Environmental Science, 2023,43(11):5974-5986.

    2023-04-06

    國家自然科學(xué)基金資助項目(52168010)

    * 責(zé)任作者, 副教授, 631000872@qq.com

    猜你喜歡
    生態(tài)服務(wù)研究
    FMS與YBT相關(guān)性的實(shí)證研究
    “生態(tài)養(yǎng)生”娛晚年
    遼代千人邑研究述論
    住進(jìn)呆萌生態(tài)房
    生態(tài)之旅
    視錯覺在平面設(shè)計中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    服務(wù)在身邊 健康每一天
    EMA伺服控制系統(tǒng)研究
    服務(wù)在身邊 健康每一天
    服務(wù)在身邊 健康每一天
    久久人人爽人人爽人人片va| 日本黄色日本黄色录像| 人妻制服诱惑在线中文字幕| 欧美激情极品国产一区二区三区 | 人妻系列 视频| 精品久久久精品久久久| 大片免费播放器 马上看| 国模一区二区三区四区视频| 日日摸夜夜添夜夜爱| 亚洲欧美日韩卡通动漫| 日本爱情动作片www.在线观看| 久久综合国产亚洲精品| 亚洲婷婷狠狠爱综合网| 欧美 亚洲 国产 日韩一| 欧美 亚洲 国产 日韩一| 亚洲欧美一区二区三区国产| 少妇人妻精品综合一区二区| 亚洲精品乱码久久久v下载方式| 男人添女人高潮全过程视频| 女人久久www免费人成看片| 高清午夜精品一区二区三区| 最近中文字幕2019免费版| 一级,二级,三级黄色视频| 日韩一本色道免费dvd| 国产精品成人在线| 少妇丰满av| 一本色道久久久久久精品综合| 午夜激情久久久久久久| 日韩强制内射视频| 亚洲,一卡二卡三卡| 国语对白做爰xxxⅹ性视频网站| 丝袜脚勾引网站| 妹子高潮喷水视频| videos熟女内射| 国产午夜精品一二区理论片| 看十八女毛片水多多多| 亚洲成色77777| 亚洲电影在线观看av| 国内少妇人妻偷人精品xxx网站| 日本午夜av视频| 亚洲一级一片aⅴ在线观看| 97在线视频观看| 国产黄片视频在线免费观看| 一级毛片我不卡| 水蜜桃什么品种好| 纯流量卡能插随身wifi吗| 免费黄频网站在线观看国产| 国产伦精品一区二区三区四那| 99视频精品全部免费 在线| 国产精品秋霞免费鲁丝片| 日日摸夜夜添夜夜爱| 特大巨黑吊av在线直播| 99re6热这里在线精品视频| 亚洲欧美精品专区久久| 我要看日韩黄色一级片| 国产在视频线精品| 国产无遮挡羞羞视频在线观看| 中文字幕人妻熟人妻熟丝袜美| 在线播放无遮挡| 欧美最新免费一区二区三区| 九九在线视频观看精品| 大话2 男鬼变身卡| 如日韩欧美国产精品一区二区三区 | 26uuu在线亚洲综合色| 国产男人的电影天堂91| 精品一区在线观看国产| 久久韩国三级中文字幕| 亚洲国产精品成人久久小说| 香蕉精品网在线| 麻豆成人av视频| 街头女战士在线观看网站| 亚洲av男天堂| 久久97久久精品| 在线观看www视频免费| 亚洲av成人精品一区久久| 高清毛片免费看| 狠狠精品人妻久久久久久综合| 国产精品一区www在线观看| 国产高清国产精品国产三级| 免费久久久久久久精品成人欧美视频 | 我要看黄色一级片免费的| 性色av一级| 国产一区有黄有色的免费视频| 国产伦在线观看视频一区| av天堂中文字幕网| 亚洲真实伦在线观看| 美女国产视频在线观看| 噜噜噜噜噜久久久久久91| 欧美日韩国产mv在线观看视频| 色婷婷久久久亚洲欧美| 久久99热6这里只有精品| 亚洲国产精品国产精品| 国产成人午夜福利电影在线观看| 成人免费观看视频高清| 久久99一区二区三区| 亚洲美女黄色视频免费看| 亚洲国产欧美在线一区| 综合色丁香网| 一本一本综合久久| 免费看日本二区| 久久狼人影院| 精品国产一区二区三区久久久樱花| 日本av免费视频播放| 久久99热这里只频精品6学生| 一级av片app| 国产av一区二区精品久久| 国产高清不卡午夜福利| 纵有疾风起免费观看全集完整版| 视频区图区小说| 丰满乱子伦码专区| 最后的刺客免费高清国语| 精品人妻熟女av久视频| 香蕉精品网在线| 国产欧美日韩精品一区二区| 欧美激情极品国产一区二区三区 | 亚洲精品乱码久久久v下载方式| 人人妻人人看人人澡| 欧美日本中文国产一区发布| 交换朋友夫妻互换小说| 欧美精品一区二区免费开放| 亚洲av男天堂| 一级毛片久久久久久久久女| 国产精品久久久久成人av| 丝袜在线中文字幕| 免费少妇av软件| 伊人久久精品亚洲午夜| 黑丝袜美女国产一区| 亚洲精品久久久久久婷婷小说| 七月丁香在线播放| 午夜免费鲁丝| 欧美日韩av久久| 97在线人人人人妻| 秋霞伦理黄片| 精品国产露脸久久av麻豆| av播播在线观看一区| 日韩中文字幕视频在线看片| 久久97久久精品| 99热这里只有是精品在线观看| 精品亚洲成a人片在线观看| 麻豆成人av视频| 我的老师免费观看完整版| 国产精品偷伦视频观看了| 免费不卡的大黄色大毛片视频在线观看| 黄色怎么调成土黄色| 亚洲国产精品成人久久小说| 乱人伦中国视频| 欧美老熟妇乱子伦牲交| 国产极品粉嫩免费观看在线 | av不卡在线播放| 在线观看免费视频网站a站| 国产精品成人在线| 嫩草影院入口| 久久精品熟女亚洲av麻豆精品| 久久久午夜欧美精品| 亚洲情色 制服丝袜| 51国产日韩欧美| 久久精品国产a三级三级三级| 少妇精品久久久久久久| 日韩熟女老妇一区二区性免费视频| 日产精品乱码卡一卡2卡三| 97超视频在线观看视频| 一区二区三区四区激情视频| 99热全是精品| tube8黄色片| 岛国毛片在线播放| 22中文网久久字幕| 极品人妻少妇av视频| 亚洲精品国产成人久久av| 国产精品久久久久成人av| 又爽又黄a免费视频| 欧美老熟妇乱子伦牲交| 国产精品久久久久久久电影| 亚洲天堂av无毛| 国产男人的电影天堂91| 丰满少妇做爰视频| 亚洲精品一区蜜桃| 亚洲精品自拍成人| 亚洲高清免费不卡视频| 久久亚洲国产成人精品v| 日韩人妻高清精品专区| 2018国产大陆天天弄谢| 久久青草综合色| 日韩 亚洲 欧美在线| 99热这里只有是精品50| 成人免费观看视频高清| 国产真实伦视频高清在线观看| 一边亲一边摸免费视频| 久久97久久精品| 免费看av在线观看网站| h视频一区二区三区| 99视频精品全部免费 在线| 国产成人精品婷婷| 免费黄网站久久成人精品| 永久网站在线| av黄色大香蕉| 亚洲精品aⅴ在线观看| 人人妻人人添人人爽欧美一区卜| 中文字幕免费在线视频6| 少妇高潮的动态图| 一区在线观看完整版| 亚洲精品国产成人久久av| 美女cb高潮喷水在线观看| 国产高清三级在线| 免费av不卡在线播放| 一区二区三区精品91| 国产精品免费大片| 亚洲综合精品二区| 亚洲av在线观看美女高潮| 国产精品一区二区三区四区免费观看| 美女主播在线视频| 国产免费视频播放在线视频| 日韩欧美一区视频在线观看 | 简卡轻食公司| 美女中出高潮动态图| 永久免费av网站大全| 精品久久国产蜜桃| 国产在线一区二区三区精| 久久久亚洲精品成人影院| 免费观看无遮挡的男女| 国产在线男女| 亚州av有码| 免费看av在线观看网站| 久久久亚洲精品成人影院| 在线观看免费高清a一片| 秋霞伦理黄片| 日韩不卡一区二区三区视频在线| 国产视频内射| 天天操日日干夜夜撸| av在线播放精品| 精品视频人人做人人爽| 国产成人91sexporn| 一级黄片播放器| 精品国产一区二区三区久久久樱花| 成人18禁高潮啪啪吃奶动态图 | 久久婷婷青草| 十分钟在线观看高清视频www | 欧美精品人与动牲交sv欧美| 一区二区三区乱码不卡18| 亚洲无线观看免费| 夫妻午夜视频| 久久精品国产鲁丝片午夜精品| 国产极品天堂在线| 丝袜脚勾引网站| 久久女婷五月综合色啪小说| 嫩草影院新地址| 青春草视频在线免费观看| av福利片在线| 热re99久久精品国产66热6| 亚洲欧美精品自产自拍| 精品国产国语对白av| 午夜免费男女啪啪视频观看| 大片电影免费在线观看免费| 一个人免费看片子| 久久精品国产自在天天线| 亚洲欧美一区二区三区国产| 国产免费福利视频在线观看| 一本久久精品| 黄色日韩在线| 久久久久久久精品精品| h视频一区二区三区| 精品一区二区免费观看| 国产欧美日韩一区二区三区在线 | 久久久精品94久久精品| 亚洲三级黄色毛片| 一本色道久久久久久精品综合| 国产精品蜜桃在线观看| 欧美日韩在线观看h| 欧美日韩av久久| 99精国产麻豆久久婷婷| 我要看日韩黄色一级片| 久久久精品免费免费高清| 日韩制服骚丝袜av| 纯流量卡能插随身wifi吗| 亚洲av中文av极速乱| 一区二区三区乱码不卡18| 国产在线视频一区二区| 久久综合国产亚洲精品| 18+在线观看网站| 亚洲精品色激情综合| 国产91av在线免费观看| 久久久午夜欧美精品| 热99国产精品久久久久久7| 看十八女毛片水多多多| 啦啦啦中文免费视频观看日本| 国产欧美日韩精品一区二区| 久久人人爽人人爽人人片va| 一级爰片在线观看| 我的女老师完整版在线观看| 高清在线视频一区二区三区| 伊人亚洲综合成人网| 成人18禁高潮啪啪吃奶动态图 | 一二三四中文在线观看免费高清| 九九爱精品视频在线观看| 国产探花极品一区二区| 少妇丰满av| 最近手机中文字幕大全| 久久久精品免费免费高清| 最近最新中文字幕免费大全7| 国产一级毛片在线| 国产极品天堂在线| 亚洲精品成人av观看孕妇| 日本黄大片高清| 欧美性感艳星| 国产欧美日韩一区二区三区在线 | 国产精品99久久久久久久久| 欧美日韩精品成人综合77777| 国精品久久久久久国模美| 十八禁高潮呻吟视频 | 国内精品宾馆在线| 亚洲人成网站在线观看播放| 这个男人来自地球电影免费观看 | 自拍偷自拍亚洲精品老妇| 亚洲欧美清纯卡通| 99视频精品全部免费 在线| 免费观看无遮挡的男女| 国产一区二区三区综合在线观看 | 又爽又黄a免费视频| 女人久久www免费人成看片| 大码成人一级视频| 99热全是精品| 伦理电影大哥的女人| 我的女老师完整版在线观看| 亚洲精品国产av成人精品| 国产伦精品一区二区三区四那| 久久国产乱子免费精品| 插阴视频在线观看视频| 精品一品国产午夜福利视频| 69精品国产乱码久久久| 插阴视频在线观看视频| 人人妻人人爽人人添夜夜欢视频 | 亚洲精品国产av成人精品| 97超碰精品成人国产| 亚洲欧美成人综合另类久久久| av女优亚洲男人天堂| 伊人久久精品亚洲午夜| 日韩av免费高清视频| 在线观看www视频免费| 不卡视频在线观看欧美| 99热国产这里只有精品6| 久久午夜综合久久蜜桃| 青春草国产在线视频| 欧美最新免费一区二区三区| 伦理电影免费视频| 久久久久久久亚洲中文字幕| 国产精品福利在线免费观看| videossex国产| a级片在线免费高清观看视频| 国产精品秋霞免费鲁丝片| 曰老女人黄片| 亚洲久久久国产精品| 看非洲黑人一级黄片| 99久国产av精品国产电影| 99国产精品免费福利视频| 在线观看三级黄色| 婷婷色av中文字幕| 你懂的网址亚洲精品在线观看| 丰满乱子伦码专区| 久久久久久人妻| 欧美xxⅹ黑人| 国产综合精华液| 久久久a久久爽久久v久久| 观看av在线不卡| 另类精品久久| 国产一区二区在线观看日韩| 极品教师在线视频| 亚洲av成人精品一二三区| 精品99又大又爽又粗少妇毛片| 少妇被粗大猛烈的视频| 99久久人妻综合| 国产伦精品一区二区三区视频9| 久久久久久久大尺度免费视频| 午夜av观看不卡| 国产精品久久久久久精品古装| 欧美精品一区二区免费开放| 18禁动态无遮挡网站| 日本色播在线视频| 纯流量卡能插随身wifi吗| 少妇被粗大猛烈的视频| 亚洲图色成人| 18禁裸乳无遮挡动漫免费视频| 国产亚洲av片在线观看秒播厂| 最近最新中文字幕免费大全7| av在线app专区| 亚洲av男天堂| 日本欧美视频一区| 国产亚洲午夜精品一区二区久久| 国产又色又爽无遮挡免| 一级毛片久久久久久久久女| 国产在线免费精品| 18禁在线播放成人免费| 成人美女网站在线观看视频| 成人毛片a级毛片在线播放| 性高湖久久久久久久久免费观看| 我的女老师完整版在线观看| 人妻系列 视频| 成人亚洲欧美一区二区av| 国产精品一区二区在线不卡| 精品亚洲成a人片在线观看| 国产精品久久久久久久电影| 亚洲欧美日韩卡通动漫| av专区在线播放| 国产一区二区在线观看av| 最近的中文字幕免费完整| 亚洲av二区三区四区| 我要看日韩黄色一级片| 欧美精品高潮呻吟av久久| 婷婷色综合大香蕉| 黄色视频在线播放观看不卡| 99精国产麻豆久久婷婷| 在线观看美女被高潮喷水网站| 欧美+日韩+精品| 久久99热这里只频精品6学生| 国产成人精品福利久久| 亚洲va在线va天堂va国产| 寂寞人妻少妇视频99o| 国产老妇伦熟女老妇高清| 欧美精品一区二区大全| 国产精品一二三区在线看| 深夜a级毛片| 特大巨黑吊av在线直播| 久久久久久伊人网av| 午夜影院在线不卡| 激情五月婷婷亚洲| 免费看光身美女| 91精品伊人久久大香线蕉| 老司机影院毛片| 美女cb高潮喷水在线观看| 日韩欧美 国产精品| 亚洲精品成人av观看孕妇| 91aial.com中文字幕在线观看| 18禁动态无遮挡网站| 精品久久久久久电影网| 亚洲三级黄色毛片| 成人漫画全彩无遮挡| 人妻制服诱惑在线中文字幕| 插阴视频在线观看视频| 美女大奶头黄色视频| 免费在线观看成人毛片| 久久精品国产自在天天线| 国产乱人偷精品视频| 精品午夜福利在线看| 国产在线男女| 黄色怎么调成土黄色| h日本视频在线播放| 国产在视频线精品| 日韩av不卡免费在线播放| 三级国产精品片| 亚洲天堂av无毛| 久久精品国产a三级三级三级| 免费不卡的大黄色大毛片视频在线观看| 精品久久久噜噜| 国产成人午夜福利电影在线观看| 黑人巨大精品欧美一区二区蜜桃 | 一级a做视频免费观看| 免费黄频网站在线观看国产| www.色视频.com| 亚洲图色成人| 看免费成人av毛片| 久久精品久久精品一区二区三区| www.av在线官网国产| av在线播放精品| 国产在线免费精品| 久久97久久精品| 国产在线免费精品| 久久 成人 亚洲| 国产精品一区二区三区四区免费观看| 日韩成人av中文字幕在线观看| 涩涩av久久男人的天堂| 少妇丰满av| 亚洲美女搞黄在线观看| 日本午夜av视频| 国产欧美日韩一区二区三区在线 | 多毛熟女@视频| 午夜免费男女啪啪视频观看| 亚洲av电影在线观看一区二区三区| 国产伦在线观看视频一区| 国产精品免费大片| h视频一区二区三区| 黑人猛操日本美女一级片| 两个人免费观看高清视频 | 中文字幕亚洲精品专区| 免费人妻精品一区二区三区视频| 一级黄片播放器| 亚洲一区二区三区欧美精品| 女性被躁到高潮视频| av.在线天堂| 亚洲精品国产av成人精品| 欧美区成人在线视频| 99九九在线精品视频 | av国产久精品久网站免费入址| 夫妻性生交免费视频一级片| 亚洲婷婷狠狠爱综合网| 最近中文字幕2019免费版| 美女福利国产在线| videos熟女内射| 看非洲黑人一级黄片| 国产精品蜜桃在线观看| 91aial.com中文字幕在线观看| 一级黄片播放器| 久久久亚洲精品成人影院| 亚洲欧美中文字幕日韩二区| 亚洲av日韩在线播放| 日本欧美国产在线视频| 黄色一级大片看看| 午夜免费鲁丝| 久久久久久久久久久免费av| 人妻少妇偷人精品九色| 日韩伦理黄色片| 美女中出高潮动态图| 校园人妻丝袜中文字幕| 夜夜看夜夜爽夜夜摸| 青春草国产在线视频| 亚洲国产精品专区欧美| 18+在线观看网站| 精品少妇内射三级| 男女边摸边吃奶| 精品99又大又爽又粗少妇毛片| 国产亚洲91精品色在线| 亚洲国产日韩一区二区| 国产老妇伦熟女老妇高清| 日本av免费视频播放| 观看免费一级毛片| 人人澡人人妻人| 欧美日韩综合久久久久久| 在线免费观看不下载黄p国产| 日韩不卡一区二区三区视频在线| 日韩欧美一区视频在线观看 | 亚洲欧美精品自产自拍| 免费看日本二区| 精品午夜福利在线看| 精品少妇黑人巨大在线播放| 亚洲va在线va天堂va国产| 免费高清在线观看视频在线观看| 涩涩av久久男人的天堂| 欧美性感艳星| 中文字幕人妻丝袜制服| 婷婷色av中文字幕| 国产男人的电影天堂91| 国产欧美亚洲国产| 国产 精品1| 最新中文字幕久久久久| 人妻 亚洲 视频| 亚洲成人手机| 亚洲精品视频女| 国产乱人偷精品视频| av天堂久久9| 99视频精品全部免费 在线| a级毛色黄片| 久久久久久久久大av| 欧美老熟妇乱子伦牲交| 中文字幕精品免费在线观看视频 | 22中文网久久字幕| 老熟女久久久| 中国国产av一级| 欧美日韩一区二区视频在线观看视频在线| 欧美丝袜亚洲另类| 亚洲国产欧美在线一区| 久久鲁丝午夜福利片| 欧美日韩av久久| 一区二区av电影网| 亚洲精品一二三| 少妇猛男粗大的猛烈进出视频| 欧美日韩一区二区视频在线观看视频在线| 成年人午夜在线观看视频| 亚洲av国产av综合av卡| 在线观看www视频免费| 久久99蜜桃精品久久| 女的被弄到高潮叫床怎么办| 国产一区二区在线观看日韩| 乱人伦中国视频| 五月玫瑰六月丁香| 久久毛片免费看一区二区三区| 国产日韩一区二区三区精品不卡 | tube8黄色片| 日本猛色少妇xxxxx猛交久久| a级一级毛片免费在线观看| 少妇裸体淫交视频免费看高清| 久久久久久久久久人人人人人人| 热99国产精品久久久久久7| 亚洲精品久久久久久婷婷小说| 热re99久久精品国产66热6| 欧美bdsm另类| 蜜桃久久精品国产亚洲av| 视频区图区小说| 免费黄频网站在线观看国产| 日韩av免费高清视频| 亚洲国产精品999| 在线观看免费高清a一片| 免费观看av网站的网址| 国产 精品1| 国产一区二区三区av在线| 丝瓜视频免费看黄片| 亚洲精品色激情综合| 在线精品无人区一区二区三| 亚洲av免费高清在线观看| 久久久久精品久久久久真实原创| 赤兔流量卡办理| 春色校园在线视频观看| av福利片在线| 欧美高清成人免费视频www| 国语对白做爰xxxⅹ性视频网站| 久久国产乱子免费精品| 日韩欧美一区视频在线观看 | 午夜日本视频在线| 精品亚洲成a人片在线观看| 国产黄片视频在线免费观看| 男女边摸边吃奶| 欧美精品高潮呻吟av久久| 国产av国产精品国产| 久久久午夜欧美精品| 国产免费又黄又爽又色| 六月丁香七月| 国产黄频视频在线观看| 日韩精品免费视频一区二区三区 | 日韩大片免费观看网站|