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

    基于MSPA-CIRCUIT的長江中游城市群熱環(huán)境網絡識別與評價研究

    2023-10-07 08:37:54方云皓趙麗元竇碧瑩王書賢
    生態(tài)環(huán)境學報 2023年7期
    關鍵詞:環(huán)境

    方云皓,趙麗元*,竇碧瑩,王書賢

    1.華中科技大學建筑與城市規(guī)劃學院,湖北 武漢 430074;2.湖北省城鎮(zhèn)化工程技術研究中心,湖北 武漢 430074;3.南京市規(guī)劃與自然資源局江寧分局,江蘇 南京 211100

    城市熱環(huán)境指與熱相關的、影響居民行為的物理環(huán)境,通?;跍囟?、熱輻射量、熱舒適度等物理參數(shù)予以定量化描述(孫鐵鋼等,2016)。在快速城市化背景下,城市規(guī)模不斷增加,由此引發(fā)城市空間環(huán)境如土地利用類型、功能設施布局、建筑空間形態(tài)等發(fā)生巨大變化,對城市熱環(huán)境產生顯著影響(Ma et al.,2010)。因此,城市熱環(huán)境也被視為城市空間環(huán)境在復雜熱力場中的綜合表現(xiàn)(辛儒鴻等,2022)。城市熱環(huán)境惡化已成為城市最難以解決的氣候問題之一,具體表征為高溫、熱浪、城市熱島等現(xiàn)象,其不僅改變了城市下墊面的生態(tài)環(huán)境,如植被生長、空氣質量、水文循環(huán)等,還影響了城鄉(xiāng)居民的生產空間、生活行為、能源消耗等(He et al.,2020)。此外,持續(xù)的熱環(huán)境惡化也會導致城鄉(xiāng)居民熱舒適性降低,高強度頻發(fā)時甚至對居民健康構成威脅,引發(fā)城鄉(xiāng)居民的發(fā)病率和死亡率增加(Lai et al.,2019),尤其對老年人與嬰幼兒來說,其更易受到熱環(huán)境惡化帶來的負面影響(Zuo et al.,2015)。因此,在應對全球氣候逐漸變暖的城市建設中,如何緩解熱環(huán)境惡化已成為中國政府、學者以及公眾等社會各界密切關注的重要議題。

    針對城市熱環(huán)境惡化問題,眾多學者基于地理學、生態(tài)學、城鄉(xiāng)規(guī)劃學等多學科視角在微觀尺度(人體尺度、房間尺度及建筑尺度)與宏觀尺度(街區(qū)尺度、城市尺度及區(qū)域尺度)將城市熱環(huán)境識別作為重點展開研究,主要原因有兩點:首先,城市熱環(huán)境識別研究涵蓋熱環(huán)境的時空格局、結構特征及演變過程等內容,是精準研判高溫災害險情、構建氣候安全框架、制定針對性氣候適應政策的主要依據(jù);其次,城市熱環(huán)境識別研究是揭示熱環(huán)境形成機制(主要基于土地利用類型、功能設施布局及建筑空間形態(tài)等層面闡述)、評估熱環(huán)境風險、模擬和預測未來熱環(huán)境發(fā)展趨勢等研究的必要前提(劉璐等,2019;付含聰?shù)龋?020;Lu et al.,2021;Chen et al.,2022;王子安等,2022),對于緩解城市熱環(huán)境惡化具有重要意義。

    當前針對城市熱環(huán)境識別研究主要包括熱環(huán)境時空格局分析與熱環(huán)境結構特征提取兩類。熱環(huán)境時空格局研究主要借助固定觀測(氣象站點、環(huán)境監(jiān)測點等形式)、流動觀測(布點觀測、探空氣球、流動車作業(yè)等形式)及衛(wèi)星遙感等技術獲取多時段氣溫或地表溫度(Land Surface Temperature,LST)數(shù)據(jù)(Unger et al.,2001;Stewart,2011;徐涵秋,2015;李召良等,2016),通過相關時空分析方法對熱環(huán)境的時空分布情況進行系統(tǒng)解析,具體來說可分為3種:第1種研究主要針對城市或城市群熱環(huán)境內部熱源中心的坐標屬性、地理區(qū)位等空間信息,通過重心遷移模型識別熱源中心并分析其時空遷移軌跡(孫芹芹等,2010;朱菲等,2022);第2種研究側重于通過剖面線方法設置道路或熱場剖面等作為城市“切割線”,比較不同方向熱環(huán)境時空差異性(徐麗華等,2007;于明洋等,2021);第三者則偏向于通過冷熱點分析模型、全局或局部自相關模型、支持向量機(Support Vector Machine,SVM)模型等對熱環(huán)境的變化趨勢與集聚特征進行測度(占文鳳等,2011;王鵬龍等,2014;陳冰倩等,2019)。熱環(huán)境結構特征研究是在時空格局分析的基礎上,進一步對城市內部各組成部分的熱環(huán)境分異規(guī)律進行挖掘的研究,涵蓋空間形態(tài)學與景觀生態(tài)學等多種理論與方法。就空間形態(tài)學而言,其主要依據(jù)土地利用類型、地表材質和人類活動等將城市下墊面劃分成不同類型區(qū)域,以評估并比較不同類型區(qū)域的熱環(huán)境空間差異,常見的分析方法包括局地氣候分區(qū)(Local Climate Zones,LCZ)框架、形態(tài)學空間格局分析(Morphological Spatial Pattern Analysis,MSPA)模型等(林中立等,2022;劉婷等,2023)。景觀生態(tài)學理論將城市熱環(huán)境視作由不同類型的斑塊組成,通過識別不同類型斑塊并計算其相應景觀格局指數(shù)(黃聚聰?shù)龋?011,2012;郭冠華等,2012),包括數(shù)量、面積、范圍、平均密度、類型比例、香濃多樣性指數(shù)等,以精準捕捉熱環(huán)境空間的演變、轉移、蔓延等特征,為改善熱環(huán)境提供精細化指導。

    盡管基于空間形態(tài)學與景觀生態(tài)學理論的熱環(huán)境結構特征研究取得豐碩成果,但其較多基于斑塊或空間單元視角展開,將熱環(huán)境斑塊或空間單元視為獨立個體而忽略了斑塊與斑塊或空間單元與空間單元之間的連通性。有證據(jù)表明,城市或城市群熱環(huán)境的連通性是影響區(qū)域內部溫度空間分布規(guī)律、城鄉(xiāng)居民熱舒適度的重要因素(Peng et al.,2020;吳健生等,2022),而識別具備連通性特征的熱環(huán)境廊道與熱環(huán)境網絡是緩解城市熱環(huán)境惡化、提升城鄉(xiāng)居民健康水平的有效手段(Yu et al.,2021;Peng et al.,2022),其主要原因為伴隨快速城市化的推進,熱環(huán)境的空間形態(tài)也隨之發(fā)生改變,尤其在城市群地區(qū),相鄰城市間的經濟、產業(yè)、交通等聯(lián)系日益增強,其熱環(huán)境會隨之連通化,由此引發(fā)熱量流通及擴散路徑也會發(fā)生變化(謝志清等,2007)。相關研究指出,熱環(huán)境的空間連通性越強,熱量流動阻力越小而擴散區(qū)域則越廣,由此易引發(fā)對更大范圍的城鄉(xiāng)居民健康構成威脅(江頌等,2022)。此外,通過電路(CIRCUIT)理論等建模方法對熱環(huán)境連通性進行模擬也被證實為一種可行的途徑,其通過模擬熱量的流動、擴散、聚集過程,為城市熱環(huán)境的監(jiān)測、預測與防范提供科學依據(jù)(喬治等,2023)。鑒于此,本研究以長江中游城市群為例,從時空動態(tài)視角探索熱環(huán)境網絡的時空演變特征,并評估其總體連通性與空間連通性,以期為區(qū)域尺度的熱環(huán)境惡化緩解方案提供氣候適應性框架。

    1 研究區(qū)域

    長江中游城市群(20°09′-33°20′N,180°21′-118°28′E)位于中國長江流域中部,承東啟西、連南接北,是長江經濟帶的重要組成部分,也是推動中部地區(qū)崛起的重要區(qū)域。根據(jù)2015年國務院批復的《長江中游城市群發(fā)展規(guī)劃》中劃定的行政區(qū)劃范圍,明確長江中游城市群是地跨湖北、湖南和江西3省共31市的特大型國家級城市群(圖1)。長江中游城市群國土面積約為3.26×105km2,地勢總體較為平坦,地形以平原與山體為主,屬亞熱帶季風性氣候,年均降水量在800-1943 mm,2020年總人口約為1.22×108人,地區(qū)生產總值達9.39×1012元。長江中游城市群的快速發(fā)展引發(fā)了城市熱環(huán)境的惡化。2000-2015年期間,城市群夏季白天平均熱島強度約0.8 ℃,夜間平均熱島強度約0.4 ℃,對城鄉(xiāng)居民的熱舒適度與身心健康構成了不同程度威脅(胡楠林等,2022)。

    圖1 研究區(qū)域Figure 1 Study area

    2 研究數(shù)據(jù)與方法

    2.1 研究數(shù)據(jù)及其預處理

    本研究的數(shù)據(jù)來源包括建成區(qū)面積數(shù)據(jù)與地表溫度數(shù)據(jù)。建成區(qū)面積數(shù)據(jù)來源于科學數(shù)據(jù)銀行提供的中國433個城市建成區(qū)邊界矢量數(shù)據(jù)(https://doi.org/10.11922/sciencedb.j00001.00332),時間分辨率為1990-2020(每隔5年)共7個年份,其主要用于城市熱島區(qū)域以及非熱島區(qū)域的劃分;地表溫度數(shù)據(jù)來源于中國科學院資源環(huán)境科學與數(shù)據(jù)中心提供的基于MODIS衛(wèi)星的中國地表溫度柵格數(shù)據(jù)集(https://www.resdc.cn/DOI/DOI.aspx?DOIID=98),其空間分辨率為1 km×1 km,時間分辨率為2000-2020年逐季度。本研究首先選取2000、2010以及2020年夏季(6-8月)的白天地表溫度數(shù)據(jù)作為源數(shù)據(jù),其次在Arcgis 10.2軟件中對源數(shù)據(jù)進行投影、裁剪等處理,最后將裁剪后的數(shù)據(jù)通過公式(1)以及柵格計算器進行單位轉換,以將熱力學溫度數(shù)據(jù)集轉換為攝氏溫度數(shù)據(jù)集,其表達式如下:

    式中:

    t1——攝氏溫度,單位為℃;

    T2——熱力學溫度,單位為K。

    2.2 研究方法

    2.2.1 基于MSPA模型的熱島斑塊類型劃分

    MSPA模型是基于數(shù)學形態(tài)學原理(如腐蝕、膨脹、開運算和閉運算等)對柵格圖像進行度量、識別以及分割的方法,其能夠依據(jù)柵格圖像的空間分異規(guī)律以及連通性視角,精準捕捉柵格圖像中的結構度量元素并提取7種類別互斥的像元斑塊,分別為核心、孤島、孔隙、邊緣、環(huán)道、橋接、支線(Soille et al.,2008)。MSPA模型由于偏向測度空間結構連接性,且具有斑塊提取準確度高以及分析結果可空間化等優(yōu)點,被廣泛運用于地理學、生態(tài)學以及遙感信息等領域的空間斑塊提取研究中。

    本研究基于MSPA模型劃分城市群熱島斑塊類型。具體來說,首先,分別統(tǒng)計2000年、2010年以及2020年夏季長江中游城市群平均地表溫度并依據(jù)其閾值進行等級劃分,由低到高依次為低溫區(qū)、次低溫區(qū)、常溫區(qū)、次高溫區(qū)、高溫區(qū);其次,高溫區(qū)域并非等同于熱島區(qū)域,綜合考慮城市群熱島的內涵及形成機制,為精準識別由城市化導致的建成區(qū)內部及其鄰接范圍的熱島區(qū)域并盡可能排除非城市化的影響,基于城市建成區(qū)邊界數(shù)據(jù)設置建成區(qū)及其外部2 km緩沖區(qū)范圍內的高溫區(qū)為熱島區(qū)域,其余區(qū)域為非熱島區(qū)域(劉婷等,2023);再次,利用Arcgis 10.2軟件的重分類工具將熱島區(qū)域與非熱島區(qū)域設置為二元柵格分析底圖,其中熱島區(qū)域為前景,非熱島區(qū)域為背景;最后,將二元柵格分析底圖導入Guidos軟件,設置鄰域轉換規(guī)則為數(shù)值8、邊緣寬度為數(shù)值2,通過數(shù)理分析方法將前景重新解譯為互不重疊的7類熱島斑塊類型,其相應分類與空間形態(tài)特征如表1與圖2所示。

    表1 基于MSPA模型的熱島斑塊分類及其空間形態(tài)特征Table 1 Classification of heat island patches and their spatial morphological characteristics based on MSPA model

    圖2 熱島斑塊類型及其空間形態(tài)特征示意圖Figure 2 Schematic diagram of the heat island patch types and their spatial morphological characteristics

    2.2.2 基于CIRCUIT理論的熱環(huán)境網絡識別

    在城鄉(xiāng)一體化加速以及建成區(qū)范圍不斷增加的雙重背景下,城市熱島斑塊也呈現(xiàn)不斷擴張趨勢(林中立等,2018;覃盟琳等,2020;喬治等,2022)。在擴張過程中,熱島斑塊與電路中的電荷類似,具有隨機游走特征(即隨機漫步理論)。CIRCUIT理論是借鑒隨機游走特征模擬相關“流體”在地表異質空間遷移和擴散運動的理論,主要用于相關廊道網絡,如生態(tài)廊道的識別(宋利利等,2016;李久林等,2020)。該理論基于歐姆定律,即電壓、電阻及電流的概念內核,可探尋源地間的最小成本路徑,將地表潛在的廊道路徑以電路網絡形式進行抽象表征,如圖3所示。在熱環(huán)境網絡中,電壓對應高低熱量源地間的熱力差,電阻對應源地間熱量流通受到的阻擋效應,電流對應源地間熱量的流通速率。本研究基于電路理論識別長江中游城市群熱環(huán)境網絡,具體來說包括3個部分:熱環(huán)境源地設置、阻力面設置及熱環(huán)境網絡識別。

    圖3 電路理論示意圖Figure 3 Schematic diagram of circuit theory

    (1)熱環(huán)境源地設置

    熱環(huán)境源地是熱環(huán)境網絡中具有一定規(guī)模的斑塊,其既是高熱量集聚空間,也是廊道交匯處的關鍵節(jié)點以及熱量傳輸?shù)摹爸修D站”,在夏季易對城鄉(xiāng)居民產生高溫威脅。本研究基于長江中游城市群熱島斑塊類型劃分結果,利用Arcgis 10.2軟件將面積大于4 km2的核心型斑塊定義為源地。

    (2)阻力面設置

    阻力面是電阻的空間形態(tài)表征,表示各熱環(huán)境源地間熱量擴散的受阻程度分布,反映地表熱量傳輸?shù)目蛇_性和趨勢。本研究選取背景以及各類型熱島斑塊,包括核心、孤島、孔隙、邊緣、環(huán)道、橋接以及支線作為阻力因子,并參照已有研究的賦值原則賦予其阻力值(喬治等,2023),以此構建長江中游城市群熱環(huán)境柵格阻力面(表2)。

    表2 各類型熱島斑塊阻力值Table 2 Resistance values for various types of heat island patches

    (3)熱環(huán)境網絡識別

    本研究通過使用Circuitscape 4.0軟件中的Build Network and Map Linkages模塊實現(xiàn)電路理論的空間可視化。首先設置初始電壓和初始電流分別為1V和1A,其次基于已設置的熱環(huán)境源地與阻力面數(shù)據(jù)的迭代遍歷運算,模擬仿真所有連接相鄰源地的熱環(huán)境廊道分布,進而對長江中游城市群熱環(huán)境網絡進行識別。

    2.2.3 熱環(huán)境網絡連通性評估

    熱環(huán)境源地與熱環(huán)境廊道間的連接程度稱為熱環(huán)境網絡連通性,其能夠反映源地間熱量的擴散過程與傳輸效能。本研究評估長江中游城市群熱環(huán)境網絡的總體連通性與空間連通性。

    (1)熱環(huán)境網絡總體連通性評估

    總體連通性反映網絡空間的總體復雜程度與連接度,映射節(jié)點與邊之間的物理和邏輯連接。在復雜網絡空間評價研究中,α指數(shù)、β指數(shù)以及γ指數(shù)是常用的定量評價拓撲指數(shù)(黃梅等,2022)。此類拓撲指數(shù)均是值越高表示熱環(huán)境網絡連通性越高,熱量傳輸效能越強,越不利于熱環(huán)境風險的降低,其不僅關注自身的網絡結構特征,更強調熱環(huán)境網絡的熱量擴散效率。因此,本研究選取上述3種指數(shù)對長江中游城市群熱環(huán)境網絡的總體連通性進行評估,其表達式為:

    式中:

    α——網絡閉合度,反映熱環(huán)境網絡中出現(xiàn)閉合回路的程度,其取值范圍為[0,1],當α=0時,表示熱環(huán)境網絡沒有回路,當α=1時,表示熱環(huán)境網絡存在最大可能數(shù)量的回路;β為網絡連接度,反映熱環(huán)境網絡中熱環(huán)境源地的平均連接數(shù),其取值范圍為[0,3],當β=0時,表示熱環(huán)境網絡不存在,當β<1時,表示熱環(huán)境網絡連接尚不完善,當β=1時,表示熱環(huán)境網絡形成單一回路,當β=2時,表示熱環(huán)境網絡連接較完善,當β=3時,表示熱環(huán)境網絡連接最完善;

    γ——網絡連通率,反映熱環(huán)境網絡中所有熱環(huán)境源地被連接的程度,其取值范圍為[0,1],當γ=0時,表示熱環(huán)境網絡中無廊道連接各源地,當γ=1時,表示熱環(huán)境網絡各源地彼此相連,熱環(huán)境廊道連接程度最高;

    N——熱環(huán)境源地(節(jié)點)數(shù)量;

    L——城市群熱環(huán)境廊道(邊)數(shù)量;

    Lmax——城市群中最大可能連接的熱環(huán)境廊道數(shù)量。

    (2)熱環(huán)境網絡空間連通性評估

    空間連通性反映網絡空間的局部區(qū)域在整個網絡空間中的重要程度,通過網絡空間中各節(jié)點和邊的影響力予以表征。在熱環(huán)境網絡中,空間連通性是衡量并識別關鍵熱環(huán)境源地及熱環(huán)境廊道的重要依據(jù)。本研究分別采取特征向量中心性(Eigenvector Centrality)與電流密度(Current Density)分別對各熱環(huán)境源地及熱環(huán)境廊道的影響程度進行測度(胡昊宇等,2022;喬治等,2023),以此來評估長江中游城市群熱環(huán)境網絡的空間連通性。

    ①特征向量中心性表示某熱環(huán)境源地在熱環(huán)境網絡中的影響力,其基本概念可概括為某熱環(huán)境源地的重要性不僅取決于自身的中心性,還與其鄰近熱環(huán)境源地的數(shù)量、中心性有關,與之相連的鄰居熱環(huán)境源地越重要,則該熱環(huán)境源地就越重要,其表達式為:

    式中:

    c(vi)——某熱環(huán)境源地vi的特征向量中心性,其取值范圍為[0,1],數(shù)值越高表示該熱環(huán)境源地在網絡中的影響力越大;

    Aij與λ——鄰接矩陣及對應特征向量值。本研究通過Gephi 0.10.0軟件的Fruchterman Reingold算法進行特征向量中心性測度及可視化。

    ②電流密度指通過單個像元的電流大小,在熱環(huán)境網絡中可表征某熱環(huán)境廊道的空間連通性,其值越高,代表該熱環(huán)境廊道的熱量傳輸效能越高。在并聯(lián)電路中,電流密度在電壓(熱力差)保持穩(wěn)定情形下,除受到電阻(阻力面)影響外,也受制于熱環(huán)境廊道的冗余度、寬度等因素,例如電流密度值會隨著熱環(huán)境廊道數(shù)量增加而降低(圖3)。本研究在Circuitscape 4.0軟件的Pinchpoint Mapper模塊中分別輸入熱環(huán)境源地及熱環(huán)境廊道空間分布數(shù)據(jù),通過選擇“All to One模式”(1個熱環(huán)境源地接地,其余熱環(huán)境源地均輸入1 V電壓與1 A電流)設置所有相鄰熱環(huán)境源地產生耦合作用,使所有連接相鄰熱環(huán)境源地的熱環(huán)境廊道路徑均產生仿真電流。

    3 結果與分析

    3.1 城市群地表溫度及熱島區(qū)域空間演變

    為探究2000-2020年長江中游城市群夏季地表溫度及熱島空間演變特征,依次選取2000、2010與2020年數(shù)據(jù)進行空間可視化,并依據(jù)地表溫度閾值區(qū)間進行等級劃分,得到由低到高5個地表溫度等級,包括低溫區(qū)(<21 ℃)、次低溫區(qū)(21-24 ℃)、常溫區(qū)(24-27 ℃)、次高溫區(qū)(27-30 ℃)、高溫區(qū)(>30 ℃),在此基礎上對各等級地表溫度面積占比進行統(tǒng)計。

    如圖4所示,2000年長江中游城市群高溫區(qū)(3%)主要呈斑塊狀分布在湖南長沙、衡陽以及江西吉安、撫州、鷹潭,少量分布在湖北武漢,次高溫區(qū)(25%)主要呈連片狀分布在湖南與江西絕大部分城市;2010年城市群高溫區(qū)(5%)主要集聚在江西,少量零散分布在湖北武漢以及湖南長沙、衡陽,次高溫區(qū)比例(32%)繼續(xù)增加,尤其在湖北地區(qū)增勢明顯;2020年城市群高溫區(qū)比例(4%)與2010年相比略有下降,在空間上集中分布在省會城市及其周邊,呈空間集聚態(tài)勢,次高溫區(qū)(36%)已擴散至城市群大部分城市。此外,2000-2020年期間,長江中游城市群常溫區(qū)、次低溫區(qū)以及低溫區(qū)比例總體呈下降趨勢,分別由2000年的56%、13%與2%下降至2020年的53%、6%與1%,其在空間上變化不明顯,主要分布在湖北宜昌、襄陽以及江西宜春??傮w而言,2000-2020年長江中游城市群夏季地表溫度的空間分布態(tài)勢逐漸由“南高北低”的兩級分化格局向以省會城市為核心的“核心-外圍”格局演變。

    圖4 2000-2020年長江中游城市群夏季地表溫度空間演變特征Figure 4 Spatial evolution characteristics of summer LST in urban agglomerations of the middle reaches of the Yangtze River from 2000 to 2020

    基于長江中游城市群地表建成區(qū)及其2 km緩沖區(qū)范圍與高溫區(qū)范圍提取熱島區(qū)域(圖5)。結果顯示,熱島區(qū)域總面積由2000年的2.80×103km2增加至2010年的6.41×103km2,后繼續(xù)增加至2020年的12.8×103km2,表明2000-2020年長江中游城市群熱島空間范圍隨建成區(qū)面積增加呈現(xiàn)持續(xù)擴張的趨勢,這與相關城市群研究結論趨同(林中立等,2018;覃盟琳等,2020)。

    圖5 2000-2020年長江中游城市群夏季熱島區(qū)域空間演變特征Figure 5 Spatial evolution characteristics of heat island area in urban agglomerations of the middle reaches of the Yangtze River from 2000 to 2020

    3.2 城市群熱島斑塊類型劃分

    如表3所示,基于MSPA模型依次識別并統(tǒng)計2000、2010與2020年長江中游城市群7種類型熱島斑塊面積及其占比,其中核心型、邊緣型熱島斑塊面積占比始終最高,其次為孤島型或支線型熱島斑塊,而孔隙型熱島斑塊面積占比始終最低。

    表3 2000-2020年長江中游城市群熱島斑塊面積及占比統(tǒng)計Table 3 Statistics on the area and proportion of heat island patches in urban agglomerations of the middle reaches of the Yangtze River from 2000 to 2020

    在2000-2010年期間,核心型熱島斑塊面積占比由31.1%增加至39.9%,增幅為8.80%,而孤島型、環(huán)道型、橋接型、支線型熱島斑塊總面積占比呈不同程度下降趨勢,由32.4%減少至23.0%,減幅為9.40%。伴隨快速城市化推進,長江中游城市群核心型熱島斑塊的面積增加包括兩方面原因,一是建成區(qū)范圍的擴大引發(fā)熱島區(qū)域面積增加,其內部核心型熱島斑塊區(qū)域也隨之擴充,二是分散、破碎的小型熱島斑塊逐漸匯聚、合并,形成了密集且均勻的核心型熱島斑塊(Du et al.,2016;Chen et al.,2021)。在2010-2020年期間,核心型熱島斑塊面積占比繼續(xù)增加至45.9%,增幅為6.00%,而其余類型熱島斑塊總面積占比呈下降趨勢。

    結合圖6可以看出,在2000-2020年期間,城市群核心型熱島斑塊在空間上逐漸向武漢、長沙以及南昌等地集聚。一方面,這些城市的人口不斷加劇,由此引發(fā)的城市活動日益增多,致使工業(yè)部門等功能片區(qū)排放出大量熱量及溫室氣體,影響大型熱島斑塊的產生。另一方面,這些城市的建成區(qū)空間不斷擴張,表征為自然環(huán)境逐漸被人工環(huán)境大規(guī)模侵占、可透性自然地表逐漸被硬質化人工地表替代,引發(fā)城市熱力與動力結構發(fā)生巨大轉變,迫使城市地表在日間儲存較多熱量,由此導致核心型熱島斑塊隨城市空間擴張而逐漸蔓延。

    圖6 2000-2020年長江中游城市群熱島斑塊空間分布Figure 6 Spatial distribution of heat island patches in urban agglomerations of the middle reaches of the Yangtze River from 2000 to 2020

    如圖7所示,依次統(tǒng)計2000、2010與2020年長江中游城市群7種類型熱島斑塊的地表溫度值。結果顯示,長江中游城市群各類型熱島斑塊的地表溫度存在差異性且隨時間推移無明顯變化。就地表溫度中位數(shù)而言,在2000-2020年期間,核心型熱島斑塊始終最高,而支線型熱島斑塊始終最低,其差值在2000、2010與2020年依次為1.50、1.16、1.30 ℃,表明長江中游城市群核心型熱島斑塊的增溫聚熱趨勢顯著。

    圖7 2000-2020年長江中游城市群熱島斑塊地表溫度箱線圖Figure 7 Box-plot of LST of heat island patches in urban agglomerations of the middle reaches of the Yangtze River from 2000 to 2020

    3.3 城市群熱環(huán)境網絡空間演變

    通過篩選2000-2020年長江中游城市群面積大于4 km2的核心型熱島斑塊,對熱環(huán)境源地進行提取。如表4和圖8所示,源地數(shù)量在2000、2010以及2020年分別為56、116、215個,隨時間推移呈不斷增加趨勢。在空間上,2000-2020年期間長江中游城市群熱環(huán)境源地逐漸向湖北武漢、湖南長沙以及江西南昌3座省會城市聚集。

    表4 2000-2020年長江中游城市群熱環(huán)境源地與熱環(huán)境廊道數(shù)量統(tǒng)計Table 4 Statistics on the number of thermal environment sources and corridors in urban agglomerations of the middle reaches of the Yangtze River from 2000 to 2020

    圖8 2000-2020年長江中游城市群熱環(huán)境網絡空間分布Figure 8 Spatial distribution of thermal environment network in urban agglomerations of the middle reaches of the Yangtze River from 2000 to 2020

    依據(jù)熱環(huán)境源地分布以及電路理論,構建2000-2020年長江中游城市群熱環(huán)境網絡。結果顯示,2000年廊道數(shù)量為89條,在城市群中有77.4%的城市覆蓋廊道分布,其中南部廊道較為密集,而北部廊道空間分布不均且較為稀疏;2010年廊道數(shù)量為192條,較2000年有所增加,在城市群中有83.9%的城市覆蓋廊道分布,熱環(huán)境網絡總體分布密集;2020年廊道數(shù)量繼續(xù)增至378條,連接區(qū)域向湖北擴散,在城市群中有100%的城市覆蓋廊道分布,表明已初步形成城市群全局熱環(huán)境網絡。

    3.4 城市群熱環(huán)境網絡連通性評估

    3.4.1 總體連通性評估

    根據(jù)長江中游城市群熱環(huán)境源地數(shù)量(N)及熱環(huán)境廊道數(shù)量(L)計算α、β、γ指數(shù),以評估城市群熱環(huán)境網絡總體連通性。如圖9所示,α指數(shù)在2000、2010與2020年分別為0.318、0.339、0.386,表明長江中游城市群閉合程度逐漸增強;β指數(shù)在2000、2010以及2020年分別為1.59、1.66、1.76,表明長江中游城市群熱環(huán)境源地的平均連接數(shù)逐漸增加,熱環(huán)境網絡已形成單一回路;γ指數(shù)在2000、2010與2020年分別為0.549、0.561、0.592,反映長江中游城市群熱環(huán)境源地的彼此相連程度逐漸提高,在2020年存在約60%的熱環(huán)境源地通過廊道彼此相連??傮w來說,α、β、γ指數(shù)隨時間推移不斷增加,長江中游城市群熱環(huán)境網絡總體連通性不斷增強,源地間熱量的傳輸效能逐步提高,對城鄉(xiāng)居民的高溫威脅程度呈上升態(tài)勢。主要原因有兩方面,一是由于城市的擴張效應與集聚效應,分別引發(fā)熱環(huán)境源地面積與數(shù)量增加,導致城市間熱量傳遞的“中轉站”增多;二是由于公路、鐵路等交通路網作為熱環(huán)境廊道的外在表征之一,在城市化建設中不斷得到延伸與完善,由此不僅改變了城市群自然下墊面類型,也為城市間熱量傳遞提供了更便捷的途徑。

    圖9 2000-2020年長江中游城市群熱環(huán)境網絡總體連通性Figure 9 Overall connectivity of the thermal environment network in urban agglomerations of the middle reaches of the Yangtze River from 2000 to 2020

    3.4.2 空間連通性評估

    本研究分別對長江中游城市群熱環(huán)境源地的特征向量中心性及熱環(huán)境廊道的電流密度進行測度,以此來評估熱環(huán)境網絡的空間連通性。

    首先,根據(jù)2000、2010與2020年熱環(huán)境源地數(shù)量分別設置56、116、215處編號,通過Fruchterman Reingold算法對各年份熱環(huán)境源地特征向量中心性測度及可視化(圖10)。結果顯示,在2000年,熱環(huán)境源地的特征向量中心性值整體較低,僅有2處高于0.8,分別位于湖南長沙、江西南昌,表明這兩處熱環(huán)境源地在2000年長江中游城市群熱環(huán)境網絡中影響力最強;在2010年,熱環(huán)境源地特征向量中心性整體得到提升,其值高于0.8的熱環(huán)境源地增至3處,分別位于湖北咸寧、湖南長沙、江西南昌;在2020年,熱環(huán)境源地特征向量中心性提升顯著,其值高于0.8的熱環(huán)境源地增至6處,其中湖北有3處(武漢1處,仙桃1處,咸寧1處),湖南有2處(長沙1處,岳陽1處),江西有1處(南昌1處)。

    圖10 2000-2020年長江中游城市群熱環(huán)境源地空間連通性Figure 10 Spatial connectivity of thermal environment sources in urban agglomerations of the middle reaches of the Yangtze River from 2000 to 2020

    其次,利用Circuitscape 4.0軟件中Pinchpoint Mapper模塊對2000年、2010年與2020年熱環(huán)境廊道的電流密度進行測度與空間可視化(圖11)。結果顯示,在2000年,江西地區(qū)熱環(huán)境廊道空間連通性較強,其平均電流密度值為0.139A,熱量傳輸效能較高,此階段城市群城市化水平較低,江西地區(qū)由于熱環(huán)境源地分布較為集中且熱環(huán)境廊道大多受自然環(huán)境隔斷,無法形成遠距離熱擴散效應,僅與臨近區(qū)域保持較高連通性。相關證據(jù)表明,源地分布越密集、距離越近,其連通性越高(劉婷等,2023)。湖北與湖南地區(qū)熱環(huán)境廊道的平均電流密度值分別為0.118A,0.127A;在2010年,電流密度高值區(qū)域逐漸向北偏移,湖北地區(qū)熱環(huán)境廊道的平均電流密度值增至0.165A,而湖南與江西地區(qū)熱環(huán)境廊道的平均電流密度值分別增至0.162A、0.153A;在2020年,湖北、湖南與江西地區(qū)熱環(huán)境廊道平均電流密度值繼續(xù)呈現(xiàn)不同程度增加趨勢,分別達到0.231A、0.202A、0.183A。綜上,伴隨快速城市化推進,在2000-2020年期間長江中游城市群熱環(huán)境網絡空間連通性逐漸增強,其中湖北地區(qū)增勢最為顯著,其總體分布趨勢由“南高北低”的空間異質性格局向“全局高連通性”的空間穩(wěn)定性格局演變。

    圖11 2000-2020年長江中游城市群熱環(huán)境廊道空間連通性Figure 11 Spatial connectivity of thermal environment corridors in urban agglomerations of the middle reaches of the Yangtze River from 2000 to 2020

    4 結論與展望

    4.1 結論

    本研究以長江中游城市群為例,基于2000、2010及2020年多源面板數(shù)據(jù)從時空動態(tài)視角探索城市群熱環(huán)境網絡的時空演變特征并評估其總體連通性與空間連通性,為長江中游城市群氣候適應發(fā)展策略提供啟示。其主要結論如下:

    (1)2000-2020年期間,長江中游城市群夏季地表溫度高溫區(qū)與次高溫區(qū)比例整體呈上升趨勢,常溫區(qū)、低溫區(qū)以及次低溫區(qū)比例整體呈下降趨勢;就各類型熱島斑塊的地表溫度中位數(shù)而言,呈現(xiàn)核心型熱島斑塊最高而支線型熱島斑塊最低的規(guī)律。

    (2)2000-2020年期間,長江中游城市群夏季熱島區(qū)域總面積由2.80×103km2增至12.8×103km2,呈現(xiàn)持續(xù)擴張趨勢,其中核心型熱島斑塊面積占比由31.1%增加至45.9%,在空間上逐漸向以武漢、長沙以及南昌等區(qū)域集聚。

    (3)2000-2020年期間,長江中游城市群熱環(huán)境源地數(shù)量由56個增至215個,熱環(huán)境廊道數(shù)量由89條增至378條,包含熱環(huán)境廊道的城市比例由77.4%增至100%,已初步形成城市群全局熱環(huán)境網絡。

    (4)2000-2020年期間,長江中游城市群熱環(huán)境網絡總體連通性逐漸增強,熱環(huán)境源地間的熱量傳輸效能呈上升態(tài)勢;熱環(huán)境網絡空間連通性隨時間推移變化顯著,其總體分布趨勢由“南高北低”的空間異質性格局向“全局高連通性”的空間穩(wěn)定性格局演變。

    4.2 展望

    在城市群尺度動態(tài)識別熱環(huán)境網絡并評估其連通性,能夠為城市熱環(huán)境提供新的研究視角,有助于在快速城市化進程中為降低城鄉(xiāng)居民熱環(huán)境風險提供新的思路。然而,本研究是利用多種網絡分析方法對熱環(huán)境空間結構內在成因機制分析的一次嘗試,在熱環(huán)境研究領域仍屬于探索性工作,后續(xù)將展開三方面完善工作:首先,本研究以熱環(huán)境源地及熱環(huán)境廊道為基礎展開熱環(huán)境網絡研究,其與相關降溫網絡(如生態(tài)網絡)之間在空間布局、土地類型、人口密度等方面存在何種聯(lián)系是下一步研究重點,這有助于為熱環(huán)境優(yōu)化策略提供建設指引;其次,本研究基于電路理論對熱環(huán)境網絡進行識別,其識別結果與其他路徑識別方法(如最小成本路徑、神經網絡等)得到的結論是否存在差異需要得到對比與驗證;最后,本研究將結合相關人口健康屬性數(shù)據(jù),針對城鄉(xiāng)居民健康水平與熱環(huán)境網絡連通性的關聯(lián)性進行測度,以增強熱環(huán)境網絡評價的實踐意義。

    猜你喜歡
    環(huán)境
    長期鍛煉創(chuàng)造體內抑癌環(huán)境
    一種用于自主學習的虛擬仿真環(huán)境
    孕期遠離容易致畸的環(huán)境
    不能改變環(huán)境,那就改變心境
    環(huán)境與保護
    環(huán)境
    孕期遠離容易致畸的環(huán)境
    高等院校環(huán)境類公選課的實踐和探討
    掌握“三個三” 兜底環(huán)境信訪百分百
    我國環(huán)境會計初探
    中國商論(2016年33期)2016-03-01 01:59:38
    在线看三级毛片| 俺也久久电影网| 麻豆国产av国片精品| 欧美高清性xxxxhd video| 在线观看舔阴道视频| 人妻夜夜爽99麻豆av| 少妇丰满av| 极品教师在线免费播放| АⅤ资源中文在线天堂| 日韩高清综合在线| 毛片女人毛片| 特大巨黑吊av在线直播| 日韩精品中文字幕看吧| 久久亚洲精品不卡| 色播亚洲综合网| 蜜桃亚洲精品一区二区三区| 中文资源天堂在线| 日日啪夜夜撸| 91久久精品国产一区二区三区| 给我免费播放毛片高清在线观看| 黄色丝袜av网址大全| 欧美人与善性xxx| 日韩一本色道免费dvd| 一a级毛片在线观看| 欧美一区二区亚洲| 中出人妻视频一区二区| 看免费成人av毛片| 欧美日韩精品成人综合77777| 丰满乱子伦码专区| 99精品久久久久人妻精品| 99热这里只有是精品50| 亚洲一区二区三区色噜噜| 一进一出好大好爽视频| 嫁个100分男人电影在线观看| 九九爱精品视频在线观看| 看免费成人av毛片| 久久久成人免费电影| 少妇猛男粗大的猛烈进出视频 | 免费大片18禁| 久久久国产成人精品二区| 色综合婷婷激情| 国产爱豆传媒在线观看| 九九热线精品视视频播放| 99在线视频只有这里精品首页| 亚洲人成网站高清观看| 精品一区二区三区人妻视频| 国产av不卡久久| 国产亚洲精品久久久com| 一本精品99久久精品77| 九九在线视频观看精品| 高清日韩中文字幕在线| 好男人在线观看高清免费视频| 一本精品99久久精品77| 国产高清视频在线播放一区| 淫秽高清视频在线观看| 亚洲一区高清亚洲精品| 黄色一级大片看看| 91麻豆精品激情在线观看国产| 精品久久久久久久末码| 99热这里只有是精品50| 91午夜精品亚洲一区二区三区 | 老司机深夜福利视频在线观看| 91午夜精品亚洲一区二区三区 | 俺也久久电影网| 最好的美女福利视频网| 神马国产精品三级电影在线观看| 婷婷精品国产亚洲av在线| 国产精品亚洲一级av第二区| 免费av观看视频| 欧美xxxx黑人xx丫x性爽| 九色成人免费人妻av| 亚洲av美国av| 99热网站在线观看| 欧美黑人欧美精品刺激| 国产精品国产三级国产av玫瑰| 精品午夜福利在线看| 国产精品一区二区三区四区久久| 丰满乱子伦码专区| 桃红色精品国产亚洲av| 午夜福利视频1000在线观看| 久久久久久九九精品二区国产| 小说图片视频综合网站| 九九热线精品视视频播放| 91久久精品国产一区二区三区| 美女cb高潮喷水在线观看| 好男人在线观看高清免费视频| 一卡2卡三卡四卡精品乱码亚洲| av专区在线播放| 黄色日韩在线| 大型黄色视频在线免费观看| 丰满的人妻完整版| 成人无遮挡网站| 亚洲av不卡在线观看| 最近最新中文字幕大全电影3| 中国美白少妇内射xxxbb| 亚洲成人免费电影在线观看| www日本黄色视频网| 亚洲美女视频黄频| 精品午夜福利在线看| 久久久成人免费电影| 亚洲精品粉嫩美女一区| 99久久中文字幕三级久久日本| 国产精品国产三级国产av玫瑰| 欧美性感艳星| 久久精品影院6| eeuss影院久久| 午夜福利欧美成人| 午夜影院日韩av| eeuss影院久久| 国产精品久久久久久精品电影| 婷婷色综合大香蕉| 国产精品无大码| 亚洲美女视频黄频| 香蕉av资源在线| 18禁在线播放成人免费| 成人无遮挡网站| 日本一本二区三区精品| 国产男靠女视频免费网站| 精品免费久久久久久久清纯| 亚洲天堂国产精品一区在线| 久久精品人妻少妇| 一夜夜www| 国产v大片淫在线免费观看| 看片在线看免费视频| 亚洲av一区综合| 日本免费a在线| 国产亚洲av嫩草精品影院| 成人av在线播放网站| 又爽又黄无遮挡网站| 欧洲精品卡2卡3卡4卡5卡区| 尾随美女入室| 成人毛片a级毛片在线播放| 日日撸夜夜添| 一夜夜www| 国产人妻一区二区三区在| 亚洲精品456在线播放app | 波多野结衣高清作品| 色哟哟·www| 色播亚洲综合网| 尤物成人国产欧美一区二区三区| 人人妻人人看人人澡| 国产精品一区二区性色av| 国内揄拍国产精品人妻在线| 最近视频中文字幕2019在线8| 三级毛片av免费| 免费av观看视频| 免费看a级黄色片| 少妇裸体淫交视频免费看高清| 99热6这里只有精品| 韩国av一区二区三区四区| 精品人妻1区二区| 嫩草影院精品99| 亚洲av免费高清在线观看| 中文字幕人妻熟人妻熟丝袜美| 欧美激情国产日韩精品一区| 亚洲自拍偷在线| 久久精品国产自在天天线| 国产成年人精品一区二区| 日韩欧美精品免费久久| 九色国产91popny在线| 国产成人aa在线观看| 天堂影院成人在线观看| 夜夜看夜夜爽夜夜摸| 亚洲精品在线观看二区| 亚洲精品色激情综合| 国内精品宾馆在线| 中文字幕免费在线视频6| 国产伦精品一区二区三区四那| 人人妻人人看人人澡| 国产精品98久久久久久宅男小说| av中文乱码字幕在线| 男女边吃奶边做爰视频| 天天躁日日操中文字幕| 国产单亲对白刺激| 亚洲欧美日韩东京热| 午夜影院日韩av| 久久精品影院6| 色精品久久人妻99蜜桃| 精品人妻视频免费看| 国产精品久久久久久久电影| 亚洲中文字幕一区二区三区有码在线看| 最好的美女福利视频网| 亚洲内射少妇av| 精品一区二区三区视频在线| 国产色婷婷99| 日本黄色视频三级网站网址| 亚洲成人中文字幕在线播放| 人人妻,人人澡人人爽秒播| 成人特级黄色片久久久久久久| 男人舔奶头视频| 少妇的逼好多水| 天堂√8在线中文| 午夜福利视频1000在线观看| 欧美另类亚洲清纯唯美| 国产精品亚洲美女久久久| 国产精品1区2区在线观看.| 久久国产乱子免费精品| 欧美一区二区亚洲| 香蕉av资源在线| 久久久久久久久大av| 久久午夜亚洲精品久久| 国产色爽女视频免费观看| 国产精品一区二区性色av| 亚洲18禁久久av| 1024手机看黄色片| 日韩欧美在线乱码| 国产男靠女视频免费网站| 精品久久国产蜜桃| 99热6这里只有精品| 免费在线观看日本一区| 不卡一级毛片| 波多野结衣高清作品| 亚洲欧美日韩高清专用| 精品国内亚洲2022精品成人| 午夜福利高清视频| 综合色av麻豆| 99在线视频只有这里精品首页| 国产综合懂色| 在线a可以看的网站| 日韩欧美一区二区三区在线观看| 欧美成人一区二区免费高清观看| 欧美高清性xxxxhd video| 国产精品乱码一区二三区的特点| 国产午夜精品久久久久久一区二区三区 | 久久久久久伊人网av| 自拍偷自拍亚洲精品老妇| 97超视频在线观看视频| 午夜福利视频1000在线观看| 成年免费大片在线观看| 99精品在免费线老司机午夜| 日本免费一区二区三区高清不卡| 国产精品人妻久久久影院| 国产精品久久久久久精品电影| 久久久久久久午夜电影| 亚洲精品影视一区二区三区av| 国产美女午夜福利| 国产高清三级在线| 999久久久精品免费观看国产| 我的女老师完整版在线观看| 我要搜黄色片| 51国产日韩欧美| 国产大屁股一区二区在线视频| 波野结衣二区三区在线| 22中文网久久字幕| 精品久久久久久久久av| 欧美成人一区二区免费高清观看| 校园春色视频在线观看| 搞女人的毛片| 村上凉子中文字幕在线| 欧美精品啪啪一区二区三区| 成人国产麻豆网| 国模一区二区三区四区视频| 久久人妻av系列| 内射极品少妇av片p| 97人妻精品一区二区三区麻豆| 国产精品久久久久久av不卡| 国产 一区 欧美 日韩| 亚洲色图av天堂| 亚洲真实伦在线观看| 美女黄网站色视频| 亚洲内射少妇av| 非洲黑人性xxxx精品又粗又长| 又黄又爽又免费观看的视频| 九九久久精品国产亚洲av麻豆| 国产男靠女视频免费网站| 长腿黑丝高跟| 99热6这里只有精品| 国产伦一二天堂av在线观看| 男女下面进入的视频免费午夜| 欧美成人a在线观看| 亚洲av中文字字幕乱码综合| 亚洲精品日韩av片在线观看| 在线观看美女被高潮喷水网站| 欧美另类亚洲清纯唯美| 免费看光身美女| 欧美bdsm另类| 一个人免费在线观看电影| 99精品久久久久人妻精品| 赤兔流量卡办理| 亚洲精品一区av在线观看| 女人十人毛片免费观看3o分钟| 国产精品久久电影中文字幕| 亚洲成a人片在线一区二区| 一级a爱片免费观看的视频| 午夜亚洲福利在线播放| 老司机午夜福利在线观看视频| 久久国产精品人妻蜜桃| 亚洲 国产 在线| 午夜亚洲福利在线播放| 欧美另类亚洲清纯唯美| 在线天堂最新版资源| 欧美日韩精品成人综合77777| 免费看美女性在线毛片视频| 亚洲国产色片| 亚洲男人的天堂狠狠| 99热精品在线国产| av在线亚洲专区| 日韩人妻高清精品专区| 久久精品国产亚洲av香蕉五月| 久久九九热精品免费| 久久久精品大字幕| 一进一出抽搐gif免费好疼| 亚洲国产欧洲综合997久久,| 久99久视频精品免费| 色综合婷婷激情| 成人三级黄色视频| 深爱激情五月婷婷| 亚洲人与动物交配视频| 日本熟妇午夜| 国产精品久久久久久精品电影| 国产美女午夜福利| 国产午夜精品论理片| 久久99热这里只有精品18| 亚洲精品久久国产高清桃花| 精品一区二区三区av网在线观看| 男人舔奶头视频| 99九九线精品视频在线观看视频| 中文字幕人妻熟人妻熟丝袜美| 国产男靠女视频免费网站| 精品久久久久久久末码| 欧美高清成人免费视频www| 少妇被粗大猛烈的视频| 男女之事视频高清在线观看| 亚洲avbb在线观看| 欧美人与善性xxx| 国产伦精品一区二区三区视频9| 欧美日韩中文字幕国产精品一区二区三区| 国产精品精品国产色婷婷| 精品欧美国产一区二区三| 日韩欧美在线二视频| 国产毛片a区久久久久| 一进一出好大好爽视频| 亚洲黑人精品在线| 精品久久久久久成人av| 国产精品1区2区在线观看.| 九色成人免费人妻av| 91精品国产九色| 国产一级毛片七仙女欲春2| 极品教师在线免费播放| 九色国产91popny在线| 国产一区二区三区av在线 | 国产视频一区二区在线看| 国产一区二区三区在线臀色熟女| 欧美极品一区二区三区四区| 日本成人三级电影网站| 欧美日韩中文字幕国产精品一区二区三区| 成熟少妇高潮喷水视频| 国产高清视频在线播放一区| а√天堂www在线а√下载| 男女那种视频在线观看| av天堂在线播放| 久久久久久久精品吃奶| 欧美潮喷喷水| 国产不卡一卡二| 乱系列少妇在线播放| 国产成人一区二区在线| 国产精品一及| 制服丝袜大香蕉在线| 最近在线观看免费完整版| 精品一区二区三区av网在线观看| 久久久久国产精品人妻aⅴ院| 亚洲精品在线观看二区| 色综合站精品国产| 免费大片18禁| 国产大屁股一区二区在线视频| 国产黄片美女视频| 国产精品av视频在线免费观看| 国产真实乱freesex| 伦理电影大哥的女人| 十八禁网站免费在线| 日本撒尿小便嘘嘘汇集6| 小说图片视频综合网站| 婷婷六月久久综合丁香| 日韩欧美精品免费久久| 日本一本二区三区精品| 国产三级中文精品| 有码 亚洲区| 久久久国产成人精品二区| 欧美日韩瑟瑟在线播放| 波多野结衣高清作品| 干丝袜人妻中文字幕| 亚洲中文字幕一区二区三区有码在线看| 制服丝袜大香蕉在线| 欧美潮喷喷水| 国产精品久久久久久精品电影| 成人av在线播放网站| 国产精品98久久久久久宅男小说| 日韩强制内射视频| 国产综合懂色| 美女免费视频网站| 无人区码免费观看不卡| 欧美最新免费一区二区三区| 久久久久久大精品| 级片在线观看| 成人三级黄色视频| 国产美女午夜福利| 日韩中文字幕欧美一区二区| 国国产精品蜜臀av免费| 精品久久久久久久久久久久久| 欧美+亚洲+日韩+国产| 嫩草影院精品99| 免费看日本二区| 一个人看的www免费观看视频| 久久亚洲精品不卡| www.www免费av| 午夜免费男女啪啪视频观看 | 亚洲精品亚洲一区二区| 日本 欧美在线| 亚洲性久久影院| 亚州av有码| 久久精品国产99精品国产亚洲性色| 1000部很黄的大片| 色精品久久人妻99蜜桃| 国产午夜精品论理片| 国产av不卡久久| 天堂影院成人在线观看| 91av网一区二区| 成人性生交大片免费视频hd| 国产高清不卡午夜福利| 观看美女的网站| 亚洲av中文字字幕乱码综合| 在线免费观看的www视频| 欧美另类亚洲清纯唯美| 免费看a级黄色片| 在线免费十八禁| 全区人妻精品视频| 亚洲av熟女| 搡老熟女国产l中国老女人| 我的女老师完整版在线观看| 久久久国产成人精品二区| 欧美日韩瑟瑟在线播放| 国产精品国产三级国产av玫瑰| 99久久精品国产国产毛片| 成人永久免费在线观看视频| 国产一区二区在线av高清观看| 热99re8久久精品国产| 国产午夜精品论理片| 久久这里只有精品中国| 国产成人av教育| 亚洲成人免费电影在线观看| 啪啪无遮挡十八禁网站| 日日摸夜夜添夜夜添av毛片 | 蜜桃久久精品国产亚洲av| 亚洲自拍偷在线| 国内精品久久久久久久电影| 国产aⅴ精品一区二区三区波| 亚洲精品一区av在线观看| 午夜精品在线福利| 亚洲最大成人中文| 欧美精品国产亚洲| 非洲黑人性xxxx精品又粗又长| 国产激情偷乱视频一区二区| 美女高潮喷水抽搐中文字幕| 亚洲avbb在线观看| 成人av一区二区三区在线看| 一边摸一边抽搐一进一小说| 日韩精品中文字幕看吧| 亚洲无线在线观看| 精品久久久久久成人av| 人人妻,人人澡人人爽秒播| 欧美成人a在线观看| 日韩国内少妇激情av| 国产精品女同一区二区软件 | 亚洲精品色激情综合| 我要搜黄色片| 免费一级毛片在线播放高清视频| 亚洲av美国av| 99久久成人亚洲精品观看| 2021天堂中文幕一二区在线观| 日韩亚洲欧美综合| 全区人妻精品视频| 国产精品久久视频播放| 国产精品乱码一区二三区的特点| 极品教师在线免费播放| 婷婷精品国产亚洲av| 好男人在线观看高清免费视频| 日日啪夜夜撸| 日本三级黄在线观看| 日本精品一区二区三区蜜桃| 亚洲精华国产精华液的使用体验 | 久久久久久伊人网av| 黄色丝袜av网址大全| 国内揄拍国产精品人妻在线| 成人鲁丝片一二三区免费| av在线蜜桃| 亚洲aⅴ乱码一区二区在线播放| 真人做人爱边吃奶动态| 美女黄网站色视频| 亚洲欧美精品综合久久99| 97碰自拍视频| 亚洲av免费高清在线观看| 国产乱人伦免费视频| 听说在线观看完整版免费高清| 午夜精品一区二区三区免费看| 久久精品久久久久久噜噜老黄 | 色播亚洲综合网| 欧美成人免费av一区二区三区| 麻豆久久精品国产亚洲av| 欧美在线一区亚洲| 国产精品不卡视频一区二区| 十八禁网站免费在线| 很黄的视频免费| 国内精品一区二区在线观看| 深夜精品福利| 99精品久久久久人妻精品| 国产精品不卡视频一区二区| 国产一区二区亚洲精品在线观看| 99久久成人亚洲精品观看| 搡女人真爽免费视频火全软件 | 听说在线观看完整版免费高清| 真人做人爱边吃奶动态| 成年女人永久免费观看视频| 免费在线观看日本一区| 国产91精品成人一区二区三区| 午夜久久久久精精品| 精品久久国产蜜桃| 成年版毛片免费区| 国产精品一区www在线观看 | 亚洲在线观看片| 久久精品影院6| 国产av麻豆久久久久久久| 亚洲熟妇熟女久久| 国产又黄又爽又无遮挡在线| 色综合亚洲欧美另类图片| 一级av片app| 我的老师免费观看完整版| 最后的刺客免费高清国语| 不卡视频在线观看欧美| 我要看日韩黄色一级片| 麻豆一二三区av精品| 色吧在线观看| 国产大屁股一区二区在线视频| 婷婷六月久久综合丁香| 亚洲精品乱码久久久v下载方式| 久久精品夜夜夜夜夜久久蜜豆| 99热这里只有精品一区| 亚洲精品日韩av片在线观看| 国产单亲对白刺激| 一级a爱片免费观看的视频| 女的被弄到高潮叫床怎么办 | 在线免费十八禁| 日本-黄色视频高清免费观看| 日本欧美国产在线视频| 一区二区三区激情视频| 91久久精品国产一区二区成人| 欧美精品啪啪一区二区三区| 在线观看美女被高潮喷水网站| 18+在线观看网站| 久久久精品欧美日韩精品| 人妻丰满熟妇av一区二区三区| 免费不卡的大黄色大毛片视频在线观看 | 国产午夜精品论理片| av在线天堂中文字幕| 国产亚洲精品综合一区在线观看| 日韩强制内射视频| 韩国av一区二区三区四区| 99九九线精品视频在线观看视频| 色5月婷婷丁香| 搡老妇女老女人老熟妇| 亚洲欧美日韩东京热| 国产精品久久久久久av不卡| 国产高清不卡午夜福利| 毛片一级片免费看久久久久 | 精品人妻视频免费看| 97超级碰碰碰精品色视频在线观看| 永久网站在线| 成人性生交大片免费视频hd| 一进一出抽搐gif免费好疼| av女优亚洲男人天堂| 国产精品综合久久久久久久免费| 狂野欧美白嫩少妇大欣赏| 国产高清视频在线播放一区| 国内精品久久久久久久电影| 亚洲av中文av极速乱 | 日本黄色视频三级网站网址| 日韩亚洲欧美综合| 国产爱豆传媒在线观看| 一进一出抽搐动态| 又粗又爽又猛毛片免费看| 麻豆精品久久久久久蜜桃| 五月伊人婷婷丁香| 直男gayav资源| 美女高潮喷水抽搐中文字幕| 联通29元200g的流量卡| 九九热线精品视视频播放| 欧美人与善性xxx| 久久午夜福利片| 欧美性猛交黑人性爽| a级一级毛片免费在线观看| 听说在线观看完整版免费高清| 国产精品乱码一区二三区的特点| 日本 欧美在线| 国产熟女欧美一区二区| 久9热在线精品视频| 亚洲国产高清在线一区二区三| 精品欧美国产一区二区三| 中文字幕av成人在线电影| 美女高潮的动态| 美女大奶头视频| 黄色丝袜av网址大全| 少妇的逼好多水| 在线播放无遮挡| 天美传媒精品一区二区| 亚洲va日本ⅴa欧美va伊人久久| 日日干狠狠操夜夜爽| 波野结衣二区三区在线| 久久久久九九精品影院| 国产精品一区二区性色av| 韩国av在线不卡| 黄片wwwwww| 老女人水多毛片| 久久久久久国产a免费观看| 欧美激情国产日韩精品一区| 色精品久久人妻99蜜桃| 两个人的视频大全免费|