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

    生態(tài)安全格局時空演變分析——以太原城市群為例

    2023-11-28 12:57:26祁強強徐占軍
    中國環(huán)境科學(xué) 2023年11期
    關(guān)鍵詞:夾點源地連通性

    祁強強,徐占軍

    生態(tài)安全格局時空演變分析——以太原城市群為例

    祁強強,徐占軍*

    (山西農(nóng)業(yè)大學(xué)資源環(huán)境學(xué)院,山西 太谷 030801)

    以太原城市群為例,利用形態(tài)學(xué)空間格局分析(MSPA)探究2000年、2010年和2020年生態(tài)源地的變化情況,并利用電路理論識別生態(tài)廊道、生態(tài)夾點和生態(tài)屏障點,從而構(gòu)建區(qū)域生態(tài)安全格局,明晰生態(tài)安全格局的演化趨勢,識別生態(tài)修復(fù)的重點區(qū)域,并制定了針對性的修復(fù)措施.研究表明,隨著時間的推移,研究區(qū)生態(tài)源地面積由2000年3083.38km2增長至2020年3144.17km2,但空間分布不均,呈西部多、中東部少的空間格局.受生態(tài)源地增加的影響,區(qū)域內(nèi)生態(tài)廊道數(shù)量由2000年69條增長至2020年90條,景觀連通性增強.其中,西部生態(tài)廊道數(shù)量多且以短距離為主,而東西部生態(tài)廊道長度長、阻力大.研究期間,生態(tài)夾點面積整體呈增加趨勢,由2000年的47.64km2增加至2010年的97.98km2,2020減少至83.77km2;障礙點面積呈減少趨勢,由2000年的451.78km2減少2020年395.95km2,這類區(qū)域應(yīng)作為生態(tài)保護和修復(fù)的重點區(qū)域.本研究可為區(qū)域生態(tài)修復(fù)和景觀規(guī)劃提供參考.

    MSPA;電路理論;生態(tài)安全格局;生態(tài)修復(fù)

    生態(tài)環(huán)境作為人類生存和發(fā)展的基礎(chǔ)[1],其質(zhì)量的好壞直接影響到人類的生存環(huán)境和社會經(jīng)濟的可持續(xù)發(fā)展[2].在過去幾十年里,全球人口的激增和經(jīng)濟的快速發(fā)展[2],嚴(yán)重影響了全球生態(tài)系統(tǒng),導(dǎo)致了棲息地喪失[3]、生物多樣性減少[4]、土壤侵蝕加劇[5]、碳儲量下降[6]等一系列生態(tài)環(huán)境問題,嚴(yán)重威脅著人類的生存和發(fā)展并直接影響著區(qū)域生態(tài)安全和可持續(xù)發(fā)展.因此,如何實現(xiàn)生態(tài)系統(tǒng)正常運行,并保證經(jīng)濟社會的協(xié)調(diào)和可持續(xù)發(fā)展,是目前亟待解決的一個重要問題.近年來,構(gòu)建生態(tài)安全格局已成保持生態(tài)安全、實現(xiàn)可持續(xù)發(fā)展的重要規(guī)劃方法[7].該方法為緩解區(qū)域生態(tài)壓力、實現(xiàn)生態(tài)保護和經(jīng)濟增長之間的均衡發(fā)展提供了一種可能的空間解決方案[7-9].通過對重點生態(tài)保護區(qū)的識別和保護,可以促進區(qū)域生態(tài)系統(tǒng)的健康發(fā)展.

    生態(tài)安全格局自提出以來,國內(nèi)外學(xué)者已開展了大量的研究.現(xiàn)相關(guān)的研究主要集中于生態(tài)安全格局的構(gòu)建[10-11]、評價與優(yōu)化[12-13]及構(gòu)建框架的優(yōu)化[14-16],研究對象涉及城市、流域、工礦區(qū)、國家公園、高原湖泊等多種類型[17-20],研究尺度也從單一尺度轉(zhuǎn)向尺度整合[21].總體來看,研究視角已從最初的理論研究逐漸向構(gòu)建方法與實證研究轉(zhuǎn)變.從構(gòu)建方法來看,目前的主流構(gòu)建范式為“源地-廊道”研究框架[22].其中源地是區(qū)域內(nèi)生境質(zhì)量較高、結(jié)構(gòu)和功能相對穩(wěn)定、對周邊區(qū)域具有重要輻射功能的斑塊[23-24].源地的識別主要包括兩種方法:一是直接選取生態(tài)系統(tǒng)服務(wù)價值較高的自然保護區(qū)、森林公園、風(fēng)景名勝區(qū)等作為生態(tài)源地;二是采用生態(tài)系統(tǒng)服務(wù)功能評價、環(huán)境敏感性測算、生態(tài)承載力評估等定量評估方法[22].但上述方法多偏重斑塊自身的屬性,而未充分考慮斑塊在景觀空間內(nèi)的整體性和連通性[10].形態(tài)學(xué)空間格局分析(MSPA)基于數(shù)學(xué)形態(tài)學(xué)原理對柵格進行度量和分割[23],在空間形態(tài)上更強調(diào)結(jié)構(gòu)連通性、在原理上更強調(diào)生態(tài)過程和生態(tài)網(wǎng)絡(luò),增強生態(tài)源地選取的科學(xué)性、客觀性[24-25].近年來,學(xué)者們將MSPA應(yīng)用到熱環(huán)境[25-26]、藍(lán)綠基礎(chǔ)設(shè)施[27-28]及生態(tài)安全格局[29-31]等源地的識別.因此本文在借鑒相關(guān)研究的基礎(chǔ)上,采用MSPA進行生態(tài)源地的識別.生態(tài)廊道是保持區(qū)域內(nèi)物質(zhì)流、生態(tài)流、生態(tài)過程及能量連續(xù)、連通的帶狀生態(tài)用地[32].廊道的識別方法主要包括最小累計阻力(MCR)及電路理論.MCR模型根據(jù) KNAAPEN 等建立的費用距離修改而來,綜合考慮距離、源地和景觀界面特征,計算物種從源運動到目的斑塊運動過程中所需要耗費的最小代價[32-33].MCR模型選擇出的生態(tài)廊道是最優(yōu)路徑,但卻忽視了生物的隨機游走性,且無法體現(xiàn)廊道的真實寬度[34].電路理論借鑒了物理學(xué)中電子隨機游走的特性[35],模擬物種個體或基因在景觀中的遷移過程,基于電流強度反映生態(tài)斑塊和廊道的相對重要性,從而預(yù)測物種擴散與遷移路徑、識別移動路徑,這一方法更加符合物種真實運動情況[36].

    城市群作為中國區(qū)域發(fā)展中最具活力和潛力的地區(qū)[37],借助高度連通的交通網(wǎng)絡(luò)和有利的區(qū)域激勵政策,推動產(chǎn)業(yè)和人口集聚,是推動產(chǎn)業(yè)調(diào)整和區(qū)域協(xié)調(diào)發(fā)展的關(guān)鍵引擎[11].構(gòu)建城市群內(nèi)生態(tài)安全格局有利于區(qū)域生態(tài)-生產(chǎn)-生活空間的整合和優(yōu)化配置,可有效減少城市化的負(fù)面影響,實現(xiàn)城市群的高質(zhì)量發(fā)展[11].然而,上述研究大多集中于珠江三角洲[38]、京津冀[39]、長江三角洲[40]等城市化水平較高的城市群,而針對太原城市群開展的研究較少.太原城市群位于山西省中心地帶,經(jīng)濟發(fā)展水平較高、發(fā)展?jié)摿^大,在全省經(jīng)濟社會發(fā)展中具有重要的戰(zhàn)略地位[41].與中部其他城市群相比,太原城市群發(fā)展緩慢,是尚在發(fā)育階段的城市群[41].近年來,太原城市群在經(jīng)濟發(fā)展過程中出現(xiàn)了一系列環(huán)境問題[42],嚴(yán)重制約城市群的高質(zhì)量發(fā)展.生態(tài)安全格局作為識別、優(yōu)化國土空間保護修復(fù)區(qū)域的主要手段,可為國土空間修復(fù)提供參考.而目前對生態(tài)安全格局的研究大多集中于生態(tài)源地和廊道的空間分布[43],很少有研究關(guān)注生態(tài)源地和斑塊的重要性.對促進區(qū)域生態(tài)進程、保障區(qū)域生態(tài)安全至關(guān)重要的斑塊、廊道和關(guān)鍵區(qū)域的空間格局和數(shù)量范圍的識別較少[44],無法實現(xiàn)對區(qū)域生態(tài)過程的有效調(diào)控,保障區(qū)域生態(tài)系統(tǒng)的可持續(xù)供給.本研究以太原城市群為例,分析20年來土地利用變化,通過MSPA進行生態(tài)源地識別,利用電路理論進行生態(tài)廊道、夾點和障礙點提取,構(gòu)建生態(tài)安全格局.從而實現(xiàn)對太原城市群生態(tài)安全格局的動態(tài)監(jiān)測、明晰生態(tài)安全格局的演化特點和趨勢、制定針對性的修復(fù)措施以防止生態(tài)功能的退化.以期為區(qū)域生態(tài)修復(fù)和景觀規(guī)劃提供參考.

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

    1.1 研究區(qū)概況

    太原城市群是國家重點建設(shè)的14個城市群之一、是以山西省省會太原為中心、以太原盆地城鎮(zhèn)密集區(qū)為主體構(gòu)成的城市群,是山西省經(jīng)濟最發(fā)達(dá)的地區(qū).依據(jù)《中華人民共和國國民經(jīng)濟和社會發(fā)展第十二個五年規(guī)劃》將太原市10個縣(市、區(qū));晉中市榆次區(qū)、太谷區(qū)、祁縣、平遙縣、介休市;呂梁市交城縣、文水縣、汾陽市、孝義市、嵐縣;忻州市靜樂縣納入太原城市群.研究區(qū)位于黃土高原生態(tài)脆弱區(qū)東部、東南部為太原盆地、西側(cè)為呂梁山區(qū)、東側(cè)為太行山區(qū),呈現(xiàn)“中間低、東西高”的地形特征,以溫帶大陸性氣候為主.研究區(qū)總面積為2.07×104km2.

    圖1 研究區(qū)位置

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

    土地利用數(shù)據(jù)來源于CLCD(China Land Cover Dataset),該數(shù)據(jù)集是由武漢大學(xué)楊杰和黃昕基于GEE(Google Earth Engine)所制作的中國年度土地覆蓋數(shù)據(jù)集(annual China Land Cover Dataset, CLCD),空間分辨率為30m,一級分類包括耕地、森林、灌木、草地、水域、雪/冰、荒地、建設(shè)用地和濕地.根據(jù)研究需要,將土地利用數(shù)據(jù)重新劃分為耕地、林地、草地、水域、未利用地和建設(shè)用地.數(shù)據(jù)時間為2000年、2010年和2020年.Dem數(shù)據(jù)來源于地理空間數(shù)據(jù)云(http://www.gscloud.cn),分辨率為30m.基于ArcGIS10.0軟件平臺,利用DEM數(shù)據(jù)計算坡度.以土地利用數(shù)據(jù)為基準(zhǔn)進行數(shù)據(jù)坐標(biāo)系和分辨率的統(tǒng)一,嚴(yán)格保證數(shù)據(jù)行列數(shù)一致.

    2 研究方法

    2.1 生態(tài)源地識別

    生態(tài)源地是區(qū)域內(nèi)生境質(zhì)量較高、結(jié)構(gòu)和功能相對穩(wěn)定、對周邊區(qū)域具有重要輻射功能的斑塊,在維持生態(tài)系統(tǒng)穩(wěn)定性、提供生態(tài)系統(tǒng)服務(wù)和促進區(qū)域生態(tài)安全等方面具有重要作用[45-46].MSPA是基于數(shù)學(xué)形態(tài)學(xué)原理而提出,能夠更為精準(zhǔn)的對柵格圖像的空間格局在功能型結(jié)構(gòu)上進行分類,從像元層面上識別出對景觀連通性具有主要作用的區(qū)域[47],從而提高生態(tài)源地選取的科學(xué)性[30].本研究采用MSPA并結(jié)合已有研究進行生態(tài)源地的選擇,將生態(tài)系統(tǒng)服務(wù)價值較高的林地作為MSPA分析的前景數(shù)據(jù),其他地類作為背景數(shù)據(jù);基于Guidos分析軟件,得到互不重疊的七類景觀要素;選取對維持景觀連通性具有重要意義的核心區(qū)作為生態(tài)源地,并剔除面積小于5km2的斑塊.

    2.2 電阻面的構(gòu)建

    在電路理論中,景觀被抽象為具有不同阻力的單元,將物種個體或生態(tài)因素視為電子,阻力面反映了物種和生態(tài)因素在景觀單元上移動時受到的阻礙程度[43,48].本文參考已有研究,選取土地利用、高程、坡度、MSPA景觀類型構(gòu)建阻力面[49-51];在阻力值賦值時,參考已有研究并結(jié)合專家意見確定,并利用層次分析法得到不同阻力因子的權(quán)重.

    2.3 生態(tài)廊道提取

    生態(tài)廊道是生態(tài)源地間物質(zhì)、能量和信息流動與交換的通道、是生態(tài)源地間最容易聯(lián)系的低阻力生態(tài)通道[43].將電路理論的連接模型與隨機游走相結(jié)合能夠較高的評價最小成本路徑[52].本研究基于電路理論,通過電流密度來識別生態(tài)廊道;若該區(qū)域電流密度較高,則表示該生態(tài)阻力較低,連接度較好.

    2.4 生態(tài)“夾點”及障礙點識別

    生態(tài)夾點是廊道中電流密度較大的區(qū)域[53]、是物種遷移交流的必經(jīng)之路或無其他可替代的區(qū)域;若此類區(qū)域生態(tài)環(huán)境遭到破壞或生態(tài)系統(tǒng)發(fā)生退化,極大可能切斷生態(tài)網(wǎng)絡(luò)之間的連通性和完整性.夾點的形成往往是由于周邊地區(qū)阻力值較大,廊道在夾點區(qū)域被壓縮在相對狹窄的范圍內(nèi).本文利用Circuitscape與Linkage Mapper工具箱中的Pinchpoint Mapper模塊來識別生態(tài)廊道中的“夾點”,即電流密度較高的區(qū)域.本文選擇“All to one”模型進行夾點識別.生態(tài)障礙點可視為阻礙生態(tài)源地間連通性的關(guān)鍵節(jié)點[54],移除后可以改善源地之間的連通性.本文利用Linkage Mapper工具箱中的Barrier Mapper模塊來識別障礙點.

    2.5 生態(tài)修復(fù)區(qū)識別

    2.5.1 生態(tài)源地重要性評估 Centrality Mapper是Linkage Mapper Toolkit的一部分,它包括Linkage Mapper和其他旨在支持區(qū)域野生動物棲息地連通性分析的模塊.利用Centrality Mapper分析得到連接網(wǎng)絡(luò),計算當(dāng)前網(wǎng)絡(luò)的生態(tài)流中心性.當(dāng)前生態(tài)流中心性是衡量一個鏈接或核心區(qū)域?qū)τ诒3终麄€網(wǎng)絡(luò)連接的重要程度.最終,在本研究中,為評估生態(tài)源地的重要性,生態(tài)源地中心性被分為低、中、高三個等級,等級越高,生態(tài)源地越重要.

    2.5.2 生態(tài)廊道阻力值評估 利用Linkage Mapper軟件輸入生態(tài)源地和綜合阻力面,計算生態(tài)源地之間的成本加權(quán)距離(CWD)和最小成本路徑(LCP_Length).用CWD/LCP_Length反映最小成本路徑的性質(zhì),當(dāng)CWD/LCP值越大時,表示物種通過此路徑時遭受的阻力較大,景觀連通性較差,反之亦然[55].

    3 結(jié)果與分析

    3.1 土地利用變化

    2000~2020年,太原城市群的土地利用類型以草地和耕地為主(圖1),面積平均占比(表1)分別為34.13%和32.80%,其次是林地,面積平均占比為25.58%,建設(shè)用地和水域面積平均占比分別為7.28%和0.20%,未利用地面積占比最少,不足0.01%. 2000~2020年林地、建設(shè)用地面積呈現(xiàn)不斷增加的趨勢,分別增加了525.20和689.95km2.耕地面積呈現(xiàn)先減少后增加趨勢,草地面積呈現(xiàn)先增加后減少,兩者總體均呈減少趨勢,分別減少885.52和360.93km2.水域和未利用地面積增加較少,分別增加了29.50和1.79km2.

    表1 2000~2020年土地利用類型占比(%)

    圖2 2000~2020年土地利用類型空間分布

    土地利用轉(zhuǎn)移矩陣圖(圖3)顯示了2000~2020年土地利用轉(zhuǎn)移模式.從凈轉(zhuǎn)出來看,耕地轉(zhuǎn)出面積最大,其次為草地,轉(zhuǎn)出面積分別為2011.5和1998.91km2,累計占凈轉(zhuǎn)出面積85.91%.其中,耕地轉(zhuǎn)為草地1129.11km2,轉(zhuǎn)為建設(shè)用地735.43km2.草地主要轉(zhuǎn)為耕地和林地,分別轉(zhuǎn)出988.94和930.51km2.從凈轉(zhuǎn)入來看,草地、耕地、林地、建設(shè)用地轉(zhuǎn)入面積都較多,分別轉(zhuǎn)入1637.98、1125.99、1055.64、和812.14km2,占總轉(zhuǎn)入面積35.09%、24.12%、22.61%、17.40%,其中林地和耕地與草地之間相互轉(zhuǎn)換,建設(shè)用地主要來源于耕地.

    圖3 2000~2020年土地利用轉(zhuǎn)移矩陣圖

    3.2 生態(tài)源地識別

    通過MSPA模型分析,選取核心區(qū)作為生態(tài)源地,并剔除面積小于5km2的斑塊,生態(tài)源地識別結(jié)果如圖4所示.三期分別識別出39、44、48塊生態(tài)源地,面積分別為3083.38km2、3140.64km2、3144.17km2,分別占研究區(qū)總面積的14.92%、15.19%和15.21%.研究區(qū)生態(tài)源地整體分布不均,呈西部多,中、東部少的空間格局.整體而言生態(tài)源地呈空間集聚分布,主要分布于呂梁的交城縣、汾陽市、文水縣及太原的古交市和陽曲縣,累計約占總生態(tài)源地面積的71.74%.小店區(qū)、尖草坪區(qū)、迎澤區(qū)、杏花嶺區(qū)及清徐縣生態(tài)源地面積較少,累計占總面積比例不足1%.原因是受研究區(qū)特殊的地形條件影響,研究區(qū)西部主要為呂梁山區(qū),以大面積的林草地為主,生態(tài)系統(tǒng)結(jié)構(gòu)與功能完整,生態(tài)服務(wù)價值較高.中部為河谷平原,受人為活動擾動較大,主要以農(nóng)田和建設(shè)用地為主.東部雖分布有生態(tài)用地,但整體面積較少,生態(tài)功能作用較弱.

    3.3 生態(tài)廊道提取

    基于Linkage mapper進行研究區(qū)生態(tài)廊道的識別(圖4).經(jīng)統(tǒng)計,2000年、2010年、2020年生態(tài)廊道分別是69、79、90條,總長度分別是1099.74km、1242.62km、1365.00km.廊道的平均長度約15.61km,且呈現(xiàn)下降趨勢.主要是由于受人類活動、退耕還林還草以及草地轉(zhuǎn)林地等的影響,研究區(qū)生態(tài)源地面積減少,數(shù)量增多,從而引起生態(tài)廊道數(shù)目不斷增多,總長度增加,平均長度減少.由圖4可知,研究區(qū)廊道空間差異明顯,西部生態(tài)源地分布相對密集、廊道數(shù)量較多、且以短距離廊道為主、景觀連通性較強.研究區(qū)西部與東部之間距離相對較長,生態(tài)廊道數(shù)量較少.并且受西部生態(tài)源地增加的影響,東西部間發(fā)散狀廊道數(shù)量增多.因此,西部和東部整體景觀連通性受文水縣,晉源區(qū)以及萬柏林區(qū)與古交市交界處斑塊的影響較大,今后應(yīng)對此類生態(tài)源地加大保護.

    3.4 生態(tài)“夾點”及障礙點識別

    利用自然斷點法將電流密度分為三級,選取電流密度最高一級作為夾點區(qū)域(圖5),夾點區(qū)域面積由2000年47.64km2增加至2010年97.98km2,增幅105.67%,2020減少至83.77km2,減少14.50%;總體呈增長趨勢.將生態(tài)夾點與生態(tài)廊道進行疊加分析,發(fā)現(xiàn)研究區(qū)幾乎所有的生態(tài)廊道都存在生態(tài)夾點,且生態(tài)夾點多位于出入生態(tài)源地的位置以及生態(tài)廊道的中間位置.因此,急需對生態(tài)夾點區(qū)域進行保護,以維持區(qū)域整體景觀連通性;一旦夾點區(qū)域斷裂,整體景觀功能極有可能發(fā)生紊亂

    圖4 生態(tài)源地及生態(tài)廊道空間分布

    圖6 生態(tài)障礙點空間分布

    研究區(qū)障礙點面積呈不斷減少的趨勢,由2000年的451.78km2減少至2010年413.94km2,至2020年達(dá)到395.95km2(圖6).疊加生態(tài)夾點和障礙點,發(fā)現(xiàn)在物種遷移的必經(jīng)之路上存在阻礙其流動的區(qū)域,對于此類區(qū)域,可優(yōu)先進行修復(fù).其中,研究區(qū)西部生態(tài)障礙點變化較小,生態(tài)障礙點變化較大的主要集中于中部地區(qū).從生態(tài)障礙點空間分布來看,西部山區(qū)生態(tài)障礙點多分布于生態(tài)廊道的中間或者生態(tài)源地邊緣位置.研究區(qū)中部生態(tài)障礙點沿生態(tài)廊道呈現(xiàn)帶狀分布,分布范圍廣泛.從縣域來看,文水縣、清徐縣、平遙縣、太谷縣、祁縣等縣域障礙點面積始終居于前列,主要原因是隨著生態(tài)源地的增加,導(dǎo)致物種遷移路徑增多,但此類地區(qū)多以農(nóng)田和建設(shè)用地為主,導(dǎo)致受到的阻力較大,極易形成生態(tài)障礙點.

    3.5 生態(tài)重要區(qū)識別

    基于中心性Mapper,計算了各生態(tài)源地的中心性(圖7).利用自然斷點法將生態(tài)源地中心性分為3個等級,中心性等級越高,生態(tài)源地越重要.圖中顯示了研究期間中心性較高的生態(tài)源地分布狀況,中高中心性生態(tài)源地主要分布在呂梁的文水縣、交城縣、汾陽市和太原市的古交市、陽曲縣,主要是由于該地區(qū)生態(tài)源地面積較大.隨著生態(tài)源地破碎化,高中心性生態(tài)源地的面積減少.因此,保護上述地區(qū)生態(tài)源地對促進整體景觀格局的穩(wěn)定具有重要意義.

    為進一步衡量生態(tài)廊道的相對阻力,本文采用CWD/LCP_Length的比值來反映生態(tài)廊道的相對阻力(圖7),比值越大,說明生態(tài)廊道的相對阻力越大,連通性越差.利用自然斷點法將生態(tài)廊道的阻力值分為高、中、低三個等級.高阻力生態(tài)廊道主要為東西向生態(tài)廊道.2020年,阻力值較高的生態(tài)廊道分布在研究區(qū)中部.該地區(qū)生態(tài)廊道的長度較長,并跨越了建筑用地和農(nóng)田等阻力較大的地區(qū).中低阻力生態(tài)廊道主要分布在西部和東部,以西部地區(qū)為主.西部地區(qū)土地利用類型為森林和草地,生態(tài)系統(tǒng)服務(wù)價值較高,阻力價值較低.

    圖7 生態(tài)源地重要性和生態(tài)廊道阻力值的空間格局

    4 討論

    目前生態(tài)廊道的識別多采用MCR模型[34]. MCR模型基于ArcGIS的成本路徑工具識別生態(tài)源地間最小耗費距離作為生態(tài)廊道,卻忽略了生物的隨機游走特性,并且不能明確生態(tài)廊道的具體范圍和關(guān)鍵節(jié)點[56].電路理論基于電子在電路中隨機游走的特性,能夠更為準(zhǔn)確的模擬物種遷移的狀況,識別待保護的“節(jié)點”,待移除的“障礙點”等國土空間生態(tài)修復(fù)的區(qū)域[57].彌補目前單點、單要素、單過程生態(tài)修復(fù)方式的不足,解決生態(tài)修復(fù)活動的空間割裂問題,為系統(tǒng)性和針對性國土空間生態(tài)修復(fù)提供重要的決策參考[20].本研究采用電路理論對生態(tài)廊道、節(jié)點和障礙點進行識別,并提出了具體的修復(fù)方向和措施.

    4.1 加強生態(tài)源地保護

    在生態(tài)安全格局構(gòu)建的基礎(chǔ)上,本文比較了2000年、2010年、2020年三個時間點生態(tài)安全格局變化.結(jié)果表明,太原城市群生態(tài)源地分布不均,呈西部多,中東部少的空間格局.生態(tài)源地集中分布于交城市、汾陽市、文水市、古交市以及陽曲縣等地.生態(tài)源地作為區(qū)域內(nèi)景觀結(jié)構(gòu)和功能較為完整的斑塊,應(yīng)著力進行保護,并逐步提高生態(tài)源地質(zhì)量[43].重點保護整體性較好且中心度較高生態(tài)斑塊,此類型斑塊對于維持區(qū)域生態(tài)系統(tǒng)的穩(wěn)定性,增強區(qū)域整體景觀連通性具有重要作用.

    4.2 生態(tài)廊道修復(fù)措施

    生態(tài)廊道是物種可以遷移的狹長地帶.它們連接棲息地、防止物種隔離,維持最低數(shù)量的種群,并有助于生物多樣性保護[58].受生態(tài)源地數(shù)增加的影響,研究期內(nèi)生態(tài)廊道數(shù)量增多,但廊道阻力值偏大,物種遷移路徑較長.對于西部相對密集的短距離生態(tài)廊道,其土地利用類型主要是林草地,今后應(yīng)以生態(tài)維護為主,以保障區(qū)域景觀生態(tài)流的暢通.研究區(qū)中部為河谷平原,以阻力值較大的農(nóng)田和建設(shè)用地為主;且該區(qū)域是城市擴張的主要區(qū)域,建設(shè)用地不斷侵占耕地,導(dǎo)致區(qū)域阻力值增大.因此,今后應(yīng)防止該地區(qū)建設(shè)用地過快增長[43].在對自然生態(tài)廊道維護的基礎(chǔ)上,加強規(guī)劃廊道的建設(shè).綜合考慮總體規(guī)劃布局以及城市資本效用,積極進行g(shù)reen corridors的建設(shè)、在建筑周邊和必要的生態(tài)區(qū)域種植植被,以加速物種的遷移和傳播[54].對生態(tài)空間已被占用,僅保護現(xiàn)有生態(tài)空間無法維持正常生態(tài)過程的生態(tài)廊道,應(yīng)進行生態(tài)廊道的重建[59].

    4.3 生態(tài)夾點及障礙點修復(fù)措施

    生態(tài)夾點作為表征景觀連通性的關(guān)鍵節(jié)點,其退化或者損失可能會切斷源地間的連通性,應(yīng)優(yōu)先保護[46].研究期內(nèi)生態(tài)夾點呈現(xiàn)不斷增加趨勢,生態(tài)夾點的增加主要有兩方面原因:一方面,受退耕還林還草以及草地轉(zhuǎn)林地等影響,導(dǎo)致經(jīng)過該地區(qū)的物種遷移路徑增多,生態(tài)夾點面積增加.另一方面,由于部分縣域本身生態(tài)源地較多,隨著生態(tài)源地的增加,縣域內(nèi)生態(tài)廊道增加,網(wǎng)絡(luò)密度加大,此類地區(qū)生態(tài)夾點的土地利用類型以林草地為主.研究期內(nèi)生態(tài)障礙點面積不斷減少,主要原因是受城市擴張的影響,導(dǎo)致物種原有遷移發(fā)生斷裂,新增長的遷移路徑受城市擴張等因素影響,導(dǎo)致物種遷移路徑變窄,障礙區(qū)面積隨之減少.由此可見,改善區(qū)域景觀連通性,僅依靠增加生態(tài)源地而不對生態(tài)廊道加以維護,會導(dǎo)致現(xiàn)有廊道阻力值不斷加大,連通性不夠,最終無法實現(xiàn)生態(tài)流的良好流通[59].針對識別出的生態(tài)夾點和障礙點,應(yīng)針對特定區(qū)域生態(tài)環(huán)境問題,應(yīng)用生態(tài)工程措施加以修復(fù),保障生態(tài)系統(tǒng)的健康發(fā)展.

    4.4 持續(xù)推進城市群內(nèi)生態(tài)補償

    城市群作為最具活力和競爭力的城市空間單元,依靠城市之間的職能分工和經(jīng)濟協(xié)調(diào)以實現(xiàn)資源的優(yōu)化配置,以最大化城市群的整體效益[60].由于自然資源稟賦、環(huán)境承載能力和社會經(jīng)濟發(fā)展水平的差異[60-61],使得各城市功能定位的不同,具有經(jīng)濟發(fā)展比較優(yōu)勢的城市往往會過度消耗生態(tài)資源甚至破壞生態(tài)環(huán)境,產(chǎn)生負(fù)外部性,而這種負(fù)外部性需要具有生態(tài)比較優(yōu)勢的城市通過生態(tài)環(huán)境的保護進行緩解,這種行為將會導(dǎo)致該類城市發(fā)展機會成本的增加[62].這種城市之間利益分配的不均若不通過政府的手段或者市場機制進行調(diào)節(jié),將會導(dǎo)致城市群內(nèi)部環(huán)境不公平問題,進而引起城市間不公平的競爭和合作.生態(tài)補償機制可以在一定程度上改善區(qū)域發(fā)展的不平衡,因此應(yīng)秉承“誰受益,誰補償,誰污染,誰付費”的原則,對城市之間的利益關(guān)系進行調(diào)整[63].對于交城、文水、陽曲等生態(tài)資源稟賦較高,能夠為周邊地區(qū)提供生態(tài)產(chǎn)品和服務(wù)的生態(tài)貢獻(xiàn)縣域,應(yīng)該從小店、尖草坪等生態(tài)受益縣域得到生態(tài)補償,并且應(yīng)持續(xù)推進生態(tài)貢獻(xiàn)縣域的生態(tài)管理工作,從而保障生態(tài)服務(wù)服務(wù)供給的可持續(xù)性[64].積極探究補償?shù)亩鄻踊?將區(qū)域的生態(tài)補償與綠色發(fā)展相結(jié)合,優(yōu)先支持綠色產(chǎn)業(yè)和生態(tài)農(nóng)業(yè),促進區(qū)域發(fā)展的良性循環(huán).

    5 結(jié)論

    5.1 研究區(qū)生態(tài)源地和廊道空間分布不均,呈西部多、中東部少的空間格局,西部生態(tài)源地分布相對密集,廊道數(shù)量較多,且以短距離廊道為主;西部與東部之間距離相對較長,生態(tài)廊道數(shù)量較少.

    5.2 源地的增加導(dǎo)致了生態(tài)廊道數(shù)量的增加,但受景觀界面阻力值增大的影響,導(dǎo)致生態(tài)遷移路徑變窄,夾點面積增加.

    5.3 本研究基于對研究區(qū)生態(tài)安全格局的動態(tài)監(jiān)測,提出了針對性的修復(fù)措施,可為區(qū)域景觀規(guī)劃與管理提供參考.

    [1] Jiang L G, Liu Y, Wu S, et al. Analyzing ecological environment change and associated driving factors in China based on NDVI time series data [J]. Ecological Indicators, 2021,129:107933.

    [2] An M, Xie P, He W J, et al. Spatiotemporal change of ecologic environment quality and human interaction factors in three gorges ecologic economic corridor, based on RSEI [J]. Ecological Indicators, 2022,141:109090.

    [3] Zhou L, Dang X W, Sun Q K, et al. Multi-scenario simulation of urban land change in Shanghai by random forest and CA-Markov model [J]. Sustainable Cities and Society, 2020,55:102045.

    [4] Cafaro P, Hansson P, & G?tmark F. Overpopulation is a major cause of biodiversity loss and smaller human populations are necessary to preserve what is left [J]. Biological Conservation, 2022,272:109646.

    [5] Liu Y X, Hou X L, Qiao J X, et al. Evaluation of soil erosion rates in the hilly-gully region of the Loess Plateau in China in the past 60 years using global fallout plutonium [J]. Catena, 2023,220:106666.

    [6] Xiang S J, Wang Y, Deng H, et al. Response and multi-scenario prediction of carbon storage to land use/cover change in the main urban area of Chongqing, China [J]. Ecological Indicators, 2022,142: 109205.

    [7] Ding M M, Liu W, Xiao L, et al. Construction and optimization strategy of ecological security pattern in a rapidly urbanizing region: A case study in central-south China [J]. Ecological Indicators, 2022,136: 108604.

    [8] 吳健生,馬洪坤,彭 建.基于“功能節(jié)點—關(guān)鍵廊道”的城市生態(tài)安全格局構(gòu)建——以深圳市為例 [J]. 地理科學(xué)進展, 2018,37(12): 1663-1671. Wu J S, Ma H K, Peng J. Improving urban ecological security pattern based on functional nodes- key corridors: A case study of Shenzhen City [J]. Progress in Geography, 2018,37(12):1663-1671.

    [9] Peng J, Yang Y, Liu Y X, et al. Linking ecosystem services and circuit theory to identify ecological security patterns [J]. Science of the total environment, 2018,644:781-790.

    [10] 陳 瑾,趙超超,趙 青,等.基于MSPA分析的福建省生態(tài)網(wǎng)絡(luò)構(gòu)建 [J]. 生態(tài)學(xué)報, 2023,43(2):603-614. Chen J, Zhao C C, Zhao Q, et al. Construction of ecological network in Fujian Province based on morphological spatial pattern anylysis [J]. Acta Ecologica Sinica, 2023,43(2):603-614.

    [11] Chen J, Wang S S, & Zou Y T. Construction of an ecological security pattern based on ecosystem sensitivity and the importance of ecological services: A case study of the Guanzhong Plain urban agglomeration, China [J]. Ecological Indicators, 2022,136:108688.

    [12] 王子琳,李志剛,方世明.基于遺傳算法和圖論法的生態(tài)安全格局構(gòu)建與優(yōu)化—以武漢市為例 [J]. 地理科學(xué), 2022,42(10):1685-1694. Wang Z L, Li Z G, Fang S M. Construction and optimization of ecological security pattern based on genetic algorithm and graph theory: A case study of Wuhan City [J]. Scientia Geographica Sinica, 2022,42(10):1685-1694.

    [13] 姚材儀,何艷梅,程建兄,等.基于MCR模型和重力模型的岷江流域生態(tài)安全格局評價與優(yōu)化建議研究 [J]. 生態(tài)學(xué)報, 2023,43(17):1-14. Yao C Y, He Y M, Cheng J X, et al. Evaluation of ecological security pattern and optimization suggestions in Minjiang River Basin based on MCR model and gravity model [J]. Acta Ecologica Sinica, 2023,43 (17):1-14.

    [14] Pan J H, Liang J, Zhao C C. Identification and optimization of ecological security pattern in arid inland basin based on ordered weighted average and ant colony algorithm: A case study of Shule River basin, NW China [J]. Ecological Indicators, 2023,154:110588.

    [15] 鄒珮雯,徐 昉.耦合ERI-MCR-PLUS模型的生態(tài)安全格局構(gòu)建及景觀生態(tài)風(fēng)險預(yù)測研究——以賽罕烏拉國家級自然保護區(qū)為例 [J]. 生態(tài)學(xué)報, 2023,43(23):1-13. Zou P W, Xu F. Ecological security pattern construction and landscape ecological risk prediction by coupling ERI-MCR-PLUS model: A case study of Saihanwula National Nature Reserve [J]. Acta Ecologica Sinica, 2023,43(23):1-13.

    [16] 安 睿,竇 超,陸硯池,等.耦合SOM-MCR模型的多特征生態(tài)安全格局構(gòu)建——以武漢城市圈為例 [J]. 生態(tài)學(xué)報, 2023,43(22):1-14. An R, Dou C, Lu Y C, et al. Construction of multi-feature ecological security patterns by coupling SOM-MCR model: A case study of Wuhan Metropolitan Area [J]. Acta Ecologica Sinica, 2023,43(22): 1-14.

    [17] Qian W Q, Zhao Y, Li X Y. Construction of ecological security pattern in coastal urban areas: A case study in Qingdao, China [J]. Ecological Indicators, 2023,154:110754.

    [18] 趙筱青,譚 琨,易 琦,等.典型高原湖泊流域生態(tài)安全格局構(gòu)建——以杞麓湖流域為例 [J]. 中國環(huán)境科學(xué), 2019,39(2):768-777. Zhao X Q, Tan K, Yi Q, et al. Construction of ecological security pattern in typical plateau lake basin——A case of the Qilu lake basin [J]. China Environmental Science, 2019,39(2):768-777.

    [19] 李洪慶,楊 瑀,張俊紅.雅魯藏布江山南寬谷流域生態(tài)安全格局構(gòu)建 [J]. 干旱區(qū)地理, 2023,46(9):1503-1513. Li H Q, Zhang Y, Zhang J H, et al. Construction of ecological security pattern in Shannan Wide Valley Basin of Yarlung Zangbo River [J]. Aris Land Geography, 2023,46(9):1503-1513.

    [20] 王回茴,李漢廷,謝苗苗,等.資源型城市工礦用地系統(tǒng)修復(fù)的生態(tài)安全格局構(gòu)建 [J]. 自然資源學(xué)報, 2020,35(1):162-173. Wang H H, Li H T, Xie M M, et al. Construction of ecological security pattern for systematic restoration of industrial and mining land in resource-based cities [J]. Journal of Natural Resources, 2020,35(1): 162-173.

    [21] 陳月嬌,李 祥,王月健,等.尺度整合視角下伊犁河谷地區(qū)生態(tài)安全格局構(gòu)建——以昭蘇縣為例 [J]. 生態(tài)學(xué)報, 2023,43(19):1-12. Chen Y J, Li X, Wang Y J, et al. Construction of ecological security pattern in Ili River Valley from the perspective of scale integration: A case of Zhaosu County [J]. Acta Ecologica Sinica, 2023,43(19):1-12.

    [22] 孫 楓,章錦河,王培家,等.城市生態(tài)安全格局構(gòu)建與評價研究:以蘇州市區(qū)為例 [J]. 地理研究, 2021,40(9):2476-2493. Sun F, Zhang J H, Wang P J, et al. Construction and evaluation of urban ecological security pattern: A case study of Suzhou city [J]. Geographical Research, 2021,40(9):2476-2493.

    [23] Mann D, Agrawal G, & Joshi P K. Spatio-temporal forest cover dynamics along road networks in the Central Himalaya [J]. Ecological Engineering, 2019,127:383-393.

    [24] 史芳寧,劉世梁,安 毅,等.城市化背景下景觀破碎化及連接度動態(tài)變化研究——以昆明市為例 [J]. 生態(tài)學(xué)報, 2020,40(10):3303- 3314. Shi F N, Liu S L, An Y, et al. Changes of landscape fragmentation and connectivity with urbanization: a case study of Kunming City [J]. Acta Ecologica Sinica, 2020,40(10):3303-3314.

    [25] 劉 婷,歐陽帥,勾蒙蒙,等.基于MSPA模型的新型城市熱景觀連通性分析 [J]. 生態(tài)學(xué)報, 2023,43(2):615-624. Liu T, Ou Y S, Gou M M, et al. Analysis connectivity of urban heat island in a new-type urbanization based on MSPA model [J]. Acta Ecologica Sinica, 2023,43(2):615-624.15.

    [26] 喬 治,陳嘉悅,王 楠,等.基于MSPA和電路理論的京津冀城市群熱環(huán)境空間網(wǎng)絡(luò) [J]. 環(huán)境科學(xué), 2023,44(6):3034-3042. Qiao Z, Chen J Y, Wang N, et al. Spatial network of urban heat environment in Beijing-Tianjin-Hebei Urban Agglomeration based on MSPA and circuit theory [J]. Environmental Science, 2023,44(6): 3034-3042.

    [27] 馬 銀,鄭敏睿,鄭新奇,等.基于CA-Markov和MSPA的綠色基礎(chǔ)設(shè)施預(yù)測與時空演變分析研究——以京津冀城市群為例 [J]. 生態(tài)學(xué)報, 2023,43(16):6785-6797. Ma Y, Zheng C R, Zheng X Q, et al. Prediction and spatiotemporal evolution analysis of green infrastructure based on CA-Markov and MSPA: A case study of Beijing-Tianjin-Hebei urban agglomeration [J]. Acta Ecologica Sinica, 2023,43(16):6785-6797.

    [28] 汪勇政,李久林,顧康康,等.基于形態(tài)學(xué)空間格局分析法的城市綠色基礎(chǔ)設(shè)施網(wǎng)絡(luò)格局優(yōu)化——以合肥市為例 [J]. 生態(tài)學(xué)報, 2022, 42(5):2022-2032. Wang Y Z, Li J L, Gu K K, et al. Optimization of urban green infrastructure network layout based on MSPA-CIRCUIT: Case of Hefei [J]. Acta Ecologica Sinica, 2022,42(5):2022-2032.

    [29] 秦子博,玄 錦,黃柳菁,等.基于MSPA和MCR模型的海島型城市生態(tài)網(wǎng)絡(luò)構(gòu)建——以福建省平潭島為例 [J]. 水土保持研究, 2023, 30(2):303-311. Qin Z B, Xuan J, Huang L J, et al. Ecological network construction of sea-island city based on MSPA and MCR Model [J]. Research of Soil and Water Conservation, 2023,30(2):303-311.

    [30] 陳竹安,馬彬彬,危小建,等.基于MSPA和MCR模型的南昌市生態(tài)網(wǎng)絡(luò)構(gòu)建與優(yōu)化 [J]. 水土保持通報, 2021,41(6):139-147. Chen Z A, Ma B B, Wei X J, et al. Construction and optimization ecological network of Nanchang city based on MSPA and MCR Model [J]. Bulletin of Soil and Water Conservation, 2021,41(6):139-147.

    [31] 王雪然,萬榮榮,潘佩佩.太湖流域生態(tài)安全格局構(gòu)建與調(diào)控——基于空間形態(tài)學(xué)-最小累積阻力模型 [J]. 生態(tài)學(xué)報, 2022,42(5): 1968-1980. Wang X R, Wan R R, Pan P P. Construction and adjustment of ecological security pattern based on MSPA-MCR Model in Taihu Lake Basin [J]. Acta Ecologica Sinica, 2022,42(5):1968-1980.

    [32] 楊 凱,曹銀貴,馮 喆,等.基于最小累積阻力模型的生態(tài)安全格局構(gòu)建研究進展 [J]. 生態(tài)與農(nóng)村環(huán)境學(xué)報, 2021,37(5):555-565. Yang K, Cao Y G, Feng Z, et al. Research progress of ecological security pattern construction based on minimum cumulative resistance model [J]. Journal of Ecology and Rural Environment, 2021,37(5): 555-565.

    [33] Adriaensen F, Chardon J P, De Blust G, et al. The application of ‘least-cost’modelling as a functional landscape model [J]. Landscape and urban planning, 2003,64(4):233-247.

    [34] 潘竟虎,王 云.基于CVOR和電路理論的討賴河流域生態(tài)安全評價及生態(tài)格局優(yōu)化 [J]. 生態(tài)學(xué)報, 2021,41(7):2582-2595. Pan J H, Wang Y. Ecological security evaluation and ecological pattern optimization in Taolai River Basin based on CVOR and circuit theory [J]. Acta Ecologica Sinica, 2021,41(7):2582-2595.

    [35] Dickson B G, Albano C M, Anantharaman R, et al. Circuit-theory applications to connectivity science and conservation [J]. Conservation Biology, 2019,33(2):239-249.

    [36] Gao M W, Hu Y C, & Bai Y P. Construction of ecological security pattern in national land space from the perspective of the community of life in mountain, water, forest, field, lake and grass: A case study in Guangxi Hechi, China [J]. Ecological Indicators, 2022,139:108867.

    [37] 魏 樂,周 亮,孫東琪,等.黃河流域城市群擴張的時空格局演化及情景模擬——以呼包鄂榆城市群為例 [J]. 地理研究, 2022,41(6): 1610-1622. Wei L, Zhou L, Sun D Q, et al. The evolution of spatio-temporal pattern and scenario simulation of urban agglomeration expansion in the Yellow River Basin: A case study in the Hohhot-Baotou-Ordos- Yulin Urban Agglomeration [J]. Geographical Research, 2022,41(6): 1610-1622.

    [38] Li Z T, Li M, & Xia B C. Spatio-temporal dynamics of ecological security pattern of the Pearl River Delta urban agglomeration based on LUCC simulation [J]. Ecological Indicators, 2020,114:106319.

    [39] Peng C C, Li B, & Nan B. An analysis framework for the ecological security of urban agglomeration: A case study of the Beijing- Tianjin-Hebei urban agglomeration [J]. Journal of Cleaner Production, 2021,315:128111.

    [40] Zhang D, Wang X R, Qu L P, et al. Land use/cover predictions incorporating ecological security for the Yangtze River Delta region, China [J]. Ecological Indicators, 2020,119:106841.

    [41] 郭彥君,郭文炯.“三生空間”視角下山西中部盆地城市群景觀生態(tài)風(fēng)險分析 [J]. 生態(tài)學(xué)雜志, 2022,41(9):1813-1824. Guo Y J, Guo W J. Landscape ecological risk analysis of urban agglomeration in the central basin of Shanxi from the perspective of “production-living-ecological space” [J]. Chinese Journal of Ecology, 2022,41(9):1813-1824.

    [42] 閆麗沙.太原城市圈城市生態(tài)系統(tǒng)健康評價 [D]. 太原:山西師范大學(xué), 2019. Yan L S. Health evaluation of urban ecosystem in Taiyuan city circle [D]. Taiyuan: Shanxi Normal University, 2019.

    [43] Li Q, Zhou Y, & Yi S Q. An integrated approach to constructing ecological security patterns and identifying ecological restoration and protection areas: a case study of Jingmen, China [J]. Ecological Indicators, 2022,137:108723.

    [44] Ran Y J, Lei D M, Li J, et al. Identification of crucial areas of territorial ecological restoration based on ecological security pattern: A case study of the central Yunnan urban agglomeration, China [J]. Ecological Indicators, 2022,143:109318.

    [45] 衛(wèi)新東,林良國,馮小龍,等.神木市生態(tài)安全格局構(gòu)建與生態(tài)問題定量診斷 [J]. 生態(tài)學(xué)報, 2023,43(1):82-94. Wei X D, Lin L G, Feng X L, et al. Construction of ecological security pattern and quantitative diagnosis of ecological problems in Shenmu City [J]. Acta Ecologica Sinica, 2023,43(1):82-94.

    [46] 周 璟,王宏衛(wèi),談 波,等.開都河流域生態(tài)安全格局構(gòu)建與生態(tài)修復(fù)分區(qū)識別 [J]. 生態(tài)學(xué)報, 2022,42(24):10127-10137. Zhou J, Wang H W, Tan B, et al. Construction of ecological security pattern and identification of ecological rehabilitation zones in Kaidu River Basin [J]. Acta Ecologica Sinica, 2022,42(24):10127-10137.

    [47] Vogt P, Riitters K H, Estreguil C, et al. Mapping spatial patterns with morphological image processing [J]. Landscape ecology, 2007,22: 171-177.

    [48] Ding M M, Liu W, Xiao L, et al. Construction and optimization strategy of ecological security pattern in a rapidly urbanizing region: A case study in central-south China [J]. Ecological Indicators, 2022,136: 108604.

    [49] Wang C X, Yu C Y, Chen T Q, et al. Can the establishment of ecological security patterns improve ecological protection? An example of Nanchang, China [J]. Science of the Total Environment, 2020,740:140051.

    [50] 曹秀鳳,劉兆順,李淑杰,等.基于生態(tài)安全格局的國土空間生態(tài)修復(fù)關(guān)鍵區(qū)域識別——以吉林省松原市為例 [J]. 中國環(huán)境科學(xué), 2022, 42(6):2779-2787. Cao X F, Liu Z S, Li S, et al. Identification of key areas of ecological protection and restoration based on the pattern of ecological security: A case of Songyuan City, Jilin province [J]. China Environmental Science, 2022, 42(6):2779-2787.

    [51] Dai L, Liu Y, Luo X. Integrating the MCR and DOI models to construct an ecological security network for the urban agglomeration around Poyang Lake, China [J]. Science of The Total Environment, 2021,754:141868.

    [52] Doyle P G, & Snell J L. Random walks and electric networks [M]. American Mathematical Soc, 1984.

    [53] 薛 強,路 路,牛 韌,等.基于地質(zhì)災(zāi)害敏感性的生態(tài)安全格局關(guān)鍵區(qū)識別與修復(fù)——以濟南市為例 [J]. 生態(tài)學(xué)報, 2021,41(22): 9050-9063. Xue Q, Lu L, Niu R, et al. Identification and restoration of key areas of ecological security pattern based on sensitivity to geological disasters: A case study of Jinan City [J]. Acta Ecologica Sinica, 2021, 41(22):9050-9063.

    [54] Fan F F, Wen X J, Feng Z M, et al. Optimizing urban ecological space based on the scenario of ecological security patterns: The case of central Wuhan, China [J]. Applied Geography, 2022,138:102619.

    [55] Li Y G, Liu W, Feng Q, et al. The role of land use change in affecting ecosystem services and the ecological security pattern of the Hexi Regions, Northwest China [J]. Science of The Total Environment, 2022,158940.

    [56] 毛誠瑞,代力民,齊 麟,等.基于生態(tài)系統(tǒng)服務(wù)的流域生態(tài)安全格局構(gòu)建——以遼寧省遼河流域為例 [J]. 生態(tài)學(xué)報, 2020,40(18): 6486-6494. Mao C R, Dai L M, Qi L, et al. Constructing ecological security pattern based on ecosystem services: A case study in Liaohe River Basin, Liaoning Province, China [J]. Acta Ecologica Sinica, 2020, 40(18):6486-6494.

    [57] 袁 媛,白中科,師學(xué)義,等.基于生態(tài)安全格局的國土空間生態(tài)保護修復(fù)優(yōu)先區(qū)確定——以河北省遵化市為例 [J]. 生態(tài)學(xué)雜志, 2022, 41(4):750-759. Yuan Y, Bai Z K, Shi X Y, et al. Determining priority areas for ecosystem preservation and restoration of territory based on ecological security pattern: A case study in Zunhua City, Hebei Province [J]. Chinese Journal of Ecology, 2022,41(4):750-759.

    [58] Wei Q Q, Halike A, Yao K X, et al. Construction and optimization of ecological security pattern in Ebinur Lake Basin based on MSPA-MCR models [J]. Ecological Indicators, 2022,138:108857.

    [59] Li P X, Cao H, Sun W, et al. Quantitative evaluation of the rebuilding costs of ecological corridors in a highly urbanized city: The perspective of land use adjustment [J]. Ecological Indicators, 2022, 141:109130.

    [60] Yang Y, Zhang Y Y, Yang H, et al. Horizontal ecological compensation as a tool for sustainable development of urban agglomerations: Exploration of the realization mechanism of Guanzhong Plain urban agglomeration in China [J]. Environmental Science & Policy, 2022, 137:301-313.

    [61] Gao X, Shen J Q, He W, et al. Spatial-temporal analysis of ecosystem services value and research on ecological compensation in Taihu Lake Basin of Jiangsu Province in China from 2005 to 2018 [J]. Journal of Cleaner Production, 2021,317:128241.

    [62] Zhou Z X, Sun X R, Zhang X T, et al. Inter-regional ecological compensation in the Yellow River Basin based on the value of ecosystem services [J]. Journal of Environmental Management, 2022, 322:116073.

    [63] Jiang Y N, Guan D J, He X J, et al. Quantification of the coupling relationship between ecological compensation and ecosystem services in the Yangtze River Economic Belt, China [J]. Land Use Policy, 2022, 114:105995.

    [64] Zhai T L, Zhang D, & Zhao C C. How to optimize ecological compensation to alleviate environmental injustice in different cities in the Yellow River Basin? A case of integrating ecosystem service supply, demand and flow [J]. Sustainable Cities and Society, 2021,75: 103341.

    Analysis on the spatio-temporal evolution of ecological security patterns: A case study over Taiyuan Urban Agglomeration.

    QIQiang-qiang, XU Zhan-jun*

    (College of Resources and Environment, Shanxi Agricultural University, Taigu 030801, China)., 2023,43(11):5987~5997

    Taking the Taiyuan urban agglomeration as an instance, Morphological Spatial Pattern Analysis (MSPA) was applied to explore the changes in ecological source areas in the years of 2000, 2010, and 2020, and circuit theory was adopted to identify ecological corridors, ecological pinch points, and ecological barrier points, so as to construct a regional ecological security pattern, clarify the evolution trend of ecological security pattern, recognize the key areas for ecological restoration, and formulate targeted restoration measures. The research showed that the ecological source area of the study area increased from 3083.38km2in 2000 to 3144.17km2in 2020 with time, which faced the challenge of uneven spatial distribution, showing a spatial structure of sufficiency in the west and inadequacy in the central and eastern regions. Influenced by the increased ecological source areas, the number of ecological corridors in the region increased from 69 in 2000 to 90 in 2020, enhancing landscape connectivity, of which, the western ecological corridors showed the features of large quantities and short distances, while the eastern and western ecological corridors were of long lengths and high resistance. During the research period, the overall area of ecological pinch points showed an increasing trend, with an increased area of 47.64km2in 2000 to 97.98km2in 2010 and a decreased area of 83.77km2; the area of barrier points showed a decreasing trend, with a decreased area of 451.78km2in 2000 to 395.95km2in 2020. Those areas should be deemed as a crucial area for ecological protection and restoration. This study was of great significance, providing guidance for regional ecological restoration and landscape planning.

    morphological spatial pattern analysis model;circuit theory;ecological security patterns;ecological restoration

    X87

    A

    1000-6923(2023)11-5987-11

    祁強強(1996-),女,甘肅平?jīng)鋈?山西農(nóng)業(yè)大學(xué)碩士研究生,主要從事土地資源管理及低碳管理方面研究.q15095538261@126.com.

    祁強強,徐占軍.生態(tài)安全格局時空演變分析——以太原城市群為例 [J]. 中國環(huán)境科學(xué), 2023,43(11):5987-5997.

    QIQiang-qiang, XU Zhan-jun. Analysis on the spatio-temporal evolution of ecological security patterns: A case study over Taiyuan Urban Agglomeration [J]. China Environmental Science, 2023,43(11):5987-5997.

    2023-03-03

    國家自然科學(xué)基金資助項目(51304130);山西省人民政府重大決策咨詢課題(ZB20211703);山西省高等學(xué)校哲學(xué)社會科學(xué)研究項目(201803010);山西省軟科學(xué)研究計劃項目(2018041060-2);山西省哲學(xué)社會科學(xué)規(guī)劃課題(2020YJ052);山西省基礎(chǔ)研究計劃項目(20210302123403)

    * 責(zé)任作者, 教授, zjxu163@126.com

    猜你喜歡
    夾點源地連通性
    ORC系統(tǒng)蒸發(fā)器夾點溫差的數(shù)值分析和實驗研究
    偏序集及其相關(guān)拓?fù)涞倪B通性?
    擬莫比烏斯映射與擬度量空間的連通性
    移民與文化認(rèn)同:土家族民歌《吳幺姑》探析
    河道-灘區(qū)系統(tǒng)連通性評價研究
    發(fā)源地
    氣候變化對渭河源地水文環(huán)境影響分析與探討
    高穩(wěn)定被動群集車聯(lián)網(wǎng)連通性研究
    不同種源地漆樹種子生物學(xué)特性研究
    機械A(chǔ)uto CAD教學(xué)中夾點操作的探究
    午夜精品久久久久久毛片777| 欧美日本中文国产一区发布| 老熟女久久久| 日韩免费高清中文字幕av| 亚洲欧洲日产国产| 麻豆国产av国片精品| 91av网站免费观看| 欧美午夜高清在线| 精品人妻1区二区| 97人妻天天添夜夜摸| 免费在线观看日本一区| 老司机福利观看| 视频在线观看一区二区三区| 国产成人精品在线电影| 久久毛片免费看一区二区三区| 99久久人妻综合| 女同久久另类99精品国产91| 99国产精品99久久久久| 久久久久精品国产欧美久久久| 日本撒尿小便嘘嘘汇集6| 国产成人影院久久av| 欧美日韩亚洲综合一区二区三区_| 精品乱码久久久久久99久播| 国产在线观看jvid| 午夜精品国产一区二区电影| av免费在线观看网站| 一区二区三区乱码不卡18| 91老司机精品| 2018国产大陆天天弄谢| 国产精品美女特级片免费视频播放器 | 精品人妻1区二区| 天堂俺去俺来也www色官网| 夜夜骑夜夜射夜夜干| 免费在线观看日本一区| 丰满人妻熟妇乱又伦精品不卡| 国产区一区二久久| 欧美日韩视频精品一区| cao死你这个sao货| 十八禁网站免费在线| 亚洲国产看品久久| 人人妻,人人澡人人爽秒播| 999久久久精品免费观看国产| 亚洲一区中文字幕在线| 精品少妇内射三级| 黄色怎么调成土黄色| 国产精品 欧美亚洲| 精品第一国产精品| 久久精品国产综合久久久| 伦理电影免费视频| 精品欧美一区二区三区在线| 一进一出抽搐动态| 久久九九热精品免费| 最近最新中文字幕大全免费视频| videosex国产| 亚洲久久久国产精品| 国产精品香港三级国产av潘金莲| 黄色 视频免费看| 亚洲精品久久午夜乱码| 麻豆国产av国片精品| 国产免费福利视频在线观看| 免费高清在线观看日韩| 在线天堂中文资源库| 亚洲成人国产一区在线观看| 亚洲精品久久成人aⅴ小说| 精品亚洲成国产av| 美女视频免费永久观看网站| 久久国产精品影院| 成年人午夜在线观看视频| 午夜91福利影院| 99国产精品99久久久久| xxxhd国产人妻xxx| 最新美女视频免费是黄的| 国产男女内射视频| 18禁裸乳无遮挡动漫免费视频| 国产成人一区二区三区免费视频网站| 国产欧美日韩一区二区精品| 桃红色精品国产亚洲av| 久久天躁狠狠躁夜夜2o2o| 亚洲av第一区精品v没综合| 黄色视频,在线免费观看| 日韩三级视频一区二区三区| 亚洲av美国av| av不卡在线播放| 人人澡人人妻人| 搡老乐熟女国产| 国产高清videossex| 亚洲avbb在线观看| 久久久久久亚洲精品国产蜜桃av| 别揉我奶头~嗯~啊~动态视频| 欧美成狂野欧美在线观看| 国产成人影院久久av| 欧美在线一区亚洲| 国产福利在线免费观看视频| 免费观看av网站的网址| 亚洲性夜色夜夜综合| 国产av精品麻豆| 午夜福利免费观看在线| 女人久久www免费人成看片| 99在线人妻在线中文字幕 | 亚洲,欧美精品.| a级片在线免费高清观看视频| 久久久精品免费免费高清| 一区二区日韩欧美中文字幕| 在线观看免费视频网站a站| 国产真人三级小视频在线观看| 亚洲欧美色中文字幕在线| 丁香六月天网| 丁香六月天网| 免费黄频网站在线观看国产| 十八禁高潮呻吟视频| 极品少妇高潮喷水抽搐| 国产日韩欧美亚洲二区| 国产日韩欧美亚洲二区| 国产野战对白在线观看| 丝袜喷水一区| 日本欧美视频一区| 国产淫语在线视频| 国产av一区二区精品久久| 午夜福利一区二区在线看| 亚洲视频免费观看视频| 国产主播在线观看一区二区| 少妇被粗大的猛进出69影院| kizo精华| 99精品欧美一区二区三区四区| 国产精品二区激情视频| 捣出白浆h1v1| 美女高潮到喷水免费观看| 国产精品.久久久| 色94色欧美一区二区| 国产一区二区三区在线臀色熟女 | www日本在线高清视频| 下体分泌物呈黄色| 亚洲欧美日韩另类电影网站| 亚洲三区欧美一区| 999久久久精品免费观看国产| 少妇粗大呻吟视频| 丝瓜视频免费看黄片| 人人妻,人人澡人人爽秒播| 丝袜喷水一区| 国产淫语在线视频| 国产不卡av网站在线观看| 国产91精品成人一区二区三区 | 黄色丝袜av网址大全| 国产在线精品亚洲第一网站| 亚洲欧美日韩高清在线视频 | 亚洲精品美女久久久久99蜜臀| www.999成人在线观看| 女性生殖器流出的白浆| 久久国产精品大桥未久av| 18禁美女被吸乳视频| 亚洲黑人精品在线| 嫁个100分男人电影在线观看| 69精品国产乱码久久久| bbb黄色大片| 最新在线观看一区二区三区| 婷婷丁香在线五月| 成人三级做爰电影| 色尼玛亚洲综合影院| 久久久久久久久免费视频了| 国产精品电影一区二区三区 | 99九九在线精品视频| 亚洲自偷自拍图片 自拍| 午夜激情av网站| 亚洲精品在线美女| 狠狠精品人妻久久久久久综合| 亚洲少妇的诱惑av| 精品一品国产午夜福利视频| 人妻一区二区av| 色94色欧美一区二区| 久久99一区二区三区| 十八禁人妻一区二区| 国产精品美女特级片免费视频播放器 | 99国产极品粉嫩在线观看| 国产又爽黄色视频| a级片在线免费高清观看视频| 国产午夜精品久久久久久| 亚洲九九香蕉| 欧美成人免费av一区二区三区 | 男女边摸边吃奶| 欧美人与性动交α欧美精品济南到| 色婷婷av一区二区三区视频| 欧美另类亚洲清纯唯美| 精品一区二区三区av网在线观看 | 日韩免费av在线播放| 在线观看人妻少妇| 久久毛片免费看一区二区三区| 黑人操中国人逼视频| 精品一区二区三区四区五区乱码| 久久久国产一区二区| 两个人看的免费小视频| 精品国产国语对白av| 水蜜桃什么品种好| 中文字幕制服av| 久久国产精品男人的天堂亚洲| 91麻豆精品激情在线观看国产 | 天天躁日日躁夜夜躁夜夜| 97人妻天天添夜夜摸| 女人被躁到高潮嗷嗷叫费观| 亚洲少妇的诱惑av| 他把我摸到了高潮在线观看 | 99精品在免费线老司机午夜| 中文字幕精品免费在线观看视频| 亚洲三区欧美一区| 女人高潮潮喷娇喘18禁视频| 久热这里只有精品99| 国产一区二区激情短视频| 黄频高清免费视频| 十八禁人妻一区二区| 日韩一卡2卡3卡4卡2021年| 高清在线国产一区| 色综合婷婷激情| 亚洲一区中文字幕在线| 欧美 日韩 精品 国产| 国产真人三级小视频在线观看| www.999成人在线观看| 狠狠狠狠99中文字幕| 国产熟女午夜一区二区三区| 欧美一级毛片孕妇| 人人妻人人添人人爽欧美一区卜| 在线观看舔阴道视频| 在线观看66精品国产| 黑人巨大精品欧美一区二区蜜桃| 80岁老熟妇乱子伦牲交| 成人免费观看视频高清| 欧美精品人与动牲交sv欧美| 99热国产这里只有精品6| 国产成人啪精品午夜网站| 一区二区三区国产精品乱码| 国产精品免费视频内射| 嫩草影视91久久| 久久精品国产亚洲av香蕉五月 | 叶爱在线成人免费视频播放| 在线观看免费日韩欧美大片| 成人国产av品久久久| 亚洲全国av大片| 另类亚洲欧美激情| 国产一区二区三区在线臀色熟女 | 亚洲国产欧美日韩在线播放| 精品一区二区三卡| 黄色视频不卡| 一级片免费观看大全| 亚洲伊人久久精品综合| 亚洲精品国产色婷婷电影| 午夜91福利影院| 妹子高潮喷水视频| 国产又爽黄色视频| 日韩熟女老妇一区二区性免费视频| 久久久久久久久久久久大奶| aaaaa片日本免费| 欧美久久黑人一区二区| www日本在线高清视频| 成人国产一区最新在线观看| 欧美老熟妇乱子伦牲交| 日本a在线网址| 精品国产一区二区久久| 亚洲情色 制服丝袜| 一边摸一边抽搐一进一出视频| 亚洲专区中文字幕在线| kizo精华| 最近最新中文字幕大全免费视频| 三级毛片av免费| 正在播放国产对白刺激| av不卡在线播放| 亚洲av成人不卡在线观看播放网| 欧美激情 高清一区二区三区| 久久人妻熟女aⅴ| 久久这里只有精品19| 久久国产精品人妻蜜桃| 国产老妇伦熟女老妇高清| 亚洲av成人不卡在线观看播放网| 午夜福利,免费看| 色精品久久人妻99蜜桃| 99久久精品国产亚洲精品| 国产成人av激情在线播放| 国产精品1区2区在线观看. | 不卡av一区二区三区| 久久人人97超碰香蕉20202| 国产91精品成人一区二区三区 | 老司机深夜福利视频在线观看| 欧美日韩黄片免| 国产免费av片在线观看野外av| 天堂中文最新版在线下载| 国产麻豆69| 一级毛片精品| 精品国产亚洲在线| 久久精品熟女亚洲av麻豆精品| 国内毛片毛片毛片毛片毛片| 99热国产这里只有精品6| 精品第一国产精品| 亚洲专区国产一区二区| 日韩 欧美 亚洲 中文字幕| 国产精品av久久久久免费| 18禁黄网站禁片午夜丰满| 日本av手机在线免费观看| 一级黄色大片毛片| 巨乳人妻的诱惑在线观看| 大陆偷拍与自拍| 欧美人与性动交α欧美精品济南到| 久久久久精品国产欧美久久久| 亚洲午夜精品一区,二区,三区| 欧美av亚洲av综合av国产av| 免费高清在线观看日韩| 亚洲人成电影免费在线| 亚洲欧美日韩高清在线视频 | 午夜福利在线观看吧| 777米奇影视久久| 欧美人与性动交α欧美软件| 人人妻,人人澡人人爽秒播| 国产一区二区激情短视频| 麻豆国产av国片精品| 国产成人av激情在线播放| 黄网站色视频无遮挡免费观看| 考比视频在线观看| 一级片免费观看大全| 国产亚洲欧美精品永久| 在线十欧美十亚洲十日本专区| 日韩成人在线观看一区二区三区| 午夜福利在线免费观看网站| 女警被强在线播放| 啦啦啦在线免费观看视频4| 精品国产亚洲在线| 1024香蕉在线观看| 最新的欧美精品一区二区| 久久久久精品国产欧美久久久| 欧美精品亚洲一区二区| 亚洲精品久久成人aⅴ小说| 亚洲精品中文字幕一二三四区 | 少妇被粗大的猛进出69影院| 18在线观看网站| 黑人操中国人逼视频| 国产一区有黄有色的免费视频| 成人黄色视频免费在线看| 国产人伦9x9x在线观看| 一级片'在线观看视频| 日日夜夜操网爽| 在线观看www视频免费| 91麻豆精品激情在线观看国产 | 在线观看免费视频日本深夜| 国产一卡二卡三卡精品| a级毛片在线看网站| 最近最新中文字幕大全电影3 | 日本a在线网址| 精品午夜福利视频在线观看一区 | 在线观看免费视频网站a站| 69av精品久久久久久 | 亚洲欧美色中文字幕在线| 香蕉丝袜av| 精品久久久久久电影网| 色综合婷婷激情| 嫩草影视91久久| 国产成人精品无人区| 桃花免费在线播放| 成年人午夜在线观看视频| 国产成人系列免费观看| 老汉色av国产亚洲站长工具| 宅男免费午夜| 精品免费久久久久久久清纯 | 国产精品国产av在线观看| 18禁美女被吸乳视频| 性少妇av在线| 在线观看免费高清a一片| 精品国产超薄肉色丝袜足j| 考比视频在线观看| 亚洲国产精品一区二区三区在线| 亚洲伊人久久精品综合| 一级毛片精品| 不卡一级毛片| 宅男免费午夜| 中亚洲国语对白在线视频| 狂野欧美激情性xxxx| 亚洲国产欧美网| 精品一区二区三区视频在线观看免费 | 国产淫语在线视频| 亚洲av日韩在线播放| 国产精品免费大片| 欧美精品亚洲一区二区| 一进一出好大好爽视频| 十八禁人妻一区二区| 91老司机精品| 精品福利观看| 国产精品麻豆人妻色哟哟久久| 99精国产麻豆久久婷婷| 精品熟女少妇八av免费久了| 亚洲精品久久午夜乱码| 丝袜在线中文字幕| 一级毛片电影观看| 精品一区二区三区视频在线观看免费 | 国产免费福利视频在线观看| 午夜激情久久久久久久| 在线观看舔阴道视频| 纯流量卡能插随身wifi吗| 欧美黑人欧美精品刺激| 国产色视频综合| 亚洲欧洲日产国产| 女性被躁到高潮视频| 久久久国产欧美日韩av| 午夜福利在线免费观看网站| 99热国产这里只有精品6| 精品国产乱子伦一区二区三区| 精品少妇一区二区三区视频日本电影| 国产又色又爽无遮挡免费看| 久久亚洲真实| 夫妻午夜视频| 一二三四在线观看免费中文在| 两个人免费观看高清视频| 天天操日日干夜夜撸| 国产激情久久老熟女| 777久久人妻少妇嫩草av网站| 一区福利在线观看| 纵有疾风起免费观看全集完整版| 欧美成人免费av一区二区三区 | 免费在线观看黄色视频的| 亚洲精品av麻豆狂野| 黑人巨大精品欧美一区二区蜜桃| 色94色欧美一区二区| 极品人妻少妇av视频| 国产亚洲精品一区二区www | 999久久久精品免费观看国产| av欧美777| 国产黄色免费在线视频| 精品第一国产精品| 欧美亚洲日本最大视频资源| 99国产精品免费福利视频| 久久久国产成人免费| 一本大道久久a久久精品| 一区二区日韩欧美中文字幕| 国产黄频视频在线观看| 中文字幕高清在线视频| 国产又爽黄色视频| 久久精品成人免费网站| 亚洲精品av麻豆狂野| 欧美在线黄色| 国产欧美日韩一区二区精品| 中文字幕精品免费在线观看视频| 国产不卡av网站在线观看| 在线观看免费午夜福利视频| 757午夜福利合集在线观看| 久久av网站| 桃花免费在线播放| 日韩一区二区三区影片| 免费高清在线观看日韩| 久久久精品94久久精品| 91字幕亚洲| a级毛片在线看网站| 男男h啪啪无遮挡| 久久久精品免费免费高清| 国产极品粉嫩免费观看在线| 老熟妇乱子伦视频在线观看| 免费看a级黄色片| 国产精品.久久久| 色综合欧美亚洲国产小说| 黄色视频在线播放观看不卡| 国产精品免费一区二区三区在线 | 亚洲九九香蕉| 青青草视频在线视频观看| 丁香欧美五月| 成人黄色视频免费在线看| 色老头精品视频在线观看| 亚洲精品粉嫩美女一区| 午夜久久久在线观看| 韩国精品一区二区三区| a级毛片在线看网站| 亚洲午夜理论影院| 免费在线观看完整版高清| av片东京热男人的天堂| 国产日韩欧美视频二区| 少妇裸体淫交视频免费看高清 | 久久热在线av| 国产成人系列免费观看| 99精品久久久久人妻精品| 国产精品1区2区在线观看. | 日本黄色视频三级网站网址 | 欧美精品一区二区大全| 丁香六月欧美| 蜜桃在线观看..| 亚洲第一欧美日韩一区二区三区 | 黄色视频不卡| 在线播放国产精品三级| 亚洲美女黄片视频| 一区在线观看完整版| 欧美精品亚洲一区二区| 极品教师在线免费播放| 黄色怎么调成土黄色| 久久精品aⅴ一区二区三区四区| 一区二区三区精品91| 午夜福利视频精品| 黑人巨大精品欧美一区二区蜜桃| 欧美亚洲 丝袜 人妻 在线| 色视频在线一区二区三区| 91字幕亚洲| 久久99一区二区三区| 国产成人免费无遮挡视频| 日韩 欧美 亚洲 中文字幕| 精品一区二区三区四区五区乱码| 久久久精品区二区三区| 交换朋友夫妻互换小说| 国产亚洲精品一区二区www | 成年人免费黄色播放视频| 777米奇影视久久| 亚洲av电影在线进入| 青青草视频在线视频观看| 国产精品香港三级国产av潘金莲| 国产成人av教育| 欧美日韩黄片免| 麻豆av在线久日| 在线十欧美十亚洲十日本专区| 国产精品一区二区在线观看99| 色视频在线一区二区三区| 少妇被粗大的猛进出69影院| 亚洲国产欧美在线一区| av线在线观看网站| 欧美人与性动交α欧美精品济南到| 又紧又爽又黄一区二区| 男女下面插进去视频免费观看| 天天躁日日躁夜夜躁夜夜| av有码第一页| 午夜免费鲁丝| 91麻豆精品激情在线观看国产 | 国产欧美日韩一区二区三| 精品久久蜜臀av无| 精品亚洲成a人片在线观看| 日韩有码中文字幕| 老汉色∧v一级毛片| 汤姆久久久久久久影院中文字幕| 在线观看舔阴道视频| av有码第一页| 免费在线观看日本一区| kizo精华| 国产片内射在线| 久久亚洲精品不卡| 老司机福利观看| 丁香欧美五月| 性少妇av在线| 99久久人妻综合| 中文字幕人妻熟女乱码| 蜜桃国产av成人99| 十八禁人妻一区二区| 水蜜桃什么品种好| 99国产精品一区二区蜜桃av | 黄色a级毛片大全视频| 在线观看一区二区三区激情| 精品国产一区二区久久| 999久久久国产精品视频| 别揉我奶头~嗯~啊~动态视频| 蜜桃国产av成人99| 欧美精品一区二区免费开放| 中文亚洲av片在线观看爽 | 视频在线观看一区二区三区| 黄色视频在线播放观看不卡| 亚洲黑人精品在线| 老熟妇乱子伦视频在线观看| 国产成人精品久久二区二区91| 欧美性长视频在线观看| 性色av乱码一区二区三区2| 久久精品亚洲精品国产色婷小说| 日韩视频一区二区在线观看| 午夜日韩欧美国产| 国产一区二区三区在线臀色熟女 | 在线看a的网站| 女性被躁到高潮视频| 亚洲一码二码三码区别大吗| 女性被躁到高潮视频| 满18在线观看网站| 一边摸一边抽搐一进一出视频| 大型av网站在线播放| 亚洲午夜精品一区,二区,三区| 国产成人精品久久二区二区91| 色老头精品视频在线观看| 天天操日日干夜夜撸| 欧美精品亚洲一区二区| 亚洲精品自拍成人| 变态另类成人亚洲欧美熟女 | 精品久久久久久电影网| 精品亚洲乱码少妇综合久久| 99热网站在线观看| 在线观看人妻少妇| 婷婷丁香在线五月| 亚洲国产毛片av蜜桃av| 欧美成人免费av一区二区三区 | 黑人巨大精品欧美一区二区mp4| 丝袜人妻中文字幕| 国产无遮挡羞羞视频在线观看| tocl精华| 伊人久久大香线蕉亚洲五| aaaaa片日本免费| 丝袜喷水一区| 人妻久久中文字幕网| 80岁老熟妇乱子伦牲交| 搡老乐熟女国产| 一二三四社区在线视频社区8| xxxhd国产人妻xxx| 视频区欧美日本亚洲| 国产亚洲精品一区二区www | 欧美在线一区亚洲| av免费在线观看网站| 这个男人来自地球电影免费观看| 午夜两性在线视频| 美女主播在线视频| 久久精品国产a三级三级三级| 桃花免费在线播放| 国产又色又爽无遮挡免费看| 黄色丝袜av网址大全| 美女高潮喷水抽搐中文字幕| 大片免费播放器 马上看| 99精品欧美一区二区三区四区| 天天躁日日躁夜夜躁夜夜| 国产成人影院久久av| 99久久人妻综合| 成人三级做爰电影| 亚洲少妇的诱惑av| av天堂久久9| 久久久国产精品麻豆| 美女国产高潮福利片在线看| 午夜精品久久久久久毛片777| 日本欧美视频一区|