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

    基于Morlet小波的呼倫湖流域降水多時(shí)間尺度分析

    2018-05-05 08:52:53韓知明賈克力
    水土保持研究 2018年1期
    關(guān)鍵詞:呼倫湖時(shí)間尺度降水量

    韓知明, 賈克力, 楊 芳, 孫 標(biāo)

    (內(nèi)蒙古農(nóng)業(yè)大學(xué) 水利與土木建筑工程學(xué)院, 呼和浩特 010018)

    濕地降水的時(shí)空分異是濕地形成、演化和分異的基礎(chǔ)[1]。隨著全球氣候變暖、經(jīng)濟(jì)高速增長及城市化水平提高,區(qū)域濕地降水時(shí)空分布發(fā)生再分配,對(duì)湖泊濕地下墊面植被群落類型和分布[2]、濕地蓄水[3]及氣候調(diào)節(jié)作用[4]等產(chǎn)生嚴(yán)重影響。尤其,在我國北方地區(qū),降水趨勢變化將改變區(qū)域生態(tài)系統(tǒng)平衡,影響自然湖泊濕地水環(huán)境質(zhì)量[5-6]?;跐竦厣鷳B(tài)保護(hù)及濕地健康監(jiān)測,有必要對(duì)區(qū)域濕地降水特征進(jìn)行時(shí)間尺度分析,為研究濕地演化及濕地管理提供重要基礎(chǔ)信息。

    小波分析是水文時(shí)間序列分析中的重要工具之一,主要用于研究水文序列的變化趨勢和周期等組成[7-8]?;谛〔ɡ碚摰亩鄷r(shí)間尺度分析可將非平穩(wěn)時(shí)間序列分解為較原始序列平穩(wěn)得多的時(shí)間序列,然后對(duì)分解后的各序列進(jìn)行分別研究[9],繼而可以清晰地提取出隱藏在水文時(shí)間序列中的多種變化周期,充分反映系統(tǒng)在多種不同時(shí)間尺度中的變化趨勢,并能對(duì)水文系統(tǒng)未來發(fā)展趨勢進(jìn)行定性化估計(jì)。例如,Souza Echer等[10]研究了巴西Pelotas百年(1894—1995年)降水量序列變化趨勢,同時(shí)使用小波變換、交叉譜方法分析了該地年降水量與ENSO、準(zhǔn)兩年振蕩和太陽黑子數(shù)的變化關(guān)系。De Jongh等[11]應(yīng)用比利時(shí)Uccle站點(diǎn)105年的降水資料,采用小波分析方法,分別對(duì)該地區(qū)年、月和季節(jié)降水變化進(jìn)行了研究。Narisma等[12]應(yīng)用小波分析方法分析了全球20世紀(jì)降水的區(qū)域突變。以上研究主要是采用長時(shí)間(百年)、單站資料或在較大區(qū)域選擇有限的站點(diǎn)進(jìn)行分析,以揭示降水變化規(guī)律特征。

    目前對(duì)呼倫湖流域降水方面的研究相對(duì)較少。李翀[13]等利用水量平衡模型來模擬呼倫湖水面蒸發(fā)與水量平衡,王志杰[14]等對(duì)呼倫湖水量的現(xiàn)狀及未來氣候下可能的存在狀態(tài)進(jìn)行了分析,萬華偉[15]等對(duì)呼倫湖水面動(dòng)態(tài)變化遙感監(jiān)測及氣候因素驅(qū)動(dòng)進(jìn)行了分析;然而,以湖泊濕地演化為背景,針對(duì)流域降水特征的研究還不完善,因此,本文利用Morlet小波理論對(duì)呼倫湖流域近55 a(1961—2014)降水資料進(jìn)行多時(shí)間尺度分析,以揭示呼倫湖流域降水的時(shí)空多尺度變化特征,為降水預(yù)測和水資源利用提供科學(xué)依據(jù)。

    1 試驗(yàn)材料與方法

    1.1 研究區(qū)概況

    呼倫湖流域位于溫帶半干旱大陸性季風(fēng)氣候區(qū),流域面積約為33 469 km2。該區(qū)域氣候主要呈現(xiàn)為春季干旱多風(fēng),夏季蒸發(fā)強(qiáng)烈降水集中,秋季氣溫降幅較大,冬季寒冷漫長。流域不同地區(qū)年均氣溫為-0.5~0.5℃,積雪期為140 d左右。全年盛行西北風(fēng),年均降雨量為200~300 mm左右,湖水封凍時(shí)間長達(dá)180 d,冰層厚可達(dá)1 m以上[16]。呼倫湖流域位于東北區(qū)域,是中國第五大淡水湖,也是東北地區(qū)第一大湖,呈西北至東南不規(guī)則斜長方形,湖長93 km,平均寬度為32 km,最大湖寬41 km,湖水面積2 339 km2。呼倫湖水系包括呼倫湖、新開河、克魯倫河、烏爾遜河、烏蘭泡、哈拉哈河等主要湖泊和支流[17]。流域內(nèi)主要的雨量站位于坤都冷水文站和阿拉坦額莫勒水文站。

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

    降水?dāng)?shù)據(jù)來源于克魯倫河阿勒坦額勒莫水文站(1961—2014),烏爾遜河坤都冷水文站(1961—2014),滿洲里氣象站(1961—2008),海拉爾壩后水文站(1961—2008)。

    1.3 研究方法

    1.3.1 Morlet小波 Morlet小波不但具有非正交性而且還是Gaussian調(diào)節(jié)的指數(shù)復(fù)值小波,在時(shí)間與頻率的局部化之間具有較好的平衡,因此本文選取復(fù)Morlet小波函數(shù)作為基函數(shù)進(jìn)行小波變換。

    Ψ0(t)=π-1/4eiw0te-t2/2

    (1)

    式中:t為時(shí)間,w0是無量綱頻率。

    1.3.2 小波變換 對(duì)于給定的水文序列和Morlet小波,其連續(xù)小波變換為:

    (2)

    式中:α為尺度因子,反映小波的周期長度;b為時(shí)間因子,反映時(shí)間上的平移; 復(fù)共軛復(fù)數(shù);Wf(a,b)稱小波變換系數(shù)。實(shí)際工作中,時(shí)間序列常常是離散的,則式(2)的離散形式為:

    (3)

    式中:Wf(a,b)能同時(shí)反映時(shí)域參數(shù)a和b頻域參數(shù),a的特性是時(shí)間序列f(t)或f(Δt)通過單位脈沖響應(yīng)的濾波器的輸出。當(dāng)a較小時(shí),對(duì)頻域的分辨率低,對(duì)時(shí)域的分辨率高;當(dāng)a增大時(shí),對(duì)頻域的分辨率高,對(duì)時(shí)域的分辨率低。因此,小波變換實(shí)現(xiàn)了窗口的大小固定,形狀可變的時(shí)頻局部化。

    1.3.3 小波方差 將時(shí)間域上的關(guān)于a的所有小波變換系數(shù)的平方進(jìn)行積分,即為小波方差:

    (4)

    小波方差a隨尺度變化過程稱為小波方差圖,反映了波動(dòng)的能量隨尺度的分布。通過小波方差圖,可以確定一個(gè)水文序列中存在的主要時(shí)間尺度,即主周期。

    2 結(jié)果與分析

    2.1 年降水序列多時(shí)間尺度分析

    根據(jù)呼倫湖流域各地區(qū)年降水小波分析結(jié)果,除了烏爾遜河,其他三地(海拉爾、滿洲里和克魯倫河)的年降水量主要存在3個(gè)明顯的特征時(shí)間尺度,見圖1,海拉爾特征時(shí)間尺度分別是6 a,11 a和27 a,滿洲里分別為10 a,16 a和27 a,克魯倫河分別出現(xiàn)在7 a,13 a和28 a,而烏爾遜河存在4個(gè)明顯的特征時(shí)間尺度,分別是5 a,13 a,20 a和29 a。其中各地能量最強(qiáng)、周期最顯著且存在于整個(gè)統(tǒng)計(jì)時(shí)間內(nèi)的特征時(shí)間尺度分別出現(xiàn)在海拉爾27 a,滿洲里27 a,克魯倫河28 a以及烏爾遜河29 a。同時(shí),各地普遍從20世紀(jì)70 a左右出現(xiàn)明顯的短周期特征,主要在海拉爾6 a,滿洲里10 a,克魯倫河7 a,烏爾遜河5 a。

    圖1呼倫湖流域各地區(qū)年降水小波圖

    2.2 年降水變化過程分析

    小波變換實(shí)部變化過程可以表達(dá)年際降水量豐枯期和突變點(diǎn)變化情況。圖2是呼倫湖流域各地區(qū)年降水量距平序列不同尺度小波系數(shù)實(shí)部變化過程圖。由圖2A可明顯看出海拉爾多年降水過程在6 a特征時(shí)間尺度上平均變化周期為5 a左右,大約經(jīng)歷了12個(gè)高低轉(zhuǎn)換期;在11 a上,平均周期為10 a左右,大約經(jīng)歷了6個(gè)高—低轉(zhuǎn)換期;而在27 a上,降水量的變化周期為20 a左右,大約經(jīng)歷了3個(gè)周期的高低變化。滿洲里多年降水變化過程(圖2B)在10 a特征時(shí)間尺度上,大約經(jīng)歷了7個(gè)高—低轉(zhuǎn)換期;在16 a特征時(shí)間尺度上,大約經(jīng)歷了4個(gè)高—低轉(zhuǎn)換期;而在27 a特征時(shí)間尺度上,大約經(jīng)歷了2個(gè)周期的高低變化。27 a左右的周期變化趨勢中包含著10 a與16 a左右的周期變化。

    克魯倫河和烏爾遜河是呼倫湖流域的主要入湖水系,也是額爾古納河水系的重要組成部分,是研究區(qū)重要的水量補(bǔ)給河流。克魯倫河多年降水變化過程(圖2C)在7 a特征時(shí)間尺度上大約經(jīng)歷了10個(gè)高—低轉(zhuǎn)換期;在11 a特征時(shí)間尺度上,大約經(jīng)歷了5個(gè)高—低轉(zhuǎn)換期;而在28 a特征時(shí)間尺度上大約經(jīng)歷了3個(gè)周期的高低變化。在大周期28 a的趨勢里包含了2個(gè)小周期的變化而且小周期重復(fù)出現(xiàn)次數(shù)有明顯增加趨勢,對(duì)呼倫湖流域湖泊濕地水量平衡變化趨勢有著重要影響。烏爾遜河多年降水變化過程(圖2D)在5 a特征時(shí)間尺度上大約經(jīng)歷了15個(gè)高—低轉(zhuǎn)換期,在13 a特征時(shí)間尺度上大約經(jīng)歷了7個(gè)高—低轉(zhuǎn)換期,在20 a特征時(shí)間尺度上大約經(jīng)歷了3個(gè)周期的高低變化,而在29 a特征時(shí)間尺度上,大約經(jīng)歷了2個(gè)周期的高低變化。其中烏爾遜河具有4個(gè)特征時(shí)間尺度,且在周期的強(qiáng)度與跨度上逐漸增大,平均周期5 a,13 a和20 a的短周期頻率有明顯的增多趨勢。強(qiáng)度最大的29 a特征時(shí)間尺度也體現(xiàn)出小周期所擁有的變化特征,所以烏爾遜河是影響呼倫湖流域最重要的因素之一,且它的降水量變化在很大程度上關(guān)系到呼倫湖未來的發(fā)展走向。

    2.3 季節(jié)性降水變化小波分析

    基于呼倫湖流域特殊的地理位置與特殊的氣候條件,降水主要集中在夏季與秋季,而克魯倫河和烏爾遜河貫穿呼倫湖流域,是呼倫貝爾草原的血液,對(duì)草原生態(tài)起著重要的調(diào)節(jié)作用,因此對(duì)兩條河流的降水在季節(jié)上做了進(jìn)一步的研究。為了更合理更鮮明地顯示出降水規(guī)律統(tǒng)計(jì)時(shí)間以4—9月的累積降水量為夏季總降水,10月—翌年3月的累積降水量為冬季總降水。

    圖2呼倫湖流域各地區(qū)年降水量距平序列不同尺度小波系數(shù)實(shí)部變化過程

    克魯倫河夏季降水量主要存在3個(gè)明顯的特征時(shí)間尺度(圖3),分別是7 a,13 a和28 a,與克魯倫河年降水量的3個(gè)特征時(shí)間尺度一樣,說明夏季降水在全年的降水中占主導(dǎo)。7 a左右的時(shí)間尺度降水周期在1970—2010年之間穩(wěn)定存在,且在1995—2005年這10 a間周期強(qiáng)度加大,與克魯倫河年降水量的變化周期一致。13 a左右的時(shí)間尺度降水周期在1995年以后穩(wěn)定存在,之前表現(xiàn)很弱。28 a左右的時(shí)間尺度的能量最強(qiáng),周期最顯著而且穩(wěn)定存在54 a中。

    圖3克魯倫河夏季降水小波圖及其對(duì)應(yīng)的不同尺度小波系數(shù)實(shí)部變化過程

    克魯倫河冬季降水量主要存在3個(gè)明顯的特征時(shí)間尺度(圖4),分別是5 a,22 a和28 a。22 a左右的時(shí)間尺度降水周期最顯著且穩(wěn)定存在54 a中。28 a左右的時(shí)間尺度的能量最強(qiáng),在1980年后震蕩趨于穩(wěn)定。與克魯倫河年降水量與夏季降水量的特征時(shí)間尺度28 a一致,說明大周期的變化相對(duì)比較穩(wěn)定。5 a左右的時(shí)間尺度降水周期在1980年之前與2000年以后穩(wěn)定存在。短周期降水變化幅度增加顯著,說明冬季降水強(qiáng)度增加,但頻率降低,是冬季極端寒冷天氣頻發(fā)及冬季平均降水過程減少的重要證據(jù)。

    通過分析可以發(fā)現(xiàn)20世紀(jì)90年代初期到21世紀(jì)初期小波系數(shù)的振幅增強(qiáng)且周期逐漸縮短,夏季降水量明顯加大,冬季降水有減小趨勢。這與克魯倫河年降水量的3個(gè)特征時(shí)間尺度上的小波系數(shù)實(shí)部變化過程一致,進(jìn)一步的表明了氣候變化下,研究區(qū)降水更加集中在了夏季(4—9月),年內(nèi)降水分配不均問題加劇。

    圖5顯示烏爾遜河夏季降水量主要存在3個(gè)明顯的特征時(shí)間尺度,分別是5 a,13 a和29 a。5 a左右的時(shí)間尺度降水周期在1975年前后5 a間穩(wěn)定存在,2010年后震蕩強(qiáng)烈。13 a左右的時(shí)間尺度降水周期在20世紀(jì)80年代后穩(wěn)定存在。29 a左右的時(shí)間尺度的能量最強(qiáng),且具有全域性。烏爾遜河夏季降水量變化過程在5 a上,降水量平均變化周期為4 a,大約經(jīng)歷了15個(gè)高—低轉(zhuǎn)換期;在13 a上,平均周變化期為8 a,經(jīng)歷了6個(gè)高—低轉(zhuǎn)換期;而在28 a上,降水量變化周期為18 a,大約經(jīng)歷了2個(gè)周期的高低變化。大周期28 a的變化里明顯包括著5 a與13 a兩個(gè)小周期變化,而且在2010年之后降水量開始顯著增強(qiáng),特征尺度變化的周期明顯縮短,震蕩劇烈。

    圖4克魯倫河冬季降水小波圖及其對(duì)應(yīng)的不同尺度小波系數(shù)實(shí)部變化過程

    圖5烏爾遜河夏季降水小波圖及其對(duì)應(yīng)的不同尺度小波系數(shù)實(shí)部變化過程

    圖6顯示烏爾遜河冬季降水量主要存在2個(gè)明顯的特征時(shí)間尺度,分別是5 a和28 a。5 a左右的時(shí)間尺度降水周期在20世紀(jì)90年代震蕩很弱,但之后開始穩(wěn)定存在。28 a左右的時(shí)間尺度能量最強(qiáng),在54 a中穩(wěn)定存在,且周期有明顯增大的趨勢。烏爾遜河冬季降水變化過程在5 a特征時(shí)間尺度上,平均變化周期為3 a,大約經(jīng)歷了16個(gè)高—低轉(zhuǎn)換期;在28 a上,降水周期為19 a,經(jīng)歷了2個(gè)周期的高低變化。

    自2006年以后烏爾遜河降水量即使在少雨的冬季也有明顯的增多,而且震蕩周期強(qiáng)烈頻繁,所以映射出夏季降水量也有顯著的同樣趨勢和強(qiáng)度。與烏爾遜河夏季降水量變化過程相比,冬季降水的長周期(28 a)變化一致,短周期進(jìn)程加快,冬夏季降水量都有增加趨勢,但降水頻次增加速度顯著低于降水幅度變化,是氣候變暖條件下,北方降水分配方式改變的重要證據(jù),也是夏季澇災(zāi)隱患的先兆。

    圖6烏爾遜河冬季降水小波圖及其對(duì)應(yīng)的不同尺度小波系數(shù)實(shí)部變化過程

    3 討 論

    由本文分析結(jié)果得出,海拉爾6 a的降水量變化周期基本在11 a的降水量變化周期中一致進(jìn)行,而6 a與11 a的降水量變化周期基本又在27 a的降水量變化周期中同步進(jìn)行,即大的周期變化趨勢中隱含著小周期的變化。這與軒偉[18]研究額爾古納河流域近50 a水文氣象要素變化分析中海拉爾存在6 a,12 a,27 a的主周期一致。

    滿洲里市地處呼倫貝爾高原西部邊緣地帶,毗鄰俄羅斯受西伯利亞寒流影響嚴(yán)重,且緯度偏高,遠(yuǎn)離海洋。海拉爾位于大興安嶺西麓的低山丘陵與呼倫貝爾高原東部邊緣的接合地帶,因大興安嶺的屏障作用,使?jié)駶櫟暮Q笮詺夂驁F(tuán)的影響減弱,大部分時(shí)間在西伯利亞高壓的控制之下,降水變化相對(duì)于滿洲里體現(xiàn)出降水量多、分布均勻、變化小等特點(diǎn)。根據(jù)趙慧穎等[19]研究,呼倫湖地區(qū)年降水量的變化是自北向南、由西向東、隨海拔的升高而增加,三者對(duì)年降水量的影響大小依次為經(jīng)度、海拔高度和緯度等特點(diǎn)得到了充分體現(xiàn)。因此,滿洲里的年降水量變化特征時(shí)間尺度較大,小尺度的周期變化雖然存在,但不顯著且不具有代表性,顯示出滿洲里降水量分配不均勻。與海拉爾相比,滿洲里年降水量變化較大,且降水周期規(guī)律沒有海拉爾的穩(wěn)定、明顯。

    呼倫貝爾草原地區(qū)降水主要來源于太平洋夏季風(fēng)帶來的濕潤空氣,由于地處內(nèi)陸又受大興安嶺山脈相隔,水汽來源不足,故年降水量較少。多年平均降水量一般為250~350 mm,降水量自東向西遞減到呼倫湖、克魯倫河一帶減至250 mm以下[20]。因此,河流補(bǔ)給是研究區(qū)湖泊濕地水量的重要保證。根據(jù)黃健[21]對(duì)克魯倫河與烏爾遜河徑流變化的研究,兩條河流約在25 a左右表現(xiàn)出明顯周期性規(guī)律。本文所選的四個(gè)研究區(qū)域的降水量特征時(shí)間大尺度都是在27 a左右,具體海拉爾與滿洲里27 a左右、克魯倫河28 a左右、烏爾遜河29 a左右,體現(xiàn)了四個(gè)研究區(qū)所具有的全域性。

    54 a來呼倫湖流域年降水總體有下降趨勢。由圖2可以看出呼倫湖流域各地區(qū)在27~29 a尺度周期振蕩最強(qiáng)最顯著,其中,海拉爾在27 a尺度上,降水變化的平均周期為10 a;滿洲里在27 a尺度上,降水變化的平均周期為9 a;克魯倫河在28 a尺度上,降水變化的平均周期為9 a;烏爾遜河在29 a尺度上,降水變化的平均周期為10 a。在2013年四個(gè)地區(qū)正處在降水的豐水期,小波實(shí)部變化還有可能繼續(xù)增大或轉(zhuǎn)而下降,但總體情況還處在豐水期。從圖1可知23~30 a的尺度在2014年還沒形成閉合,還有發(fā)展的趨勢,所以在2014年后的一段時(shí)間還將處于豐水期,過后9 a左右的時(shí)間段里將處于枯水期。

    呼倫湖主要水源補(bǔ)給河流克魯倫河與烏爾遜河在季節(jié)性降水變化小波分析中可以十分明顯的看出研究區(qū)降水的季節(jié)性變化:夏秋多雨,冬春少雨雪。降水多集中在6—9月,占常年的80%左右,冬春降雨雪一般在30~50 mm,約占全年降水量的20%左右。其原因是該地區(qū)天氣系統(tǒng)受蒙古高壓、季風(fēng)環(huán)流和地形因素的影響。呼倫湖地區(qū)冬季降水量的增加主要集中在20世紀(jì)90年代以來,而春季降水量的增加則進(jìn)入21世紀(jì)后更為明顯,這正是呼倫湖地區(qū)雪災(zāi)在20世紀(jì)90年代后期的冬季和21世紀(jì)初春季多有發(fā)生的原因所在[22]。

    伴隨氣候變化,湖泊濕地的生態(tài)環(huán)境演變過程發(fā)生了變化,降水量、降水日數(shù)、日極端降水等的變化趨勢[23-24]是氣候變化的重要參數(shù),同時(shí)對(duì)湖泊濕地生態(tài)系統(tǒng)的健康起著非常重要的作用?;谛〔ǚ治鎏峁┑淖赃m應(yīng)時(shí)域和頻域同時(shí)局部化的分析方法,研究降水的時(shí)間變化周期[25-27],能自動(dòng)調(diào)節(jié)時(shí)—頻窗,以適應(yīng)實(shí)際分析需要,從而可克服傳統(tǒng)的滑動(dòng)平均、濾波等方法不具有局部化且診斷突變點(diǎn)存在困難的不足,從多年降水變化過程提取反映呼倫湖流域生態(tài)環(huán)境的變化的有效信息[28],可為呼倫湖流域生態(tài)環(huán)境綜合治理提供參考依據(jù)。

    4 結(jié) 論

    (1) 海拉爾的年降水量周期變化明顯較其他3地穩(wěn)定,而且周期顯著,變化趨勢小,但在20世紀(jì)末開始,周期變化出現(xiàn)顯著浮動(dòng),而且降水量有所減少,周期有所縮短。

    (2) 克魯倫河的年降水量變化周期顯著,而且在大的周期28 a左右下包括著13 a左右周期變化,在13 a周期變化中又包括著7 a左后周期變化。說明降水周期有一定的層次性與規(guī)律性,而且在1995—2010年雖然降水量劇增,但還在小尺度7 a左右變化,所以克魯倫河降水量周期性所受干擾較小,受夏季(4—9月)降水量的直接影響。

    (3) 烏爾遜河的年降水量特征時(shí)間尺度較多于其他3地,但周期規(guī)律并沒有其他3地明顯,且相較于克魯倫河而言,同樣是在降水尺度大周期28 a左右下,小周期規(guī)律缺少一定的層次性。但同樣是受夏季(4—9月)的降水量的影響,與克魯倫河一致。在2008年以后出現(xiàn)了頻繁的高強(qiáng)度降水,這一變化同時(shí)出現(xiàn)在了年、夏季和冬季的周期變化中,所以體現(xiàn)出了在年際方面季節(jié)的時(shí)效性與同步性。

    (4) 根據(jù)周期的變化趨勢可以預(yù)測出四個(gè)研究區(qū)現(xiàn)在正處于豐水期,海拉爾與滿洲里將在2017年迎來突變年,由豐水期轉(zhuǎn)為枯水期;克魯倫河與烏爾遜河將在2019年迎來突變年,由豐水期轉(zhuǎn)為枯水期。利用Morlet小波預(yù)測2018—2035年的降水得出:2018—2026年,呼倫湖流域降水量減少,進(jìn)枯入水期;反之,2027—2035年,呼倫湖流域降水量充沛,進(jìn)入豐水期。

    參考文獻(xiàn):

    [1] 鄧偉,胡金明.濕地水文學(xué)研究進(jìn)展及科學(xué)前沿問題[J].濕地科學(xué),2003,1(1):12-20.

    [2] Price J S, Waddington J M. Advances in Canadian wetland hydrology an biogeochemistry[J]. Hydrological Processes, 2000,14(9):1579-1589.

    [3] 徐衛(wèi)紅,葛德祥,李娜,等.1960—2011年洞庭湖流域降水量變化特征[J].濕地科學(xué),2016,14(1):108-112.

    [4] Hamlin L, Pietroniro A, Prowse T, et al. Application of indexed snowmelt algorithms in a northern wetland regime[J]. Hydrological Processes, 1998,12(10/11):1641-1657.

    [5] 賈文雄,何元慶,李宗省,等.近50年來河西走廊平原區(qū)氣候變化的區(qū)域特征及突變分析[J].地理科學(xué),2008,28(4):525-531.

    [6] 王傳輝,周順武,唐曉萍,等.近48年青藏高原強(qiáng)降水量的時(shí)空分布特征[J].地理科學(xué),2011,31(4):470-477.

    [7] Weng H, Lau K M. Wavelets, period doubling, and time-frequency localization with application to organization of convection over the tropical western Pacific[J]. Journal of the Atmospheric Sciences, 1994,51(17):2523-2541.

    [8] Babu GLS, Srivastava A, Sahana V. Analysis of stability of earthen dams in kachchh region, Gujarat, India[J]. Engineering Geology, 2007,94(3):123-136.

    [9] Cho S E. Effects of spatial variability of soil properties on slope stability[J]. Engineering Geology, 2007,92(3):97-109.

    [10] Souza Echer M P, Echer E, Nordemann D J, et al. Wavelet analysis of a centennial (1895—1994) southern Brazil rainfall series (Pelotas, 31°46′ 19″S 52°20′33″W)[J]. Climatic Change, 2008,87(3):489-497.

    [11] De Jongh I L M, Verhoest N E C, De Troch F P. Analysisof a105-year time series of precipitation observed at Uccle, Belgium[J]. International Journal of Climatology, 2006,26(14):2023-2039.

    [12] Narisma G T, Foley J A, Licker R, et al. Abrupt changes in rainfall during the twentieth century[J]. Geophysical Research Letters, 2007,34(6): L06710, doi:10.1029/2006 GL028628.

    [13] 李翀,馬巍,葉柏生,等.呼倫湖水面蒸發(fā)及水量平衡估計(jì)[J].水文,2006,26(5):41-44.

    [14] 王志杰,李暢游,賈克力,等.呼倫湖生態(tài)系統(tǒng)健康狀況變化趨勢分析[C].世界湖泊大會(huì),2009.

    [15] 萬華偉,康峻,高帥,等.呼倫湖水面動(dòng)態(tài)變化遙感監(jiān)測及氣候因素驅(qū)動(dòng)分析[J].中國環(huán)境科學(xué),2016,36(3):894-898.

    [16] 張娜,烏力吉,劉松濤,等.呼倫湖地區(qū)氣候變化特征及其對(duì)湖泊面積的影響[J].干旱區(qū)資源與環(huán)境,2015(7):192-197.

    [17] Li C, Sun B, Jia K, et al. Multi-band remote sensing based retrieval model and3 D analysis of water depth in Hulun Lake, China[J]. Mathematical & Computer Modelling, 2013,58(3/4):771-781.

    [18] 軒瑋,李翀,趙慧穎,等.額爾古納河流域近50年水文氣象要素變化分析[J].水文,2011,31(5):80-87.

    [19] 趙慧穎,李成才,趙恒和,等.呼倫湖濕地氣候變化及其對(duì)水環(huán)境的影響[J].冰川凍土,2007,29(5):795-801.

    [20] 王文華,朱添.呼倫貝爾草原地區(qū)水文特性分析[J].東北水利水電,2011(7):41-43.

    [21] 黃健.基于小波理論的呼倫湖流域水文序列隨機(jī)分析[D].呼和浩特:內(nèi)蒙古農(nóng)業(yè)大學(xué),2011.

    [22] 趙慧穎,烏力吉,郝文俊.氣候變化對(duì)呼倫湖濕地及其周邊地區(qū)生態(tài)環(huán)境演變的影響[J].生態(tài)學(xué)報(bào),2008,28(3):1064-1071.

    [23] Moberg A, Jones P D. Trends in indices for extremes in daily temperature and precipitation in central and western Europe, 1901—99[J]. International Journal of Climatology, 2005,25(9):1149-1171.

    [24] Baigorria G A, Jones J W, Obrien J J. Understanding rainfall spatial variability in southeast USA at different timescales[J]. International Journal of Climatology, 2007,27(6):749-760.

    [25] Nolin A W, Hall-McKim E A. Frequency modes of monsoon precipitation in Arizona and New Mexico[J]. Monthly Weather Review, 2006,134(12):3774-3781.

    [26] Partal T, ?zgür K. Wavelet and neuro-fuzzy conjunction model for precipitation forecasting[J]. Journal of Hydrology, 2007,342(1/2):199-212.

    [27] Gan T Y, Gobena A K, Wang Q. Precipitation of southwestern Canada: Wavelet, scaling, multifractal analysis, and teleconnection to climate anomalies[J]. Journal of Geophysical Research, 2007,112, D10110, doi:10.1029/2006 J D007157.

    [28] Bhattacharya G, Jana D, Ojha S, et al. Direct search for minimum reliability index of earth slopes[J]. Computers & Geotechnics, 2003,30(6):455-462.

    猜你喜歡
    呼倫湖時(shí)間尺度降水量
    繪制和閱讀降水量柱狀圖
    時(shí)間尺度上非完整系統(tǒng)的Noether準(zhǔn)對(duì)稱性與守恒量
    坐在呼倫湖畔看星星
    遼河(2022年1期)2022-02-14 19:45:39
    坐在呼倫湖畔看星星
    遼河(2022年1期)2022-02-14 05:15:04
    時(shí)間尺度上Lagrange 系統(tǒng)的Hojman 守恒量1)
    交直流混合微電網(wǎng)多時(shí)間尺度協(xié)同控制
    能源工程(2021年1期)2021-04-13 02:06:12
    新巴爾虎寫意(組詩)
    北極光(2019年5期)2019-11-13 04:54:18
    降水量是怎么算出來的
    啟蒙(3-7歲)(2019年8期)2019-09-10 03:09:08
    1988—2017年呼和浩特市降水演變特征分析
    大連市暴雨多時(shí)間尺度研究分析
    亚洲avbb在线观看| 日本欧美视频一区| 久久人妻福利社区极品人妻图片| 久久这里只有精品19| 黄片播放在线免费| 999精品在线视频| 亚洲精品久久成人aⅴ小说| 香蕉丝袜av| 极品人妻少妇av视频| 一级毛片精品| 午夜视频精品福利| 成人免费观看视频高清| 三级毛片av免费| 久久久久久久久免费视频了| 国产成人欧美| 国产三级黄色录像| 高潮久久久久久久久久久不卡| 成年版毛片免费区| avwww免费| 亚洲七黄色美女视频| 免费在线观看亚洲国产| 麻豆一二三区av精品| 精品一区二区三区视频在线观看免费 | 丝袜美腿诱惑在线| 桃色一区二区三区在线观看| 真人做人爱边吃奶动态| 婷婷精品国产亚洲av在线| 亚洲精品国产精品久久久不卡| 欧美乱妇无乱码| 国产av一区在线观看免费| 激情在线观看视频在线高清| 精品国产一区二区三区四区第35| 亚洲精品中文字幕一二三四区| aaaaa片日本免费| 美女福利国产在线| 久久人人精品亚洲av| 成在线人永久免费视频| 老司机午夜福利在线观看视频| 无限看片的www在线观看| 亚洲精品一区av在线观看| 国产免费现黄频在线看| 国产国语露脸激情在线看| 国产激情欧美一区二区| 99国产极品粉嫩在线观看| 国产精品日韩av在线免费观看 | 看片在线看免费视频| 一进一出抽搐gif免费好疼 | 国产精品99久久99久久久不卡| 国产不卡一卡二| 亚洲欧美一区二区三区黑人| 法律面前人人平等表现在哪些方面| 又黄又爽又免费观看的视频| 日韩成人在线观看一区二区三区| 亚洲黑人精品在线| 欧美久久黑人一区二区| 午夜视频精品福利| 久久人妻福利社区极品人妻图片| 精品久久久久久成人av| 亚洲中文字幕日韩| 人人妻人人爽人人添夜夜欢视频| 国产极品粉嫩免费观看在线| 亚洲一区二区三区不卡视频| 一区福利在线观看| 男人的好看免费观看在线视频 | 99国产综合亚洲精品| 色婷婷av一区二区三区视频| 女性被躁到高潮视频| 美女午夜性视频免费| 亚洲一码二码三码区别大吗| 国产精品永久免费网站| 激情视频va一区二区三区| 日韩欧美三级三区| 国产欧美日韩综合在线一区二区| 欧美乱色亚洲激情| 精品人妻在线不人妻| 中国美女看黄片| 亚洲欧洲精品一区二区精品久久久| 女同久久另类99精品国产91| 久久久水蜜桃国产精品网| 欧美日韩一级在线毛片| 色综合婷婷激情| 亚洲少妇的诱惑av| 欧美日韩乱码在线| 成人18禁在线播放| 身体一侧抽搐| 国产一区二区激情短视频| 在线观看免费午夜福利视频| 美女福利国产在线| 午夜福利欧美成人| 最近最新免费中文字幕在线| 一级毛片女人18水好多| 久久人人精品亚洲av| 51午夜福利影视在线观看| 成熟少妇高潮喷水视频| 夜夜看夜夜爽夜夜摸 | 丰满人妻熟妇乱又伦精品不卡| 亚洲精品一卡2卡三卡4卡5卡| 国产一卡二卡三卡精品| 99热国产这里只有精品6| 99久久99久久久精品蜜桃| 亚洲精品一区av在线观看| 亚洲欧美一区二区三区黑人| 欧洲精品卡2卡3卡4卡5卡区| 亚洲五月色婷婷综合| 人人妻人人爽人人添夜夜欢视频| 色婷婷av一区二区三区视频| 亚洲专区国产一区二区| 成年人免费黄色播放视频| 国产成人一区二区三区免费视频网站| av电影中文网址| 巨乳人妻的诱惑在线观看| 精品福利永久在线观看| av国产精品久久久久影院| 韩国精品一区二区三区| 欧美在线一区亚洲| 欧美久久黑人一区二区| 午夜影院日韩av| 欧美日韩亚洲综合一区二区三区_| 可以免费在线观看a视频的电影网站| 曰老女人黄片| 一二三四社区在线视频社区8| bbb黄色大片| 国产伦一二天堂av在线观看| 亚洲精品一二三| 男人操女人黄网站| 不卡av一区二区三区| 国产蜜桃级精品一区二区三区| 国产精华一区二区三区| 91精品国产国语对白视频| 久久久国产欧美日韩av| 桃红色精品国产亚洲av| netflix在线观看网站| 亚洲国产精品sss在线观看 | 亚洲专区国产一区二区| 亚洲av成人av| 亚洲欧美精品综合久久99| 亚洲精品一卡2卡三卡4卡5卡| 中文字幕人妻丝袜一区二区| 狂野欧美激情性xxxx| 侵犯人妻中文字幕一二三四区| 91精品三级在线观看| 亚洲人成伊人成综合网2020| 久久国产精品影院| 美女高潮到喷水免费观看| 日韩一卡2卡3卡4卡2021年| 在线十欧美十亚洲十日本专区| 国产免费av片在线观看野外av| tocl精华| 女性被躁到高潮视频| 久久中文看片网| 国产极品粉嫩免费观看在线| 黄色成人免费大全| 欧美日韩中文字幕国产精品一区二区三区 | 1024视频免费在线观看| 在线观看舔阴道视频| 国产又色又爽无遮挡免费看| 女警被强在线播放| 岛国在线观看网站| 国产精品二区激情视频| 色综合欧美亚洲国产小说| 欧美成人午夜精品| 香蕉丝袜av| 法律面前人人平等表现在哪些方面| 不卡一级毛片| 久久精品国产亚洲av高清一级| 级片在线观看| 亚洲av成人不卡在线观看播放网| 大型av网站在线播放| 亚洲情色 制服丝袜| 99久久人妻综合| 亚洲av熟女| 国产精品一区二区免费欧美| 99re在线观看精品视频| 久久香蕉激情| 在线观看66精品国产| 19禁男女啪啪无遮挡网站| 可以免费在线观看a视频的电影网站| 亚洲欧美日韩另类电影网站| 免费女性裸体啪啪无遮挡网站| 久久午夜综合久久蜜桃| 日本欧美视频一区| 久久午夜亚洲精品久久| 亚洲色图av天堂| 国产午夜精品久久久久久| 久久草成人影院| 中文亚洲av片在线观看爽| 90打野战视频偷拍视频| 久久人人爽av亚洲精品天堂| 国产真人三级小视频在线观看| 亚洲五月婷婷丁香| 美女国产高潮福利片在线看| 久久精品人人爽人人爽视色| 制服人妻中文乱码| 日本欧美视频一区| 极品人妻少妇av视频| 麻豆成人av在线观看| 淫秽高清视频在线观看| 国产精品一区二区免费欧美| 在线观看免费视频网站a站| 亚洲人成电影免费在线| 黄网站色视频无遮挡免费观看| 精品久久久久久久毛片微露脸| √禁漫天堂资源中文www| 999久久久国产精品视频| av网站免费在线观看视频| 国产精品电影一区二区三区| 国产欧美日韩一区二区三| 国产麻豆69| 侵犯人妻中文字幕一二三四区| 亚洲av电影在线进入| 男女高潮啪啪啪动态图| 嫩草影视91久久| 一边摸一边抽搐一进一小说| 亚洲va日本ⅴa欧美va伊人久久| 99精品欧美一区二区三区四区| 欧美中文日本在线观看视频| 亚洲精品在线美女| a级毛片在线看网站| 亚洲国产欧美网| 国产精品国产av在线观看| x7x7x7水蜜桃| 国产成人欧美| 十八禁人妻一区二区| 香蕉久久夜色| 99在线人妻在线中文字幕| 熟女少妇亚洲综合色aaa.| 999久久久国产精品视频| 中亚洲国语对白在线视频| 精品高清国产在线一区| 少妇被粗大的猛进出69影院| www国产在线视频色| 12—13女人毛片做爰片一| 69av精品久久久久久| 黄色毛片三级朝国网站| 亚洲av美国av| 级片在线观看| avwww免费| 精品无人区乱码1区二区| 国产精品亚洲一级av第二区| cao死你这个sao货| 高清欧美精品videossex| 十八禁网站免费在线| 成人三级做爰电影| 村上凉子中文字幕在线| 国产成人一区二区三区免费视频网站| 999精品在线视频| 亚洲黑人精品在线| 黑人操中国人逼视频| 美女高潮喷水抽搐中文字幕| www.熟女人妻精品国产| 国产xxxxx性猛交| 99久久国产精品久久久| 亚洲国产欧美日韩在线播放| 超色免费av| 热99国产精品久久久久久7| 久久香蕉国产精品| 国产人伦9x9x在线观看| 日韩欧美一区视频在线观看| 亚洲欧美精品综合久久99| 亚洲专区中文字幕在线| 国产av又大| 成年版毛片免费区| 国产成人精品久久二区二区91| 真人做人爱边吃奶动态| 成年人免费黄色播放视频| 高清在线国产一区| 美女大奶头视频| 视频区图区小说| 激情视频va一区二区三区| 欧美在线一区亚洲| 欧美日韩瑟瑟在线播放| 久久人妻熟女aⅴ| 日日爽夜夜爽网站| 日本黄色日本黄色录像| 精品久久久久久久久久免费视频 | 女人被狂操c到高潮| 不卡一级毛片| 久久九九热精品免费| 亚洲午夜精品一区,二区,三区| 亚洲性夜色夜夜综合| 精品国内亚洲2022精品成人| 人人妻人人爽人人添夜夜欢视频| 久久影院123| 亚洲片人在线观看| 人妻久久中文字幕网| 国产av在哪里看| 久9热在线精品视频| 九色亚洲精品在线播放| 亚洲黑人精品在线| 老熟妇乱子伦视频在线观看| 高清黄色对白视频在线免费看| 9191精品国产免费久久| 免费高清在线观看日韩| 看黄色毛片网站| 天堂动漫精品| 精品国产乱子伦一区二区三区| 国产精品国产av在线观看| 成人三级做爰电影| 国产一区在线观看成人免费| 国产97色在线日韩免费| 午夜免费成人在线视频| 日本vs欧美在线观看视频| 成人av一区二区三区在线看| 中文欧美无线码| 侵犯人妻中文字幕一二三四区| 日日干狠狠操夜夜爽| 国产熟女xx| 视频区欧美日本亚洲| 怎么达到女性高潮| 老汉色av国产亚洲站长工具| 人妻丰满熟妇av一区二区三区| 精品国产美女av久久久久小说| 校园春色视频在线观看| 日韩 欧美 亚洲 中文字幕| 欧美丝袜亚洲另类 | 丝袜美腿诱惑在线| 香蕉久久夜色| 一边摸一边抽搐一进一小说| 一级a爱视频在线免费观看| а√天堂www在线а√下载| 亚洲av成人一区二区三| 午夜福利免费观看在线| 午夜日韩欧美国产| 99精国产麻豆久久婷婷| 国产亚洲精品一区二区www| 欧美日韩中文字幕国产精品一区二区三区 | 久久久国产一区二区| 波多野结衣av一区二区av| aaaaa片日本免费| www.自偷自拍.com| 热99re8久久精品国产| 欧美一区二区精品小视频在线| 嫁个100分男人电影在线观看| 亚洲色图av天堂| 久久久久久久精品吃奶| 亚洲av五月六月丁香网| 国产精品 国内视频| 曰老女人黄片| 1024视频免费在线观看| 午夜两性在线视频| 超色免费av| 国产精品亚洲一级av第二区| 久久性视频一级片| 老司机靠b影院| 两性夫妻黄色片| 一区二区三区激情视频| a级毛片在线看网站| 9热在线视频观看99| 日日干狠狠操夜夜爽| 超碰成人久久| 香蕉丝袜av| 美国免费a级毛片| 黑丝袜美女国产一区| 国产亚洲精品综合一区在线观看 | 亚洲av五月六月丁香网| 最近最新免费中文字幕在线| 国产亚洲精品久久久久久毛片| 国产精品久久视频播放| 亚洲三区欧美一区| 国产av又大| 丝袜在线中文字幕| 国产成人av教育| 97碰自拍视频| 久久久精品国产亚洲av高清涩受| 涩涩av久久男人的天堂| 成年女人毛片免费观看观看9| 亚洲自拍偷在线| 欧美黄色淫秽网站| 国产不卡一卡二| 亚洲男人的天堂狠狠| 精品国产一区二区三区四区第35| svipshipincom国产片| 久热这里只有精品99| 成人18禁在线播放| 1024香蕉在线观看| 怎么达到女性高潮| 国产精品九九99| 日日爽夜夜爽网站| 极品教师在线免费播放| 欧美精品啪啪一区二区三区| 色精品久久人妻99蜜桃| 国产精品98久久久久久宅男小说| 亚洲人成电影观看| 丝袜在线中文字幕| 很黄的视频免费| 男女午夜视频在线观看| 欧美一级毛片孕妇| 两个人免费观看高清视频| 高清在线国产一区| 91精品三级在线观看| 亚洲色图综合在线观看| 午夜福利在线免费观看网站| 久久精品亚洲精品国产色婷小说| 日韩av在线大香蕉| 亚洲精品成人av观看孕妇| 在线观看免费日韩欧美大片| 国产欧美日韩一区二区三区在线| 亚洲成a人片在线一区二区| 日韩成人在线观看一区二区三区| 国产激情久久老熟女| 国产又色又爽无遮挡免费看| 国产人伦9x9x在线观看| 久久精品影院6| 一级毛片高清免费大全| 中文字幕另类日韩欧美亚洲嫩草| 欧美日韩精品网址| 青草久久国产| 久久人人97超碰香蕉20202| 三级毛片av免费| 男女午夜视频在线观看| 亚洲精品国产一区二区精华液| 80岁老熟妇乱子伦牲交| 亚洲欧美日韩高清在线视频| 色综合站精品国产| 久久精品亚洲熟妇少妇任你| 一级,二级,三级黄色视频| 女同久久另类99精品国产91| 99香蕉大伊视频| 久久人人97超碰香蕉20202| 操美女的视频在线观看| 午夜福利欧美成人| 久久人妻av系列| 免费高清视频大片| 黑人操中国人逼视频| 午夜亚洲福利在线播放| 精品福利观看| 如日韩欧美国产精品一区二区三区| 成人av一区二区三区在线看| 桃色一区二区三区在线观看| 久久亚洲精品不卡| 黑人欧美特级aaaaaa片| 淫秽高清视频在线观看| 欧美日韩瑟瑟在线播放| 嫁个100分男人电影在线观看| 亚洲成人久久性| 成人手机av| 国产精品98久久久久久宅男小说| 欧美日韩乱码在线| 啦啦啦在线免费观看视频4| 丰满的人妻完整版| 老司机福利观看| 老汉色∧v一级毛片| 99国产精品一区二区蜜桃av| 亚洲午夜理论影院| 日韩三级视频一区二区三区| 99国产精品99久久久久| 亚洲欧美一区二区三区黑人| 午夜两性在线视频| 国产成人一区二区三区免费视频网站| 欧美在线黄色| 免费av中文字幕在线| 深夜精品福利| 亚洲中文日韩欧美视频| 午夜亚洲福利在线播放| 国产精品免费视频内射| 如日韩欧美国产精品一区二区三区| 好男人电影高清在线观看| 国产单亲对白刺激| 精品人妻在线不人妻| 50天的宝宝边吃奶边哭怎么回事| 久久久久久大精品| 免费日韩欧美在线观看| 视频在线观看一区二区三区| 不卡一级毛片| 日本撒尿小便嘘嘘汇集6| 美女午夜性视频免费| 最新美女视频免费是黄的| 国产一区二区在线av高清观看| 欧美日韩国产mv在线观看视频| 99热国产这里只有精品6| 91字幕亚洲| 精品欧美一区二区三区在线| 岛国视频午夜一区免费看| 亚洲中文日韩欧美视频| 精品福利永久在线观看| 老司机靠b影院| 国产区一区二久久| 国产欧美日韩一区二区三| 亚洲精品av麻豆狂野| 免费女性裸体啪啪无遮挡网站| 亚洲专区国产一区二区| 国产极品粉嫩免费观看在线| 午夜福利欧美成人| 看免费av毛片| 中文字幕另类日韩欧美亚洲嫩草| 女同久久另类99精品国产91| 日本wwww免费看| 午夜a级毛片| 人人妻人人添人人爽欧美一区卜| 日韩欧美免费精品| 夜夜看夜夜爽夜夜摸 | 天堂俺去俺来也www色官网| 免费在线观看黄色视频的| 亚洲少妇的诱惑av| 国产人伦9x9x在线观看| 国产精品爽爽va在线观看网站 | 久久久久精品国产欧美久久久| 亚洲精品美女久久av网站| 国产成+人综合+亚洲专区| 这个男人来自地球电影免费观看| 欧美+亚洲+日韩+国产| 欧美日韩亚洲国产一区二区在线观看| 天堂√8在线中文| 无限看片的www在线观看| 久久久久国内视频| 波多野结衣高清无吗| 国产精品久久电影中文字幕| 日韩三级视频一区二区三区| 岛国在线观看网站| 一级作爱视频免费观看| 亚洲国产看品久久| 日韩大尺度精品在线看网址 | 老司机福利观看| 亚洲久久久国产精品| 亚洲一区二区三区不卡视频| 亚洲五月婷婷丁香| 亚洲自拍偷在线| 丰满迷人的少妇在线观看| 亚洲av电影在线进入| 一二三四在线观看免费中文在| 免费在线观看黄色视频的| 欧美日韩视频精品一区| 美女大奶头视频| 淫妇啪啪啪对白视频| 亚洲九九香蕉| 亚洲中文字幕日韩| 首页视频小说图片口味搜索| 精品乱码久久久久久99久播| 亚洲专区字幕在线| 制服诱惑二区| 美女午夜性视频免费| 脱女人内裤的视频| cao死你这个sao货| 91字幕亚洲| 亚洲第一欧美日韩一区二区三区| 日本精品一区二区三区蜜桃| 免费日韩欧美在线观看| av福利片在线| 精品人妻在线不人妻| 看黄色毛片网站| 18美女黄网站色大片免费观看| 久久欧美精品欧美久久欧美| www.熟女人妻精品国产| 日韩视频一区二区在线观看| 亚洲成人免费电影在线观看| 性欧美人与动物交配| 99精品欧美一区二区三区四区| 欧美 亚洲 国产 日韩一| 亚洲五月婷婷丁香| 国产在线观看jvid| 精品国内亚洲2022精品成人| 少妇裸体淫交视频免费看高清 | 久久天躁狠狠躁夜夜2o2o| 日本免费一区二区三区高清不卡 | 欧美成人午夜精品| 最新在线观看一区二区三区| 日韩欧美三级三区| 久久精品亚洲熟妇少妇任你| 99精品久久久久人妻精品| 亚洲精品久久成人aⅴ小说| 精品高清国产在线一区| www.www免费av| 色在线成人网| xxxhd国产人妻xxx| 夜夜爽天天搞| 精品人妻1区二区| 国产精品一区二区免费欧美| 亚洲国产精品一区二区三区在线| 欧美一区二区精品小视频在线| 午夜成年电影在线免费观看| 天天添夜夜摸| 亚洲精品在线美女| 校园春色视频在线观看| 国产精品爽爽va在线观看网站 | 脱女人内裤的视频| 日韩人妻精品一区2区三区| 99久久人妻综合| 久久国产精品人妻蜜桃| 久久精品国产亚洲av高清一级| 纯流量卡能插随身wifi吗| 在线播放国产精品三级| 久久人妻福利社区极品人妻图片| 免费在线观看影片大全网站| 日本三级黄在线观看| 免费高清在线观看日韩| 亚洲精品在线美女| 亚洲avbb在线观看| 嫩草影院精品99| 麻豆成人av在线观看| 性欧美人与动物交配| 欧美激情极品国产一区二区三区| 777久久人妻少妇嫩草av网站| 操美女的视频在线观看| 亚洲av成人av| 国产精品免费视频内射| 国产成人欧美在线观看| 亚洲精品国产精品久久久不卡| 一级黄色大片毛片| 成年女人毛片免费观看观看9| 亚洲中文av在线| 精品少妇一区二区三区视频日本电影| 亚洲一区高清亚洲精品| 国产黄a三级三级三级人| 在线观看免费视频日本深夜| www国产在线视频色| 女性生殖器流出的白浆| 日本黄色视频三级网站网址| 亚洲avbb在线观看| 俄罗斯特黄特色一大片| 国产高清视频在线播放一区| av天堂久久9| 亚洲七黄色美女视频| 神马国产精品三级电影在线观看 | 午夜亚洲福利在线播放| 男人舔女人的私密视频|