• <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í)間尺度研究分析
    精品一区在线观看国产| 大陆偷拍与自拍| 国产av精品麻豆| 欧美精品人与动牲交sv欧美| 成人国语在线视频| 免费观看性生交大片5| 欧美人与性动交α欧美软件 | 岛国毛片在线播放| 成人毛片60女人毛片免费| 2022亚洲国产成人精品| 亚洲内射少妇av| 亚洲美女视频黄频| 亚洲欧美成人综合另类久久久| 超色免费av| 国产福利在线免费观看视频| 插逼视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 成人亚洲欧美一区二区av| 欧美激情 高清一区二区三区| 久久久a久久爽久久v久久| 女的被弄到高潮叫床怎么办| 丝袜脚勾引网站| 一本久久精品| 欧美国产精品va在线观看不卡| 国产精品久久久久久av不卡| 国产成人免费无遮挡视频| 欧美国产精品一级二级三级| 丰满乱子伦码专区| 伦理电影免费视频| 美女xxoo啪啪120秒动态图| 免费久久久久久久精品成人欧美视频 | 免费黄频网站在线观看国产| 捣出白浆h1v1| 欧美日韩精品成人综合77777| 夫妻午夜视频| 曰老女人黄片| 亚洲第一av免费看| 亚洲欧美一区二区三区黑人 | 中文字幕av电影在线播放| 最近的中文字幕免费完整| 欧美精品高潮呻吟av久久| 国产白丝娇喘喷水9色精品| 我的女老师完整版在线观看| 高清欧美精品videossex| 男人舔女人的私密视频| videosex国产| 亚洲欧洲国产日韩| 国产高清三级在线| 久久久国产精品麻豆| 午夜av观看不卡| 免费久久久久久久精品成人欧美视频 | 亚洲成人一二三区av| 国产淫语在线视频| 人体艺术视频欧美日本| 国产极品粉嫩免费观看在线| 嫩草影院入口| 九九爱精品视频在线观看| 亚洲久久久国产精品| 成人毛片a级毛片在线播放| 日日撸夜夜添| 男男h啪啪无遮挡| 国产无遮挡羞羞视频在线观看| 久久毛片免费看一区二区三区| av网站免费在线观看视频| 日本vs欧美在线观看视频| 91在线精品国自产拍蜜月| 国产亚洲一区二区精品| 十八禁网站网址无遮挡| 日本色播在线视频| 看免费av毛片| 久久精品国产亚洲av天美| 亚洲国产日韩一区二区| 草草在线视频免费看| 日韩成人av中文字幕在线观看| 一边摸一边做爽爽视频免费| 免费在线观看完整版高清| 亚洲精品乱码久久久久久按摩| 高清av免费在线| 亚洲四区av| 又黄又爽又刺激的免费视频.| 国产有黄有色有爽视频| 有码 亚洲区| 国产片特级美女逼逼视频| 乱码一卡2卡4卡精品| videosex国产| 各种免费的搞黄视频| av国产久精品久网站免费入址| 精品99又大又爽又粗少妇毛片| 成年人免费黄色播放视频| 一区在线观看完整版| 青青草视频在线视频观看| 国产亚洲一区二区精品| 99国产精品免费福利视频| 黄片无遮挡物在线观看| 国产成人免费无遮挡视频| 精品99又大又爽又粗少妇毛片| 大香蕉久久网| av免费在线看不卡| 亚洲av免费高清在线观看| 欧美精品高潮呻吟av久久| 成年美女黄网站色视频大全免费| 99热网站在线观看| 免费日韩欧美在线观看| 色5月婷婷丁香| 亚洲精品,欧美精品| 国产激情久久老熟女| 国产视频首页在线观看| 最后的刺客免费高清国语| av女优亚洲男人天堂| 国产精品国产三级国产专区5o| 18禁动态无遮挡网站| 亚洲国产精品国产精品| 国产成人一区二区在线| 免费观看av网站的网址| 精品国产乱码久久久久久小说| 日韩大片免费观看网站| 久久精品久久精品一区二区三区| a级毛片黄视频| 如何舔出高潮| 亚洲精品一二三| 成人亚洲精品一区在线观看| 免费高清在线观看视频在线观看| 日韩中文字幕视频在线看片| 欧美精品亚洲一区二区| 高清视频免费观看一区二区| 七月丁香在线播放| 欧美xxxx性猛交bbbb| 80岁老熟妇乱子伦牲交| 亚洲熟女精品中文字幕| 中文字幕精品免费在线观看视频 | 97在线视频观看| 水蜜桃什么品种好| 欧美日本中文国产一区发布| 母亲3免费完整高清在线观看 | 国产成人一区二区在线| 日本黄大片高清| 国产熟女午夜一区二区三区| 午夜激情久久久久久久| 男女国产视频网站| 热re99久久国产66热| 国产乱来视频区| 日本欧美视频一区| 久久鲁丝午夜福利片| 男女啪啪激烈高潮av片| av又黄又爽大尺度在线免费看| 久久99热这里只频精品6学生| 99热6这里只有精品| tube8黄色片| 大香蕉久久网| 亚洲国产精品一区二区三区在线| 五月伊人婷婷丁香| 有码 亚洲区| 性色avwww在线观看| 狠狠精品人妻久久久久久综合| 精品亚洲成国产av| 最近中文字幕2019免费版| 日本午夜av视频| 国产精品一区二区在线不卡| 一本大道久久a久久精品| 亚洲av在线观看美女高潮| 成年人免费黄色播放视频| 视频在线观看一区二区三区| 搡女人真爽免费视频火全软件| 久久99热这里只频精品6学生| 久久人人爽人人片av| 精品一区二区三卡| 五月开心婷婷网| 一级片免费观看大全| 免费不卡的大黄色大毛片视频在线观看| 国产精品女同一区二区软件| 日韩 亚洲 欧美在线| 欧美少妇被猛烈插入视频| 国产精品熟女久久久久浪| 超碰97精品在线观看| 少妇被粗大的猛进出69影院 | 超色免费av| 精品一区在线观看国产| 亚洲欧美一区二区三区国产| 国产日韩欧美在线精品| 色网站视频免费| 天天操日日干夜夜撸| 极品人妻少妇av视频| 亚洲精品av麻豆狂野| 秋霞在线观看毛片| 在线观看免费视频网站a站| 亚洲精品美女久久av网站| 蜜臀久久99精品久久宅男| 成人国语在线视频| 精品少妇内射三级| 免费黄网站久久成人精品| 不卡视频在线观看欧美| 欧美日韩视频精品一区| 亚洲欧洲精品一区二区精品久久久 | 少妇熟女欧美另类| 美女中出高潮动态图| 久久久国产精品麻豆| 国产精品无大码| 男女边吃奶边做爰视频| 国产精品久久久久成人av| 中文字幕免费在线视频6| 精品酒店卫生间| 国产熟女午夜一区二区三区| 伦理电影免费视频| 高清欧美精品videossex| 啦啦啦啦在线视频资源| videos熟女内射| 久久久久精品人妻al黑| 亚洲精品视频女| 亚洲精品色激情综合| 国产免费视频播放在线视频| 久久久精品区二区三区| 日韩一区二区视频免费看| 肉色欧美久久久久久久蜜桃| 亚洲国产av新网站| av在线播放精品| 国产精品国产三级国产专区5o| 成年女人在线观看亚洲视频| 亚洲成人av在线免费| 国产成人免费观看mmmm| 日韩欧美精品免费久久| 婷婷色av中文字幕| 日本免费在线观看一区| 日日爽夜夜爽网站| 老熟女久久久| av.在线天堂| 精品第一国产精品| 一个人免费看片子| 美女主播在线视频| 久久国内精品自在自线图片| 日韩视频在线欧美| 欧美人与性动交α欧美精品济南到 | 久久免费观看电影| 9191精品国产免费久久| 亚洲国产精品999| 99热网站在线观看| 国产欧美日韩一区二区三区在线| 国产永久视频网站| 性色avwww在线观看| 日韩在线高清观看一区二区三区| 香蕉丝袜av| 亚洲精品色激情综合| 日日爽夜夜爽网站| 亚洲欧美日韩另类电影网站| 一区二区日韩欧美中文字幕 | av天堂久久9| 久久久久国产精品人妻一区二区| 熟女人妻精品中文字幕| 精品一区二区三卡| 一级片免费观看大全| 夫妻午夜视频| 韩国av在线不卡| 热99久久久久精品小说推荐| a级毛片在线看网站| 一区二区av电影网| 高清欧美精品videossex| 精品视频人人做人人爽| 亚洲av日韩在线播放| 成人国产麻豆网| 中文字幕av电影在线播放| 婷婷成人精品国产| 在线精品无人区一区二区三| 国产女主播在线喷水免费视频网站| 制服诱惑二区| 91精品国产国语对白视频| 欧美成人午夜免费资源| 在线天堂最新版资源| 一边亲一边摸免费视频| 性色av一级| 亚洲精品美女久久久久99蜜臀 | 免费在线观看黄色视频的| 各种免费的搞黄视频| 春色校园在线视频观看| 精品亚洲乱码少妇综合久久| 十八禁网站网址无遮挡| 久久精品久久久久久噜噜老黄| 日日摸夜夜添夜夜爱| 久久久久久久国产电影| 全区人妻精品视频| 日韩一本色道免费dvd| 日韩精品有码人妻一区| 草草在线视频免费看| 大片电影免费在线观看免费| 免费日韩欧美在线观看| 亚洲国产看品久久| 午夜日本视频在线| 久久午夜综合久久蜜桃| 亚洲欧美日韩另类电影网站| 亚洲欧美成人精品一区二区| 69精品国产乱码久久久| 久久精品国产亚洲av天美| 美女中出高潮动态图| 久久久久精品人妻al黑| 欧美日韩一区二区视频在线观看视频在线| 国产精品久久久久久久电影| 中文天堂在线官网| 亚洲色图 男人天堂 中文字幕 | 亚洲 欧美一区二区三区| 另类精品久久| 韩国av在线不卡| 欧美国产精品va在线观看不卡| 天美传媒精品一区二区| 男女高潮啪啪啪动态图| 自拍欧美九色日韩亚洲蝌蚪91| 欧美人与性动交α欧美精品济南到 | 日韩中字成人| 1024视频免费在线观看| 午夜福利在线观看免费完整高清在| www.av在线官网国产| 亚洲精品久久成人aⅴ小说| 97在线人人人人妻| 超碰97精品在线观看| 狂野欧美激情性bbbbbb| 国产在视频线精品| 精品少妇内射三级| 免费观看性生交大片5| 久久影院123| 又粗又硬又长又爽又黄的视频| 一区二区av电影网| 欧美日韩视频精品一区| 午夜av观看不卡| 国产精品一国产av| 制服诱惑二区| 成人黄色视频免费在线看| 久久99精品国语久久久| 日韩av不卡免费在线播放| 国产永久视频网站| 一级,二级,三级黄色视频| 国产精品 国内视频| 久久久亚洲精品成人影院| 22中文网久久字幕| 日韩人妻精品一区2区三区| 久久久久久久久久久久大奶| 边亲边吃奶的免费视频| 寂寞人妻少妇视频99o| 交换朋友夫妻互换小说| 人成视频在线观看免费观看| 亚洲一码二码三码区别大吗| 狠狠婷婷综合久久久久久88av| 黑丝袜美女国产一区| 色94色欧美一区二区| 国产成人精品在线电影| 久久久久久久久久久免费av| 国产免费现黄频在线看| 久久女婷五月综合色啪小说| 日韩中字成人| 日韩不卡一区二区三区视频在线| 日韩制服骚丝袜av| 亚洲成国产人片在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品人妻久久久影院| 亚洲精品色激情综合| 九色成人免费人妻av| 国产精品欧美亚洲77777| 有码 亚洲区| 欧美另类一区| 黑人猛操日本美女一级片| 日韩三级伦理在线观看| 七月丁香在线播放| 最新中文字幕久久久久| 国产精品偷伦视频观看了| 亚洲一码二码三码区别大吗| 制服丝袜香蕉在线| 久久久久久久久久人人人人人人| 国产精品人妻久久久久久| 三上悠亚av全集在线观看| 一级黄片播放器| 一区二区三区乱码不卡18| 久久久久久伊人网av| 在线观看人妻少妇| 麻豆精品久久久久久蜜桃| 亚洲国产最新在线播放| 国产精品国产三级专区第一集| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 精品一区二区免费观看| 高清黄色对白视频在线免费看| 中国美白少妇内射xxxbb| 丁香六月天网| 老司机亚洲免费影院| 一级片'在线观看视频| 精品99又大又爽又粗少妇毛片| 男人爽女人下面视频在线观看| 国产高清国产精品国产三级| 免费观看a级毛片全部| 国产免费福利视频在线观看| 人妻一区二区av| 久久久久久人妻| 秋霞在线观看毛片| av卡一久久| 黄网站色视频无遮挡免费观看| 日韩在线高清观看一区二区三区| 汤姆久久久久久久影院中文字幕| 亚洲在久久综合| h视频一区二区三区| 新久久久久国产一级毛片| 热re99久久国产66热| 久久97久久精品| av天堂久久9| 丰满乱子伦码专区| 狂野欧美激情性bbbbbb| 秋霞伦理黄片| 欧美激情极品国产一区二区三区 | 久久精品国产自在天天线| www.色视频.com| 黄色怎么调成土黄色| av国产精品久久久久影院| 日日摸夜夜添夜夜爱| 亚洲欧美日韩卡通动漫| 国产综合精华液| 一本色道久久久久久精品综合| 国产国语露脸激情在线看| 最黄视频免费看| 欧美人与性动交α欧美精品济南到 | 午夜免费男女啪啪视频观看| 久热这里只有精品99| 国产熟女欧美一区二区| 一级片'在线观看视频| 在线 av 中文字幕| 高清毛片免费看| videos熟女内射| 国产精品.久久久| 日韩人妻精品一区2区三区| av片东京热男人的天堂| av一本久久久久| 久久综合国产亚洲精品| 美女xxoo啪啪120秒动态图| 亚洲精品久久午夜乱码| 欧美性感艳星| 少妇被粗大猛烈的视频| 日韩一本色道免费dvd| 日韩制服骚丝袜av| 99热网站在线观看| 美女国产视频在线观看| 亚洲国产精品专区欧美| 一本—道久久a久久精品蜜桃钙片| 久久久久视频综合| 久久久久视频综合| 777米奇影视久久| 69精品国产乱码久久久| 最新的欧美精品一区二区| 日产精品乱码卡一卡2卡三| 免费高清在线观看视频在线观看| 91aial.com中文字幕在线观看| 亚洲av国产av综合av卡| 伦理电影大哥的女人| 国产男女超爽视频在线观看| 最近中文字幕2019免费版| 成人毛片60女人毛片免费| 美女大奶头黄色视频| 99精国产麻豆久久婷婷| 在线看a的网站| 亚洲精品自拍成人| 国产一区二区在线观看日韩| 视频区图区小说| 国产成人精品婷婷| 久久午夜综合久久蜜桃| 久久人人爽av亚洲精品天堂| 99国产精品免费福利视频| 黄色配什么色好看| 亚洲美女黄色视频免费看| 最黄视频免费看| 制服人妻中文乱码| 久久久久久人妻| 国产片特级美女逼逼视频| 日韩制服丝袜自拍偷拍| 亚洲国产看品久久| 尾随美女入室| 黑丝袜美女国产一区| 一区在线观看完整版| av福利片在线| 欧美精品国产亚洲| 国产日韩一区二区三区精品不卡| 午夜福利视频精品| 丝袜脚勾引网站| 热re99久久精品国产66热6| 婷婷色综合www| 亚洲人与动物交配视频| 香蕉精品网在线| 又黄又爽又刺激的免费视频.| 精品一区二区三卡| 精品视频人人做人人爽| 欧美日本中文国产一区发布| videosex国产| 亚洲av中文av极速乱| 久久这里只有精品19| 多毛熟女@视频| av卡一久久| 水蜜桃什么品种好| 九九爱精品视频在线观看| 国产免费又黄又爽又色| 在线观看免费视频网站a站| 国产精品不卡视频一区二区| 岛国毛片在线播放| 亚洲精品aⅴ在线观看| 国产一区二区三区av在线| 一区二区av电影网| 日韩三级伦理在线观看| 成人18禁高潮啪啪吃奶动态图| 9191精品国产免费久久| av女优亚洲男人天堂| 国产乱来视频区| 涩涩av久久男人的天堂| 欧美变态另类bdsm刘玥| 五月开心婷婷网| 日韩av在线免费看完整版不卡| 国产欧美日韩综合在线一区二区| 国产亚洲午夜精品一区二区久久| 日日摸夜夜添夜夜爱| 国产精品99久久99久久久不卡 | 国产一区亚洲一区在线观看| 成人漫画全彩无遮挡| 久久国产亚洲av麻豆专区| 2021少妇久久久久久久久久久| 建设人人有责人人尽责人人享有的| av黄色大香蕉| 亚洲久久久国产精品| 伊人亚洲综合成人网| 又大又黄又爽视频免费| 高清av免费在线| 91在线精品国自产拍蜜月| 人人妻人人添人人爽欧美一区卜| 美国免费a级毛片| 男女午夜视频在线观看 | 亚洲av成人精品一二三区| 视频中文字幕在线观看| 黄色 视频免费看| 国产色婷婷99| 成年美女黄网站色视频大全免费| 亚洲av福利一区| 嫩草影院入口| 久久久精品区二区三区| 一级黄片播放器| 日韩欧美一区视频在线观看| videosex国产| 国产 一区精品| 欧美国产精品一级二级三级| 亚洲综合精品二区| 天天躁夜夜躁狠狠躁躁| 亚洲人成网站在线观看播放| 国产av精品麻豆| 成人手机av| av一本久久久久| 欧美日韩av久久| 在线观看人妻少妇| 久久久久久久久久久免费av| 99热网站在线观看| a级毛片黄视频| 亚洲精华国产精华液的使用体验| 午夜日本视频在线| 青春草亚洲视频在线观看| 少妇熟女欧美另类| 黑丝袜美女国产一区| 又粗又硬又长又爽又黄的视频| 久久久精品区二区三区| 欧美激情 高清一区二区三区| 男女下面插进去视频免费观看 | 91国产中文字幕| 国产精品三级大全| 97人妻天天添夜夜摸| 精品亚洲成国产av| 亚洲中文av在线| 制服诱惑二区| 精品第一国产精品| 青春草国产在线视频| 亚洲精华国产精华液的使用体验| 高清黄色对白视频在线免费看| 七月丁香在线播放| 午夜福利影视在线免费观看| 日日啪夜夜爽| 免费女性裸体啪啪无遮挡网站| 婷婷色麻豆天堂久久| 国产精品国产三级国产av玫瑰| 亚洲av.av天堂| 9191精品国产免费久久| 高清在线视频一区二区三区| av线在线观看网站| 国产有黄有色有爽视频| av免费在线看不卡| 亚洲国产最新在线播放| 亚洲图色成人| 欧美激情极品国产一区二区三区 | 免费女性裸体啪啪无遮挡网站| 美女国产视频在线观看| 精品久久久精品久久久| 亚洲国产毛片av蜜桃av| 高清黄色对白视频在线免费看| 亚洲综合色网址| 热99国产精品久久久久久7| 精品国产一区二区三区四区第35| 成人毛片60女人毛片免费| 九色成人免费人妻av| 久久久久精品性色| 久久精品aⅴ一区二区三区四区 | 精品国产一区二区三区久久久樱花| 亚洲av成人精品一二三区| 最后的刺客免费高清国语| 又黄又粗又硬又大视频| 国产精品久久久av美女十八| 成人毛片60女人毛片免费| www.熟女人妻精品国产 | 一本久久精品| 一个人免费看片子| 最近中文字幕2019免费版| 国产精品成人在线| 久久影院123| 夫妻午夜视频| 黄色视频在线播放观看不卡| 亚洲,欧美精品.| 男女边摸边吃奶| 一级毛片电影观看| 99九九在线精品视频| 成人漫画全彩无遮挡| 久久久久久久大尺度免费视频| 国产白丝娇喘喷水9色精品| 看免费成人av毛片| 日韩制服骚丝袜av|