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

    陜西省土壤濕度時空變化特征及其驅(qū)動力分析

    2023-11-02 11:47:08王征
    城市勘測 2023年5期
    關(guān)鍵詞:土壤濕度陜西省風(fēng)速

    王征

    (中煤航測遙感集團(tuán)有限公司,陜西 西安 710018)

    0 引 言

    土壤濕度(Soil Moisture)是表示土壤干旱程度的重要指標(biāo),是陸-氣系統(tǒng)水熱過程的重要紐帶,同時也是植被生長的主要水分來源[1]。目前主要通過兩種手段獲得土壤濕度信息[2],第一種是田間實測法,即利用野外觀測站或者采樣點直接測量不同深度的土壤濕度值,該方法雖然可以準(zhǔn)確地獲得土壤濕度信息,但是獲得的基于站點尺度數(shù)據(jù),不能夠刻畫大尺度、大范圍內(nèi)的土壤水分時空變化。第二種是采用遙感數(shù)據(jù)對土壤濕度信息進(jìn)行反演,是當(dāng)前較為常見的方法。通常結(jié)合熱慣量法、植被指數(shù)法等方法反演土壤濕度。馬春峰等人[3]以MODIS為數(shù)據(jù)源,采用熱慣量法對黑河流域土壤濕度信息進(jìn)行估測,但該方法對于日較差數(shù)據(jù)精度有較高要求,無法滿足遙感定量研究的需要。2003年姚春生研究發(fā)現(xiàn)[4],土壤水分決定著植被冠層溫度,而冠層溫度間接決定土壤供水情況,基于此,提出溫度植被干旱指數(shù)(TVDI)。至此,國內(nèi)外學(xué)者基于TVDI監(jiān)測土壤濕度開展大量研究。Cao等人[5]利用AVHRR和MODIS數(shù)據(jù)構(gòu)建了Ts-NDVI特征空間監(jiān)測蒙古高原土壤濕度并分析干旱的時空變化特征,研究發(fā)現(xiàn)Ts-NDVI特征空間是監(jiān)測區(qū)域干旱的有效方法。Zhang等人[6]利用NOAA/AVHRR數(shù)據(jù)提取TVDI對2006年夏季四川盆地土壤濕度進(jìn)行監(jiān)測,研究發(fā)現(xiàn)TVDI是監(jiān)測區(qū)域干旱的有效手段。程夢園等人[7]基于LST和NDVI數(shù)據(jù)構(gòu)建Ts-NDVI特征空間對2005—2020年香日德-柴達(dá)木河流域的土壤濕度信息進(jìn)行提取,研究發(fā)現(xiàn)Ts-NDVI特征空間反演該流域土壤濕度具有較好的適用性,同時TVDI與海拔表現(xiàn)負(fù)相關(guān),而與平均氣溫表現(xiàn)正相關(guān)趨勢??祱虻热薣8]基于MODIS數(shù)據(jù)中的NDVI和LST數(shù)據(jù)構(gòu)建NDVI-LST特征空間對2000—2019年蒙古高原干旱程度進(jìn)行監(jiān)測,研究發(fā)現(xiàn)該地區(qū)的干旱情況小幅度加重,其增長速率為 0.000 1/a。眾多研究表明應(yīng)用TVDI可以實現(xiàn)土壤濕度的反演,但已有研究在長時序土壤濕度動態(tài)監(jiān)測和驅(qū)動因素方面仍然存在不足,理論情況下TVDI特征空間能夠達(dá)到裸土到全植被覆蓋的地表情況,但在實際應(yīng)用中,研究區(qū)域地表覆蓋情況難以滿足,得到的干濕邊實際上只是特征空間內(nèi)部存在的邊,并非理論的干濕邊。陜西省是重要的糧食作物產(chǎn)區(qū),而土壤濕度作為農(nóng)作物生產(chǎn)過程中的核心要素,需要長時序動態(tài)土壤濕度觀測數(shù)據(jù),但是陜西省氣象站點分布較少,難以獲取土壤濕度信息。

    因此,本研究基于2000—2020年NDVI和LST數(shù)據(jù)構(gòu)建改進(jìn)的雙拋物線NDVI-LST特征空間,通過對TVDI指數(shù)進(jìn)行計算從而監(jiān)測陜西省土壤濕度情況,并結(jié)合氣候因素(降水、氣溫、日照時數(shù)和風(fēng)速)探討與其相關(guān)關(guān)系,以期為陜西省旱情監(jiān)測提供決策支持。

    1 研究區(qū)概況

    陜西省位于我國西北部地區(qū),地理范圍為105°29′~111°15′E,31°42′~39°35′N之間[9],總面積約為2.0×105km2(圖1)。陜西省區(qū)域內(nèi)地形地貌豐富多樣,地勢表現(xiàn)為南北高,中間低,平均海拔約為 1 100 m。

    圖1 研究區(qū)行政區(qū)劃圖

    陜西省地處濕潤地區(qū)和干旱區(qū)的交界帶,陜南為濕潤區(qū),主要有秦嶺和大巴山系,海拔高度在200~1 200 m之間,水熱條件相對較好,年均氣溫約為14℃,年降水量達(dá)到 800 mm;關(guān)中位于中間,屬于半濕潤區(qū),主要是由渭河干流及兩岸支流沉積形成,地勢較低,海拔高度僅在300~800 m內(nèi),年均氣溫約為12℃,年降水量達(dá)到 600 mm;陜北則為半干旱區(qū),包括榆林和延安地區(qū),海拔高度在900~1 600 m之間,年均氣溫約為9℃,年降水量約為 500 mm。與此同時,陜西省地處長江流域中上游地區(qū),屬于國家生態(tài)環(huán)境重點建設(shè)區(qū)域,但是陜西省也是中國水土流失問題較為嚴(yán)重的地區(qū)之一,因此,監(jiān)測區(qū)域土壤濕度時空變化對于陸地表面能量平衡具有重要作用[10]。

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

    2.1 數(shù)據(jù)源及預(yù)處理

    本研究采用的遙感數(shù)據(jù)為美國航天航空局的MODIS數(shù)據(jù)產(chǎn)品(https://modis.gsfc.nasa.gov/data/),獲取數(shù)據(jù)為2000—2020年,空間分辨率為 1 km×1 km,時間分辨率為 8 d的MOD11A2的合成產(chǎn)品和 16 d的MOD13A2合成產(chǎn)品。利用MRT軟件進(jìn)行一系列預(yù)處理操作,如拼接、投影轉(zhuǎn)換和格式轉(zhuǎn)換等[11],為消除云、霧以及大氣氣溶膠等外界因素的影響,本研究采用最大值合成法(MVC)計算得到月最大NDVI和月最大值LST,最后采用Python進(jìn)行裁剪處理。

    氣象數(shù)據(jù)采用中國氣象數(shù)據(jù)網(wǎng)(http://data.cma.cn/)2000—2020年陜西省及相鄰31個氣象站點的逐日氣溫、降水量、平均風(fēng)速和日照時數(shù)數(shù)據(jù),采用MATLAB軟件對該數(shù)據(jù)集進(jìn)行清洗和異常值處理,對于空缺數(shù)據(jù)采用雙線性內(nèi)插方法對其進(jìn)行填補[12],最后采用ANUSPLIN插值方法對其進(jìn)行空間插值,得到2000—2020年逐年 1 km分辨率的年均氣溫(Tmean)、總降水量(PRE)、年均平均風(fēng)速(WIN)和年日照總時數(shù)(SSD)數(shù)據(jù)。

    2.2 溫度植被干旱指數(shù)法(Ts-NDVI特征空間)

    研究發(fā)現(xiàn)[13],當(dāng)研究區(qū)植被覆蓋達(dá)到從裸土到全植被覆蓋、土壤含水量達(dá)到從完全干旱到濕潤的條件下,NDVI和地表溫度的分布圖表現(xiàn)出類四邊形形狀(如圖2所示)。四邊形四個頂點分別代表4種極限情況,AD為干邊,土壤含水量有效性最小,地表水分蒸發(fā)量也較小;BC為濕邊,所包含的土壤含水量最大,此時地表蒸散量屬于潛在蒸散量。從A到B點,地表土壤水分蒸發(fā)急劇增大,從B到C點,植被覆蓋度由裸土達(dá)到全植被覆蓋。

    圖2 Ts-NDVI特征空間原理

    一般來說,當(dāng)監(jiān)測區(qū)完全達(dá)到濕潤狀態(tài)時,地表土壤內(nèi)的水分也隨之偏大,單一時段下特征空間濕邊比理論值高;當(dāng)監(jiān)測區(qū)完全達(dá)到干旱狀態(tài)時,地表土壤內(nèi)的水分也隨之偏小,其特征空間干邊比理論值低;為彌補這一問題,本研究采用長時序遙感數(shù)據(jù)對陜西省土壤濕度進(jìn)行監(jiān)測,以此彌補單一時段下Ts-NDVI特征空間的弊端。序列數(shù)據(jù)越完整,其合成的特征空間就更符合理論特征空間。計算公式如下[14]:

    (1)

    式中:LST表示像元地表溫度;LSTmin和LSTmax分別表示相同NDVI值所對應(yīng)地表溫度的最小值和最大值;TVDI在0~1之間,其值越大表示越干旱。

    本研究采用對干濕邊方程改進(jìn)后的特征空間來計算TVDI,其計算公式如下[15]:

    LSTmin=a1+b1×NDVI+c1×NDVI2

    (2)

    LSTmin=a2+b2×NDVI+c2×NDVI2

    (3)

    其中:a1、b1、c1、a2、b2和c2分別為干濕邊方程的系數(shù)。

    本文將TVDI作為衡量干旱程度的指標(biāo),為消除量綱,對NDVI和LST均歸一化處理,步長為0.01,編寫代碼獲取同一NDVI下全部單元格對應(yīng)的LST最大值和最小值,從而確定Ts-NDVI特征空間圖,將不同時期特征空間的干濕邊進(jìn)行線性擬合,最終得到逐年TVDI值。

    2.3 Theil-Sen Median趨勢分析和Mann-Kendall檢驗

    基于Theil-Sen Median趨勢分析結(jié)合Mann-Kendall趨勢檢驗方法能夠分析長時間序列數(shù)據(jù),該方法相較于傳統(tǒng)的一元線性回歸分析方法而言,具有無可比擬的優(yōu)點,如無數(shù)據(jù)情況下能夠服從特定分布、可抵抗一定程度的數(shù)據(jù)誤差,其計算公式為[16]:

    (4)

    式中:Median為中值函數(shù);TVDIi、TVDIj分別為該像元在i和j年的TVDI數(shù)據(jù)。β>0,即TVDI呈上升狀態(tài);β=0,即TVDI呈穩(wěn)定狀態(tài);β>0,即TVDI呈下降狀態(tài)。

    Mann-Kendall是一種非參數(shù)突變檢驗方法,它能夠判斷趨勢的顯著程度,當(dāng)n>10時,標(biāo)準(zhǔn)的正態(tài)系統(tǒng)變量通過下式計算[17]:

    (5)

    2.4 Hurst指數(shù)

    Hurst指數(shù)是用于表征長時序變量在未來變化的演變趨勢,該方法已在水文、植被、氣象等領(lǐng)域得到廣泛應(yīng)用[18]。目前,Hurst指數(shù)常使用重標(biāo)極差(R/S),以此確定陜西省未來的干旱變化情況。有研究發(fā)現(xiàn)[6],采用R/S分析法得到的Hurst指數(shù)更真實,計算公式為

    均值序列:

    (6)

    累計離差:

    (7)

    極差:

    (8)

    標(biāo)準(zhǔn)差:

    (9)

    根據(jù)Hurst指數(shù)的H值能夠確定陜西省的TVDI時間序列的持續(xù)性。Hurst指數(shù)主要分為三種情況:當(dāng)0.5

    3 結(jié)果分析

    3.1 基于TVDI的陜西省干旱與氣候因素的時間變化特征

    2000—2020年陜西省年均TVDI介于0.48~0.6之間,總體呈現(xiàn)波動下降趨勢,其增長速率為-0.002/a(p<0.01),這表明陜西省土壤濕度整體呈現(xiàn)濕潤狀態(tài)。其中,2010年TVDI均值最小,為0.492;2005年TVDI均值最大,為0.585。研究時段可以2010年為界限,劃分2000—2010年的由干往濕發(fā)展和2010—2020年的由濕往干發(fā)展。2000—2010年年均TVDI呈現(xiàn)下降趨勢,其增長速率為-0.004/a;而2010—2020年年均TVDI呈現(xiàn)上升趨勢,其增長速率為0.002/a。

    對于氣候因素而言,平均氣溫、總降水、平均風(fēng)速和日照總時數(shù)總體呈現(xiàn)波動上升的趨勢,其增長速率分別為0.004/a、1.732/a、0.005/a和4.527/a(圖3)。

    圖3 TVDI與氣溫、降水、風(fēng)速和日照時數(shù)的時間變化趨勢

    3.2 基于TVDI的陜西省干旱空間變化特征

    由于陜西省地區(qū)氣候類型多樣,各地區(qū)植被覆蓋度也存在明顯不同,因此年均TVDI值不足以表示干旱空間變化特征。本研究結(jié)合Theil-Sen Median趨勢分析與Mann-Kendall檢驗方法進(jìn)行空間疊加,對陜西省內(nèi)像元尺度下的旱情趨勢進(jìn)行分析(圖4左)。首先對所有像元的年均TVDI值進(jìn)行趨勢分析,獲取到趨勢值β的空間分布,再利用檢驗方法得到TVDI變化顯著性空間分布,將計算得到趨勢值和顯著性進(jìn)行疊加分類,得到陜西省年均TVDI的5種類型變化趨勢,如顯著改善、輕微改善、基本穩(wěn)定、輕微退化和嚴(yán)重退化。結(jié)果發(fā)現(xiàn),陜西省的干旱變化中,面積占比最大的類型是輕微改善(50.91%),占比最小的類型是嚴(yán)重退化(0.95%),這表明2000—2020年陜西省土壤濕度在空間上往濕潤方向發(fā)展為主,其中,呈現(xiàn)持續(xù)改善的地區(qū)主要分布在延安市的吳起縣、安塞縣、延川縣、甘泉縣,渭南市的合陽縣、蒲城縣,寶雞市的隴縣、千陽縣,西安市的藍(lán)田縣、戶縣,漢中市的寧強(qiáng)縣等地區(qū),主要占陜西省總面積的16.31%。呈現(xiàn)退化趨勢的面積占陜西省總面積的18.13%,其中呈現(xiàn)嚴(yán)重退化的面積占比0.95%,主要分布在榆林市的定邊縣,寶雞市的太白縣等地區(qū)。

    圖4 2000—2020年陜西省TVDI時空變化趨勢及未來趨勢分析

    為揭示陜西省土壤濕度未來變化的持續(xù)性特征,本研究進(jìn)一步將趨勢值與Hurst指數(shù)進(jìn)行耦合,得到陜西省年均TVDI的5種類型未來變化趨勢,包括持續(xù)下降、下降、上升、持續(xù)上升、穩(wěn)定不變(圖4右)。陜西省年均TVDI的Hurst指數(shù)值為0.441,大于0.5的像元數(shù)所占比例為25.70%,小于0.5的像元數(shù)所占比例為74.30%。結(jié)果表明陜西省大部分地區(qū)TVDI未來變化趨勢與過去相反。在空間分布上,陜西省有52.37%的區(qū)域未來干旱情況有好轉(zhuǎn)趨勢。陜北地區(qū)植被以低矮灌木為主,生態(tài)環(huán)境相對較為脆弱,所以干旱的未來趨勢無法確定。蒸發(fā)量、降水量和CO2濃度等因素都是影響植被生理效應(yīng)的主要因素,但是這些因素在未來存在較強(qiáng)的不確定性,這就使得未來干旱演變軌跡很難監(jiān)測。但是本文通過對TVDI的趨勢值和Hurst指數(shù)進(jìn)行耦合綜合分析,實現(xiàn)了陜西省未來干旱情況的監(jiān)測。

    3.3 陜西省TVDI與氣候因素相關(guān)性分析

    氣溫、降水、日照時數(shù)和風(fēng)速變化都是影響土壤濕度的重要影響因子,因此,本研究基于陜西省的氣象數(shù)據(jù)進(jìn)一步結(jié)合陜西省TVDI時序數(shù)據(jù)分析其與降水、氣溫、日照時數(shù)和風(fēng)速的相關(guān)性(圖5)。結(jié)果發(fā)現(xiàn),陜西省TVDI與總降水之間以負(fù)相關(guān)關(guān)系為主(-0.109),這表面在降水量增加的同時,TVDI值變小,土壤濕度呈現(xiàn)增加趨勢。根據(jù)相關(guān)系數(shù)可看出,負(fù)相關(guān)占陜西省總面積的70.86%,主要分布在榆林市的橫山縣、神木縣、佳縣、米脂縣,漢中市的城固縣、寧強(qiáng)縣等地區(qū),降水是土壤水分的重要補給,能夠緩解干旱的加重,當(dāng)降水量逐漸減少時,土壤中的水分含量也在逐漸減少,特別是淺層土壤。正相關(guān)占陜西省總面積的29.14%,主要分布在咸陽市的乾縣、彬縣、永壽縣、禮泉縣,渭南市的白水縣等地區(qū),隨著降水量的增加,植被蒸騰作用導(dǎo)致植物需水量大幅度增加,也間接性導(dǎo)致這些地區(qū)土壤水分消耗加快,使得TVDI與降水出現(xiàn)了正相關(guān)關(guān)系。

    圖5 2000—2020年陜西省TVDI與降水、氣溫、風(fēng)速和日照時數(shù)之間的空間相關(guān)性

    陜西省TVDI與平均氣溫主要呈正相關(guān)(0.003),這表明當(dāng)氣溫升高時,TVDI值變大,則土壤濕度逐漸降低。根據(jù)相關(guān)系數(shù)結(jié)果可知,正相關(guān)占陜西省總面積的50.85%,主要分布在榆林市的子洲縣、橫山縣、定邊縣,延安市的甘泉縣、安塞縣、子長縣、延長縣,渭南市的大荔縣、澄城縣等地區(qū),當(dāng)溫度高于植被適宜值時,植被葉面蒸騰作用以及土壤水分蒸散量加快,使得該地區(qū)具有干旱化趨勢。負(fù)相關(guān)占陜西省總面積的49.15%,主要分布在漢中市的佛坪縣,西安市的周至縣、戶縣等地區(qū),適宜的溫度使得地表土壤水分不易于蒸發(fā),植被鎖水能力增強(qiáng),使得該地區(qū)愈發(fā)濕潤。

    陜西省TVDI與平均風(fēng)速之間主要呈正相關(guān)(0.006),這表明當(dāng)風(fēng)速增加時,TVDI值變大,則土壤濕度逐漸降低。根據(jù)相關(guān)系數(shù)結(jié)果可知,正相關(guān)占陜西省總面積的52.92%,主要分布于榆林市的子洲縣、神木縣、定邊縣,延安市的甘泉縣、延長縣、子長縣,渭南市的大荔縣、澄城縣等地區(qū)。大風(fēng)會帶走土壤中的水分,使得該地區(qū)具有干旱化趨勢。負(fù)相關(guān)占陜西省總面積的47.08%,主要分布于漢中市佛坪縣,西安市戶縣等地區(qū),適當(dāng)?shù)娘L(fēng)速不僅帶來植被生長所需要的二氧化碳,而且加快了植被葉面蒸騰促進(jìn)了植物根部吸收土壤水分。

    陜西省TVDI與日照總時數(shù)主要呈正相關(guān)(0.118),這說明隨著日照時數(shù)的增加,TVDI值變大,則土壤濕度逐漸降低。根據(jù)相關(guān)系數(shù)結(jié)果可知,正相關(guān)占陜西省總面積的71.58%,主要分布于榆林市的橫山縣、定邊縣、子洲縣,延安市和渭南市的大部分地區(qū)。隨著日照時間的增長,加快土壤水分的增加,再加之這些地區(qū)本身干旱缺水。負(fù)相關(guān)占陜西省總面積的28.42%,主要分布于漢中市的佛坪縣,西安市的戶縣等地區(qū),適宜的日照可以促進(jìn)植被光合作用,有利于植被的生長,植被鎖水能力增強(qiáng),使得該地區(qū)愈發(fā)濕潤。

    4 結(jié) 論

    (1)從時間演變特征上分析1,陜西省21年的年均TVDI在0.48~0.6波動,增速為-0.002/a(p<0.01),整體呈現(xiàn)波動下降趨勢;由Theil-Sen Median趨勢分析與Mann-Kendall檢驗疊加分析表明,空間上陜西省主要以輕微改善為主,說明土壤濕度在空間上往濕潤方向發(fā)展為主。

    (2)在未來演變特征上,陜西省年均TVDI表現(xiàn)出較強(qiáng)的反向持續(xù)性,其中Hurst指數(shù)小于0.5的地區(qū)占74.30%,表明陜西省大部分地區(qū)TVDI具有相反的變化趨勢。通過采用Theil-Sen Median趨勢分析與Hurst指數(shù)結(jié)果疊加分析表明,陜西省大部分地區(qū)未來干旱情況總體呈現(xiàn)好轉(zhuǎn)趨勢。

    (3)從TVDI與氣候因素相關(guān)性分析結(jié)果來看,TVDI與降水呈負(fù)相關(guān),與氣溫、風(fēng)速和日照時數(shù)呈正相關(guān),并且TVDI與降水、氣溫、風(fēng)速的相關(guān)系數(shù)分別為-0.109、0.003、0.006和0.118。

    猜你喜歡
    土壤濕度陜西省風(fēng)速
    陜西省自強(qiáng)中等專業(yè)學(xué)校簡介
    基于Kmeans-VMD-LSTM的短期風(fēng)速預(yù)測
    基于最優(yōu)TS評分和頻率匹配的江蘇近海風(fēng)速訂正
    海洋通報(2020年5期)2021-01-14 09:26:54
    土壤濕度傳感器在園林綠化灌溉上的應(yīng)用初探
    基于51單片機(jī)控制花盆土壤濕度
    電子制作(2019年15期)2019-08-27 01:12:12
    陜西省抓黨建促脫貧攻堅的實踐與思考
    聚焦兩會
    陜西畫報(2018年1期)2018-11-17 19:33:14
    陜西省閱讀文化節(jié)
    基于GARCH的短時風(fēng)速預(yù)測方法
    四川盆地土壤濕度時空分布及影響因子分析
    久久毛片免费看一区二区三区| 精品熟女少妇av免费看| 欧美亚洲 丝袜 人妻 在线| 欧美激情国产日韩精品一区| av有码第一页| av播播在线观看一区| 日本午夜av视频| 色吧在线观看| 国产日韩一区二区三区精品不卡| 看免费av毛片| 又黄又粗又硬又大视频| 香蕉精品网在线| 国产精品人妻久久久久久| 丝袜脚勾引网站| 在线看a的网站| 亚洲国产精品专区欧美| 国产免费视频播放在线视频| 亚洲av综合色区一区| 午夜免费观看性视频| 免费女性裸体啪啪无遮挡网站| 搡女人真爽免费视频火全软件| 久久99热这里只频精品6学生| 国产精品一二三区在线看| av视频免费观看在线观看| 男女免费视频国产| 国产精品人妻久久久影院| 国产精品久久久久久av不卡| 中文字幕制服av| 亚洲成av片中文字幕在线观看 | 免费av不卡在线播放| 观看av在线不卡| 成人18禁高潮啪啪吃奶动态图| 国产 精品1| 国产色婷婷99| 永久网站在线| 日韩三级伦理在线观看| 日韩 亚洲 欧美在线| 日韩视频在线欧美| 日本黄大片高清| 亚洲精品乱久久久久久| videossex国产| 久久久久久久大尺度免费视频| 人人妻人人爽人人添夜夜欢视频| 91aial.com中文字幕在线观看| 国产片内射在线| 人妻少妇偷人精品九色| 久久狼人影院| 精品酒店卫生间| 国产精品国产三级专区第一集| 色哟哟·www| 色婷婷久久久亚洲欧美| 一区二区日韩欧美中文字幕 | 人人妻人人添人人爽欧美一区卜| 婷婷色av中文字幕| 久久久精品免费免费高清| 亚洲精品aⅴ在线观看| 亚洲国产欧美在线一区| 欧美xxxx性猛交bbbb| 两个人免费观看高清视频| 国产免费一级a男人的天堂| 亚洲成国产人片在线观看| 亚洲精品国产av成人精品| 亚洲国产精品国产精品| 国产成人精品久久久久久| 18禁裸乳无遮挡动漫免费视频| 久久精品国产综合久久久 | 午夜激情久久久久久久| av不卡在线播放| 丝袜喷水一区| 亚洲国产精品一区二区三区在线| 国产亚洲欧美精品永久| 亚洲精品第二区| 一区二区日韩欧美中文字幕 | 我要看黄色一级片免费的| 新久久久久国产一级毛片| 欧美日韩综合久久久久久| 精品一区二区免费观看| 2021少妇久久久久久久久久久| 老女人水多毛片| 色网站视频免费| 亚洲综合色惰| 午夜老司机福利剧场| 十八禁高潮呻吟视频| 99国产综合亚洲精品| 七月丁香在线播放| 亚洲精品国产av成人精品| 欧美精品高潮呻吟av久久| 国产免费福利视频在线观看| 激情五月婷婷亚洲| 婷婷成人精品国产| 九草在线视频观看| 成人国语在线视频| 狂野欧美激情性bbbbbb| 亚洲精品美女久久久久99蜜臀 | 亚洲美女搞黄在线观看| 久久人妻熟女aⅴ| 777米奇影视久久| 久久久久网色| 日韩一区二区视频免费看| 亚洲精品视频女| 欧美精品av麻豆av| 久久久久久人人人人人| 人体艺术视频欧美日本| 日日撸夜夜添| 男人爽女人下面视频在线观看| 日韩视频在线欧美| 香蕉国产在线看| 免费观看性生交大片5| 内地一区二区视频在线| 成人毛片60女人毛片免费| av一本久久久久| 欧美丝袜亚洲另类| 99精国产麻豆久久婷婷| 又粗又硬又长又爽又黄的视频| 精品一区二区三区视频在线| 免费在线观看完整版高清| 69精品国产乱码久久久| 高清欧美精品videossex| 男女无遮挡免费网站观看| 伊人亚洲综合成人网| av天堂久久9| 亚洲精品中文字幕在线视频| 波野结衣二区三区在线| 久久久精品94久久精品| 亚洲国产成人一精品久久久| 日本-黄色视频高清免费观看| 国产精品免费大片| 国产精品久久久久久精品古装| 性高湖久久久久久久久免费观看| 成人国语在线视频| 天堂俺去俺来也www色官网| 日韩,欧美,国产一区二区三区| 国产欧美日韩综合在线一区二区| av不卡在线播放| 考比视频在线观看| 搡女人真爽免费视频火全软件| 夜夜骑夜夜射夜夜干| 狂野欧美激情性bbbbbb| 国产白丝娇喘喷水9色精品| 亚洲一级一片aⅴ在线观看| 欧美日韩av久久| 蜜桃在线观看..| 男人操女人黄网站| 久久国内精品自在自线图片| 亚洲精品久久成人aⅴ小说| 日产精品乱码卡一卡2卡三| 久久97久久精品| 亚洲精品一区蜜桃| 日韩一本色道免费dvd| 伦理电影免费视频| 乱人伦中国视频| 久久人妻熟女aⅴ| 激情五月婷婷亚洲| 高清在线视频一区二区三区| 精品福利永久在线观看| 国产1区2区3区精品| 精品亚洲成国产av| 日产精品乱码卡一卡2卡三| 一二三四中文在线观看免费高清| 久久久久久久亚洲中文字幕| 国产精品偷伦视频观看了| 国产1区2区3区精品| 少妇精品久久久久久久| 午夜av观看不卡| 好男人视频免费观看在线| 国产色婷婷99| 91aial.com中文字幕在线观看| 亚洲婷婷狠狠爱综合网| 日韩电影二区| av在线观看视频网站免费| av在线播放精品| 精品亚洲成国产av| 两性夫妻黄色片 | 国产伦理片在线播放av一区| 国产成人免费观看mmmm| 热re99久久精品国产66热6| 亚洲成人av在线免费| 精品久久国产蜜桃| 日韩制服丝袜自拍偷拍| 国产精品国产三级专区第一集| 黄色毛片三级朝国网站| 国产午夜精品一二区理论片| 成人国产麻豆网| 国产老妇伦熟女老妇高清| 97超碰精品成人国产| 亚洲av成人精品一二三区| 国产精品嫩草影院av在线观看| 久久久欧美国产精品| 国产一区二区在线观看日韩| 中文乱码字字幕精品一区二区三区| 日韩成人伦理影院| 国产男人的电影天堂91| 成人二区视频| 午夜福利在线观看免费完整高清在| 男女高潮啪啪啪动态图| 乱码一卡2卡4卡精品| 亚洲精品视频女| 国产免费视频播放在线视频| 考比视频在线观看| 一本—道久久a久久精品蜜桃钙片| 51国产日韩欧美| 成人免费观看视频高清| 国产亚洲最大av| 免费黄网站久久成人精品| 人成视频在线观看免费观看| 国产精品久久久久久av不卡| 97超碰精品成人国产| 精品久久久久久电影网| 婷婷色av中文字幕| 国产永久视频网站| 亚洲欧美清纯卡通| 欧美激情极品国产一区二区三区 | 精品亚洲成a人片在线观看| 色94色欧美一区二区| 日本av免费视频播放| 久久青草综合色| 久久久国产欧美日韩av| 午夜av观看不卡| 在现免费观看毛片| 人人妻人人澡人人爽人人夜夜| 色哟哟·www| 久久人人97超碰香蕉20202| 天天操日日干夜夜撸| 久久久久久久精品精品| 日韩视频在线欧美| 人人妻人人澡人人看| 黄色视频在线播放观看不卡| 色94色欧美一区二区| 国产视频首页在线观看| 国产一区有黄有色的免费视频| 超色免费av| 又黄又粗又硬又大视频| 久久久国产精品麻豆| 国产精品三级大全| 99久久人妻综合| 国产极品天堂在线| 久久久久精品性色| 亚洲成人av在线免费| 青春草亚洲视频在线观看| 老女人水多毛片| 国产一区二区三区av在线| 亚洲av免费高清在线观看| 在线免费观看不下载黄p国产| 少妇人妻精品综合一区二区| 美女中出高潮动态图| 国产女主播在线喷水免费视频网站| 各种免费的搞黄视频| 中文字幕精品免费在线观看视频 | 亚洲欧美成人综合另类久久久| 午夜福利视频精品| 飞空精品影院首页| 九草在线视频观看| 国产极品粉嫩免费观看在线| 久久97久久精品| 99久久人妻综合| 全区人妻精品视频| 国产一区亚洲一区在线观看| 国产成人aa在线观看| 日韩熟女老妇一区二区性免费视频| 97在线视频观看| 极品少妇高潮喷水抽搐| 亚洲熟女精品中文字幕| 建设人人有责人人尽责人人享有的| 国产精品久久久久久久电影| 日日爽夜夜爽网站| 亚洲成av片中文字幕在线观看 | 久久人人爽av亚洲精品天堂| av视频免费观看在线观看| av片东京热男人的天堂| 天天躁夜夜躁狠狠躁躁| 只有这里有精品99| 国产成人aa在线观看| 99视频精品全部免费 在线| 黄片播放在线免费| 9色porny在线观看| 久久这里有精品视频免费| 国产成人91sexporn| 丰满饥渴人妻一区二区三| 国产免费一级a男人的天堂| 嫩草影院入口| 捣出白浆h1v1| 中文字幕最新亚洲高清| 欧美变态另类bdsm刘玥| 日韩av免费高清视频| 高清毛片免费看| 久久99蜜桃精品久久| 欧美日韩av久久| 99视频精品全部免费 在线| xxxhd国产人妻xxx| av视频免费观看在线观看| 国产欧美另类精品又又久久亚洲欧美| 乱人伦中国视频| 亚洲熟女精品中文字幕| 欧美丝袜亚洲另类| 亚洲精品,欧美精品| 亚洲国产精品成人久久小说| 肉色欧美久久久久久久蜜桃| av卡一久久| 99re6热这里在线精品视频| 精品久久国产蜜桃| 夜夜爽夜夜爽视频| 各种免费的搞黄视频| 欧美成人午夜免费资源| 亚洲美女黄色视频免费看| 国产精品一区www在线观看| 免费观看av网站的网址| 91在线精品国自产拍蜜月| 在线观看一区二区三区激情| 成人黄色视频免费在线看| 国产免费福利视频在线观看| 国产一区二区激情短视频 | 亚洲成av片中文字幕在线观看 | 国产黄色免费在线视频| 在线观看国产h片| 久久久久久人妻| 看免费成人av毛片| 亚洲av男天堂| 久久国内精品自在自线图片| 成年女人在线观看亚洲视频| 亚洲精品久久久久久婷婷小说| 亚洲成人手机| 国产色婷婷99| 亚洲精品乱码久久久久久按摩| av在线播放精品| 又黄又爽又刺激的免费视频.| 1024视频免费在线观看| 久久97久久精品| 国产一区有黄有色的免费视频| 国产极品天堂在线| 美女内射精品一级片tv| 日韩视频在线欧美| 夫妻性生交免费视频一级片| 国产精品免费大片| 晚上一个人看的免费电影| 亚洲成色77777| 亚洲内射少妇av| 久久久久久久大尺度免费视频| 超色免费av| 久久精品夜色国产| 1024视频免费在线观看| 国产 精品1| 纵有疾风起免费观看全集完整版| 9热在线视频观看99| 日本-黄色视频高清免费观看| 日韩欧美精品免费久久| 成年动漫av网址| 交换朋友夫妻互换小说| 亚洲精品美女久久av网站| 国产成人91sexporn| 国产精品嫩草影院av在线观看| 人妻 亚洲 视频| 午夜日本视频在线| 9色porny在线观看| 熟妇人妻不卡中文字幕| 成人毛片60女人毛片免费| 国产精品欧美亚洲77777| 纯流量卡能插随身wifi吗| 永久网站在线| 综合色丁香网| 一级爰片在线观看| av在线观看视频网站免费| 一区二区三区乱码不卡18| 蜜桃在线观看..| 久久久久久伊人网av| 欧美成人午夜精品| 亚洲人成77777在线视频| 丰满乱子伦码专区| 纯流量卡能插随身wifi吗| 母亲3免费完整高清在线观看 | 欧美人与性动交α欧美精品济南到 | 我的女老师完整版在线观看| 丁香六月天网| 大码成人一级视频| 免费黄频网站在线观看国产| 我的女老师完整版在线观看| 天天影视国产精品| 国产亚洲午夜精品一区二区久久| 久久久欧美国产精品| av天堂久久9| 日韩成人伦理影院| 看免费av毛片| 女性生殖器流出的白浆| 精品国产国语对白av| 老司机影院毛片| 最后的刺客免费高清国语| 欧美日韩国产mv在线观看视频| 久久久国产精品麻豆| 亚洲高清免费不卡视频| 亚洲第一区二区三区不卡| 亚洲精品国产av成人精品| 日本免费在线观看一区| 亚洲欧美中文字幕日韩二区| 最新的欧美精品一区二区| 人妻系列 视频| 街头女战士在线观看网站| 婷婷成人精品国产| 亚洲av成人精品一二三区| 内地一区二区视频在线| 日韩大片免费观看网站| 日韩精品免费视频一区二区三区 | 97在线人人人人妻| 少妇的丰满在线观看| 久久国产精品男人的天堂亚洲 | 伦精品一区二区三区| 久久久久久久精品精品| 80岁老熟妇乱子伦牲交| 免费播放大片免费观看视频在线观看| 欧美国产精品va在线观看不卡| 新久久久久国产一级毛片| 精品熟女少妇av免费看| 国产视频首页在线观看| 欧美成人午夜免费资源| 国产福利在线免费观看视频| 久久99热这里只频精品6学生| 老熟女久久久| 丝瓜视频免费看黄片| 亚洲欧美一区二区三区黑人 | 又粗又硬又长又爽又黄的视频| 美女福利国产在线| 日韩在线高清观看一区二区三区| 18禁裸乳无遮挡动漫免费视频| 国产精品蜜桃在线观看| 最新的欧美精品一区二区| 日韩欧美精品免费久久| 91精品三级在线观看| 日韩av免费高清视频| 自线自在国产av| 九草在线视频观看| 国产亚洲最大av| av线在线观看网站| av免费观看日本| 国产精品国产三级国产专区5o| 99国产综合亚洲精品| 日韩成人av中文字幕在线观看| 久久狼人影院| 精品一区二区免费观看| 久热久热在线精品观看| 亚洲一区二区三区欧美精品| 亚洲精品日本国产第一区| 五月开心婷婷网| 国产精品女同一区二区软件| 精品人妻在线不人妻| 成人国语在线视频| 精品99又大又爽又粗少妇毛片| 男女啪啪激烈高潮av片| 欧美人与性动交α欧美软件 | 少妇的逼好多水| 国产日韩欧美亚洲二区| 欧美精品一区二区免费开放| 丝袜脚勾引网站| 亚洲第一区二区三区不卡| 18禁在线无遮挡免费观看视频| 国产探花极品一区二区| 久久久精品94久久精品| 国产免费又黄又爽又色| 99热网站在线观看| 18禁动态无遮挡网站| 亚洲精品自拍成人| www日本在线高清视频| 美女脱内裤让男人舔精品视频| 久久久久人妻精品一区果冻| 国产色婷婷99| 亚洲美女视频黄频| 国产精品久久久久久久久免| 亚洲精品乱码久久久久久按摩| 亚洲精品日本国产第一区| 人体艺术视频欧美日本| 综合色丁香网| 亚洲欧美成人精品一区二区| 亚洲国产精品专区欧美| 国产欧美日韩一区二区三区在线| 久久人人爽人人爽人人片va| 久久久久久久精品精品| 久久精品熟女亚洲av麻豆精品| 搡女人真爽免费视频火全软件| 中文字幕制服av| av免费观看日本| 亚洲精品456在线播放app| av在线观看视频网站免费| 国产精品99久久99久久久不卡 | 亚洲国产最新在线播放| 十八禁网站网址无遮挡| 亚洲欧洲精品一区二区精品久久久 | 午夜视频国产福利| 国产精品人妻久久久影院| 涩涩av久久男人的天堂| 久久久久久久精品精品| 在线观看人妻少妇| 亚洲色图 男人天堂 中文字幕 | 精品一区二区三卡| 国产一区亚洲一区在线观看| 国产日韩欧美亚洲二区| 国产免费福利视频在线观看| 最近手机中文字幕大全| 精品熟女少妇av免费看| 国内精品宾馆在线| 中文字幕人妻丝袜制服| 中文字幕精品免费在线观看视频 | 成人18禁高潮啪啪吃奶动态图| 午夜视频国产福利| 男女国产视频网站| 亚洲丝袜综合中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 涩涩av久久男人的天堂| 少妇人妻 视频| 久久久久久久大尺度免费视频| av播播在线观看一区| 免费看av在线观看网站| 两性夫妻黄色片 | 亚洲久久久国产精品| 美女xxoo啪啪120秒动态图| 午夜精品国产一区二区电影| av在线观看视频网站免费| 最近的中文字幕免费完整| 狠狠婷婷综合久久久久久88av| 热re99久久国产66热| 岛国毛片在线播放| xxxhd国产人妻xxx| 亚洲久久久国产精品| 美女xxoo啪啪120秒动态图| 宅男免费午夜| 制服人妻中文乱码| 尾随美女入室| 亚洲欧洲精品一区二区精品久久久 | 久久久久精品人妻al黑| 欧美精品av麻豆av| 午夜视频国产福利| 大片电影免费在线观看免费| 天天影视国产精品| 一本久久精品| 插逼视频在线观看| 9191精品国产免费久久| 国产男女超爽视频在线观看| av电影中文网址| 韩国av在线不卡| 十八禁网站网址无遮挡| 丰满乱子伦码专区| 夜夜爽夜夜爽视频| 国产一区二区在线观看av| 蜜桃国产av成人99| 亚洲精品自拍成人| 国产老妇伦熟女老妇高清| 亚洲av电影在线进入| 久久 成人 亚洲| 熟女av电影| 夫妻午夜视频| 亚洲第一av免费看| 色吧在线观看| 麻豆精品久久久久久蜜桃| 国产有黄有色有爽视频| 黄色怎么调成土黄色| 乱码一卡2卡4卡精品| 亚洲精品,欧美精品| av卡一久久| 久热这里只有精品99| 精品久久久久久电影网| 国产一级毛片在线| 免费高清在线观看日韩| 王馨瑶露胸无遮挡在线观看| 国产av码专区亚洲av| 国产精品久久久久成人av| 五月玫瑰六月丁香| 人妻 亚洲 视频| 黄片播放在线免费| 大香蕉久久网| 亚洲精品国产av成人精品| 亚洲精品aⅴ在线观看| 国产亚洲一区二区精品| 少妇被粗大猛烈的视频| 亚洲国产成人一精品久久久| 国产成人91sexporn| 夜夜爽夜夜爽视频| 丰满迷人的少妇在线观看| 深夜精品福利| 亚洲高清免费不卡视频| 色网站视频免费| 我要看黄色一级片免费的| 日本-黄色视频高清免费观看| a 毛片基地| 大香蕉久久成人网| 亚洲精品国产av蜜桃| 日韩三级伦理在线观看| 天天影视国产精品| 伦理电影免费视频| 久久久久人妻精品一区果冻| 国产毛片在线视频| 国产日韩一区二区三区精品不卡| 久久久久久久久久成人| 最近2019中文字幕mv第一页| 日本与韩国留学比较| 三上悠亚av全集在线观看| 亚洲av成人精品一二三区| 涩涩av久久男人的天堂| 午夜久久久在线观看| 亚洲欧美成人精品一区二区| 一区二区三区四区激情视频| 大陆偷拍与自拍| 9热在线视频观看99| 日本与韩国留学比较| 在线亚洲精品国产二区图片欧美| 菩萨蛮人人尽说江南好唐韦庄| 久久久久久久亚洲中文字幕| 日韩av免费高清视频| 亚洲精品一二三| 少妇猛男粗大的猛烈进出视频| 日本免费在线观看一区| 咕卡用的链子| 男女下面插进去视频免费观看 | 日韩,欧美,国产一区二区三区| 好男人视频免费观看在线| 国产成人精品久久久久久| 97精品久久久久久久久久精品| 香蕉国产在线看| 男女午夜视频在线观看 | 久久久国产欧美日韩av|