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

    南水北調(diào)中線工程水源區(qū)“三生”空間演化特征與生態(tài)效應(yīng)

    2024-02-10 00:00:00段亞明付景保周翼
    中國人口·資源與環(huán)境 2024年11期

    摘要 “三生”空間演化及其生態(tài)效應(yīng)研究可為國土空間規(guī)劃與生態(tài)環(huán)境保護提供決策參考。南水北調(diào)中線工程是緩解華北地區(qū)水資源短缺的重大戰(zhàn)略性基礎(chǔ)工程,水源區(qū)生態(tài)環(huán)境穩(wěn)定是保障水質(zhì)的基礎(chǔ)。為此,該研究立足生態(tài)系統(tǒng)服務(wù)視角,在開展南水北調(diào)中線工程水源區(qū)“三生”空間識別與生態(tài)系統(tǒng)服務(wù)評估的基礎(chǔ)上,分析水源區(qū)“三生”空間演化特征及其生態(tài)效應(yīng)。結(jié)果表明:①南水北調(diào)中線工程水源區(qū)內(nèi)各類“三生”空間的空間格局基本穩(wěn)定,除農(nóng)業(yè)生產(chǎn)空間、草地空間外,其他空間的面積均有增長,其中林地空間的增量最多,非農(nóng)生產(chǎn)空間增幅最大。不同時期內(nèi),農(nóng)業(yè)生產(chǎn)空間、林地空間與草地空間的轉(zhuǎn)換面積均遠高于其他空間。②水源區(qū)整體及各類空間內(nèi)的生態(tài)系統(tǒng)服務(wù)相對穩(wěn)定,林地空間、草地空間與農(nóng)業(yè)生產(chǎn)空間的生態(tài)系統(tǒng)服務(wù)水平顯著高于其他空間。③“三生”空間轉(zhuǎn)換對水源區(qū)生態(tài)系統(tǒng)服務(wù)存在明顯的正向與負向效應(yīng),其中,農(nóng)業(yè)生產(chǎn)空間轉(zhuǎn)為林地空間、草地空間以及草地空間轉(zhuǎn)為林地空間對水源涵養(yǎng)與水土保持服務(wù)增長的貢獻率較高。④南水北調(diào)中線工程水源區(qū)“三生”空間演化的生態(tài)效應(yīng)與人類活動密切相關(guān),農(nóng)業(yè)生產(chǎn)、城鄉(xiāng)建設(shè)、水利建設(shè)以及生態(tài)保護等人類活動通過改變空間利用方式推動生態(tài)系統(tǒng)服務(wù)演化發(fā)展。其中,農(nóng)業(yè)生產(chǎn)、城鄉(xiāng)建設(shè)與水利建設(shè)活動對水源區(qū)生態(tài)系統(tǒng)服務(wù)具有負向效應(yīng),生態(tài)保護活動則發(fā)揮正向效應(yīng)。未來,政府應(yīng)進一步推進南水北調(diào)中線工程水源區(qū)全域土地綜合整治,綜合統(tǒng)籌水源區(qū)內(nèi)農(nóng)業(yè)生產(chǎn)與生態(tài)保護的整體格局。

    關(guān)鍵詞 “三生”空間;演化特征;生態(tài)效應(yīng);生態(tài)系統(tǒng)服務(wù)評估;南水北調(diào)中線工程水源區(qū)

    中圖分類號 F205;X321 文獻標志碼 A 文章編號 1002-2104(2024)11-0138-13 DOI:10. 12062/cpre. 20240902

    國土空間是人類生存與發(fā)展的重要載體[1]。改革開放后,中國經(jīng)濟社會快速發(fā)展,工業(yè)化、城鎮(zhèn)化進程持續(xù)推進,國土空間的開發(fā)與利用程度也不斷增強,并引發(fā)交通擁堵、環(huán)境污染、生態(tài)退化等諸多負面問題[2-4]。為指導(dǎo)國土空間有序發(fā)展、緩解人類活動對國土空間的擾動,中國立足本國國情提出了“生產(chǎn)空間集約高效、生活空間宜居適度、生態(tài)空間山清水秀”的國土空間發(fā)展理念[5]。在現(xiàn)實需求與政策引領(lǐng)的共同作用下,“三生”空間研究作為探明國土空間發(fā)展規(guī)律、優(yōu)化國土空間開發(fā)格局與實現(xiàn)經(jīng)濟-社會協(xié)調(diào)發(fā)展的重要手段,成為地理學(xué)、規(guī)劃學(xué)等學(xué)科探討的研究熱點之一[6-8]。

    縱觀已有研究,“三生”空間研究主要涉及“三生”空間基礎(chǔ)理論[9-10]、分類與識別[11-12]、演化與模擬[13-14]、沖突與優(yōu)化[15-16]等方面。隨著生態(tài)文明戰(zhàn)略的持續(xù)推進,“三生”空間演化的生態(tài)效應(yīng)研究也逐漸受到關(guān)注[17-18]。現(xiàn)有研究主要基于不同空間類型的生態(tài)環(huán)境質(zhì)量指數(shù),通過計算不同區(qū)域國土空間的生態(tài)質(zhì)量變化與各類空間轉(zhuǎn)換的生態(tài)貢獻率,分析流域、市、縣等不同尺度內(nèi)“三生”空間演化的生態(tài)效應(yīng)[19-21]。例如,黃晶等[22]、董冬等[23]均基于生態(tài)質(zhì)量指數(shù)法分別測度并分析了典型干旱綠洲農(nóng)業(yè)區(qū)以及長三角城市群“三生”空間演化的生態(tài)效應(yīng)與驅(qū)動力。雖然該方法具有直觀、可操作性強等特點,但是由于各類空間的生態(tài)環(huán)境質(zhì)量指數(shù)多基于專家打分法等主觀方式獲得,難以客觀反映區(qū)域內(nèi)及同類型“三生”空間內(nèi)部生態(tài)質(zhì)量的空間異質(zhì)性,降低了研究結(jié)果的說服力。生態(tài)系統(tǒng)服務(wù)評估是基于生態(tài)學(xué)原理量化生態(tài)系統(tǒng)對人類提供的各種惠益,是反映生態(tài)環(huán)境質(zhì)量的重要參考,在制定生態(tài)補償策略、環(huán)境評估和生態(tài)保護區(qū)劃定中應(yīng)用廣泛[24-25]。相對于生態(tài)環(huán)境質(zhì)量指數(shù)法,生態(tài)系統(tǒng)服務(wù)評估可以定量測度不同空間內(nèi)生態(tài)系統(tǒng)服務(wù)的高低,進而準確反映“三生”空間內(nèi)部的生態(tài)質(zhì)量,彌補了生態(tài)環(huán)境質(zhì)量指數(shù)法的不足,為“三生”空間演化的生態(tài)效應(yīng)研究提供了新的方法與手段[26-27]。因此,基于生態(tài)系統(tǒng)服務(wù)視角分析“三生”空間演化的生態(tài)效應(yīng),更有助于準確揭示國土空間開發(fā)對生態(tài)環(huán)境的影響,對于經(jīng)濟社會協(xié)調(diào)發(fā)展、改善生態(tài)環(huán)境與實施國土空間規(guī)劃具有重要參考意義,但目前基于生態(tài)系統(tǒng)服務(wù)視角探究“三生”空間演化的生態(tài)效應(yīng)的研究相對較少[28-29]。

    干凈、安全的水資源是人類生存與發(fā)展的基礎(chǔ),也是聯(lián)合國《2030年可持續(xù)發(fā)展議程》制定的可持續(xù)發(fā)展目標之一。水源區(qū)生態(tài)環(huán)境質(zhì)量是維持、凈化水資源的基本保障,而分析水源區(qū)“三生”空間演化及其生態(tài)效應(yīng)對認識水源區(qū)生態(tài)環(huán)境現(xiàn)狀、制定國土空間規(guī)劃與開展生態(tài)環(huán)境保護具有重要借鑒意義,但現(xiàn)有研究對水源區(qū)“三生”空間演化的生態(tài)效應(yīng)及其驅(qū)動機制缺乏關(guān)注[30]。南水北調(diào)中線工程是中國重要的跨流域、跨區(qū)域調(diào)水工程,旨在緩解中國華北地區(qū)水資源嚴重短缺問題,自2014年通水以來累計調(diào)水600億m3,對保障京津冀等工程沿線地區(qū)的生態(tài)安全和水安全具有重要意義[31]。南水北調(diào)中線工程水源區(qū)包括丹江口水庫庫區(qū)以及漢江、丹江的部分支流,整體位于伏牛山南麓與大巴山之間,具有重要的水源涵養(yǎng)與水土保持功能,區(qū)域內(nèi)生態(tài)環(huán)境狀況對水質(zhì)具有重要支撐作用[32]。自中線工程開工以來,各級政府不斷通過防護林建設(shè)、濕地修復(fù)等手段改善水源區(qū)生態(tài)環(huán)境,但是隨著中部崛起戰(zhàn)略的實施,水源區(qū)內(nèi)城鎮(zhèn)化、工業(yè)化進程逐步加快,生產(chǎn)、生活、生態(tài)空間的轉(zhuǎn)換日趨激烈并對區(qū)域內(nèi)生態(tài)環(huán)境產(chǎn)生巨大影響。

    為此,本研究立足生態(tài)系統(tǒng)服務(wù)視角,基于遙感、氣象等多源數(shù)據(jù)與GIS空間分析、InVEST模型等多學(xué)科方法,在開展南水北調(diào)中線工程水源區(qū)“三生”空間識別與生態(tài)系統(tǒng)服務(wù)評估的基礎(chǔ)上,系統(tǒng)分析水源區(qū)“三生”空間演化特征及其生態(tài)環(huán)境效應(yīng),以期為南水北調(diào)中線工程水源區(qū)“三生”空間協(xié)調(diào)發(fā)展與國土空間可持續(xù)利用提供科學(xué)參考。

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

    1. 1 研究區(qū)概況

    南水北調(diào)中線工程水源區(qū)(105°28'48″ E~112°22'12″E、31°27'36″ N~34°25'48″ N)指丹江口庫區(qū)及其上游地區(qū),位于中國第二、三階梯的過渡地帶,地勢由南、北向中部傾斜,具有平原、山地等多種地貌類型,北部以秦嶺支脈伏牛山南麓為界,南部以大巴山—武當山為界,中部為漢水谷地與南陽盆地,海拔高度在80~3 056 m之間,屬北亞熱帶向暖溫帶過渡氣候區(qū),涉及河南、湖北、陜西等6個省份的49個區(qū)(縣)??紤]研究數(shù)據(jù)的可得性與水源區(qū)覆蓋的完整性,本研究選取漢中、安康、商洛、十堰以及南陽市下屬4個區(qū)(縣)為研究區(qū),占水源區(qū)總面積的81. 55%(圖1)。2000年后,南水北調(diào)中線工程開工建設(shè),丹江口大壩也由162 m加高至176. 6 m,庫區(qū)水位進一步抬升。為提高水源區(qū)內(nèi)水源涵養(yǎng)能力、降低水土流失風(fēng)險,各級政府積極在水源區(qū)內(nèi)開展植樹造林、退耕還林還草、護坡草帶建設(shè)等多種生態(tài)工程。與此同時,水源區(qū)內(nèi)城鄉(xiāng)建設(shè)與經(jīng)濟發(fā)展同步推進,截至2020年,第一、二、三產(chǎn)業(yè)產(chǎn)值分別達到954. 72億元、2 565. 11億元與3 012. 91億元,區(qū)域總?cè)丝? 553. 47萬人。經(jīng)濟發(fā)展、城鄉(xiāng)建設(shè)與南水北調(diào)中線工程持續(xù)推進使水源區(qū)“三生”空間發(fā)生劇烈變化,也使水源區(qū)生態(tài)環(huán)境發(fā)生改變。因此,揭示南水北調(diào)中線工程水源區(qū)“三生”空間演化特征及其生態(tài)效應(yīng)對保護水源區(qū)生態(tài)環(huán)境、確保南水北調(diào)中線工程穩(wěn)定運行具有重要意義。

    1. 2 數(shù)據(jù)來源

    本研究所使用的基礎(chǔ)數(shù)據(jù)包括南水北調(diào)中線工程水源區(qū)降雨數(shù)據(jù)、遙感影像、數(shù)字高程模型(DEM)、行政邊界數(shù)據(jù)、城鎮(zhèn)建成區(qū)邊界數(shù)據(jù)、社會經(jīng)濟統(tǒng)計數(shù)據(jù)等。其中,降雨數(shù)據(jù)來自國家地球系統(tǒng)科學(xué)數(shù)據(jù)中心(http://www. geodata. cn/)。2000、2010與2020年3期南水北調(diào)中線工程水源區(qū)Landsat遙感影像(均通過軟件ENVI5. 3進行輻射定標、大氣校正等預(yù)處理工作,空間分辨率30 m)與DEM高程數(shù)據(jù),均來自地理空間數(shù)據(jù)云(https://www.gscloud. cn/)。2000、2010與2020年的城鎮(zhèn)建成區(qū)邊界數(shù)據(jù)來自國家對地觀測科學(xué)數(shù)據(jù)中心(https://www. chinageoss.cn/)。漢中、安康、商洛、十堰與南陽的社會經(jīng)濟統(tǒng)計數(shù)據(jù)均來自各市統(tǒng)計局。

    2 研究方法

    2. 1 “三生”空間分類體系與識別方法

    “三生”空間分類與識別是開展“三生”空間演化研究的基礎(chǔ)?!叭笨臻g作為一種綜合性的國土空間分區(qū)方式,包含生產(chǎn)、生活、生態(tài)3種大類空間[33]。其中,生產(chǎn)空間指為人類提供各類產(chǎn)品的空間;生活空間指為人類提供各類生活性服務(wù)的空間;生態(tài)空間指對維持生物生存環(huán)境穩(wěn)定與生態(tài)系統(tǒng)平衡具有重要作用的空間[34-36]。由于各種大類空間內(nèi)的生產(chǎn)活動、土地權(quán)屬或自然景觀存在差異,導(dǎo)致各大類空間內(nèi)存在不同的子空間[37]。例如,城鎮(zhèn)生活空間、鄉(xiāng)村生活空間等。而各種“三生”空間子空間往往具有不同的演化特征,分析各類子空間的演化特征可以為城鄉(xiāng)空間、各類產(chǎn)業(yè)空間與生態(tài)空間的管控提供更精準的決策支持[37]。因此,在參考已有研究[38-40]的基礎(chǔ)上,遵循系統(tǒng)性、可行性與實用性等原則,構(gòu)建“三生”空間分類體系(表1)。首先,根據(jù)空間的主導(dǎo)功能將“三生”空間劃分為生產(chǎn)、生活與生態(tài)空間。然后,根據(jù)生產(chǎn)活動、土地權(quán)屬與自然景觀差異,分別對生產(chǎn)、生活與生態(tài)空間進行細化。根據(jù)內(nèi)部生產(chǎn)活動不同,生產(chǎn)空間被劃分為農(nóng)業(yè)生產(chǎn)空間與非農(nóng)生產(chǎn)空間,分別支撐農(nóng)業(yè)與非農(nóng)業(yè)生產(chǎn)活動。根據(jù)土地權(quán)屬差異,生活空間被劃分為城鎮(zhèn)生活空間與鄉(xiāng)村生活空間,其土地權(quán)屬分別由國家與農(nóng)村集體所有。根據(jù)自然景觀差異,生態(tài)空間被劃分為林地空間、草地空間、水域空間與裸地空間,分別以林地、草地、水域與裸土為主導(dǎo)景觀。

    在識別方法上,首先,以2000、2010與2020年3期南水北調(diào)中線工程水源區(qū)遙感影像為基礎(chǔ),基于隨機森林算法,提取水源區(qū)農(nóng)業(yè)生產(chǎn)空間、林地空間、草地空間、水域空間、裸地空間、其他空間(包括:非農(nóng)生產(chǎn)空間與城、鄉(xiāng)生活空間),總體識別精度分別為82. 51%、83. 62%、85. 52%。然后,根據(jù)Google Earth高分辨率遙感影像,對空間識別結(jié)果進行人工目視修正。最后,結(jié)合GoogleEarth高分辨率遙感影像、城鎮(zhèn)建成區(qū)邊界數(shù)據(jù)與實地調(diào)研,對3個時期的非農(nóng)生產(chǎn)空間、城鎮(zhèn)生活空間與鄉(xiāng)村生活空間進行目視識別(圖2)。

    2. 2 “三生”空間轉(zhuǎn)移矩陣

    “三生”空間轉(zhuǎn)移矩陣能夠精確刻畫不同類型“三生”空間轉(zhuǎn)換的數(shù)量、方向與空間位置,從而揭示“三生”空間的演化格局[41-42]。本研究利用軟件ArcGIS10. 2中的柵格計算器工具得到2000—2010年與2010—2020年南水北調(diào)中線工程水源區(qū)“ 三生”空間轉(zhuǎn)移矩陣。計算公式如下:

    2. 3 生態(tài)系統(tǒng)服務(wù)評估

    生態(tài)系統(tǒng)服務(wù)指生態(tài)系統(tǒng)與生態(tài)過程對維持人類賴以生存的自然環(huán)境條件中發(fā)揮的作用[43]。生態(tài)系統(tǒng)服務(wù)評估則是基于生態(tài)系統(tǒng)的過程、功能機制以及與人類的互動關(guān)系對生態(tài)系統(tǒng)服務(wù)的價值量或物理量進行定量測度[44]。根據(jù)南水北調(diào)中線工程的功能定位與《丹江口庫區(qū)及上游水污染防治和水土保持“十四五”規(guī)劃》,水源區(qū)具有重要的水源涵養(yǎng)與水土保持功能。InVEST模型由斯坦福大學(xué)研發(fā),可對陸地、海洋等多種生態(tài)系統(tǒng)開展生態(tài)系統(tǒng)服務(wù)評估,具有運算速度快、數(shù)據(jù)兼容性強與操作簡單等特點,在生態(tài)系統(tǒng)服務(wù)評估中受到廣泛應(yīng)用[45-47]。因此,本研究主要采用InVEST模型對南水北調(diào)中線工程水源區(qū)的水源涵養(yǎng)與水土保持服務(wù)進行測度。

    2. 3. 1 水源涵養(yǎng)服務(wù)

    水源涵養(yǎng)服務(wù)評估旨在測度不同空間對降雨進行截留、滲透與蓄積的能力,通常以水源涵養(yǎng)量衡量不同區(qū)域的水源涵養(yǎng)能力[48]。本研究基于水量平衡原理對水源區(qū)的水源涵養(yǎng)服務(wù)進行測算[46]。首先,基于InVEST 模型Water Yield模塊計算區(qū)域的產(chǎn)水量(WY),然后根據(jù)流速系數(shù)、地形指數(shù)、土壤飽和導(dǎo)水率等參數(shù)計算水源涵養(yǎng)量。具體公式如下:

    Water Yield模塊運行的關(guān)鍵參數(shù)包括:年總降雨量、年潛在蒸散量、土壤的最大根系埋藏深度、植物可利用水量、生物物理系數(shù)表、Z系數(shù)。其中,年總降雨量、年潛在蒸散量、土壤的最大根系埋藏深度分別通過氣象數(shù)據(jù)與土壤數(shù)據(jù)直接獲得。植物可利用水量(PAWC)的計算見公式4。在參考現(xiàn)有研究[49-50]、《InVEST模型使用手冊》以及開展專家咨詢的基礎(chǔ)上,將Z 系數(shù)設(shè)置為8. 6,并制作用于Water Yield模塊運行的生物物理系數(shù)表(表2)。

    2. 3. 2 水土保持服務(wù)

    水土保持服務(wù)評價旨在測算各類空間對于由水蝕導(dǎo)致的土壤侵蝕作用的減弱能力,通常以土壤保持量衡量區(qū)域的水土保持能力[51]。本研究采用InVEST 模型的SDR模塊對水源區(qū)2000、2010與2020年的水土保持服務(wù)進行計算。模型運行的關(guān)鍵參數(shù)包括:降雨侵蝕率因子、土壤可蝕性因子、生物物理系數(shù)表(表3)。其中,降雨侵蝕率與土壤可蝕性分別根據(jù)Wischmeier[52]與《生態(tài)保護紅線劃定技術(shù)指南》提供的公式計算,詳見式(5)—式(10)。在參考現(xiàn)有研究[53-54]、《InVEST模型使用手冊》以及開展專家咨詢的基礎(chǔ)上,制作用于SDR模塊運行的生物物理系數(shù)表。

    2. 4 生態(tài)貢獻率

    生態(tài)貢獻率指區(qū)域內(nèi)“三生”空間轉(zhuǎn)換對生態(tài)環(huán)境變化的貢獻程度,為識別驅(qū)動生態(tài)環(huán)境變化的主導(dǎo)因素提供科學(xué)依據(jù)。本研究在參考現(xiàn)有研究[55-57]的基礎(chǔ)上,確定貢獻率計算方法并定量測度不同空間轉(zhuǎn)換類型對生態(tài)系統(tǒng)服務(wù)變化的影響。公式如下:

    3 “三生”空間與生態(tài)系統(tǒng)服務(wù)演化分析

    3. 1 “三生”空間演化分析

    3. 1. 1 “三生”空間格局與數(shù)量變化

    根據(jù)圖2、表4,2000—2020年南水北調(diào)中線工程水源區(qū)“三生”空間格局基本穩(wěn)定,林地空間、草地空間與農(nóng)業(yè)種植空間的占比較高,其余空間的占比均未超過2%。農(nóng)業(yè)生產(chǎn)空間主要分布于漢水谷地、南陽盆地以及伏牛山與大巴山區(qū)的山谷地帶,具有集中連片的分布特征,總面積減少1 074. 32 km2。林地空間、草地空間均主要分布于伏牛山與大巴山區(qū),也具有集中連片的分布特征,但兩類空間的數(shù)量變化相反,其中林地空間增長1 648. 75 km2,草地空間減少1 651. 29 km2。水域空間以漢江水系為主體,呈線狀分布于水源區(qū)中部,總面積增長37. 45%。城鎮(zhèn)生活空間、非農(nóng)生產(chǎn)空間主要圍繞各市、區(qū)(縣)政府所在地布局,具有明顯的集聚分布特征,分別增長105. 98% 與1 303. 98%。鄉(xiāng)村生活空間主要圍繞各村級行政機構(gòu)擴張發(fā)展,空間布局相對分散,總面積增長38. 28%。裸地空間主要呈零星狀分布于水源區(qū)內(nèi),總面積僅增長8. 87 km2。

    3. 1. 2 “三生”空間轉(zhuǎn)換

    根據(jù)表5、表6,2000—2020年南水北調(diào)中線工程水源區(qū)“三生”空間發(fā)生顯著轉(zhuǎn)換,2010—2020年的總轉(zhuǎn)換面積約為2000—2010年的兩倍,農(nóng)業(yè)生產(chǎn)、城鄉(xiāng)建設(shè)、水利建設(shè)與生態(tài)保護活動是水源區(qū)“三生”空間演化的主要驅(qū)動力。不同時期內(nèi),農(nóng)業(yè)生產(chǎn)空間、林地空間與草地空間的轉(zhuǎn)換面積均遠高于其他空間。農(nóng)業(yè)生產(chǎn)空間、草地空間的轉(zhuǎn)換均受農(nóng)業(yè)生產(chǎn)與生態(tài)保護活動的雙重影響。其中,農(nóng)業(yè)生產(chǎn)空間的轉(zhuǎn)出與轉(zhuǎn)入規(guī)模相近,轉(zhuǎn)出面積略大于轉(zhuǎn)入面積,轉(zhuǎn)換類型均以林地空間、草地空間為主。而草地空間的轉(zhuǎn)出面積約為轉(zhuǎn)入面積的兩倍,轉(zhuǎn)出類型以農(nóng)業(yè)生產(chǎn)空間、林地空間為主,轉(zhuǎn)入類型則以農(nóng)業(yè)生產(chǎn)空間為主。林地空間轉(zhuǎn)換主要受生態(tài)保護活動影響,以轉(zhuǎn)入為主,且轉(zhuǎn)入面積逐階段擴大,轉(zhuǎn)入類型均以農(nóng)業(yè)生產(chǎn)空間、草地空間為主。農(nóng)業(yè)生產(chǎn)空間、林地空間與草地空間的空間轉(zhuǎn)換在水源區(qū)廣泛分布。其余空間在不同時期的轉(zhuǎn)換面積相對較小。丹江口水庫壩體加高工程帶動水源區(qū)水面提升,使水域空間面積顯著增長,轉(zhuǎn)入類型以農(nóng)業(yè)生產(chǎn)空間、林地空間與草地空間為主,主要分布于水源區(qū)中部與東部地區(qū)。城鎮(zhèn)生活空間、鄉(xiāng)村生活空間與非農(nóng)生產(chǎn)空間的轉(zhuǎn)換則主要受城鄉(xiāng)建設(shè)活動的驅(qū)動。其中,城鎮(zhèn)生活空間與鄉(xiāng)村生活空間均以農(nóng)業(yè)生產(chǎn)空間的轉(zhuǎn)入為主,空間轉(zhuǎn)換主要分布于城鎮(zhèn)與鄉(xiāng)村生活空間的外圍,城鎮(zhèn)生活空間的轉(zhuǎn)入規(guī)模遠小于鄉(xiāng)村生活空間。非農(nóng)生產(chǎn)空間也以轉(zhuǎn)入為主,且轉(zhuǎn)入規(guī)模逐階段擴張,轉(zhuǎn)入類型以農(nóng)業(yè)種植空間、林地空間為主,主要分布于城鎮(zhèn)生活空間外圍,其中位于水源區(qū)東部與西部地區(qū)的空間轉(zhuǎn)換較為集中。裸地空間的轉(zhuǎn)換規(guī)模在所有空間中最小,轉(zhuǎn)入面積均大于轉(zhuǎn)出面積,空間轉(zhuǎn)換零星散布于水源區(qū)內(nèi)。

    3. 2 生態(tài)系統(tǒng)服務(wù)演化分析

    3. 2. 1 不同類型“三生”空間生態(tài)系統(tǒng)服務(wù)變化

    根據(jù)表7、表8,2000—2020年南水北調(diào)中線工程水源區(qū)內(nèi)各類“三生”空間的生態(tài)系統(tǒng)服務(wù)總量均存在不同程度變化,但各類空間的單位面積生態(tài)系統(tǒng)服務(wù)量與生態(tài)系統(tǒng)服務(wù)總量的比重相對穩(wěn)定。2000—2020年,林地空間、草地空間與農(nóng)業(yè)生產(chǎn)空間的平均水源涵養(yǎng)深度、水源涵養(yǎng)總量以及平均土壤保持量與土壤保持總量均明顯高于其他空間,分別在210 mm、5×109 t 與3. 9×104 t/km2、1×109 t以上。不同時期內(nèi),上述3類空間的水源涵養(yǎng)總量占比與土壤保持總量占比也均在99%以上。除農(nóng)業(yè)生產(chǎn)空間、草地空間外,其余空間的水源涵養(yǎng)總量與土壤保持總量均有增長。其中,非農(nóng)生產(chǎn)空間、水域空間、城鎮(zhèn)生活空間與裸地空間的增幅較大,水源涵養(yǎng)與土壤保持總量的增幅分別在60%與40%以上,而林地空間的增幅最小,但其水源涵養(yǎng)與土壤保持服務(wù)的增長總量均遠高于其他空間,分別為5. 81×108 t與1. 37×108 t。

    3. 2. 2 生態(tài)系統(tǒng)服務(wù)空間分異特征與變化

    根據(jù)圖3,2000—2020年,南水北調(diào)中線工程水源區(qū)內(nèi)水源涵養(yǎng)與水土保持服務(wù)的空間格局基本穩(wěn)定,并具有不同的空間分異特征。水源涵養(yǎng)服務(wù)主要呈“南高北低”的分布特征。其中,高值區(qū)主要分布于漢江以北的大巴山、武當山區(qū),水源涵養(yǎng)深度大多在400~500 mm以及500 mm以上;中值區(qū)主要分布于大巴山、武當山區(qū)的低海拔地區(qū)以及伏牛山的高海拔區(qū)域,水源涵養(yǎng)深度約為300~400 mm;低值區(qū)則主要位于水源區(qū)中部及北部的低海拔區(qū)域,水源涵養(yǎng)深度大多在200~300 mm以及200 mm以下。水土保持服務(wù)則呈“高值分散,低值連片”的分布特征。其中,高、中值區(qū)主要分布于坡度較陡的區(qū)域,單位面積的土壤保持量大多在1. 5×105 t/km2~2×105 t/km2及2×105 t/km2以上。低值區(qū)分布廣泛且集中連片,主要分布于坡度平緩的區(qū)域,單位面積的土壤保持總量均在5×104t/km2以下。在變化趨勢上,水源涵養(yǎng)與水土保持服務(wù)也保持相對穩(wěn)定。根據(jù)表7、表9,2000—2020年平均水源涵養(yǎng)深度與水源涵養(yǎng)總量僅降低0. 27%,單位面積土壤保持量與土壤保持總量則均增長0. 04%。

    4 “三生”空間演化的生態(tài)效應(yīng)與驅(qū)動機制

    4. 1 “三生”空間演化的生態(tài)效應(yīng)

    根據(jù)圖4、表10,在2000—2010年與2010—2020年2個時期,南水北調(diào)中線工程水源區(qū)內(nèi)“三生”空間演化對其生態(tài)系統(tǒng)服務(wù)產(chǎn)生顯著影響,可分為正向與負向效應(yīng)。2000—2020年,水源涵養(yǎng)服務(wù)的增長與減少區(qū)域在水源區(qū)分布廣泛。不同時期內(nèi),農(nóng)業(yè)生產(chǎn)空間轉(zhuǎn)為林地空間、草地空間以及草地空間轉(zhuǎn)為林地空間對水源涵養(yǎng)服務(wù)增長的貢獻率較高,3種空間轉(zhuǎn)換類型的貢獻率總和分別為86. 32%與86. 95%。2010—2020年,由于漢中、安康與商洛地區(qū)草地空間大規(guī)模轉(zhuǎn)為林地空間,使該區(qū)域水源涵養(yǎng)服務(wù)的增長區(qū)域相對集中并具有較大規(guī)模。而農(nóng)業(yè)生產(chǎn)空間、林地空間轉(zhuǎn)為水域空間、林地空間、草地空間轉(zhuǎn)為農(nóng)業(yè)生產(chǎn)空間、農(nóng)業(yè)生產(chǎn)空間轉(zhuǎn)為鄉(xiāng)村生活空間、林地空間轉(zhuǎn)變?yōu)榉寝r(nóng)生產(chǎn)空間對水源區(qū)水源涵養(yǎng)服務(wù)減少的貢獻率較高。其中,農(nóng)業(yè)生產(chǎn)空間轉(zhuǎn)為鄉(xiāng)村生活空間的貢獻主要體現(xiàn)在2000—2010年,主要集中在南陽、漢中與商洛地區(qū);林地空間轉(zhuǎn)變?yōu)榉寝r(nóng)生產(chǎn)空間的貢獻則主要體現(xiàn)在2010—2020 年,主要集中在十堰市。在水土保持方面,2000—2020年,水土保持服務(wù)的增長與減少區(qū)域也廣泛分布于水源區(qū)內(nèi)。不同時期內(nèi),對水土保持服務(wù)增長貢獻較高的空間轉(zhuǎn)換類型與水源涵養(yǎng)服務(wù)相同,均為農(nóng)業(yè)生產(chǎn)空間轉(zhuǎn)為草地空間、林地空間以及草地空間轉(zhuǎn)為林地空間,3 類空間轉(zhuǎn)換類型的貢獻率總和分別為87. 67% 與89. 56%。而林地空間、草地空間轉(zhuǎn)為農(nóng)業(yè)生產(chǎn)空間以及林地空間、草地空間、農(nóng)業(yè)生產(chǎn)空間轉(zhuǎn)為裸地空間對水土保持服務(wù)減少的貢獻率較高,特別是林地空間、草地空間轉(zhuǎn)為農(nóng)業(yè)生產(chǎn)空間,兩種空間轉(zhuǎn)換類型的貢獻率總和分別為74. 58%與82. 34%。

    基于南水北調(diào)中線工程水源區(qū)“三生”空間演化與生態(tài)效應(yīng)分析發(fā)現(xiàn),水源區(qū)農(nóng)業(yè)生產(chǎn)空間、林地空間與草地空間的空間沖突日趨激烈,使當?shù)厣鷳B(tài)系統(tǒng)服務(wù)產(chǎn)生劇烈的空間波動。根據(jù)表5、表6、表10,在2000—2010 年與2010—2020年,分別有848. 77 km2與1 987. 8 km2的農(nóng)業(yè)生產(chǎn)空間轉(zhuǎn)換為林地空間與草地空間,同時也有991. 15km2與1 553. 21 km2的林地空間、草地空間轉(zhuǎn)換為農(nóng)業(yè)生產(chǎn)空間,而上述空間轉(zhuǎn)換在水源區(qū)水源涵養(yǎng)與水土保持服務(wù)總量增減變化中具有較高的貢獻率,其中農(nóng)業(yè)生產(chǎn)空間轉(zhuǎn)為林地空間、草地空間的貢獻率總和分別為68. 14%、62. 09%與71. 25%、61. 48%,林地空間、草地空間轉(zhuǎn)換為農(nóng)業(yè)生產(chǎn)空間的貢獻總和則分別為37. 58%、48. 13% 與74. 58%、82. 34%。因此,建議通過全域土地綜合整治,統(tǒng)籌南水北調(diào)中線工程水源區(qū)內(nèi)農(nóng)業(yè)生產(chǎn)與生態(tài)保護的整體格局,在劃定永久基本農(nóng)田與生態(tài)保護紅線的基礎(chǔ)上,結(jié)合農(nóng)戶意愿與生態(tài)保護需求,合理確定其他農(nóng)業(yè)生產(chǎn)空間與各類生態(tài)空間內(nèi)的生態(tài)保護與農(nóng)業(yè)開發(fā)格局,進而減少農(nóng)業(yè)生產(chǎn)與生態(tài)保護活動的空間競爭。

    4. 2 生態(tài)效應(yīng)的驅(qū)動機制分析

    基于南水北調(diào)中線工程水源區(qū)“三生”空間演化及其生態(tài)效應(yīng)分析結(jié)果以及水源區(qū)實地調(diào)研,發(fā)現(xiàn)農(nóng)業(yè)生產(chǎn)、城鄉(xiāng)建設(shè)、水利建設(shè)與生態(tài)保護等人類活動是水源區(qū)“三生”空間生態(tài)效應(yīng)的直接驅(qū)動力,主要通過改變國土空間利用方式推動生態(tài)系統(tǒng)服務(wù)演化發(fā)展(圖5)。其中,農(nóng)業(yè)生產(chǎn)活動、城鄉(xiāng)建設(shè)活動與水利建設(shè)活動是南水北調(diào)中線工程水源區(qū)生態(tài)系統(tǒng)服務(wù)流失的主要原因。農(nóng)業(yè)是水源區(qū)內(nèi)鄉(xiāng)村居民的主要收入來源,為增加務(wù)農(nóng)收入,鄉(xiāng)村居民不斷增加農(nóng)業(yè)生產(chǎn)投入并通過侵占林地空間、草地空間擴大農(nóng)業(yè)生產(chǎn)規(guī)模。根據(jù)前文,林地空間、草地空間轉(zhuǎn)變?yōu)檗r(nóng)業(yè)生產(chǎn)空間對水源涵養(yǎng)與水土保持服務(wù)的減少的貢獻度均較高。地方政府為發(fā)展經(jīng)濟、改善城鄉(xiāng)居民生活環(huán)境不斷引導(dǎo)水源區(qū)內(nèi)工業(yè)、服務(wù)業(yè)發(fā)展并開展城鄉(xiāng)建設(shè)活動。2000—2020年,工業(yè)、服務(wù)業(yè)總產(chǎn)值從349. 89億元增長到5 258. 6億元。非農(nóng)生產(chǎn)空間與城、鄉(xiāng)生活空間作為工業(yè)、服務(wù)業(yè)與城鄉(xiāng)居民生活的載體,各類空間面積也分別增長1 303. 98%與47. 11%,但也進一步導(dǎo)致農(nóng)業(yè)生產(chǎn)空間、林地空間的流失,進而使水源區(qū)生態(tài)系統(tǒng)服務(wù)減少。為配合南水北調(diào)中線工程“自流引水”的建設(shè)目標,丹江口大壩加高加寬續(xù)建工程陸續(xù)開工,水庫蓄水位從2005年的157 m提高到2013年的170 m,使得水面上漲并淹沒周圍農(nóng)業(yè)生產(chǎn)空間、林地空間等空間類型。根據(jù)表10,在2000—2010年與2010—2020年,農(nóng)業(yè)生產(chǎn)空間與林地空間轉(zhuǎn)變?yōu)樗蚩臻g對水源涵養(yǎng)服務(wù)減少的貢獻率總和分別為30. 87%與25. 39%。生態(tài)保護活動則是水源區(qū)生態(tài)系統(tǒng)服務(wù)增長的主要原因。2000年以來,為提升水源區(qū)的水源涵養(yǎng)與水土保持服務(wù),各級政府積極在水源區(qū)內(nèi)開展植樹造林、退耕還林、還草等生態(tài)工程,促進大量農(nóng)業(yè)種植空間向林地空間、草地空間以及草地空間向林地空間轉(zhuǎn)換。根據(jù)表10,在2000—2010年與2010—2020年,上述空間轉(zhuǎn)換對水源區(qū)水源涵養(yǎng)與水土保持服務(wù)增長的貢獻率總和分別為86. 32%、86. 95%與87. 67%、89. 56%。

    5 結(jié)論與討論

    5. 1 結(jié)論

    本研究基于多源數(shù)據(jù)與多學(xué)科方法,在識別南水北調(diào)中線工程水源區(qū)“三生”空間并評估水源涵養(yǎng)與水土保持兩類生態(tài)系統(tǒng)服務(wù)的基礎(chǔ)上,分析水源區(qū)“三生”空間與生態(tài)系統(tǒng)服務(wù)演化特征,并基于生態(tài)系統(tǒng)服務(wù)視角解析水源區(qū)“三生”空間演化的生態(tài)效應(yīng)及其驅(qū)動力,得出如下結(jié)論。

    (1)南水北調(diào)中線工程水源區(qū)內(nèi)各類“三生”空間的空間格局基本穩(wěn)定,除農(nóng)業(yè)生產(chǎn)空間、草地空間外,其他空間的面積均有增長,其中林地空間的增量最多,非農(nóng)生產(chǎn)空間增幅最大。不同時期內(nèi),農(nóng)業(yè)生產(chǎn)空間、林地空間與草地空間的轉(zhuǎn)換面積均遠高于其他空間。

    (2)水源區(qū)整體及各類空間內(nèi)的生態(tài)系統(tǒng)服務(wù)相對穩(wěn)定,林地空間、草地空間與農(nóng)業(yè)生產(chǎn)空間的生態(tài)系統(tǒng)服務(wù)水平顯著高于其他空間。“三生”空間轉(zhuǎn)換對水源區(qū)生態(tài)系統(tǒng)服務(wù)存在明顯的正向與負向效應(yīng),其中,農(nóng)業(yè)生產(chǎn)空間轉(zhuǎn)為林地空間、草地空間以及草地空間轉(zhuǎn)為林地空間對水源涵養(yǎng)與水土保持服務(wù)增長的貢獻率較高。

    (3)南水北調(diào)中線工程水源區(qū)“三生”空間演化的生態(tài)效應(yīng)與人類活動密切相關(guān),農(nóng)業(yè)生產(chǎn)、城鄉(xiāng)建設(shè)、水利建設(shè)以及生態(tài)保護等人類活動通過改變空間利用方式推動生態(tài)系統(tǒng)服務(wù)演化發(fā)展。其中,農(nóng)業(yè)生產(chǎn)、城鄉(xiāng)建設(shè)與水利建設(shè)活動對水源區(qū)生態(tài)系統(tǒng)服務(wù)具有負向效應(yīng),生態(tài)保護活動則發(fā)揮正向效應(yīng)。

    (4)南水北調(diào)中線工程水源區(qū)內(nèi)農(nóng)業(yè)生產(chǎn)空間、林地空間與草地空間的空間競爭日趨激烈,未來,政府應(yīng)進一步推進南水北調(diào)中線工程水源區(qū)全域土地綜合整治,綜合統(tǒng)籌水源區(qū)內(nèi)農(nóng)業(yè)生產(chǎn)與生態(tài)保護的整體格局。

    5. 2 討論

    本研究立足生態(tài)系統(tǒng)服務(wù)視角,在開展南水北調(diào)中線工程水源區(qū)“三生”空間識別與生態(tài)系統(tǒng)服務(wù)評估的基礎(chǔ)上,系統(tǒng)分析水源區(qū)“三生”空間演化特征及其生態(tài)環(huán)境效應(yīng),豐富“三生”空間演化與生態(tài)效應(yīng)研究的內(nèi)容與方法體系?;贗nVEST模型開展“三生”空間生態(tài)系統(tǒng)服務(wù)評估,可以定量刻畫不同“三生”空間及各類空間內(nèi)部生態(tài)質(zhì)量的高低,有利于從數(shù)量和空間兩個方面解析“三生”空間演化對生態(tài)環(huán)境的影響。相比于傳統(tǒng)的生態(tài)環(huán)境質(zhì)量指數(shù)法,本方法基于生態(tài)學(xué)原理評估“三生”空間的生態(tài)系統(tǒng)服務(wù)能力,具有良好的科學(xué)性與客觀性,可以準確反映“三生”空間內(nèi)部生態(tài)質(zhì)量,彌補了生態(tài)環(huán)境質(zhì)量指數(shù)法的不足,從而為解釋國土空間演化的生態(tài)效應(yīng)與制定國土空間規(guī)劃策略提供決策支持。但受方法限制,本研究對水域空間水源涵養(yǎng)功能的測算相對片面。今后可基于生態(tài)學(xué)、水文學(xué)等多學(xué)科知識完善水源涵養(yǎng)功能評價方法,并對“三生”空間演化的生態(tài)效應(yīng)進行深入分析。

    參考文獻

    [1] 劉繼來,劉彥隨,李裕瑞. 中國“三生空間” 分類評價與時空格局

    分析[J]. 地理學(xué)報,2017,72(7):1290-1304.

    [2] 董寅,金貴,鄧祥征. 中國國土空間布局優(yōu)化研究[J]. 地理學(xué)

    報,2024,79(3):672-687.

    [3] 孫傳諄,李鵬,鄧羽,等. 基于生態(tài)系統(tǒng)空間異質(zhì)性的美麗中國

    生態(tài)建設(shè)分區(qū)[J]. 地理學(xué)報,2022,77(11):2902-2919.

    [4] QU Y B,ZHAO W Y,ZHAO L J,et al. Research on hollow village

    governance based on action network:mode,mechanism and countermeasures——

    comparison of different patterns in plain agricultural

    areas of China[J]. Land,2022,11(6):792.

    [5] 李廣東,方創(chuàng)琳. 城市生態(tài)-生產(chǎn)-生活空間功能定量識別與分

    析[J]. 地理學(xué)報,2016,71(1):49-65.

    [6] 徐凱,房艷剛,周宇航. 基于“三生” 視角的鎮(zhèn)域鄉(xiāng)村功能空間分

    化研究:以山東省為例[J]. 地理科學(xué)進展,2023,42(4):

    644-656.

    [7] 劉超,張祚,紀旭,等. 基于土地功能強弱對比的三生空間識別

    與優(yōu)化:以襄陽市為例[J]. 資源科學(xué),2023,45(7):1366-1379.

    [8] JIN J F,YIN S Y,YIN H M,et al. Eco?environmental effects of‘ production-

    living-ecological’ space land use changes and recommendations

    for ecological restoration:a case study of the Weibei dryland

    in Shaanxi Province[J]. Land,2023,12(5):1060.

    [9] 黃安,許月卿,盧龍輝,等.“生產(chǎn)-生活-生態(tài)” 空間識別與優(yōu)化

    研究進展[J]. 地理科學(xué)進展,2020,39(3):503-518.

    [10] LIN G,JIANG D,F(xiàn)U J Y,et al.A review on the overall optimization

    of production-living-ecological space:theoretical basis and

    conceptual framework[J].Land,2022,11(3):345.

    [11] 武子豪,祖健,史云揚,等. 城市職能視角下“三生” 空間的識別

    與評價:以京津冀城市群為例[J]. 資源科學(xué),2022,44(11):

    2247-2259.

    [12] JI Z X,LIU C,XU Y Q,et al.Quantitative identification and the

    evolution characteristics of production-living-ecological space in

    the mountainous area:from the perspective of multifunctional land

    [J].Journal of geographical sciences,2023,33(4):779-800.

    [13] 陳鴻基,楊慶媛,彭立嫻,等. 三峽庫區(qū)縣域“三生空間”時空演

    變特征與情景模擬[J]. 農(nóng)業(yè)工程學(xué)報,2022,38(13):285-294.

    [14] LI H Z,F(xiàn)ANG C Y,XIA Y,et al.Multi?scenario simulation of production-

    living-ecological space in the Poyang Lake area based on

    remote sensing and RF?Markov?FLUS model[J].Remote sensing,

    2022,14(12):2830.

    [15] 潘方杰,萬慶,曾菊新,等. 城鎮(zhèn)化進程中湖北省“生產(chǎn)-生活-

    生態(tài)” 空間沖突演化及其影響因素[J]. 經(jīng)濟地理,2023,43

    (2):80-92.

    [16] 李紅潤,劉慧芳,王瑾,等. 基于Markov-FLUS-MCR模型的晉

    中市“ 三生” 空間優(yōu)化[J]. 農(nóng)業(yè)工程學(xué)報,2022,38(10):

    267-276.

    [17] 韓美,孔祥倫,李云龍,等. 黃河三角洲“三生”用地轉(zhuǎn)型的生態(tài)

    環(huán)境效應(yīng)及其空間分異機制[J]. 地理科學(xué),2021,41(6):1009-

    1018.

    [18] 姜昕彤,翟石艷,王錚,等. 基于未來土地利用模擬模型的鄭州

    市“三生空間”模擬及生態(tài)環(huán)境效應(yīng)分析[J]. 生態(tài)學(xué)報,2023,

    43(15):6225-6242.

    [19] 馮雨柯,翟石艷,姜昕彤,等. 中原城市群“三生”空間模擬及生

    態(tài)效應(yīng)評估:基于未來土地利用模擬模型和共享社會經(jīng)濟路徑

    情景[J]. 生態(tài)學(xué)報,2023,43(20):8292-8308.

    [20] 孔冬艷,陳會廣,吳孔森. 中國“ 三生空間”演變特征、生態(tài)環(huán)境

    效應(yīng)及其影響因素[J]. 自然資源學(xué)報,2021,36(5):1116-

    1135.

    [21] YANG Y Y,BAO W K,LI Y H,et al. Land use transition and its

    eco?environmental effects in the Beijing?Tianjin?Hebei urban agglomeration:

    a production-living-ecological perspective[J]. Land,

    2020,9(9):285.

    [22] 黃晶,薛東前,董朝陽,等. 干旱綠洲農(nóng)業(yè)區(qū)土地利用轉(zhuǎn)型生態(tài)

    環(huán)境效應(yīng)及分異機制:基于三生空間主導(dǎo)功能判別視角[J]. 地

    理科學(xué)進展,2022,41(11):2044-2060.

    [23] 董冬,羅毅,顧康康.“三生空間”視角下長三角城市群土地利用

    轉(zhuǎn)型生態(tài)環(huán)境效應(yīng)時空分異及驅(qū)動力[J]. 長江流域資源與環(huán)

    境,2023,32(8):1664-1676.

    [24] COSTANZA R,D'ARGE R,DE GROOT R,et al. The value of the

    world's ecosystem services and natural capital[J]. Ecological economics,

    1998,25(1):3-15.

    [25] DE GROOT R.Function?analysis and valuation as a tool to assess

    land use conflicts in planning for sustainable,multi?functional landscapes

    [J].Landscape and urban planning,2006,75(3/4):175-186.

    [26] RYALLS J M W,LANGFORD B,MULLINGER N J,et al. Anthropogenic

    air pollutants reduce insect?mediated pollination services

    [J]. Environmental pollution,2022,297:118847.

    [27] WANG A Y,LIAO X Y,TONG Z J,et al.Spatial?temporal dynamic

    evaluation of the ecosystem service value from the perspective of‘ production-

    living-ecological’ spaces:a case study in Dongliao River Basin,

    China[J].Journal of cleaner production,2022,333:130218.

    [28] ZHAO F F,LIU X X,ZHAO X,et al.Effects of production-livingecological

    space changes on the ecosystem service value of the

    Yangtze River Delta urban agglomeration in China[J].Environmental

    monitoring and assessment,2023,195(9):1133.

    [29] PAN F J,SHU N N,WAN Q,et al. Land use function transition and

    associated ecosystem service value effects based on production-living-

    ecological space:a case study in the Three Gorges Reservoir area

    [J]. Land,2023, 12(2):391.

    [30] WANG L T,WANG S X,ZHOU Y,et al.Landscape pattern variation,

    protection measures,and land use/land cover changes in drinking water

    source protection areas:a case study in Danjiangkou Reservoir,

    China[J]. Global ecology and conservation,2020,21:e00827.

    [31] 汪倩,陳軍飛. 南水北調(diào)工程通水對受水區(qū)水資源集約利用效

    率的影響:基于河南省市級層面的實證[J]. 中國人口·資源與

    環(huán)境,2022,32(6):155-164.

    [32] 吳夢. 南水北調(diào)中線工程受水區(qū)灰水足跡荷載系數(shù)時空變化

    及驅(qū)動因素研究[J]. 長江流域資源與環(huán)境,2024,33(5):

    971-981.

    [33] DUAN Y M,WANG H,HUANG A,et al.Identification and spatialtemporal

    evolution of rural‘ production-living-ecological’ space from

    the perspective of villagers’ behavior:A case study of Ertai Town,

    Zhangjiakou City[J].Land use policy,2021,106:105457.

    [34] FU J Y,BU Z Q,JIANG D,et al. Sustainable land use diagnosis

    based on the perspective of production-living-ecological spaces in

    China[J]. Land use policy,2022,122:106386.

    [35] 凌子燕,李延順,蔣衛(wèi)國,等. 山江海交錯帶城市群國土三生空

    間動態(tài)變化特征:以廣西北部灣城市群為例[J]. 經(jīng)濟地理,

    2022,42(2):18-24.

    [36] 段亞明,黃安,盧龍輝,等.“生產(chǎn)-生活-生態(tài)” 空間的概念與理

    論研究[J]. 中國農(nóng)業(yè)大學(xué)學(xué)報,2023,28(4):170-182.

    [37] 許月卿,盧龍輝,冀正欣,等. 國土空間精細化管理視角下“三

    生”空間細分類型的演化特征與驅(qū)動機制研究:以河北省張家

    口市為例[J]. 北京大學(xué)學(xué)報(自然科學(xué)版),2024,60(4):

    732-744.

    [38] SONG Y Y,XIA S Y,XUE D Q,et al. Land space change process

    and its eco?environmental effects in the Guanzhong Plain urban agglomeration

    of China[J]. Land,2022,11(9):1547.

    [39] 楊燕,胡書玲,王坤,等. 少數(shù)民族旅游村寨的三生空間演變及

    驅(qū)動機制:以肇興侗寨為例[J]. 經(jīng)濟地理,2023,43(12):

    222-230.

    [40] 付晶瑩,郜強,江東,等. 黑土保護與糧食安全背景下齊齊哈爾

    市國土空間優(yōu)化調(diào)控路徑[J]. 地理學(xué)報,2022,77(7):1662-

    1680.

    [41] CHEN Z P,LIU Y Q,TU S S.Comprehensive eco?environmental effects

    caused by land use transition from the perspective of productionliving-

    ecological spaces in a typical region:a case study of the Guangxi

    Zhuang autonomous region,China[J].Land,2022,11(12):2160.

    [42] 袁雪松,周俊,胡蓓蓓,等. 基于FLUS模型粵港澳大灣區(qū)“三生

    空間”多情景模擬預(yù)測[J]. 地理科學(xué),2023,43(3):564-574.

    [43] LOVELL S T,DESANTIS S,NATHAN C A,et al. Integrating agroecology

    and landscape multifunctionality in Vermont:an evolving

    framework to evaluate the design of agroecosystems[J]. Agricultural

    systems,2010, 103(5):327-341.

    [44] YAO J,LIU M,CHEN N N,et al. The spatial variation of air purification

    benefit provided by street tree assemblages in Shenyang,

    China[J]. Urban ecosystems,2022,25(3):725-732.

    [45] JIANG W. Ecosystem services research in China:a critical review

    [J]. Ecosystem services,2017,26:10-16.

    [46] LIU R,NIU X,WANG B,et al. InVEST model?based spatiotemporal

    analysis of water supply services in the Zhangcheng district[J].

    Forests,2021,12(8):1082.

    [47] 劉美娟,仲俊濤,王蓓,等. 基于InVEST模型的青海湖流域產(chǎn)

    水功能時空變化及驅(qū)動因素分析[J]. 地理科學(xué),2023,43(3):

    411-422.

    [48] YANG X,CHEN R S,MEADOWS M E,et al. Modelling water

    yield with the InVEST model in a data scarce region of Northwest

    China[J]. Water supply,2020,20(3):1035-1045.

    [49] LIU L H,ZHENG L,WANG Y,et al.Land use and ecosystem services

    evolution in Danjiangkou Reservoir area,China:implications

    for sustainable management of national projects[J].Land,2023,12

    (4):788.

    [50] 胡硯霞,于興修,廖雯,等. 漢江流域產(chǎn)水量時空格局及影響因

    素研究[J]. 長江流域資源與環(huán)境,2022,31(1):73-82.

    [51] YIN C C,ZHAO W W,PEREIRA P.Soil conservation service underpins

    sustainable development goals[J]. Global ecology and conservation,

    2022,33:e01974.

    [52] WISCHMEIER W H, JOHNSON C B,CROSS B V. A soil erodibility

    monograph for farmland and construction sites[J]. Journal of

    soil and water conservation,1971,26:189-193.

    [53] 廖雯,胡硯霞,于興修,等. 丹江口庫區(qū)2010—2020年土壤保持

    功能時空特征及其影響因素[J]. 水土保持通報,2021,41(6):

    288-294.

    [54] 李亦秋,馮仲科,韓烈保,等. 丹江口庫區(qū)及上游生態(tài)系統(tǒng)土壤

    保持效益價值評估[J]. 中國人口·資源與環(huán)境,2010,20(5):

    64-69.

    [55] LI K,ZHANG B Y,XIAO W D,et al. Land use transformation

    based on production-living-ecological space and associated ecoenvironment

    effects:a case study in the Yangtze River Delta urban

    agglomeration[J].Land, 2022,11(7):1076.

    [56] LI C X,WU J Y. Land use transformation and eco?environmental

    effects based on production-living-ecological spatial synergy:evidence

    from Shaanxi Province,China[J]. Environmental science

    and pollution research international,2022,29(27):41492-41504.

    [57] 董建紅,張志斌,笪曉軍,等.“三生” 空間視角下土地利用轉(zhuǎn)型

    的生態(tài)環(huán)境效應(yīng)及驅(qū)動力:以甘肅省為例[J]. 生態(tài)學(xué)報,2021,

    41(15):5919-5928.

    (責(zé)任編輯:劉呈慶)

    亚洲成人免费电影在线观看| 成人精品一区二区免费| 国产69精品久久久久777片| 国产淫片久久久久久久久 | 男人舔奶头视频| 大型黄色视频在线免费观看| 免费看a级黄色片| 久久香蕉精品热| 好看av亚洲va欧美ⅴa在| 高清在线国产一区| 亚洲男人的天堂狠狠| 国产成人a区在线观看| av福利片在线观看| 又黄又爽又免费观看的视频| 小说图片视频综合网站| 一区二区三区激情视频| 午夜免费成人在线视频| 动漫黄色视频在线观看| 人妻夜夜爽99麻豆av| 全区人妻精品视频| 国产精品亚洲一级av第二区| 久久国产精品影院| 在线观看免费视频日本深夜| 婷婷精品国产亚洲av| av福利片在线观看| 夜夜看夜夜爽夜夜摸| 91久久精品电影网| 国产三级在线视频| 丁香欧美五月| 国产精品香港三级国产av潘金莲| av欧美777| 91av网一区二区| avwww免费| 少妇丰满av| 久久婷婷人人爽人人干人人爱| 久久久久久久久大av| 99在线人妻在线中文字幕| 免费大片18禁| 国内毛片毛片毛片毛片毛片| а√天堂www在线а√下载| 一个人观看的视频www高清免费观看| 国产亚洲av嫩草精品影院| 男女视频在线观看网站免费| 男女视频在线观看网站免费| 老司机午夜十八禁免费视频| 久久久久性生活片| 久久伊人香网站| 99热这里只有精品一区| 国产又黄又爽又无遮挡在线| 中文字幕久久专区| 成人精品一区二区免费| 成人永久免费在线观看视频| 国产精品99久久久久久久久| 婷婷亚洲欧美| 亚洲专区国产一区二区| 国产99白浆流出| 国产乱人视频| 男人舔奶头视频| 国产野战对白在线观看| 日本与韩国留学比较| 无遮挡黄片免费观看| 草草在线视频免费看| 日韩欧美三级三区| 亚洲精品粉嫩美女一区| 中文字幕熟女人妻在线| 精品久久久久久久末码| 真人一进一出gif抽搐免费| 国内精品久久久久久久电影| 国产 一区 欧美 日韩| 亚洲av成人精品一区久久| 91在线观看av| 亚洲va日本ⅴa欧美va伊人久久| 在线观看免费视频日本深夜| 男女午夜视频在线观看| 久久精品国产清高在天天线| 淫秽高清视频在线观看| 日本 欧美在线| 麻豆成人午夜福利视频| 亚洲精品456在线播放app | 草草在线视频免费看| 国产一区二区在线av高清观看| 一个人免费在线观看电影| www.www免费av| 少妇人妻精品综合一区二区 | 国产精品久久久久久精品电影| 91麻豆精品激情在线观看国产| a级毛片a级免费在线| 亚洲性夜色夜夜综合| 成人午夜高清在线视频| 国产欧美日韩一区二区精品| 精品99又大又爽又粗少妇毛片 | 成人av在线播放网站| 欧洲精品卡2卡3卡4卡5卡区| 亚洲天堂国产精品一区在线| 啦啦啦观看免费观看视频高清| 欧美日韩精品网址| 又紧又爽又黄一区二区| 久久精品国产亚洲av香蕉五月| 嫩草影院精品99| 俺也久久电影网| 国产成人aa在线观看| 尤物成人国产欧美一区二区三区| 一本综合久久免费| 黄色日韩在线| 亚洲一区高清亚洲精品| 一区二区三区国产精品乱码| 一进一出抽搐gif免费好疼| 2021天堂中文幕一二区在线观| 亚洲av免费在线观看| 99久久无色码亚洲精品果冻| 在线国产一区二区在线| 亚洲 欧美 日韩 在线 免费| 亚洲成人中文字幕在线播放| 亚洲国产欧美人成| 女人十人毛片免费观看3o分钟| 精品欧美国产一区二区三| 国产美女午夜福利| 欧美3d第一页| 亚洲人成网站高清观看| 成人永久免费在线观看视频| 国产乱人视频| 他把我摸到了高潮在线观看| 欧美性感艳星| 真人一进一出gif抽搐免费| 精品99又大又爽又粗少妇毛片 | 国产精品日韩av在线免费观看| 久久精品91无色码中文字幕| 成年免费大片在线观看| 成人国产综合亚洲| 亚洲片人在线观看| 久久九九热精品免费| 精品国内亚洲2022精品成人| 97人妻精品一区二区三区麻豆| 亚洲国产精品成人综合色| av黄色大香蕉| 国产av在哪里看| 亚洲国产精品成人综合色| 久久人妻av系列| 久久久久久人人人人人| 久久6这里有精品| 制服人妻中文乱码| 国产午夜精品久久久久久一区二区三区 | 亚洲18禁久久av| 美女 人体艺术 gogo| 欧美日本视频| 老汉色∧v一级毛片| 波多野结衣高清作品| 欧美色视频一区免费| 可以在线观看毛片的网站| 又粗又爽又猛毛片免费看| 欧美性猛交╳xxx乱大交人| 免费在线观看成人毛片| 久久人妻av系列| 中文亚洲av片在线观看爽| 九九久久精品国产亚洲av麻豆| 国模一区二区三区四区视频| 国产又黄又爽又无遮挡在线| 动漫黄色视频在线观看| 亚洲精品粉嫩美女一区| 成人av在线播放网站| 亚洲中文日韩欧美视频| 欧美+亚洲+日韩+国产| 哪里可以看免费的av片| 色精品久久人妻99蜜桃| 久久草成人影院| 国产精品99久久久久久久久| 日本免费一区二区三区高清不卡| 午夜福利欧美成人| 九九热线精品视视频播放| 午夜免费激情av| 欧美极品一区二区三区四区| 一本综合久久免费| 男人舔奶头视频| 国产三级黄色录像| 麻豆成人午夜福利视频| 一本精品99久久精品77| 又黄又爽又免费观看的视频| 人人妻,人人澡人人爽秒播| 国产三级黄色录像| 欧美国产日韩亚洲一区| 99国产极品粉嫩在线观看| 天天躁日日操中文字幕| 午夜精品一区二区三区免费看| av中文乱码字幕在线| 欧美xxxx黑人xx丫x性爽| 亚洲在线观看片| 国产男靠女视频免费网站| 国产精品久久视频播放| 精品久久久久久久人妻蜜臀av| 五月伊人婷婷丁香| 熟女人妻精品中文字幕| 日本成人三级电影网站| 精品免费久久久久久久清纯| 小蜜桃在线观看免费完整版高清| 日韩欧美三级三区| 精品99又大又爽又粗少妇毛片 | 国产黄片美女视频| 国内毛片毛片毛片毛片毛片| 美女被艹到高潮喷水动态| 国产一区二区在线观看日韩 | 精品不卡国产一区二区三区| 国产伦人伦偷精品视频| 免费观看精品视频网站| 99久久精品热视频| 免费看日本二区| 色哟哟哟哟哟哟| 一夜夜www| 国内精品久久久久精免费| 2021天堂中文幕一二区在线观| 日韩高清综合在线| 国产精品乱码一区二三区的特点| 亚洲av成人av| 欧美日韩亚洲国产一区二区在线观看| 日本免费一区二区三区高清不卡| 日韩成人在线观看一区二区三区| 久久久国产精品麻豆| 欧美一区二区国产精品久久精品| 岛国在线免费视频观看| 很黄的视频免费| 丁香六月欧美| 亚洲中文字幕一区二区三区有码在线看| 国内揄拍国产精品人妻在线| 国产高清有码在线观看视频| 欧美极品一区二区三区四区| 亚洲色图av天堂| 国产精品久久久久久精品电影| 国产成人影院久久av| 国产精品一区二区免费欧美| 国产视频一区二区在线看| 亚洲美女视频黄频| 欧美xxxx黑人xx丫x性爽| 高清在线国产一区| 国产精品98久久久久久宅男小说| 最新中文字幕久久久久| 女人十人毛片免费观看3o分钟| 国产精品亚洲av一区麻豆| 国产伦在线观看视频一区| 亚洲精品一区av在线观看| 成人av一区二区三区在线看| 丰满人妻一区二区三区视频av | 9191精品国产免费久久| 国产91精品成人一区二区三区| 久久午夜亚洲精品久久| 亚洲av中文字字幕乱码综合| 婷婷精品国产亚洲av在线| 免费大片18禁| 国产成人欧美在线观看| 好看av亚洲va欧美ⅴa在| 婷婷精品国产亚洲av在线| a级毛片a级免费在线| 亚洲成a人片在线一区二区| 亚洲一区二区三区不卡视频| 国产真实伦视频高清在线观看 | 全区人妻精品视频| 岛国视频午夜一区免费看| 久久午夜亚洲精品久久| 欧美一级a爱片免费观看看| 好男人电影高清在线观看| 亚洲性夜色夜夜综合| 午夜久久久久精精品| 久久国产乱子伦精品免费另类| 精品一区二区三区视频在线 | 久久久国产成人免费| 国产精品99久久99久久久不卡| or卡值多少钱| 一级毛片高清免费大全| 19禁男女啪啪无遮挡网站| 中文字幕人妻熟人妻熟丝袜美 | 久久久久久久亚洲中文字幕 | 欧美黑人巨大hd| 国产熟女xx| 亚洲国产欧洲综合997久久,| 国产aⅴ精品一区二区三区波| 国产真实伦视频高清在线观看 | 日本一二三区视频观看| 精品不卡国产一区二区三区| 丝袜美腿在线中文| 少妇熟女aⅴ在线视频| 国产精品日韩av在线免费观看| 国产综合懂色| 国产伦精品一区二区三区视频9 | 嫩草影视91久久| 亚洲成av人片免费观看| 亚洲成人免费电影在线观看| 午夜老司机福利剧场| 91在线观看av| 一区福利在线观看| 最近在线观看免费完整版| 熟妇人妻久久中文字幕3abv| 久久久久久久午夜电影| 成人一区二区视频在线观看| 十八禁网站免费在线| 级片在线观看| 久久久久久久久大av| 欧美另类亚洲清纯唯美| 精华霜和精华液先用哪个| 久久久久性生活片| 午夜福利在线在线| 亚洲欧美激情综合另类| 久久精品人妻少妇| 女警被强在线播放| 91字幕亚洲| 1024手机看黄色片| 人妻夜夜爽99麻豆av| 午夜老司机福利剧场| 99久久九九国产精品国产免费| 俄罗斯特黄特色一大片| 国产亚洲精品一区二区www| 美女高潮喷水抽搐中文字幕| 日本一本二区三区精品| 丰满人妻一区二区三区视频av | 亚洲精品日韩av片在线观看 | 国产高清视频在线播放一区| 男人和女人高潮做爰伦理| 在线观看免费视频日本深夜| 性色avwww在线观看| 亚洲国产精品合色在线| 久久久久久久精品吃奶| 亚洲国产精品合色在线| 欧美激情在线99| 欧美+日韩+精品| 综合色av麻豆| 精品久久久久久久人妻蜜臀av| 十八禁网站免费在线| 国产 一区 欧美 日韩| 黄色女人牲交| 日韩欧美国产在线观看| 欧美午夜高清在线| 最新美女视频免费是黄的| 久久精品国产亚洲av香蕉五月| 成熟少妇高潮喷水视频| 一a级毛片在线观看| 中出人妻视频一区二区| 国产一区二区在线观看日韩 | 此物有八面人人有两片| 国产成人av教育| ponron亚洲| 久9热在线精品视频| 国产一级毛片七仙女欲春2| 久久久久国内视频| 亚洲国产高清在线一区二区三| 国产精品久久视频播放| 久久香蕉精品热| 国模一区二区三区四区视频| 亚洲中文字幕日韩| 男人的好看免费观看在线视频| 欧美日韩综合久久久久久 | 国产精品久久久人人做人人爽| 99热这里只有精品一区| 国产蜜桃级精品一区二区三区| 亚洲男人的天堂狠狠| 中文字幕人成人乱码亚洲影| 欧美另类亚洲清纯唯美| 久久99热这里只有精品18| 国产精品 欧美亚洲| 3wmmmm亚洲av在线观看| 亚洲久久久久久中文字幕| 国产不卡一卡二| 国产精品久久久久久久久免 | 亚洲av美国av| 日本在线视频免费播放| 神马国产精品三级电影在线观看| 免费人成视频x8x8入口观看| 久久人妻av系列| 成人av一区二区三区在线看| 亚洲av熟女| 国产欧美日韩精品一区二区| 国产精品香港三级国产av潘金莲| 别揉我奶头~嗯~啊~动态视频| 久久精品亚洲精品国产色婷小说| 国产91精品成人一区二区三区| 首页视频小说图片口味搜索| 日韩中文字幕欧美一区二区| 日本免费a在线| 免费在线观看日本一区| 村上凉子中文字幕在线| 青草久久国产| 天天一区二区日本电影三级| 深夜精品福利| 精品午夜福利视频在线观看一区| av在线蜜桃| 一夜夜www| 中文字幕人妻熟人妻熟丝袜美 | 午夜福利在线观看吧| 久久久精品大字幕| 中文字幕人成人乱码亚洲影| 国产麻豆成人av免费视频| 又紧又爽又黄一区二区| 免费大片18禁| 国产探花在线观看一区二区| 美女大奶头视频| 99精品久久久久人妻精品| 黄色日韩在线| 叶爱在线成人免费视频播放| 亚洲在线观看片| 国产精品av视频在线免费观看| 一本一本综合久久| 伊人久久精品亚洲午夜| 欧美一区二区国产精品久久精品| 国产在视频线在精品| 国内揄拍国产精品人妻在线| 国产又黄又爽又无遮挡在线| 国产国拍精品亚洲av在线观看 | 99热只有精品国产| 国产av麻豆久久久久久久| 丰满乱子伦码专区| 两性午夜刺激爽爽歪歪视频在线观看| 老鸭窝网址在线观看| 国产精品1区2区在线观看.| 久久久久国内视频| 12—13女人毛片做爰片一| 无人区码免费观看不卡| 日韩欧美精品v在线| 欧美成人免费av一区二区三区| av欧美777| 午夜久久久久精精品| 日本精品一区二区三区蜜桃| 午夜免费男女啪啪视频观看 | 在线天堂最新版资源| 亚洲最大成人手机在线| 波多野结衣高清无吗| 国产成人欧美在线观看| 国产蜜桃级精品一区二区三区| 九色国产91popny在线| 天堂动漫精品| 18+在线观看网站| 一个人免费在线观看电影| 白带黄色成豆腐渣| 欧美乱妇无乱码| 中文字幕av成人在线电影| 夜夜夜夜夜久久久久| eeuss影院久久| 99久久成人亚洲精品观看| 国产成年人精品一区二区| 丁香六月欧美| 一夜夜www| av福利片在线观看| 成人av在线播放网站| 人人妻,人人澡人人爽秒播| 久久精品国产99精品国产亚洲性色| 午夜福利欧美成人| 日本一本二区三区精品| 国产三级黄色录像| 99精品在免费线老司机午夜| 一进一出好大好爽视频| 久久伊人香网站| 欧美在线黄色| 免费在线观看亚洲国产| 免费在线观看影片大全网站| 久久久久国内视频| 欧美不卡视频在线免费观看| 综合色av麻豆| 国产精品野战在线观看| 岛国视频午夜一区免费看| 欧美日韩一级在线毛片| 久久久久久九九精品二区国产| 亚洲成a人片在线一区二区| 少妇丰满av| 男女那种视频在线观看| 久久久国产精品麻豆| 99国产精品一区二区三区| 最近最新免费中文字幕在线| 一进一出抽搐gif免费好疼| 色视频www国产| 午夜福利欧美成人| 天堂网av新在线| 国产三级黄色录像| 天天一区二区日本电影三级| 午夜福利视频1000在线观看| 日韩精品青青久久久久久| 特大巨黑吊av在线直播| 色吧在线观看| 亚洲成av人片在线播放无| 丝袜美腿在线中文| 国产亚洲精品久久久久久毛片| av女优亚洲男人天堂| 日韩精品青青久久久久久| 色av中文字幕| 国产精品一区二区三区四区久久| 欧美日韩国产亚洲二区| 精品免费久久久久久久清纯| 特级一级黄色大片| 99国产综合亚洲精品| 波多野结衣巨乳人妻| 亚洲精品国产精品久久久不卡| 亚洲av成人精品一区久久| 丰满的人妻完整版| 亚洲久久久久久中文字幕| 97人妻精品一区二区三区麻豆| 成年免费大片在线观看| 国产精品影院久久| av女优亚洲男人天堂| 999久久久精品免费观看国产| 十八禁人妻一区二区| 一个人看的www免费观看视频| 亚洲无线在线观看| 看片在线看免费视频| 国产精品三级大全| a级一级毛片免费在线观看| 男女床上黄色一级片免费看| 日本在线视频免费播放| 国产主播在线观看一区二区| 亚洲avbb在线观看| 国产精品av视频在线免费观看| 久久伊人香网站| 18禁国产床啪视频网站| 国产精品一区二区免费欧美| 亚洲狠狠婷婷综合久久图片| 国产精品99久久99久久久不卡| av女优亚洲男人天堂| 亚洲自拍偷在线| 亚洲天堂国产精品一区在线| 综合色av麻豆| 在线a可以看的网站| 美女 人体艺术 gogo| 一a级毛片在线观看| 看免费av毛片| 久久6这里有精品| 欧美日韩黄片免| 久久精品国产99精品国产亚洲性色| 国产精品免费一区二区三区在线| 亚洲成av人片在线播放无| 国产精品久久视频播放| 久久欧美精品欧美久久欧美| 精品久久久久久成人av| 国产三级黄色录像| 欧美不卡视频在线免费观看| 国产美女午夜福利| 久久婷婷人人爽人人干人人爱| 国产av麻豆久久久久久久| 亚洲aⅴ乱码一区二区在线播放| 悠悠久久av| 国产v大片淫在线免费观看| 国产精品一区二区三区四区久久| 天天一区二区日本电影三级| tocl精华| 欧美bdsm另类| 看片在线看免费视频| 欧美在线一区亚洲| 国产爱豆传媒在线观看| 国产精品久久电影中文字幕| 亚洲精品粉嫩美女一区| h日本视频在线播放| 熟女人妻精品中文字幕| 男女做爰动态图高潮gif福利片| 少妇的逼水好多| 日本 欧美在线| 人人妻,人人澡人人爽秒播| av黄色大香蕉| 90打野战视频偷拍视频| 熟女电影av网| 欧美3d第一页| 少妇丰满av| 在线播放国产精品三级| 精品乱码久久久久久99久播| 精品久久久久久久久久久久久| 国产精品久久久人人做人人爽| 亚洲人与动物交配视频| 免费看光身美女| 欧美成人性av电影在线观看| 国产精品一区二区三区四区免费观看 | 亚洲精品美女久久久久99蜜臀| 看免费av毛片| 亚洲国产欧美人成| 在线观看舔阴道视频| 午夜亚洲福利在线播放| 99国产精品一区二区三区| 色播亚洲综合网| 亚洲专区国产一区二区| 国产亚洲av嫩草精品影院| 人妻夜夜爽99麻豆av| 中文字幕精品亚洲无线码一区| 人人妻,人人澡人人爽秒播| 久久久国产成人精品二区| 精华霜和精华液先用哪个| av黄色大香蕉| 久久久久久久久久黄片| 麻豆一二三区av精品| 在线观看午夜福利视频| 在线视频色国产色| 看免费av毛片| 精品久久久久久久末码| 国产午夜精品久久久久久一区二区三区 | 久久香蕉精品热| 亚洲一区高清亚洲精品| 亚洲熟妇熟女久久| av天堂中文字幕网| av片东京热男人的天堂| 宅男免费午夜| 岛国在线观看网站| 亚洲人成伊人成综合网2020| 深夜精品福利| 欧美日韩精品网址| 老司机午夜十八禁免费视频| 国产精品一区二区三区四区久久| 亚洲色图av天堂| 免费在线观看成人毛片| 天堂动漫精品| 免费av不卡在线播放| 亚洲一区高清亚洲精品| 长腿黑丝高跟| 在线a可以看的网站| 久久精品国产99精品国产亚洲性色| 一二三四社区在线视频社区8| 国产精品98久久久久久宅男小说| 亚洲一区高清亚洲精品| 少妇的逼好多水| 免费av不卡在线播放| 狂野欧美激情性xxxx| 女生性感内裤真人,穿戴方法视频| 国产视频一区二区在线看| 精品午夜福利视频在线观看一区| 成年女人永久免费观看视频| 在线观看66精品国产| 国产单亲对白刺激|