• <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)決定
    成年女人毛片免费观看观看9| www.自偷自拍.com| av欧美777| 国产精品精品国产色婷婷| 精品欧美国产一区二区三| 国产亚洲欧美精品永久| 久久精品亚洲精品国产色婷小说| 日韩成人在线观看一区二区三区| 国产高清有码在线观看视频 | 欧美在线一区亚洲| 日韩一卡2卡3卡4卡2021年| 丰满人妻熟妇乱又伦精品不卡| 男男h啪啪无遮挡| 不卡一级毛片| 国产av不卡久久| 亚洲成av人片免费观看| 久久午夜综合久久蜜桃| 欧美日韩精品网址| 中文字幕人成人乱码亚洲影| 久久久精品欧美日韩精品| 99精品久久久久人妻精品| 久久人妻av系列| 一级黄色大片毛片| 日韩 欧美 亚洲 中文字幕| 国产精品av久久久久免费| 亚洲人成电影免费在线| 一级毛片女人18水好多| 午夜两性在线视频| 国产真人三级小视频在线观看| avwww免费| 午夜两性在线视频| 欧美乱色亚洲激情| 免费在线观看视频国产中文字幕亚洲| 精品国产国语对白av| 国产久久久一区二区三区| 黄色丝袜av网址大全| 99久久精品国产亚洲精品| 国产精华一区二区三区| 久久精品91无色码中文字幕| 欧美日韩乱码在线| 国产黄a三级三级三级人| 99re在线观看精品视频| av中文乱码字幕在线| 国产精品免费视频内射| 亚洲午夜精品一区,二区,三区| 美国免费a级毛片| 免费在线观看影片大全网站| 一二三四社区在线视频社区8| 男人的好看免费观看在线视频 | 成人国产综合亚洲| 美女午夜性视频免费| or卡值多少钱| 午夜精品在线福利| 精品第一国产精品| 亚洲中文字幕一区二区三区有码在线看 | 色尼玛亚洲综合影院| 亚洲av五月六月丁香网| 久久狼人影院| 亚洲黑人精品在线| svipshipincom国产片| 免费无遮挡裸体视频| 午夜福利欧美成人| 99久久无色码亚洲精品果冻| 久久人妻av系列| 精品国产亚洲在线| 亚洲av成人av| 后天国语完整版免费观看| 午夜免费鲁丝| 欧美日本亚洲视频在线播放| 色播亚洲综合网| 午夜成年电影在线免费观看| 日韩大码丰满熟妇| 国产v大片淫在线免费观看| 老司机深夜福利视频在线观看| 亚洲欧美日韩无卡精品| 男女之事视频高清在线观看| 久久久久国产一级毛片高清牌| 欧美精品啪啪一区二区三区| 国产私拍福利视频在线观看| 一进一出抽搐动态| 人人妻人人看人人澡| 韩国精品一区二区三区| or卡值多少钱| 美女扒开内裤让男人捅视频| 欧美成人午夜精品| 国产精品,欧美在线| 两性夫妻黄色片| 脱女人内裤的视频| 此物有八面人人有两片| 叶爱在线成人免费视频播放| svipshipincom国产片| 亚洲人成网站在线播放欧美日韩| 亚洲熟妇熟女久久| 久久久精品国产亚洲av高清涩受| 国产亚洲精品av在线| 久久国产乱子伦精品免费另类| 精品免费久久久久久久清纯| 丰满人妻熟妇乱又伦精品不卡| 黄网站色视频无遮挡免费观看| 国产精品免费一区二区三区在线| 日韩有码中文字幕| 国产亚洲精品综合一区在线观看 | 91av网站免费观看| 特大巨黑吊av在线直播 | 精品欧美国产一区二区三| 日日爽夜夜爽网站| 午夜福利在线在线| 亚洲av第一区精品v没综合| 成人亚洲精品av一区二区| 亚洲欧美精品综合一区二区三区| 长腿黑丝高跟| 桃色一区二区三区在线观看| 一级黄色大片毛片| 国产av不卡久久| 久久久精品国产亚洲av高清涩受| 777久久人妻少妇嫩草av网站| 亚洲激情在线av| 久久九九热精品免费| 国产成人av激情在线播放| 少妇粗大呻吟视频| 99国产综合亚洲精品| 99久久综合精品五月天人人| 亚洲色图av天堂| 日本 av在线| 一边摸一边抽搐一进一小说| 一级毛片高清免费大全| 亚洲国产精品sss在线观看| 国产伦人伦偷精品视频| 男人操女人黄网站| 免费观看人在逋| 久久精品国产亚洲av高清一级| 国产麻豆成人av免费视频| 国产aⅴ精品一区二区三区波| 国产伦在线观看视频一区| 午夜福利高清视频| 色综合站精品国产| 久久精品成人免费网站| 很黄的视频免费| 黑人操中国人逼视频| 999久久久国产精品视频| 国产一区二区激情短视频| 久99久视频精品免费| 国产1区2区3区精品| 国产一级毛片七仙女欲春2 | 欧美日韩福利视频一区二区| 黄片大片在线免费观看| 在线看三级毛片| 亚洲男人的天堂狠狠| 精品久久久久久久人妻蜜臀av| 日韩中文字幕欧美一区二区| 免费看a级黄色片| 在线免费观看的www视频| 久久久久久久午夜电影| 人成视频在线观看免费观看| 国产成人欧美在线观看| 国产v大片淫在线免费观看| 色在线成人网| 久久久水蜜桃国产精品网| 久9热在线精品视频| 中文字幕最新亚洲高清| 美国免费a级毛片| 老汉色∧v一级毛片| 男女午夜视频在线观看| 欧美激情极品国产一区二区三区| 一本久久中文字幕| 精品久久久久久久久久久久久 | 此物有八面人人有两片| 国产av不卡久久| 国产精华一区二区三区| 黄色丝袜av网址大全| 白带黄色成豆腐渣| 国产成人系列免费观看| 亚洲性夜色夜夜综合| 非洲黑人性xxxx精品又粗又长| 啦啦啦观看免费观看视频高清| 少妇 在线观看| 日韩成人在线观看一区二区三区| 亚洲无线在线观看| 久久香蕉激情| 国产精品一区二区免费欧美| 婷婷精品国产亚洲av| 国内揄拍国产精品人妻在线 | 亚洲av电影不卡..在线观看| 国产熟女xx| 久久欧美精品欧美久久欧美| 国产黄a三级三级三级人| 18美女黄网站色大片免费观看| 免费在线观看黄色视频的| 欧美成人免费av一区二区三区| 精品卡一卡二卡四卡免费| 大型黄色视频在线免费观看| 亚洲精品美女久久久久99蜜臀| 老鸭窝网址在线观看| 色尼玛亚洲综合影院| 欧美日韩瑟瑟在线播放| 免费看a级黄色片| 亚洲精品国产精品久久久不卡| 精品不卡国产一区二区三区| 国产一区在线观看成人免费| 国产亚洲欧美在线一区二区| 日韩高清综合在线| 国产日本99.免费观看| 国产aⅴ精品一区二区三区波| 成人国产综合亚洲| 97人妻精品一区二区三区麻豆 | 亚洲精品中文字幕一二三四区| 麻豆国产av国片精品| 最近最新中文字幕大全电影3 | 国产蜜桃级精品一区二区三区| or卡值多少钱| 91国产中文字幕| 成人18禁在线播放| 国产精品国产高清国产av| 最好的美女福利视频网| 国产精品一区二区三区四区久久 | 日本一本二区三区精品| 在线观看66精品国产| 久久精品aⅴ一区二区三区四区| 制服丝袜大香蕉在线| 国产视频内射| 国产1区2区3区精品| 国产99白浆流出| 亚洲av美国av| 禁无遮挡网站| 精品人妻1区二区| 免费av毛片视频| 亚洲七黄色美女视频| 免费看日本二区| 欧美绝顶高潮抽搐喷水| 免费看a级黄色片| 欧美性猛交黑人性爽| 一进一出抽搐动态| 人人妻人人看人人澡| 人人妻,人人澡人人爽秒播| 视频区欧美日本亚洲| xxxwww97欧美| 男男h啪啪无遮挡| 国产av在哪里看| 免费看美女性在线毛片视频| aaaaa片日本免费| 欧美最黄视频在线播放免费| 精品国产超薄肉色丝袜足j| 免费av毛片视频| 精品久久久久久久末码| 国产成人啪精品午夜网站| 欧美激情久久久久久爽电影| 国产男靠女视频免费网站| 一个人免费在线观看的高清视频| 久久久国产成人免费| 国产激情偷乱视频一区二区| 日本黄色视频三级网站网址| 久久婷婷人人爽人人干人人爱| 精品免费久久久久久久清纯| 久久草成人影院| 久久久久国产精品人妻aⅴ院| 男人舔女人的私密视频| 久久久久久久午夜电影| 国产成人一区二区三区免费视频网站| 中文字幕久久专区| 成人av一区二区三区在线看| 国产av不卡久久| 精品一区二区三区四区五区乱码| 麻豆一二三区av精品| 18禁观看日本| 免费高清在线观看日韩| 日韩有码中文字幕| 麻豆一二三区av精品| 日韩欧美国产一区二区入口| 久久久水蜜桃国产精品网| 一本一本综合久久| 精品一区二区三区av网在线观看| 亚洲精品国产区一区二| 国产亚洲欧美在线一区二区| 午夜老司机福利片| 欧美精品啪啪一区二区三区| 欧美在线黄色| 欧美在线一区亚洲| 国产精品 欧美亚洲| 欧美三级亚洲精品| 亚洲精品美女久久av网站| 在线看三级毛片| 一本综合久久免费| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品自产拍在线观看55亚洲| 久久久久国产一级毛片高清牌| 搡老熟女国产l中国老女人| 好男人电影高清在线观看| 免费女性裸体啪啪无遮挡网站| 女人被狂操c到高潮| 亚洲av片天天在线观看| 97人妻精品一区二区三区麻豆 | aaaaa片日本免费| 欧美日本视频| 国产亚洲精品第一综合不卡| 欧美性猛交黑人性爽| aaaaa片日本免费| 国产1区2区3区精品| 国产成人精品久久二区二区免费| 黑人操中国人逼视频| 国产成人欧美| 欧美三级亚洲精品| 无遮挡黄片免费观看| 美女高潮喷水抽搐中文字幕| 中国美女看黄片| 国产单亲对白刺激| 中文字幕最新亚洲高清| 亚洲午夜精品一区,二区,三区| 中文字幕人成人乱码亚洲影| 国产精品 国内视频| 亚洲片人在线观看| 亚洲成av片中文字幕在线观看| 一级作爱视频免费观看| 免费观看精品视频网站| 午夜激情av网站| 久久久久久免费高清国产稀缺| 亚洲第一青青草原| 可以在线观看的亚洲视频| 欧美 亚洲 国产 日韩一| 美女午夜性视频免费| 亚洲国产中文字幕在线视频| 欧美成人午夜精品| 亚洲av电影在线进入| 可以免费在线观看a视频的电影网站| 黑人巨大精品欧美一区二区mp4| av有码第一页| 色av中文字幕| 夜夜看夜夜爽夜夜摸| 亚洲男人天堂网一区| 亚洲人成伊人成综合网2020| 又黄又爽又免费观看的视频| 99在线人妻在线中文字幕| 性欧美人与动物交配| 亚洲中文字幕日韩| 亚洲熟妇熟女久久| 欧美丝袜亚洲另类 | 国内揄拍国产精品人妻在线 | 亚洲欧美日韩高清在线视频| 伊人久久大香线蕉亚洲五| 久久久久久久精品吃奶| 黑人欧美特级aaaaaa片| 婷婷精品国产亚洲av在线| 99在线人妻在线中文字幕| 日本撒尿小便嘘嘘汇集6| 国产精品98久久久久久宅男小说| АⅤ资源中文在线天堂| 欧美日韩一级在线毛片| 亚洲自偷自拍图片 自拍| 亚洲男人天堂网一区| 成人国语在线视频| 久久久久久免费高清国产稀缺| 制服丝袜大香蕉在线| 中文字幕最新亚洲高清| 一区二区三区高清视频在线| 久久久久久国产a免费观看| 久久这里只有精品19| 日本熟妇午夜| 香蕉av资源在线| 99久久精品国产亚洲精品| 一二三四在线观看免费中文在| 亚洲久久久国产精品| 国产亚洲欧美98| 午夜激情av网站| 成人亚洲精品一区在线观看| www.精华液| 午夜福利高清视频| 夜夜爽天天搞| 精品国产乱子伦一区二区三区| 精品国产一区二区三区四区第35| av福利片在线| 欧美日韩精品网址| 国产精品日韩av在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久免费高清国产稀缺| 搡老妇女老女人老熟妇| 中文字幕久久专区| 国产精品香港三级国产av潘金莲| 亚洲人成电影免费在线| 成年版毛片免费区| 国产一区二区在线av高清观看| 久久精品国产99精品国产亚洲性色| 侵犯人妻中文字幕一二三四区| a级毛片a级免费在线| 丰满人妻熟妇乱又伦精品不卡| 日韩免费av在线播放| 亚洲精华国产精华精| 欧美乱码精品一区二区三区| 国产欧美日韩精品亚洲av| 中出人妻视频一区二区| 成人特级黄色片久久久久久久| 国产精品精品国产色婷婷| 国产av一区在线观看免费| 最好的美女福利视频网| 国产亚洲av嫩草精品影院| 成年版毛片免费区| 亚洲中文av在线| 午夜久久久在线观看| 97碰自拍视频| 男人的好看免费观看在线视频 | 女性被躁到高潮视频| 久久久水蜜桃国产精品网| 99热这里只有精品一区 | 亚洲无线在线观看| 欧美另类亚洲清纯唯美| 999久久久国产精品视频| 黄频高清免费视频| svipshipincom国产片| 国产高清视频在线播放一区| or卡值多少钱| svipshipincom国产片| 日韩精品青青久久久久久| 18美女黄网站色大片免费观看| 99久久99久久久精品蜜桃| 黄片小视频在线播放| e午夜精品久久久久久久| 亚洲五月色婷婷综合| 亚洲九九香蕉| 国产成人av激情在线播放| 日本 av在线| av片东京热男人的天堂| 黄网站色视频无遮挡免费观看| 搡老妇女老女人老熟妇| 黄色毛片三级朝国网站| 18禁黄网站禁片免费观看直播| 国产极品粉嫩免费观看在线| 亚洲一区中文字幕在线| 成人免费观看视频高清| 国产久久久一区二区三区| 精品人妻1区二区| 国产精华一区二区三区| 女性被躁到高潮视频| 亚洲狠狠婷婷综合久久图片| 午夜久久久久精精品| 亚洲美女黄片视频| 免费在线观看黄色视频的| 一级a爱视频在线免费观看| 欧美日韩亚洲综合一区二区三区_| 黑人欧美特级aaaaaa片| av欧美777| 亚洲免费av在线视频| 亚洲国产欧美日韩在线播放| 人成视频在线观看免费观看| 99国产精品99久久久久| 2021天堂中文幕一二区在线观 | 黑人操中国人逼视频| 精品乱码久久久久久99久播| 国产激情久久老熟女| 国产亚洲精品第一综合不卡| 日本一区二区免费在线视频| 亚洲电影在线观看av| 亚洲黑人精品在线| 91九色精品人成在线观看| 十八禁网站免费在线| 丰满人妻熟妇乱又伦精品不卡| 欧美av亚洲av综合av国产av| 操出白浆在线播放| 女同久久另类99精品国产91| netflix在线观看网站| 欧美激情 高清一区二区三区| 国产久久久一区二区三区| 男女下面进入的视频免费午夜 | 亚洲专区字幕在线| 女性生殖器流出的白浆| 50天的宝宝边吃奶边哭怎么回事| 日本 av在线| 亚洲第一欧美日韩一区二区三区| 伊人久久大香线蕉亚洲五| av福利片在线| 国产99白浆流出| 亚洲精品久久国产高清桃花| 免费无遮挡裸体视频| 男女之事视频高清在线观看| 亚洲av片天天在线观看| 黄色视频不卡| 两个人看的免费小视频| 久久99热这里只有精品18| 老司机午夜福利在线观看视频| 嫩草影视91久久| 在线观看免费日韩欧美大片| 黄色成人免费大全| 国产伦人伦偷精品视频| 9191精品国产免费久久| av有码第一页| 高潮久久久久久久久久久不卡| 香蕉av资源在线| 国产成人啪精品午夜网站| 国产真人三级小视频在线观看| 桃色一区二区三区在线观看| 日韩免费av在线播放| 男人操女人黄网站| 亚洲熟女毛片儿| 99国产精品一区二区蜜桃av| 国产成人精品久久二区二区免费| 在线天堂中文资源库| 精品一区二区三区av网在线观看| 精品国产美女av久久久久小说| 成熟少妇高潮喷水视频| 欧美av亚洲av综合av国产av| 99久久精品国产亚洲精品| 国产av在哪里看| 99久久精品国产亚洲精品| 嫩草影视91久久| bbb黄色大片| 国产高清视频在线播放一区| 香蕉国产在线看| 欧美日韩黄片免| 99riav亚洲国产免费| 久久国产精品影院| 欧美成人一区二区免费高清观看 | x7x7x7水蜜桃| 12—13女人毛片做爰片一| 俺也久久电影网| 久热爱精品视频在线9| 国产精品 国内视频| 亚洲午夜精品一区,二区,三区| 久久久久久久精品吃奶| 欧美日韩瑟瑟在线播放| 欧美zozozo另类| 99国产精品一区二区三区| 一夜夜www| 国产精品乱码一区二三区的特点| 日本a在线网址| 色婷婷久久久亚洲欧美| av视频在线观看入口| 美女国产高潮福利片在线看| 亚洲七黄色美女视频| 亚洲午夜精品一区,二区,三区| 最近最新中文字幕大全免费视频| 欧美国产精品va在线观看不卡| 最近最新中文字幕大全免费视频| 岛国在线观看网站| 色播亚洲综合网| 亚洲精品一卡2卡三卡4卡5卡| av超薄肉色丝袜交足视频| 国产精品98久久久久久宅男小说| 色精品久久人妻99蜜桃| 麻豆成人av在线观看| 高清毛片免费观看视频网站| 国产亚洲av高清不卡| 欧美zozozo另类| 国产乱人伦免费视频| 久久精品91无色码中文字幕| 精品不卡国产一区二区三区| a在线观看视频网站| 精品国产超薄肉色丝袜足j| 国产av又大| 亚洲专区中文字幕在线| 18禁黄网站禁片午夜丰满| 国产精品影院久久| 非洲黑人性xxxx精品又粗又长| av中文乱码字幕在线| 成人国语在线视频| 国产真人三级小视频在线观看| 精品日产1卡2卡| 99re在线观看精品视频| 高清在线国产一区| 亚洲av美国av| 日日夜夜操网爽| 亚洲av电影在线进入| 中文字幕精品免费在线观看视频| 51午夜福利影视在线观看| 婷婷亚洲欧美| 久久精品成人免费网站| 亚洲 欧美 日韩 在线 免费| 亚洲色图 男人天堂 中文字幕| 伊人久久大香线蕉亚洲五| 97超级碰碰碰精品色视频在线观看| 国产精品影院久久| 夜夜夜夜夜久久久久| 国产真人三级小视频在线观看| 国内久久婷婷六月综合欲色啪| 日韩欧美一区二区三区在线观看| 亚洲天堂国产精品一区在线| 十八禁网站免费在线| 在线国产一区二区在线| 日韩欧美免费精品| 手机成人av网站| 18美女黄网站色大片免费观看| 亚洲三区欧美一区| 久久香蕉激情| 2021天堂中文幕一二区在线观 | 无遮挡黄片免费观看| 一卡2卡三卡四卡精品乱码亚洲| 午夜福利一区二区在线看| 黄片播放在线免费| 日韩三级视频一区二区三区| 精品卡一卡二卡四卡免费| а√天堂www在线а√下载| 一二三四在线观看免费中文在| 欧美乱妇无乱码| 久久久久亚洲av毛片大全| 精品日产1卡2卡| 国产91精品成人一区二区三区| 一级毛片女人18水好多| 听说在线观看完整版免费高清| 白带黄色成豆腐渣| 午夜免费激情av| 天天躁夜夜躁狠狠躁躁| 99久久精品国产亚洲精品| 男女下面进入的视频免费午夜 | 悠悠久久av| 欧美黄色淫秽网站| 男女做爰动态图高潮gif福利片| 黄片播放在线免费| 国产精品九九99| 麻豆一二三区av精品| 中文字幕另类日韩欧美亚洲嫩草| 亚洲色图av天堂| a在线观看视频网站| 999久久久国产精品视频| 最新在线观看一区二区三区| 久久精品aⅴ一区二区三区四区| 久久久久久久久免费视频了| 女性生殖器流出的白浆| 757午夜福利合集在线观看| 一本久久中文字幕|