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

    內(nèi)蒙古地表凍融指數(shù)動態(tài)變化與驅(qū)動因素分析

    2022-02-03 07:24:22張昊琛薩楚拉孟凡浩王牧蘭高紅豆ADIYASaruulzaya
    干旱區(qū)研究 2022年6期
    關(guān)鍵詞:多年凍土凍土凍融

    張昊琛, 薩楚拉, 孟凡浩, 羅 敏, 王牧蘭,高紅豆, ADIYA Saruulzaya

    (1.內(nèi)蒙古師范大學(xué)地理科學(xué)學(xué)院,內(nèi)蒙古 呼和浩特 010022;2.內(nèi)蒙古自治區(qū)遙感與地理信息系統(tǒng)重點實驗室,內(nèi)蒙古 呼和浩特 010022;3.蒙古科學(xué)院地理與地球生態(tài)研究所,蒙古 烏蘭巴托 15170)

    凍土一般是指溫度在0 ℃或0 ℃以下,并含有冰的各種巖土和土壤,包括永久凍土(PF)和季節(jié)性凍土(SFG)。我國是世界凍土大區(qū),多年凍土和季節(jié)凍土分別占國土面積的22%和53.5%[1]。IPCC 第5次評估顯示,全球地表溫度在1880—2012年期間升高了0.85 ℃[2]。全球變暖在高緯度和高海拔區(qū)域加劇,對凍土的形成有較大的影響,尤其是在我國北方高緯度,溫度升高幅度大于全國平均升溫水平[3-5]。隨著氣候變化,我國凍土大致體現(xiàn)在:最大凍土深度下降、凍結(jié)日期推遲、消融日期提前、凍土下限上升[6-7]。20 世紀(jì)70 年代至21 世紀(jì)初,我國凍土面積減少約18.6%,面積由2.15×106km2減少到1.75×106km2[8]。按照對未來的天氣情況預(yù)報,未來冰凍圈各組成部分面積將持續(xù)減少,而多年凍土的活動層厚度將會進一步增加,季節(jié)凍土的凍結(jié)深度也將進一步下降[9-10]。

    凍融指數(shù)是研究凍土特性的重要參數(shù),也是判定其存在和長期變化的一個重要指標(biāo)[11]。凍土研究中凍融指數(shù)有助于繪制和預(yù)測永久凍土分布,估算土壤凍結(jié)深度,并計算活躍層厚度,凍結(jié)指數(shù)的下降和融化指數(shù)的上升在一定程度上說明凍土環(huán)境的退化[12]。利用季節(jié)凍土區(qū)的凍結(jié)指數(shù),可以預(yù)測季節(jié)凍融深度[13];利用多年凍土區(qū)的凍融指數(shù),不僅能夠繪制出凍土的分布圖,還能對多年凍土的活動層厚度進行預(yù)測[14-15];Peng 等[12]的研究結(jié)果顯示,北半球高緯度和高海拔地區(qū)的凍結(jié)指數(shù)下降幅度最大,凍結(jié)指數(shù)的下降速率要比融化指數(shù)的上升速率快得多;Shi 等[16]研究表明,從1901—2015 年,環(huán)北極地區(qū)的凍結(jié)和融化指數(shù)分別呈現(xiàn)降低、上升趨勢,在1988年之后,其變化趨勢更為明顯;Liu等[17]基于地面溫度觀測資料和再分析產(chǎn)品計算了雅魯藏布江上游凍融指數(shù),并基于凍融指數(shù)模擬了凍土面積的變化,效果良好;Wu等[18]對青藏高原凍融指數(shù)的變化規(guī)律進行了分析,發(fā)現(xiàn)多年凍土區(qū)的凍結(jié)指數(shù)比季節(jié)凍土區(qū)的凍結(jié)指數(shù)降低速度要快,而融化指數(shù)上升速率更慢;Wang 等[19]研究認(rèn)為,黃河源區(qū)海拔高度對凍融指數(shù)的空間分布有一定的影響,地表凍融指數(shù)的變化比空氣凍融指數(shù)的變化要大。

    內(nèi)蒙古地區(qū)的氣象觀測站數(shù)量較少,有些站點建設(shè)時間較晚,總體分布稀疏,個別站點存在數(shù)據(jù)缺失問題,氣候變化條件下的凍融指數(shù)變化研究較為欠缺。僅有的研究主要聚焦在典型草原或多年凍土區(qū),且研究對象為空氣凍融指數(shù)[20-21],不能完全代表內(nèi)蒙古全區(qū)的地表凍融特征,并且研究區(qū)地表凍融指數(shù)變化的驅(qū)動因素有待探究。鑒于此,本文基于內(nèi)蒙古氣象站1980—2019 年日均地表溫度數(shù)據(jù)、結(jié)合中國第一代全球陸面再分析產(chǎn)品(CRA)數(shù)據(jù)以及NDVI 數(shù)據(jù),對內(nèi)蒙古近40 a 地表凍融指數(shù)時空變化特征及驅(qū)動因素進行分析,該研究能夠為下一步估算內(nèi)蒙古季節(jié)凍土凍結(jié)深度、多年凍土活動層厚度、模擬凍土面積變化等提供重要參數(shù)數(shù)據(jù),以及為內(nèi)蒙古農(nóng)牧業(yè)生產(chǎn)布局、草場生態(tài)評價及凍害防治等方面提供一定的科學(xué)依據(jù)。

    1 數(shù)據(jù)與方法

    1.1 研究區(qū)概況

    內(nèi)蒙古位于126°04′~97°12′E,37°24′~53°23′N之間(圖1),地處北部邊疆,總面積118.3×104km2,占全國總面積的12.3%。地貌以高原為主,兼有平原、沙地、丘陵與山地,大部分地區(qū)海拔在1000 m以上[22]。氣候?qū)贉貛Т箨懶约撅L(fēng)氣候,處于干旱與半干旱氣候帶,年平均氣溫0~8 ℃,且由東北向西南遞增;全區(qū)降水量少且分布不均,年總降水量50~450 mm,由東北向西南遞減[23]。氣候分異明顯使得研究區(qū)域從東北向西南呈現(xiàn)出寒溫帶針葉林、落葉闊葉林、草甸草原、典型草原、荒漠草原等植被分布[22]。研究區(qū)位于環(huán)北極地區(qū)多年凍土區(qū)歐亞大陸多年凍土南界。冬季嚴(yán)寒,夏季高溫,且局地日夜溫差大,春冬季土壤凍融現(xiàn)象突出,夏秋季土壤干濕交替明顯[24]。多年凍土主要分布于呼倫貝爾地區(qū),其余地區(qū)廣泛分布季節(jié)凍土,淺層表層土壤冬季結(jié)冰,夏季完全解凍。

    圖1 內(nèi)蒙古氣象站點/凍土區(qū)分布及土地覆蓋類型Fig.1 Distribution of meteorological stations/permafrost regions and land cover types in Inner Mongolia

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

    使用的45 個氣象臺站1980—2019 年日平均0 cm 地表溫度觀測數(shù)據(jù)來自于國家氣象信息中心的中國地面氣候資料日值數(shù)據(jù)集(http://data.cma.cn/)。為更好地體現(xiàn)研究區(qū)的凍融特征,盡量采用過去40 a來不同觀測站的共同時間數(shù)據(jù)來求取地表凍融指數(shù)。對缺測數(shù)據(jù)采用如下方式內(nèi)插[25]:缺測1 d,對其前后各1 d數(shù)據(jù)的平均值內(nèi)插;缺測2 d,在缺測當(dāng)天選擇該日前2 d 的數(shù)值取平均,缺測第2 d選取該日后2 d的數(shù)值取平均;如果有連續(xù)2 d以上的缺測,可以利用對相鄰站同一月份的溫度平均值做線性相關(guān)性檢查,再選取溫度相關(guān)性最佳的站點用線性回歸的溫度估計值內(nèi)插。采用反距離插值法(inverse distance weighted)對內(nèi)蒙古45個氣象站點的地表凍融指數(shù)進行空間插值,插值分辨率為34 km。

    凍土類型分區(qū)依據(jù)青藏高原科學(xué)數(shù)據(jù)中心(http://data.tpdc.ac.cn)提供的1:1000凍土分區(qū)和類型地圖資料進行劃分,通過對凍土類型整合,整體劃分為多年凍土區(qū)和季節(jié)凍土區(qū),用來對地表凍融指數(shù)和影響因素分區(qū)討論。0.1 m和0.4 m表層土壤含水量、降水量、雪深數(shù)據(jù)均源自我國國家氣象信息中心的中國第一代全球大氣和陸面再分析產(chǎn)品CRA/Land-逐月陸面產(chǎn)品(http://data.cma.cn/),空間分辨率為34 km,用來分析內(nèi)蒙古地表凍融指數(shù)影響因素時空變化。

    GIMMS NDVI3g 數(shù)據(jù)由美國馬里蘭大學(xué)GLCF(Global Land Cover Facility)研究組生產(chǎn)的15 d 合成產(chǎn) 品 數(shù) 據(jù) (https://ecocast.arc.nasa.gov/data/pub/gimms/),MODIS NDVI 數(shù)據(jù)為美國國家航天航空局(NASA)提供的MODIS 16 d合成產(chǎn)品數(shù)據(jù)(https://mo?dis.gsfc.nasa.gov/data/)。在GIMMS NDVI 和MODIS NDVI 的基礎(chǔ)上,采用線性回歸分析,對2 種數(shù)據(jù)進行時間序列的重建,從而使數(shù)據(jù)重采樣到同一個時空尺度上。使用重合時間的GIMMS NDVI 和MO?DIS NDVI 進行建模,根據(jù)建模結(jié)果對GIMMS 數(shù)據(jù)作延長,獲取1982—2019年NDVI數(shù)據(jù)。對于NDVI數(shù)據(jù)的預(yù)處理,主要分為坐標(biāo)變換、數(shù)值裁剪、數(shù)據(jù)轉(zhuǎn)換、異常值刪除、月最大值獲取等過程,使參與計算的所有NDVI 值在0~1 范圍內(nèi),根據(jù)研究區(qū)土地覆蓋分類圖(圖1),結(jié)合研究區(qū)實際情況,認(rèn)為ND?VI<0.15 的地區(qū)無植被覆蓋。考慮到數(shù)據(jù)空間分辨率的一致性,對NDVI數(shù)據(jù)重采樣為34 km。

    1.3 研究方法

    (1)地表凍融指數(shù)計算?;贙lene 等[26]和Frauenfeld等[15]對凍結(jié)期和融化期的限制,界定凍結(jié)期介于每年7 月1 日至次年6 月30 日,融化期介于每年1月1日至12月31日,以最大可能統(tǒng)計一年中寒冷期(溫度<0 ℃)和溫暖期(溫度>0 ℃)的溫度累加值,并確保負(fù)溫和正溫計算分別在整個寒冷期和溫暖期內(nèi)。年凍結(jié)指數(shù)計算如下:

    式中:FI和TI分別為凍結(jié)指數(shù)和融化指數(shù)(℃·d);Ti和Tj分別為負(fù)值逐日地表溫度和正值逐日地表溫度,單位為℃;Nf和NT為年內(nèi)地表溫度小于和大于0 ℃的天數(shù),單位為d。

    (2)趨勢分析。利用一元線性回歸方法計算凍融指數(shù)和影響因素的年際變化趨勢,并通過最小二乘法擬合的斜率反映,當(dāng)斜率大于0時,說明有增大的傾向,當(dāng)斜率小于0時,說明有下降的傾向。其變化趨勢的計算公式為:

    式中:Slope表示地表凍融指數(shù)或影響因素的變化趨勢;n為時間序列長度;為研究對象X第i年的平均值。

    (3)相關(guān)分析法。通過計算地表凍融指數(shù)與對應(yīng)的影響因素相關(guān)系數(shù),表示其與地表凍融指數(shù)的相關(guān)性。計算公式如下:

    式中:n表示研究時間段的累計年數(shù);rxy表示x和y的相關(guān)系數(shù);xi和yi分別表示第i年x和y的值;和表示2 個要素樣本值的平均值;rxy在-1~1 之間,隨著絕對值的增大,地表凍融指數(shù)與影響因素的變化關(guān)系也相應(yīng)較好;反之越差,當(dāng)rxy>0 表示相關(guān)性為正,rxy<0,表示相關(guān)性為負(fù)[27]。最后采用T檢驗法進行顯著性檢驗,將地表凍融指數(shù)與影響因素的相關(guān)程度劃分為顯著相關(guān)(0.01≤P<0.05)和不顯著相關(guān)(P≥0.05)。

    (4)灰色關(guān)聯(lián)分析。采用灰色關(guān)聯(lián)分析法,研究影響地表凍融指數(shù)的主要因子及重要性?;疑P(guān)聯(lián)度可以測量2 個系統(tǒng)之間的關(guān)系,并確定相關(guān)因子中的主控因子[28]。地表凍融指數(shù)的時間序列作為目標(biāo)變量,將0.1 m土壤含水量、0.4 m土壤含水量、降水量等影響因素作為影響變量?;疑P(guān)聯(lián)系數(shù)表達為:

    2 結(jié)果與分析

    2.1 內(nèi)蒙古地表凍融指數(shù)時空變化

    2.1.1 地表凍融指數(shù)的時間變化特征 通過對內(nèi)蒙古地區(qū)1980—2019 年的45 個觀測站SFI(地表凍結(jié)指數(shù))和STI(地表融化指數(shù))進行分析,得到研究區(qū)40 a 整體和不同凍土區(qū)SFI 與STI 的時間變化特征。如圖2a和圖2b所示,整體上SFI近40 a呈現(xiàn)顯著下降趨勢,多年凍土區(qū)SFI 下降速率高于季節(jié)凍土區(qū),內(nèi)蒙古SFI 多年變化范圍介于956.1~1848.3 ℃·d,最低值和最高值分別發(fā)生于2016年和1984年,年均值為1457.0 ℃·d,變化率為-156.4 ℃·d·(10a)-1(P<0.01),表明土壤的季節(jié)凍結(jié)潛力有連續(xù)降低的趨勢[30]。進一步分析各站點SFI 變化特征,近40 a 所有站點的SFI 均呈下降趨勢,但各站SFI變化特征差異也較大,變化最小值和最大值分別出現(xiàn)于阿拉善盟阿拉善右旗站點[-20.7 ℃·d·(10a)-1]和興安盟阿爾山站點[-711.1 ℃·d·(10a)-1]。如圖2c 和圖2d 所示,整體上STI 近40 a 呈現(xiàn)顯著上升趨勢,多年凍土區(qū)STI上升速率高于季節(jié)凍土區(qū),內(nèi)蒙古STI 多年變化范圍介于3717.6~4442.3 ℃·d,最低值和最高值分別發(fā)生于1984 年和2007 年,年均值為4040.9 ℃·d,上升速率為152.4 ℃·d·(10a)-1(P<0.01),說明土壤融化期間最大融化深度增加[30]。1980—2019年間,所有站點的STI均呈上升趨勢,各站STI 變化特征存在差異,變化特征最小值和最大值分別出現(xiàn)于赤峰市赤峰站點[43.1 ℃·d·(10a)-1]和赤峰市巴林左旗站點[221.1 ℃·d·(10a)-1]。

    圖2 1980—2019年內(nèi)蒙古整體和不同凍土區(qū)地表凍融指數(shù)年際變化趨勢Fig.2 Interannual variation trend of surface freezing-thawing indices in Inner Mongolia and different permafrost regions from 1980 to 2019

    2.1.2 地表凍融指數(shù)空間分布與變化 內(nèi)蒙古SFI多年平均值的空間分布特征如圖3a 所示,受經(jīng)度、緯度和海拔的綜合影響,SFI 均值總體呈現(xiàn)西低東高、由西南向東北逐漸增加的特征,其中SFI受緯度的影響最大(表1),相關(guān)系數(shù)達到0.92(P<0.01),每上升1°,SFI 增加165.2 ℃·d,該分布特征同環(huán)北極凍融指數(shù)分布變化特征一致[12,16]。SFI 高值區(qū)位于呼倫貝爾北部地區(qū),SFI 低值區(qū)位于阿拉善西部地區(qū),與凍土區(qū)分布狀況類似,多年凍土區(qū)SFI 普遍大于季節(jié)凍土區(qū)。內(nèi)蒙古SFI 年際變化率空間分布特征見圖3b,SFI 年際變化率的變化范圍介于-10.3~-710.6 ℃·d·(10a)-1,SFI 變化率在空間上自西南向東北呈增加趨勢,表明內(nèi)蒙古東北地區(qū)SFI 下降速率較西南地區(qū)更快。

    內(nèi)蒙古STI 多年平均值的空間分布特征見圖3c,STI 均值的空間分布特征呈現(xiàn)西部高,東部低,由西南向東北逐漸減小的趨勢,多年凍土區(qū)的STI普遍小于季節(jié)凍土區(qū),在呼和浩特、赤峰等地出現(xiàn)了點狀高值區(qū),這與近幾年農(nóng)業(yè)開發(fā)、城市擴張、工業(yè)發(fā)展等社會經(jīng)濟活動密切相關(guān)。同時,上述地區(qū)處于生態(tài)環(huán)境較為脆弱,人為活動對環(huán)境的影響更加敏感。STI同緯度的相關(guān)系數(shù)最高(表1),為-0.73(P<0.01),每上升1°STI 減少141.1 ℃·d,STI 高值區(qū)位于阿拉善西北部地區(qū),STI 低值區(qū)位于呼倫貝爾中部地區(qū)。內(nèi)蒙古STI年際變化率空間分布特征見圖3d,STI上升速率較高的地區(qū)(Slope>150 ℃·d·(10a)-1)主要分布在研究區(qū)中部、中東部和東北部分地區(qū),可能是因為STI 上升速率較高的地區(qū)受人類活動等因素的影響,其次,呼倫貝爾高原大部分地區(qū)的凍融指數(shù)變化率最大,這也證實了該區(qū)域(原島狀凍土區(qū))目前處于季節(jié)性凍土的狀態(tài)[31]。

    表1 地表凍融指數(shù)與地理因子的相關(guān)系數(shù)Tab.1 Correlation coefficient between surface freezingthawing index and geographical factors

    圖3 1980—2019年內(nèi)蒙古地表凍融指數(shù)均值和變化趨勢Fig.3 Mean value and variation trend of surface freezing-thawing index in Inner Mongolia from 1980 to 2019

    2.2 內(nèi)蒙古地表凍融指數(shù)影響因素時空變化

    2.2.1 內(nèi)蒙古地表凍融指數(shù)影響因素時間變化 如圖4 所示,近40 a 研究區(qū)0.1 m 和0.4 m 土壤含水量年均值分別介于0.1401~0.1673 m3·m-3和0.1388~0.1779 m3·m-3之間,0.4 m 土壤含水量略大于0.1 m土壤含水量,不同深度土壤含水量的總體變化趨勢均顯示出向上態(tài)勢,平均年際變化率為0.001 m3·m-3·(10a)-1,而0.1 m 和0.4 m 土壤含水量并不是呈現(xiàn)一致性的向上態(tài)勢,而是隨氣候波動經(jīng)歷了“上升-下降-上升”3 個波動性變化階段,且0.1 m 和0.4 m 深度的土壤含水量年際變化趨勢大體相同:1980—1998年為波動上升階段,1998—2007年呈波動下降趨勢,2007—2019 年又呈波動上升趨勢;降水量以10.0 mm·(10a)-1的速率呈不顯著增加趨勢(P>0.05),年均值介于240.1~415.9 mm 之間,降水量與0.1 m和0.4 m土壤含水量的年際波動變化趨勢極為相似,同樣經(jīng)歷了“上升-下降-上升”3 個波動變化階段,1980—1998 年、2007—2019 年屬于降水偏豐的年代;1998—2007 年屬于降水偏枯的年代;研究區(qū)雪深以11.4~32.8 mm 的起伏趨勢為主,多年積雪厚度為22.1 mm,年平均雪深以-0.82 mm·(10a)-1的速率呈微弱減少趨勢(P>0.05);基于1982—2019年內(nèi)蒙古NDVI 數(shù)據(jù),發(fā)現(xiàn)內(nèi)蒙古NDVI 總體以0.007·(10a)-1的速率呈顯著增加趨勢(P<0.05),年均值介于0.41~0.47之間。

    圖4 1980—2019年影響因素年際變化趨勢Fig.4 Interannual variation trend of impact factors from 1980 to 2019

    2.2.2 內(nèi)蒙古地表凍融指數(shù)影響因素空間分布與變化 如圖5 所示,研究區(qū)內(nèi)各影響因素空間分布整體上存在明顯的緯度地帶性和區(qū)域性規(guī)律,且均呈現(xiàn)出由西南向東北遞減的發(fā)展趨勢,與SFI 空間分布特征一致,與STI 空間分布特征相反。但上述要素在部分地區(qū)空間分布和空間變化存在一定差異:0.1 m和0.4 m土壤含水量呈現(xiàn)出東北部和中西部土壤水分相對較高、西部相對較低的空間格局,可能受到不同土地覆蓋類型、積雪分布及人為活動等多重因素的共同作用;由于大興安嶺、陰山、賀蘭山等地勢因素的阻礙,東南地區(qū)的水汽抬升,容易造成降雨,西北山區(qū)的水汽含量偏低,因此東南地區(qū)的降水量比西北地區(qū)大;雪深與降水量分布特征基本一致,但與溫度相反,研究區(qū)東北部降水量較多、年均地表溫度較低且森林覆蓋率高,有利于積雪的形成和積累,中部和西部地區(qū)降水較少,平均地表溫度較高,植被覆蓋率較低,不利于積雪的形成與累積,使得雪深較淺且分布均勻;內(nèi)蒙古NDVI由西南向東北呈增加趨勢,部分地區(qū)NDVI 存在區(qū)域性特征和空間異質(zhì)性較明顯,可能與灌溉和開墾農(nóng)田等人類活動有關(guān)。

    圖5 1980—2019年地表凍融指數(shù)影響因素多年平均Fig.5 Annual average of surface freezing-thawing index influencing factors from 1980 to 2019

    如圖6所示,近40 a來內(nèi)蒙古大部分地區(qū)0.1 m和0.4 m 土壤含水量年際變化呈上升態(tài)勢,0.1 m 和0.4 m 土壤含水量呈減少態(tài)勢的區(qū)域分別僅占研究區(qū)整體的19.9%和36.1%,且大部分集中在研究區(qū)西北部和東北部的多年凍土區(qū);內(nèi)蒙古中部和西部的年降水量有所增加,特別是中偏南地區(qū)降水量增加較多,東部降水量則有所減少,特別是在呼倫貝爾西北部、赤峰中部以及通遼中部地區(qū)降水量減少情況較為突出;研究區(qū)38.5%的區(qū)域雪深呈下降態(tài)勢,且10.8%的區(qū)域雪深平均下降速度在0.5 mm·(10a)-1以上,主要分布在呼倫貝爾東北部,并位于多年凍土區(qū);61.5%的區(qū)域雪深呈現(xiàn)增加趨勢,主要分布于西南部的沙漠及中部和東部的部分草原,多位于季節(jié)凍土區(qū),且以上區(qū)域的雪深平均增加速率僅為0.4 mm·(10a)-1,這一變化特征可能與研究區(qū)近年來溫度和降水變化有關(guān);研究時段內(nèi)蒙古有75.2%區(qū)域的NDVI 呈現(xiàn)增長的趨勢,僅24.8%的區(qū)域呈現(xiàn)下降的趨勢,下降區(qū)域以錫林郭勒盟大部、赤峰市北部、通遼市北部、興安盟南部、呼倫貝爾市林區(qū)為主,這一變化特征可能與上述地區(qū)近年來溫度升高、降水減少及人類活動有關(guān)。

    圖6 1980—2019年地表凍融指數(shù)影響因素年際變化Fig.6 Interannual variation of surface freezing-thawing index influencing factors from 1980 to 2019

    2.3 地表凍融指數(shù)與影響因素相關(guān)分析

    地表凍融指數(shù)的變化受多種因素的影響,具有較強的隨機性,適于進行灰色關(guān)聯(lián)度分析。如表2所示,內(nèi)蒙古地表凍融指數(shù)年均值與0.1 m 土壤含水量、0.4 m土壤含水量、降水量、雪深、NDVI的灰色關(guān)聯(lián)度系數(shù)皆大于0.65,說明地表凍融指數(shù)變化與各影響因素存在較大的相關(guān)性。SFI與影響因素的灰色關(guān)聯(lián)序為:0.4 m 土壤含水量>0.1 m 土壤含水量>NDVI>降水量>雪深,說明了0.4 m 土壤含水量對SFI變化的影響較大,是影響內(nèi)蒙古整體SFI變化的主控因素,0.1 m 土壤含水量、NDVI、降水量次之,最后是雪深;STI 與影響因素的灰色關(guān)聯(lián)序為:NDVI>0.1 m 土壤含水量>0.4 m 土壤含水量>降水量>雪深,表明NDVI 對STI 變化的影響較大,是影響內(nèi)蒙古整體STI變化的主控因素,0.1 m土壤含水量、0.4 m土壤含水量、降水量次之,最后是雪深。

    表2 地表凍融指數(shù)與影響因素灰色關(guān)聯(lián)系分析Tab.2 Grey correlation analysis between surface freezing-thawing index and influencing factors

    為探討地表凍融指數(shù)對影響因素變化的動態(tài)響應(yīng),基于像元尺度逐像元計算了內(nèi)蒙古地表凍融指數(shù)與影響因素的相關(guān)關(guān)系。SFI 與各影響因素(除雪深外)相關(guān)關(guān)系有超過65%的地區(qū)呈負(fù)相關(guān)(表3),如圖7所示,從空間分布來看,在多年凍土區(qū)大部分地區(qū)呈現(xiàn)正相關(guān)關(guān)系,但在季節(jié)凍土區(qū)大部分地區(qū)卻呈現(xiàn)負(fù)相關(guān)關(guān)系,造成不同凍土區(qū)SFI 與影響因素相關(guān)性差異的原因可能與影響因素的時空變化差異有關(guān)。其中,0.1 m土壤含水量與SFI顯著負(fù)相關(guān)(P<0.05)的區(qū)域大部分在季節(jié)凍土區(qū)的中部和東南部,顯著正相關(guān)(P<0.05)區(qū)域集中分布在多年凍土區(qū)北部小部分地區(qū),0.4 m 土壤含水量與SFI 顯著負(fù)相關(guān)(P<0.05)的地區(qū)在季節(jié)凍土區(qū)的南部和東南部均有分布,顯著正相關(guān)(P<0.05)的地區(qū)則在多年凍土區(qū)北部呈現(xiàn)片狀分布;降水量與SFI顯著負(fù)相關(guān)(P<0.05)的地區(qū)主要在鄂爾多斯市以東、呼和浩特市和錫林郭勒盟以南,呈正相關(guān)關(guān)系的區(qū)域主要在多年凍土區(qū)北部和阿拉善中西部;NDVI 與SFI顯著負(fù)相關(guān)(P<0.05)的區(qū)域主要在鄂爾多斯東部、呼和浩特市南部和多年凍土區(qū)部分林地,顯著正相關(guān)(P<0.05)區(qū)域主要在多年凍土區(qū)部分林地和濕地地區(qū);雪深與SFI顯著負(fù)相關(guān)(P<0.05)的區(qū)域主要在阿拉善西南部和季節(jié)凍土區(qū)中部地區(qū)零星分布,顯著正相關(guān)(P<0.05)區(qū)域集中在多年凍土區(qū)北部和阿拉善西北部小片地區(qū)。

    圖7 SFI與影響因素相關(guān)性T檢驗Fig.7 Correlation T test between SFI and influencing factors

    表3 SFI與影響因素相關(guān)系數(shù)面積占比Tab.3 Area ratio of correlation coefficient between SFI and influencing factors

    STI 與影響因素的相關(guān)關(guān)系以負(fù)相關(guān)為主(表4)。如圖8所示,雖然STI與各影響因素在多年凍土區(qū)及大部分地區(qū)呈現(xiàn)負(fù)相關(guān)關(guān)系,但在季節(jié)凍土區(qū)部分地區(qū)卻呈現(xiàn)正相關(guān)關(guān)系,造成不同凍土區(qū)STI與影響因素相關(guān)性差異的原因可能與影響因素的時空變化差異有關(guān)。其中,0.1 m 土壤含水量與STI顯著負(fù)相關(guān)(P<0.05)區(qū)域集中在多年凍土區(qū)東部和多年凍土區(qū)周圍的盟市,顯著正相關(guān)(P<0.05)的區(qū)域在黃河流域內(nèi)蒙古地區(qū)有所分布;0.4 m土壤含水量與STI 顯著負(fù)相關(guān)(P<0.05)區(qū)域主要在多年凍土區(qū)北部和阿拉善西部零星地區(qū),顯著正相關(guān)(P<0.05)的區(qū)域主要在阿拉善東部、巴彥淖爾市大部和鄂爾多斯市周邊;降水量與STI顯著負(fù)相關(guān)(P<0.05)區(qū)域主要在錫林郭勒盟東部、興安盟東部及多年凍土區(qū)東北部,呈正相關(guān)關(guān)系的區(qū)域主要在季節(jié)凍土區(qū)中部和西部;雪深與STI顯著負(fù)相關(guān)(P<0.05)區(qū)域主要在錫林郭勒盟西北部、赤峰市北部及多年凍土區(qū)東北部,顯著正相關(guān)(P<0.05)的區(qū)域集中在巴彥淖爾市中部。NDVI 與SFI 顯著負(fù)相關(guān)(P<0.05)的區(qū)域主要在赤峰市北部、興安盟、多年凍土區(qū)部分林地和錫林郭勒盟荒漠草原地區(qū),顯著正相關(guān)(P<0.05)的區(qū)域主要在黃河流域內(nèi)蒙古地區(qū)、通遼市和多年凍土區(qū)北部。

    圖8 STI與影響因素相關(guān)性T檢驗Fig.8 Correlation T test between STI and influencing factors

    表4 STI與影響因素相關(guān)系數(shù)面積占比Tab.4 Area ratio of correlation coefficient between STI and influencing factors

    3 討論

    氣候變暖背景下,內(nèi)蒙古1980—2019年間地表凍結(jié)和融化指數(shù)分別顯著上升和下降,多年凍土區(qū)地表凍融指數(shù)變化速率較季節(jié)凍土區(qū)更快。研究結(jié)果基本與Peng等[12]、Wu等[18]、李昀赟等[32]分別在北半球、青藏高原、東北地區(qū)的結(jié)論相吻合,說明在內(nèi)蒙古整體凍土環(huán)境退化明顯,多年凍土區(qū)冬季溫度變暖趨勢更為顯著,其凍土退化更嚴(yán)重,該區(qū)地表凍融指數(shù)的顯著變化可能導(dǎo)致多年凍土區(qū)面積縮減,凍結(jié)深度減少,活動層厚度增加等一系列凍土退化現(xiàn)象的發(fā)生[12];季節(jié)凍土區(qū)的地表凍融指數(shù)變化速率同樣不容忽視,該區(qū)是內(nèi)蒙古農(nóng)牧業(yè)生產(chǎn)主要布局地區(qū),凍融指數(shù)的變化同農(nóng)作物的種植,畜牧業(yè)的發(fā)展息息相關(guān)。氣候變暖雖然是造成該區(qū)多年凍土區(qū)凍土退化顯著的根本性因素,但是多種人為經(jīng)濟活動對研究區(qū)多年凍土退化起到了加速促進作用[24]。近40 a研究區(qū)表層土壤含水量、降水量的增加趨勢表明,內(nèi)蒙古地區(qū)的凍土環(huán)境正向暖濕化發(fā)展,與馬梓策等[33]得出的研究結(jié)論相一致。但是基于本文研究發(fā)現(xiàn),不同凍土區(qū)的發(fā)展情況存在較大差異,多年凍土區(qū)正向暖干化發(fā)展,與研究區(qū)整體和季節(jié)凍土區(qū)大部分地區(qū)變化趨勢相反。地面溫度升高、降水減少、潛在蒸發(fā)量增加是導(dǎo)致土壤含水量下降的主要因素[34]。地表凍融指數(shù)是根據(jù)地表溫度計算所得,與氣溫和地表溫度有密切的聯(lián)系,但表層土壤含水量、降水量、雪深和NDVI 同樣對凍融指數(shù)的變化起著至關(guān)重要的作用,總體上,表層土壤含水量和NDVI對地表凍融指數(shù)的變化起主導(dǎo)作用。此外,凍土環(huán)境退化對生物多樣性的作用是負(fù)面的。溫度增高、雨量減少的情況對草原植物生長發(fā)育不利,對森林、濕地植物有促進作用[35];季節(jié)凍土地帶冬季增溫明顯,農(nóng)作物凍害損失減少,從而促進了糧食作物單產(chǎn)的增加和質(zhì)量提高,利于家畜健康過冬度春,同時也會使病蟲害及蟲卵易于過冬,容易造成病蟲害的防治困難[36]。

    20世紀(jì)80年代前研究區(qū)氣象數(shù)據(jù)缺測較多,本研究只采用了1980—2019年的氣象資料,下一步還需依據(jù)更長時間序列和更多的站點數(shù)據(jù)開展研究。在對地表凍融指數(shù)處理后,需要根據(jù)凍融指數(shù),進一步確定研究區(qū)多年凍土空間分布范圍和活動層厚度及季節(jié)凍土凍結(jié)深度的時空變化特征,對內(nèi)蒙古地區(qū)凍土對氣候變化的響應(yīng)作出較為全面、精確的預(yù)測。目前,由于中國全球陸面再分析產(chǎn)品(CRA/Land)數(shù)據(jù)在凍土研究中仍具有空間分辨率略小的缺陷,在今后凍土變化研究中,擬將與空間分辨率更高的再分析資料及實際觀測資料相結(jié)合,以盡可能降低研究結(jié)論的不確定性。

    4 結(jié)論

    本文基于內(nèi)蒙古氣象站1980—2019 年日均地表溫度數(shù)據(jù)、結(jié)合中國第一代全球陸面再分析產(chǎn)品(CRA)數(shù)據(jù)以及NDVI 數(shù)據(jù),借助GIS 空間分析功能,通過RStudio編程處理分析了內(nèi)蒙古近40 a地表凍融指數(shù)時空變化格局及關(guān)鍵影響要素,主要結(jié)論如下:

    (1)1980—2019年間,內(nèi)蒙古SFI與STI分別呈現(xiàn)明顯減少與增加的態(tài)勢,緯度是制約內(nèi)蒙古地表凍融指數(shù)空間分布的重要因素,與SFI 呈現(xiàn)極明顯正相關(guān)的關(guān)系,則與STI 相反。多年凍土區(qū)地表凍融指數(shù)變化速率較季節(jié)凍土區(qū)更為劇烈,對氣候變暖的響應(yīng)更為敏感。

    (2)近40 a 研究區(qū)0.1 m 和0.4 m 土壤含水量、降水量、NDVI 均呈增加趨勢,僅雪深呈減少趨勢。其中不同深度土壤含水量平均年際變化率為0.001 m3·m-3·(10a)-1;降水量以10.0 mm·(10a)-1的速率呈不顯著增加趨勢;NDVI總體以0.007·(10a)-1的速率呈顯著增加趨勢;年平均雪深以-0.82 mm·(10a)-1的速率呈微弱減少趨勢。

    (3)多年凍土區(qū)0.1 m 和0.4 m 土壤含水量、降水量、NDVI 及雪深的多年均值普遍大于季節(jié)凍土區(qū),同樣影響因素在年際變化上也表現(xiàn)出不同的空間差異性,多年凍土區(qū)呈暖干化發(fā)展趨勢,季節(jié)凍土區(qū)呈暖濕化發(fā)展趨勢。

    (4)地表凍融指數(shù)與影響因素以負(fù)相關(guān)關(guān)系為主,SFI 與影響因素在多年凍土區(qū)大部呈正相關(guān)關(guān)系,在季節(jié)凍土區(qū)大部呈負(fù)相關(guān)關(guān)系,STI反之。內(nèi)蒙古地表凍融指數(shù)變化受影響因素共同驅(qū)動,而0.4 m土壤含水量是影響SFI變化的主導(dǎo)因素,NDVI是影響STI變化的主導(dǎo)因素。

    猜你喜歡
    多年凍土凍土凍融
    中國東北多年凍土退化對植被季節(jié)NDVI 的影響研究
    北極凍土在求救
    凍土下的猛犸墳場
    太陽能制冷在多年凍土熱穩(wěn)定維護中的傳熱效果研究
    間苯三酚在凍融胚胎移植中的應(yīng)用
    反復(fù)凍融作用下巖橋破壞的試驗研究
    多年凍土地基隔熱保溫技術(shù)研究綜述
    多年凍土區(qū)鐵路路堤臨界高度研究
    26
    降調(diào)節(jié)方案在凍融胚胎移植周期中的應(yīng)用
    成人av一区二区三区在线看| 又紧又爽又黄一区二区| 中文字幕精品免费在线观看视频| 99国产精品一区二区三区| 国产精品免费视频内射| 国产精品亚洲一级av第二区| 12—13女人毛片做爰片一| 极品少妇高潮喷水抽搐| avwww免费| 91麻豆精品激情在线观看国产 | 久久久国产一区二区| 两个人看的免费小视频| 一级毛片精品| 中文字幕精品免费在线观看视频| 女同久久另类99精品国产91| 91成人精品电影| 久久久国产欧美日韩av| av在线播放免费不卡| 国产一区二区三区综合在线观看| 在线观看免费午夜福利视频| 热re99久久国产66热| 老司机靠b影院| 亚洲av美国av| 亚洲五月婷婷丁香| 国产日韩一区二区三区精品不卡| 黑人操中国人逼视频| av免费在线观看网站| 日韩视频在线欧美| cao死你这个sao货| 国产男女内射视频| 9热在线视频观看99| 亚洲少妇的诱惑av| 国产又色又爽无遮挡免费看| 他把我摸到了高潮在线观看 | 精品一区二区三区av网在线观看 | 亚洲专区中文字幕在线| 国产男靠女视频免费网站| av欧美777| avwww免费| 成年人免费黄色播放视频| 香蕉丝袜av| 99riav亚洲国产免费| 久久久久视频综合| 我的亚洲天堂| 中文字幕精品免费在线观看视频| 最黄视频免费看| 久久久精品国产亚洲av高清涩受| 美女扒开内裤让男人捅视频| 色视频在线一区二区三区| 一本综合久久免费| 9热在线视频观看99| 黄色 视频免费看| 色老头精品视频在线观看| 女人爽到高潮嗷嗷叫在线视频| 99精品在免费线老司机午夜| 亚洲av片天天在线观看| 国产精品免费视频内射| 日韩制服丝袜自拍偷拍| 午夜福利一区二区在线看| 精品乱码久久久久久99久播| 国产精品.久久久| 亚洲一码二码三码区别大吗| 后天国语完整版免费观看| av又黄又爽大尺度在线免费看| 人人澡人人妻人| 亚洲av片天天在线观看| 水蜜桃什么品种好| 在线亚洲精品国产二区图片欧美| 欧美激情极品国产一区二区三区| 两人在一起打扑克的视频| 丰满少妇做爰视频| 天天影视国产精品| 国内毛片毛片毛片毛片毛片| 天天操日日干夜夜撸| av网站在线播放免费| 成年人午夜在线观看视频| av在线播放免费不卡| 夜夜爽天天搞| 高清在线国产一区| 悠悠久久av| 国产午夜精品久久久久久| 欧美日韩福利视频一区二区| 中文字幕最新亚洲高清| 热99国产精品久久久久久7| 好男人电影高清在线观看| 亚洲,欧美精品.| 黄片大片在线免费观看| 欧美精品av麻豆av| 色老头精品视频在线观看| 亚洲少妇的诱惑av| 国产精品.久久久| a在线观看视频网站| 午夜激情久久久久久久| 国产精品99久久99久久久不卡| 久久精品aⅴ一区二区三区四区| 亚洲国产欧美一区二区综合| 久久午夜综合久久蜜桃| 少妇粗大呻吟视频| 久久中文字幕一级| 法律面前人人平等表现在哪些方面| 亚洲精品一卡2卡三卡4卡5卡| 一进一出好大好爽视频| 最新的欧美精品一区二区| 国产麻豆69| 香蕉国产在线看| 免费看十八禁软件| 青草久久国产| 亚洲熟女精品中文字幕| av网站在线播放免费| 99国产精品99久久久久| 女同久久另类99精品国产91| 亚洲精品久久成人aⅴ小说| 免费高清在线观看日韩| a在线观看视频网站| 亚洲第一欧美日韩一区二区三区 | 999久久久精品免费观看国产| 在线观看一区二区三区激情| 国产一区二区激情短视频| 日韩人妻精品一区2区三区| 成人国语在线视频| 亚洲精品av麻豆狂野| 日韩一区二区三区影片| 新久久久久国产一级毛片| 国产欧美日韩精品亚洲av| 一级片免费观看大全| 精品午夜福利视频在线观看一区 | 丝袜喷水一区| 十八禁网站网址无遮挡| 男女免费视频国产| 国产色视频综合| 国产一区二区在线观看av| 夫妻午夜视频| 色播在线永久视频| 动漫黄色视频在线观看| 国产又色又爽无遮挡免费看| 国产精品av久久久久免费| 69精品国产乱码久久久| 人人澡人人妻人| 黄色成人免费大全| 亚洲欧美色中文字幕在线| 国产男女超爽视频在线观看| 国产熟女午夜一区二区三区| 伦理电影免费视频| 手机成人av网站| videos熟女内射| 国产午夜精品久久久久久| 久久精品亚洲av国产电影网| 丰满人妻熟妇乱又伦精品不卡| 成年人午夜在线观看视频| 国产精品麻豆人妻色哟哟久久| 中文字幕另类日韩欧美亚洲嫩草| 黄片播放在线免费| 欧美在线黄色| 亚洲人成电影免费在线| 欧美国产精品va在线观看不卡| 考比视频在线观看| 国产成人一区二区三区免费视频网站| 一级a爱视频在线免费观看| 中文字幕人妻丝袜一区二区| 欧美一级毛片孕妇| 纯流量卡能插随身wifi吗| 免费观看av网站的网址| av天堂久久9| 欧美+亚洲+日韩+国产| 51午夜福利影视在线观看| av一本久久久久| 欧美日韩亚洲高清精品| 精品福利永久在线观看| 99精品在免费线老司机午夜| 国产又色又爽无遮挡免费看| 女警被强在线播放| 又紧又爽又黄一区二区| 欧美激情高清一区二区三区| 午夜免费鲁丝| 天天躁夜夜躁狠狠躁躁| 国产在线观看jvid| 老汉色av国产亚洲站长工具| 最近最新免费中文字幕在线| 亚洲人成电影免费在线| 91麻豆av在线| 黄色毛片三级朝国网站| 婷婷成人精品国产| 久久av网站| 亚洲人成电影观看| 正在播放国产对白刺激| 中亚洲国语对白在线视频| 高清av免费在线| 在线十欧美十亚洲十日本专区| 亚洲美女黄片视频| 国产一区二区三区在线臀色熟女 | 欧美变态另类bdsm刘玥| 成人手机av| 欧美一级毛片孕妇| 热99国产精品久久久久久7| 成人特级黄色片久久久久久久 | 亚洲男人天堂网一区| 久久久久久久精品吃奶| 一级片'在线观看视频| 欧美日韩精品网址| 亚洲精品一二三| 国产色视频综合| 久久99热这里只频精品6学生| 久久久精品区二区三区| 精品少妇黑人巨大在线播放| 欧美精品一区二区免费开放| 男男h啪啪无遮挡| 乱人伦中国视频| 免费高清在线观看日韩| 一本—道久久a久久精品蜜桃钙片| 高清欧美精品videossex| 波多野结衣av一区二区av| 男女免费视频国产| 欧美成人免费av一区二区三区 | 成人国语在线视频| 999精品在线视频| 日韩欧美三级三区| 美女主播在线视频| 一区二区三区国产精品乱码| 亚洲精品在线观看二区| 国产免费福利视频在线观看| 老司机午夜福利在线观看视频 | 在线观看免费视频日本深夜| 久久久欧美国产精品| 亚洲国产精品一区二区三区在线| 十分钟在线观看高清视频www| 亚洲av片天天在线观看| 纯流量卡能插随身wifi吗| 精品少妇黑人巨大在线播放| 麻豆国产av国片精品| 啦啦啦中文免费视频观看日本| 亚洲国产成人一精品久久久| 岛国毛片在线播放| 一进一出抽搐动态| 亚洲国产欧美网| 丝袜喷水一区| 国产成人精品久久二区二区91| 免费久久久久久久精品成人欧美视频| av免费在线观看网站| 免费在线观看完整版高清| 亚洲色图 男人天堂 中文字幕| 精品少妇久久久久久888优播| 嫩草影视91久久| av有码第一页| 女警被强在线播放| 99九九在线精品视频| 午夜免费成人在线视频| 日韩欧美免费精品| 精品国产一区二区久久| 窝窝影院91人妻| 好男人电影高清在线观看| 亚洲七黄色美女视频| www.精华液| 岛国在线观看网站| 久久亚洲精品不卡| 免费看十八禁软件| 一级片'在线观看视频| 成人永久免费在线观看视频 | 一级,二级,三级黄色视频| 狂野欧美激情性xxxx| 在线十欧美十亚洲十日本专区| av欧美777| 十八禁人妻一区二区| 亚洲精品久久午夜乱码| 大型黄色视频在线免费观看| 成在线人永久免费视频| 欧美一级毛片孕妇| 麻豆成人av在线观看| 99在线人妻在线中文字幕 | 亚洲精品国产一区二区精华液| 女同久久另类99精品国产91| 亚洲精品美女久久av网站| 国产麻豆69| 岛国在线观看网站| 蜜桃在线观看..| 别揉我奶头~嗯~啊~动态视频| 美女福利国产在线| 丝瓜视频免费看黄片| 国产有黄有色有爽视频| 亚洲欧洲日产国产| 亚洲精品国产区一区二| 亚洲第一av免费看| 韩国精品一区二区三区| 亚洲av第一区精品v没综合| 少妇的丰满在线观看| 精品国产国语对白av| 狠狠精品人妻久久久久久综合| 另类亚洲欧美激情| 18在线观看网站| 亚洲色图综合在线观看| 嫩草影视91久久| 亚洲午夜理论影院| 国产精品亚洲一级av第二区| 宅男免费午夜| 久久精品亚洲av国产电影网| 久久国产亚洲av麻豆专区| 美女午夜性视频免费| 精品免费久久久久久久清纯 | 香蕉丝袜av| 男人操女人黄网站| 青青草视频在线视频观看| 手机成人av网站| 一进一出好大好爽视频| 欧美亚洲 丝袜 人妻 在线| 成人18禁高潮啪啪吃奶动态图| 在线亚洲精品国产二区图片欧美| 国产有黄有色有爽视频| 亚洲色图av天堂| 久久久久久久国产电影| 久久中文字幕人妻熟女| 多毛熟女@视频| 999精品在线视频| 国产亚洲欧美在线一区二区| 熟女少妇亚洲综合色aaa.| 丰满饥渴人妻一区二区三| 人人妻人人澡人人爽人人夜夜| 精品午夜福利视频在线观看一区 | 伊人久久大香线蕉亚洲五| 久久人妻熟女aⅴ| 一级毛片电影观看| 老熟妇仑乱视频hdxx| 最近最新中文字幕大全免费视频| 91老司机精品| videos熟女内射| 久久久久国内视频| 精品卡一卡二卡四卡免费| 久久精品aⅴ一区二区三区四区| 国产av又大| 一二三四在线观看免费中文在| 亚洲国产精品一区二区三区在线| 亚洲va日本ⅴa欧美va伊人久久| 成在线人永久免费视频| 女人爽到高潮嗷嗷叫在线视频| 视频区欧美日本亚洲| 一级毛片女人18水好多| 十八禁网站免费在线| 9色porny在线观看| 天堂动漫精品| 窝窝影院91人妻| 99国产精品免费福利视频| 老司机在亚洲福利影院| 亚洲欧美一区二区三区黑人| 在线 av 中文字幕| 欧美在线一区亚洲| 亚洲色图av天堂| 亚洲av日韩在线播放| 国产精品久久久久成人av| 欧美日韩亚洲综合一区二区三区_| 91精品国产国语对白视频| 一二三四社区在线视频社区8| 国产精品一区二区在线不卡| 国产精品免费一区二区三区在线 | 中文字幕高清在线视频| 免费一级毛片在线播放高清视频 | 热99国产精品久久久久久7| 91av网站免费观看| 亚洲专区字幕在线| 久久人妻福利社区极品人妻图片| 亚洲专区字幕在线| 在线看a的网站| 欧美国产精品va在线观看不卡| 国产亚洲精品久久久久5区| 欧美乱码精品一区二区三区| 男女之事视频高清在线观看| 女性被躁到高潮视频| 在线看a的网站| 女性被躁到高潮视频| 最近最新中文字幕大全免费视频| 国产精品99久久99久久久不卡| 老熟女久久久| 日本wwww免费看| 男人舔女人的私密视频| 99久久精品国产亚洲精品| 日本五十路高清| 操美女的视频在线观看| 在线观看人妻少妇| 日韩中文字幕视频在线看片| 国产又色又爽无遮挡免费看| www.熟女人妻精品国产| 深夜精品福利| 中文字幕av电影在线播放| 亚洲精品国产精品久久久不卡| 国产99久久九九免费精品| 亚洲熟妇熟女久久| 亚洲av片天天在线观看| 中文字幕人妻丝袜制服| www日本在线高清视频| cao死你这个sao货| 精品少妇一区二区三区视频日本电影| 欧美乱妇无乱码| 美女视频免费永久观看网站| 一区二区av电影网| 涩涩av久久男人的天堂| 婷婷丁香在线五月| 考比视频在线观看| 天堂8中文在线网| 国产一区二区三区综合在线观看| 亚洲成人手机| 男女之事视频高清在线观看| 欧美精品亚洲一区二区| 大香蕉久久成人网| 亚洲精品中文字幕在线视频| 久久久精品国产亚洲av高清涩受| 亚洲国产毛片av蜜桃av| 国产男女内射视频| 久久九九热精品免费| 免费在线观看影片大全网站| 日本wwww免费看| 精品亚洲成国产av| av网站在线播放免费| 曰老女人黄片| 欧美日本中文国产一区发布| 日日爽夜夜爽网站| 不卡av一区二区三区| 一个人免费在线观看的高清视频| 91老司机精品| 精品人妻1区二区| 女警被强在线播放| 老熟妇仑乱视频hdxx| 18在线观看网站| 亚洲一区中文字幕在线| 久久中文字幕一级| 啦啦啦视频在线资源免费观看| a级片在线免费高清观看视频| 精品国产乱码久久久久久小说| 日本黄色视频三级网站网址 | 电影成人av| 亚洲自偷自拍图片 自拍| 久久亚洲精品不卡| 亚洲精品乱久久久久久| 精品国产超薄肉色丝袜足j| 欧美老熟妇乱子伦牲交| 亚洲国产av新网站| 999久久久精品免费观看国产| 国产成人精品无人区| 丰满少妇做爰视频| 99re6热这里在线精品视频| 成人影院久久| tube8黄色片| 性少妇av在线| 老汉色∧v一级毛片| 久久天躁狠狠躁夜夜2o2o| 亚洲专区字幕在线| 91av网站免费观看| 无限看片的www在线观看| 露出奶头的视频| 国产精品电影一区二区三区 | 777久久人妻少妇嫩草av网站| 丝袜美腿诱惑在线| 精品高清国产在线一区| 高清视频免费观看一区二区| 黄色丝袜av网址大全| 久久亚洲精品不卡| 成人国产av品久久久| 搡老岳熟女国产| 欧美中文综合在线视频| 亚洲成av片中文字幕在线观看| 成人精品一区二区免费| 国产免费福利视频在线观看| 在线观看免费视频日本深夜| 女警被强在线播放| 巨乳人妻的诱惑在线观看| 国产高清国产精品国产三级| 久久久久国内视频| 不卡一级毛片| 操美女的视频在线观看| 欧美黑人精品巨大| 精品国内亚洲2022精品成人 | 精品国产乱码久久久久久男人| 欧美日韩av久久| 嫁个100分男人电影在线观看| 亚洲国产成人一精品久久久| 日韩中文字幕视频在线看片| 久久中文看片网| 日日摸夜夜添夜夜添小说| 国产亚洲精品久久久久5区| 一级a爱视频在线免费观看| 欧美精品一区二区免费开放| 新久久久久国产一级毛片| 精品午夜福利视频在线观看一区 | 精品久久久久久电影网| 亚洲精品久久成人aⅴ小说| 夫妻午夜视频| 两个人免费观看高清视频| 久久久国产一区二区| av在线播放免费不卡| 午夜激情久久久久久久| 97人妻天天添夜夜摸| 满18在线观看网站| 99九九在线精品视频| 十八禁高潮呻吟视频| 久久久国产成人免费| 9热在线视频观看99| 亚洲av日韩在线播放| 国产成人精品久久二区二区免费| 国产三级黄色录像| 黄色怎么调成土黄色| 国产一卡二卡三卡精品| 免费av中文字幕在线| 成年动漫av网址| 在线av久久热| 午夜福利视频在线观看免费| 久久午夜综合久久蜜桃| 成在线人永久免费视频| 亚洲精品成人av观看孕妇| 99久久国产精品久久久| 日韩欧美一区视频在线观看| 久久 成人 亚洲| 狠狠狠狠99中文字幕| 国产亚洲一区二区精品| 99热网站在线观看| 1024香蕉在线观看| 欧美日韩国产mv在线观看视频| 电影成人av| 纯流量卡能插随身wifi吗| 亚洲综合色网址| 亚洲av美国av| 亚洲专区字幕在线| 国产色视频综合| 久久久国产欧美日韩av| 国产人伦9x9x在线观看| 国产免费现黄频在线看| 手机成人av网站| av福利片在线| av在线播放免费不卡| av国产精品久久久久影院| 男女之事视频高清在线观看| 国产精品久久久久久精品古装| 国产成人欧美在线观看 | 国产成人影院久久av| 亚洲专区中文字幕在线| 夫妻午夜视频| 亚洲av电影在线进入| 九色亚洲精品在线播放| 午夜老司机福利片| 考比视频在线观看| 精品卡一卡二卡四卡免费| 亚洲国产欧美一区二区综合| 91字幕亚洲| 成人18禁在线播放| 久久国产精品大桥未久av| 少妇 在线观看| 老司机午夜十八禁免费视频| 国产精品久久久人人做人人爽| 99国产精品99久久久久| 麻豆av在线久日| 日本五十路高清| svipshipincom国产片| a级片在线免费高清观看视频| 精品一区二区三区av网在线观看 | 2018国产大陆天天弄谢| 欧美变态另类bdsm刘玥| 精品一区二区三区av网在线观看 | 99精品在免费线老司机午夜| 欧美精品一区二区大全| 午夜精品国产一区二区电影| 午夜久久久在线观看| 久久久国产精品麻豆| 黑人巨大精品欧美一区二区mp4| 国产精品秋霞免费鲁丝片| 亚洲视频免费观看视频| 久久精品亚洲精品国产色婷小说| a在线观看视频网站| 国产男靠女视频免费网站| 黑人欧美特级aaaaaa片| 在线看a的网站| 大型黄色视频在线免费观看| 可以免费在线观看a视频的电影网站| 免费在线观看完整版高清| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美日韩精品网址| 视频区图区小说| 欧美午夜高清在线| 成人永久免费在线观看视频 | 一个人免费看片子| 免费在线观看日本一区| 一级片免费观看大全| 精品亚洲乱码少妇综合久久| 久久久久久久久免费视频了| 亚洲一码二码三码区别大吗| 欧美+亚洲+日韩+国产| 中文字幕最新亚洲高清| 免费女性裸体啪啪无遮挡网站| 悠悠久久av| 欧美国产精品一级二级三级| 大陆偷拍与自拍| 亚洲精品国产区一区二| 国产精品免费大片| 啦啦啦 在线观看视频| 国产精品免费视频内射| 十分钟在线观看高清视频www| 十八禁网站免费在线| 97人妻天天添夜夜摸| e午夜精品久久久久久久| 国产精品二区激情视频| 黄色视频在线播放观看不卡| 国产成人av教育| 国产91精品成人一区二区三区 | 十八禁人妻一区二区| 国产日韩欧美在线精品| 在线观看免费高清a一片| 一级,二级,三级黄色视频| 免费黄频网站在线观看国产| 中文字幕av电影在线播放| 又大又爽又粗| 巨乳人妻的诱惑在线观看| 热re99久久国产66热| 国产色视频综合| 久久亚洲精品不卡| 美国免费a级毛片| 国产熟女午夜一区二区三区| 国产免费av片在线观看野外av| 国产成+人综合+亚洲专区| 精品国产乱码久久久久久小说| 宅男免费午夜| 亚洲午夜精品一区,二区,三区|