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

    長(zhǎng)三角地區(qū)景觀生態(tài)風(fēng)險(xiǎn)與人類活動(dòng)強(qiáng)度的響應(yīng)及其時(shí)空關(guān)聯(lián)

    2022-03-09 03:27:32劉艷曉李燦鋒
    關(guān)鍵詞:長(zhǎng)三角時(shí)空尺度

    方 林,方 斌,劉艷曉,蔡 俊,李燦鋒

    長(zhǎng)三角地區(qū)景觀生態(tài)風(fēng)險(xiǎn)與人類活動(dòng)強(qiáng)度的響應(yīng)及其時(shí)空關(guān)聯(lián)

    方 林1,方 斌2※,劉艷曉3,蔡 俊1,李燦鋒4

    (1. 安徽農(nóng)業(yè)大學(xué)經(jīng)濟(jì)管理學(xué)院,合肥 230036;2. 南京師范大學(xué)地理科學(xué)學(xué)院,南京 210023;3. 東北師范大學(xué)地理科學(xué)學(xué)院,長(zhǎng)春 130024;4. 中國(guó)地質(zhì)調(diào)查局昆明自然資源綜合調(diào)查中心,昆明 650100)

    人類活動(dòng)是誘發(fā)區(qū)域生態(tài)風(fēng)險(xiǎn)的重要因素,但生態(tài)風(fēng)險(xiǎn)與人類活動(dòng)在不同尺度下的表征效果不一,響應(yīng)程度不同,從多尺度科學(xué)地研究區(qū)域景觀生態(tài)風(fēng)險(xiǎn)與人類活動(dòng)強(qiáng)度的關(guān)聯(lián)特征對(duì)于協(xié)調(diào)區(qū)域人地關(guān)系地域系統(tǒng)矛盾,推動(dòng)地區(qū)可持續(xù)發(fā)展具有重要意義。該研究以長(zhǎng)三角地區(qū)為研究對(duì)象,基于1990—2020年4期土地利用、夜間燈光、人口空間分布數(shù)據(jù),建立市域、縣域、格網(wǎng)三級(jí)尺度,構(gòu)建景觀生態(tài)風(fēng)險(xiǎn)及人類活動(dòng)強(qiáng)度評(píng)估模型,刻畫不同尺度下二者的時(shí)空響應(yīng)特征,并基于Copula函數(shù)、雙變量空間自相關(guān)及耦合協(xié)調(diào)度模型,揭示景觀生態(tài)風(fēng)險(xiǎn)與人類活動(dòng)強(qiáng)度的時(shí)空關(guān)聯(lián)性。結(jié)果表明:1)1990—2020年間不同尺度下長(zhǎng)三角地區(qū)景觀生態(tài)風(fēng)險(xiǎn)均呈現(xiàn)北高南低的空間分布格局,且高風(fēng)險(xiǎn)地區(qū)持續(xù)減少,低風(fēng)險(xiǎn)區(qū)均呈增加趨勢(shì),其中市域尺度下高風(fēng)險(xiǎn)地區(qū)的減勢(shì)最明顯,而縣域尺度下,低及中低風(fēng)險(xiǎn)區(qū)總體增勢(shì)最為顯著;2)1990—2020年間不同尺度下長(zhǎng)三角地區(qū)人類活動(dòng)強(qiáng)度均呈現(xiàn)東北高西南低的空間分布格局且高強(qiáng)度區(qū)域均呈現(xiàn)顯著上升趨勢(shì),而低強(qiáng)度區(qū)域均明顯下降,其中格網(wǎng)尺度下高強(qiáng)度區(qū)域增幅最大,達(dá)13.42個(gè)百分點(diǎn),市域尺度下低強(qiáng)度區(qū)域的減幅最為明顯,達(dá)9.76個(gè)百分點(diǎn);3)1990—2020年間長(zhǎng)三角地區(qū)人類活動(dòng)強(qiáng)度與景觀生態(tài)風(fēng)險(xiǎn)呈現(xiàn)正向相關(guān)性,但人類活動(dòng)強(qiáng)度對(duì)于景觀生態(tài)風(fēng)險(xiǎn)的影響正在削弱??臻g上高-高(H-H)集聚區(qū)域呈現(xiàn)持續(xù)減少趨勢(shì),而低-低(L-L)呈現(xiàn)先增后減的變化趨勢(shì),低-高(L-H)及高-低(H-L)在2010年后增勢(shì)明顯。研究結(jié)果可為長(zhǎng)三角地區(qū)國(guó)土空間優(yōu)化及生態(tài)安全維護(hù)提供理論依據(jù),為區(qū)域生態(tài)風(fēng)險(xiǎn)防范和生態(tài)修復(fù)提供參考。

    土地利用;景觀;生態(tài);人類活動(dòng);尺度響應(yīng);Copula函數(shù);長(zhǎng)三角地區(qū)

    0 引 言

    自工業(yè)文明誕生以來(lái),人類憑借先進(jìn)生產(chǎn)力得到了極大的物質(zhì)滿足,但膨脹的原材料需求,卻給自然生態(tài)系統(tǒng)帶來(lái)了沉重的負(fù)擔(dān)[1]。已有研究表明,受氣候變化及人類活動(dòng)影響,全球陸地生態(tài)系統(tǒng)出現(xiàn)了不同程度的退化,且這一趨勢(shì)今后或?qū)⑦M(jìn)一步加劇[2]。事實(shí)上,區(qū)域尺度下較短時(shí)限內(nèi)影響并誘發(fā)區(qū)域生態(tài)風(fēng)險(xiǎn)演變的自然因素往往趨于穩(wěn)定,難以發(fā)生大的突變,但人類活動(dòng)變化卻異常激烈[3],其顯性作用于陸地表面間接改變土地利用的類型、格局、強(qiáng)度,并對(duì)區(qū)域景觀格局及生態(tài)過程產(chǎn)生深刻影響,最終誘導(dǎo)區(qū)域景觀生態(tài)風(fēng)險(xiǎn)(Landscape Ecological Risk,LER)發(fā)生演變。全面客觀且精細(xì)地刻畫人類活動(dòng)與LER的時(shí)空關(guān)聯(lián)特征是協(xié)調(diào)人地關(guān)系地域系統(tǒng)矛盾、探索區(qū)域可持續(xù)發(fā)展道路的重要科學(xué)問題[4]。

    LER是指人類活動(dòng)影響下景觀格局與生態(tài)過程相互作用所產(chǎn)生的不利后果[5]。LER評(píng)價(jià)起源于20世紀(jì)90年代,較之于傳統(tǒng)生態(tài)風(fēng)險(xiǎn)評(píng)價(jià),其著重強(qiáng)調(diào)風(fēng)險(xiǎn)空間異質(zhì)性的定量表征,其研究結(jié)果是支撐區(qū)域生態(tài)網(wǎng)絡(luò)建設(shè)和生態(tài)可持續(xù)管理的重要依據(jù)[6]。歷經(jīng)30多年的研究發(fā)展,LER研究已發(fā)展成為了地理學(xué)及生態(tài)學(xué)的研究熱點(diǎn)。目前,基于單一尺度進(jìn)行LER評(píng)價(jià)的研究日益成熟[7],其研究對(duì)象的選擇既包括傳統(tǒng)的流域[8]、行政區(qū)[9]及城市地域單元[10],同時(shí)也包括工礦開采區(qū)[11]、海岸帶[12]、農(nóng)牧交錯(cuò)帶[13]等重點(diǎn)風(fēng)險(xiǎn)控制區(qū)。但尺度變化及尺度差異增加了LER評(píng)價(jià)結(jié)果的不確定性,一定程度上造成了LER研究成果的割裂與脫節(jié),極大地削弱了研究成果向?qū)嵺`的轉(zhuǎn)化與應(yīng)用。近年來(lái),部分學(xué)者愈發(fā)關(guān)注景觀生態(tài)研究過程中的尺度問題,如基于適宜空間粒度[14]選取研究尺度,但卻鮮有研究關(guān)注LER評(píng)價(jià)過程中尺度變化的響應(yīng)特征。

    人類活動(dòng)可以理解為人類為滿足自身生存發(fā)展需要對(duì)自然生態(tài)環(huán)境所進(jìn)行的各種開發(fā)、利用與保護(hù)行為的總稱[15-16]。人類活動(dòng)強(qiáng)度(Human Activity Intensity,HAI)是人類行為對(duì)自然環(huán)境所產(chǎn)生干擾作用程度的客觀表征。定量化及可視化表征區(qū)域HAI是揭示其與生態(tài)環(huán)境變化間密切關(guān)系的基礎(chǔ)。盡管伴隨著3S技術(shù)的快速發(fā)展,土地利用、生態(tài)系統(tǒng)服務(wù)等可視化空間制圖技術(shù)愈發(fā)完備[17],但受限于人類活動(dòng)的多元特征及復(fù)雜趨向,與人類活動(dòng)相關(guān)的空間數(shù)據(jù)制備發(fā)展較為緩慢。高質(zhì)量HAI空間數(shù)據(jù)的缺失極大地限制了研究人員對(duì)人類活動(dòng)與生態(tài)風(fēng)險(xiǎn)間復(fù)雜作用關(guān)系的揭露,這也造成了二者研究的長(zhǎng)期剝離與孤立。因此,長(zhǎng)期以來(lái)相關(guān)研究多將人類活動(dòng)視為輔助信息變量,鮮有關(guān)于HAI-LER關(guān)系的定量分析[18]。而HAI和LER評(píng)估的綜合運(yùn)用,不僅能深化人—地間復(fù)雜互動(dòng)關(guān)系的研究,還能將人類福祉與生態(tài)環(huán)境變化相聯(lián)系,更好地為區(qū)域生態(tài)環(huán)境保護(hù)政策制定提供理論支持。

    長(zhǎng)三角地區(qū)是人活動(dòng)密集、經(jīng)濟(jì)發(fā)展顯著且生態(tài)地位突出的典型區(qū)域之一??焖俪鞘谢率褂騼?nèi)耕地保護(hù)、經(jīng)濟(jì)發(fā)展及生態(tài)維護(hù)之間的矛盾激化,其所產(chǎn)生的次生環(huán)境問題嚴(yán)重遲滯阻礙了長(zhǎng)三角地區(qū)的可持續(xù)發(fā)展。為緩解區(qū)域人—地間矛盾關(guān)系,區(qū)域政府實(shí)行了退耕還林、生態(tài)功能區(qū)劃定等一系列政策法規(guī)。因此,本研究旨在以長(zhǎng)三角地區(qū)為研究對(duì)象,探討區(qū)域LER與HAI的尺度響應(yīng)及時(shí)空關(guān)聯(lián)特征,揭示域內(nèi)人地關(guān)系的演變狀況,以期為地區(qū)國(guó)土空間優(yōu)化與生態(tài)安全維護(hù)以及區(qū)域生態(tài)風(fēng)險(xiǎn)防范和生態(tài)修復(fù)提供參考。

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

    1.1 研究區(qū)概況

    長(zhǎng)三角地區(qū)包括上海、浙江、江蘇、安徽在內(nèi)的3省和1直轄市共41個(gè)城市,是中國(guó)人口集聚和經(jīng)濟(jì)活躍程度最高的地區(qū)之一,據(jù)國(guó)家統(tǒng)計(jì)局資料顯示(http://www.stats.gov.cn),截至2020年年底,域內(nèi)人口總量約為2.35億,占中國(guó)總?cè)丝诘?6.66%;地區(qū)經(jīng)濟(jì)生產(chǎn)總值達(dá)24.47萬(wàn)億元,占全國(guó)GDP比值的24.08%。此外,地處長(zhǎng)江下游平原,地形相對(duì)平坦,水熱條件優(yōu)越,年平均氣溫約為17 ℃。年降水量約1 368 mm,自然資源十分富足。受人類活動(dòng)影響,長(zhǎng)三角地區(qū)生態(tài)環(huán)境受到嚴(yán)重挑戰(zhàn),生態(tài)環(huán)境保護(hù)與經(jīng)濟(jì)社會(huì)發(fā)展間矛盾沖突愈發(fā)明顯。系統(tǒng)梳理不同尺度下LER與區(qū)域HAI的時(shí)空關(guān)系演變規(guī)律,對(duì)于推進(jìn)長(zhǎng)三角生態(tài)綠色一體化發(fā)展示范區(qū)建設(shè),加快構(gòu)建區(qū)域生態(tài)環(huán)境共保聯(lián)治機(jī)制具有重要意義。

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

    以長(zhǎng)三角地區(qū)為研究對(duì)象,基于1990—2020年4期土地利用數(shù)據(jù)、夜間燈光數(shù)據(jù)、人口空間分布數(shù)據(jù),對(duì)長(zhǎng)三角地區(qū)LER與HAI的時(shí)空關(guān)聯(lián)特征展開分析。其中土地利用數(shù)據(jù)來(lái)源于中國(guó)科學(xué)院資源與環(huán)境科學(xué)數(shù)據(jù)中心(http://www.resdc.cn/),該數(shù)據(jù)空間分辨率為30 m′30 m,數(shù)據(jù)精度達(dá)90%;1992—2013年DMSP/OLS穩(wěn)定夜間燈光數(shù)據(jù)的像元灰度值范圍為0~63。2012—2018年NPP/VIIRS影像為月尺度數(shù)據(jù)。2種夜間燈光數(shù)據(jù)均來(lái)源于NOAA網(wǎng)站NGDC數(shù)據(jù)中心(https://www.ngdc. noaa.gov/eog/download.html),空間分辨率為1 km′1 km。人口空間分布數(shù)據(jù)來(lái)源于資源環(huán)境科學(xué)數(shù)據(jù)注冊(cè)與出版系統(tǒng)的中國(guó)人口空間分布公里網(wǎng)格數(shù)據(jù)集(http://www.resdc.cn/),空間分辨率為1 km′1 km。其中部分年份的夜光數(shù)據(jù)及人口空間分布數(shù)據(jù)缺乏,本研究采用相近年份的夜間燈光遙感數(shù)據(jù)及人口空間分布數(shù)據(jù)進(jìn)行替代和線性插值處理。為便于研究運(yùn)算,借助ArcGIS 10.5軟件將所有柵格數(shù)據(jù)重采樣至100 m′100 m,投影統(tǒng)一為Krasovsky_1940_Albers。

    2 研究方法

    2.1 研究尺度劃分

    本研究以長(zhǎng)三角地區(qū)為研究對(duì)象,探討區(qū)域LER與HAI的時(shí)空演變規(guī)律,借助ArcGIS10.5軟件,創(chuàng)建3種不同類型的評(píng)價(jià)單元尺度,即市域尺度、縣域尺度、格網(wǎng)尺度。在格網(wǎng)尺度下,參照已有研究[19],采用10 km×10 km格網(wǎng)尺度既能保證研究的精度需求,且運(yùn)算量相對(duì)較低。其中縣域及市域尺度,參照2020年行政區(qū)劃標(biāo)準(zhǔn)進(jìn)行劃分。分別計(jì)算3種不同尺度下的LER及HAI變化,對(duì)比研究區(qū)不同尺度下的景觀生態(tài)風(fēng)險(xiǎn)及人類活動(dòng)強(qiáng)度的空間分布及其異質(zhì)性,揭示長(zhǎng)三角地區(qū)LER與HAI對(duì)尺度變化的響應(yīng)特征。

    2.2 LER評(píng)估

    景觀格局變化引起的生態(tài)風(fēng)險(xiǎn)受外部因素和內(nèi)部因素的影響,其中外部因素是指生態(tài)系統(tǒng)在受到脅迫時(shí)的干擾程度,用景觀干擾度進(jìn)行表征;內(nèi)部因素是指生態(tài)系統(tǒng)抵御壓力時(shí)的脆弱程度,借助景觀脆弱度進(jìn)行表征。研究基于Fragstats4.2軟件,基于研究區(qū)實(shí)際情況及已有研究[20],圍繞景觀干擾度指數(shù)、景觀脆弱度指數(shù)構(gòu)建長(zhǎng)三角LER評(píng)價(jià)模型,具體計(jì)算式為

    式中D第景觀類型的景觀干擾指數(shù);、、分別代表相應(yīng)的各景觀指數(shù)的權(quán)重,++=1,參照尉芳等[21]研究成果,、、分別定為0.5、0.3、0.2;F表示景觀破碎度,S表示景觀分離度,F(xiàn)D表示景觀分形維數(shù);N表示景觀類型的斑塊數(shù)量,A表示景觀類型的面積,為景觀類型的總面積(hm2),P為景觀類型的周長(zhǎng)(m)。R表示第類景觀類型的生態(tài)損失度指數(shù);V第景觀類型的景觀脆弱度指數(shù),其中景觀脆弱度指數(shù)參照已有研究成果[22],采取專家打分法,對(duì)不同景觀類型進(jìn)行賦值,其中未利用地、水域、耕地、草地、林地、建設(shè)用地分別賦值為6、5、4、3、2、1,經(jīng)歸一化處理得到各景觀類型自身的脆弱度指數(shù):未利用地0.29、水域0.24、耕地0.19、草地0.14、林地0.10、建設(shè)用地0.05;ERI為第個(gè)風(fēng)險(xiǎn)小區(qū)生態(tài)風(fēng)險(xiǎn)指數(shù),A為第個(gè)風(fēng)險(xiǎn)小區(qū)第類景觀的面積,A為第個(gè)風(fēng)險(xiǎn)小區(qū)的面積。

    2.3 HAI表征

    HAI是度量人類改造自然力度的有效指標(biāo)[18],現(xiàn)有關(guān)于HAI的表征方法可概括為直接或間接空間化兩類空間表征的方法。其中土地利用變化是人類經(jīng)濟(jì)社會(huì)活動(dòng)的綜合反映,同時(shí)可看作人類活動(dòng)作用于陸地表面的顯性后果,是用以表征HAI的理想媒介。但單純以土地利用變化等替評(píng)價(jià)HAI的方法,既難以反映多元、復(fù)雜和綜合的人類活動(dòng),同時(shí)也會(huì)與同樣以土地利用數(shù)據(jù)為基礎(chǔ)的景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)產(chǎn)生一定程度的內(nèi)生效應(yīng)。為此,本研究在綜合前人研究的基礎(chǔ)上,引入夜間燈光及人口空間分布數(shù)據(jù),并借鑒榮益等[23]的研究成果,綜合Lohani清單法、Leopold矩陣法及Delphi法3種方法,并通過求取平均值的方法減少系數(shù)誤差,最后針對(duì)不同土地利用類型賦予不同的人類活動(dòng)影響強(qiáng)度系數(shù),進(jìn)而構(gòu)建長(zhǎng)三角地區(qū)HAI綜合評(píng)價(jià)模型,其表達(dá)式如下:

    式中、、分別表示歸一化后的夜間燈光指數(shù)NTL、人口分布數(shù)據(jù)PD、以及土地利用數(shù)據(jù)LU的指標(biāo)權(quán)重。參照陳泓瑾等[4]研究成果,、、分別取0.3、0.3、0.4。

    2.4 LER和HAI分級(jí)標(biāo)準(zhǔn)

    有鑒于研究涉及的時(shí)間跨度較大,尺度較多,為便于直觀展示1990—2020年間不同尺度下LER及HAI時(shí)空演化特征,參照已有研究成果[24],采用Jenks最佳自然斷裂法,將LER與HAI由低向高劃分為5個(gè)等級(jí)。

    2.5 時(shí)空關(guān)聯(lián)分析

    2.5.1 Copula函數(shù)建模

    Copula理論是由Sklar于1959年首次提出,旨在揭示不同變量間線性及非線性的復(fù)雜關(guān)系,傳統(tǒng)皮爾遜相關(guān)分析僅能夠度量線性變化下保持穩(wěn)定的指標(biāo),其數(shù)據(jù)要求嚴(yán)格的正態(tài)分布,難以捕捉變量間的非線性關(guān)系[25]。Copula函數(shù)可有效避免以上問題,現(xiàn)已在經(jīng)濟(jì)學(xué)、生態(tài)學(xué)、氣象學(xué)等復(fù)雜科學(xué)領(lǐng)域得到廣泛的應(yīng)用[26]。常見的Copula參數(shù)族包括橢球Copula(二元Gaussian Copula、二元t-Copula)、阿基米德Copula(二元Gumbel Copula、Clayton Copula、Frank Copula)。不同Copula函數(shù)的選擇直接關(guān)乎不同變量間關(guān)聯(lián)測(cè)度的結(jié)果,因此Copula建模的一個(gè)核心問題便是Copula函數(shù)的選擇問題。本研究基于Copula理論,采用非參數(shù)估計(jì),借助核分布估計(jì)法調(diào)用合適的Copula函數(shù),構(gòu)建LER-HAI互動(dòng)關(guān)系分析模型,采用平方歐氏距離2(其值越小,模型擬合度越高)對(duì)模型進(jìn)行驗(yàn)證,完成Copula函數(shù)的尋優(yōu)計(jì)算,分析LER-HAI的時(shí)序演變規(guī)律。研究選用Kendall以及Spearman秩相關(guān)系數(shù)對(duì)LER與HAI的互動(dòng)關(guān)系進(jìn)行測(cè)度,二者在數(shù)據(jù)服從正態(tài)分布時(shí)與Spearman相關(guān)測(cè)度結(jié)果一致,相較之下,Kendall和Spearman秩相關(guān)系數(shù)能夠同時(shí)測(cè)度變量間線性及非線性的復(fù)雜作用關(guān)系,其相關(guān)性測(cè)度結(jié)果更加可靠[27]。建模過程通過MATLAB 19b軟件完成,不同函數(shù)的表達(dá)式見文獻(xiàn)[28]。

    2.5.2 雙變量局部空間自相關(guān)

    為進(jìn)一步刻畫長(zhǎng)三角地區(qū)HAI與LER在空間上的互動(dòng)關(guān)系變化,本研究借助雙變量局部空間自相關(guān)分析方法,探究LER與HAI的空間關(guān)聯(lián)模式,并通過GeoDa 1.18軟件繪制LISA聚類圖。模型的具體表達(dá)式參照文獻(xiàn)[29]。

    2.5.3 耦合測(cè)度

    借助物理學(xué)中的耦合模型,構(gòu)建LER與HAI的耦合協(xié)調(diào)度模型,定量分析長(zhǎng)三角地區(qū)LER與HAI耦合過程與演進(jìn)趨勢(shì),計(jì)算式為

    式中表示長(zhǎng)三角地區(qū)LER與HAI的耦合系數(shù),反映二者間相互影響的程度;表示二者間的綜合調(diào)和指數(shù)(本研究參照前人[30]研究成果,待定系數(shù)取0.5);表示LER與HAI的耦合協(xié)調(diào)程度。

    3 結(jié)果與分析

    3.1 長(zhǎng)三角地區(qū)LER時(shí)空特征分析

    1990—2020年間不同尺度下長(zhǎng)三角地區(qū)景觀生態(tài)時(shí)空變化特征,如圖1、圖2所示。整體來(lái)看,1990—2020年間不同尺度下長(zhǎng)三角LER高風(fēng)險(xiǎn)面積占比均呈現(xiàn)下降趨勢(shì),低風(fēng)險(xiǎn)地區(qū)均呈現(xiàn)上升趨勢(shì),中風(fēng)險(xiǎn)地區(qū)擴(kuò)大趨勢(shì)明顯;3種尺度下的LER空間分布存在明顯的共性規(guī)律,即北高南低的空間分布格局,高風(fēng)險(xiǎn)區(qū)主要分布在長(zhǎng)江以北地區(qū),低風(fēng)險(xiǎn)區(qū)主要分布于皖西、皖南、浙西南地區(qū)。不同尺度下LER呈現(xiàn)不同的尺度響應(yīng)特征,其中市域尺度下,高風(fēng)險(xiǎn)區(qū)較其他尺度減少最為明顯,1990 —2020年間下降了19.51個(gè)百分點(diǎn),低風(fēng)險(xiǎn)區(qū)增長(zhǎng)較其他尺度亦最為明顯,但中低風(fēng)險(xiǎn)區(qū)域呈現(xiàn)減少趨勢(shì),因此低及中低風(fēng)險(xiǎn)區(qū)總體增長(zhǎng)趨勢(shì)較低,僅為4.88%??h域尺度層面,高風(fēng)險(xiǎn)及中高風(fēng)險(xiǎn)區(qū)域面積占比均呈現(xiàn)減少趨勢(shì),30 a間分別減少17.60和3.60個(gè)百分點(diǎn),低及中低風(fēng)險(xiǎn)區(qū)總體增幅較其他尺度最為明顯,1990 —2020年間共增加了16.29個(gè)百分點(diǎn)。格網(wǎng)尺度層面,高風(fēng)險(xiǎn)區(qū)占比較低,減勢(shì)平緩,中高風(fēng)險(xiǎn)區(qū)占比1990—2020年間減少了8.79個(gè)百分點(diǎn),中風(fēng)險(xiǎn)區(qū)域整體穩(wěn)定在31.00%左右,低及中低風(fēng)險(xiǎn)區(qū)總體占比較大,但總體增幅較低,30 a間總體增加了10.39個(gè)百分點(diǎn)。

    3.2 長(zhǎng)三角地區(qū)HAI時(shí)空特征分析

    1990—2020年間不同尺度下長(zhǎng)三角地區(qū)HAI的時(shí)空變化特征,如圖3所示。3種尺度下,區(qū)域HAI時(shí)空變化整體保持一致,其中高強(qiáng)度及中高強(qiáng)度區(qū)域面積占比持續(xù)上升,中強(qiáng)度Ⅲ區(qū)域面積占比相對(duì)穩(wěn)定,中低強(qiáng)度及低強(qiáng)度區(qū)域面積占比均呈現(xiàn)下降趨勢(shì),這表明1990—2020年間長(zhǎng)三角區(qū)域人類活動(dòng)呈現(xiàn)低強(qiáng)度向高強(qiáng)度遞進(jìn)的時(shí)序變化特征??臻g上長(zhǎng)三角地區(qū)HAI呈現(xiàn)東北高西南低的空間分布格局,其中高強(qiáng)度地區(qū)主要圍繞以上海為核心的長(zhǎng)三角城市群附近,而低強(qiáng)度區(qū)則主要分布在皖西、皖南及浙西地區(qū)的丘陵山地附近。具體到不同的尺度條件下長(zhǎng)三角地區(qū)HAI的時(shí)空特征而言,市域尺度向格網(wǎng)尺度的遞減變化過程中,高強(qiáng)度區(qū)域占比不斷擴(kuò)大,在格網(wǎng)尺度下,高強(qiáng)度區(qū)占比由1990年的0.91%增長(zhǎng)到2020年的14.33%,增幅最為明顯,達(dá)13.42個(gè)百分點(diǎn),縣域尺度下,中高強(qiáng)度區(qū)的增幅變化最為明顯,30 a間增幅達(dá)24.10個(gè)百分點(diǎn),中等強(qiáng)度區(qū)區(qū)域在格網(wǎng)尺度下減幅最為明顯,1990—2020年間減少了14.54個(gè)百分點(diǎn),中低強(qiáng)度區(qū)在縣域尺度下變化最為明顯,達(dá)16.29個(gè)百分點(diǎn),市域尺度下低強(qiáng)度區(qū)域的減幅最為明顯,30 a間減少達(dá)9.76個(gè)百分點(diǎn)。

    圖1 1990—2020年不同尺度下景觀生態(tài)風(fēng)險(xiǎn)(LER)等級(jí)占比變化特征

    圖2 1990—2020年不同尺度下LER時(shí)空變化特征

    3.3 長(zhǎng)三角地區(qū)LER與HAI時(shí)空關(guān)聯(lián)分析

    3.3.1 特征尺度選擇

    有鑒于多尺度所帶來(lái)的冗雜性,本研究參照黃木易等[31]的研究方法,基于半變異系數(shù),選取具有最佳空間響應(yīng)特征的尺度,不同尺度下LER和HAI的半變異系數(shù)如表1所示,其中半變異系數(shù)越大,表明其空間差異的信息量越豐富,反之信息量越少。綜合考慮LER和HAI的半變異系數(shù)值及樣本量,本研究將格網(wǎng)尺度作為L(zhǎng)ER及HAI時(shí)空關(guān)聯(lián)研究的特征尺度。

    3.3.2 LER與HAI時(shí)間關(guān)聯(lián)特征

    Copula函數(shù)相關(guān)性測(cè)度結(jié)果如圖4所示。選用Clayton Copula函數(shù)時(shí),平方歐氏距離2相對(duì)較小,這說明Clayton Copula函數(shù)下的模型擬合優(yōu)度最好,因此本研究主要圍繞Clayton Copula函數(shù)運(yùn)算結(jié)果對(duì)1990—2020年間HAI和LER相關(guān)性進(jìn)行分析。綜合Kendall以及Spearman的相關(guān)性測(cè)度結(jié)果可知,長(zhǎng)三角地區(qū)LER與HAI呈現(xiàn)明顯的正向相關(guān)性,即伴隨著HAI的增加,區(qū)域LER將不斷增長(zhǎng)。但值得注意的是1990—2020年間Kendall以及Spearman相關(guān)系數(shù)均呈現(xiàn)由高向低的下降趨勢(shì),其中Kendall由1990年的0.458 5下降到2020年的0.309 8;而Spearman由1990年的0.634 8下降到2020年的0.447 3。綜合以上結(jié)果,本研究認(rèn)為1990—2020年間HAI與LER始終呈現(xiàn)正向相關(guān)關(guān)系,且相關(guān)性減弱趨勢(shì)明顯,但在未來(lái)較長(zhǎng)一段時(shí)間仍將呈現(xiàn)正向相關(guān)性。這意味著長(zhǎng)三角地區(qū)HAI的增加對(duì)于區(qū)域LER的影響正在削弱,即二者的沖突性降低,協(xié)調(diào)性增加。

    圖3 1990—2020年不同尺度下人類活動(dòng)(HAI)強(qiáng)度時(shí)空變化特征

    表1 長(zhǎng)三角地區(qū)LER與HAI半變異系數(shù)

    3.3.3 LER與HAI空間關(guān)聯(lián)特征

    雙變量局部空間自相關(guān)分析結(jié)果如圖5所示,H-H區(qū)表示高HAI與高LER集聚;L-L表示低HAI與低LER集聚;H-L表示高HAI與低LER集聚;L-H表示低HAI與高LER集聚。H-H區(qū)域早期(1990年)主要分布在安徽中部、江蘇東南部以及上海周邊地區(qū),但2000—2020年間,H-H區(qū)域呈現(xiàn)明顯減少趨勢(shì),并向江蘇中部、安徽中北部轉(zhuǎn)移。L-L地區(qū)主要分布于皖西大別山、皖南丘陵、浙西南地區(qū),1990—2000年間L-L區(qū)域面積占比穩(wěn)定,2010—2020年間L-L占比呈現(xiàn)一定的下降趨勢(shì)。L-H地區(qū)與H-L地區(qū)則呈零星狀的空間分布特征,整體面積占比較少,L-H與H-L地區(qū)占比均呈現(xiàn)增長(zhǎng)趨勢(shì),其中L-H主要增長(zhǎng)區(qū)域?yàn)榘不?、江蘇中部及皖江城市帶地區(qū);而H-L地區(qū)的增長(zhǎng)區(qū)域則分布于以上海為核心的城市群周圍。

    3.3.4 LER與HAI耦合度分析

    為進(jìn)一步探明LER與HAI之間互動(dòng)關(guān)系,本研究對(duì)二者間耦合協(xié)調(diào)關(guān)系的演進(jìn)趨勢(shì)進(jìn)行分析,結(jié)果如圖6所示。長(zhǎng)三角地區(qū)1990—2020年間的耦合協(xié)調(diào)度“質(zhì)心”值為0.305 8。按照此前研究對(duì)于耦合協(xié)調(diào)度的范圍定義[30],0.2<≤0.4,表示輕度失調(diào)。這表明LER與HAI處于沖突階段,且輕度失調(diào)。而1990—2020年間LER與HAI的耦合協(xié)調(diào)度整體上呈現(xiàn)增長(zhǎng)趨勢(shì),由1990年的0.303 1上升至2020年的0.311 2。這表明LER與HAI間關(guān)系逐漸由沖突向協(xié)調(diào)轉(zhuǎn)變。耦合度分析結(jié)果印證了此前關(guān)于長(zhǎng)三角地區(qū)LER與HAI時(shí)空關(guān)聯(lián)特征的結(jié)論。

    注:Kendall τ,Spearman ρ表示HAI與LER的相關(guān)系數(shù);d2表示不同Copula函數(shù)的平方歐氏距離。

    注:H-H區(qū)表示高高集聚;L-L表示低低集聚;H-L表示高低集聚;L-H表示低高集聚。

    圖6 1990—2020年長(zhǎng)三角地區(qū)耦合協(xié)調(diào)度變化

    4 討 論

    4.1 LER與HAI多尺度時(shí)空變化特征

    1990—2020年間長(zhǎng)三角地區(qū)不同尺度下LER等級(jí)較高地區(qū)均呈現(xiàn)持續(xù)減少趨勢(shì)且減勢(shì)在近年來(lái)有著明顯的加劇,而低風(fēng)險(xiǎn)等級(jí)地區(qū)呈現(xiàn)明顯的上升趨勢(shì)且近年來(lái)增勢(shì)愈發(fā)明顯。這與Ran等[32]此前的研究結(jié)論一致。呈現(xiàn)如此變化趨勢(shì)得益于2000年來(lái)中國(guó)政府針對(duì)長(zhǎng)三角地區(qū)生態(tài)環(huán)境治理所提出的諸如:“退耕還林還草(2001年)”“大別山水土保持生態(tài)功能區(qū)(2010年)”“新安江流域生態(tài)補(bǔ)償試點(diǎn)(2011年)”等一系列支持性政策有關(guān)。長(zhǎng)三角地區(qū)LER呈現(xiàn)北高南低的空間分布格局,區(qū)域東北部大面積的高風(fēng)險(xiǎn)地區(qū)減少及西南部大面積低風(fēng)險(xiǎn)區(qū)域的增長(zhǎng)預(yù)示著長(zhǎng)三角地區(qū)生態(tài)環(huán)境整體向好的趨勢(shì)特征,同時(shí)也表明區(qū)域政府的生態(tài)治理保護(hù)措施取得了顯著成效。

    1990—2020年間,不同尺度下長(zhǎng)三角地區(qū)HAI均呈現(xiàn)明顯的上升趨勢(shì),特別是高強(qiáng)度地區(qū)面積占比持續(xù)增長(zhǎng),低強(qiáng)度地區(qū)面積占比呈現(xiàn)明顯的降低趨勢(shì),這與長(zhǎng)三角地區(qū)良好的自然稟賦、得天獨(dú)厚的區(qū)位優(yōu)勢(shì)、發(fā)達(dá)的經(jīng)濟(jì)基礎(chǔ)有著密切的關(guān)系,三者顯著影響并重塑著區(qū)域人口空間分布格局,使得長(zhǎng)三角地區(qū)成為全國(guó)性的人口集聚中心。人口集聚與HAI攀升是區(qū)域社會(huì)經(jīng)濟(jì)發(fā)展的必然結(jié)果,該趨勢(shì)或在未來(lái)進(jìn)一步加劇。預(yù)防長(zhǎng)三角東北部HAI較高地區(qū)的生態(tài)環(huán)境質(zhì)量下降及生態(tài)風(fēng)險(xiǎn)升高,是區(qū)域生態(tài)治理與保護(hù)的重點(diǎn)。

    尺度對(duì)于揭示LER與HAI的時(shí)空演化特征十分關(guān)鍵,由于生態(tài)系統(tǒng)及人類活動(dòng)的高度復(fù)雜性,不同尺度下LER與HAI的時(shí)空特征存在一定的差異。但研究發(fā)現(xiàn),不同尺度下,無(wú)論是長(zhǎng)三角地區(qū)LER或是區(qū)域HAI的時(shí)空特征均存在明顯的共性,誠(chéng)然亦存在一定的差異特征,這與張弛等[33]此前的研究結(jié)論較為一致。無(wú)論是LER還是HAI,二者的時(shí)空特征均會(huì)受到取樣尺度變化的影響,但這種影響并不是顛覆性的。由于較大尺度的空間幅度及粒度均高于小尺度,這就使得較大尺度下的研究分辨率偏低,局部范圍的景觀或人類活動(dòng)信息易被模糊與忽略,但這并不意味著較大尺度下研究結(jié)論的失真,相反不同尺度有其對(duì)應(yīng)的實(shí)踐目標(biāo)。單一尺度下的LER或是HAI研究會(huì)忽略甚至扭曲其時(shí)空演變的某種規(guī)律,淡化或削弱研究結(jié)果的適用性或?qū)嵺`性,區(qū)域政府在開展當(dāng)?shù)厣鷳B(tài)保護(hù)與治理實(shí)踐時(shí),應(yīng)依據(jù)不同的治理范圍及治理目標(biāo)選擇合適的尺度研究結(jié)果。

    4.2 LER與HAI的時(shí)空關(guān)聯(lián)

    研究發(fā)現(xiàn),1990—2020年間HAI與LER呈現(xiàn)較強(qiáng)的正向相關(guān)性,即HAI的增加會(huì)激發(fā)區(qū)域LER的產(chǎn)生,這與前人[34-35]的研究結(jié)論一致。但與之不同的是,長(zhǎng)三角地區(qū)人類活動(dòng)與生態(tài)風(fēng)險(xiǎn)間的相關(guān)系數(shù)呈現(xiàn)持續(xù)衰減態(tài)勢(shì),這表明區(qū)域HAI的增加對(duì)于LER的影響正在削弱,二者關(guān)系正經(jīng)歷沖突向協(xié)調(diào)的轉(zhuǎn)變。耦合度分析結(jié)果表明長(zhǎng)三角地區(qū)LER與HAI處于輕度失調(diào)水平,但區(qū)域LER與HAI的耦合協(xié)調(diào)度卻呈現(xiàn)明顯的上升趨勢(shì),這也印證了“長(zhǎng)三角地區(qū)人類活動(dòng)與生態(tài)風(fēng)險(xiǎn)間關(guān)系由沖突向協(xié)調(diào)轉(zhuǎn)變”這一結(jié)論。

    傳統(tǒng)觀點(diǎn)認(rèn)為人類活動(dòng)加劇會(huì)給自然生態(tài)系統(tǒng)帶來(lái)沉重的負(fù)擔(dān),激化區(qū)域人地關(guān)系地域系統(tǒng)矛盾。但近年來(lái)部分學(xué)者[36-37]研究認(rèn)為,人類社會(huì)發(fā)展與自然生態(tài)環(huán)境的關(guān)系會(huì)呈現(xiàn)先惡化后改善的倒U形發(fā)展趨勢(shì),即經(jīng)典的“環(huán)境庫(kù)茲涅茨曲線”假說。受經(jīng)濟(jì)社會(huì)發(fā)展影響,人口會(huì)向發(fā)達(dá)的城市地區(qū)集聚,高強(qiáng)度的人口負(fù)荷,會(huì)對(duì)區(qū)域自然生態(tài)產(chǎn)生一定的影響,其影響程度與區(qū)域生態(tài)承載力有關(guān),因此長(zhǎng)三角地區(qū)HAI與LER呈現(xiàn)正向相關(guān)關(guān)系,但伴隨著經(jīng)濟(jì)發(fā)展至一定水平,人們生態(tài)保護(hù)意識(shí)的提高加之區(qū)域政府對(duì)于生態(tài)保護(hù)與治理力度的提升,HAI的增加對(duì)于區(qū)域LER的影響會(huì)逐步降低,甚至出現(xiàn)“拐點(diǎn)”。中國(guó)政府近20多年的生態(tài)保護(hù)實(shí)踐也表明:通過政策干預(yù)手段,能夠有效降低區(qū)域生態(tài)風(fēng)險(xiǎn),及時(shí)修復(fù)受創(chuàng)的自然生態(tài)系統(tǒng)。

    盡管HAI對(duì)LER的激化效應(yīng)明顯淡化,但“沖突”向“協(xié)調(diào)”關(guān)系轉(zhuǎn)換的過程是漫長(zhǎng)的,未來(lái)區(qū)域政府仍需加強(qiáng)域內(nèi)的生態(tài)環(huán)境保護(hù)。當(dāng)?shù)卣山Y(jié)合HAI與LER的空間關(guān)聯(lián)特征,因地制宜分區(qū)治理。其中,高HAI和高LER地區(qū)主要分布于皖中、蘇東南及上海周邊地區(qū),但伴隨著長(zhǎng)三角核心城市群產(chǎn)業(yè)升級(jí),經(jīng)濟(jì)轉(zhuǎn)型,外來(lái)務(wù)工人員呈現(xiàn)一定的回流趨勢(shì)[38-39],區(qū)域發(fā)展對(duì)于自然生態(tài)系統(tǒng)的擾動(dòng)逐漸減少,因此在近年來(lái)高強(qiáng)度-高風(fēng)險(xiǎn)地區(qū)整體減少,且逐漸向江蘇中部及安徽中北部擴(kuò)散,而這些地區(qū)人口稠密,經(jīng)濟(jì)起步晚,發(fā)展迅速,城市擴(kuò)張明顯,對(duì)生態(tài)的擾動(dòng)程度明顯增強(qiáng)。因此2010年后江蘇中部及安徽中北部地區(qū)低HAI-高LER區(qū)域明顯增多,而上海等長(zhǎng)三角核心城市群附近地區(qū)的高HAI-低LER區(qū)域也呈現(xiàn)了增長(zhǎng)趨勢(shì)。此類現(xiàn)象的產(chǎn)生,本質(zhì)是長(zhǎng)三角內(nèi)部不同地區(qū)發(fā)展階段的差異所造成的[40]。未來(lái)需警惕長(zhǎng)三角腹地經(jīng)濟(jì)騰飛過程中所產(chǎn)生的次生環(huán)境問題,加快長(zhǎng)三角核心區(qū)域的產(chǎn)業(yè)升級(jí)及經(jīng)濟(jì)轉(zhuǎn)型節(jié)奏,擴(kuò)大高HAI-低LER區(qū)域面積,協(xié)調(diào)二者間的關(guān)系;同時(shí)著重加強(qiáng)皖西、皖南、浙西南等低HAI-低LER地區(qū)的生態(tài)保護(hù),防止區(qū)域經(jīng)濟(jì)發(fā)展侵吞原有的自然生態(tài)空間。

    Copula函數(shù)的應(yīng)用能夠很好的彌補(bǔ)現(xiàn)有研究對(duì)變量間可能存在的非線性關(guān)系關(guān)注不足的問題,但受限于模型及連續(xù)數(shù)據(jù)缺失的影響,研究對(duì)于HAI與LER線性及非線性關(guān)系轉(zhuǎn)換的閾值揭示不足,后續(xù)需要探究更為適當(dāng)方法揭示二者關(guān)系轉(zhuǎn)換的閾值,并對(duì)二者的關(guān)聯(lián)機(jī)制進(jìn)行剖析。

    5 結(jié) 論

    本研究基于多期土地利用、夜間燈光及人口空間分布數(shù)據(jù),分別探討了景觀生態(tài)風(fēng)險(xiǎn)(LER)與人類活動(dòng)強(qiáng)度(HAI)的尺度響應(yīng)特征,并基于Copula函數(shù)、雙變量空間自相關(guān)及耦合協(xié)調(diào)度模型,對(duì)二者的空間關(guān)聯(lián)特征進(jìn)行細(xì)致分析,研究結(jié)論如下:

    1)1990—2020年間不同尺度下長(zhǎng)三角地區(qū)LER均呈現(xiàn)北高南低的空間分布格局,且高風(fēng)險(xiǎn)地區(qū)呈現(xiàn)減少趨勢(shì),低風(fēng)險(xiǎn)地區(qū)占比明顯增多;市域尺度下,高風(fēng)險(xiǎn)地區(qū)減勢(shì)最明顯,30 a減少了19.51個(gè)百分點(diǎn);縣域尺度下,低及中低風(fēng)險(xiǎn)區(qū)總體增勢(shì)最為明顯,30 a間增加了16.29個(gè)百分點(diǎn)。

    2)不同尺度下長(zhǎng)三角地區(qū)HAI均呈現(xiàn)東北高西南低的空間分布格局,高強(qiáng)度地區(qū)主要分布在以上海為核心的長(zhǎng)三角城市群附近,而低強(qiáng)度區(qū)則主要分布在皖西、皖南及浙西南地區(qū);1990—2020年間不同尺度下高強(qiáng)度區(qū)域均呈現(xiàn)顯著上升趨勢(shì),而低強(qiáng)度區(qū)域明顯下降,其中格網(wǎng)尺度下高強(qiáng)度區(qū)域增幅最大,達(dá)13.42個(gè)百分點(diǎn),市域尺度下低強(qiáng)度區(qū)域的減幅最為明顯,30 a間減少了9.76個(gè)百分點(diǎn)。

    3)格網(wǎng)尺度為L(zhǎng)ER及HAI時(shí)空關(guān)聯(lián)研究的特征尺度。1990—2020年間長(zhǎng)三角地區(qū)HAI與LER呈現(xiàn)正向相關(guān)性,但相關(guān)性呈現(xiàn)減弱趨勢(shì),而二者耦合協(xié)調(diào)度呈現(xiàn)上升趨勢(shì),從1990年的0.303 1上升到2020年的0.311 2,這說明區(qū)域人類活動(dòng)強(qiáng)度對(duì)景觀生態(tài)風(fēng)險(xiǎn)的影響在減弱,兩者之間的關(guān)系正逐步從沖突轉(zhuǎn)向協(xié)調(diào)。空間上,H-H集聚區(qū)域呈現(xiàn)持續(xù)減少趨勢(shì),而L-L呈現(xiàn)先增后減,整體減少的變化趨勢(shì),L-H及H-L在2010年后增勢(shì)明顯,未來(lái)應(yīng)結(jié)合區(qū)域LER與HAI空間關(guān)聯(lián)特征,對(duì)區(qū)域生態(tài)環(huán)境進(jìn)行分區(qū)治理管控。

    [1] 段群滔,羅立輝. 人類活動(dòng)強(qiáng)度空間化方法綜述與展望:以青藏高原為例[J]. 冰川凍土,2021,43(5):1582-1593.

    Duan Quntao, Luo Lihui. Summary and prospect of spatialization method of human activity intensity: Taking the Qinghai-Tibet Plateau as an example[J]. Journal of Glaciology and Geocryology, 2021, 43(5): 1582-1593. (in Chinese with English abstract)

    [2] Burkhard B, Kroll F, Nedkov S, et. al. Mapping ecosystem service supply, demand and budgets[J]. Ecological Indicators, 2012, 21: 17-29.

    [3] 劉坤,于賜剛,張藝凡,等. 青藏高原自然保護(hù)區(qū)人類活動(dòng)及其影響研究現(xiàn)狀與熱點(diǎn)[J]. 應(yīng)用與環(huán)境生物學(xué)報(bào),2022,28(2):508-516.

    Liu Kun, Yu Cigang, Zhang Yifan, et al. Research status and current hotspots on the human impact on natural reserves in the Qinghai-Tibetan Plateau[J]. Chinese Journal of Applied and Environmental, 2022, 28(2): 508-516. (in Chinese with English abstract)

    [4] 陳泓瑾,劉琳,張正勇,等. 天山北坡人類活動(dòng)強(qiáng)度與地表溫度的時(shí)空關(guān)聯(lián)性[J]. 地理學(xué)報(bào),2022,77(5):1244-1259.

    Chen Hongjin, Liu Lin, Zhang Zhengyong, et al. Spatiotemporal correlation between human activity intensity and surface temperature on the north slope of Tianshan Mountains[J]. Acta Geographica Sinica, 2022, 77(5): 1244-1259. (in Chinese with English abstract)

    [5] 彭建,黨威雄,劉焱序,等. 景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)研究進(jìn)展與展望[J]. 地理學(xué)報(bào),2015,70(4):664-677.

    Peng Jian, Dang Weixiong, Liu Yanxu, et al. Review on landscape ecological risk assessment[J]. Acta Geographica Sinica, 2015, 70(4): 664-677. (in Chinese with English abstract)

    [6] 劉希朝,李效順,蔣冬梅. 基于土地利用變化的黃河流域景觀格局及生態(tài)風(fēng)險(xiǎn)評(píng)估[J]. 農(nóng)業(yè)工程學(xué)報(bào),2021,37(4):265-274.

    Liu Xizhao, Li Xiaoshun, Jiang Dongmei. Landscape pattern identification and ecological risk assessment using land-use change in the Yellow River Basin[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(4): 265-274. (in Chinese with English abstract)

    [7] 馬勝,梁小英,劉迪,等. 生態(tài)脆弱區(qū)多尺度景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià):以陜西省米脂縣高渠鄉(xiāng)為例[J]. 生態(tài)學(xué)雜志,2018,37(10):3171-3178.

    Ma Sheng, Liang Xiaoying, Liu Di, et al. Multi-scale landscape ecological risk assessment in ecologically fragile regions: A case study in Gaoqu Town in Mizhi County, Shaanxi Province[J]. Chinese Journal of Ecology 2018, 37(10): 3171-3178. (in Chinese with English abstract)

    [8] 楊伶,鄧敏,王金龍,等. 近40年來(lái)洞庭湖流域土地利用及生態(tài)風(fēng)險(xiǎn)時(shí)空演變分析[J]. 生態(tài)學(xué)報(bào),2021,41(10):3929-3939.

    Yang Lin, Deng Min, Wang Jinlong, et al. Spatiotemporal evolution of landuseand ecological risk in Dongting Lake Basin during 1980-2018[J]. Acta Ecologica Sinica, 2021, 41(10): 3929-3939. (in Chinese with English abstract)

    [9] Zeng C, He J, He Q, et al. Assessment of land use pattern and landscape ecological risk in the Chengdu-Chongqing Economic Circle, Southwestern China[J]. Land, 2022, 11(5): 659-676.

    [10] Hou Y, Ding W, Liu C, et al. Influences of impervious surfaces on ecological risks and controlling strategies in rapidly urbanizing regions[J]. Science of the Total Environment, 2022, 825: 153823.

    [11] 李少玲,謝苗苗,李漢廷,等. 資源型城市景觀生態(tài)風(fēng)險(xiǎn)的時(shí)空分異:以烏海市為例[J]. 地學(xué)前緣,2021,28(4):100-109.

    Li Shaoling, Xie Miaomiao, Li Hanting, et al. Spatiotemporaldynamics of landscape ecological risk in resource-based cities: A case study of Wuhai[J]. Earth Science Frontiers, 2021, 28(4): 100-109. (in Chinese with English abstract)

    [12] 陳心怡,謝跟蹤,張金萍. 海口市海岸帶近30年土地利用變化的景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)[J]. 生態(tài)學(xué)報(bào),2021,41(3):975-986.

    Chen Xinyi, Xie Gengzong, Zhang Jinping. Landscape ecological risk assessment of land use changes in the coastal area of Haikou City in the past 30 years[J]. Acta Ecologica Sinica, 2021, 41(3): 975-986. (in Chinese with English abstract)

    [13] 徐蘭,羅維,周寶同. 基于土地利用變化的農(nóng)牧交錯(cuò)帶典型流域生態(tài)風(fēng)險(xiǎn)評(píng)價(jià):以洋河為例[J]. 自然資源學(xué)報(bào),2015,30(4):580-590.

    Xu Lan, Luo Wei, Zhou Baotong. Landscape ecological risk assessment of farming-pastoral ecozone based on land use change: A case study of the Yanghe Watershed, China[J]. Journal of Natural Resources, 2015, 30(4): 580-590. (in Chinese with English abstract)

    [14] 楊馗,信桂新,蔣好雨,等基于最佳尺度的景觀生態(tài)風(fēng)險(xiǎn)時(shí)空變化研究:以重慶市江津區(qū)為例[J]. 生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2021,37(5):576-586.

    Yang Kui, Xin Guixin, Jiang Haoyu, et al. Study on spatiotemporal changes of landscape ecological risk based on the optimal spatial scale: A case study of Jiangjin District, Chongqing City[J]. Journal of Ecology and Rural Environment, 2021, 37(5): 576-586. (in Chinese with English abstract)

    [15] Yang Z, Zhan J, Wang C, et al. Coupling coordination analysis and spatiotemporal heterogeneity between sustainable development and ecosystem services in Shanxi Province, China[J]. Science of the Total Environment, 2022, 836: 155625.

    [16] Guo S, Wu C, Wang Y, et al. Threshold effect of ecosystem services in response to climate change, human activity and landscape pattern in the upper and middle Yellow River of China[J]. Ecological Indicators, 2022, 136: 108603.

    [17] Gong J, Liu D, Zhang J, et al. Tradeoffs/synergies of multiple ecosystem services based on land use simulation in a mountain-basin area, western China[J]. Ecological Indicators, 2019, 99: 283-293.

    [18] 黃孟勤,李陽(yáng)兵,李明珍,等. 三峽庫(kù)區(qū)人類活動(dòng)強(qiáng)度與景觀格局的耦合響應(yīng)[J]. 生態(tài)學(xué)報(bào),2022,42(10):3959-3973.

    Huang Mengqin, Li Yangbing, Li Mingzhen, et al. Coupling response of Human activity intensity and landscape pattern in the Three Gorges Reservoir Area[J]. Acta Ecologica Sinica, 2022, 42(10): 3959-3973. (in Chinese with English abstract)

    [19] 方林,蔡俊,劉艷曉,等. 長(zhǎng)三角地區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值動(dòng)態(tài)演化及驅(qū)動(dòng)力分析[J]. 生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2022,38(5):556-565.

    Fang Lin, Cai Jun, Liu Yanxiao, et al. Dynamic evolution of ecosystem service value in Yangtze River Delta and analysis on the driving factors[J]. Journal of Ecology and Rural Environment, 2022, 38(5): 556-565. (in Chinese with English abstract)

    [20] Liao G, He P, Gao X, et al. Land use optimization of rural production-living-ecological space at different scales based on the BP-ANN and CLUE-S models[J]. Ecological Indicators, 2022, 137: 108710.

    [21] 尉芳,劉京,夏利恒,等. 基于LUCC的陜西渭北旱塬區(qū)景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)[J]. 中國(guó)環(huán)境科學(xué),2022,42(4):1963-1974.

    Wei Fang, Liu Jing, Xia Liheng, et al. Landscape ecological risk assessment in Weibei dryland region of Shaanxi Province based on LUCC[J]. China Environmental Science,2022, 42(4): 1963-1974. (in Chinese with English abstract)

    [22] Huang X, Wang X, Zhang X, et al. Ecological risk assessment and identification of risk control priority areas based on degradation of ecosystem services: A case study in the Tibetan Plateau[J]. Ecological Indicators, 2022, 141: 109078.

    [23] 榮益,李超,許策,等. 城鎮(zhèn)化過程中生態(tài)系統(tǒng)服務(wù)價(jià)值變化及人類活動(dòng)影響的空間分異:以黃驊市為例[J]. 生態(tài)學(xué)雜志,2017,36(5):1374-1381.

    Rong Yi, Li Chao, Xu Ce, et al. Ecosystem service values and spatial differentiation changes during urbanization: A case study of Huanghua City[J]. Chinese Journal of Ecology, 2017, 36(5): 1374-1381. (in Chinese with English abstract)

    [24] 于淑會(huì),康園園,鄧偉,等. 太行山東部縣域“三生”用地轉(zhuǎn)型與景觀生態(tài)風(fēng)險(xiǎn)分析:以河北省平山縣為例[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)(中英文),2022,30(7):1113-1122.

    Yu Shuhui, Kang Yuanyuan, Deng Wei, et al. Analysis of “production-living-ecological” land transformation and landscape ecological risk in the eastern counties of the Taihang Mountain: A case study in Pingshan County[J]. Chinese Journal of Eco-Agriculture, 2022, 30(7): 1113-1122. (in Chinese with English abstract)

    [25] 董立寬,方斌,王晨歌. 基于Copula函數(shù)的茶園土壤銅鋅空間協(xié)同效應(yīng)研究[J]. 自然資源學(xué)報(bào),2018,33(5):867-878.

    Dong Likuan, Fang Bin, Wang Chen’ge. Study on the spatial synergistic effect of copper and Zincin tea garden soil based on Copula function[J]. Journal of Natural Resources, 2018, 33(5): 867-878. (in Chinese with English abstract)

    [26] Abakah E J A, Addo E, Gil-Alana L A, et al. Re-examination of international bond market dependence: Evidence from a pair copula approach[J]. International Review of Financial Analysis, 2021, 74: 101678.

    [27] 張齡月. 一類有關(guān)隨機(jī)向量間非線性相依度量的研究及其應(yīng)用[D]. 大連:大連理工大學(xué),2021.

    Zhang Linyue. The Measures of Nonlinear Dependence for Random Vectors and Their Application[D]. Dalian: Dalian University of Technology, 2021. (in Chinese with English abstract)

    [28] 魏立力,李賀. 基于Copula函數(shù)的CPI與PPI相關(guān)性分析[J]. 數(shù)學(xué)的實(shí)踐與認(rèn)識(shí),2020,50(6):1-7.

    Wei Lili, Li He. Dependence analysis of CPI and PPI based on Copula function[J]. Mathematics in Practice and Theory, 2020, 50(6): 1-7. (in Chinese with English abstract)

    [29] 朱從謀,李武艷,杜瑩瑩,等. 浙江省耕地多功能價(jià)值時(shí)空變化與權(quán)衡—協(xié)同關(guān)系[J]. 農(nóng)業(yè)工程學(xué)報(bào),2020,36(14):263-272.

    Zhu Congmou, Li Wuyan, Du Yingying, et al. Spatial-temporal change, trade-off and synergy relationships of cropland multifunctional value in Zhejiang Province, China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(14): 263-272. (in Chinese with English abstract)

    [30] 朱慶瑩,涂濤,陳銀蓉,等. 長(zhǎng)江經(jīng)濟(jì)帶產(chǎn)業(yè)結(jié)構(gòu)與土地利用效率協(xié)調(diào)度的時(shí)空格局分析[J]. 農(nóng)業(yè)工程學(xué)報(bào),2021,37(23):235-246.

    Zhu Qingying, Tu Tao, Chen Yinrong, et al. Spatio-temporal pattern for the coordination degree between industrial structure and land use efficiency of Yangtze River Economic Zone[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(23): 235-246. (in Chinese with English abstract)

    [31] 黃木易,岳文澤,方斌,等. 1970-2015年大別山區(qū)生態(tài)服務(wù)價(jià)值尺度響應(yīng)特征及地理探測(cè)機(jī)制[J]. 地理學(xué)報(bào),2019,74(9):1904-1920.

    Huang Muyi, Yue Wenze, Fang Bin, et al. Scale response characteristics and geographic exploration mechanism of spatial differentiation of ecosystem service values in Dabie Mountain area, central China from 1970 to 2015[J]. Acta Geographica Sinica, 2019, 74(9): 1904-1920. (in Chinese with English abstract)

    [32] Ran P, Hu S, Frazier A E E, et al. Exploring changes in landscape ecological risk in the Yangtze River Economic Belt from a spatiotemporal perspective[J]. Ecological Indicators, 2022, 137: 108744.

    [33] 張弛,李枝堅(jiān),曾輝. 基于小波變換的生態(tài)系統(tǒng)服務(wù)權(quán)衡尺度效應(yīng)與影響因素研究:以珠江三角洲為例[J]. 地理研究,2022,41(5):1279-1297.

    Zhang Chi, Li Zhijian, Zeng Hui. Scale effects on ecosystem service trade-off and its influencing factors based on wavelet transform: A case study in the Pearl River Delta, China[J]. Geographical Research, 2022, 41(5): 1279-1297. (in Chinese with English abstract)

    [34] 周婷,陳萬(wàn)旭,李江風(fēng),等. 神農(nóng)架林區(qū)人類活動(dòng)與生境質(zhì)量的空間關(guān)系[J]. 生態(tài)學(xué)報(bào),2021,41(15):6134-6145.

    Zhou Ting, Chen Wanxu, Li Jangfeng, et al. Spatial relationship between human activities and habitat quality in Shennongjia Forest Region from 1995 to 2015[J]. Acta Ecologica Sinica, 2021, 41(15): 6134-6145. (in Chinese with English abstract)

    [35] 李輝,紅英,鄧國(guó)榮,等. 1982—2015年氣候變化和人類活動(dòng)對(duì)內(nèi)蒙古草地凈初級(jí)生產(chǎn)力的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào),2021,32(2):415-424.

    Li Hui, Hong Ying, Deng Guorong, et al. Impacts of climate change and human activities on net primary productivity of grasslands in Inner Mongolia, China during 1982-2015[J]. Chinese Journal of Applied Ecology, 2021, 32(2): 415-424. (in Chinese with English abstract)

    [36] 陳改君,呂培亮. “生態(tài)正義”何以實(shí)現(xiàn)?——基于環(huán)境庫(kù)茲涅茨曲線的檢驗(yàn)性分析[J]. 湖南社會(huì)科學(xué),2022(3):45-51.

    Chen Gaijun, Lu Peiliang. “How to achieve ecological justice”? —Testability analysis based on environmental Kuznets curve[J]. Hunan Social Sciences, 2022(3): 45-51. (in Chinese with English abstract)

    [37] 劉曉紅,江可申. 基于靜態(tài)與動(dòng)態(tài)空間面板模型分析城鎮(zhèn)化對(duì)霧霾的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2017,33(20):218-225.

    Liu Xiaohong, Jiang Keshen. Analyzing effect of urbanization on haze pollution based on static and dynamic spatial panel model[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(20): 218-225. (in Chinese with English abstract)

    [38] 王子成,趙忠. 農(nóng)民工遷移模式的動(dòng)態(tài)選擇:外出、回流還是再遷移[J]. 管理世界,2013(1):78-88.

    Wang Zicheng, Zhao Zhong. The dynamic choice of migrant workers' migration mode: Going out, returning or re-migration[J]. Journal of Management World, 2013(1): 78-88. (in Chinese with English abstract)

    [39] 朱存泉. 長(zhǎng)三角地區(qū)勞動(dòng)力回流的動(dòng)因與對(duì)策研究[D]. 蘇州:蘇州大學(xué),2019.

    Zhu Cunquan. Research on the Causes and Countermeasures of Labor Return in Yangtze River Delta Region[D]. Suzhou: Soochow University, 2019. (in Chinese with English abstract)

    [40] 孫秋鵬. 經(jīng)濟(jì)高質(zhì)量發(fā)展對(duì)環(huán)境保護(hù)和生態(tài)文明建設(shè)的推動(dòng)作用[J]. 當(dāng)代經(jīng)濟(jì)管理,2019,41(11):9-14.

    Sun Qiupeng. The promoting role of high-quality economic development in environmental protection and ecological civilization construction[J]. Contemporary Economic Management, 2019, 41(11): 9-14. (in Chinese with English abstract)

    Scale response and spatiotemporal correlations between landscape ecological risk and human activity intensity in the Yangtze River Delta region

    Fang Lin1, Fang Bin2※, Liu Yanxiao3, Cai Jun1, Li Canfeng4

    (1.,,230036,; 2.,,210023,; 3.,,130024,; 4.,,650100,)

    Human activity has been one of the most important factors inducing regional ecological risks in recent years. But, there are different visualization and response degree of ecological risks and human activities in different land scales. It is a high demand to scientifically analyze the correlation characteristics of regional Landscape Ecological Risk (LER) and Human Activity Intensity (HAI), in order to coordinate the territorial system of human-environment interaction for better regional high-quality development. The ever-increasing conflict can be found between ecological protection and economic and social development, due mainly to rapid urbanization. The Yangtze River Delta (YRD) is one of the typical regions with intensive human activities, remarkable economic development, and outstanding ecological status. Taking the YRD as the subject, a three-level scale of city, county, and grid was established to construct the LER and HAI assessment models, in order to characterize the spatial and temporal response at different scales. The data was also collected from the land use, nighttime lighting, and population spatial distribution in the four periods from 1990 to 2020. The Copula function, bivariate spatial autocorrelation, and coupled coordination degree models were used to reveal the spatial and temporal correlation between the LER and HAI. The results show that (1) the spatial distribution pattern of LER was high in the north of the study area, and low in the south under different scales. The high risk areas continued to decrease during 1990-2020, whereas, the low risk areas showed an increasing trend. There was the most outstanding reduction in the high risk areas at the city scale, with a decrease of 19.51 percentage points over 30 years. By contrast, there was the most significant increase in the low-medium risk and low risk areas at the county scale, indicating an increase of 16.29 percentage points over 30 years. (2) The spatial distribution pattern of HAI was high in the northeast of the study area during 1990-2020, while low in the southwest under different scales. All the regions with high HAI showed a significant increasing trend. There was a significant decreasing trend in the regions with low HAI. The high intensity region presented the largest increase at 13.42 percentage points on the grid scale, whereas, the low intensity region was the most significant decrease at 9.76 percentage points on the city scale. (3) The correlation between HAI and LER shared a positive correlation between 1990 and 2020, but the correlation between them was weakening, indicating the weak influence of regional HAI on LER. By contrast, the coupling and coordination between them showed an increasing trend from 0.3031 in 1990 to 0.3112 in 2020. It infers that the relationship between them was gradually shifting from conflict to coordination. Spatially, there was a continuous decreasing trend in the H-H clustering area. Furthermore, the L-L clustering area showed an increasing and then decreasing trend with an overall decrease. Both L-H and H-L increased significantly after 2010. The spatial correlation characteristics of regional LER and HAI were combined to manage and control the regional ecological environment. The finding can provide a theoretical basis for land management and landscape planning, together with the spatial guidelines for the regional ecological risk prevention and ecological restoration.

    land use; landscape; ecology; human activity; scale response; copula function; Yangtze River Delta

    10.11975/j.issn.1002-6819.2022.22.023

    X321

    A

    1002-6819(2022)-22-0210-10

    方林,方斌,劉艷曉,等. 長(zhǎng)三角地區(qū)景觀生態(tài)風(fēng)險(xiǎn)與人類活動(dòng)強(qiáng)度的響應(yīng)及其時(shí)空關(guān)聯(lián)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2022,38(22):210-219.doi:10.11975/j.issn.1002-6819.2022.22.023 http://www.tcsae.org

    Fang Lin, Fang Bin, Liu Yanxiao, et al. Scale response and spatiotemporal correlations between landscape ecological risk and human activity intensity in the Yangtze River Delta region[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(22): 210-219. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2022.22.023 http://www.tcsae.org

    2022-08-27

    2022-09-27

    國(guó)家自然科學(xué)基金云南聯(lián)合基金重點(diǎn)項(xiàng)目(U2102209);國(guó)家自然科學(xué)基金項(xiàng)目(42071229,41671174)

    方林,研究方向?yàn)橥恋乩门c生態(tài)系統(tǒng)服務(wù)。Email:fl342427@126.com

    方斌,教授,博士生導(dǎo)師,研究方向?yàn)楦乇Wo(hù)、土地利用及土地生態(tài)。Email:wenyanfang731@163.com

    猜你喜歡
    長(zhǎng)三角時(shí)空尺度
    跨越時(shí)空的相遇
    緊扣一體化 全面融入長(zhǎng)三角
    鏡中的時(shí)空穿梭
    財(cái)產(chǎn)的五大尺度和五重應(yīng)對(duì)
    “首屆長(zhǎng)三角新青年改稿會(huì)”作品選
    玩一次時(shí)空大“穿越”
    宇宙的尺度
    太空探索(2016年5期)2016-07-12 15:17:55
    時(shí)空之門
    長(zhǎng)三角瞭望
    長(zhǎng)三角瞭望
    亚洲av第一区精品v没综合| 久久亚洲国产成人精品v| 欧美国产日韩亚洲一区| 亚洲av不卡在线观看| 国产大屁股一区二区在线视频| 欧美人与善性xxx| 你懂的网址亚洲精品在线观看 | 亚洲美女黄片视频| 免费搜索国产男女视频| 亚洲国产色片| 国产伦精品一区二区三区四那| 中出人妻视频一区二区| 成人性生交大片免费视频hd| 禁无遮挡网站| 插阴视频在线观看视频| 丝袜美腿在线中文| 久久6这里有精品| 美女xxoo啪啪120秒动态图| 看非洲黑人一级黄片| 国产白丝娇喘喷水9色精品| 啦啦啦啦在线视频资源| 免费观看人在逋| 亚洲精品一区av在线观看| 国产色婷婷99| 国产伦一二天堂av在线观看| 亚洲av中文av极速乱| 69av精品久久久久久| 国产精品嫩草影院av在线观看| 久久久a久久爽久久v久久| 婷婷精品国产亚洲av| 精品午夜福利在线看| 一边摸一边抽搐一进一小说| 级片在线观看| 22中文网久久字幕| 最新中文字幕久久久久| 九九久久精品国产亚洲av麻豆| 3wmmmm亚洲av在线观看| 淫秽高清视频在线观看| 久久精品夜色国产| 久久精品国产鲁丝片午夜精品| 可以在线观看毛片的网站| av在线亚洲专区| 亚洲经典国产精华液单| 欧美激情久久久久久爽电影| 丝袜美腿在线中文| 亚洲图色成人| 无遮挡黄片免费观看| 深夜a级毛片| 一级a爱片免费观看的视频| 日本五十路高清| 男人舔女人下体高潮全视频| 亚洲国产欧洲综合997久久,| 听说在线观看完整版免费高清| 一级毛片久久久久久久久女| 亚洲七黄色美女视频| 女人十人毛片免费观看3o分钟| 校园人妻丝袜中文字幕| 久久久久久久午夜电影| 俄罗斯特黄特色一大片| 高清毛片免费观看视频网站| 久久久久国内视频| 国产精品一二三区在线看| 国产大屁股一区二区在线视频| 国产伦精品一区二区三区四那| 国产亚洲精品av在线| 亚洲人成网站在线观看播放| 亚洲欧美成人精品一区二区| 国产精品国产高清国产av| 小说图片视频综合网站| 22中文网久久字幕| 亚洲成人中文字幕在线播放| 男人舔女人下体高潮全视频| 精品久久国产蜜桃| 搡老妇女老女人老熟妇| 春色校园在线视频观看| 少妇人妻精品综合一区二区 | 国内精品宾馆在线| 精品一区二区三区视频在线| 人妻久久中文字幕网| 黄色欧美视频在线观看| 日韩在线高清观看一区二区三区| 欧美最新免费一区二区三区| 男女啪啪激烈高潮av片| 亚洲电影在线观看av| 夜夜看夜夜爽夜夜摸| 久久久午夜欧美精品| 久久亚洲精品不卡| 久久久久久大精品| 69av精品久久久久久| 男女之事视频高清在线观看| 99国产精品一区二区蜜桃av| 国产成人影院久久av| 亚洲精品国产成人久久av| 国产精品一区二区三区四区免费观看 | 亚洲精品久久国产高清桃花| 丝袜美腿在线中文| 国产高清有码在线观看视频| 亚洲无线在线观看| 亚洲精品色激情综合| 久久久成人免费电影| 国产熟女欧美一区二区| 成人一区二区视频在线观看| 亚洲人与动物交配视频| 观看美女的网站| 欧美日韩在线观看h| 久久久久国产网址| 国产精品久久电影中文字幕| 国产精华一区二区三区| 最好的美女福利视频网| 波野结衣二区三区在线| 午夜精品一区二区三区免费看| 麻豆精品久久久久久蜜桃| 乱人视频在线观看| 国产精品av视频在线免费观看| 亚洲激情五月婷婷啪啪| 日本成人三级电影网站| 天堂影院成人在线观看| 99久久成人亚洲精品观看| 亚洲精品乱码久久久v下载方式| 观看免费一级毛片| 特级一级黄色大片| 欧美一区二区精品小视频在线| 国产精品一及| 一区二区三区免费毛片| 又爽又黄无遮挡网站| 天堂影院成人在线观看| 精品久久久久久久久久久久久| 亚洲在线自拍视频| 精品久久久噜噜| 亚洲天堂国产精品一区在线| 老女人水多毛片| 日韩在线高清观看一区二区三区| 精品熟女少妇av免费看| 久久国产乱子免费精品| 国产av不卡久久| 久久6这里有精品| 一a级毛片在线观看| 少妇被粗大猛烈的视频| 国产精品一及| 日本a在线网址| 99久国产av精品国产电影| 麻豆av噜噜一区二区三区| 中文字幕久久专区| 国产精品1区2区在线观看.| 国产亚洲精品久久久久久毛片| 又黄又爽又刺激的免费视频.| 校园春色视频在线观看| 能在线免费观看的黄片| 少妇裸体淫交视频免费看高清| 欧美性猛交黑人性爽| 深夜精品福利| 国产私拍福利视频在线观看| 搡老岳熟女国产| 一区二区三区高清视频在线| 亚洲性久久影院| 国产探花在线观看一区二区| 日本在线视频免费播放| 美女 人体艺术 gogo| 在现免费观看毛片| 卡戴珊不雅视频在线播放| 亚洲最大成人av| 亚洲18禁久久av| 最好的美女福利视频网| 日本-黄色视频高清免费观看| 国语自产精品视频在线第100页| 少妇的逼水好多| 啦啦啦韩国在线观看视频| 亚洲欧美日韩东京热| 国产高清视频在线播放一区| 日韩中字成人| 久久精品综合一区二区三区| 日韩欧美精品免费久久| 极品教师在线视频| 国产在线男女| 日本黄色视频三级网站网址| 搡老岳熟女国产| 亚洲无线在线观看| 丝袜美腿在线中文| 99在线视频只有这里精品首页| 最近的中文字幕免费完整| 特大巨黑吊av在线直播| 国产高清三级在线| 欧美最黄视频在线播放免费| 久久久久久久久中文| 国产亚洲精品综合一区在线观看| 日韩欧美在线乱码| 尾随美女入室| 色5月婷婷丁香| av福利片在线观看| 亚洲成人久久性| 欧美一区二区精品小视频在线| 国产一级毛片七仙女欲春2| 国产极品精品免费视频能看的| 三级国产精品欧美在线观看| 悠悠久久av| 亚洲av中文av极速乱| 精品人妻视频免费看| 久久久成人免费电影| 日本精品一区二区三区蜜桃| 97在线视频观看| 噜噜噜噜噜久久久久久91| or卡值多少钱| 两个人视频免费观看高清| 国产精品久久久久久av不卡| av女优亚洲男人天堂| 亚洲人与动物交配视频| 在线看三级毛片| 天堂网av新在线| 99精品在免费线老司机午夜| 国内精品一区二区在线观看| 插逼视频在线观看| 蜜臀久久99精品久久宅男| 一区福利在线观看| 国产一区二区激情短视频| 我要搜黄色片| 俺也久久电影网| 久久人人爽人人片av| 老司机福利观看| 久久人妻av系列| 亚洲中文日韩欧美视频| 国产精品一区二区三区四区免费观看 | 国产极品精品免费视频能看的| 欧美+亚洲+日韩+国产| 婷婷色综合大香蕉| 久久中文看片网| 91久久精品国产一区二区三区| 女人十人毛片免费观看3o分钟| 久久九九热精品免费| 韩国av在线不卡| 欧美另类亚洲清纯唯美| 欧美一级a爱片免费观看看| 欧美日本视频| 网址你懂的国产日韩在线| 久久精品91蜜桃| 亚洲美女黄片视频| 精品乱码久久久久久99久播| 亚洲最大成人av| 欧美成人a在线观看| 亚洲欧美日韩高清专用| 极品教师在线视频| 国语自产精品视频在线第100页| 人人妻人人澡人人爽人人夜夜 | 久久久久久久久久成人| 国产午夜精品论理片| 国产精品乱码一区二三区的特点| 网址你懂的国产日韩在线| av.在线天堂| 免费在线观看影片大全网站| 亚洲av中文av极速乱| 91av网一区二区| 久久这里只有精品中国| 国产一区亚洲一区在线观看| 久久精品综合一区二区三区| 少妇裸体淫交视频免费看高清| 69av精品久久久久久| 久久久久九九精品影院| av女优亚洲男人天堂| 在线看三级毛片| 欧美另类亚洲清纯唯美| 又黄又爽又刺激的免费视频.| 国产精品电影一区二区三区| av在线亚洲专区| 黄色日韩在线| 我的老师免费观看完整版| 国产精品久久久久久亚洲av鲁大| 美女大奶头视频| 99久久精品热视频| 国产欧美日韩精品亚洲av| 欧美性感艳星| 亚洲婷婷狠狠爱综合网| 99在线视频只有这里精品首页| a级毛色黄片| 欧洲精品卡2卡3卡4卡5卡区| 久久久精品大字幕| 国产爱豆传媒在线观看| 国产男人的电影天堂91| 少妇裸体淫交视频免费看高清| 国产黄色视频一区二区在线观看 | av视频在线观看入口| 老司机福利观看| 一个人看的www免费观看视频| 此物有八面人人有两片| 男女边吃奶边做爰视频| 国产片特级美女逼逼视频| 成人一区二区视频在线观看| 22中文网久久字幕| 久久久久久国产a免费观看| 国产麻豆成人av免费视频| 亚洲中文日韩欧美视频| 国产精品无大码| 欧美zozozo另类| 亚洲第一电影网av| 国产精品久久久久久久电影| 国产欧美日韩精品亚洲av| 51国产日韩欧美| 欧美激情国产日韩精品一区| 男女下面进入的视频免费午夜| 在线观看午夜福利视频| 青春草视频在线免费观看| 亚洲中文字幕一区二区三区有码在线看| 99热只有精品国产| 日韩亚洲欧美综合| 青春草视频在线免费观看| 亚洲电影在线观看av| 亚洲18禁久久av| 久久国内精品自在自线图片| 国产大屁股一区二区在线视频| 免费观看在线日韩| 亚洲内射少妇av| 午夜爱爱视频在线播放| 日韩成人伦理影院| avwww免费| 最近在线观看免费完整版| 亚洲第一电影网av| 免费观看在线日韩| 亚洲一级一片aⅴ在线观看| 日本欧美国产在线视频| 99热这里只有是精品50| 在线观看av片永久免费下载| 久久久精品94久久精品| 国产精品一二三区在线看| 91狼人影院| 国内揄拍国产精品人妻在线| 我要搜黄色片| 国产精品一区www在线观看| 别揉我奶头 嗯啊视频| 亚洲欧美日韩东京热| 久久欧美精品欧美久久欧美| 日韩三级伦理在线观看| 六月丁香七月| 午夜久久久久精精品| 在线天堂最新版资源| 一区福利在线观看| 国产精品精品国产色婷婷| 三级毛片av免费| 亚洲在线观看片| 成人毛片a级毛片在线播放| 免费搜索国产男女视频| 成人高潮视频无遮挡免费网站| 国产一区二区激情短视频| 国产 一区 欧美 日韩| 欧美日韩精品成人综合77777| 国产麻豆成人av免费视频| 精品熟女少妇av免费看| 久久韩国三级中文字幕| 天天一区二区日本电影三级| 欧美一区二区亚洲| 在线国产一区二区在线| 国产成人福利小说| 人妻丰满熟妇av一区二区三区| 日韩欧美一区二区三区在线观看| 精品日产1卡2卡| 人人妻,人人澡人人爽秒播| 日本欧美国产在线视频| 免费观看在线日韩| 亚洲精品国产av成人精品 | 亚洲内射少妇av| 丝袜喷水一区| 国产精品一及| 乱码一卡2卡4卡精品| 亚洲精品成人久久久久久| aaaaa片日本免费| 欧美一区二区精品小视频在线| 最近手机中文字幕大全| 国产熟女欧美一区二区| 性色avwww在线观看| 黄色配什么色好看| 日本a在线网址| 亚洲国产日韩欧美精品在线观看| 国产人妻一区二区三区在| 大香蕉久久网| av视频在线观看入口| 国产精品亚洲美女久久久| a级毛片免费高清观看在线播放| 亚洲欧美日韩高清在线视频| 亚洲av美国av| 国产在线精品亚洲第一网站| 国产蜜桃级精品一区二区三区| 久久精品国产99精品国产亚洲性色| 亚洲精品亚洲一区二区| 国产麻豆成人av免费视频| 麻豆久久精品国产亚洲av| 美女内射精品一级片tv| 亚洲经典国产精华液单| 亚洲高清免费不卡视频| 亚洲av二区三区四区| 久久6这里有精品| 内射极品少妇av片p| 亚洲最大成人av| 真人做人爱边吃奶动态| 成人特级av手机在线观看| 欧美一级a爱片免费观看看| 日日撸夜夜添| 久久人妻av系列| 香蕉av资源在线| 直男gayav资源| 男女下面进入的视频免费午夜| 久久草成人影院| 国产一区二区三区av在线 | 欧美人与善性xxx| 老女人水多毛片| 色综合亚洲欧美另类图片| 精品久久久久久成人av| 搡女人真爽免费视频火全软件 | 久久久午夜欧美精品| 狂野欧美白嫩少妇大欣赏| 22中文网久久字幕| 久久精品夜色国产| 午夜福利在线在线| 国产日本99.免费观看| 男女啪啪激烈高潮av片| 久久人人爽人人片av| 欧美+亚洲+日韩+国产| 能在线免费观看的黄片| 97碰自拍视频| 在线看三级毛片| 婷婷精品国产亚洲av在线| 少妇人妻精品综合一区二区 | 男人舔女人下体高潮全视频| 亚洲人与动物交配视频| 在线观看免费视频日本深夜| 噜噜噜噜噜久久久久久91| 久久人妻av系列| 国产一区二区三区av在线 | 两性午夜刺激爽爽歪歪视频在线观看| а√天堂www在线а√下载| 国产成人一区二区在线| 最近中文字幕高清免费大全6| 一个人看视频在线观看www免费| 一个人免费在线观看电影| 两性午夜刺激爽爽歪歪视频在线观看| 两个人视频免费观看高清| 国产精品乱码一区二三区的特点| 91久久精品国产一区二区成人| 婷婷亚洲欧美| 中国美女看黄片| 国产高清不卡午夜福利| 禁无遮挡网站| 成人av一区二区三区在线看| 国产人妻一区二区三区在| 午夜福利高清视频| aaaaa片日本免费| 国产精品福利在线免费观看| 成人鲁丝片一二三区免费| 国产一区二区三区av在线 | 少妇高潮的动态图| 18禁在线无遮挡免费观看视频 | 51国产日韩欧美| 嫩草影视91久久| aaaaa片日本免费| 亚洲美女黄片视频| 亚洲av电影不卡..在线观看| 国产伦精品一区二区三区视频9| 内地一区二区视频在线| www.色视频.com| 最近最新中文字幕大全电影3| 久久久久久久久久黄片| 精品国内亚洲2022精品成人| a级毛片免费高清观看在线播放| 午夜日韩欧美国产| 久久综合国产亚洲精品| 亚洲七黄色美女视频| 小蜜桃在线观看免费完整版高清| 99九九线精品视频在线观看视频| 久久精品国产亚洲av涩爱 | 亚洲美女搞黄在线观看 | 少妇的逼水好多| 亚洲最大成人av| 国产伦在线观看视频一区| 男女之事视频高清在线观看| 亚洲经典国产精华液单| 99久久精品国产国产毛片| 国产女主播在线喷水免费视频网站 | 免费看光身美女| 少妇猛男粗大的猛烈进出视频 | 日韩欧美三级三区| 亚洲在线观看片| 亚洲国产欧美人成| 伦精品一区二区三区| 亚洲av熟女| 三级男女做爰猛烈吃奶摸视频| 一本久久中文字幕| 国产高清视频在线播放一区| a级毛色黄片| 亚洲精品成人久久久久久| 欧洲精品卡2卡3卡4卡5卡区| 久久久久久久久久黄片| 精品久久久久久久人妻蜜臀av| av在线天堂中文字幕| 亚洲欧美日韩东京热| 亚洲高清免费不卡视频| 日本黄大片高清| 最近2019中文字幕mv第一页| 无遮挡黄片免费观看| 国产精品久久电影中文字幕| 国产精品免费一区二区三区在线| 亚洲丝袜综合中文字幕| 国产精品女同一区二区软件| 黄色一级大片看看| 午夜福利视频1000在线观看| 亚洲激情五月婷婷啪啪| 99久国产av精品| 久久久成人免费电影| 精品久久久噜噜| 床上黄色一级片| 午夜精品一区二区三区免费看| 国产精品久久久久久精品电影| 最近2019中文字幕mv第一页| 欧美激情国产日韩精品一区| 国国产精品蜜臀av免费| 在线播放国产精品三级| 无遮挡黄片免费观看| 久久九九热精品免费| 久久久久国产精品人妻aⅴ院| a级毛片免费高清观看在线播放| 国产精品一区二区免费欧美| 国产精品av视频在线免费观看| 91精品国产九色| 国产三级在线视频| 狂野欧美白嫩少妇大欣赏| 麻豆久久精品国产亚洲av| 亚洲精品久久国产高清桃花| 国产精品国产高清国产av| 偷拍熟女少妇极品色| 国产亚洲精品综合一区在线观看| 长腿黑丝高跟| 久久久久久久久久久丰满| 日本色播在线视频| 国产亚洲欧美98| 91午夜精品亚洲一区二区三区| av女优亚洲男人天堂| 欧美日韩一区二区视频在线观看视频在线 | 久久久久久大精品| 人妻制服诱惑在线中文字幕| 俺也久久电影网| 男女啪啪激烈高潮av片| 亚洲av不卡在线观看| 日韩欧美三级三区| 日韩在线高清观看一区二区三区| 久久韩国三级中文字幕| 国产免费一级a男人的天堂| 成年免费大片在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 国产黄色视频一区二区在线观看 | 国产精品久久久久久亚洲av鲁大| 欧美另类亚洲清纯唯美| 国产一区亚洲一区在线观看| 国产亚洲精品av在线| 99久久无色码亚洲精品果冻| 国产女主播在线喷水免费视频网站 | 99久久久亚洲精品蜜臀av| 久久综合国产亚洲精品| 天堂av国产一区二区熟女人妻| 国产高清视频在线观看网站| 精品久久久久久久久久久久久| 亚洲欧美成人综合另类久久久 | 丰满的人妻完整版| 成年免费大片在线观看| ponron亚洲| 男人和女人高潮做爰伦理| 欧美不卡视频在线免费观看| 我要搜黄色片| 日韩欧美在线乱码| 色综合站精品国产| 看片在线看免费视频| АⅤ资源中文在线天堂| 欧美最黄视频在线播放免费| 国产男人的电影天堂91| 精品一区二区三区av网在线观看| 黑人高潮一二区| 亚洲在线自拍视频| 国产高清视频在线观看网站| 成人美女网站在线观看视频| 欧美精品国产亚洲| 啦啦啦韩国在线观看视频| 国产精品三级大全| 小蜜桃在线观看免费完整版高清| 成人特级av手机在线观看| 精品久久久久久久人妻蜜臀av| 成人特级av手机在线观看| 午夜日韩欧美国产| 免费观看人在逋| 亚洲成人av在线免费| 狂野欧美白嫩少妇大欣赏| 热99在线观看视频| 亚洲av中文字字幕乱码综合| 99久久无色码亚洲精品果冻| 午夜福利高清视频| 成年av动漫网址| 亚洲av.av天堂| 可以在线观看毛片的网站| 日本a在线网址| 身体一侧抽搐| 老司机午夜福利在线观看视频| 五月玫瑰六月丁香| 真人做人爱边吃奶动态| 久久久精品94久久精品| 国产精品人妻久久久久久| 欧美性猛交╳xxx乱大交人| 欧美高清成人免费视频www| 中文字幕熟女人妻在线| 中文资源天堂在线| 日本免费a在线| 日本熟妇午夜| 卡戴珊不雅视频在线播放| 国语自产精品视频在线第100页| 18禁在线无遮挡免费观看视频 | 日本免费a在线| 少妇高潮的动态图| 中国美女看黄片| 国产一级毛片七仙女欲春2| 国产片特级美女逼逼视频| 又粗又爽又猛毛片免费看| 成年免费大片在线观看| 五月伊人婷婷丁香|