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

    基于溫度植被干旱指數(shù)的華北平原干旱監(jiān)測及時空變化分析

    2023-05-16 02:47:52張兆旭崔津茍文濤肖月
    農業(yè)大數(shù)據(jù)學報 2023年1期
    關鍵詞:華北平原植被河北省

    張兆旭,崔津,茍文濤,肖月

    研究論文

    基于溫度植被干旱指數(shù)的華北平原干旱監(jiān)測及時空變化分析

    張兆旭*,崔津,茍文濤,肖月

    天津工業(yè)大學環(huán)境科學與工程學院,天津 300387

    干旱是地區(qū)水資源供需失衡和水資源短缺的主要原因,干旱的發(fā)生是一個長期、連續(xù)且復雜的過程,是大氣、土壤以及農作物綜合作用的結果,干旱對全國人民的生產和生活都有很大的影響,如何更好地監(jiān)測干旱,對于工農業(yè)、人民的日常生活都有著重要的意義。植被指數(shù)和地表溫度可以描述植被對干旱脅迫的響應,從而反映出土壤水分狀況。文章利用植被指數(shù)和地表溫度的二維特征空間建立了溫度植被干旱指數(shù),將兩個指標的獨特生態(tài)生理學意義結合起來,選取2010—2019年(共10年)的溫度植被干旱指數(shù),對華北平原的干旱進行監(jiān)測與分析,同時將干旱劃分為5個等級(重度濕潤、輕度濕潤、正常、輕度干旱和重度干旱),從年變化、季節(jié)變化、月變化分析華北平原干旱時空變化情況。結果表明,總體上華北平原2010—2019年的干旱以輕度干旱為主,且干旱情況呈現(xiàn)上升的趨勢,干旱指數(shù)值從最低點0.580上升到了最高點0.602,分別出現(xiàn)在2011年與2017年,表明在近10年來,2011年的干旱程度最低,而2017年的干旱程度最嚴重;從季節(jié)角度分析,華北平原的干旱多發(fā)在夏季,夏季的干旱指數(shù)值為0.714;從月度角度分析,華北平原溫度植被干旱指數(shù)最高點與最低點分別出現(xiàn)在7月(0.736)與1月(0.446)。本文對華北平原的干旱研究為該地區(qū)的抗旱防旱提供了可靠的研究支撐。

    干旱監(jiān)測;溫度植被干旱指數(shù);華北平原;遙感

    1 引言

    干旱是由于土壤缺水,植株根部汲取不了充足的水分去補充蒸騰的消耗而引起的自然災害,是人類從古至今面臨的主要災害之一[1]。干旱類型通常包括氣象干旱、農業(yè)干旱、水文干旱以及社會經濟干旱[2-3]。氣象干旱通常主要以降水短缺作為指標,指某一時段內,由于蒸發(fā)量和降水量的收支不平衡,水分支出大于水分收入而造成的水分短缺現(xiàn)象;農業(yè)干旱使土壤的總水分嚴重虧缺,通過研究分析各種植物根系在其生長與發(fā)育變化的整個過程中土壤水分特性變化規(guī)律與形態(tài)特征,反映土壤濕度遠小于正常植物生長最低供水量的程度;當河流流量遠遠小于正常值或含水層水位突然下降時,就會發(fā)生水文干旱;而社會經濟干旱則是指人類工業(yè)生產與消費過程遭受大自然系統(tǒng)和人類社會經濟體系嚴重缺水影響的現(xiàn)象[4-6]。

    干旱的常規(guī)探測方式主要是通過對氣象站點的數(shù)據(jù)進行統(tǒng)計和分析。但是,氣象站的數(shù)據(jù)采集主要依靠氣象站的觀察,而且由于氣象站的位置不可能覆蓋全部地區(qū),因此無法及時、全面地獲得資料,使得對旱情的實時和精確監(jiān)測都會下降[7]。相對于以氣象站為基礎的觀測方式,利用遙感技術進行干旱監(jiān)測可以獲得更大范圍的數(shù)據(jù),更適合地區(qū)和長期的監(jiān)測和評價[8]。近年來,衛(wèi)星遙感技術的客觀、宏觀、真實和高性價比等諸多特征被科學研究工作者所接受,同時其快速發(fā)展也為干旱監(jiān)測開拓了一個嶄新的道路[9]。遙感技術也給干旱監(jiān)測工作帶來了大量的信息,主要觀測手段包括光學、熱紅外以及微波遙感等[10]。

    Sandholt等人綜合考慮地表溫度(Land Surface Temperature,LST)和歸一化植被指數(shù)(Normalized Difference Vegetation Index,NDVI)構建了溫度植被干旱指數(shù)(Temperature Vegetation Dryness Index,TVDI),利用TVDI進行全國旱情監(jiān)測研究能夠較好地反映表層土壤水分變化趨勢,TVDI已被廣泛應用于旱情監(jiān)測中[11]。在TVDI的特征空間中,濕邊代表地表在各植被覆蓋水平上有效的土壤水分最高值,此時水分不是植被生長的限制因素,地表實際蒸散等于潛在蒸散,被認為特征空間的濕邊(wet edge);干邊代表各種植被覆蓋水平上有效土壤水分最低值,地表蒸散也最小,被認為特征空間的干邊(dry edge)。王鵬新等人(2003)提出了條件植被溫度指數(shù)(Vegetation Temperature Condition Index,VTCI)[12],VTCI也被廣泛地應用到干旱監(jiān)測中。干旱是一個緩慢而長期的過程,大多數(shù)農業(yè)干旱指數(shù)是基于線性模型或者單一數(shù)據(jù)源來構建干旱模型(如NDVI、VCI、TCI、SMCI和PCI等),但農業(yè)干旱與不同時間、不同地區(qū)的VI、LST異常不存在線性關系。TVDI整合了農作物生長狀況、地表水熱平衡等信息,發(fā)揮了多源遙感數(shù)據(jù)的優(yōu)勢。目前,TVDI指數(shù)已經被廣泛地應用到了干旱監(jiān)測中[13-15]。

    本文以華北平原為研究區(qū)域,基于2010—2019年該區(qū)域的TVDI進行研究分析,以探索2010—2019華北平原干旱情況的變化過程與趨勢,尋找變化規(guī)律。同時結合華北平原統(tǒng)計數(shù)據(jù),分析華北平原農業(yè)變化情況,同時結合干旱指數(shù),對不同年份、不同月份以及不同季節(jié)的平均分布數(shù)值進行分析,概括總結干旱時空分布變化特征,并且從定性與定量兩個角度分析華北平原時空變化特征。

    2 研究區(qū)概況

    華北平原是中國三大平原之一,是我國人口最集中的區(qū)域之一,橫貫京、津、冀、魯、豫、皖和蘇7個省市(圖1)。華北平原的年平均雨量在500~900 mm之間,年輻射總量為4605~5860 萬億焦(米·年),年平均日照時間為2800 h,主要種植的農作物為小麥和玉米[16]。

    注:基于自然資源部標準地圖服務網站審圖號為GS(2019)1823號的標準地圖制作,底圖邊界無修改,圖4、5、7、9同

    3 研究區(qū)數(shù)據(jù)及方法

    3.1 研究數(shù)據(jù)

    3.1.1 溫度植被干旱指數(shù)數(shù)據(jù)

    TVDI是一種基于光學與熱紅外遙感通道數(shù)據(jù)進行植被覆蓋區(qū)域表層土壤水分反演的方法,TVDI綜合考慮了植被指數(shù)和地表溫度對干旱監(jiān)測的影響,提高了區(qū)域干旱監(jiān)測的準確性。同時由于TVDI的物理意義明確、獲取方便,因此在特定年份的某個時段內,特別是反映整個地區(qū)的干旱程度和狀況,可以用TVDI來研究干旱水平的空間變化特征。

    3.1.2 遙感數(shù)據(jù)預處理

    本文的數(shù)據(jù)來源為國家地球系統(tǒng)科學數(shù)據(jù)中心下載的TVDI數(shù)據(jù)集(http://www.geodata.cn/),數(shù)據(jù)空間投影為WGS-1984地理坐標系,時間為2010—2019年共10年的TVDI數(shù)據(jù)。通過在ArcGIS與ENVI上對數(shù)據(jù)進行處理,得到華北平原10年間各個年份的年分布、月分布以及各個季節(jié)的分布結果。之后再將所需數(shù)據(jù)分別繪制出圖,得到匹配的圖表以便后續(xù)進行分析。

    3.1.3 統(tǒng)計數(shù)據(jù)

    依據(jù)國家統(tǒng)計局網站(http://www.stats.gov.cn/)下載了北京、天津、河北、河南、山東、江蘇和安徽的小麥產量和玉米產量數(shù)據(jù),通過對不同省份和直轄市的產量相加,最終得到華北平原2010—2019年產量統(tǒng)計數(shù)據(jù)。

    4 分析與討論

    4.1 華北平原統(tǒng)計數(shù)據(jù)分析

    華北平原區(qū)土壤豐富,土質肥沃,盛產多種經濟農產品,是中國十分關鍵的糧食產區(qū),最主要的經濟作物是小麥和玉米。而華北平原區(qū)則屬于旱作居多的農耕區(qū)。長江北部區(qū)域以二年三熟的居多,經濟作物一般以小麥、玉米居多,但隨著農業(yè)灌溉事業(yè)發(fā)展,每年二熟的播種面積也在不斷擴大。在長江以南區(qū)域多數(shù)地方均為每年二熟,其中以二年三熟和三年五熟的居多,復種指數(shù)也位居華北地區(qū)前列。

    以小麥為例,通過收集小麥的產量數(shù)據(jù)計算得到華北平原2010—2019年的小麥產量增速圖(圖2)。在圖2中我們可以很明顯地看出小麥產量的增速發(fā)生了多次改變,農業(yè)干旱的發(fā)生是影響農作物生長的主要因素,干旱發(fā)生后如果沒有及時灌溉,最直接的結果是導致農作物產量的降低。從圖2中可以發(fā)現(xiàn),華北平原的小麥產量增速以2017—2018年最為明顯,在這兩年間小麥產量增速大幅下降,甚至在2018年出現(xiàn)了逆增長,增速為負數(shù)的情況,這表明在2017年華北平原有發(fā)生干旱的可能。

    同時以玉米為例,通過計算后同樣得到2010— 2019年間華北平原玉米產量增速圖(圖3)。通過圖3我們可以發(fā)現(xiàn)玉米產量也隨著時間逐年升高,且最高的產量發(fā)生在2015年。玉米產量的增速依舊發(fā)生了多次轉變,在2014年增速突然大幅下降,不過在2015年略有回升。但在2017年之后的幾年增速再一次下降,甚至跌至負數(shù)。

    4.2 華北平原多年平均干旱空間分布

    4.2.1 干旱等級劃分

    依照TVDI干旱等級劃分準則,本文將干旱狀況分為5個等級(見表1),并分別在圖中由綠至紅以不同的顏色表示。

    4.2.2 華北平原2010—2019年多年平均干旱空間分布

    基于每月的TVDI數(shù)據(jù),裁剪得到華北平原每月TVDI變化圖,然后對2010—2019年的全部TVDI數(shù)據(jù)進行區(qū)域平均后得到的研究時間內年平均干旱圖(圖4)。

    圖2 2010—2019年小麥產量增速圖

    圖3 2010—2019年中國玉米產量增速圖

    表1 干旱劃分等級

    從圖4中可以看出,華北平原總體的干旱情況偏向正常與輕度干旱之間。從地域的角度分析,河南、安徽全省以及河北南部、山東南部、江蘇南部地區(qū)均出現(xiàn)了輕度干旱的現(xiàn)象,而北京、天津、河北北部、山東北部以及江蘇北部相對正常。由圖4可以發(fā)現(xiàn):

    圖4 2010—2019年華北平原年均干旱圖(TVDI)

    (1)華北平原2010—2019年的干旱以輕度干旱為主;

    (2)輕度干旱主要發(fā)生在河南省、安徽省大部分區(qū)域以及河北省南部、山東南部和江蘇南部等部分區(qū)域。

    4.3 華北平原2010—2019 年干旱時空變化分析

    4.3.1 年變化分析

    在對2010—2019年每年所有的TVDI進行年平均后,可以得到各年間的平均干旱數(shù)值,將每年的數(shù)據(jù)綜合統(tǒng)計后,得到了圖5(a)~(j)每年的干旱圖。之后將各個年份的區(qū)域平均值提取計算得到每年的干旱局域平均值(圖6)。

    空間上,從圖5中可以發(fā)現(xiàn),2010年華北平原的干旱主要發(fā)生在河南省、安徽省的大部分區(qū)域以及河北的南部區(qū)域;2011年華北平原的干旱則出現(xiàn)在安徽省大部分區(qū)域與河南省的南部區(qū)域;2012年華北平原河南省的干旱程度逐漸上升,河北省南部地區(qū)出現(xiàn)輕度干旱,但安徽省的干旱情況得到一定程度的緩解;2013年華北平原的干旱高度集中在河南省與安徽省,基本上涵蓋了省內全部區(qū)域;2014年華北平原干旱情況愈加嚴重,河南省、安徽省、山東省大部分區(qū)域與河北省、江蘇省南部均出現(xiàn)輕度干旱;2015年華北平原干旱情況與前一年相比變化不大,僅僅河北省與江蘇省的情況略有好轉;2016年華北平原的干旱得到抑制,只有河南省大部分區(qū)域與安徽省、山東省的小部分區(qū)域出現(xiàn)干旱情況;2017年華北平原的干旱達到了近幾年的最高值,全部平原只有河北省北部區(qū)域保持正常水平;2018年華北平原干旱危害輕微下降,山東省、江蘇省部分區(qū)域也恢復正常;2019年華北平原河北省、山東省與江蘇省的北部區(qū)域均恢復正常,只剩下河南省、安徽省大部分區(qū)域與河北省、山東省和江蘇省南部區(qū)域仍出現(xiàn)輕度干旱情況。

    時間上,從圖6中我們可以直觀地發(fā)現(xiàn)在10年間華北平原的干旱情況整體呈現(xiàn)出的是一個上升的趨勢。TVDI越大,則代表干旱情況越嚴重,TVDI數(shù)值從最低點0.580上升到了最高點0.602,分別出現(xiàn)在2011年與2017年,代表了在研究時間段內,2011年的干旱程度最低,而2017年的干旱程度最高,整個華北平原將近3/4的區(qū)域都處于輕度干旱范圍。這10年間的TVDI數(shù)值在大趨勢是增加的基礎上經歷了4次轉折:首先2011—2014年干旱程度逐年上升,而2014年之后干旱情況又有所好轉,這說明2014年處在一個轉折點的位置。下降持續(xù)到2016年結束,在2017年干旱程度達到了頂峰,由此之后干旱變得越來越嚴重,雖然相較于頂峰有所下降,但每年基本上都會有超過2/3的區(qū)域仍處在輕度干旱范圍。

    圖6 2010—2019年TVDI散點圖

    4.3.2 季節(jié)變化分析

    在對2010—2019年所有的TVDI進行季節(jié)平均后(春季為3—5月,夏季為6—8月,秋季為9—11月,冬季為12—2月),可以得到各季節(jié)華北平原的季節(jié)干旱空間分布圖(圖7),然后將每個季節(jié)的空間分布圖進行區(qū)域平均,得到了各季節(jié)TVDI區(qū)域平均數(shù)值的柱狀圖(圖8)。

    對圖7進行空間分析,可以發(fā)現(xiàn),華北平原的春季干旱以正常區(qū)域與輕度干旱為主,其中輕度干旱主要分布在河北省南部、安徽省南部、山東省東部與河南省西部區(qū)域,而其他區(qū)域則處于正常狀況;夏季的華北平原屬于干旱頻發(fā)時期,幾乎全部平原均處于輕度干旱;華北平原的秋季主要以正常為主,輕度干旱主要集中在河南省與安徽省的大部分區(qū)域;華北平原在冬季會出現(xiàn)輕度濕潤的情況,尤其濕潤主要集中在河北省北部區(qū)域。對圖8進行定量分析,可以發(fā)現(xiàn)華北平原的夏季是TVDI數(shù)值最高的,而冬季則最低。春、秋兩季情況比較來說相差較小。

    (a)春季,(b)夏季,(c)秋季,(d)冬季

    圖8 2010—2019年各季節(jié)TVDI柱狀圖

    正如一般情況,華北平原的夏季干旱現(xiàn)象十分明顯,基本上華北平原內的所有地區(qū)均處于輕度干旱范圍;而冬季正好與之相反,平原內無明顯干旱地區(qū),甚至部分北方地區(qū)出現(xiàn)輕度濕潤情況。同時需要注意的是,華北平原的春季也較明顯地發(fā)生了部分干旱現(xiàn)象,原因則是華北平原的春季素來有“春旱”的說法,這是指華北平原的春季氣溫回升較快,蒸發(fā)較強,而夏季風弱,雨季未到,降水較少,導致了春季的土壤缺水較為嚴重,極易引發(fā)干旱。

    4.3.3 月變化分析

    在對2010—2019年所有的溫度植被干旱指數(shù)進行月平均后,可以得到各月華北平原的月份干旱空間分布圖(圖9),然后將各個月的數(shù)據(jù)綜合統(tǒng)計,計算出平均數(shù)值后,得到了各月份TVDI的柱狀圖(圖10)。

    對華北平原1—12月(圖9)干旱圖進行空間分析可以發(fā)現(xiàn),1月華北平原主要包括正常情況與輕度濕潤,其中輕度濕潤情況主要分布在河北省與北京市的北部區(qū)域,其他地區(qū)則為正常;2月華北平原的濕潤地區(qū)逐漸減少;3月華北平原輕度濕潤基本上已經消失不見,且在河北省北部與河南省西部均出現(xiàn)零星的輕度干旱情況;4月華北平原干旱情況愈加嚴重,河北省、山東省與河南省部分區(qū)域均出現(xiàn)輕度干旱情況;隨著氣溫的上升,地表溫度也隨之上升,再加上雨季未到,降水較少,5—8月的華北平原全部區(qū)域持續(xù)出現(xiàn)輕度干旱現(xiàn)象;9月份后,降水開始增加,9月的華北平原雖然大部分地區(qū)仍處于輕度干旱情況,但河北省的北部區(qū)域已經恢復正常;10月的華北平原氣溫逐漸下降,輕度干旱情況變?yōu)榧性诤幽鲜 不帐〈蟛糠值貐^(qū)與河北省、山東省以及江蘇省的南部區(qū)域,其余地區(qū)已恢復正常;11月華北平原臨近冬季,全部區(qū)域均恢復正常,河北省北部區(qū)域甚至出現(xiàn)輕度濕潤現(xiàn)象;12月華北平原仍以正常情況為主,但河北省北部大部分區(qū)域全部變?yōu)檩p度濕潤,干旱情況全部消失不見。

    (a)一月,(b)二月,(c)三月,(d)四月,(e)五月,(f)六月 (g)七月,(h)八月,(i)九月,(j)十月,(k)十一月,(l)十二月

    圖9 2010—2019年各月份華北平原平均干旱圖

    Fig. 9 Average drought map of North China Plain from 2010 to 2019

    對圖10進行定量分析,可以發(fā)現(xiàn),TVDI的最高點與最低點分別出現(xiàn)在7月與1月,其中5—8月的干旱程度較為嚴重,基本上干旱已經覆蓋了整片華北平原,均屬于輕度干旱的范圍;而12月與1月則干旱程度較輕,在北方的部分區(qū)域均一定程度上出現(xiàn)了輕度濕潤;其他月份則相對來說比較正常,干旱情況處在于可控范圍內。

    圖10 2010—2019年各月份TVDI柱狀圖

    5 結論

    華北平原是全國人口最多的平原之一,而其土地質量的情況也由此十分重要,更是關系到糧食等作物的產量與經濟情況。本文以華北平原為研究區(qū),以遙感數(shù)據(jù),地表溫度-植被指數(shù)特征空間原理為基礎,構建地表溫度-植被指數(shù)特征空間,計算TVDI指數(shù),獲取華北平原年、月、季節(jié)干旱分布圖,分析華北平原干旱情況。主要結論如下。

    (1)將計算得到的TVDI分為重度濕潤(0~0.2)、輕度濕潤(0.2~0.4)、正常(0.4~0.6)、輕度干旱(0.6~0.8)和重度干旱(0.8~1)5個等級,根據(jù)干旱等級劃分整個華北平原,并依此繪制華北平原平均干旱分布圖。華北平原2010—2019年的干旱程度以輕度干旱為主,輕度干旱主要發(fā)生在河南省、安徽省大部分區(qū)域以及河北省的南部,山東的南部、江蘇南部等部分區(qū)域。

    (2)華北平原從年變化角度分析主要輕度干旱年份在2017—2019年。整體干旱情況呈現(xiàn)上升的趨勢。TVDI最高與最低分別出現(xiàn)在2011年與2017年,且在這之間的變化又經歷了多次轉折。

    (3)從季節(jié)角度分析,華北平原的干旱多發(fā)在夏季,同時春季也有干旱現(xiàn)象的發(fā)生,華北平原一直以來的“春旱”情況,研究結果表明春季的干旱情況較為嚴重。

    (4)通過月干旱分析,華北平原一年間TVDI最高點與最低點分別出現(xiàn)在7月與1月,其中5—8月的干旱情況因為溫度上升導致最為嚴重,而12月至1月的土地在氣溫回落之后則以濕潤為主,其他月份的干旱均較為正常,僅僅部分出現(xiàn)干旱情況。

    干旱是一個包括氣象、水文、農業(yè)等綜合作用的過程,隨著遙感技術的快速發(fā)展,不同遙感技術提供了不同的手段來監(jiān)測農業(yè)干旱。本文基于TVDI指數(shù)監(jiān)測了華北平原的干旱情況,為抗旱對策的制定提供依據(jù)。通過對旱情的監(jiān)測與分析,為農業(yè)生產的結構調整、農業(yè)技術的實施、農業(yè)人員的工作、農業(yè)生產的可持續(xù)發(fā)展提供了堅實可信的依據(jù)。

    [1] 張兆旭, 秦其明. 基于多源遙感數(shù)據(jù)的黃淮海平原農業(yè)干旱監(jiān)測研究[J]. 農業(yè)與技術, 2021, 41(2): 84-92.

    Zhang Z X, Qin Q M. Study on agricultural drought monitoring in Huang-Huai-Hai plain based on multi-source remote sensing data[J]. Agriculture and Technology, 2021, 41(2): 84-92.

    [2] Mladenova I E, Bolten J D, Crow W T, et al. Evaluating the operational application of SMAP for global agricultural drought monitoring[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2019, 12(9): 3387-3397.

    [3] Wilhite D A, Glantz M H. Understanding: the drought phenomenon: the role of definitions[J]. Water International, 1985, 10(3): 111- 120.

    [4] Mishra A K, Desai V R. Drought forecasting using stochastic models[J]. Stochastic Environmental Research and Risk Assessment, 2005, 19(5): 326-339.

    [5] Mishra A K Singh V P. A review of drought concepts[J]. Journal of hydrology, 2010, 391(1-2): 202-216.

    [6] Wilhite D A, Svoboda M D, Hayes M J. Understanding the complex impacts of drought: A key to enhancing drought mitigation and preparedness[J]. Water Resources Management, 2007, 21(5): 763-774.

    [7] Maes W, Steppe K. Estimating evapotranspiration and drought stress with ground-based thermal remote sensing in agriculture: a review[J]. Journal of Experimental Botany, 2012, 63(13): 4671-4712.

    [8] KanellouE, DomenikiotisC, TsirosE, et al. Satellite- based drought estimation in Thessaly[J]. European Water, 2008, 23(24): 111-122.

    [9] QinQ, WuZ, ZhangT, et al. Optical and thermal remote sensing for monitoring agricultural drought[J]. Remote Sensing, 2021, 13(24): 1-34.

    [10] 宋煒. 遙感監(jiān)測干旱的方法綜述[J]. 中國西部科技, 2011, 10(13): 42-44.

    Song W. Summary of drought monitoring methods by remote sensing[J]. Science and Technology of West China, 2011, 10(13): 42-44.

    [11] SandholtI, RasmussenK, AndersenJ. A simple interpretation of the surface temperature/vegetation index space for assessment of surface moisture status[J]. Remote Sensing of Environment, 2002, 79(2-3): 213-224.

    [12] 王鵬新, 龔健雅, 李小文, 等. 基于植被指數(shù)和土地表面溫度的干旱監(jiān)測模型[J]. 地球科學進展, 2003, 18(4):40-46.

    Wang P X, Gong J Y, Li X W, et al. Drought monitoring model based on vegetation index and land surface temperature[J]. Advances in Earth Science, 2003, 18(4):40-46.

    [13] BaiJ, YuanY, DiL. Comparison between TVDI and CWSI for drought monitoring in the Guanzhong Plain, China[J]. Journal of Integrative Agriculture, 2017, 16(2): 389-397.

    [14] Liang L, Zhao S, Qin Z, et al. Drought change trend using MODIS TVDI and its relationship with climate factors in China from 2001 to 2010[J]. Journal of Integrative Agriculture, 2014, 13(7): 1501- 1508.

    [15] Shashikant V, Mohamed Shariff A R, Wayayok A, et al. Utilizing TVDI and NDWI to classify severity of agricultural drought in Chuping, Malaysia[J]. Agronomy, 2021, 11(6): 1-18.

    [16] Zhang Z, Xu W, Shi Z, et al. Establishment of a comprehensive drought monitoring index based on multisource remote sensing data and agricultural drought monitoring[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2021, 14: 2113-2126.

    Drought Monitoring and Spatiotemporal Changes Analysis in North China Plain Based on Temperature Vegetation Dryness Index

    ZHANG Zhaoxu*, CUI Jin, GOU Wentao, XIAO Yue

    School of Environmental Science and Engineering, Tiangong University, Tianjin 300387, China

    Drought is the main reason for the imbalance between supply and demand of water resources and the shortage of water resources in the region. The occurrence of drought is a long-term, continuous and complex process, which is the result of the comprehensive action of atmosphere, soil and crops. Drought has a great impact on the production and life of the people all over the country. How to better monitor drought is of great significance to industry, agriculture and people's daily life. Vegetation index and temperature can describe the response of vegetation to drought stress, thus reflecting the soil water status. In this paper, the temperature vegetation dryness index is established by using the two-dimensional characteristic space of vegetation index and surface temperature. The unique ecological and physiological significance of the two indexes is combined. The temperature vegetation dryness index from 2010 to 2019 (10 years) is selected to monitor and analyze the drought in the North China Plain. Drought is divided into five grades (heavy humidity, light humidity, normal, light drought, heavy drought), from the annual change seasonal and monthly changes analyze the temporal and spatial changes of drought. The results show that the drought is mainly mild drought, and the drought situation shows an upward trend. From the perspective of season, drought often occurs in summer. From the month, the highest and lowest temperature vegetation dryness index appeared in July and January respectively. The study of drought in North China Plain provides a reliable research support for drought resistance and drought prevention in this region.

    drought monitoring; temperature vegetation dryness index; north China plain; remote sensing

    張兆旭,崔津,茍文濤,等. 基于溫度植被干旱指數(shù)的華北平原干旱監(jiān)測及時空變化分析[J]. 農業(yè)大數(shù)據(jù)學報, 2023,5(1):95-107.

    ZHANG Zhaoxu, CUI Jin, GOU Wentao, et al. Drought monitoring and spatiotemporal changes analysis in North China Plain based on temperature vegetation dryness index[J]. Journal of Agricultural Big Data, 2023, 5(1):95-107.

    10.19788/j.issn.2096-6369.230118

    2023-

    海水資源利用技術發(fā)展研究與報告編制(22-02-01018A-0023), 北京未名福科技有限公司開放基金(22-02-01018A-0017)

    第一作者張兆旭兼通信作者,男,博士研究生,研究方向:定量遙感與干旱監(jiān)測;Email:zhangzhaoxu@tiangong. edu. cn。

    猜你喜歡
    華北平原植被河北省
    基于植被復綠技術的孔植試驗及應用
    河北地質(2022年2期)2022-08-22 06:24:04
    2020年河北省對外經濟
    國家開發(fā)銀行河北省分行
    河北省
    河北省
    綠色植被在溯溪旅游中的應用
    清晨
    詩潮(2017年2期)2017-03-16 11:04:01
    基于原生植被的長山群島植被退化分析
    基于NDVI的鶴壁市植被覆蓋動態(tài)分析研究
    河南科技(2014年4期)2014-02-27 14:07:25
    華北平原淺層地下水污染嚴重
    城市地質(2013年1期)2013-03-11 15:20:07
    亚洲欧美精品综合一区二区三区| 日日夜夜操网爽| 在线观看美女被高潮喷水网站 | 国产精品亚洲美女久久久| 蜜桃久久精品国产亚洲av| 熟女少妇亚洲综合色aaa.| 亚洲人成伊人成综合网2020| 精品一区二区三区四区五区乱码| 99久久国产精品久久久| 成人国语在线视频| 国产日本99.免费观看| 成人三级做爰电影| 给我免费播放毛片高清在线观看| 19禁男女啪啪无遮挡网站| 欧美黄色片欧美黄色片| 90打野战视频偷拍视频| 最近视频中文字幕2019在线8| 少妇粗大呻吟视频| svipshipincom国产片| 欧美乱妇无乱码| 露出奶头的视频| 国产三级黄色录像| 人人妻人人澡欧美一区二区| 在线观看舔阴道视频| 午夜日韩欧美国产| 亚洲成人免费电影在线观看| 欧美黑人欧美精品刺激| 在线观看美女被高潮喷水网站 | 99国产精品一区二区三区| 久久久久久大精品| 18美女黄网站色大片免费观看| 亚洲精品一区av在线观看| 亚洲成av人片在线播放无| 久久久久久久久免费视频了| 久久九九热精品免费| 日韩欧美 国产精品| 精品一区二区三区av网在线观看| 老司机午夜福利在线观看视频| 香蕉丝袜av| 国产精品久久久久久人妻精品电影| 香蕉丝袜av| 18禁裸乳无遮挡免费网站照片| 欧美zozozo另类| 我要搜黄色片| 十八禁人妻一区二区| 免费看十八禁软件| 麻豆成人午夜福利视频| 免费看十八禁软件| 亚洲国产精品成人综合色| 免费人成视频x8x8入口观看| 深夜精品福利| 欧美又色又爽又黄视频| 亚洲人成网站在线播放欧美日韩| 久久久久免费精品人妻一区二区| 小说图片视频综合网站| 亚洲精品国产精品久久久不卡| 久久久精品欧美日韩精品| 午夜老司机福利片| 窝窝影院91人妻| 两个人的视频大全免费| 在线免费观看的www视频| 18禁裸乳无遮挡免费网站照片| 日本一二三区视频观看| 18禁黄网站禁片午夜丰满| 啦啦啦免费观看视频1| 18禁美女被吸乳视频| 美女午夜性视频免费| 麻豆久久精品国产亚洲av| 曰老女人黄片| 精品国产乱子伦一区二区三区| 精品免费久久久久久久清纯| av国产免费在线观看| 亚洲专区字幕在线| 一级作爱视频免费观看| x7x7x7水蜜桃| 最近最新中文字幕大全电影3| 成人特级黄色片久久久久久久| 国产爱豆传媒在线观看 | 色精品久久人妻99蜜桃| 在线播放国产精品三级| 精品国产美女av久久久久小说| 午夜a级毛片| 男女床上黄色一级片免费看| 国模一区二区三区四区视频 | 国产在线精品亚洲第一网站| 激情在线观看视频在线高清| 啦啦啦免费观看视频1| 狠狠狠狠99中文字幕| 狂野欧美白嫩少妇大欣赏| 亚洲av中文字字幕乱码综合| 免费在线观看影片大全网站| 变态另类丝袜制服| 午夜福利视频1000在线观看| 亚洲人成网站高清观看| 老司机在亚洲福利影院| 国产精品免费一区二区三区在线| 欧美在线一区亚洲| 欧美国产日韩亚洲一区| 亚洲av片天天在线观看| 女人高潮潮喷娇喘18禁视频| 欧美一区二区国产精品久久精品 | 亚洲欧美日韩高清在线视频| 午夜福利在线在线| 1024香蕉在线观看| 无人区码免费观看不卡| 国产精品永久免费网站| 99久久久亚洲精品蜜臀av| 五月伊人婷婷丁香| 大型av网站在线播放| 又紧又爽又黄一区二区| 午夜福利在线在线| 免费人成视频x8x8入口观看| 亚洲精品国产一区二区精华液| 曰老女人黄片| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人精品无人区| 国产成人aa在线观看| 久久精品国产综合久久久| 国产高清视频在线播放一区| 51午夜福利影视在线观看| 可以在线观看的亚洲视频| 午夜福利成人在线免费观看| 天天一区二区日本电影三级| 天天躁夜夜躁狠狠躁躁| 久久精品影院6| 午夜福利视频1000在线观看| 两个人看的免费小视频| 人人妻,人人澡人人爽秒播| 国产精品,欧美在线| 最近在线观看免费完整版| 国产久久久一区二区三区| 国产伦在线观看视频一区| 国产在线观看jvid| 免费看美女性在线毛片视频| 日韩欧美国产在线观看| 夜夜看夜夜爽夜夜摸| 欧美av亚洲av综合av国产av| 久久性视频一级片| 法律面前人人平等表现在哪些方面| 日韩欧美国产在线观看| 国产亚洲av高清不卡| av在线播放免费不卡| 亚洲性夜色夜夜综合| 老司机午夜福利在线观看视频| 真人做人爱边吃奶动态| 不卡一级毛片| 一夜夜www| 久久这里只有精品19| www日本在线高清视频| 国产亚洲精品第一综合不卡| 亚洲黑人精品在线| 久久久久精品国产欧美久久久| 精品少妇一区二区三区视频日本电影| 久久中文看片网| 国产欧美日韩精品亚洲av| 97超级碰碰碰精品色视频在线观看| 99精品在免费线老司机午夜| 村上凉子中文字幕在线| 俄罗斯特黄特色一大片| x7x7x7水蜜桃| 757午夜福利合集在线观看| ponron亚洲| 国内毛片毛片毛片毛片毛片| 国产精品亚洲av一区麻豆| 午夜精品一区二区三区免费看| 一二三四在线观看免费中文在| 99久久久亚洲精品蜜臀av| 亚洲精品美女久久久久99蜜臀| 两性夫妻黄色片| 精品久久久久久久久久久久久| 12—13女人毛片做爰片一| 久久中文看片网| 操出白浆在线播放| 久久精品国产综合久久久| 每晚都被弄得嗷嗷叫到高潮| 91av网站免费观看| 亚洲av电影不卡..在线观看| 此物有八面人人有两片| 欧美一区二区国产精品久久精品 | 精品午夜福利视频在线观看一区| 欧美不卡视频在线免费观看 | 制服人妻中文乱码| 欧美三级亚洲精品| 国产精品av久久久久免费| 亚洲人成伊人成综合网2020| 真人做人爱边吃奶动态| 国产精品,欧美在线| 男女视频在线观看网站免费 | 不卡一级毛片| 国产伦人伦偷精品视频| 少妇的丰满在线观看| 亚洲中文字幕日韩| 亚洲av美国av| 久久中文字幕一级| 亚洲av成人精品一区久久| 成年版毛片免费区| 亚洲九九香蕉| 欧美丝袜亚洲另类 | 国产精品 欧美亚洲| 亚洲午夜精品一区,二区,三区| 波多野结衣高清作品| 日本一二三区视频观看| 少妇熟女aⅴ在线视频| 国产精品久久久av美女十八| 国产黄色小视频在线观看| 欧美久久黑人一区二区| 少妇粗大呻吟视频| 18禁裸乳无遮挡免费网站照片| 欧美大码av| 亚洲精品在线观看二区| 国产精品久久视频播放| 91老司机精品| 日韩成人在线观看一区二区三区| 国产亚洲精品久久久久5区| 午夜激情福利司机影院| 熟女电影av网| 国产精品一及| 亚洲欧美激情综合另类| 激情在线观看视频在线高清| 国产精品av久久久久免费| 首页视频小说图片口味搜索| 久热爱精品视频在线9| 久久精品国产清高在天天线| 黄频高清免费视频| 一区二区三区国产精品乱码| 精品久久久久久,| 男人的好看免费观看在线视频 | 波多野结衣高清作品| 在线观看日韩欧美| 欧美不卡视频在线免费观看 | 黄色a级毛片大全视频| 亚洲电影在线观看av| 免费看美女性在线毛片视频| √禁漫天堂资源中文www| 欧美 亚洲 国产 日韩一| 亚洲熟妇熟女久久| 99久久精品国产亚洲精品| 欧美 亚洲 国产 日韩一| 日本 av在线| 这个男人来自地球电影免费观看| 超碰成人久久| 日本精品一区二区三区蜜桃| 精品日产1卡2卡| 国产爱豆传媒在线观看 | 好看av亚洲va欧美ⅴa在| 国产精品久久久人人做人人爽| 久久久久国内视频| 精品欧美一区二区三区在线| 国产亚洲精品第一综合不卡| 91成年电影在线观看| 久久久久精品国产欧美久久久| 久久精品91蜜桃| 搡老妇女老女人老熟妇| 国产人伦9x9x在线观看| 欧美在线黄色| 欧美成狂野欧美在线观看| 天天一区二区日本电影三级| 岛国视频午夜一区免费看| 亚洲成a人片在线一区二区| 色老头精品视频在线观看| 国产高清激情床上av| 非洲黑人性xxxx精品又粗又长| 在线观看舔阴道视频| 亚洲精品美女久久久久99蜜臀| 狂野欧美白嫩少妇大欣赏| 日本成人三级电影网站| videosex国产| 91老司机精品| 色在线成人网| 女警被强在线播放| 亚洲一区中文字幕在线| 听说在线观看完整版免费高清| 成人亚洲精品av一区二区| 一本大道久久a久久精品| 美女黄网站色视频| 国产麻豆成人av免费视频| 麻豆成人av在线观看| 久久精品影院6| 高清在线国产一区| 男人舔奶头视频| 制服诱惑二区| 国产熟女xx| 亚洲熟妇中文字幕五十中出| 97超级碰碰碰精品色视频在线观看| 亚洲电影在线观看av| 一个人免费在线观看的高清视频| 国产高清videossex| 两个人的视频大全免费| 美女黄网站色视频| 国产aⅴ精品一区二区三区波| 国产精品久久久久久久电影 | 国产午夜精品论理片| 欧美av亚洲av综合av国产av| 男人的好看免费观看在线视频 | 男女午夜视频在线观看| 久久久久九九精品影院| 又黄又爽又免费观看的视频| 国产成年人精品一区二区| 久久伊人香网站| 三级国产精品欧美在线观看 | 少妇裸体淫交视频免费看高清 | 97碰自拍视频| 亚洲国产精品成人综合色| 精品第一国产精品| 村上凉子中文字幕在线| 免费在线观看亚洲国产| 亚洲电影在线观看av| 最近最新免费中文字幕在线| 99热这里只有是精品50| 女人高潮潮喷娇喘18禁视频| 狂野欧美白嫩少妇大欣赏| 99国产精品一区二区三区| 国产精品自产拍在线观看55亚洲| 午夜久久久久精精品| 午夜a级毛片| av福利片在线| 亚洲精品中文字幕在线视频| 51午夜福利影视在线观看| 国产熟女xx| 久久精品aⅴ一区二区三区四区| 成年女人毛片免费观看观看9| 亚洲欧美精品综合一区二区三区| 久久精品国产亚洲av高清一级| 久久国产精品人妻蜜桃| 午夜福利高清视频| 可以在线观看的亚洲视频| 亚洲国产日韩欧美精品在线观看 | 久久久久久九九精品二区国产 | 婷婷精品国产亚洲av| 免费高清视频大片| 可以在线观看毛片的网站| 国产亚洲精品第一综合不卡| 又爽又黄无遮挡网站| 日本成人三级电影网站| 桃红色精品国产亚洲av| 一本综合久久免费| 制服丝袜大香蕉在线| 久久久久精品国产欧美久久久| 亚洲片人在线观看| 久久久国产成人免费| 亚洲一区二区三区色噜噜| 亚洲九九香蕉| 欧美乱码精品一区二区三区| 18禁国产床啪视频网站| 午夜久久久久精精品| 波多野结衣高清无吗| 亚洲人与动物交配视频| 黑人巨大精品欧美一区二区mp4| 亚洲aⅴ乱码一区二区在线播放 | 国内毛片毛片毛片毛片毛片| 国产精品一区二区免费欧美| 免费电影在线观看免费观看| 欧美3d第一页| 97人妻精品一区二区三区麻豆| 久久性视频一级片| 黄色 视频免费看| 午夜亚洲福利在线播放| 亚洲欧美日韩东京热| 日韩欧美三级三区| 国产成人av激情在线播放| 日韩 欧美 亚洲 中文字幕| 欧美成人免费av一区二区三区| 我的老师免费观看完整版| 国产精华一区二区三区| 99国产精品一区二区三区| 一个人免费在线观看电影 | 国产单亲对白刺激| 色在线成人网| 男人的好看免费观看在线视频 | 欧美日韩一级在线毛片| 曰老女人黄片| 亚洲aⅴ乱码一区二区在线播放 | 制服丝袜大香蕉在线| 久久国产精品人妻蜜桃| 国产高清视频在线观看网站| 欧美色欧美亚洲另类二区| 99国产极品粉嫩在线观看| 18禁裸乳无遮挡免费网站照片| 在线观看免费视频日本深夜| 制服丝袜大香蕉在线| 欧美 亚洲 国产 日韩一| 一二三四社区在线视频社区8| 精品国产乱子伦一区二区三区| 国产成人系列免费观看| 色av中文字幕| 一区二区三区高清视频在线| 国产成人aa在线观看| 婷婷精品国产亚洲av| 麻豆av在线久日| 国产激情欧美一区二区| 床上黄色一级片| 老司机午夜福利在线观看视频| 巨乳人妻的诱惑在线观看| 久久久国产成人免费| 精品电影一区二区在线| 在线观看午夜福利视频| 18禁黄网站禁片午夜丰满| 不卡一级毛片| 一进一出抽搐动态| 男人舔奶头视频| 国产视频内射| www日本黄色视频网| 成人av在线播放网站| 亚洲专区国产一区二区| 制服诱惑二区| 在线a可以看的网站| 久久久久国内视频| 又黄又爽又免费观看的视频| 亚洲五月天丁香| 亚洲 欧美一区二区三区| 久久香蕉国产精品| 正在播放国产对白刺激| 天堂√8在线中文| 欧美日韩亚洲国产一区二区在线观看| 久久这里只有精品中国| 欧美乱妇无乱码| 不卡一级毛片| 日本 欧美在线| 亚洲专区字幕在线| 亚洲国产精品成人综合色| 亚洲av成人一区二区三| 欧美日韩亚洲综合一区二区三区_| 高潮久久久久久久久久久不卡| 亚洲国产欧美人成| 久久久精品大字幕| 床上黄色一级片| 级片在线观看| 国产午夜精品久久久久久| 国产人伦9x9x在线观看| 黄色视频不卡| 午夜福利成人在线免费观看| 免费电影在线观看免费观看| 蜜桃久久精品国产亚洲av| 精华霜和精华液先用哪个| 亚洲成av人片免费观看| 欧美成人午夜精品| tocl精华| 熟女电影av网| 久99久视频精品免费| 久热爱精品视频在线9| 久久九九热精品免费| 国产精品1区2区在线观看.| 午夜福利在线观看吧| 国产亚洲精品av在线| 老汉色av国产亚洲站长工具| 淫妇啪啪啪对白视频| 国内揄拍国产精品人妻在线| 国产在线精品亚洲第一网站| 亚洲精品久久成人aⅴ小说| 精品第一国产精品| 久久亚洲精品不卡| 99热这里只有是精品50| 90打野战视频偷拍视频| 久久久久国内视频| 欧美激情久久久久久爽电影| 国产黄a三级三级三级人| 国产亚洲精品综合一区在线观看 | 日本撒尿小便嘘嘘汇集6| 我的老师免费观看完整版| 很黄的视频免费| 国产高清视频在线播放一区| 久久精品国产99精品国产亚洲性色| ponron亚洲| 黄片大片在线免费观看| 亚洲专区国产一区二区| 黄频高清免费视频| 欧美日韩国产亚洲二区| 国产一区二区三区在线臀色熟女| 麻豆一二三区av精品| 国产免费av片在线观看野外av| 校园春色视频在线观看| 1024香蕉在线观看| 国产久久久一区二区三区| 亚洲国产欧洲综合997久久,| 中文字幕最新亚洲高清| 国产精品九九99| 高潮久久久久久久久久久不卡| 久久久精品欧美日韩精品| 老司机福利观看| 亚洲av美国av| 国产精华一区二区三区| 国产精品影院久久| 国产精品久久电影中文字幕| 午夜福利欧美成人| 最近最新中文字幕大全电影3| 国产av一区二区精品久久| 国产成人av激情在线播放| 国产精品免费一区二区三区在线| 日韩欧美在线乱码| 日韩三级视频一区二区三区| 老司机深夜福利视频在线观看| 国产三级黄色录像| 亚洲av中文字字幕乱码综合| 国产高清有码在线观看视频 | 五月玫瑰六月丁香| 国产av一区在线观看免费| 久久久水蜜桃国产精品网| 国产在线精品亚洲第一网站| 俄罗斯特黄特色一大片| 国产高清videossex| 久久久久久久午夜电影| 成人高潮视频无遮挡免费网站| cao死你这个sao货| 久久久久久人人人人人| 欧美日韩福利视频一区二区| 国产在线观看jvid| 久久久久久亚洲精品国产蜜桃av| 18禁国产床啪视频网站| 精品高清国产在线一区| 国产精品自产拍在线观看55亚洲| 麻豆成人av在线观看| 18美女黄网站色大片免费观看| 一级黄色大片毛片| 日日干狠狠操夜夜爽| 国内精品一区二区在线观看| 哪里可以看免费的av片| 人妻丰满熟妇av一区二区三区| 一区福利在线观看| 在线观看免费视频日本深夜| 国产一区二区激情短视频| 午夜激情福利司机影院| 亚洲欧美精品综合一区二区三区| 一个人免费在线观看电影 | 最好的美女福利视频网| 欧美色视频一区免费| 欧美性长视频在线观看| 午夜亚洲福利在线播放| 一级毛片高清免费大全| 亚洲狠狠婷婷综合久久图片| 精品午夜福利视频在线观看一区| 99久久久亚洲精品蜜臀av| 国产成人一区二区三区免费视频网站| 亚洲va日本ⅴa欧美va伊人久久| 舔av片在线| 黄色a级毛片大全视频| 亚洲18禁久久av| www.熟女人妻精品国产| 国产精品久久久av美女十八| 欧美日韩一级在线毛片| 成年版毛片免费区| 这个男人来自地球电影免费观看| 国产精品电影一区二区三区| 波多野结衣高清作品| 久久久久性生活片| 白带黄色成豆腐渣| 黄频高清免费视频| 国产久久久一区二区三区| 好男人电影高清在线观看| 精品久久久久久久末码| 脱女人内裤的视频| 女人高潮潮喷娇喘18禁视频| 久久天躁狠狠躁夜夜2o2o| 久久国产精品影院| 99久久无色码亚洲精品果冻| 伊人久久大香线蕉亚洲五| 国产精华一区二区三区| 美女免费视频网站| 精品久久久久久成人av| 琪琪午夜伦伦电影理论片6080| 人人妻人人看人人澡| av有码第一页| 一本综合久久免费| 欧美性猛交╳xxx乱大交人| 亚洲 欧美 日韩 在线 免费| 青草久久国产| 欧美av亚洲av综合av国产av| 精品电影一区二区在线| 1024手机看黄色片| 亚洲片人在线观看| 一区二区三区国产精品乱码| 欧美日韩乱码在线| 一区二区三区高清视频在线| 欧美日韩福利视频一区二区| 国产蜜桃级精品一区二区三区| 久久 成人 亚洲| 久久中文字幕人妻熟女| 久久婷婷成人综合色麻豆| 日本一区二区免费在线视频| 国产真人三级小视频在线观看| 人妻久久中文字幕网| 麻豆国产97在线/欧美 | 美女免费视频网站| 国产亚洲欧美98| 一二三四社区在线视频社区8| 亚洲av五月六月丁香网| 十八禁人妻一区二区| 国产在线观看jvid| 亚洲成av人片免费观看| 国产精品香港三级国产av潘金莲| 99精品在免费线老司机午夜| 天堂√8在线中文| av免费在线观看网站| 丝袜人妻中文字幕| 一级a爱片免费观看的视频| 美女 人体艺术 gogo| 精品乱码久久久久久99久播| 熟妇人妻久久中文字幕3abv| 两个人视频免费观看高清| 天天躁夜夜躁狠狠躁躁| 老司机午夜十八禁免费视频| 日韩欧美三级三区| 熟妇人妻久久中文字幕3abv| 午夜福利在线观看吧| 亚洲国产精品sss在线观看| 桃色一区二区三区在线观看| 亚洲18禁久久av| 国语自产精品视频在线第100页| 十八禁人妻一区二区| 国产一区在线观看成人免费| 首页视频小说图片口味搜索| 91麻豆精品激情在线观看国产| 欧美乱码精品一区二区三区| 亚洲九九香蕉| 日韩国内少妇激情av| av在线天堂中文字幕|