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

    1961—2017年黑龍江省蒸發(fā)量演變特征及其與氣候因子的關(guān)系*

    2021-07-26 09:23:36李秀芬姜麗霞李險(xiǎn)峰朱海霞宮麗娟趙慧穎
    氣象 2021年6期
    關(guān)鍵詞:氣候因子距平蒸發(fā)量

    李秀芬 姜麗霞 李險(xiǎn)峰 趙 放 朱海霞 王 萍 宮麗娟 趙慧穎

    1 黑龍江省氣象科學(xué)研究所,哈爾濱 150030 2 中國(guó)氣象局東北地區(qū)生態(tài)氣象創(chuàng)新開放實(shí)驗(yàn)室,哈爾濱 150030 3 黑龍江省水利科學(xué)研究院,哈爾濱150080 4 哈爾濱市氣象局,哈爾濱 150028

    提 要: 利用黑龍江省80個(gè)站1961—2017年器測(cè)蒸發(fā)量觀測(cè)資料及常規(guī)氣象觀測(cè)資料,采用線性傾向估計(jì)、累積距平、Mann-Kendal突變分析、數(shù)理統(tǒng)計(jì)和Mexican hat小波分析等方法,分析了黑龍江省年和四季器測(cè)蒸發(fā)量的時(shí)空演變特征,并探討了其與氣候因子的關(guān)系。結(jié)果表明:黑龍江省年蒸發(fā)量的空間分布的地理特征明顯,其值隨緯度、經(jīng)度、海拔高度的增加而遞減,遞減率分別為55.4 mm/°N、45.2 mm/°E、88.8 mm/(100 m)。1961—2017年,黑龍江省年蒸發(fā)量呈顯著下降趨勢(shì),降幅達(dá)13.7 mm/(10 a),存在8 a和24 a周期,全省下降趨勢(shì)站點(diǎn)比例達(dá)70.0%,其中62.5%站點(diǎn)的下降趨勢(shì)通過(guò)0.05的顯著性水平檢驗(yàn),遠(yuǎn)超上升站點(diǎn)比例,總體存在“蒸發(fā)悖論”。季節(jié)間對(duì)照發(fā)現(xiàn),春季蒸發(fā)量降幅較大且趨勢(shì)極顯著,存在24 a、準(zhǔn)2 a 周期,有67個(gè)站點(diǎn)表現(xiàn)為下降趨勢(shì),其中44個(gè)站呈顯著下降趨勢(shì)(P<0.05);夏季、秋季的降幅較小且變化不顯著,均存在7 a周期;冬季則表現(xiàn)為小幅不顯著的增加趨勢(shì),存在24、11、2 a周期,有23站冬季蒸發(fā)量呈顯著上升趨勢(shì)。突變檢驗(yàn)發(fā)現(xiàn),年、春季、冬季蒸發(fā)量存在明顯的突變時(shí)間,夏季和冬季則無(wú)明顯突變。年、季節(jié)蒸發(fā)量與平均溫度、風(fēng)速存在正相關(guān)關(guān)系,與相對(duì)濕度存在負(fù)相關(guān)關(guān)系。風(fēng)速顯著下降是導(dǎo)致年蒸發(fā)量顯著減少的主導(dǎo)因素,風(fēng)速顯著下降及增濕明顯的疊加作用,致使春季蒸發(fā)量的下降趨勢(shì)更顯著,而氣候的暖干化使得冬季蒸發(fā)量呈較弱的上升趨勢(shì)。

    引 言

    蒸發(fā)作為陸-氣水循環(huán)系統(tǒng)的紐帶,在地表水資源重新分配中扮演著重要的角色,平均而言,陸地蒸發(fā)產(chǎn)生的降水占全球陸地總降水量的65%(Chahine,1992)。21世紀(jì)以來(lái),由蒸散發(fā)再凝結(jié)形成的降水量增加,大氣內(nèi)循環(huán)活躍程度加大(姜彤等,2020)。蒸發(fā)量的變化影響著全球水分循環(huán)和能量平衡,作為重要的指標(biāo),在水文監(jiān)測(cè)預(yù)測(cè)、水利工程設(shè)計(jì)、水資源評(píng)價(jià)及氣候區(qū)劃等領(lǐng)域應(yīng)用廣泛(莊曉翠等,2010;榮艷淑等,2012;程海濤和慕彩蕓,2009)。蒸發(fā)皿蒸發(fā)量(即器測(cè)蒸發(fā)量)雖然不能代表實(shí)際蒸發(fā)量,但是作為區(qū)域蒸發(fā)能力的一種體現(xiàn),其值的大小及變化趨勢(shì)可以代表某地區(qū)大氣蒸發(fā)能力,同時(shí)也能反映該區(qū)域干濕狀況的變化以及對(duì)區(qū)域氣候變化的影響。因此,研究器測(cè)蒸發(fā)量的變化對(duì)深入了解氣候變化規(guī)律、探討氣候變化成因具有十分重要的意義。器測(cè)蒸發(fā)量資料具有覆蓋范圍廣、累積序列長(zhǎng)、可比性好等優(yōu)勢(shì),為此,許多學(xué)者側(cè)重于器測(cè)蒸發(fā)量開展研究(Peterson et al,1995;邱新法等,2003;Liu et al,2009;2011;石明星等,2018),獲得了大量的成果。眾多研究發(fā)現(xiàn),氣候變暖背景下的全球大部分地區(qū)的器測(cè)蒸發(fā)量隨溫度升高不升反降趨勢(shì)明顯(Peterson et al,1995; Chattopadhyay and Hulme,1997;Cohen et al,2002;Burn and Hesch,2007;Liu et al,2009),即存在“蒸發(fā)悖論”(Roderick and Farquhar,2002)。導(dǎo)致“蒸發(fā)悖論”現(xiàn)象產(chǎn)生的氣候機(jī)制因地域不同而有所差異。在中國(guó),器測(cè)蒸發(fā)量表現(xiàn)為顯著下降趨勢(shì),但其主要影響因子略有差異,曾燕等(2007)認(rèn)為主要影響因子為輻射和氣溫日較差;申雙和和盛瓊(2008)將其歸因于日照時(shí)數(shù)和風(fēng)速的變化;祁添垚等(2015)的研究發(fā)現(xiàn),相對(duì)濕度是器測(cè)蒸發(fā)量變化的關(guān)鍵影響因子。在黃河流域,“蒸發(fā)悖論”具有空間上和時(shí)間上的不一致性,其主導(dǎo)因素是風(fēng)速的明顯減小(馬雪寧等,2012)。在中國(guó)北方地區(qū),則多將蒸發(fā)量的減小歸因于氣溫日較差的減小、風(fēng)速的下降、太陽(yáng)總輻射的減少及水汽壓的降低(Liu et al, 2011;朱紅蕊等,2013;岳元等,2017;楊璐等,2019)。京津冀地區(qū),平均風(fēng)速是其平原區(qū)蒸發(fā)量變化的主導(dǎo)因子,而在山區(qū)和高原地區(qū),日照時(shí)數(shù)則是主導(dǎo)因子(于占江和楊鵬,2018)。在同樣的氣候背景下,部分地區(qū)的蒸發(fā)皿蒸發(fā)量也有增多現(xiàn)象(石明星等,2018)。大量的研究成果表明,無(wú)論蒸發(fā)量的增多還是減少,重要的是在氣候變暖的背景下存在蒸發(fā)量發(fā)生了明顯變化這一事實(shí),但由于研究區(qū)域、研究時(shí)段、站點(diǎn)數(shù)等諸多因素的不同,導(dǎo)致結(jié)論存在一定的偏差,從而凸顯了開展長(zhǎng)時(shí)間序列、精細(xì)化的小區(qū)域的動(dòng)態(tài)研究的重要性。

    黑龍江省作為中國(guó)糧倉(cāng)和東北亞生態(tài)屏障,擔(dān)負(fù)著保障國(guó)家糧食安全和生態(tài)安全的雙重重任。在氣候變化背景下,黑龍江省氣候變暖顯著(于梅等,2009;陳晶,2013;肖冰霜等,2016)的同時(shí),水資源供給條件、分配方式也發(fā)生了相應(yīng)的變化,特別是20世紀(jì)90年代以來(lái),全省旱情自南向北逐漸加強(qiáng)(Zhao et al, 2019),同時(shí)極端降水事件頻繁發(fā)生(朱海霞等,2019)。氣候干濕的劇烈交替,較大地影響著該地區(qū)生態(tài)環(huán)境和社會(huì)經(jīng)濟(jì)的發(fā)展。為加深對(duì)黑龍江省蒸發(fā)能力隨氣候變化演變特征的理解,驗(yàn)證黑龍江省是否存在“蒸發(fā)悖論”,本文綜合考慮熱力、動(dòng)力、濕度等3類因子,建立蒸發(fā)量的多元線性回歸模型,以探討其與氣候因子的關(guān)系。研究結(jié)果有助于揭示區(qū)域水旱災(zāi)害成因,為準(zhǔn)確判識(shí)、預(yù)估旱澇趨勢(shì),合理制定水資源高效利用途徑提供科學(xué)參考依據(jù)。

    1 資料與方法

    1.1 資 料

    本文所用資料來(lái)自黑龍江省80個(gè)常規(guī)氣象站(圖1)1961—2017年的逐日觀測(cè)資料,包括器測(cè)蒸發(fā)量(小型、E-601型)、平均氣溫、最高氣溫、最低氣溫、平均風(fēng)速、日照時(shí)數(shù)、降水量、平均相對(duì)濕度和平均水汽壓等。針對(duì)2002—2017年5—9月E-601B型蒸發(fā)皿蒸發(fā)量資料,按照換算系數(shù)(任芝花等,2002)統(tǒng)一折算成小型蒸發(fā)皿蒸發(fā)量后使用;對(duì)于缺測(cè)漏測(cè)值則通過(guò)建立其與本站氣溫、降水、日照、相對(duì)濕度及風(fēng)速等多元線性回歸模型進(jìn)行插補(bǔ)(祁添垚等,2015)。季節(jié)劃分采用氣象季節(jié),即:春季為3—5月、夏季為6—8月、秋季為9—11月、冬季為12月—次年2月。各站年、季節(jié)蒸發(fā)量通過(guò)相應(yīng)時(shí)段的逐日蒸發(fā)量觀測(cè)值累加獲得;黑龍江省年、季節(jié)蒸發(fā)量為80個(gè)站年、季節(jié)蒸發(fā)量的算術(shù)平均。

    圖1 黑龍江省80個(gè)氣象觀測(cè)站分布 (填色為海拔高度)Fig.1 Locations of 80 meteorological stations in Heilongjiang Province (colored: altitude)

    1.2 研究方法

    1.2.1 線性傾向分析法

    本文采用線性傾向估計(jì)法(魏鳳英,2007;趙慧穎等,2017),通過(guò)建立蒸發(fā)量(EPi)與對(duì)應(yīng)年(ti)建立的一元線性回歸方程,分析年、季節(jié)蒸發(fā)量的變化趨勢(shì),公式如下:

    EPi=a+bti

    (1)

    式中:a為回歸常數(shù),b為回歸系數(shù),通過(guò)最小二乘法確定。當(dāng)b>0(b<0)時(shí),說(shuō)明蒸發(fā)量隨時(shí)間呈上升(下降)趨勢(shì),以10b值的大小表示蒸發(fā)量上升(下降)的傾向程度。其年、季節(jié)變化的顯著性水平通過(guò)t檢驗(yàn)進(jìn)行判斷:α=0.1,變化趨勢(shì)較顯著;α=0.05,變化趨勢(shì)顯著;α=0.01,變化趨勢(shì)極顯著。

    1.2.2 累積距平分析法

    本文通過(guò)累積距平(魏鳳英,2007)曲線的起伏,

    直觀判斷年、季節(jié)蒸發(fā)量的變化趨勢(shì)及持續(xù)性變化,對(duì)于氣候要素(xi)序列,其某一時(shí)刻t的累積距平(Xt)表示為:

    (2)

    (3)

    1.2.3 Mann-Kendall突變檢驗(yàn)法

    采用非參數(shù)Mann-Kendall法(魏鳳英,2007)對(duì)黑龍江省年、季節(jié)器測(cè)蒸發(fā)量進(jìn)行突變檢驗(yàn)。對(duì)具有n個(gè)樣本的時(shí)間序列x的順序x1,x2,…,xn構(gòu)建秩序列Sk,即

    (4)

    式中:xi和xj分別為第i年和第j年年、季節(jié)蒸發(fā)量。

    在時(shí)間序列隨機(jī)獨(dú)立的假定下,定義統(tǒng)計(jì)量

    (5)

    式中:UF1=0;E(Sk)、var(Sk)分別為Sk的均值和方差,計(jì)算公式如下:

    (6)

    重復(fù)上述計(jì)算過(guò)程,按時(shí)間序列x的逆序xn,xn-1,…,x1計(jì)算統(tǒng)計(jì)量UBk,UB1=0。繪制UFk和UBk曲線,若UFk>0或UFk<0,則表明序列呈上升或下降趨勢(shì);若UFk超過(guò)臨界值時(shí),則表明上升或下降趨勢(shì)顯著;若UFk和UBk兩條曲線在臨界線之間有交點(diǎn),一般來(lái)說(shuō)交點(diǎn)對(duì)應(yīng)的時(shí)刻即突變開始的時(shí)間。對(duì)于存在的虛假突變點(diǎn),結(jié)合累積距平曲線剔除。本研究采用α=0.05顯著性水平,其臨界值為±1.96。

    1.2.4 Mexican hat小波分析方法

    本文采用小波分析法來(lái)識(shí)別年、季節(jié)蒸發(fā)量序列在一定時(shí)間尺度下的周期變化特征、小波方差分析該序列的周期尺度。連續(xù)小波變換可以定義為(Zhao et al,2016;馮禹昊和朱江玲,2019):

    (7)

    -∞

    (8)

    2 結(jié)果與分析

    2.1 蒸發(fā)量的地理分布特征

    利用黑龍江省80個(gè)站多年平均年蒸發(fā)量(PE)與地理因子——緯度(N)、經(jīng)度(E)、海拔高度(H)建立多元回歸方程,以描述年蒸發(fā)量的地理分布特征。建立的回歸方程如下:

    PE=9 859.670-55.360N-45.182E-0.888H

    (9)

    式中:復(fù)相關(guān)系數(shù)R=0.801,通過(guò)了0.001的顯著性水平檢驗(yàn)。分析圖2可知,年蒸發(fā)量總體表現(xiàn)為松嫩、三江兩大平原區(qū)高,大、小興安嶺山地及牡丹江半山區(qū)低。其大值區(qū)位于松嫩平原,絕大部分地區(qū)的年蒸發(fā)量大于1 200 mm,最大值出現(xiàn)在其南部地區(qū)(大于1 600 mm);三江平原次之,大部分地區(qū)年蒸發(fā)量為1 200~1 600 mm;大、小興安嶺山地及牡丹江半山區(qū)大部最少(小于1 000 mm)。分析式(9)可知,年蒸發(fā)量隨地理因子的變化增減明顯,在

    圖2 1961—2017年黑龍江省多年 平均年蒸發(fā)量的空間分布Fig.2 Spatial distribution of annual average pan evaporation in Heilongjiang from 1961 to 2017

    同經(jīng)度、同高度下,緯度每增加1°,年蒸發(fā)量減少55.4 mm;在同緯度、同高度下,經(jīng)度每增加1°,年蒸發(fā)量減少45.2 mm;在同經(jīng)度、同緯度下,海拔高度每增高100 m,年蒸發(fā)量減少88.8 mm。

    2.2 蒸發(fā)量時(shí)間變化

    1961—2017年,黑龍江省年蒸發(fā)量變化總體呈波動(dòng)下降的趨勢(shì)(圖3a),降幅為13.7 mm/(10 a),通過(guò)了0.05顯著性水平檢驗(yàn);年蒸發(fā)量的多年均值為1 296.0 mm,最大值為1 507.8 mm,出現(xiàn)在1982年,最小值僅為1 147.3 mm,出現(xiàn)在2013年,差值達(dá)360.5 mm。突變檢驗(yàn)結(jié)果顯示(圖4a),年蒸發(fā)量和曲線在0.05顯著性水平下存在多個(gè)交點(diǎn),其中,1977—1978年的曲線超出了0.05顯著性水平上限,1993—1997年超出了0.05顯著性水平下限,表明年蒸發(fā)量存在明顯的時(shí)間突變,結(jié)合累積距平曲線(圖4b),確定1983年、2001年為突變年。綜合分析圖3和圖4,近57年,年蒸發(fā)量經(jīng)歷了3個(gè)階段變化:1961—1982年為偏多期,68.2%的年份為正距平年;1983—2000年為持續(xù)偏少期,94.4%的年份為負(fù)距平年;2001—2017年為先增后減期,58.8%的年份為負(fù)距平年。

    圖3 1961—2017年黑龍江省年蒸發(fā)皿蒸發(fā)量時(shí)間變化曲線(a)和距平直方圖(b)Fig.3 Time variation curve (a) and anomaly histogram (b) of annual evaporation in Heilongjiang from 1961 to 2017

    圖4 1961—2017年黑龍江省年蒸發(fā)量Mann-Kendall突變檢驗(yàn)(a)和累積距平曲線(b)Fig.4 Mann-Kendall test (a) and accumulated anomaly curve (b) of annual pan evaporation in Heilongjiang from 1961 to 2017

    春季蒸發(fā)量以13.1 mm/(10 a)的幅度呈顯著下降趨勢(shì)(圖5a),通過(guò)0.01的顯著性水平檢驗(yàn)。突變檢驗(yàn)發(fā)現(xiàn)(圖6a),UF和UB曲線在0.05顯著性水平下有1個(gè)交點(diǎn),且20世紀(jì)90年代以來(lái),UF值超過(guò)顯著性水平0.05顯著性水平下限,甚至超過(guò)0.001顯著性水平下限(臨界值為±2.56),表明春季蒸發(fā)量存在明顯的時(shí)間突變,突變開始時(shí)間為1986年,突變前后其均值分別為477.1、430.0 mm,相差47.1 mm。綜合分析圖5a和圖6a,1961—1985年為春季蒸發(fā)量偏多階段,80.0%的年份為正距平年,1986—2017年為偏少期,71.9%的年份為負(fù)距平年。

    夏季、秋季蒸發(fā)量均呈小幅度不顯著的減少趨勢(shì)(圖5b、5c),累積距平曲線特征與年蒸發(fā)量累積距平曲線相似。突變檢驗(yàn)發(fā)現(xiàn)(圖6b、6c),UF和UB線在0.05顯著性水平下均存在多個(gè)交點(diǎn),但UF值均未超過(guò)0.05顯著性水平,表明夏季、秋季蒸發(fā)量不存在明顯的時(shí)間突變。

    圖5 1961—2017年黑龍江省春(a),夏(b),秋(c)和冬(d)季蒸發(fā)皿蒸發(fā)量特征曲線Fig.5 Characteristic curves of seasonal pan evaporation in Heilongjiang from 1961 to 2017 (a) spring, (b) summer, (c) autumn, (d) winter

    圖6 1961—2017年黑龍江省春(a),夏(b),秋(c)和冬(d)季蒸發(fā)量突變檢驗(yàn)Fig.6 Mann-Kendall test curves of spring (a), summer (b), autumn (c) and winter (d) pan evaporation in Heilongjiang 1961—2017年

    冬季蒸發(fā)量呈小幅度不顯著上升趨勢(shì)(圖5d),UF和UB線在0.05顯著性水平下有1個(gè)交點(diǎn)(圖6d),其中,1964年達(dá)顯著性水平下限,1997—2009年,UF值在顯著性水平線上波動(dòng),且多次超過(guò)其上限值,表明冬季蒸發(fā)量存在明顯的時(shí)間突變,突變開始時(shí)間為1986年,突變前后其均值分別為39.2、44.2 mm,相差6.0 mm。其累積距平曲線表現(xiàn)為以1986年為轉(zhuǎn)折,呈現(xiàn)先降后升特征。

    2.3 蒸發(fā)量變化趨勢(shì)的空間分布

    從圖7可以看出,1961—2017年,黑龍江省80個(gè)站年蒸發(fā)量下降趨勢(shì)的范圍較大,其中35個(gè)站的年總蒸發(fā)量顯著下降趨勢(shì),站點(diǎn)主要分布在松嫩平原、三江平原和牡丹江山區(qū);9個(gè)站呈顯著上升趨勢(shì),有4個(gè)站分布在黑河市,其他5個(gè)站分別分布在大興安嶺、哈爾濱、佳木斯、雞西和牡丹江等地。計(jì)算年蒸發(fā)量的傾向率發(fā)現(xiàn),遜克站增幅最大,傾向率達(dá)55.2 mm/(10 a),烏伊嶺站增幅最小,僅為0.3 mm/(10 a);降幅最大和最小站均出現(xiàn)在綏化市,其中,蘭西站降幅達(dá)90.6 mm/(10 a),望奎站僅為0.3 mm/(10 a)。

    器測(cè)蒸發(fā)量的變化趨勢(shì)在不同季節(jié)表現(xiàn)為不同的空間分布特征(圖8)。從圖8a可知,春季蒸發(fā)量下降趨勢(shì)的空間分布特征與圖7基本相似,但范圍略大,44個(gè)站春季蒸發(fā)量呈顯著下降趨勢(shì),集中在松嫩平原、三江平原及牡丹江半山區(qū);僅4個(gè)站蒸發(fā)量呈顯著上升趨勢(shì),分布在大興安嶺北部及黑河中部地區(qū)。夏季蒸發(fā)量下降的區(qū)域明顯縮小(圖8b),僅14個(gè)站蒸發(fā)量呈顯著下降趨勢(shì),集中分布在哈爾濱西部、伊春南部、綏化南部;顯著增加站點(diǎn)有所增加,達(dá)9個(gè)站,分布在黑南部(4站)及呼瑪、克山、富錦、尚志、穆棱等地。秋季蒸發(fā)量下降區(qū)域進(jìn)一步縮小(圖8c),僅9個(gè)站蒸發(fā)量呈顯著下降趨勢(shì),零星分布在中、東部地區(qū);顯著增加站點(diǎn)增加至17個(gè)站,除伊春外各地市均有分布。冬季蒸發(fā)量上升、下降趨勢(shì)站點(diǎn)相當(dāng)(圖8d),上升趨勢(shì)趨強(qiáng),有23個(gè)站蒸發(fā)量呈顯著增加趨勢(shì),主要分布在大、小興安嶺地區(qū)及三江平原北部,僅龍江、明水、鐵力、虎林等4個(gè)站蒸發(fā)量呈明顯下降趨勢(shì)。

    圖7 1961—2017年黑龍江省年總蒸發(fā)量 變化趨勢(shì)的空間分布Fig.7 Spatial distribution of the trends of annual evaporation in Heilongjiang from 1961 to 2017

    圖8 1961—2017年黑龍江省春(a)、夏(b)、秋(c)、冬(d)季器測(cè)蒸發(fā)量變化趨勢(shì)的空間分布Fig.8 Spatial distribution of the trends of spring (a), summer (b), autumn (c) and winter (d) evaporation in Heilongjiang from 1961 to 2017

    2.4 蒸發(fā)量周期特征

    1961—2017年Mexican hat小波分析發(fā)現(xiàn),年蒸發(fā)量在5~10 a尺度上周期規(guī)律明顯(圖9a),期間經(jīng)歷了“偏少—偏多—偏少—偏多—偏少”5個(gè)循環(huán)交替,到2017年等值線也沒(méi)有閉合,說(shuō)明在未來(lái)一段時(shí)間內(nèi)年蒸發(fā)量仍將處于偏少期;且隨著時(shí)間尺度的增加,在21~26 a尺度上,在20世紀(jì)80年代中期存在較明顯的由偏多向偏少的轉(zhuǎn)折。

    分析各季節(jié)蒸發(fā)量的周期特征表現(xiàn)。春季蒸發(fā)量在5 a以下尺度上雖波動(dòng)劇烈,但周期規(guī)律不明顯;在21~27 a尺度上,在20世紀(jì)80年代中期存在明顯的由偏多向偏少轉(zhuǎn)折(圖9b)。從圖9c、9d可以看出,夏季、秋季蒸發(fā)量具有相似的周期規(guī)律,均在6~9 a尺度上經(jīng)歷了“偏少—偏多—偏少—偏多—偏少”5個(gè)循環(huán)交替,且到2017年等值線均未閉合,說(shuō)明在未來(lái)一段時(shí)間內(nèi)夏季、秋季蒸發(fā)量仍將處于偏少期。冬季蒸發(fā)量的周期變化規(guī)律不同于其他三個(gè)季節(jié)(圖9e),在2 a尺度上存在多個(gè)短期周期振蕩;在11 a左右尺度上則經(jīng)歷了“偏少—偏多—偏少”3個(gè)循環(huán)交替;在20~27 a尺度上存在與春季相反的振蕩趨勢(shì),在20世紀(jì)80年代中后期存在較明顯由偏少向偏多的轉(zhuǎn)折。

    圖9 1961—2017年黑龍江省年(a)、春(b)、夏(c)、秋(d)、冬(e)季蒸發(fā)量小波分析Fig.9 Wavelet analysis of pan evaporation annual (a) and spring (b), summer (c), autumn (d), winter (e) in Heilongjiang from 1961 to 2017

    利用小波系數(shù)方差對(duì)黑龍江省蒸發(fā)量的周期進(jìn)行分析發(fā)現(xiàn)(圖10),1961—2017年黑龍江省年蒸發(fā)量存在8 a、24 a周期;春季存在24 a、準(zhǔn)2 a周期;夏季、秋季均存在7 a周期;冬季存在24 a、11 a、2 a周期。

    圖10 同圖9,但為小波方差Fig.10 Same asFig.9, but for wavelet variance

    2.5 蒸發(fā)量與氣候因子的關(guān)系

    綜合分析年、季節(jié)蒸發(fā)量與相應(yīng)時(shí)段各氣候因子的相關(guān)關(guān)系及各氣候因子的時(shí)間變化趨勢(shì)(表1),結(jié)合各氣候因子間的相關(guān)關(guān)系(表略),篩選了與蒸發(fā)量相關(guān)密切且獨(dú)立性較好的平均溫度(T)、相對(duì)濕度(RH)、風(fēng)速(W)等3個(gè)因子,構(gòu)建黑龍江省年、季節(jié)蒸發(fā)量的氣候影響多元線性回歸方程。

    表1 黑龍江省平均各氣候因子的氣候變化趨勢(shì)及與相應(yīng)時(shí)段蒸發(fā)量相關(guān)系數(shù)Table 1 Climate change trend of the average climate factors in Heilongjiang and the correlation coefficient with the evaporation in the corresponding period

    由于各氣候因子與蒸發(fā)量的單位、量綱不盡相同,為準(zhǔn)確描述氣候因子對(duì)蒸發(fā)量的影響程度,先對(duì)其進(jìn)行標(biāo)準(zhǔn)化處理,然后建立多元回歸方程,并計(jì)算各氣候因子對(duì)蒸發(fā)量變化的相對(duì)貢獻(xiàn)率(張嘉琪和任志遠(yuǎn),2014)。結(jié)果見(jiàn)表2。式(10)~(14)中,EPy、EPSp、EPSu、EPA、EPW分別為年和春、夏、秋、冬季全省平均蒸發(fā)量;Ty、TSp、TSu、TA、TW為相應(yīng)時(shí)段的全省平均溫度;RHy、RHSp、RHSu、RHA、RHW為平均相對(duì)濕度;Wy、WSp、WSu、WA、WW為平均風(fēng)速。

    由式(10)~式(14)可知,年、季節(jié)蒸發(fā)量對(duì)平均溫度、風(fēng)速存在正相關(guān)關(guān)系,對(duì)相對(duì)濕度存在負(fù)相關(guān)關(guān)系,且各氣候因子與蒸發(fā)量的偏相關(guān)關(guān)系均通過(guò)了0.05及以上顯著性水平檢驗(yàn)。以相對(duì)貢獻(xiàn)率來(lái)判斷各氣候因子對(duì)蒸發(fā)量的影響程度(表2)。年尺度上,平均溫度、相對(duì)濕度、風(fēng)速對(duì)蒸發(fā)量的影響程度大小排序?yàn)橄鄬?duì)濕度>風(fēng)速>平均溫度,春季為風(fēng)速>相對(duì)濕度>平均溫度,夏季和秋季為相對(duì)濕度>平均溫度>風(fēng)速,冬季為平均溫度>相對(duì)濕度>風(fēng)速。結(jié)合表1分析可得,風(fēng)速顯著下降對(duì)年蒸發(fā)量的負(fù)效應(yīng)遠(yuǎn)大于顯著增溫和不顯著干化產(chǎn)生的正效應(yīng),導(dǎo)致年蒸發(fā)量呈顯著下降趨勢(shì);春季風(fēng)速顯著下降、增濕明顯對(duì)其蒸發(fā)量產(chǎn)生的疊加負(fù)效應(yīng)遠(yuǎn)大于顯著增溫產(chǎn)生的正效應(yīng),致使春季蒸發(fā)量較年蒸發(fā)量的下降趨勢(shì)更明顯;夏季顯著增溫對(duì)其蒸發(fā)量的正效應(yīng)較大程度削弱了風(fēng)速顯著降低產(chǎn)生的負(fù)效應(yīng),一定程度減緩了夏季蒸發(fā)量的減少趨勢(shì);秋季增溫顯著、干旱化明顯對(duì)其蒸發(fā)量產(chǎn)生的正效應(yīng)較大程度削弱了風(fēng)速顯著下降產(chǎn)生的負(fù)效應(yīng),因此秋季蒸發(fā)量的下降趨勢(shì)更不明顯;而冬季則因氣候暖干對(duì)其蒸發(fā)量的正效應(yīng)略大于風(fēng)速顯著下降產(chǎn)生的負(fù)效應(yīng),使得冬季蒸發(fā)量呈較弱的上升趨勢(shì)。

    表2 蒸發(fā)量氣候影響多元回歸方程及各氣候因子的相對(duì)貢獻(xiàn)率Table 2 Climate response equation of evaporation and relative contribution rates of climate factors

    3 結(jié)論與討論

    利用黑龍江省80個(gè)站1961—2017年器測(cè)蒸發(fā)量觀測(cè)資料,運(yùn)用線性傾向估計(jì)、累積距平、Mann-Kendall突變檢驗(yàn)、數(shù)理統(tǒng)計(jì)及Mexican hat小波分析等方法,對(duì)黑龍江省年、季節(jié)器測(cè)蒸發(fā)量的空間分布特征、趨勢(shì)變化、周期規(guī)律進(jìn)行了分析,并探討了其與氣候因子的關(guān)系,結(jié)果表明:

    (1)1961—2017年,黑龍江省年、春季蒸發(fā)皿蒸發(fā)量分別以13.7、13.1 mm/(10 a)幅度呈顯著下降趨勢(shì)。夏季和秋季蒸發(fā)量的降幅分別為2.3 mm/(10 a)和0.3 mm/(10 a),冬季蒸發(fā)量的增幅為0.9 mm/(10 a),此三個(gè)季節(jié)的增減趨勢(shì)均不明顯。年蒸發(fā)量下降趨勢(shì)站點(diǎn)比例達(dá)70.0%,其中,62.5%站點(diǎn)的下降趨勢(shì)通過(guò)0.05的顯著性水平檢驗(yàn),遠(yuǎn)超上升站點(diǎn)比例,驗(yàn)證了增暖顯著的氣候背景下黑龍江省總體存在“蒸發(fā)悖論”現(xiàn)象。季節(jié)蒸發(fā)量的增減趨勢(shì)具有不同的空間分布特征。春季蒸發(fā)量顯著下降趨勢(shì)區(qū)域最大,冬季顯著增加趨勢(shì)區(qū)域最大。

    (2)突變檢驗(yàn)發(fā)現(xiàn),黑龍江省年、春季、冬季蒸發(fā)量存在明顯的突變時(shí)間,夏季和秋季蒸發(fā)量不存在明顯的時(shí)間突變。年蒸發(fā)量存在兩個(gè)突變點(diǎn),分別為1983年和2001年。春季和冬季蒸發(fā)量均存在一個(gè)突變點(diǎn),且突變開始時(shí)間相同,為1986年,但變化趨勢(shì)相反,春季自突變年起蒸發(fā)量呈顯著減少趨勢(shì),冬季則呈顯著增加趨勢(shì)。

    (3)黑龍江省年、季節(jié)蒸發(fā)量的周期規(guī)律明顯。年蒸發(fā)量存在8 a和24 a周期;春季存在24 a、準(zhǔn)2 a 周期;夏季、秋季均存在7 a周期;冬季存在24、11、2 a周期。在未來(lái)一段時(shí)間內(nèi),黑龍江省年和春、夏、秋三季的蒸發(fā)量將繼續(xù)處于偏少階段,而冬季蒸發(fā)量則仍將保持稍多態(tài)勢(shì)。

    (4)年蒸發(fā)量的地理特征明顯,其值隨緯度、經(jīng)度、海拔高度的增減(增高)而遞減。年蒸發(fā)量的最大值出現(xiàn)在松嫩平原南部(大于1 600 mm)。

    (5)年、季節(jié)蒸發(fā)量對(duì)平均溫度、風(fēng)速等的變化存在正相關(guān)關(guān)系,對(duì)相對(duì)濕度的存在負(fù)相關(guān)關(guān)系。年、夏季、秋季相對(duì)濕度對(duì)其蒸發(fā)量變化的貢獻(xiàn)最大,相對(duì)貢獻(xiàn)率分別為49.5%、55.0%、61.2%;風(fēng)速對(duì)春季蒸發(fā)量變化的貢獻(xiàn)最大,相對(duì)貢獻(xiàn)率為38.8%,相對(duì)濕度次之,為35.0%;平均溫度對(duì)冬季蒸發(fā)量變化貢獻(xiàn)最大,相對(duì)貢獻(xiàn)率為57.3%。

    研究表明,蒸發(fā)量的變化是多環(huán)境因子共同作用的結(jié)果(左洪超等,2006;楊司琪等,2019),其相關(guān)關(guān)系存在一定的地區(qū)差異(邱新法等,2003;曾燕等,2007;馬雪寧等,2012;岳元等,2017;于占江和楊鵬,2018)。本研究表明,年、四季蒸發(fā)量的變化趨勢(shì)與風(fēng)速的變化趨勢(shì)一致,這與已有的研究結(jié)論基本吻合(申雙和和盛瓊,2008;馬雪寧等,2012)。年、季節(jié)蒸發(fā)量變化趨勢(shì)的空間分布特征與朱紅蕊等(2013)的研究結(jié)論大體一致;但與曾燕等(2007)、申雙和和盛瓊(2008)、祁添垚等(2015)的研究結(jié)論存在較大的差異,差異產(chǎn)生的可能原因是研究資料的時(shí)間序列不同,上述文獻(xiàn)所用資料的時(shí)段分別為1960—2000年、1957—2001年、1960—2005年,本研究的資料時(shí)段為1961—2017年,因此,導(dǎo)致上述文獻(xiàn)與本研究結(jié)論在要素變化趨勢(shì)、周期特征上不一致。本研究系統(tǒng)分析了器測(cè)蒸發(fā)量的周期規(guī)律,并定量評(píng)價(jià)了其與氣候因子的關(guān)系,這有別于已有的研究成果,對(duì)已有成果是有益的補(bǔ)充。

    猜你喜歡
    氣候因子距平蒸發(fā)量
    青海省果洛州高原鼠兔宜生區(qū)分布及其對(duì)氣候變化的響應(yīng)
    颶風(fēng)Edouard(2014)暖心結(jié)構(gòu)的多資料對(duì)比分析
    近40年阿里地區(qū)云量和氣溫的年際變化
    西藏科技(2018年9期)2018-10-17 05:51:30
    1958—2013年沽源縣蒸發(fā)量變化特征分析
    1981—2010年菏澤市定陶區(qū)蒸發(fā)量變化特征分析
    新疆民豐縣地表水面蒸發(fā)量分析
    甘肅省降水和冰雹天氣氣候分析
    達(dá)孜縣夏秋季大小型蒸發(fā)量特征、影響因子與差異分析
    地球(2016年7期)2016-08-23 03:01:35
    四川水稻產(chǎn)量及其構(gòu)成要素對(duì)不同生育期氣候因子的響應(yīng)分析
    基于距平百分率的萬(wàn)源市氣象干旱灰色預(yù)測(cè)
    午夜激情久久久久久久| 日韩大码丰满熟妇| 亚洲av美国av| 日本91视频免费播放| 一本大道久久a久久精品| 一区二区三区四区激情视频| 国产片内射在线| 桃花免费在线播放| av福利片在线| 久久精品成人免费网站| 美女脱内裤让男人舔精品视频| 国产欧美日韩精品亚洲av| 欧美激情 高清一区二区三区| 午夜免费观看性视频| 久久国产精品男人的天堂亚洲| 久久久国产欧美日韩av| 国产精品久久久av美女十八| av线在线观看网站| 亚洲精品一区蜜桃| 伊人久久大香线蕉亚洲五| 亚洲欧美激情在线| 欧美性长视频在线观看| 在线十欧美十亚洲十日本专区| 91九色精品人成在线观看| 色播在线永久视频| 久久精品国产综合久久久| 亚洲中文日韩欧美视频| 99久久99久久久精品蜜桃| 精品少妇内射三级| 咕卡用的链子| 五月天丁香电影| 下体分泌物呈黄色| 黑人巨大精品欧美一区二区mp4| 男人爽女人下面视频在线观看| 国产亚洲一区二区精品| www.自偷自拍.com| 亚洲欧美日韩高清在线视频 | 国产欧美日韩精品亚洲av| 久久人妻熟女aⅴ| 大香蕉久久成人网| bbb黄色大片| 亚洲欧美精品自产自拍| 超碰成人久久| 人人澡人人妻人| 午夜免费鲁丝| 精品国产乱子伦一区二区三区 | 女人被躁到高潮嗷嗷叫费观| 12—13女人毛片做爰片一| 黄色 视频免费看| 在线亚洲精品国产二区图片欧美| 女人被躁到高潮嗷嗷叫费观| 亚洲欧美激情在线| 老司机靠b影院| 久久人人97超碰香蕉20202| 亚洲伊人色综图| 精品一区二区三卡| 免费在线观看完整版高清| 视频在线观看一区二区三区| 1024香蕉在线观看| 午夜视频精品福利| 黑人巨大精品欧美一区二区mp4| 欧美另类一区| 国产精品偷伦视频观看了| 亚洲激情五月婷婷啪啪| 亚洲精品久久久久久婷婷小说| 两性夫妻黄色片| 亚洲美女黄色视频免费看| 欧美精品一区二区免费开放| 国产精品秋霞免费鲁丝片| 免费不卡黄色视频| 久久久久国产一级毛片高清牌| 国产精品免费大片| 国产人伦9x9x在线观看| 久久久久久久久免费视频了| xxxhd国产人妻xxx| 久久99热这里只频精品6学生| 脱女人内裤的视频| 91精品国产国语对白视频| 一区二区三区乱码不卡18| 欧美日韩av久久| 日韩中文字幕欧美一区二区| 久久久久久久大尺度免费视频| 日韩制服丝袜自拍偷拍| 国产精品影院久久| 国产真人三级小视频在线观看| 午夜免费观看性视频| 亚洲男人天堂网一区| 国产一卡二卡三卡精品| 99re6热这里在线精品视频| 国产欧美亚洲国产| 日本a在线网址| 一级,二级,三级黄色视频| 在线观看免费日韩欧美大片| 欧美亚洲 丝袜 人妻 在线| 日本一区二区免费在线视频| 日韩制服骚丝袜av| 久久久国产一区二区| 青春草亚洲视频在线观看| 亚洲av成人一区二区三| 老司机亚洲免费影院| 亚洲专区国产一区二区| 婷婷色av中文字幕| 欧美日韩亚洲高清精品| 狠狠精品人妻久久久久久综合| 女人高潮潮喷娇喘18禁视频| 日本欧美视频一区| 国产精品国产av在线观看| 亚洲成av片中文字幕在线观看| 99国产综合亚洲精品| 狠狠狠狠99中文字幕| 一级,二级,三级黄色视频| 中文字幕人妻熟女乱码| av片东京热男人的天堂| 国产激情久久老熟女| 一级毛片女人18水好多| 91精品国产国语对白视频| 亚洲精品国产av成人精品| 久久久久网色| 在线观看www视频免费| 久久久久久人人人人人| 精品亚洲成国产av| 精品一区二区三卡| 黄片小视频在线播放| 黄色视频,在线免费观看| 日韩欧美免费精品| 超色免费av| 中亚洲国语对白在线视频| 国产av又大| 欧美日韩福利视频一区二区| 日韩大片免费观看网站| 国产亚洲欧美精品永久| 日本猛色少妇xxxxx猛交久久| www.熟女人妻精品国产| 国产一区二区 视频在线| 人人妻人人添人人爽欧美一区卜| 国产成人精品久久二区二区免费| 91九色精品人成在线观看| 一二三四在线观看免费中文在| 999精品在线视频| av在线app专区| 人妻 亚洲 视频| 岛国毛片在线播放| 满18在线观看网站| 交换朋友夫妻互换小说| av超薄肉色丝袜交足视频| 99精品久久久久人妻精品| av网站在线播放免费| 91精品伊人久久大香线蕉| 免费在线观看视频国产中文字幕亚洲 | 操出白浆在线播放| 黄频高清免费视频| 又大又爽又粗| 国产精品久久久久久人妻精品电影 | 国产高清视频在线播放一区 | 美女午夜性视频免费| h视频一区二区三区| 欧美黑人欧美精品刺激| 成在线人永久免费视频| 一级,二级,三级黄色视频| 国产成人av教育| 91大片在线观看| 99九九在线精品视频| 下体分泌物呈黄色| 亚洲精品在线美女| 国产精品麻豆人妻色哟哟久久| 美女中出高潮动态图| 日本精品一区二区三区蜜桃| 少妇裸体淫交视频免费看高清 | 国产成人a∨麻豆精品| a级毛片黄视频| 另类精品久久| 在线观看人妻少妇| 大香蕉久久网| 交换朋友夫妻互换小说| 热99re8久久精品国产| 欧美日韩亚洲国产一区二区在线观看 | 国产片内射在线| 国产1区2区3区精品| 五月天丁香电影| 十八禁高潮呻吟视频| 亚洲精品美女久久av网站| 欧美精品一区二区免费开放| 一级片'在线观看视频| 精品国产国语对白av| 国产日韩欧美视频二区| 国产日韩欧美亚洲二区| 免费在线观看完整版高清| 永久免费av网站大全| 亚洲伊人久久精品综合| 久久热在线av| 日韩免费高清中文字幕av| 黄色视频在线播放观看不卡| 亚洲伊人色综图| 国产精品国产av在线观看| 免费女性裸体啪啪无遮挡网站| 国产精品秋霞免费鲁丝片| 久久ye,这里只有精品| 免费在线观看视频国产中文字幕亚洲 | 精品人妻在线不人妻| 亚洲九九香蕉| 欧美精品一区二区免费开放| 日韩大片免费观看网站| 黄片播放在线免费| 高清欧美精品videossex| 老熟妇仑乱视频hdxx| 日韩欧美一区视频在线观看| 午夜福利免费观看在线| 中亚洲国语对白在线视频| 国产精品亚洲av一区麻豆| 一级毛片女人18水好多| 黑人巨大精品欧美一区二区mp4| 久久狼人影院| 色94色欧美一区二区| 99国产精品一区二区蜜桃av | 少妇猛男粗大的猛烈进出视频| 黑人巨大精品欧美一区二区mp4| 波多野结衣av一区二区av| 国产欧美亚洲国产| 桃红色精品国产亚洲av| 亚洲精品国产一区二区精华液| 久久精品国产综合久久久| 丰满饥渴人妻一区二区三| av在线老鸭窝| 免费观看av网站的网址| 国产国语露脸激情在线看| 三级毛片av免费| 天堂8中文在线网| 久久久久网色| 亚洲第一av免费看| 俄罗斯特黄特色一大片| 免费在线观看黄色视频的| 亚洲 欧美一区二区三区| 欧美日韩黄片免| tube8黄色片| 亚洲精品美女久久久久99蜜臀| 1024视频免费在线观看| 亚洲av男天堂| 亚洲色图综合在线观看| 超碰97精品在线观看| 国产精品久久久人人做人人爽| 狠狠狠狠99中文字幕| 18禁裸乳无遮挡动漫免费视频| 免费一级毛片在线播放高清视频 | 操美女的视频在线观看| av国产精品久久久久影院| 国产日韩欧美视频二区| 精品人妻在线不人妻| 中文字幕人妻丝袜制服| 精品一区二区三区四区五区乱码| 午夜91福利影院| 国产欧美日韩一区二区三区在线| 黄片小视频在线播放| a在线观看视频网站| 黑人猛操日本美女一级片| 777久久人妻少妇嫩草av网站| 丝瓜视频免费看黄片| 欧美黑人欧美精品刺激| 麻豆av在线久日| 欧美日韩亚洲国产一区二区在线观看 | 狂野欧美激情性xxxx| 亚洲国产欧美网| 91字幕亚洲| 精品久久久久久久毛片微露脸 | av电影中文网址| 精品人妻一区二区三区麻豆| 国产一区二区三区在线臀色熟女 | 法律面前人人平等表现在哪些方面 | 国产亚洲一区二区精品| 欧美午夜高清在线| 麻豆av在线久日| 法律面前人人平等表现在哪些方面 | 亚洲第一av免费看| 精品免费久久久久久久清纯 | 国产免费av片在线观看野外av| 一级,二级,三级黄色视频| 男女国产视频网站| 九色亚洲精品在线播放| 两个人免费观看高清视频| 国产一区二区三区在线臀色熟女 | 正在播放国产对白刺激| 久久久久国产精品人妻一区二区| av视频免费观看在线观看| 天天操日日干夜夜撸| 法律面前人人平等表现在哪些方面 | 婷婷色av中文字幕| 在线永久观看黄色视频| 大型av网站在线播放| 91麻豆av在线| 热99国产精品久久久久久7| 国产亚洲精品第一综合不卡| 大陆偷拍与自拍| e午夜精品久久久久久久| 美女福利国产在线| 丝瓜视频免费看黄片| 国产精品一区二区在线观看99| 欧美黄色淫秽网站| 在线 av 中文字幕| 考比视频在线观看| 这个男人来自地球电影免费观看| 91大片在线观看| 婷婷色av中文字幕| 久久久水蜜桃国产精品网| 99香蕉大伊视频| 色播在线永久视频| 丝袜喷水一区| 久久精品国产亚洲av高清一级| 女人精品久久久久毛片| 免费在线观看视频国产中文字幕亚洲 | 久久亚洲精品不卡| 亚洲av日韩在线播放| 午夜成年电影在线免费观看| 国产亚洲av高清不卡| 中文字幕人妻丝袜一区二区| 国产极品粉嫩免费观看在线| 久久 成人 亚洲| 丁香六月天网| 满18在线观看网站| 青春草亚洲视频在线观看| 久久久久久免费高清国产稀缺| 日韩,欧美,国产一区二区三区| 国产成人av教育| 岛国在线观看网站| 他把我摸到了高潮在线观看 | 免费高清在线观看视频在线观看| 久久久久久久久免费视频了| 午夜精品国产一区二区电影| 国产欧美亚洲国产| 国产成人av教育| 欧美国产精品va在线观看不卡| 国产精品国产av在线观看| 男人操女人黄网站| 大片免费播放器 马上看| 精品亚洲成国产av| 亚洲精华国产精华精| www.精华液| 国产精品久久久久久精品古装| 国产三级黄色录像| 亚洲免费av在线视频| 欧美国产精品va在线观看不卡| 考比视频在线观看| 日本av免费视频播放| 国产免费一区二区三区四区乱码| 热99久久久久精品小说推荐| 制服诱惑二区| 在线观看舔阴道视频| 成人av一区二区三区在线看 | 国产成人一区二区三区免费视频网站| 亚洲欧美精品自产自拍| 久久久精品免费免费高清| 午夜两性在线视频| 精品久久久精品久久久| 多毛熟女@视频| 极品人妻少妇av视频| 精品国产国语对白av| 无遮挡黄片免费观看| 大型av网站在线播放| 丰满迷人的少妇在线观看| 久久久久久久国产电影| 免费女性裸体啪啪无遮挡网站| 各种免费的搞黄视频| 97精品久久久久久久久久精品| 欧美日韩中文字幕国产精品一区二区三区 | 一级毛片女人18水好多| 乱人伦中国视频| 国产精品 国内视频| 脱女人内裤的视频| 两人在一起打扑克的视频| 在线观看免费高清a一片| 日韩 亚洲 欧美在线| 考比视频在线观看| 一区二区三区精品91| 各种免费的搞黄视频| 午夜激情av网站| 成年人免费黄色播放视频| 国产区一区二久久| 亚洲精品成人av观看孕妇| 精品久久久久久久毛片微露脸 | 一本综合久久免费| 人妻 亚洲 视频| 黄色片一级片一级黄色片| 五月天丁香电影| 99re6热这里在线精品视频| 美女主播在线视频| 欧美性长视频在线观看| 欧美黑人精品巨大| 操出白浆在线播放| 热99国产精品久久久久久7| 在线永久观看黄色视频| 精品一区二区三区四区五区乱码| 欧美日韩精品网址| 19禁男女啪啪无遮挡网站| 亚洲国产欧美在线一区| 激情视频va一区二区三区| 成在线人永久免费视频| 伦理电影免费视频| 日韩欧美国产一区二区入口| 日韩人妻精品一区2区三区| 中文字幕高清在线视频| 巨乳人妻的诱惑在线观看| 成人黄色视频免费在线看| 欧美国产精品va在线观看不卡| 久久久久精品人妻al黑| 国产区一区二久久| 建设人人有责人人尽责人人享有的| 国产在线免费精品| 少妇的丰满在线观看| 一级a爱视频在线免费观看| 国产精品.久久久| 亚洲国产中文字幕在线视频| 999精品在线视频| 亚洲av日韩在线播放| 日日摸夜夜添夜夜添小说| 国产av一区二区精品久久| e午夜精品久久久久久久| 黄网站色视频无遮挡免费观看| 女人久久www免费人成看片| 国产片内射在线| 国产亚洲精品第一综合不卡| 一区二区三区精品91| 亚洲专区字幕在线| 宅男免费午夜| 国产区一区二久久| 999久久久精品免费观看国产| 肉色欧美久久久久久久蜜桃| 国产高清国产精品国产三级| 18禁观看日本| 视频在线观看一区二区三区| 9191精品国产免费久久| www.999成人在线观看| √禁漫天堂资源中文www| 韩国精品一区二区三区| 亚洲国产欧美一区二区综合| 国产成人a∨麻豆精品| 久久天堂一区二区三区四区| 久久这里只有精品19| 丰满饥渴人妻一区二区三| 国产免费福利视频在线观看| 国产精品一区二区免费欧美 | 国产麻豆69| 亚洲精品国产av成人精品| av欧美777| 精品国产一区二区久久| 亚洲成人免费电影在线观看| 欧美日韩亚洲综合一区二区三区_| 91老司机精品| 巨乳人妻的诱惑在线观看| 三级毛片av免费| 99热全是精品| 国产有黄有色有爽视频| 久久99热这里只频精品6学生| 午夜福利视频在线观看免费| 久久久国产成人免费| 亚洲五月色婷婷综合| 欧美xxⅹ黑人| 99久久99久久久精品蜜桃| 母亲3免费完整高清在线观看| 亚洲七黄色美女视频| 黄色怎么调成土黄色| videos熟女内射| a在线观看视频网站| 王馨瑶露胸无遮挡在线观看| 久久av网站| 一区二区av电影网| 99热全是精品| 天堂中文最新版在线下载| 久久99热这里只频精品6学生| 80岁老熟妇乱子伦牲交| 狂野欧美激情性xxxx| 三上悠亚av全集在线观看| 又大又爽又粗| 女性被躁到高潮视频| 精品久久久久久久毛片微露脸 | 精品少妇内射三级| 宅男免费午夜| 国产成人影院久久av| 亚洲自偷自拍图片 自拍| 亚洲国产日韩一区二区| 自拍欧美九色日韩亚洲蝌蚪91| av天堂在线播放| 99精品久久久久人妻精品| 久久国产精品影院| 亚洲精品中文字幕一二三四区 | 亚洲av片天天在线观看| 美女高潮喷水抽搐中文字幕| 法律面前人人平等表现在哪些方面 | 免费在线观看完整版高清| 2018国产大陆天天弄谢| 国产一区有黄有色的免费视频| 国产精品秋霞免费鲁丝片| 婷婷成人精品国产| 日韩人妻精品一区2区三区| 国产精品欧美亚洲77777| 91九色精品人成在线观看| 女人久久www免费人成看片| 国产精品一区二区在线不卡| 国产在线观看jvid| 国产成人免费观看mmmm| 三上悠亚av全集在线观看| 日本一区二区免费在线视频| 一级黄色大片毛片| 丝袜美腿诱惑在线| 久久久久久久精品精品| 99精品欧美一区二区三区四区| 亚洲欧洲日产国产| 精品一区二区三卡| 国产成人一区二区三区免费视频网站| 亚洲精品在线美女| 中亚洲国语对白在线视频| 国产成+人综合+亚洲专区| videosex国产| 久久青草综合色| 十八禁网站网址无遮挡| 国产高清国产精品国产三级| 久久精品aⅴ一区二区三区四区| av又黄又爽大尺度在线免费看| 亚洲av成人不卡在线观看播放网 | 黄色a级毛片大全视频| 亚洲,欧美精品.| 搡老乐熟女国产| 久热这里只有精品99| 2018国产大陆天天弄谢| 午夜影院在线不卡| 桃花免费在线播放| 悠悠久久av| 中文欧美无线码| 亚洲中文字幕日韩| 久久久久久久久久久久大奶| 精品卡一卡二卡四卡免费| www.av在线官网国产| 一级毛片女人18水好多| 精品少妇黑人巨大在线播放| 天堂8中文在线网| 久久精品aⅴ一区二区三区四区| 一本色道久久久久久精品综合| 精品亚洲成国产av| 亚洲国产av影院在线观看| 欧美 日韩 精品 国产| 久久性视频一级片| 欧美日韩黄片免| 久久香蕉激情| 女人被躁到高潮嗷嗷叫费观| 欧美日韩一级在线毛片| 久久久久久久大尺度免费视频| 亚洲色图 男人天堂 中文字幕| 12—13女人毛片做爰片一| 69精品国产乱码久久久| 午夜视频精品福利| 精品熟女少妇八av免费久了| 亚洲专区中文字幕在线| 亚洲欧美成人综合另类久久久| 老熟女久久久| 日本vs欧美在线观看视频| 精品亚洲成国产av| 各种免费的搞黄视频| 在线观看人妻少妇| 女人精品久久久久毛片| 无遮挡黄片免费观看| 18禁黄网站禁片午夜丰满| 免费在线观看影片大全网站| 亚洲一卡2卡3卡4卡5卡精品中文| 纯流量卡能插随身wifi吗| 免费观看a级毛片全部| 久久久久久久国产电影| 欧美日韩福利视频一区二区| 成年人午夜在线观看视频| 狠狠婷婷综合久久久久久88av| 亚洲自偷自拍图片 自拍| 狠狠狠狠99中文字幕| 欧美激情久久久久久爽电影 | 欧美在线黄色| 亚洲国产av影院在线观看| 伊人久久大香线蕉亚洲五| 亚洲国产av新网站| 纯流量卡能插随身wifi吗| 操出白浆在线播放| 亚洲全国av大片| 国产区一区二久久| 国产xxxxx性猛交| www.熟女人妻精品国产| 午夜福利一区二区在线看| 黑人猛操日本美女一级片| 国产欧美日韩一区二区三 | 老司机在亚洲福利影院| 亚洲激情五月婷婷啪啪| 国产精品一二三区在线看| 久久性视频一级片| 中文字幕制服av| 亚洲伊人久久精品综合| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品国产一区二区精华液| 午夜福利一区二区在线看| 亚洲欧洲精品一区二区精品久久久| 午夜福利,免费看| 国精品久久久久久国模美| 久久久国产精品麻豆| 真人做人爱边吃奶动态| 国产精品自产拍在线观看55亚洲 | 亚洲国产欧美一区二区综合| 天堂俺去俺来也www色官网| 在线观看免费日韩欧美大片| 丝袜在线中文字幕| 亚洲伊人久久精品综合| 搡老岳熟女国产| 夫妻午夜视频| 亚洲久久久国产精品| xxxhd国产人妻xxx| 人人妻人人澡人人爽人人夜夜| 又黄又粗又硬又大视频| 久久免费观看电影| 欧美亚洲日本最大视频资源| 色婷婷av一区二区三区视频| 巨乳人妻的诱惑在线观看| 国产激情久久老熟女| 宅男免费午夜| 午夜福利视频精品|