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

    莫莫格濕地冷島效應(yīng)遙感分析

    2019-12-30 09:44項小云李鳳秀杜嘉劉雁宋開山趙博宇周昊昊馬學(xué)垚
    安徽農(nóng)業(yè)科學(xué) 2019年23期
    關(guān)鍵詞:冷島反演土地利用

    項小云 李鳳秀 杜嘉 劉雁 宋開山 趙博宇 周昊昊 馬學(xué)垚

    摘要?以2015年春季、夏季、秋季共3期莫莫格濕地的Landsat 8 遙感影像為基礎(chǔ),利用大氣校正法反演研究區(qū)的地表溫度;在地表溫度反演結(jié)果和土地利用數(shù)據(jù)的基礎(chǔ)上,利用ArcGIS軟件的空間分析工具進(jìn)行緩沖區(qū)分析、疊加分析和地統(tǒng)計分析,從而完成莫莫格濕地的冷島效應(yīng)分析。結(jié)果表明,濕地對周邊不同土地利用類型的冷島效應(yīng)(降溫作用)隨著離濕地距離的增加而逐漸減弱。濕地對不同土地利用類型的降溫作用并不相同,對居住地類型的影響范圍在240~540 m,對旱地類型的影響范圍在480~660 m;濕地對鹽堿地的降溫作用比較顯著,影響范圍在600~660 m。濕地在不同季節(jié)對其他土地利用類型降溫作用并不相同,在春夏秋3個季節(jié)中降溫作用由強到弱依次為夏季、秋季、春季。

    關(guān)鍵詞?冷島效應(yīng);大氣校正法;Landsat 8影像;莫莫格濕地

    中圖分類號?P?208;P?237文獻(xiàn)標(biāo)識碼?A

    文章編號?0517-6611(2019)23-0082-06

    doi:10.3969/j.issn.0517-6611.2019.23.026

    開放科學(xué)(資源服務(wù))標(biāo)識碼(OSID):

    Remote Sensing Analysis of Cold Island Effect in Momoge Wetland

    XIANG Xiao?yun1, LI Feng?xiu2 ,DU Jia3 et al

    (1.College of Tourist and Geoscience, Jilin Normal University, Siping,Jilin 136000;2.Henan Meteorological Center,Zhengzhou,Henan 450003;3.Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences, Changchun, Jilin 130102)

    Abstract?Based on the Landsat 8 image of the Momoge Wetland in the spring, summer and autumn of 2015, land surface temperature was estimated by using atmospheric correction method. We used buffering analysis, overlay analysis and geostatistical analysis in spatial analysis tools of ArcGIS software to complete the cold island effect study of Momoge Wetland based on the remote sensing inversion results and land use data. The results showed the cold island effect of the wetland on different land use types was gradually weakened when the distance from the wetland increased.Wetland had different cooling effects on different land use types. The impact on residence, dryland and saline?alkali land was roughly 240-540 m, 480-660 m and 600-660 m, respectively. Wetland has different cooling effects on other land use types in different seasons. In the three seasons of spring, summer and autumn, the cooling effect from strong to weak is summer, autumn and spring.

    Key words?Cold island effect;Atmospheric correction method;Landsat 8 image;Momoge Wetland

    濕地是位于陸生生態(tài)系統(tǒng)和水生生態(tài)系統(tǒng)之間的過渡性地帶,是水陸相互作用形成的一種獨特的生態(tài)系統(tǒng),濕地具有較大的比熱容量和較強的水分蒸散度,其自身的“碳匯”作用使其對周邊的大氣溫度和濕度都起到顯著的調(diào)節(jié)作用[1]。由于濕地長期和季節(jié)性積水,其水熱容量較大,氣候較周邊區(qū)域冷。正是由于這種冷島作用,所以濕地能夠在一定程度上對周邊不同土地利用類型具有冷島效應(yīng)。

    目前關(guān)于濕地的冷島效應(yīng)國內(nèi)有很多學(xué)者進(jìn)行了研究,楊一鵬等[2]提出氣候調(diào)節(jié)功能在干旱區(qū)的濕地表現(xiàn)尤為顯著,每年氣候調(diào)節(jié)所帶來的生態(tài)服務(wù)價值在濕地總生態(tài)服務(wù)價值中占比可達(dá)到32%,居于各項生態(tài)價值的首位;聶曉等[3]對沼澤濕地的局地小氣候效應(yīng)進(jìn)行研究,結(jié)果表明,沼澤濕地具有一定的冷濕效應(yīng),且在不同高度上冷濕作用表現(xiàn)有所不同,距離地面越近作用越顯著;拱秀麗等[4]以洪河自然保護(hù)區(qū)為研究區(qū),對沼澤濕地冷濕小氣候效應(yīng)的時空分布特征進(jìn)行模擬研究,發(fā)現(xiàn)濕地在垂直方向和水平方向上其冷濕效應(yīng)的效果存在差異。以往對濕地冷島效應(yīng)的研究多關(guān)注于濕地具有冷濕作用的功能研究,對于其冷島效應(yīng)影響程度和影響范圍的研究較少。

    遙感技術(shù)因其具備觀測范圍廣、速度快、周期短等優(yōu)勢,被廣泛地用于不同尺度的氣候變化研究。目前關(guān)于地表溫度的遙感反演算法有很多,包括輻射傳輸方程法[5-6]、單窗算法[7-9]、劈窗算法[10]以及多通道多角度算法[11-13],根據(jù)不同的研究內(nèi)容和不同的遙感傳感器數(shù)選擇相應(yīng)的算法[14-15]。1984年開始發(fā)展的Landsat遙感數(shù)據(jù),還有從2000年NASA開始提供MODIS熱紅外遙感影像,雖然空間分辨率較小,但在時間分辨率上具有很大的優(yōu)勢(同一地區(qū)一天可過境4次),因此采用熱紅外遙感數(shù)據(jù)可以大范圍定量地反演地表溫度?;谶b感技術(shù)的冷島效應(yīng)的研究正逐漸成為研究濕地冷島效應(yīng)的重要手段。

    長期以來,由于濕地本身的特性(低反照率、高發(fā)射率等),濕地的地表溫度遙感反演一直是定量遙感的熱點。葉智威等[16]選擇以 Landsat TM6 數(shù)據(jù)為基礎(chǔ),利用單窗算法對洪澤湖地區(qū)地表溫度進(jìn)行反演并分析,結(jié)果表明洪澤湖區(qū)陸地與水體溫度空間差異明顯,濕地在水體溫度中較低;杜嘉等[17]以Landsat/TM圖像、實測地面數(shù)據(jù)和MODIS地表發(fā)射數(shù)據(jù)為基礎(chǔ),運用3種不同算法對洪河濕地的地表溫度進(jìn)行估算,對比不同算法之間的差異,結(jié)果顯示使用覃志豪的單窗算法進(jìn)行估算的結(jié)果較為理想;杜培軍等[18]使用單窗算法和Landsat TM /ETM + 影像對江蘇省鹽城濱海濕地進(jìn)行地表溫度的反演,對1992—2009年濕地的土地利用變化和地表溫度變化之間的相關(guān)性進(jìn)行分析,結(jié)果顯示濕地土地利用變化會帶來地表溫度的波動。以往溫度反演研究多傾向于沿海地區(qū)和湖泊地區(qū)的地表溫度反演,對濕地類型的地表溫度反演的研究相對較少。

    該研究使用大氣校正法和Landsat 8 遙感數(shù)據(jù)為基礎(chǔ),以莫莫格濕地為研究區(qū)進(jìn)行地表溫度反演,利用地面實測數(shù)據(jù)和MOD11A1產(chǎn)品數(shù)據(jù)對其反演結(jié)果精度進(jìn)行驗證;利用ArcGIS軟件中強大的空間分析工具,包括緩沖區(qū)分析、疊加分析和地統(tǒng)計分析,對濕地冷島效應(yīng)對周邊不同土地利用類型的影響程度和影響范圍進(jìn)行進(jìn)一步的分析,目的在于通過定量的方法直觀地展示濕地的冷島效應(yīng),旨在加強人們對干旱-半干旱地區(qū)濕地功能的認(rèn)識。

    1?資料與方法

    1.1?研究區(qū)概況

    莫莫格濕地地理坐標(biāo)位于45°42′25″~46°18′00″N、123°27′00″~124°04′33″E,地處松嫩平原西部(圖1)。研究區(qū)內(nèi)主要以保護(hù)濕地生態(tài)系統(tǒng)和瀕危野生動物物種及其棲息環(huán)境為主,總面積14.4萬hm2,自然濕地面積在研究區(qū)總面積中占比80%以上,是吉林省最大的濕地保留地,也是松嫩平原保存最為完整的濕地之一[19]。松嫩平原區(qū)域內(nèi)有扎龍、莫莫格、查干湖、向海和大蘇布5個國家級濕地自然保護(hù)區(qū),其中莫莫格濕地作為典型的內(nèi)陸鹽堿化沼澤濕地且地處于內(nèi)陸干旱、半干旱過渡帶對當(dāng)?shù)厣鷳B(tài)環(huán)境具有重要意義。

    研究區(qū)的氣候類型是屬于溫帶大陸性季風(fēng)氣候。全年平均氣溫約為4.2 ℃,整年降水量約為391.8 mm,各季節(jié)分配不均,降水主要集中在6—8月份,占全年降水量的766%。研究區(qū)內(nèi)的地勢起伏較小且較為平坦,平均海拔約為142 m,總體呈現(xiàn)出西北高、東南低的趨勢,坡度一般為5°~15°,區(qū)域內(nèi)的相對海拔高差也在35 m以下。

    1.2?數(shù)據(jù)來源

    該研究中地表溫度反演所選用的遙感數(shù)據(jù)包括Landsat 8遙感影像和MOD11A1數(shù)據(jù)。Landsat 8遙感數(shù)據(jù)來源于美國地質(zhì)調(diào)查局USGS 網(wǎng)站(http://glovis.usgs.gov/)。Landsat 8衛(wèi)星于2013年2月11日在美國加州范登堡空軍基地發(fā)射成功,搭載的傳感器分別是OLI陸地成像儀和TIRS熱紅外傳感器,并于2013年5月提供全球免費下載。OLI傳感器在可見光、近紅外和短波紅外區(qū)域接收9個光譜波段的數(shù)據(jù),空間分辨率為30 m;TIRS傳感器將原來的熱紅外波段一分為二,設(shè)置成2個熱紅外通道(Band10:10.6~119 μm,Band11:11.50~12.51 μm),空間分辨率為100 m[20]。當(dāng)前熱紅外遙感技術(shù)可以獲取到地球表面熱紅外波譜的輻射能量,根據(jù)其地球表面不同物體的發(fā)射率特征可以反演出其熱力學(xué)溫度,所以在反演地表溫度時,熱紅外遙感影像可以提供長期的遙感影像支持。選擇影像的原則以能覆蓋研究區(qū)范圍且天氣晴朗、影像上基本無云或云覆蓋度在5%以下的影像。研究中使用的遙感影像如表1所示。

    MODSI LST數(shù)據(jù)主要來源為從LAAADS Web(http://ladsweb.modaps.eosdis.nasa.gov)選擇并下載。研究中所要使用到的MODIS LST地表溫度值是MODIS LST經(jīng)過處理的產(chǎn)品數(shù)據(jù),其原因在于可以更加便捷地通過處理獲取到相應(yīng)范圍內(nèi)的地表溫度值,并且可以提高所獲取到的地表溫度值的準(zhǔn)確性。MOD11A1數(shù)據(jù)獲取后,利用ERDAS 9.2和ArcGIS 10.3軟件進(jìn)行相應(yīng)的處理。MOD11A1產(chǎn)品數(shù)據(jù)在該研究中具有2個作用,第1個作用是可以驗證Landsat 8熱紅外遙感影像的反演精度,第2個作用是為研究中確定對Landsat 8反演結(jié)果分析的范圍提供一定的依據(jù)。

    在研究濕地冷島效應(yīng)過程中,除了使用Landsat 8遙感數(shù)據(jù)和MODSI LST產(chǎn)品(MOD11A1)外,還使用了土地利用數(shù)據(jù)做輔助數(shù)據(jù)。土地利用數(shù)據(jù)主要是采用人機交互的遙感解譯方法所獲取的。通過對遙感影像進(jìn)行輻射定標(biāo)、大氣校正及幾何校正等一系列處理,在ArcGIS軟件中進(jìn)行土地利用的解譯及完成矢量化,得到研究所需的土地利用數(shù)據(jù)。

    根據(jù)研究需求,將研究區(qū)內(nèi)土地利用類型主要劃分為居民地、旱地、林地、草地、鹽堿地、水田、水庫坑塘、河流和草本沼澤,共計分為9種類型。通過解譯和矢量化等一系列處理后,得到研究區(qū)2015年土地利用示意圖(圖2)。研究所得到的2015年各土地利用分類的Kappa系數(shù)為0.84,對應(yīng)的分類精度級別非常好,滿足進(jìn)一步研究需要。

    1.3?研究方法

    1.3.1?大氣校正法。

    輻射傳輸方程又稱大氣校正法,是地表溫度反演的基本方法之一,在所有熱紅外遙感波段中有著巨大的優(yōu)勢。在研究中使用Landsat 8 進(jìn)行地表溫度的反演,因為Landsat 8衛(wèi)星發(fā)射運行時間較短,TIRS第11波段暫時存在定標(biāo)不穩(wěn)定性[21],為了避免第10波段與第11波段所產(chǎn)生的較大誤差,因此采用大氣校正法和Landsat 8數(shù)據(jù)的第10波段對研究區(qū)地表溫度進(jìn)行反演。

    大氣校正法反演地表溫度的原理具體表現(xiàn)為:衛(wèi)星傳感器接收到的熱紅外輻射亮度值Lλ由大氣向上輻射亮度L↑、地面的真實輻射亮度經(jīng)過大氣層之后到達(dá)衛(wèi)星傳感器的能量以及大氣向下輻射到達(dá)地面后反射的能量共3個部分組成。衛(wèi)星傳感器接收到的熱紅外輻射亮度值Lλ的表達(dá)式如下所示:

    Lλ=[εB(TS)+(1-ε)L↓]τ+L↑(1)

    式中,ε為地表比輻射率,TS為地表真實溫度(K),B(TS)為黑體熱輻射亮度,τ為大氣在熱紅外波段的透過率。溫度為T的黑體在熱紅外波段的輻射亮度B(TS)為:

    B(TS)=[Lλ-L↑-τ(1-ε)L↓]/τε(2)

    TS可以用普朗克公式的函數(shù)獲?。?/p>

    TS=K2/ln(K1B(TS)+1)(3)

    對于TIRS Band10而言,K1 = 774.89 W/(m2·sr·μm),K2 = 1 321.08 K。

    1.3.2?驗證方法。

    根據(jù)已知的Landsat 8影像第10波段的空間分辨率為100 m,再以選取的研究區(qū)邊界為準(zhǔn),在ArcGIS軟件中生成長和寬均為100 m的正方形矢量格網(wǎng)圖層,如圖3所示。

    將得到的驗證格網(wǎng)矢量圖層分別疊加到對應(yīng)的Landsat 8反演結(jié)果遙感影像和MOD11A1影像中,使用ArcGIS中的分區(qū)統(tǒng)計功能將各個網(wǎng)格內(nèi)的溫度平均值提取出來,再通過Landsat 8反演結(jié)果與MOD11A1的地表溫度建立相關(guān)性并進(jìn)行對比分析,以達(dá)到驗證Landsat 8遙感影像反演地表溫度結(jié)果的精度。驗證過程中計算LST平均絕對誤差所采用的公式如下:

    MAD=|x-|n(4)

    式中,MAD為平均絕對誤差,x為數(shù)值,為數(shù)值平均值,n為數(shù)據(jù)個數(shù)。

    1.4?MODIS LST數(shù)據(jù)緩沖區(qū)分區(qū)統(tǒng)計分析結(jié)果

    獲取與3期研究區(qū)熱紅外遙感影像所對應(yīng)的MOD11A1地表溫度產(chǎn)品影像。MOD11A1是MODIS LST經(jīng)過處理的產(chǎn)品數(shù)據(jù),可以通過一系列處理獲取到相應(yīng)范圍內(nèi)的地表溫度值,并且可以提高所獲取到的地表溫度值的準(zhǔn)確性。MODIS過境時間為當(dāng)?shù)貢r間11:30,Landsat 8過境時刻為當(dāng)?shù)貢r間10:50左右,二者之間相差40 min,可以近似認(rèn)為是相同時刻過境,因此可以用MODIS數(shù)據(jù)對Landsat 8地表溫度反演結(jié)果進(jìn)行驗證。利用已獲取的MOD11A1數(shù)據(jù)進(jìn)行地表溫度緩沖區(qū)分區(qū)統(tǒng)計計算。

    在研究中,將研究區(qū)的所有濕地類型(草本沼澤、水田、河流、水庫坑塘)進(jìn)行合并,向外每1 km制作該區(qū)域3 km范圍的緩沖區(qū)。經(jīng)過分區(qū)統(tǒng)計后3期遙感影像對應(yīng)日期的結(jié)果如圖4所示。從圖4可以看出,緩沖區(qū)內(nèi)平均地表溫度從0 m開始不斷升高,在1~2 km平均地表溫度仍在不斷升高,但升溫幅度有所減緩,自3 km處開始溫度變化趨于平緩。但是因為MODIS LST產(chǎn)品的空間分辨率為1 km,所以對1 km范圍內(nèi)的溫度變化無法做詳細(xì)的分析,無法與中等分辨率影像匹配,這也是MODIS 遙感數(shù)據(jù)的局限性所在。由此可見,濕地冷島效應(yīng)在大約1 km范圍內(nèi)最為劇烈,為后續(xù)分析濕地冷島效應(yīng)的影響范圍提供一定依據(jù)。

    2?結(jié)果與分析

    2.1?Landsat 8反演地表溫度

    2.1.1?地表溫度反演結(jié)果驗證。

    在對Landsat 8地表溫度反演結(jié)果進(jìn)行精度驗證時,使用與遙感影像對應(yīng)日期的白城、大安、鎮(zhèn)賚、泰來和杜蒙5個氣象站點逐時的地面實測數(shù)據(jù)進(jìn)行進(jìn)一步驗證。驗證結(jié)果如圖5所示,通過驗證可以看出,Landsat 8的反演結(jié)果與地面實測數(shù)據(jù)溫度值的相關(guān)性R2值為0.948 9,溫度誤差在0.56~4.79 ℃,平均絕對誤差為2.5 ℃。造成這種誤差的原因在于Landsat 8的溫度反演結(jié)果是反演衛(wèi)星過境時刻的地表溫度值,但是驗證的氣象站點的溫度數(shù)據(jù)無法精確到與過境時刻一致,并且氣象站點的逐時溫度數(shù)據(jù)為這個時間段內(nèi)的溫度平均值,不能完全代表像元尺度的地表溫度,因此會有一定的溫度誤差。經(jīng)驗證,該反演結(jié)果可以達(dá)到研究的精度要求,可以進(jìn)行后續(xù)的分析。

    根據(jù)MODIS LST 產(chǎn)品具有每個像素地表溫度和發(fā)射率

    值的全球數(shù)據(jù)產(chǎn)品的準(zhǔn)確性,使用與遙感影像對應(yīng)日期的地表溫度實測數(shù)據(jù)以及MOD11A1數(shù)據(jù)進(jìn)行進(jìn)一步驗證。通過對研究區(qū)內(nèi)格網(wǎng)做地表溫度反演值的分區(qū)統(tǒng)計后,可以看出Landsat 8的反演結(jié)果與MOD11A1數(shù)據(jù)提取平均溫度值的相關(guān)性R2值均在0.6以上,且其平均絕對誤差也均小于2.4 K。產(chǎn)生這種誤差的原因在于Landsat 8與MOD11A1影像的空間分辨率不同,MOD11A1的空間分辨率為1 km,遠(yuǎn)高于Landsat系列傳感器的空間分辨率,研究區(qū)內(nèi)草地、林地和居住地等土地利用類型中存在面積較小或分布較為分散的現(xiàn)象,所以在1 km的MODIS產(chǎn)品數(shù)據(jù)上容易與其他土地利用方式形成混合像元由此帶來一定的誤差[22]。

    2.1.2?地表溫度反演結(jié)果分析。

    通過以3期Landsat 8遙感影像為基礎(chǔ),利用大氣校正法對3期影像進(jìn)行地表溫度反演,得到2015年春季、夏季、秋季3幅地表溫度反演圖(圖6)。從圖6可以看出,研究區(qū)內(nèi)溫度較低的區(qū)域均為濕地類型區(qū),以濕地類型為中心向外,周邊不同土地利用類型地表溫度逐漸升高。利用3期溫度反演結(jié)果和土地利用圖疊加分析可以發(fā)現(xiàn),不同土地利用類型的地表溫度并不一致。進(jìn)一步分析使用ArcGIS軟件計算不同土地利用類型的地表溫度,結(jié)果顯示各土地利用類型地表溫度由高到低排序依次為居民地類型、鹽堿地類型、旱地類型、濕地類型、林地、草地。

    根據(jù)上述,研究區(qū)范圍內(nèi)所得到的反演結(jié)果均在正常值以內(nèi),莫莫格濕地研究區(qū)內(nèi)所反演的最低地表溫度>-4 ℃,最高地表溫度<44 ℃,這與之前在獲取反演參數(shù)的計算及選取上較為精準(zhǔn),反演結(jié)果符合自然規(guī)律。

    2.2?莫莫格濕地地表溫度的空間分布

    MODIS LST數(shù)據(jù)緩沖區(qū)分區(qū)統(tǒng)計分析結(jié)果顯示,在1 km范圍內(nèi)濕地對周邊區(qū)域降溫作用較為明顯。所以利用Landsat 8遙感影像反演地表溫度的結(jié)果和土地利用數(shù)據(jù)進(jìn)行疊加分析,對莫莫格濕地進(jìn)行緩沖區(qū)分區(qū)統(tǒng)計計算。由于研究區(qū)內(nèi)林地和草地的面

    積所占比例很?。值卣急?.06%,草地占比0.31%),所以在對不同土地利用類型溫度反演分析中將其去除。

    從研究區(qū)各土地利用LST分區(qū)統(tǒng)計(圖7)可以看出,3月研究區(qū)內(nèi)地表溫度由高到低依次為旱地類型、居住地類型、鹽堿地類型,在6月和10月結(jié)果中研究區(qū)內(nèi)地表溫度由高到低依次為居住地類型、旱地類型、鹽堿地類型。綜合來看,除3月份以外,居住地類型的地表溫度均值均高于旱地類型,這可能是由于3月份研究區(qū)旱地的農(nóng)業(yè)生產(chǎn)并未展開,其實際上為裸土地,6月和10月地表有農(nóng)作物覆蓋,而植被覆蓋能夠有效降低其地表溫度[23],所以3月旱地類型的地表溫度均值較高。

    從圖7可以看出,在莫莫格濕地緩沖區(qū)范圍內(nèi),濕地類型的地表溫度較低,周邊各地物類型的地表溫度均高于其地表溫度。這主要是因為濕地內(nèi)部長期和季節(jié)性積水的特性,使其具有較大的比熱容量,因此其地表溫度相對周邊區(qū)域溫度更低。周邊各土地利用類型中,居住地類型和旱地類型的地表溫度一般要高于鹽堿地類型的地表溫度。居住地類型主要是由于其獨特的下墊面性質(zhì)(居民地一般是水泥混凝土或者磚瓦構(gòu)造),因此其地表溫度較高。

    總體而言,莫莫格濕地研究區(qū)中的居住地、旱地和鹽堿地土地利用類型的地表溫度總體趨勢表現(xiàn)為:隨著距離濕地類型距離的逐漸增加,其地表溫度隨之開始升高。這主要是因為濕地地區(qū)獨特的水熱條件使其地表溫度較低,對各周邊地物來說具有冷島效應(yīng),因此對周邊不同土地利用類型有著降低其地表溫度的作用。

    2.3?濕地冷島效應(yīng)影響范圍分析

    在分析研究區(qū)濕地類型對周邊各土地利用類型地表溫度的影響范圍中,根據(jù)之前緩沖區(qū)分區(qū)統(tǒng)計分析結(jié)果,計算研究區(qū)范圍內(nèi)各緩沖距離之間的不同土地利用類型的地表溫度斜率。然后選擇各距離斜率接近于0值的距離,作為濕地類型影響周邊不同土地利用類型地表溫度的范圍。在上述濕地地表溫度的空間分布分析中由于林地和草地所占面積比例較小將其排除在外,所以在影響范圍的研究中同樣將二者排除在外,主要研究居住地、旱地以及鹽堿地3種土地利用類型受濕地冷島效應(yīng)的影響。

    通過“2.2”分析可以發(fā)現(xiàn),濕地類型對周邊各地物具有降溫作用,但降溫作用具有一定的范圍。因此,在各土地利用類型地表溫度變化中,對研究區(qū)內(nèi)濕地類型對不同土地利用類型地表溫度影響范圍(降溫范圍)進(jìn)行進(jìn)一步分析,結(jié)果如表2所示。從表2可以看出,在莫莫格濕地研究區(qū)內(nèi),濕地類型對居住地類型地表溫度的影響邊界范圍在240~540 m,在0~540 m居住地類型的地表溫度隨著距濕地類型越遠(yuǎn),其地表溫度開始上升,在540~990 m居住地類型的地表溫度隨著距濕地類型越遠(yuǎn),其地表溫度受濕地類型的影響開始減弱;濕地類型對旱地類型地表溫度的影響邊界范圍在480~660 m,在0~660 m旱地類型的地表溫度隨著距濕地類型越遠(yuǎn),其地表溫度逐漸升高,在660~990 m旱地類型的地表溫度受濕地類型的影響開始減弱;濕地類型對鹽堿地類型地表溫度的影響邊界范圍在600~660 m,在0~660 m鹽堿地類型的地表溫度距濕地類型越遠(yuǎn),其地表溫度逐漸升高,在660~990 m受濕地類型的影響開始減弱。從不同季節(jié)角度分析可以發(fā)現(xiàn),不同季節(jié)濕地類型對周邊不同土地利用類型的影響范圍各不相同,春季影響范圍在330~600 m;夏季影響范圍在540~660 m;秋季影響范圍在240~660 m。

    綜上所述,可以發(fā)現(xiàn)濕地類型對周邊地區(qū)不同土地利用類型的降溫效果各不相同,并且在不同季節(jié)內(nèi)濕地類型對不同土地利用類型的影響范圍也不相同??傮w來看,濕地類型對鹽堿地類型的降溫效果最為明顯,其次是旱地類型和居住地類型,這主要是因為在同樣氣象條件下,鹽堿地由于比熱更小地表溫度更高,更易與周邊濕地發(fā)生能量交換;鹽堿地比居民地和旱地的植被更少,導(dǎo)致其相應(yīng)的空氣動力學(xué)阻抗更小,也導(dǎo)致其更易受到濕地的影響。在春夏秋3個季節(jié)中,夏季的濕地冷島效應(yīng)最強,對周邊地區(qū)的降溫范圍最大。這主要是由于濕地冷島效應(yīng)影響范圍與濕地通量貢獻(xiàn)區(qū)面積直接相關(guān),而通量貢獻(xiàn)區(qū)面積大小受大氣穩(wěn)定度、觀測高度以及下墊面粗糙度等因素的影響[24-25]。在大氣穩(wěn)定狀態(tài)下,大氣湍流運動較弱,氣體垂直擴散速度較緩慢,所以通量信息可以代表較遠(yuǎn)的區(qū)域;而在不穩(wěn)定狀態(tài)下,大氣湍流運動劇烈,氣體垂直擴散速度較快,因此通量信息代表的區(qū)域較近[26]。該研究區(qū)中,夏季相對春秋兩季,大氣狀態(tài)更加穩(wěn)定,通量貢獻(xiàn)區(qū)面積更大,因此夏季的濕地冷島效應(yīng)更強。

    3?結(jié)論與討論

    利用大氣校正法對2015年3期Landsat 8影像進(jìn)行地表溫度反演工作,通過后續(xù)對反演結(jié)果進(jìn)行空間分析等進(jìn)一步研究,發(fā)現(xiàn)濕地的冷島效應(yīng)對周邊不同土地利用類型有著一定程度的降溫作用,旨在人們對濕地冷島效應(yīng)的重要功能進(jìn)行更深層次的認(rèn)識。

    首先,通過對3期地表溫度反演結(jié)果進(jìn)行地表溫度空間分布分析,發(fā)現(xiàn)在研究區(qū)范圍內(nèi),除濕地類型外(濕地地表溫度最低),總體來看居住地和旱地類型地表溫度均值較高,鹽堿地類型地表溫度均值較低。其次,通過對Landsat 8遙感影像反演地表溫度的結(jié)果和土地利用數(shù)據(jù)進(jìn)行疊加分析,結(jié)果表明濕地對周邊各土地利用類型的影響受距離遠(yuǎn)近的影響,在一定范圍內(nèi)隨著距離濕地距離的增加其地表溫度不斷升高。

    最后,對研究區(qū)內(nèi)各緩沖距離之間的各土地利用類型的地表溫度之間計算出斜率,進(jìn)一步分析研究區(qū)內(nèi)濕地類型對周邊不同土地利用類型地表溫度的影響范圍,結(jié)果表明,濕地對不同土地利用類型的降溫作用和降溫范圍各不相同,降溫作用由強到弱依次為居住地類型、旱地類型、鹽堿地類型;在不同季節(jié)中,濕地對周邊各土地利用類型的降溫作用也不相同,降溫作用由強到弱依次為夏季、春季、秋季。

    盡管遙感反演溫度的精度驗證結(jié)果可以達(dá)到研究的精度要求,但是其反演精度仍存在進(jìn)步的空間,在今后的研究中通過利用不同算法和高分辨率遙感影像來進(jìn)一步提高反演精度。在研究中,利用大氣校正法和Landsat 8遙感數(shù)據(jù)對莫莫格濕地冷島效應(yīng)進(jìn)行分析,得到較為理想的結(jié)果。

    但濕地冷島效應(yīng)在空間上是連續(xù)的,遙感影像在空間上是間斷、離散的,不能完美地展示其過程,需要用更高分辨率的影像和氣象模型相結(jié)合,對濕地冷島效應(yīng)過程進(jìn)一步分析。

    參考文獻(xiàn)

    [1] SUN R H,CHEN A L,CHEN L D,et al.Cooling effects of wetlands in an urban region:The case of Beijing[J].Ecological indicators,2012,20:57-64.

    [2] 楊一鵬,曹廣真,侯鵬,等.城市濕地氣候調(diào)節(jié)功能遙感監(jiān)測評估[J].地理研究,2013,32(1):73-80.

    [3] 聶曉,王毅勇.沼澤濕地局地小氣候 “冷濕島” 效應(yīng)[J].生態(tài)與農(nóng)村環(huán)境學(xué)報,2010,26(2):189-192.

    [4] 拱秀麗,王毅勇,聶曉,等.沼澤濕地與周邊旱田氣溫,相對濕度差異分析[J].東北林業(yè)大學(xué)學(xué)報,2011,39(11):93-96.

    [5] SOBRINO J A,JIMNEZ?MUOZA J C,PAOLINIB L.Land surface temperature retrieval from Landsat TM5[J].Remote sensing of environment,2004,90(4):434-440.

    [6] LI Z L,TANG B H,WU H,et al.Satellite?derived land surface temperature:Current status and perspectives[J].Remote sensing of environment,2013,131:14-37.

    [7] 覃志豪,ZHANG M H,ARNON KARNIELI,et al.用陸地衛(wèi)星 TM6 數(shù)據(jù)演算地表溫度的單窗算法[J].地理學(xué)報,2001,56(4):456-466.

    [8] JIMNEZ?MUOZ J C,SOBRINO J A.A generalized single?channel method for retrieving land surface temperature from remote sensing data[J].Journal of geophysical research,2003,108(D22):1-9.

    [9] 周紀(jì),李京,趙祥,等.用HJ-1B 衛(wèi)星數(shù)據(jù)反演地表溫度的修正單通道算法[J].紅外與毫米波學(xué)報,2011,30(1):61-67.

    [10] 楊槐.從 Landsat8 影像反演地表溫度的劈窗算法研究[J].測繪地理信息,2014,39(4):73-77.

    [11] SOBRINO J A,LI Z L,STOLL M P,et al.Multi?channel and multi?angle algorithms for estimating sea and land surface temperature with ASTR data[J].International journal of remote sensing,1996,17(11):2089-2114.

    [12] GILLESPIE A,ROKUGAWA S,MATSUNAGA T,et al.A temperature and emissivity separation algorithm for Advanced Spaceborne Thermal Emission and Reflection Radiometer(ASTER)images[J].IEEE Transactions on Geoscience and Remote Sensing,1998,36(4):1113-1126.

    [13] 毛克彪,施建成,覃志豪,等.一個針對 ASTER 數(shù)據(jù)同時反演地表溫度和比輻射率的四通道算法[J].遙感學(xué)報,2006,10(4):593-599.

    [14]DU C,REN H Z,QIN Q M,et al.A practical split?window algorithm for estimating land surface temperature from Landsat 8 data[J].Remote sensing,2015,7(1):647-665.

    [15] 胡德勇,喬琨,王興玲,等.單窗算法結(jié)合 Landsat8 熱紅外數(shù)據(jù)反演地表溫度[J].遙感學(xué)報,2015,19(6):964-976.

    [16] 葉智威,覃志豪,宮輝力.洪澤湖區(qū)的Landsat TM6 地表溫度遙感反演和空間差異分析[J].首都師范大學(xué)學(xué)報(自然科學(xué)版),2009(1):88-95.

    [17] 杜嘉,張柏,宋開山,等.基于 Landsat-5TM 的洪河濕地地表溫度估算方法對比研究[J].遙感技術(shù)與應(yīng)用,2009,24(3):312-319.

    [18] 杜培軍,陳宇,譚琨.江蘇濱海濕地土地利用/覆蓋變化與地表溫度響應(yīng)遙感監(jiān)測[J].國土資源遙感,2014,26(2):112-120.

    [19] 于秀麗.松嫩平原莫莫格濕地土壤不同形態(tài)碳轉(zhuǎn)化過程及其與土壤鐵耦合關(guān)系研究[D].長春:東北師范大學(xué),2016.

    [20] IRONS J R,DWYER J L,BARSI J A.The next Landsat satellite:The Landsat Data Continuity Mission[J].Remote sensing of environment,2012,122:11-21.

    [21] USGS.2014.January 29,2014?Landsat 8 Reprocessing to Begin February 3,2014[EB/OL].[2015-04-26].http://landsat.usgs.gov /calibration_notices.

    [22] MA J,ZHOU J,LIU S M,et al.Review on validation of remotely sensed land surface temperature[J].Advances in earth science,2017,32(6):615-629.

    [23] 楊士弘.城市生態(tài)環(huán)境學(xué)[M].北京:科學(xué)出版社,2003:114-119.

    [24] 米娜,于貴瑞,溫學(xué)發(fā),等.中國通量觀測網(wǎng)絡(luò)(ChinaFLUX)通量觀測空間代表性初步研究[J].中國科學(xué)(D輯),2006,36(S1):22-33.

    [25] 金瑩,張志強,方顯瑞,等.楊樹人工林生態(tài)系統(tǒng)通量貢獻(xiàn)區(qū)分析[J].生態(tài)學(xué)報,2012,32(12):3966-3974.

    [26] 王江濤,周劍虹,歐強,等.崇明東灘濱海圍墾濕地CO2通量貢獻(xiàn)區(qū)分析[J].生態(tài)與農(nóng)村環(huán)境學(xué)報,2014,30(5):588-594.

    猜你喜歡
    冷島反演土地利用
    反演對稱變換在解決平面幾何問題中的應(yīng)用
    極熱條件下空冷島永磁電動機散熱能力研究
    高寒條件下空冷島風(fēng)機集群永磁直驅(qū)系統(tǒng)批量改造分析
    高寒條件下空冷島變頻永磁直驅(qū)改造及應(yīng)用
    空冷島的熱態(tài)沖洗及防凍措施
    基于低頻軟約束的疊前AVA稀疏層反演
    土地利用生態(tài)系統(tǒng)服務(wù)研究進(jìn)展及啟示
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    濱??h土地利用挖潛方向在哪里
    疊前同步反演在港中油田的應(yīng)用
    av又黄又爽大尺度在线免费看 | 亚洲欧美日韩高清专用| 精品免费久久久久久久清纯| 色哟哟·www| 日韩制服骚丝袜av| 亚洲精品亚洲一区二区| 在线观看66精品国产| 日韩国内少妇激情av| 国产伦精品一区二区三区四那| 99热这里只有是精品50| 高清毛片免费看| 欧美极品一区二区三区四区| 国产精品久久久久久精品电影小说 | 国产精品久久久久久av不卡| 午夜亚洲福利在线播放| 亚洲欧洲日产国产| АⅤ资源中文在线天堂| 国产精品一区二区三区四区久久| 亚洲成人中文字幕在线播放| 欧美日本亚洲视频在线播放| 国产爱豆传媒在线观看| 日本-黄色视频高清免费观看| 成年版毛片免费区| 黄色配什么色好看| 亚洲无线观看免费| 美女被艹到高潮喷水动态| 精品久久久噜噜| 国产av一区在线观看免费| 国产精品蜜桃在线观看 | 联通29元200g的流量卡| 国产又黄又爽又无遮挡在线| 亚洲久久久久久中文字幕| 51国产日韩欧美| 一进一出抽搐动态| 久久久久网色| 一本一本综合久久| 免费观看的影片在线观看| 精品久久久久久久末码| 欧美在线一区亚洲| 免费人成视频x8x8入口观看| 成人午夜高清在线视频| 一夜夜www| 亚洲欧美日韩卡通动漫| 免费观看人在逋| 亚洲av第一区精品v没综合| 亚洲欧洲日产国产| 国产成人a区在线观看| 国产亚洲av片在线观看秒播厂 | 国产精品99久久久久久久久| 国产精品人妻久久久久久| 麻豆av噜噜一区二区三区| 性色avwww在线观看| 欧美精品国产亚洲| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 色综合色国产| 国产精品国产高清国产av| 国产精品乱码一区二三区的特点| 亚洲欧美日韩东京热| 免费看光身美女| 欧美激情国产日韩精品一区| 日本熟妇午夜| 岛国毛片在线播放| 国产毛片a区久久久久| 国语自产精品视频在线第100页| 日本熟妇午夜| 国产精品乱码一区二三区的特点| 国产精品一区二区三区四区免费观看| 看片在线看免费视频| 赤兔流量卡办理| 国产成人午夜福利电影在线观看| 国产欧美日韩精品一区二区| 国产精品电影一区二区三区| 国产色婷婷99| 日韩大尺度精品在线看网址| 天堂av国产一区二区熟女人妻| 免费一级毛片在线播放高清视频| 免费观看精品视频网站| 午夜视频国产福利| 国产探花在线观看一区二区| a级一级毛片免费在线观看| 人人妻人人看人人澡| 在线播放无遮挡| 国产伦在线观看视频一区| 欧美激情国产日韩精品一区| 久久久久国产网址| 国产精品一区二区三区四区免费观看| 久久久久久久久久久免费av| 成人亚洲欧美一区二区av| 亚洲经典国产精华液单| 一级黄色大片毛片| 成人特级黄色片久久久久久久| 一个人看的www免费观看视频| 欧美不卡视频在线免费观看| 女同久久另类99精品国产91| 在线观看免费视频日本深夜| 亚洲第一区二区三区不卡| 成人毛片60女人毛片免费| 22中文网久久字幕| 欧美日本亚洲视频在线播放| eeuss影院久久| 国产蜜桃级精品一区二区三区| 国产成年人精品一区二区| 99在线人妻在线中文字幕| 成年女人永久免费观看视频| or卡值多少钱| av在线播放精品| 天天躁日日操中文字幕| 亚洲va在线va天堂va国产| 亚洲欧美中文字幕日韩二区| 少妇人妻精品综合一区二区 | 看片在线看免费视频| 中文资源天堂在线| 久久久久久久久久久免费av| 日本-黄色视频高清免费观看| 久久热精品热| 亚洲精品成人久久久久久| avwww免费| 国产伦精品一区二区三区视频9| 精品一区二区免费观看| 深夜精品福利| 一级毛片电影观看 | 亚洲精品久久久久久婷婷小说 | 亚洲国产日韩欧美精品在线观看| 免费看光身美女| 高清午夜精品一区二区三区 | 欧美xxxx黑人xx丫x性爽| 男人的好看免费观看在线视频| 国产真实乱freesex| 非洲黑人性xxxx精品又粗又长| 女的被弄到高潮叫床怎么办| 亚洲图色成人| 99精品在免费线老司机午夜| 欧美+日韩+精品| 夫妻性生交免费视频一级片| 欧美精品一区二区大全| 18禁黄网站禁片免费观看直播| 免费电影在线观看免费观看| 久久精品影院6| 欧美性猛交╳xxx乱大交人| 91在线精品国自产拍蜜月| 精品熟女少妇av免费看| 国产成人精品久久久久久| 一区二区三区高清视频在线| 人人妻人人澡欧美一区二区| 18禁黄网站禁片免费观看直播| 在线免费观看的www视频| 国产一区二区亚洲精品在线观看| 舔av片在线| 两个人视频免费观看高清| 亚洲精品国产av成人精品| 亚洲人与动物交配视频| 蜜桃亚洲精品一区二区三区| 老女人水多毛片| av在线亚洲专区| 久久久久久久久中文| 日韩国内少妇激情av| 日本成人三级电影网站| 久久久久久九九精品二区国产| 亚洲成人中文字幕在线播放| 美女黄网站色视频| 九九热线精品视视频播放| 亚洲五月天丁香| 麻豆成人午夜福利视频| 中出人妻视频一区二区| 国产 一区 欧美 日韩| 美女被艹到高潮喷水动态| 久久久久久久久久成人| 你懂的网址亚洲精品在线观看 | 91狼人影院| 在线观看66精品国产| 国产乱人视频| 热99在线观看视频| 天堂√8在线中文| 成人午夜高清在线视频| 精品一区二区三区视频在线| 性欧美人与动物交配| 尤物成人国产欧美一区二区三区| 精品国内亚洲2022精品成人| 伦精品一区二区三区| 国产亚洲精品久久久久久毛片| 国内少妇人妻偷人精品xxx网站| 国产 一区精品| av免费在线看不卡| 欧美性猛交黑人性爽| 美女内射精品一级片tv| 99在线视频只有这里精品首页| 一卡2卡三卡四卡精品乱码亚洲| 日韩视频在线欧美| 两性午夜刺激爽爽歪歪视频在线观看| 日本黄色视频三级网站网址| 国产一区二区三区av在线 | 免费看av在线观看网站| 老女人水多毛片| 精品少妇黑人巨大在线播放 | 精品少妇黑人巨大在线播放 | 欧美最新免费一区二区三区| 国产免费一级a男人的天堂| 97人妻精品一区二区三区麻豆| 精品一区二区免费观看| or卡值多少钱| 久久国内精品自在自线图片| 日本黄色片子视频| 精品久久久噜噜| 日韩,欧美,国产一区二区三区 | 蜜桃久久精品国产亚洲av| 欧美一级a爱片免费观看看| 99久久人妻综合| 国产激情偷乱视频一区二区| 欧美不卡视频在线免费观看| 国产极品天堂在线| 97超视频在线观看视频| 丝袜美腿在线中文| 午夜激情福利司机影院| 亚洲图色成人| 国产精品99久久久久久久久| 亚洲精品亚洲一区二区| 国产精品久久久久久精品电影| 激情 狠狠 欧美| 免费观看的影片在线观看| 欧美极品一区二区三区四区| 搞女人的毛片| 久久欧美精品欧美久久欧美| 欧美精品国产亚洲| 99国产精品一区二区蜜桃av| 午夜激情欧美在线| 又黄又爽又刺激的免费视频.| 老熟妇乱子伦视频在线观看| 国产三级在线视频| 日韩三级伦理在线观看| 美女内射精品一级片tv| 亚洲欧美精品综合久久99| 国产成人午夜福利电影在线观看| 小说图片视频综合网站| av免费观看日本| 欧美精品国产亚洲| 97热精品久久久久久| 中文字幕免费在线视频6| av专区在线播放| 午夜福利在线观看免费完整高清在 | 国产高清视频在线观看网站| 69av精品久久久久久| 少妇熟女欧美另类| 国内精品久久久久精免费| .国产精品久久| av免费观看日本| 乱人视频在线观看| 亚洲成人久久爱视频| 久久精品国产亚洲av天美| 内射极品少妇av片p| 国产精品一区www在线观看| 国产日韩欧美在线精品| 插逼视频在线观看| 婷婷六月久久综合丁香| 国产精品女同一区二区软件| 在线免费十八禁| 国产乱人偷精品视频| 国产精品国产高清国产av| 2021天堂中文幕一二区在线观| 国产一区二区在线av高清观看| 亚洲国产精品成人综合色| 欧美3d第一页| 在线国产一区二区在线| 身体一侧抽搐| 九九久久精品国产亚洲av麻豆| 综合色丁香网| 久99久视频精品免费| 又黄又爽又刺激的免费视频.| av福利片在线观看| 亚洲欧美成人综合另类久久久 | 麻豆国产97在线/欧美| 插阴视频在线观看视频| 国产精品蜜桃在线观看 | 两个人的视频大全免费| 久久这里只有精品中国| 天天躁日日操中文字幕| 欧美3d第一页| 国产成人一区二区在线| 亚洲av中文av极速乱| 久久精品影院6| 精品久久久久久久末码| 国产精品蜜桃在线观看 | 国产极品精品免费视频能看的| 国产高清激情床上av| 亚洲国产精品成人综合色| 成人二区视频| 亚洲欧美成人精品一区二区| 久久热精品热| 黑人高潮一二区| 高清毛片免费观看视频网站| 欧洲精品卡2卡3卡4卡5卡区| .国产精品久久| 热99在线观看视频| 一个人观看的视频www高清免费观看| 国产精品伦人一区二区| 亚洲无线观看免费| 中文欧美无线码| 噜噜噜噜噜久久久久久91| 国产精品一区www在线观看| 欧美最黄视频在线播放免费| 国产老妇女一区| 久久亚洲精品不卡| 亚洲精品影视一区二区三区av| 一卡2卡三卡四卡精品乱码亚洲| 在线播放无遮挡| 99国产极品粉嫩在线观看| a级毛片免费高清观看在线播放| 99久国产av精品国产电影| 欧美一级a爱片免费观看看| 欧美日韩综合久久久久久| 午夜a级毛片| 亚洲精品影视一区二区三区av| 精品久久久久久成人av| 日韩,欧美,国产一区二区三区 | 国产极品天堂在线| 国国产精品蜜臀av免费| 久久人人爽人人片av| 午夜免费男女啪啪视频观看| 久久久久久九九精品二区国产| 国产精品三级大全| 午夜爱爱视频在线播放| 99久久精品一区二区三区| 1000部很黄的大片| 国产一区二区在线av高清观看| 国产黄色视频一区二区在线观看 | 国产日韩欧美在线精品| 全区人妻精品视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 寂寞人妻少妇视频99o| 日韩欧美三级三区| 一本久久精品| 国产综合懂色| 午夜福利在线观看免费完整高清在 | 淫秽高清视频在线观看| 午夜激情福利司机影院| 久久久久久久久中文| 成人毛片60女人毛片免费| 欧美精品国产亚洲| 国产人妻一区二区三区在| 中文字幕久久专区| 国产精品综合久久久久久久免费| 国产亚洲91精品色在线| 人妻少妇偷人精品九色| 91精品一卡2卡3卡4卡| 午夜福利在线在线| 99九九线精品视频在线观看视频| 久久99热6这里只有精品| 看黄色毛片网站| 男人舔奶头视频| 听说在线观看完整版免费高清| 啦啦啦啦在线视频资源| 日本熟妇午夜| 亚洲欧美成人综合另类久久久 | 精品午夜福利在线看| 国产片特级美女逼逼视频| 午夜视频国产福利| 免费看a级黄色片| 我要搜黄色片| 高清毛片免费看| 一级黄色大片毛片| 久久久久久久久久黄片| 亚洲欧美中文字幕日韩二区| 国产成人午夜福利电影在线观看| 日日啪夜夜撸| 日韩成人av中文字幕在线观看| 免费黄网站久久成人精品| 国产高清不卡午夜福利| 国产一区亚洲一区在线观看| 亚洲国产色片| 人人妻人人看人人澡| av国产免费在线观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲图色成人| 免费观看在线日韩| 床上黄色一级片| 久久人人爽人人片av| 内射极品少妇av片p| 日韩精品有码人妻一区| 成人无遮挡网站| 国产色爽女视频免费观看| 三级国产精品欧美在线观看| 国产一区二区激情短视频| 麻豆av噜噜一区二区三区| 中文字幕av在线有码专区| 婷婷色综合大香蕉| 欧美成人一区二区免费高清观看| 乱码一卡2卡4卡精品| 日韩一本色道免费dvd| 青春草亚洲视频在线观看| 精品99又大又爽又粗少妇毛片| 日本成人三级电影网站| 男女边吃奶边做爰视频| 久久精品国产99精品国产亚洲性色| 国语自产精品视频在线第100页| 深夜精品福利| 麻豆精品久久久久久蜜桃| 国产精品.久久久| 国产成年人精品一区二区| 热99re8久久精品国产| 九草在线视频观看| 亚洲18禁久久av| 国产69精品久久久久777片| 久久久久久久久久久免费av| 亚洲美女视频黄频| 欧美色视频一区免费| 色5月婷婷丁香| 国产成人91sexporn| 可以在线观看的亚洲视频| 久久鲁丝午夜福利片| 亚洲av.av天堂| 午夜激情欧美在线| 岛国毛片在线播放| av又黄又爽大尺度在线免费看 | 婷婷六月久久综合丁香| 免费观看的影片在线观看| 91精品一卡2卡3卡4卡| 亚洲欧美成人综合另类久久久 | 嫩草影院精品99| 亚洲精品乱码久久久v下载方式| 高清午夜精品一区二区三区 | 99热这里只有是精品50| 国产三级中文精品| av视频在线观看入口| 日韩视频在线欧美| 亚洲欧美日韩无卡精品| 亚洲五月天丁香| 亚洲av中文字字幕乱码综合| 久久精品久久久久久噜噜老黄 | 久久午夜福利片| 欧美变态另类bdsm刘玥| 欧美日韩精品成人综合77777| 最后的刺客免费高清国语| 人人妻人人看人人澡| 日韩制服骚丝袜av| 亚洲国产精品成人久久小说 | 哪里可以看免费的av片| 亚洲精品456在线播放app| 精品不卡国产一区二区三区| 波多野结衣高清作品| 欧美性感艳星| 亚洲国产精品久久男人天堂| 一区二区三区高清视频在线| 欧美在线一区亚洲| 啦啦啦观看免费观看视频高清| 亚洲真实伦在线观看| 能在线免费看毛片的网站| 国产精品人妻久久久久久| 丰满的人妻完整版| 国产日本99.免费观看| 国产白丝娇喘喷水9色精品| 亚洲精品粉嫩美女一区| 欧美成人精品欧美一级黄| 久久久久久久久久久免费av| 国产国拍精品亚洲av在线观看| 99久国产av精品| 我的老师免费观看完整版| 日日干狠狠操夜夜爽| 色尼玛亚洲综合影院| 一进一出抽搐gif免费好疼| 欧美在线一区亚洲| 精品久久久久久久人妻蜜臀av| 黄片wwwwww| 欧美色视频一区免费| 婷婷亚洲欧美| 国产私拍福利视频在线观看| 亚洲人与动物交配视频| 欧美日韩综合久久久久久| 99九九线精品视频在线观看视频| 国产精品永久免费网站| 国产亚洲5aaaaa淫片| 99久久精品国产国产毛片| 人妻制服诱惑在线中文字幕| 精品午夜福利在线看| 中文字幕免费在线视频6| 久久欧美精品欧美久久欧美| 亚洲va在线va天堂va国产| 国产成人一区二区在线| 长腿黑丝高跟| 欧洲精品卡2卡3卡4卡5卡区| 此物有八面人人有两片| 国产亚洲精品久久久com| 九草在线视频观看| 男女边吃奶边做爰视频| 日韩欧美国产在线观看| 国产精品福利在线免费观看| 日韩亚洲欧美综合| 白带黄色成豆腐渣| 99九九线精品视频在线观看视频| 级片在线观看| 亚洲精品日韩av片在线观看| 国产老妇伦熟女老妇高清| 老熟妇乱子伦视频在线观看| 激情 狠狠 欧美| 啦啦啦观看免费观看视频高清| 看非洲黑人一级黄片| 少妇人妻精品综合一区二区 | 久久欧美精品欧美久久欧美| 国内精品宾馆在线| 国内精品一区二区在线观看| 91久久精品国产一区二区三区| 久久精品国产鲁丝片午夜精品| 日韩欧美精品v在线| 亚洲欧美中文字幕日韩二区| 少妇高潮的动态图| 深爱激情五月婷婷| 国产av麻豆久久久久久久| 国产精品,欧美在线| 国产单亲对白刺激| 亚洲国产高清在线一区二区三| 亚洲av第一区精品v没综合| 国产探花在线观看一区二区| 男女做爰动态图高潮gif福利片| 麻豆乱淫一区二区| 一区二区三区高清视频在线| 男女视频在线观看网站免费| 免费看光身美女| 欧美性猛交黑人性爽| 十八禁国产超污无遮挡网站| 色噜噜av男人的天堂激情| 中文字幕精品亚洲无线码一区| 国产午夜福利久久久久久| 青青草视频在线视频观看| 精品99又大又爽又粗少妇毛片| 22中文网久久字幕| 观看美女的网站| 精品久久久久久久久久久久久| 亚洲精品久久久久久婷婷小说 | av免费观看日本| 国产精品一二三区在线看| 三级国产精品欧美在线观看| 在线观看午夜福利视频| 一个人看视频在线观看www免费| 欧美性猛交黑人性爽| 美女被艹到高潮喷水动态| 久久中文看片网| 国产av一区在线观看免费| 成人三级黄色视频| 少妇熟女欧美另类| ponron亚洲| 久久午夜亚洲精品久久| 免费av不卡在线播放| 日日摸夜夜添夜夜爱| 国产免费男女视频| 中国国产av一级| 看黄色毛片网站| 又粗又硬又长又爽又黄的视频 | 岛国毛片在线播放| 99九九线精品视频在线观看视频| 久久这里有精品视频免费| 国产一区亚洲一区在线观看| 国产大屁股一区二区在线视频| 久久精品国产自在天天线| 又爽又黄无遮挡网站| 日本与韩国留学比较| 亚洲人与动物交配视频| 我的老师免费观看完整版| 春色校园在线视频观看| 在线a可以看的网站| 国产精品久久视频播放| 免费无遮挡裸体视频| 性插视频无遮挡在线免费观看| 色综合色国产| 亚洲aⅴ乱码一区二区在线播放| 国产探花极品一区二区| 青青草视频在线视频观看| 亚洲精华国产精华液的使用体验 | 一本精品99久久精品77| 丰满人妻一区二区三区视频av| 尤物成人国产欧美一区二区三区| 九九久久精品国产亚洲av麻豆| 日本三级黄在线观看| 亚洲av.av天堂| 乱码一卡2卡4卡精品| 精品一区二区免费观看| 高清在线视频一区二区三区 | 亚洲欧美成人精品一区二区| 国产一区二区三区在线臀色熟女| 成人午夜精彩视频在线观看| 少妇的逼好多水| 国产探花在线观看一区二区| 在线免费观看的www视频| 观看免费一级毛片| 色吧在线观看| 又黄又爽又刺激的免费视频.| 日韩欧美三级三区| 又黄又爽又刺激的免费视频.| 国产一区亚洲一区在线观看| 天天躁日日操中文字幕| 三级男女做爰猛烈吃奶摸视频| 午夜亚洲福利在线播放| 日韩国内少妇激情av| 亚洲国产精品成人综合色| 看十八女毛片水多多多| 卡戴珊不雅视频在线播放| 国产白丝娇喘喷水9色精品| 亚洲精品久久国产高清桃花| 国产亚洲精品久久久久久毛片| 99久久精品国产国产毛片| 国产成人精品久久久久久| 精品一区二区三区视频在线| 一进一出抽搐gif免费好疼| 国产美女午夜福利| 波多野结衣高清无吗| 午夜激情欧美在线| 最好的美女福利视频网| 成人鲁丝片一二三区免费| 婷婷色综合大香蕉| 校园春色视频在线观看| 搡女人真爽免费视频火全软件| 别揉我奶头 嗯啊视频| 99热网站在线观看| 久久久a久久爽久久v久久| 国产精品一区二区三区四区久久| 啦啦啦韩国在线观看视频| 欧美性猛交黑人性爽| 国产伦理片在线播放av一区 |