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

    1964~2013年青海省不同生態(tài)功能區(qū)蒸發(fā)皿蒸發(fā)量的變化特征

    2016-06-01 03:30:27祁棟林肖宏斌李曉東周萬福蘇文將胡愛軍
    干旱氣象 2016年2期
    關(guān)鍵詞:青海省

    祁棟林,肖宏斌,李曉東,3,周萬福,蘇文將,胡愛軍

    (1.青海省氣象科學(xué)研究所,青?!∥鲗帯?10001;2.青海省防災(zāi)減災(zāi)重點(diǎn)實(shí)驗(yàn)室,青?!∥鲗帯?10001;3.中國氣象局成都高原氣象研究所,四川 成都 610072)

    ?

    1964~2013年青海省不同生態(tài)功能區(qū)蒸發(fā)皿蒸發(fā)量的變化特征

    祁棟林1,2,肖宏斌1,2,李曉東1,2,3,周萬福1,2,蘇文將1,2,胡愛軍1,2

    (1.青海省氣象科學(xué)研究所,青海西寧810001;2.青海省防災(zāi)減災(zāi)重點(diǎn)實(shí)驗(yàn)室,青海西寧810001;3.中國氣象局成都高原氣象研究所,四川成都610072)

    摘要:選用1964~2013年青海省43個(gè)氣象站常規(guī)觀測資料,分析4個(gè)生態(tài)功能區(qū)蒸發(fā)皿蒸發(fā)量的時(shí)空變化特征,并采用完全相關(guān)法進(jìn)行蒸發(fā)皿蒸發(fā)量變化趨勢的成因分析。結(jié)果表明:近50 a來,青海省4個(gè)生態(tài)功能區(qū)蒸發(fā)皿蒸發(fā)量的年和季節(jié)變化特征明顯,柴達(dá)木盆地和三江源區(qū)年蒸發(fā)量整體分別呈顯著下降和上升趨勢,而環(huán)青海湖區(qū)和東部農(nóng)業(yè)區(qū)變化趨勢不明顯;4個(gè)生態(tài)功能區(qū)年蒸發(fā)量均呈先降后升的階段性變化特征,但轉(zhuǎn)換時(shí)間不盡一致。青海省春、夏、秋季和年蒸發(fā)量均從西北向東南逐漸減小,冬季從南向北逐漸減小。柴達(dá)木盆地和東部農(nóng)業(yè)區(qū)年蒸發(fā)量突變時(shí)間分別是1998年和2002年;環(huán)青海湖區(qū)和三江源區(qū)突變不明顯。影響青海省4個(gè)不同生態(tài)功能區(qū)年蒸發(fā)量的主要?dú)庀笠蜃硬煌翰襁_(dá)木盆地為風(fēng)速、日照時(shí)數(shù)、平均氣溫和氣溫日較差,環(huán)青海湖區(qū)為平均氣溫、風(fēng)速和氣溫日較差,東部農(nóng)業(yè)區(qū)為風(fēng)速和相對濕度,平均氣溫、相對濕度和氣溫日較差是三江源區(qū)蒸發(fā)量上升的主要因素。

    關(guān)鍵詞:蒸發(fā)皿蒸發(fā)量;氣候突變;不同生態(tài)功能區(qū);青海省

    引言

    20世紀(jì)是全球近千年來增暖幅度最顯著的時(shí)期,IPCC第5次評估報(bào)告指出,21世紀(jì)全球氣候?qū)⒗^續(xù)變暖[1]。隨著全球氣溫升高,水循環(huán)加快,勢必將改變?nèi)蛩Y源的時(shí)空分布,進(jìn)而影響生態(tài)環(huán)境和社會(huì)經(jīng)濟(jì)發(fā)展[2]。從1950年代開始,青藏高原冬季氣溫表現(xiàn)出逐步變暖趨勢。相比全球和北半球平均氣溫變暖趨勢,青藏高原變暖趨勢出現(xiàn)更早[3]。蒸散(發(fā))作為潛熱通量是決定天氣與氣候的重要因子,是水循環(huán)中最直接受土地利用和氣候變化影響的一項(xiàng),全球性蒸散(發(fā))對大氣環(huán)流和降水均有重要影響[4]。實(shí)際蒸散(發(fā))的測定非常困難,蒸發(fā)皿蒸發(fā)量雖不能直接代表水面蒸發(fā),但與水面蒸發(fā)之間存在很好的相關(guān)關(guān)系,是水文、氣象臺(tái)站常規(guī)觀測項(xiàng)目之一,具有觀測資料累積序列長、可比性好等優(yōu)勢。研究認(rèn)為[5-13],全球范圍內(nèi)云量或氣溶膠增加所引起的輻射量下降以及氣溫日較差和風(fēng)速的下降是蒸發(fā)皿蒸發(fā)量與潛在蒸發(fā)量下降的主要原因,但不同地區(qū)表現(xiàn)有所差異。

    青海省地處青藏高原東北部,全省平均海拔在3 000 m以上,境內(nèi)地貌特征復(fù)雜多樣,各地氣候條件和植被狀況相差較大。根據(jù)地理位置和地貌特征將青海省劃分為4個(gè)生態(tài)功能區(qū)[14]。其中,環(huán)青海湖區(qū)和三江源區(qū)為天然草場區(qū),主要分布著高寒草甸和草原草場,牧草生長季為4~9月,總體上三江源區(qū)的植被狀況稍好于環(huán)青海湖區(qū);東部農(nóng)業(yè)區(qū)主要以農(nóng)作物為主,作物生長季地面植物覆蓋度較好,冬季和春季多為裸地;柴達(dá)木盆地以荒漠區(qū)為主,氣候干燥,年降水量較少,沙化嚴(yán)重,主要分布著一些灌木類植物。由于4個(gè)不同生態(tài)功能區(qū)地貌和氣候特征差異,導(dǎo)致近年來各地氣候變化趨勢有所不同[14-15]。劉蓓[16]利用1961~2003年44個(gè)臺(tái)站資料分析發(fā)現(xiàn),青海省整體及其東北黃河流域區(qū)、柴達(dá)木盆地區(qū)和南部高原區(qū)年蒸發(fā)量呈逐年下降趨勢;申紅艷等[17]研究表明1961~2010年青海省整體及其西北部和東部農(nóng)業(yè)區(qū)年蒸發(fā)量呈明顯下降趨勢,而地處高原南部的三江源區(qū)顯示輕微增加趨勢;李景鑫等[18]研究指出,青海西寧市蒸發(fā)量的年際變化呈顯著下降趨勢。目前,針對青海省整體及其不同生態(tài)功能區(qū)的蒸發(fā)量變化特征以及影響因子的研究較少,且數(shù)據(jù)大多截止于2010年,時(shí)效性和綜合性明顯不足。本文針對1964~2013年青海省4個(gè)不同生態(tài)功能區(qū)蒸發(fā)皿蒸發(fā)量的氣候變化特征和影響氣象因子開展研究,找出其變化差異,以期為今后更好地認(rèn)識(shí)不同生態(tài)功能區(qū)蒸發(fā)皿蒸發(fā)量變化特征提供基礎(chǔ)。

    1資料與方法

    1.1資料

    考慮到氣象臺(tái)站的搬遷,為確保記錄的連續(xù)性、均一性和可靠性,在保證站點(diǎn)數(shù)及反映4個(gè)生態(tài)功能區(qū)差異的情況下,選取了青海省43個(gè)氣象站1964~2013年逐月平均氣溫、氣溫日較差、日照時(shí)數(shù)、平均風(fēng)速、平均相對濕度、降水量和20 cm蒸發(fā)皿蒸發(fā)量等觀測資料。其中,1964~2003年、2004~2013年1~4月和10~12月為20 cm口徑小型蒸發(fā)皿蒸發(fā)量,而2004~2013年5~9月為E601B型蒸發(fā)器蒸發(fā)量。根據(jù)4個(gè)不同生態(tài)功能區(qū)1998~1999年5~9月20 cm口徑小型蒸發(fā)皿月蒸發(fā)量(Y)和E601B型蒸發(fā)器月蒸發(fā)量(x)對比觀測資料,得到兩者線性回歸方程(表1),回歸方程均通過F檢驗(yàn)。

    根據(jù)表1各生態(tài)功能區(qū)的線性回歸方程,將2004~2013年5~9月E601B型蒸發(fā)器蒸發(fā)量換算成20 cm口徑小型蒸發(fā)皿蒸發(fā)量(以下簡稱蒸發(fā)量),從而得到1964~2013年統(tǒng)一標(biāo)準(zhǔn)的蒸發(fā)量時(shí)間序列。青海省4個(gè)生態(tài)功能區(qū):柴達(dá)木盆地(包括9個(gè)站)、環(huán)青海湖區(qū)(包括8個(gè)站)、東部農(nóng)業(yè)區(qū)(包括12個(gè)站)和三江源區(qū)(包括14個(gè)站),站點(diǎn)分布見圖1。季節(jié)劃分:3~5月為春季,6~8月為夏季,9~11月為秋季,12月至次年2月為冬季。

    1.2方法

    采用線性趨勢法、Mann-Kendall檢驗(yàn)法[19]分別進(jìn)行蒸發(fā)量時(shí)間序列的趨勢及突變分析。

    采用完全相關(guān)系數(shù)法[10,20]分析蒸發(fā)量變化的主要?dú)庀笥绊懸蜃印F溆?jì)算公式:

    (1)

    其中,|QC|為完全相關(guān)系數(shù),C1為氣候要素與蒸發(fā)量的相關(guān)系數(shù),C2為氣候要素與時(shí)間的相關(guān)系數(shù)。只有當(dāng)某氣候要素與蒸發(fā)量相關(guān)性較強(qiáng),且隨時(shí)間發(fā)生顯著變化時(shí),它與蒸發(fā)量的完全相關(guān)系數(shù)才比較大,且C1和C2必須通過95%的信度檢驗(yàn)。

    圖1 青海省43個(gè)氣象站分布及生態(tài)功能區(qū)劃分

    表1 青海省4個(gè)生態(tài)功能區(qū)20 cm口徑小型蒸發(fā)皿和E601B型蒸發(fā)器蒸發(fā)量訂正方程

    2結(jié)果分析

    2.1空間分布特征

    圖2給出1964~2013年青海省全年、冬季蒸發(fā)量和年蒸發(fā)量氣候傾向率的空間分布。由圖2a可看出,青海省年蒸發(fā)量呈現(xiàn)出從西北向東南逐漸減小的分布特征,最大值在柴達(dá)木盆地的冷湖站,達(dá)3 082.5 mm,最小值在環(huán)青海湖區(qū)的門源站,為1 143.3 mm,4個(gè)生態(tài)功能區(qū)的年蒸發(fā)量由高到低依次是柴達(dá)木盆地、東部農(nóng)業(yè)區(qū)、環(huán)青海湖地區(qū)和三江源區(qū)。春、夏、秋季蒸發(fā)量的空間分布特征與年蒸發(fā)量相似(圖略),而冬季蒸發(fā)量的空間分布特征明顯不同,表現(xiàn)為由南向北逐漸減小的特征(圖2b)。近50 a來,青海省各區(qū)域年蒸發(fā)量的變化趨勢不同,氣候傾向率大體自西北向東南逐漸增大,東南部表現(xiàn)為上升趨勢,中部、西北部表現(xiàn)為下降趨勢,其中玉樹站上升趨勢最大,為81.1 mm/10 a,其次是久治站,氣候傾向率為71.2 mm/10 a,而茫崖站下降趨勢最大,為-145.7 mm/10 a,中部一帶有一下降大值中心(圖2c),為諾木洪站,氣候傾向率為-140.7 mm/10 a;四季蒸發(fā)量的氣候傾向率分布特征與年蒸發(fā)量較相似(圖略)。

    各氣象站與各自生態(tài)功能區(qū)域的氣候傾向率并不完全一致。就各生態(tài)功能區(qū)內(nèi)氣象站點(diǎn)年蒸發(fā)量的變化趨勢而言(表2),柴達(dá)木盆地9個(gè)代表站中有8個(gè)呈下降趨勢,均通過α=0.05及以上的顯著性檢驗(yàn),下降幅度最大的是茫崖站,氣候傾向率為-145.7 mm/10 a;環(huán)青海湖區(qū)8個(gè)代表站中有3個(gè)呈下降趨勢,其中2個(gè)臺(tái)站通過α=0.05及以上的顯著性檢驗(yàn),下降幅度最大的是野牛溝站,氣候傾向率為-35.4 mm/10 a,而其余5個(gè)上升臺(tái)站中有3個(gè)通過α=0.05及以上的顯著性檢驗(yàn),增幅最大的是剛察站,氣候傾向率為30.5 mm/10 a;東部農(nóng)業(yè)區(qū)12個(gè)代表站中有5個(gè)為下降趨勢,其中3個(gè)通過α=0.01的顯著性檢驗(yàn),下降幅度最大的是西寧站,氣候傾向率為-95.8 mm/10 a,而上升的7個(gè)臺(tái)站中僅有1個(gè)通過α=0.01的顯著性檢驗(yàn),增幅最大的是同仁站,氣候傾向率為53.1 mm/10 a;三江源區(qū)14個(gè)代表站中有11個(gè)為上升趨勢,其中8個(gè)通過α=0.01的顯著性檢驗(yàn),升幅最大的是玉樹站,氣候傾向率為81.1 mm/10 a,下降的3個(gè)臺(tái)站中有2個(gè)通過α=0.05及以上的顯著性檢驗(yàn),降幅最大的是五道梁,氣候傾向率為-31.9 mm/10 a。

    圖2 1964~2013年青海省年(a)、冬季(b)平均蒸發(fā)量的空間分布(a,b,單位:mm)

    地區(qū)站點(diǎn)總數(shù)下降站點(diǎn)數(shù)量信度0.05信度0.01上升站點(diǎn)數(shù)量信度0.05信度0.01柴達(dá)木盆地9835100環(huán)青海湖區(qū)8311521東部農(nóng)業(yè)區(qū)12503701三江源區(qū)143111108

    2.2時(shí)間變化特征

    2.2.1月和季節(jié)變化

    青海省4個(gè)生態(tài)功能區(qū)平均蒸發(fā)量的逐月變化顯示(圖3a),柴達(dá)木盆地、環(huán)青海湖區(qū)和東部農(nóng)業(yè)區(qū)平均蒸發(fā)量的月變化表現(xiàn)為單峰型分布,而三江源區(qū)表現(xiàn)為弱雙峰型分布。1~5月蒸發(fā)量持續(xù)增加,8月后持續(xù)下降,5~8月是各生態(tài)功能區(qū)全年蒸發(fā)量中最大的4個(gè)月,其中,5~8月蒸發(fā)量占全年比例柴達(dá)木盆地最大(55.9%),三江源區(qū)最小(47.4%)。三江源區(qū)、環(huán)青海湖區(qū)和東部農(nóng)業(yè)區(qū)月蒸發(fā)量最大值出現(xiàn)在5月,柴達(dá)木盆地出現(xiàn)在7月;柴達(dá)木盆地、環(huán)青海湖區(qū)和三江源區(qū)最小值出現(xiàn)在1月,東部農(nóng)業(yè)區(qū)出現(xiàn)在12月;其中柴達(dá)木盆地月蒸發(fā)量最大值是最小值的7倍左右,其他3個(gè)生態(tài)功能區(qū)在3~5倍左右。4個(gè)生態(tài)功能區(qū)各氣象站間蒸發(fā)量月變化規(guī)律略有不同,主要表現(xiàn)在最大值和最小值出現(xiàn)的月份不完全一致。

    青海省4個(gè)生態(tài)功能區(qū)蒸發(fā)量存在明顯的季節(jié)變化(圖3b),夏季蒸發(fā)量最大,其次為春季和秋季,冬季蒸發(fā)量最小,可見春、夏季蒸發(fā)量對這4個(gè)生態(tài)功能區(qū)的水循環(huán)起重要作用。另外發(fā)現(xiàn),除冬季外,其他季節(jié)柴達(dá)木盆地的蒸發(fā)量遠(yuǎn)大于其他3個(gè)功能區(qū),這可能與柴達(dá)木盆地干旱氣候有關(guān)。

    2.2.2年際變化

    圖4給出1964~2013年青海省4個(gè)生態(tài)功能區(qū)平均蒸發(fā)量的年變化趨勢。整體上看出,近50 a青海省柴達(dá)木盆地和三江源區(qū)蒸發(fā)量變化趨勢均顯著(通過α=0.001的顯著性檢驗(yàn)),前者呈下降趨勢,后者呈上升趨勢,氣候傾向率分別為-77.9 mm/10 a和30.1 mm/10 a,而環(huán)青海湖區(qū)和東部農(nóng)業(yè)區(qū)蒸發(fā)量沒有明顯變化趨勢,這與劉蓓[16]和申紅艷[17]等的研究結(jié)果不一致。表明青海省平均蒸發(fā)量整體呈下降趨勢,但局部區(qū)域與整體變化形式并不完全同步。

    另外發(fā)現(xiàn),4個(gè)生態(tài)功能區(qū)年蒸發(fā)量多項(xiàng)式擬合曲線均呈先降后升的特征,但轉(zhuǎn)換時(shí)間不完全一致。其中,柴達(dá)木盆地在1990年代中期以前蒸發(fā)量呈明顯下降趨勢,以后進(jìn)入緩慢波動(dòng)狀態(tài)(圖4a);環(huán)青海湖區(qū)在1980年代末以前蒸發(fā)量波動(dòng)下降,進(jìn)入1990年代后轉(zhuǎn)為上升趨勢(圖4b);東部農(nóng)業(yè)區(qū)蒸發(fā)量在1990年代中期以前呈顯著下降趨勢,而后呈顯著上升趨勢(圖4c);三江源區(qū)蒸發(fā)量在1990年代以前波動(dòng)下降,以后至2000年代中期呈顯著上升趨勢,隨后出現(xiàn)快速下降趨勢(圖4d)。

    圖3 1964~2013年青海省4個(gè)生態(tài)功能區(qū)平均蒸發(fā)量的月(a)和季節(jié)(b)變化

    圖4 1964~2013年青海省4個(gè)生態(tài)功能區(qū)蒸發(fā)量的年變化趨勢

    表3是青海省4個(gè)生態(tài)功能區(qū)四季蒸發(fā)量的氣候傾向率和相關(guān)系數(shù)。可以看出,近50 a來,柴達(dá)木盆地春、夏、秋季蒸發(fā)量均呈顯著下降趨勢,冬季呈顯著上升趨勢,其中春、夏季蒸發(fā)量的大幅下降造成其年蒸發(fā)量顯著減?。蝗磪^(qū)四季蒸發(fā)量均呈上升趨勢,其中夏、秋、冬季通過α=0.05及以上的顯著性檢驗(yàn),夏、秋季蒸發(fā)量的上升造成該區(qū)年蒸發(fā)量的顯著增大;環(huán)青海湖區(qū)夏、秋、冬季蒸發(fā)量均呈上升趨勢,而春季表現(xiàn)為下降趨勢,其中僅秋季通過α=0.05的顯著性檢驗(yàn);東部農(nóng)業(yè)區(qū)除夏季外,其他季節(jié)均表現(xiàn)為不顯著下降趨勢。

    從青海省4個(gè)生態(tài)功能區(qū)各月蒸發(fā)量變化趨勢來看(表4),近50 a來,柴達(dá)木盆地各月蒸發(fā)量均呈下降趨勢,線性趨勢率在-12.5~-1.0 mm/10 a之間,除2月外,其他月份均通過α=0.05及以上的顯著性檢驗(yàn),減幅最大的是7月,最小的是2月;三江源區(qū)各月蒸發(fā)量均為上升趨勢,僅7~9月和11月通過α=0.05及以上的顯著性檢驗(yàn),9月增幅最大,為8.0 mm/10 a;環(huán)青海湖區(qū)春季到夏初為下降趨勢,其他月份為上升趨勢,僅9月通過α=0.05的顯著性檢驗(yàn);東部農(nóng)業(yè)區(qū)夏季到秋初為上升趨勢,其他月份均呈下降趨勢,且均未通過顯著性檢驗(yàn)。

    表3 1964~2013年青海省4個(gè)生態(tài)功能區(qū)四季蒸發(fā)量的氣候傾向率(單位:mm/10 a)和相關(guān)系數(shù)

    注:*、**和***分別表示通過α=0.05、0.01和0.001的顯著性檢驗(yàn),下同

    表4 1964~2013年青海省4個(gè)生態(tài)功能區(qū)逐月蒸發(fā)量的氣候傾向率(單位:mm/10 a)和相關(guān)系數(shù)

    2.2.3年代際變化

    從青海省4個(gè)生態(tài)功能區(qū)蒸發(fā)量的年代際距平變化(表5)來看,1970年代除三江源區(qū)外,其他3個(gè)功能區(qū)為正距平,其中柴達(dá)木盆地最大;1980~1990年代,4個(gè)功能區(qū)均為負(fù)距平;2000年代除柴達(dá)木盆地外,其他3個(gè)功能區(qū)為正距平,其中三江源區(qū)最大。可見,在全球變暖背景下,青海省不同生態(tài)功能區(qū)蒸發(fā)量的年代際變化特征差異較大。

    表5 1964~2013年青海省4個(gè)生態(tài)功能區(qū)

    2.3突變特征

    從1964~2013年青海省4個(gè)生態(tài)功能區(qū)年蒸發(fā)量的Mann-Kendall檢驗(yàn)(圖5)可以看出:柴達(dá)木盆地年蒸發(fā)量1964年開始波動(dòng)下降,1975年后有短期上升趨勢,1980年代蒸發(fā)量開始持續(xù)下降,1990年代末下降趨勢超過α=0.05的顯著性水平臨界線,表明蒸發(fā)量減小趨勢明顯(圖5a)。根據(jù)UF和UB曲線的交點(diǎn),確定突變開始時(shí)間是1998年。環(huán)青海湖區(qū)在2005年出現(xiàn)交點(diǎn),但下降趨勢沒有超過α=0.05的顯著性水平臨界線,表明下降不明顯(圖5b)。東部農(nóng)業(yè)區(qū)自1964年以來一直處于波動(dòng)下降趨勢,2004年以后下降趨勢明顯,突變時(shí)間為2002年(圖5c)。三江源區(qū)自1964年開始波動(dòng)上升,1979年以后上升趨勢較顯著,1980年代后期開始轉(zhuǎn)入下降趨勢,在2005~2010年出現(xiàn)3個(gè)交點(diǎn),但UF線未超過α=0.05的顯著性水平臨界線,表明突變不明顯(圖5d)。

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

    蒸發(fā)量作為大氣蒸發(fā)潛力的一個(gè)重要指標(biāo),其影響因素眾多,既有氣象要素也有非氣象要素。本文主要選取濕度(相對濕度、降水量和水汽壓)、熱力(空氣溫度、日照時(shí)數(shù)和氣溫日較差)和動(dòng)力(風(fēng)速)3類氣象因子來分析對蒸發(fā)量的影響。

    表6給出1964~2013年青海省4個(gè)生態(tài)功能區(qū)上述各氣象因子的氣候傾向率及其與時(shí)間、蒸發(fā)量的相關(guān)系數(shù)和完全相關(guān)系數(shù)。可知:近50 a來,柴達(dá)木盆地日照時(shí)數(shù)、風(fēng)速和氣溫日較差有顯著減小趨勢,平均溫度、水汽壓和降水量為顯著上升趨勢,而相對濕度變化不明顯,且除相對濕度外,其他氣象因子與蒸發(fā)量的相關(guān)性較強(qiáng),致使平均風(fēng)速與蒸發(fā)量的完全相關(guān)系數(shù)最大,其次為日照時(shí)數(shù)、平均氣溫和氣溫日較差,水汽壓和降水量的完全相關(guān)系數(shù)較小,由此可知,影響柴達(dá)木盆地蒸發(fā)量下降的主要因素是風(fēng)速、日照時(shí)數(shù)、平均氣溫和氣溫日較差。環(huán)青海湖區(qū)平均氣溫顯著上升,風(fēng)速和氣溫日較差下降趨勢明顯,其他氣象因子變化趨勢不明顯,且除水汽壓、日照時(shí)數(shù)外,其他氣象因子均與蒸發(fā)量顯著相關(guān),致使平均氣溫與蒸發(fā)量的完全相關(guān)系數(shù)最大,其次為平均風(fēng)速,氣溫日較差的完全相關(guān)系數(shù)較小,可見平均氣溫、風(fēng)速和氣溫日較差是環(huán)青海湖區(qū)蒸發(fā)量上升的主要因素。東部農(nóng)業(yè)區(qū)除平均氣溫外,其他氣象因子與蒸發(fā)量相關(guān)性顯著,且風(fēng)速、相對濕度和日照時(shí)數(shù)下降趨勢明顯,而平均氣溫上升趨勢明顯,致使平均風(fēng)速與蒸發(fā)量的完全相關(guān)系數(shù)最大,其次為相對濕度,日照時(shí)數(shù)的完全相關(guān)系數(shù)較小,故風(fēng)速和相對濕度是東部農(nóng)業(yè)區(qū)蒸發(fā)量下降的主要因素。三江源區(qū)平均氣溫和水汽壓上升趨勢明顯,相對濕度、風(fēng)速和氣溫日較差下降趨勢明顯,但水汽壓、風(fēng)速、降水量和日照時(shí)數(shù)與蒸發(fā)量相關(guān)性不顯著,致使平均氣溫與蒸發(fā)量的完全相關(guān)系數(shù)最大,其次為相對濕度,氣溫日較差的完全相關(guān)系數(shù)較小,因此,平均氣溫、相對濕度和氣溫日較差是三江源區(qū)蒸發(fā)量上升的主要因素。綜上所述,青海省不同生態(tài)功能區(qū)影響蒸發(fā)量的主要因子不同,這與申紅艷等[17]的研究結(jié)論“近50 a來影響整體青海高原蒸發(fā)量的主要因子為平均氣溫和相對濕度”不一致,而與左洪超等[7]的“中國66%站點(diǎn)蒸發(fā)量下降是由多環(huán)境因子共同作用,且不同地區(qū)有不同的表現(xiàn)”結(jié)果一致。

    圖5 1964~2013年青海省4個(gè)生態(tài)功能區(qū)年蒸發(fā)量M-K突變檢驗(yàn)

    平均氣溫氣溫日較差水汽壓相對濕度降水量風(fēng)速日照時(shí)數(shù)柴達(dá)木盆地氣候傾向率0.5-0.20.1-0.36.0-0.3-41.4與時(shí)間相關(guān)系數(shù)0.87***0.64***0.51***0.270.36*0.83***0.64***與蒸發(fā)量相關(guān)系數(shù)-0.47***0.62***-0.60***-0.16-0.56***0.74***0.68***完全相關(guān)系數(shù)0.410.400.31-0.200.620.43環(huán)青海湖區(qū)氣候傾向率0.4-0.10.0-0.36.7-0.1-11.6與時(shí)間相關(guān)系數(shù)0.80***0.33*0.160.210.210.69***0.24與蒸發(fā)量相關(guān)系數(shù)0.39**0.52***-0.03-0.70***-0.40**0.38**0.16完全相關(guān)系數(shù)0.320.17---0.26-東部農(nóng)業(yè)區(qū)氣候傾向率0.40.00.0-0.71.8-0.1-26.9與時(shí)間相關(guān)系數(shù)0.82***0.110.270.39**0.040.73***0.43**與蒸發(fā)量相關(guān)系數(shù)0.250.46***-0.67***-0.75***-0.38**0.54***0.39**完全相關(guān)系數(shù)---0.29-0.390.17三江源區(qū)氣候傾向率0.39-0.130.05-0.78.6-0.11-4.3與時(shí)間相關(guān)系數(shù)0.76***-0.35**0.40**-0.42**0.26-0.56***-0.08與蒸發(fā)量相關(guān)系數(shù)0.77***-0.37**0.04-0.79***0.170.130.09完全相關(guān)系數(shù)0.590.13-0.33---

    注:-表示完全相關(guān)系數(shù)沒有通過顯著性檢驗(yàn)

    3結(jié)論

    (1)青海省年蒸發(fā)量呈從西北向東南逐漸減小的空間分布特征,春、夏、秋季蒸發(fā)量的分布特征與年蒸發(fā)量相似,而冬季的分布特征為從南向北逐漸減小。

    (2)近50 a來,年蒸發(fā)量的氣候傾向率大體自西北向東南逐漸增加,東南部為蒸發(fā)量的上升區(qū)域,中部、西北部為下降區(qū)域;四季蒸發(fā)量的氣候傾向率分布特征與年蒸發(fā)量相似。就不同生態(tài)功能區(qū)而言,柴達(dá)木盆地和三江源區(qū)年蒸發(fā)量整體分別呈顯著下降和上升趨勢,其氣候傾向率分別為-77.9 mm/10 a和30.1 mm/10 a,而環(huán)青海湖區(qū)和東部農(nóng)業(yè)區(qū)的變化趨勢不明顯。4個(gè)功能區(qū)年蒸發(fā)量多項(xiàng)式擬合曲線均呈現(xiàn)先降后升的特征,但轉(zhuǎn)換時(shí)間不完全一致。柴達(dá)木盆地春、夏、秋季蒸發(fā)量均呈顯著下降趨勢,冬季呈顯著上升趨勢;三江源區(qū)蒸發(fā)量除春季外均呈顯著上升趨勢;環(huán)青海湖區(qū)秋季蒸發(fā)量呈顯著上升趨勢,其他季節(jié)變化趨勢不明顯;東部農(nóng)業(yè)區(qū)四季蒸發(fā)量變化趨勢均不明顯。

    (3)青海省4個(gè)生態(tài)功能區(qū)蒸發(fā)量的月和季節(jié)變化特征明顯。月和四季蒸發(fā)量由大至小依次為柴達(dá)木盆地、東部農(nóng)業(yè)區(qū)、環(huán)青海湖區(qū)和三江源區(qū)。柴達(dá)木盆地、環(huán)青海湖區(qū)和東部農(nóng)業(yè)區(qū)蒸發(fā)量的月變化表現(xiàn)為單峰型分布,而三江源區(qū)表現(xiàn)為弱雙峰型分布。夏季蒸發(fā)量最大,其次為春季和秋季,冬季蒸發(fā)量最少。

    (4)柴達(dá)木盆地年蒸發(fā)量在1990年代末減小趨勢明顯,突變時(shí)間是1998年;東部農(nóng)業(yè)區(qū)自1964年一直處于波動(dòng)下降趨勢,2004年以后下降趨勢明顯,突變時(shí)間為2002年;環(huán)青海湖區(qū)和三江源區(qū)突變不明顯。

    (5)影響青海省4個(gè)不同生態(tài)功能區(qū)蒸發(fā)量的主要因子不同:柴達(dá)木盆地為風(fēng)速、日照時(shí)數(shù)、平均氣溫和氣溫日較差;環(huán)青海湖區(qū)為平均氣溫、風(fēng)速和氣溫日較差;東部農(nóng)業(yè)區(qū)為風(fēng)速和相對濕度;平均氣溫、相對濕度和氣溫日較差是三江源區(qū)蒸發(fā)量上升的主要因素。

    本文對青海省4個(gè)生態(tài)功能區(qū)蒸發(fā)量變化特征及其影響氣象因素做了初步分析,進(jìn)一步的物理成因還需探討。由于影響蒸發(fā)量變化的因素眾多,既有氣象要素也有非氣象要素,物理機(jī)制復(fù)雜,它是多個(gè)環(huán)境因子共同非線性相互作用的結(jié)果,將蒸發(fā)量的變化歸結(jié)為任何單一環(huán)境因子的變化所致都會(huì)產(chǎn)生較大的偏頗。

    參考文獻(xiàn):

    [1] Houghton J T, Ding Y, Griggs D J, et al. IPCC, 2001 :Climate Change 2001: The Scientific Basis[M]. Contribution of Working Group Ι to the Third Assessment Report of the Intergovermental Panel on Climate Change. New York:Cambridge University Press,2001,881.

    [2] 江濤,陳永勤,陳俊合,等. 未來氣候變化對我國水文水資源影響的研究[J]. 中山大學(xué)學(xué)報(bào)(自然科學(xué)版),2000,(增刊3):151-157.

    [3] Liu X D, Chen B D. Climatic warming in the Tibetan Plateau during recent decades[J]. International Journal of Climatology, 2000,20(14):1729-1742.

    [4] 施能,陳家其,屠其璞. 中國近100年來4個(gè)年代際的氣候變化特征[J]. 氣象學(xué)報(bào),1995,53(4):431-439.

    [5] Peterson T C, Golubev V S, Groisman P Y. Evaporation losing its strength[J]. Nature, 1995,377:687-688.

    [6] 申雙和,盛瓊. 45年來中國蒸發(fā)皿蒸發(fā)量的變化特征及其成因[J]. 氣象學(xué)報(bào),2008,66(3):452-460.

    [7] 左洪超,李棟梁,胡隱樵,等. 近40 a中國氣候變化趨勢及其同蒸發(fā)皿觀測的蒸發(fā)量變化的關(guān)系[J]. 科學(xué)通報(bào),2005,50(11):1125-1130.

    [8] Liu X M, Zheng H X, Zhang M H, et al. Identification of dominant climate factor for pan evaporation trend in the Tibetan Plateau[J]. Journal of Geographical Sciences. 2011,21(4):594-608.

    [9] Zhang Y Q, Liu C M, Tang Y H, et al.Trends in pan evaporation and reference and actual evapotranspiration across the Tibetan Plateau[J]. Journal of Geographysical Research, 2006, 112, D12110, doi:10.1029/2006 JD008161.

    [10] 劉波,馬柱國,丁裕國,等. 中國北方近45年蒸發(fā)變化的特征及與環(huán)境的關(guān)系[J]. 高原氣象,2006,25(5):840-848.

    [11] 王艷君,姜彤,許崇育. 長江流域蒸發(fā)皿蒸發(fā)量及影響因素變化趨勢[J]. 自然資源學(xué)報(bào),2005,20(6):864-870.

    [12] 張艷玲,王海艷,熊坤,等. 商丘近40 a蒸發(fā)量變化特征及原因分析[J]. 氣象與環(huán)境科學(xué),2007,30(2):72-75.

    [13] 苗運(yùn)玲,卓世新,楊艷玲,等. 新疆哈密市近50 a蒸發(fā)量變化特征及影響因子[J]. 干旱氣象,2013,31(1):95-106.

    [14] 李紅梅,李作偉,王振宇,等. 青藏高原不同生態(tài)功能區(qū)氣候突變時(shí)間的比較分析[J]. 冰川凍土,2012,34(6):1388-1393.

    [15] 李林,朱西德,周陸生,等. 三江源地區(qū)氣候變化及其對生態(tài)環(huán)境的影響[J]. 氣象,2004,30(8):18-22.

    [16] 劉蓓. 43年來青海省蒸發(fā)皿蒸發(fā)量變化及其影響因子[J]. 干旱區(qū)研究,2010,27(6):892-897.

    [17] 申紅艷,馬明亮,汪青春,等. 1961~2010年青海高原蒸發(fā)皿蒸發(fā)量變化及其對水資源的影響[J]. 氣象與環(huán)境學(xué)報(bào),2013,29(6):87-94.

    [18] 李景鑫,王式功,李艷,等. 西寧市蒸發(fā)量變化特征及影響因素[J]. 干旱氣象,2013,31(3):497-503.

    [19] 魏鳳英. 現(xiàn)代氣候統(tǒng)計(jì)診斷與預(yù)測技術(shù)(第2版)[M]. 北京:氣象出版社,2007.36-69.

    [20] 韓軍彩,張秉祥,高祺,等. 石家莊市蒸發(fā)皿蒸發(fā)量的變化特征及其影響因子分析[J]. 干旱氣象,2009,27(4):340-345.

    Variation Characteristics of Pan Evaporation in Different Ecological Function Areas of Qinghai Province During 1964-2013

    QI Donglin1,2,XIAO Hongbin1,2,LI Xiaodong1,2,3,ZHOU Wanfu1,2,SU Wenjiang1,2,HU Aijun1,2

    (1.InstituteofQinghaiMeteorologicalScienceResearch,Xi’ning810001,China;2.QinghaiKeyLaboratoryofDisasterPreventingandReducing,Xi’ning810001,China;3.InstituteofChenduPlateauMeteorologyScienceResearch,CMA,Chengdu610072,China)

    Abstract:Based on the conventional meteorological data of 43 weather stations in different ecological function areas of Qinghai Province during 1964-2013, the temporal and spatial variation characteristics of pan evaporation were analyzed in the past 50 years, and their change causes of pan evaporation in four ecological function areas were analyzed by using the full correlation analysis. The results showed that the annual and seasonal changes of pan evaporation were significant in different ecological function areas of Qinghai Province. In the recent 50 years, the pan evaporation in Qaidam Basin and three river source area showed a significant decreasing and increasing trends, respectively, while that in Qinghai Lake area and eastern agriculture area hadn’t obvious change trend. The annual evaporation in four ecological function areas decreased firstly and then increased during 1964-2013, but the change time was not completely consistent. The spatial distributions of annual pan evaporation decreased from northwest to southeast, and that in spring, summer, fall was similar to the annual, but that in winter decreased gradually from south to north. The annual pan evaporation mutated in 1998 in Qaidam Basin and in 2002 in eastern agriculture area, while that didn’t mutate in three river source area and Qinghai Lake area during 1964-2013. The influences of climate factors on evaporation were different in four ecological function areas of Qinghai Province.The main climate factors affecting evaporation in Qaidam Basin were wind speed, sunshine hours, average temperature and daily range of temperature, those were average temperature, wind speed and daily range of temperature in Qinghai Lake area, and in eastern agriculture area those were wind speed and relative humidity. However, the effects of average temperature and relative humidity on the increasing of pan evaporation in three river source area were critical.

    Key words:pan evaporation; climate mutation; different ecological function areas; Qinghai Province

    中圖分類號:P332

    文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號:1006-7639(2016)-02-0234-09

    doi:10.11755/j.issn.1006-7639(2016)-02-0234

    作者簡介:祁棟林(1967- ),男,高級工程師,主要從事氣候變化和大氣成分研究. E-mail:qidl007@163.com

    基金項(xiàng)目:國家自然科學(xué)基金地區(qū)科學(xué)基金項(xiàng)目(41161009)、公益性行業(yè)(氣象)專項(xiàng)(GYHY201306054、GYHY201506001)和中國氣象局成都高原氣象研究所基金項(xiàng)目(LPM2014004)共同資助

    收稿日期:2015-05-26;改回日期:2015-08-01

    祁棟林,肖宏斌,李曉東,等.1964~2013年青海省不同生態(tài)功能區(qū)蒸發(fā)皿蒸發(fā)量的變化特征[J].干旱氣象,2016,34(2):234-242, [QI Donglin, XIAO Hongbin, LI Xiaodong, et al. Variation Characteristics of Pan Evaporation in Different Ecological Function Areas of Qinghai Province During 1964-2013[J]. Journal of Arid Meteorology, 2016, 34(2):234-242], doi:10.11755/j.issn.1006-7639(2016)-02-0234

    猜你喜歡
    青海省
    青海省高速公路集中監(jiān)控運(yùn)營管理模式探索與實(shí)踐
    人物志
    中國公路(2020年19期)2020-11-09 06:02:06
    2017年青海省一般公共預(yù)算收支決算總表
    青海省非物質(zhì)文化遺產(chǎn)保護(hù)辦法
    青海省國土綠化提速三年行動(dòng)計(jì)劃
    落實(shí)“四個(gè)扎扎實(shí)實(shí)”展現(xiàn)開行使命擔(dān)當(dāng)——國家開發(fā)銀行青海省分行
    青海省交通運(yùn)輸行業(yè)數(shù)據(jù)中心節(jié)能探索
    2015年青海省一般公共預(yù)算收支決算總表
    青海省人民政府辦公廳關(guān)于轉(zhuǎn)發(fā)省科協(xié)青海省自然科學(xué)優(yōu)犁匆圈獨(dú)砰審和勤加祛的通知
    青海省人民政府關(guān)于第二屆青海省質(zhì)量獎(jiǎng)的授獎(jiǎng)決定
    精品无人区乱码1区二区| 精品国内亚洲2022精品成人| 欧美又色又爽又黄视频| 别揉我奶头 嗯啊视频| 琪琪午夜伦伦电影理论片6080| .国产精品久久| 成人无遮挡网站| 在线a可以看的网站| 精品人妻1区二区| 天堂网av新在线| 可以在线观看的亚洲视频| 国产91精品成人一区二区三区| 精品人妻偷拍中文字幕| 真人做人爱边吃奶动态| 国产精品爽爽va在线观看网站| av在线蜜桃| 欧美一级a爱片免费观看看| 无遮挡黄片免费观看| 大又大粗又爽又黄少妇毛片口| 日本五十路高清| 女人被狂操c到高潮| 琪琪午夜伦伦电影理论片6080| 成人国产综合亚洲| 欧美极品一区二区三区四区| 国模一区二区三区四区视频| 国产av在哪里看| 欧美日本视频| 老司机午夜福利在线观看视频| 少妇裸体淫交视频免费看高清| 久久精品夜夜夜夜夜久久蜜豆| 国产伦在线观看视频一区| 国产单亲对白刺激| 能在线免费观看的黄片| 88av欧美| 最近视频中文字幕2019在线8| 窝窝影院91人妻| 国产麻豆成人av免费视频| 91在线精品国自产拍蜜月| 久久婷婷人人爽人人干人人爱| 国产高清视频在线播放一区| 亚洲无线观看免费| 直男gayav资源| 精品人妻熟女av久视频| 欧美成人a在线观看| 免费高清视频大片| 97碰自拍视频| 久久久久久久久久成人| 亚洲国产精品久久男人天堂| 嫩草影院入口| 亚洲一级一片aⅴ在线观看| 三级国产精品欧美在线观看| 午夜免费男女啪啪视频观看 | 中文字幕人妻熟人妻熟丝袜美| 国产精品自产拍在线观看55亚洲| 亚洲国产精品sss在线观看| 日本熟妇午夜| 女人被狂操c到高潮| 日韩av在线大香蕉| av国产免费在线观看| 又黄又爽又免费观看的视频| 高清在线国产一区| 在线观看美女被高潮喷水网站| 久久久午夜欧美精品| 亚洲不卡免费看| 亚洲熟妇熟女久久| 日韩一本色道免费dvd| 中国美白少妇内射xxxbb| 成人美女网站在线观看视频| 97热精品久久久久久| 国产精品av视频在线免费观看| 亚洲欧美日韩东京热| 男人狂女人下面高潮的视频| 狠狠狠狠99中文字幕| 国产成年人精品一区二区| x7x7x7水蜜桃| av视频在线观看入口| 日日撸夜夜添| 欧美又色又爽又黄视频| 国产精品一区二区性色av| 又黄又爽又免费观看的视频| 九九爱精品视频在线观看| 国产精品一区二区三区四区久久| 亚洲专区中文字幕在线| 特级一级黄色大片| 欧美中文日本在线观看视频| 99久久无色码亚洲精品果冻| 99热只有精品国产| 国产精品免费一区二区三区在线| 亚洲精品日韩av片在线观看| 成人国产综合亚洲| 可以在线观看毛片的网站| 免费观看人在逋| 欧美+亚洲+日韩+国产| 亚洲精华国产精华液的使用体验 | 成年女人毛片免费观看观看9| a在线观看视频网站| 国产伦人伦偷精品视频| 亚洲精品色激情综合| 亚洲美女黄片视频| 精品福利观看| 成人欧美大片| 日韩 亚洲 欧美在线| 五月玫瑰六月丁香| 午夜影院日韩av| 国产精品1区2区在线观看.| 国产av麻豆久久久久久久| 久久精品国产99精品国产亚洲性色| 日韩精品中文字幕看吧| 午夜福利成人在线免费观看| 俺也久久电影网| 九色成人免费人妻av| 婷婷六月久久综合丁香| 99热精品在线国产| 国产人妻一区二区三区在| 久久亚洲真实| 制服丝袜大香蕉在线| 久久久久久久久大av| 欧美一区二区精品小视频在线| 国产av麻豆久久久久久久| 最新在线观看一区二区三区| 欧美成人a在线观看| 丰满的人妻完整版| 在线天堂最新版资源| 真实男女啪啪啪动态图| 精品久久久久久久久久久久久| 精品久久久久久成人av| 午夜免费激情av| 黄色视频,在线免费观看| 亚洲真实伦在线观看| 精品久久久噜噜| 亚洲国产欧美人成| 国产伦在线观看视频一区| 亚洲美女搞黄在线观看 | 日本熟妇午夜| 国产综合懂色| 热99re8久久精品国产| av天堂在线播放| 少妇高潮的动态图| 久久久久免费精品人妻一区二区| 成年免费大片在线观看| 午夜免费激情av| 俄罗斯特黄特色一大片| 欧美bdsm另类| 久久精品人妻少妇| 三级男女做爰猛烈吃奶摸视频| 成人综合一区亚洲| 欧美成人免费av一区二区三区| 成年人黄色毛片网站| 亚洲自偷自拍三级| 午夜免费成人在线视频| 亚洲一区二区三区色噜噜| 欧美不卡视频在线免费观看| 久久久久久久久久成人| 亚洲午夜理论影院| 国产精品永久免费网站| x7x7x7水蜜桃| 久9热在线精品视频| 亚洲最大成人手机在线| 毛片一级片免费看久久久久 | 久久精品夜夜夜夜夜久久蜜豆| 亚洲成人久久性| 国产精品国产高清国产av| 91在线观看av| 99久久成人亚洲精品观看| 最新在线观看一区二区三区| 久久久国产成人免费| 免费观看的影片在线观看| 国产成人福利小说| 成人国产麻豆网| 人人妻,人人澡人人爽秒播| 国产精品久久久久久久电影| 国产亚洲精品av在线| 午夜久久久久精精品| 色av中文字幕| 狂野欧美白嫩少妇大欣赏| 禁无遮挡网站| 亚洲欧美日韩卡通动漫| 一级毛片久久久久久久久女| 亚洲av一区综合| 国产一区二区在线av高清观看| 日韩欧美一区二区三区在线观看| 午夜福利欧美成人| 最好的美女福利视频网| 日本免费一区二区三区高清不卡| 婷婷精品国产亚洲av| 亚洲成人免费电影在线观看| 婷婷亚洲欧美| 欧美+日韩+精品| 国产精品野战在线观看| 男插女下体视频免费在线播放| 成年女人永久免费观看视频| 日韩强制内射视频| 成年女人永久免费观看视频| 色哟哟·www| 变态另类丝袜制服| 99久久中文字幕三级久久日本| 99热只有精品国产| 国产精品爽爽va在线观看网站| 亚洲av成人av| 最新在线观看一区二区三区| 国产真实乱freesex| 国内精品宾馆在线| 99riav亚洲国产免费| 一本久久中文字幕| 嫩草影院入口| 男女视频在线观看网站免费| 日本a在线网址| 亚洲狠狠婷婷综合久久图片| 久久亚洲精品不卡| 国内精品美女久久久久久| 全区人妻精品视频| 亚洲欧美日韩东京热| 国产精品乱码一区二三区的特点| 国产精品99久久久久久久久| 欧美日韩国产亚洲二区| 亚洲午夜理论影院| 成人精品一区二区免费| 99精品在免费线老司机午夜| 午夜免费男女啪啪视频观看 | 久久久久久久久久黄片| 国产亚洲91精品色在线| 国产一区二区在线av高清观看| 免费不卡的大黄色大毛片视频在线观看 | 天堂av国产一区二区熟女人妻| 色av中文字幕| 亚洲不卡免费看| 国产久久久一区二区三区| 亚洲av成人精品一区久久| 国产在线男女| 黄色日韩在线| 国产高潮美女av| 免费看a级黄色片| 美女大奶头视频| 97碰自拍视频| eeuss影院久久| 99热网站在线观看| 亚洲av美国av| 免费高清视频大片| 亚洲综合色惰| 97碰自拍视频| 精华霜和精华液先用哪个| 色噜噜av男人的天堂激情| 国产精品av视频在线免费观看| 美女黄网站色视频| 成人国产麻豆网| 床上黄色一级片| 毛片一级片免费看久久久久 | 69av精品久久久久久| 性插视频无遮挡在线免费观看| 欧美成人a在线观看| 国产亚洲av嫩草精品影院| 精品人妻熟女av久视频| 国产精品美女特级片免费视频播放器| 久久精品国产亚洲网站| 一级av片app| 欧洲精品卡2卡3卡4卡5卡区| 中文在线观看免费www的网站| 人人妻人人看人人澡| 免费av不卡在线播放| 国产精品久久久久久久电影| 免费在线观看影片大全网站| 日韩欧美国产一区二区入口| 久久国产乱子免费精品| 色综合亚洲欧美另类图片| 国产精品无大码| 深夜a级毛片| 高清日韩中文字幕在线| 少妇的逼好多水| 狂野欧美激情性xxxx在线观看| 国产高清有码在线观看视频| 99riav亚洲国产免费| 欧美激情在线99| 午夜福利视频1000在线观看| 久久精品国产亚洲网站| 成熟少妇高潮喷水视频| 国产一区二区三区在线臀色熟女| 春色校园在线视频观看| 久久精品国产99精品国产亚洲性色| 久久精品国产鲁丝片午夜精品 | 少妇高潮的动态图| 黄色日韩在线| av在线观看视频网站免费| 美女黄网站色视频| 黄色视频,在线免费观看| 午夜久久久久精精品| 亚洲黑人精品在线| 欧美色视频一区免费| 国产亚洲欧美98| 人妻久久中文字幕网| 欧美一级a爱片免费观看看| 中文资源天堂在线| 国产精品久久久久久精品电影| 欧美性猛交黑人性爽| 亚洲人与动物交配视频| 亚洲久久久久久中文字幕| 欧美日韩亚洲国产一区二区在线观看| 看片在线看免费视频| 动漫黄色视频在线观看| 国产精品免费一区二区三区在线| 国产人妻一区二区三区在| 中国美白少妇内射xxxbb| 国产精品免费一区二区三区在线| 91精品国产九色| 三级国产精品欧美在线观看| 日韩亚洲欧美综合| av福利片在线观看| 午夜亚洲福利在线播放| 国产精品不卡视频一区二区| 深夜精品福利| 日韩欧美在线二视频| 国产精品无大码| 亚洲欧美日韩高清专用| 波野结衣二区三区在线| 日本免费a在线| 99久久成人亚洲精品观看| 亚洲精品一区av在线观看| 男插女下体视频免费在线播放| 国产伦一二天堂av在线观看| 97超视频在线观看视频| 尤物成人国产欧美一区二区三区| 国产 一区精品| 最近在线观看免费完整版| 成人永久免费在线观看视频| 精品乱码久久久久久99久播| 中文字幕免费在线视频6| 22中文网久久字幕| 午夜老司机福利剧场| 国产在线男女| 久久精品影院6| 一卡2卡三卡四卡精品乱码亚洲| 日韩中文字幕欧美一区二区| 亚洲四区av| 又黄又爽又免费观看的视频| 我要看日韩黄色一级片| 精品一区二区三区av网在线观看| 97超视频在线观看视频| 欧美高清性xxxxhd video| 人妻夜夜爽99麻豆av| 欧美日韩乱码在线| 中文资源天堂在线| 免费不卡的大黄色大毛片视频在线观看 | 成人av在线播放网站| 永久网站在线| 色播亚洲综合网| 国产激情偷乱视频一区二区| 91麻豆精品激情在线观看国产| 亚洲一区高清亚洲精品| 国产极品精品免费视频能看的| 天美传媒精品一区二区| 亚洲美女搞黄在线观看 | av在线蜜桃| 亚洲国产精品成人综合色| 少妇的逼水好多| 久久久久久伊人网av| 国产精品伦人一区二区| 欧美不卡视频在线免费观看| 色精品久久人妻99蜜桃| 欧美精品啪啪一区二区三区| or卡值多少钱| 天堂网av新在线| 99热只有精品国产| 免费黄网站久久成人精品| av中文乱码字幕在线| 久久午夜亚洲精品久久| 日韩中文字幕欧美一区二区| 欧美性猛交╳xxx乱大交人| 午夜日韩欧美国产| 亚洲av一区综合| 国产精品99久久久久久久久| 欧美激情在线99| 欧美不卡视频在线免费观看| 在线观看午夜福利视频| 国产午夜福利久久久久久| 给我免费播放毛片高清在线观看| 国产一区二区三区视频了| 午夜福利视频1000在线观看| 两个人的视频大全免费| 能在线免费观看的黄片| 日本五十路高清| 日本黄色片子视频| 黄片wwwwww| 特大巨黑吊av在线直播| 国产免费av片在线观看野外av| 蜜桃久久精品国产亚洲av| 国产精品,欧美在线| 精品99又大又爽又粗少妇毛片 | 成年女人毛片免费观看观看9| 深夜a级毛片| 女人被狂操c到高潮| 嫩草影视91久久| 人妻少妇偷人精品九色| 久久香蕉精品热| 国产精品精品国产色婷婷| 亚洲avbb在线观看| 久久热精品热| 热99re8久久精品国产| 日本黄大片高清| 亚洲七黄色美女视频| 亚洲性久久影院| 精品久久久久久久久av| 日韩大尺度精品在线看网址| 日本精品一区二区三区蜜桃| 久久国产精品人妻蜜桃| 国产精品人妻久久久影院| 亚洲第一电影网av| 国产主播在线观看一区二区| 日本在线视频免费播放| 少妇丰满av| 亚洲五月天丁香| 久久香蕉精品热| 免费搜索国产男女视频| 老司机深夜福利视频在线观看| 久久婷婷人人爽人人干人人爱| 欧美另类亚洲清纯唯美| 亚洲国产精品sss在线观看| 午夜爱爱视频在线播放| 午夜日韩欧美国产| 成人特级av手机在线观看| 亚洲三级黄色毛片| 国产蜜桃级精品一区二区三区| 乱码一卡2卡4卡精品| 在线观看一区二区三区| 欧美另类亚洲清纯唯美| 久久久成人免费电影| 精品久久久久久成人av| 91av网一区二区| 能在线免费观看的黄片| 欧美+亚洲+日韩+国产| 国产淫片久久久久久久久| 国产高潮美女av| 亚洲最大成人手机在线| 永久网站在线| 无人区码免费观看不卡| 精品人妻1区二区| 三级毛片av免费| 精品久久久久久,| 99热这里只有是精品在线观看| 亚洲最大成人手机在线| 国产精品一区二区三区四区免费观看 | 国产男靠女视频免费网站| 国产成人av教育| 国产高清不卡午夜福利| 午夜老司机福利剧场| 波野结衣二区三区在线| 色综合站精品国产| 久久久久免费精品人妻一区二区| 国产精品爽爽va在线观看网站| 亚洲在线自拍视频| 久久久精品欧美日韩精品| 免费看光身美女| 亚洲精品久久国产高清桃花| 成人av一区二区三区在线看| www.www免费av| 91久久精品电影网| 99在线视频只有这里精品首页| 欧美+亚洲+日韩+国产| bbb黄色大片| 成人av一区二区三区在线看| 国产中年淑女户外野战色| 亚洲精品亚洲一区二区| 韩国av在线不卡| 91av网一区二区| 一级av片app| 国产极品精品免费视频能看的| 无人区码免费观看不卡| 久久久午夜欧美精品| 成人亚洲精品av一区二区| eeuss影院久久| 中文字幕人妻熟人妻熟丝袜美| 国国产精品蜜臀av免费| 男女那种视频在线观看| 女生性感内裤真人,穿戴方法视频| 成人一区二区视频在线观看| 床上黄色一级片| 国产精品人妻久久久久久| 日韩欧美免费精品| 日韩,欧美,国产一区二区三区 | 日韩一区二区视频免费看| 香蕉av资源在线| 联通29元200g的流量卡| 国产一区二区亚洲精品在线观看| 在线国产一区二区在线| 日日撸夜夜添| 大又大粗又爽又黄少妇毛片口| АⅤ资源中文在线天堂| 一本精品99久久精品77| 国国产精品蜜臀av免费| 国产精品不卡视频一区二区| 欧美性猛交╳xxx乱大交人| 亚洲经典国产精华液单| 99riav亚洲国产免费| 伊人久久精品亚洲午夜| 久久午夜福利片| 欧美区成人在线视频| 午夜老司机福利剧场| 在线观看av片永久免费下载| 韩国av在线不卡| 国产乱人视频| 特大巨黑吊av在线直播| 婷婷六月久久综合丁香| 国产高清三级在线| 亚洲熟妇熟女久久| 人人妻人人看人人澡| 有码 亚洲区| 在线国产一区二区在线| 日韩欧美免费精品| 在线观看免费视频日本深夜| 精品欧美国产一区二区三| 色哟哟哟哟哟哟| 亚洲成a人片在线一区二区| 深夜a级毛片| 校园春色视频在线观看| 国产精品一及| 成人一区二区视频在线观看| 日韩强制内射视频| 白带黄色成豆腐渣| 日韩欧美精品免费久久| 国产老妇女一区| 狂野欧美白嫩少妇大欣赏| 精品人妻1区二区| 色视频www国产| 国产精品一及| 国产精品电影一区二区三区| 国产一区二区在线av高清观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 欧美精品啪啪一区二区三区| 久9热在线精品视频| 亚洲精品亚洲一区二区| 又爽又黄无遮挡网站| 日韩欧美一区二区三区在线观看| 啦啦啦观看免费观看视频高清| 国产一区二区在线av高清观看| 亚洲最大成人中文| 午夜福利在线观看吧| 久久精品91蜜桃| 国产亚洲欧美98| 日本 av在线| 又爽又黄a免费视频| 久久这里只有精品中国| 午夜精品一区二区三区免费看| 色哟哟·www| 国产毛片a区久久久久| 日本一二三区视频观看| 国产精品1区2区在线观看.| 一区二区三区免费毛片| 两个人的视频大全免费| 日本黄色片子视频| eeuss影院久久| 日日摸夜夜添夜夜添小说| 国产精品精品国产色婷婷| 精品不卡国产一区二区三区| 久99久视频精品免费| av天堂在线播放| 久久国产精品人妻蜜桃| 国产一区二区亚洲精品在线观看| 成年版毛片免费区| 午夜免费激情av| 久久热精品热| 内地一区二区视频在线| 亚洲18禁久久av| 久久久精品大字幕| 一级毛片久久久久久久久女| h日本视频在线播放| 亚洲精品影视一区二区三区av| 欧美又色又爽又黄视频| 欧美极品一区二区三区四区| 久久精品夜夜夜夜夜久久蜜豆| 国产高清三级在线| 久久久久久大精品| 国产av麻豆久久久久久久| 亚洲五月天丁香| 色哟哟哟哟哟哟| 狂野欧美白嫩少妇大欣赏| 午夜a级毛片| 在线看三级毛片| 免费电影在线观看免费观看| 看片在线看免费视频| 国产伦人伦偷精品视频| 国产一级毛片七仙女欲春2| 国产亚洲欧美98| 男女做爰动态图高潮gif福利片| 国产精品电影一区二区三区| 欧美性感艳星| 国产精品一区二区三区四区免费观看 | 亚洲国产欧美人成| 日韩国内少妇激情av| 中亚洲国语对白在线视频| 婷婷精品国产亚洲av在线| 亚洲avbb在线观看| a级毛片免费高清观看在线播放| 亚洲最大成人中文| 午夜免费激情av| 亚洲精品日韩av片在线观看| 乱码一卡2卡4卡精品| 久久久久久伊人网av| 国产探花在线观看一区二区| 一夜夜www| 亚洲成人久久性| 婷婷丁香在线五月| 热99在线观看视频| 国产伦一二天堂av在线观看| 一边摸一边抽搐一进一小说| 国产黄片美女视频| 一本一本综合久久| 亚洲乱码一区二区免费版| 又爽又黄无遮挡网站| 亚洲av美国av| 尾随美女入室| 1024手机看黄色片| 亚洲精品日韩av片在线观看| 又爽又黄a免费视频| 69人妻影院| 久久精品91蜜桃| 欧美极品一区二区三区四区|