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

    秦嶺南北地區(qū)絕對濕度的時空變化及其與潛在蒸發(fā)量的關系

    2015-03-10 08:40:54穆興民馬文勇于新洋劉憲鋒李建國劉思潔
    生態(tài)學報 2015年2期
    關鍵詞:子區(qū)谷地蒸發(fā)量

    蔣 沖,穆興民,馬文勇,于新洋,劉憲鋒,李建國,劉思潔,王 飛,*

    1 西北農(nóng)林科技大學水土保持研究所,楊凌 712100 2 北京師范大學全球變化與地球系統(tǒng)科學研究院地表過程與資源生態(tài)國家重點實驗室,北京 100875 3 中國科學院地理科學與資源研究所,北京 100101 4 北京師范大學資源學院,北京 100875 5 中國科學院寒區(qū)旱區(qū)環(huán)境與工程研究所,蘭州 730000 6 北京大學遙感與地理信息研究所,北京 100871

    秦嶺南北地區(qū)絕對濕度的時空變化及其與潛在蒸發(fā)量的關系

    蔣 沖1, 2,穆興民1,馬文勇3,于新洋3,劉憲鋒4,李建國5,劉思潔6,王 飛1,*

    1 西北農(nóng)林科技大學水土保持研究所,楊凌 712100 2 北京師范大學全球變化與地球系統(tǒng)科學研究院地表過程與資源生態(tài)國家重點實驗室,北京 100875 3 中國科學院地理科學與資源研究所,北京 100101 4 北京師范大學資源學院,北京 100875 5 中國科學院寒區(qū)旱區(qū)環(huán)境與工程研究所,蘭州 730000 6 北京大學遙感與地理信息研究所,北京 100871

    基于秦嶺南北地區(qū)1960—2011年的氣象觀測資料,分別利用絕對濕度計算公式和Penman-Monteith公式估算大氣絕對濕度和潛在蒸發(fā)量,并進一步分析了絕對濕度的空間分布規(guī)律、時空變化特征及其與潛在蒸發(fā)量變化的響應關系。結果表明:(1)秦嶺南北地區(qū)絕對濕度由南向北順次遞減,具有較好的海拔地帶性和緯度地帶性。季節(jié)平均濕度以夏季為最大,冬季最小。(2)近52年絕對濕度除巴巫谷地表現(xiàn)出不顯著下降趨勢外,其它子區(qū)均呈現(xiàn)出增加趨勢。1986年和1998年是濕度變化的轉折點,1960—1986年以微弱下降為主,此后直至1998年震蕩上升,1998年以后呈下降趨勢。(3)年尺度和春、秋兩季,除巴巫谷地外,濕度和潛在蒸發(fā)量均表現(xiàn)出顯著的負相關關系;而在夏季和冬季,除漢水流域和巴巫谷地不顯著正相關以外,其它區(qū)域也均為負相關。年度和春、秋兩季兩個指標負相關的緊密程度隨著區(qū)域的南移而逐漸減弱。1960—2011年間,年度和季節(jié)尺度潛在蒸發(fā)和絕對濕度呈反向變化趨勢;而在1960—1989年間,兩者同向變化;1990—2011年間,年度和春、冬兩季潛在蒸發(fā)上升,而同期的絕對濕度則是先上升再下降。(4)實際蒸發(fā)量的增加導致空氣中水汽含量(濕度)增加,反過來抑制了水面蒸發(fā)(潛在蒸發(fā)量)。秦嶺以北、秦嶺南坡、漢水流域和巴巫谷地部分地區(qū)絕對濕度和潛在蒸發(fā)量為互補關系,由北向南隨著水分限制作用的不斷減弱兩者逐漸轉化為不顯著的正相關關系。

    秦嶺南北;絕對濕度;時空變化;潛在蒸發(fā)量;互補關系

    全球變暖已經(jīng)成為不爭的事實[1]。近年來,眾多學者使用歷史氣象數(shù)據(jù)對全國和區(qū)域尺度氣溫的變化情況進行了廣泛而深入的研究,研究區(qū)域包括中國全境[2]、青藏高原[3]、黃土高原[4]、華北平原[5]、黃河流域[6]、長江流域[7]、渭河流域[8]、漢江流域[9]、陜西關中[10]、和陜南地區(qū)[10]等。研究結果表明我國大部分地區(qū)存在顯著升溫趨勢,但升溫的快慢程度存在明顯區(qū)域差異。根據(jù)氣象學有關理論,氣溫升高會加速空氣中水汽分子運動速度從而加速自然氣候的蒸發(fā)能力[11- 13]。但事實卻不盡如此,在全國尺度[2,14]、干旱半干旱地區(qū)[15]、黃土高原[4]、黃河流域[6]、長江流域[7]、渭河流域[8]、漢江流域[9]、陜西關中[10]等地,無論是蒸發(fā)皿蒸發(fā)量還是根據(jù)Penman-Monteith公式計算得到的潛在蒸發(fā)量都表現(xiàn)出了波動下降的趨勢即所謂的“蒸發(fā)悖論”現(xiàn)象。部分學者還對蒸發(fā)下降的可能成因進行了探討,他們發(fā)現(xiàn)與蒸發(fā)密切相關的氣象要素的變化(風速減小和太陽輻射減少[2- 16]等)和人類活動(農(nóng)業(yè)灌溉和下墊面變化[9- 16]等)可以導致蒸發(fā)量的變化。目前,關于“蒸發(fā)悖論”的機理尚無定論,對于蒸發(fā)下降的原因歸納起來主要包括:云量和氣溶膠等污染物的增加引起的太陽輻射下降[2- 16],空氣濕度的增加導致的水汽壓差減小[5- 16],夏季風變化引起的風速下降等[2- 16]。但從全球范圍來看,云量和氣溶膠增加引起的輻射量(潛熱通量)下降是蒸發(fā)減少的主要原因[12- 13]。

    蒸發(fā)是大氣中水汽的重要來源,而絕對濕度則是衡量大氣干濕程度(水汽含量)的物理量,具體是指單位體積空氣中所含水蒸汽的質(zhì)量,即水汽密度(g/m3),絕對濕度的變化可以從側面反映蒸發(fā)能力的變化。土壤水分蒸發(fā)理論表明,地表的蒸發(fā)能力與蒸發(fā)皿的蒸發(fā)能力成正比,蒸發(fā)皿蒸發(fā)量下降的情況下,地表的蒸發(fā)量也應相應減小,而地表蒸發(fā)又是大氣水汽的重要來源[17],進而造成大氣絕對水汽含量的下降。秦嶺南北地區(qū)主體位于我國東部季風區(qū)西緣,與青藏高原氣候區(qū)和西北干旱區(qū)交匯,氣候條件較為復雜。該區(qū)整體和各子區(qū)在1960—2011年間經(jīng)歷了明顯的升溫過程[18],而潛在蒸發(fā)量卻表現(xiàn)出了明顯的下降趨勢,即“蒸發(fā)悖論”現(xiàn)象[10]。此外,該區(qū)干旱事件發(fā)生頻率較高,氣候、農(nóng)業(yè)和水文干旱也日趨嚴重。過去52年間降水量微弱減少,以干旱指數(shù)(降水量比上蒸發(fā)量)度量的干旱程度有加重趨勢,區(qū)域整體特別是秦嶺以北地區(qū)有向暖干化發(fā)展的態(tài)勢[19]。到目前為止,關于秦嶺南北地區(qū)蒸發(fā)量的已有研究[7, 9- 10, 16,19]更多地是考慮多個氣象要素對蒸發(fā)過程的共同作用。通過相關分析、偏相關分析、敏感性分析、歸因分析等方法探究蒸發(fā)下降的主導因子,所選用的指標也大多是降水、相對濕度和水汽壓等,較少涉及到絕對濕度。絕對濕度、相對濕度和水汽壓3個指標并不相同,所代表的含義也不一致。相對濕度與絕對濕度相比,前者只是表征空氣中水汽的相對含量并不能反映出水汽密度的具體值,而且受溫度影響較大。而絕對濕度則不然,其不受溫度影響,可以較好地量化水汽密度的變化。基于上述認識,本研究利用秦嶺南北地區(qū)的歷史氣象觀測資料,分別計算了絕對濕度和潛在蒸發(fā)量,并分析了絕對濕度的變化趨勢。重點探討氣候變暖和蒸發(fā)下降(“蒸發(fā)悖論”)背景下大氣水汽資源(絕對濕度)的變化趨勢,絕對濕度和潛在蒸發(fā)究竟是同向變化還是逆向變化,細節(jié)特征又如何,并嘗試利用蒸發(fā)互補理論對其變化的可能成因給出了初步的說明和解釋。

    1 研究區(qū)概況

    圖1 研究區(qū)范圍及氣象觀測站點分布 Fig.1 Location of the study area and distribution of meteorological stations

    秦嶺橫亙于于我國中部,東西走向,山脈以北屬黃河流域,氣候較為干旱;山脈以南為長江流域,氣候溫暖濕潤。在氣候、植被和土壤分布上,山脈南北存在明顯差異,因而,秦嶺一直被認為是我國境內(nèi)的一條重要地理分界線。秦嶺的范圍有廣義和狹義之分,狹義的秦嶺只限于陜西南部、渭河與漢江之間的山區(qū);而廣義的秦嶺東西長約1500km,西接昆侖山,起自岷山以北,向東延伸經(jīng)由甘肅境內(nèi)、陜西南部到達河南境內(nèi)伏牛山一帶。本研究參考周旗等[18]的研究成果,將研究區(qū)域劃分成4個子區(qū),以研究絕對濕度變化的區(qū)域分異特點。分別是秦嶺北坡及其以北的暖溫帶地區(qū)(以下簡稱“秦嶺以北”);伏牛山及其以東平原(因大部分區(qū)域屬秦嶺南坡山地,以下簡稱“秦嶺南坡”);秦嶺以南的漢水谷地、巴山、涢水谷地及淮河上游北亞熱帶地區(qū)(因大部分區(qū)域屬漢水流域,下簡稱“漢水流域”);巴山南麓、巫山谷地及江漢平原西北部(下簡稱“巴巫谷地”),研究區(qū)域和氣象站點的空間分布如圖1。秦嶺地區(qū)具有大陸性季風氣候特點,年均氣溫12—17℃,≥10℃的年積溫為3700—4900℃,年均降水量600—1200mm,降水變率大,季節(jié)分配不均勻,無霜期200—250d[18- 19]。

    2 數(shù)據(jù)及方法

    基于氣象觀測資料的可獲取性(比濕和露點溫度資料不可得),選用逐日平均氣溫和水汽壓實測資料推算絕對濕度。數(shù)據(jù)獲取自中國氣象科學數(shù)據(jù)共享服務網(wǎng)(http://cdc.cma.gov.cn/home.do),主要包括1960—2011年間秦嶺南北地區(qū)45個氣象臺站的逐日最高溫度、最低溫度、平均氣溫、相對濕度、水汽壓、日照時數(shù)和平均風速數(shù)據(jù)。1960—2001年20cm口徑蒸發(fā)皿蒸發(fā)量逐日觀測資料用于檢驗Penman-Monteith公式估算結果的準確性,獲取自黃土高原地區(qū)氣候數(shù)據(jù)庫(http://www.loess.csdb.cn/pdmp/foreground/viewMetadata.action?id=8)。上述觀測資料在整編發(fā)布過程中經(jīng)過了嚴格的質(zhì)量控制,本研究在數(shù)據(jù)預處理過程中剔除了缺測年份大于10a的站點,并用線性回歸法和臨近站點插補法對漏測數(shù)據(jù)進行了完善。

    依據(jù)上述方法分別計算各臺站日均絕對濕度,采用算術平均值法分別統(tǒng)計各站點絕對濕度在不同時間尺度的平均值,并進一步求均值獲得不同子區(qū)和區(qū)域整體絕對濕度的長期變化序列。季節(jié)劃分按照慣例,以3—5月為春季,6—8月為夏季,9—11月為秋季,12月至次年2月為冬季。為了揭示絕對濕度和蒸發(fā)量的時變特征和響應關系,采用目前被廣泛認可,準確性較高,由世界糧農(nóng)組織(FAO)推薦的Penman-Monteith公式估算潛在蒸發(fā)量,估算結果與蒸發(fā)皿蒸發(fā)量觀測數(shù)據(jù)進行了對比,以檢驗其估算準確性。Penman-Monteith公式的具體計算過程和結果檢驗參見文獻[10, 19]。

    3 結果與分析

    3.1 絕對濕度的空間分布特征

    秦嶺南北地區(qū)年尺度和不同季節(jié)絕對濕度呈現(xiàn)南高北低的空間分布格局,具有較好的海拔地帶性和緯度地帶性,即隨著海拔和緯度的上升(下降)而相應減少(增加)(圖2,圖3)。各子區(qū)濕度大小排序為巴巫谷地(12.0 g/m3)>漢水流域(10.6 g/m3)>秦嶺南北(10.0 g/m3)>秦嶺南坡(9.2 g/m3)>秦嶺以北(8.3 g/m3),各子區(qū)濕度相對較高的站點依次為萬州(13.3 g/m3)>鐘祥(11.8 g/m3)>西華(10.4 g/m3)>開封(9.7 g/m3),相對較小的站點依次為廣元(10.3 g/m3)>欒川(8.3 g/m3)>武都(8.1 g/m3)>洛川(6.7 g/m3);季節(jié)尺度上濕度的分布規(guī)律與年尺度基本相同(圖略),排序為夏季(17.3 g/m3)>秋季(10.1 g/m3)>春季(8.8 g/m3)>冬季(4.1 g/m3)。上述空間分布格局主要是由于秦嶺南北的大部分地區(qū)受季風影響,冬季風來自高緯度大陸,水汽較少,而夏季風來自低緯度海洋,高溫而潮濕,造成了濕度的年內(nèi)變化[21- 23]。

    圖2 年平均絕對濕度的海拔地帶性和緯度地帶性Fig.2 The altitudinal and latitudinal zonality of yearly averaged absolute humidity

    圖3 年度和夏季平均絕對濕度空間分布Fig.3 Spatial distribution of yearly and summer averaged absolute humidity

    圖4 年平均絕對濕度和累積距平變化時間序列Fig.4 Time series of yearly averaged absolute humidity and its accumulative anomaly

    3.2 絕對濕度的時空變化

    3.2.1 年平均絕對濕度年際變化

    秦嶺南北整體和各子區(qū)濕度變化趨勢基本一致,除巴巫谷地(-0.04 g m-310a-1)不顯著下降外,其它子區(qū)均呈現(xiàn)出增加趨勢(圖4),增加速率排序為秦嶺南坡(0.080 g m-310a-1)>秦嶺以北(0.028 g m-310a-1)>漢水流域(0.026 g m-310a-1)>秦嶺南北(0.023 g m-310a-1)。由時間變化曲線(圖4)和累積距平曲線(圖4)可知,研究區(qū)整體和各子區(qū)濕度呈現(xiàn)波動變化,各站點變差系數(shù)介于0.05(長武)和0.02(梁平)之間,區(qū)域平均變差系數(shù)分別為秦嶺以北0.033、秦嶺南坡0.034、漢水流域0.025、巴巫谷地0.019和秦嶺南北0.024,其中以巴巫谷地的波動幅度最小。近52年間,絕對濕度最大值出現(xiàn)在1998年,最小值出現(xiàn)在1966年,1986年和1998年是變化的轉折點,1960—1986年以微弱下降為主,此后直至1998年震蕩上升,1998年以后各子區(qū)變化步調(diào)基本保持一致,呈同步下降趨勢。

    年尺度上,分別有51%、30%和13%的站點達到95%、99%和99.9%的顯著水平(圖5)。秦嶺以北上升速率為0.03 g m-310a-1,上升和下降的站點分別占到58%和42%,以武功的上升速度最快;秦嶺南坡上升速率為0.08 g m-310a-1,所有站點濕度都上升,以鎮(zhèn)安的上升速度最快;漢水流域上升速率為0.03 g m-310a-1,上升和下降的站點分別占到86%和14%,以石泉的上升速度最快;巴巫谷地的下降速率為-0.04 g m-310a-1,上升和下降的站點分別占到30%和70%,以閬中的上升速度最快,巴中的下降速度最快。

    3.2.2 季平均絕對濕度的年際變化

    春季,區(qū)域整體和秦嶺南坡(0.041 g m-310a-1)以外的各子區(qū)絕對濕度均表現(xiàn)出下降趨勢,下降速率排序為巴巫谷地(-0.070 g m-310a-1)>秦嶺以北(-0.034 g m-310a-1)>漢水流域(-0.026 g m-310a-1)>秦嶺南北(-0.022 g m-310a-1),其中巴巫谷地的下降趨勢通過了95%的顯著性檢驗。各子區(qū)濕度下降站點所占比例排序為巴巫谷地(80%)>漢水流域(71%)>秦嶺南北(62%)>秦嶺以北(58%)>秦嶺南坡(33%),其中16%站點的下降趨勢達到95%及以上顯著水平(圖5);

    夏季,區(qū)域整體和巴巫谷地(-0.112 g m-310a-1)以外的各子區(qū)絕對濕度均表現(xiàn)出上升趨勢,上升速率排序為秦嶺南坡(0.151 g m-310a-1)>漢水流域(0.037 g m-310a-1)>秦嶺南北(0.020 g m-310a-1)>秦嶺以北(0.005 g m-310a-1),其中秦嶺南坡的上升趨勢通過了95%的顯著性檢驗。各子區(qū)濕度上升站點所占比例排序為秦嶺南坡(100%)>漢水流域(79%)>秦嶺南北(62%)>秦嶺以北(58%)>巴巫谷地(10%),其中22%站點的上升趨勢達到95%及以上顯著水平(圖5);

    秋季,各子區(qū)絕對濕度呈上升趨勢,上升速率排序為秦嶺南坡(0.115 g m-310a-1)>漢水流域(0.059 g m-310a-1)>秦嶺南北(0.051 g m-310a-1)>秦嶺以北(0.017 g m-310a-1)>巴巫谷地(0.015 g m-310a-1),其中秦嶺南坡的上升趨勢通過了95%的顯著性檢驗。各子區(qū)濕度上升站點所占比例排序為秦嶺南坡(89%)>漢水流域(79%)>秦嶺南北(64%)>巴巫谷地(50%)=秦嶺以北(50%),其中18%站點的上升趨勢達到95%及以上顯著水平(圖5);

    冬季,各子區(qū)絕對濕度均表現(xiàn)出上升趨勢,上升速率排序為秦嶺南坡(0.093 g m-310a-1)>秦嶺南北(0.060 g m-310a-1)>漢水流域(0.056 g m-310a-1)>秦嶺以北(0.048 g m-310a-1)>巴巫谷地(0.042 g m-310a-1),其中秦嶺南坡和秦嶺南北整體的上升趨勢通過了95%的顯著性檢驗。各子區(qū)濕度上升站點所占比例排序為秦嶺以北(100%)=秦嶺南坡(100%)>秦嶺南北(93%)>巴巫谷地(90%)>漢水流域(86%),其中44%站點的上升趨勢達到95%及以上顯著水平(圖5)。

    圖5 近52年年度和季節(jié)平均絕對濕度變化的空間分布Fig.5 Spatial distribution of yearly and seasonal averaged absolute humidity variation in recent 52a

    3.3 絕對濕度和潛在蒸發(fā)量變化的比較分析

    蒸發(fā)是大氣中水汽的重要來源,而絕對濕度則是衡量大氣干濕程度(水汽含量)的物理量,因此,絕對濕度的變化可以從側面反映蒸發(fā)能力的變化。在存在明顯“蒸發(fā)悖論”現(xiàn)象的背景下,本研究試圖從絕對濕度角度分析大氣水汽資源的變化趨勢,以及絕對濕度和潛在蒸發(fā)變化的響應關系。參考3.2.1小節(jié)的分析結果將絕對濕度變化劃分為3個時間段,即1960—1986年、1987—1998年和1999—2011年。同理根據(jù)潛在蒸發(fā)量累積距平曲線將其劃分為1960—1989年和1990—2011年兩個時間段進行分析。區(qū)域整體和各子區(qū)潛在蒸發(fā)量和絕對濕度相關系數(shù)見表1,兩者不同時段氣候傾向率及其顯著性檢驗結果見表2。由表1可知,年尺度上,除巴巫谷地外,絕對濕度和潛在蒸發(fā)量均表現(xiàn)出顯著的負相關關系,即反向變化趨勢;季節(jié)尺度上,春季和秋季的相關系數(shù)都是負值,各子區(qū)無一例外,且大部分達到顯著水平;夏季和冬季的規(guī)律較為接近,漢水流域和巴巫谷地呈不顯著正相關,其它區(qū)域均為負相關。此外,年度、春季和秋季兩個指標的相關系數(shù)都表現(xiàn)出由北向南順次遞減的變化規(guī)律,說明濕度和蒸發(fā)的緊密程度有明顯地域差異,隨著子區(qū)的南移而順次減弱。

    表1 潛在蒸發(fā)量與絕對濕度相關系數(shù)Table 1 Correlation coefficients between potential evaporation and absolute humidity

    **0.01顯著水平0.01 significant level; *0.1顯著水平0.1 significant level

    NRQ: northern regions of Qinling Mountains; SSQ: southern slope of Qinling Mountains; HRB: Han River Basin; BWV: Bawu Valley; NSQ: northern and southern regions of Qinling Mountains

    由表2可知,1960—2011年間,年度、春季和冬季的潛在蒸發(fā)量經(jīng)歷了先降后升的變化階段,1960—1989年下降,1990—2011年上升。年尺度上,前30a蒸發(fā)量的下降速率要明顯大于后22a的上升速率;春季則剛好相反,后22a的速率大于前30a;冬季蒸發(fā)變化在這兩個時段的差異相對較小;夏季和秋季的蒸發(fā)量變化與年尺度不同,兩個時段和52a整體都表現(xiàn)出了下降趨勢。對于夏季而言,前30a的下降速率要大于52a整體和后22a;秋季,不同時段間速率變化差異不大。絕對濕度方面,近52年經(jīng)歷了“降—升—降”的變化過程,除春季外,其他季節(jié)和年尺度的傾向率均為正值,1960—1986年和1999—2011年下降,而1986—1998年增加。其中,后13a的下降速度要明顯大于前27a,各季節(jié)和子區(qū)的變化規(guī)律基本一致。

    表2 年度和季節(jié)潛在蒸發(fā)量和絕對濕度變化傾向率Table 2 The slope ratio of annual and seasonal potential evaporation and absolute humidity

    ***0.001顯著水平; **0.01顯著水平; *0.1顯著水平

    近52年年度和季節(jié)潛在蒸發(fā)和絕對濕度呈反向變化趨勢,區(qū)域整體和各子區(qū)基本一致。其中,年度和春季以外的其它季節(jié)濕度上升,蒸發(fā)下降,而春季表現(xiàn)出來的規(guī)律則相反,濕度下降而蒸發(fā)上升;1960—1989年間,年度和季節(jié)尺度的蒸發(fā)都在下降,同期的濕度也在下降(個別區(qū)域上升),兩者基本保持同向變化趨勢;1990—2011年間,年度和春、冬兩季的蒸發(fā)都在上升,而同期的濕度則是先同向變化(上升)再反向變化(下降),區(qū)域整體和各子區(qū)的規(guī)律基本一致??臻g分布特征方面(圖6),52a整體兩者呈反向變化趨勢的站點所占比例分別為年度69%、春季80%、夏季62%、秋季67%和冬季78%。春季、秋季和冬季的濕度和蒸發(fā)反向變化站點均勻分布,沒有明顯的區(qū)域聚集特征;而在年尺度上,同向變化的站點主要集中于巴巫谷地;在夏季,濕度和蒸發(fā)同時下降的站點也集中于這一區(qū)域。

    3.4 絕對濕度和潛在蒸發(fā)量反向變化的成因分析

    1960—2011年間,潛在蒸發(fā)量呈下降趨勢,而絕對濕度表現(xiàn)出上升趨勢。潛在蒸發(fā)與絕對濕度負相關關系明顯的現(xiàn)象,反映出了蒸發(fā)皿蒸發(fā)量和潛在蒸發(fā)量的代表性問題。事實上,蒸發(fā)皿蒸發(fā)只是有限水面的自由蒸發(fā),嚴格意義上它只代表一個地區(qū)接受太陽能量多少,而不能代表實際蒸發(fā)量,這一點在干旱和半干旱地區(qū)體現(xiàn)得尤為明顯。絕對濕度反映了空氣中的實際水汽含量,而水汽含量的高低又直接取決于地表實際蒸發(fā)量,換言之,絕對濕度反映了地面實際蒸發(fā)量。蒸發(fā)互補理論假定(圖7),在給定的輻射條件下,當充分供水時實際蒸散發(fā)量與潛在蒸散發(fā)量相等;當下墊面供水量減少時,實際蒸散發(fā)量會減少,從而釋放出更多的能量成為顯熱,從而導致潛在蒸散發(fā)增加[24- 30]。實際蒸發(fā)量受水分與能量條件的控制。如果能量條件固定,互補關系成立;如果能量條件變化不大,互補關系仍然成立;如果能量條件變化較大,互補關系可能發(fā)生改變[28- 30]。

    圖6 近52年年度和季節(jié)絕對濕度和潛在蒸發(fā)量變化趨勢的空間分布Fig.6 Spatial distribution of absolute humidity and potential evaporation change trend in recent 52a and seasons

    就全國而言,在以干旱為主的區(qū)域決定實際蒸發(fā)最重要的因子是水分,而在以濕潤為主的區(qū)域,則是能量[24]。王艷君等[7]分析了長江流域實際蒸發(fā)量和潛在蒸發(fā)量的關系,并利用干燥度指數(shù)R(潛在蒸發(fā)量與降水的比值)來判定研究區(qū)域的干濕條件。當R<0.8時,實際蒸發(fā)量與潛在蒸發(fā)量為明顯的正相關關系,此時氣候為濕潤條件,實際蒸發(fā)量主要受制于其它氣象要素的變化,與降水關系不大;當R>1.0時,即干旱環(huán)境,實際蒸發(fā)量與潛在蒸發(fā)量為明顯的互補關系,實際蒸發(fā)量主要受降水控制。Golubev等[25]通過對美國和前蘇聯(lián)8個區(qū)域的蒸發(fā)皿蒸發(fā)量和實際蒸發(fā)量的變化研究得到類似的結論,認為當R<0.7時,蒸發(fā)皿蒸發(fā)量與實際蒸發(fā)量的變化表現(xiàn)為一致趨勢;當R≥0.8時,兩者為互補關系。Cong等[26]對中國十大流域水文變化趨勢的研究結果表明,在我國北方干旱區(qū)實際蒸發(fā)量的下降主要由降水變化引起,南方濕潤區(qū)實際蒸發(fā)量的下降主要受潛在蒸發(fā)量的影響。Teuling等[27]和Roderick等[11]指出,在能量受限條件下,下降的蒸發(fā)皿蒸發(fā)量通常意味著實際蒸發(fā)量的下降;在水分受限條件下,即干旱缺水地區(qū),下降的蒸發(fā)皿蒸發(fā)量常常意味著實際蒸發(fā)量的增加。

    圖7 蒸發(fā)互補關系曲線Fig.7 Curve of evaporation complementary relationship

    本研究所關注的秦嶺南北地區(qū)地跨兩個氣候帶,秦嶺以北的黃河流域(黃土高原)屬于典型的干旱半干旱地帶,而秦嶺以南和漢水流域為半干旱半濕潤地區(qū),巴巫谷地嚴格意義上講屬于濕潤易旱地區(qū)(降水充沛但季節(jié)性干旱頻發(fā))。參考蔣沖等[19]的研究成果和王艷君等[7]對區(qū)域干濕條件的劃分標準,界定當R<0.8時為濕潤條件,而R>1.0時為干旱環(huán)境,由此可知巴巫谷地(R=0.86)為濕潤區(qū),而漢水流域(1.06)、秦嶺南坡(1.25)和秦嶺以北(1.64)為干旱環(huán)境,各子區(qū)由北向南潛在蒸發(fā)量受水分的限制作用不斷減弱,能量限制不斷加強。結合干旱指數(shù)和表1可知,秦嶺以北、秦嶺南坡、漢水流域和巴巫谷地部分站點絕對濕度(實際蒸發(fā)量)和潛在蒸發(fā)量為互補關系(顯著負相關),由北向南隨著水分限制作用的不斷減弱兩者的負相關關系也逐漸減弱,直至巴巫谷地的正相關,各子區(qū)潛在蒸發(fā)和絕對濕度變化趨勢相反的站點個數(shù)也隨著區(qū)域的南移而逐漸減少;季節(jié)尺度上,干燥度指數(shù)由北向南順次遞減,蒸發(fā)和濕度的相關關系也基本上由(顯著)負相關轉變?yōu)檎嚓P或不顯著負相關。上述結果符合蒸發(fā)互補理論,也和王艷君等[7]和Cong等[26]的研究結果趨勢上基本一致。

    4 結論

    (1)秦嶺南北地區(qū)年尺度和不同季節(jié)絕對濕度由南向北順次遞減,具有較好的海拔地帶性和緯度地帶性,各子區(qū)按濕度大小排序為巴巫谷地>漢水流域>秦嶺南北>秦嶺南坡>秦嶺以北。季節(jié)平均濕度以夏季為最大,冬季最小。

    (2)秦嶺南北整體和各子區(qū)絕對濕度變化趨勢基本一致,除巴巫谷地表現(xiàn)出不顯著的下降趨勢外,其它子區(qū)均呈現(xiàn)出增加趨勢,增加速率排序為秦嶺南坡>秦嶺以北>漢水流域>秦嶺南北。1986年和1998年是絕對濕度變化的轉折點,1960—1986年以微弱下降為主,此后直至1998年震蕩上升,1998年以后呈下降趨勢。

    (3)年尺度和春、秋兩季,除巴巫谷地外,絕對濕度和潛在蒸發(fā)量均表現(xiàn)出顯著的負相關關系,而在夏季和冬季,除漢水流域和巴巫谷地呈不顯著正相關以外,其它區(qū)域均為負相關。年度和春、秋兩季兩個指標負相關的緊密程度隨著區(qū)域的南移而逐漸減弱。

    (4)1960—2011年間,年度和季節(jié)尺度潛在蒸發(fā)和絕對濕度反向變化;而在1960—1989年間,兩者同向變化;1990—2011年間,年度和春、冬兩季潛在蒸發(fā)上升,而同期的絕對濕度則是先上升再下降。

    (5)潛在蒸發(fā)和絕對濕度變化的相反趨勢驗證了“蒸發(fā)悖論”的普遍存在,可以用蒸發(fā)互補理論來解釋。實際蒸發(fā)量的增加導致空氣中水汽含量增加,反過來抑制了水面蒸發(fā)(潛在蒸發(fā)量)。秦嶺以北、秦嶺南坡、漢水流域和巴巫谷地部分地區(qū)絕對濕度和潛在蒸發(fā)量為互補關系,由北向南隨著水分限制作用的不斷減弱兩者逐漸轉化為不顯著的正相關關系;季節(jié)尺度上,潛在蒸發(fā)和絕對濕度的相關關系也基本上由(顯著)負相關轉變?yōu)檎嚓P或不顯著負相關。

    [1] Solomon S, Qin D, Manning M, Chen Z, Marquis M, Averyt K B, Tignor M, Miller H L. Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment of the Intergovernmental Panel on Climate Change. New York: Cambridge University Press, 2007: 10- 10.

    [2] Yin Y H, Wu S H, Chen G, Dai E F. Attribution analyses of potential evapotranspiration changes in China since the 1960s. Theoretical and Applied Climatology, 2010, 101(1/2): 19- 28.

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

    [4] Li Z, Zheng F L, Liu W Z. Spatiotemporal characteristics of reference evapotranspiration during 1961- 2009 and its projected changes during 2011- 2099 on the Loess Plateau of China. Agricultural and Forest Meteorology, 2012, 154- 155: 147- 155.

    [5] Gao G, Xu C Y, Chen D L, Singh V P. Spatial and temporal characteristics of actual evapotranspiration over Haihe River basin in China. Stochastic Environmental Research and Risk Assessment, 2012, 26(5): 655- 669.

    [6] Wang W G, Shao Q X, Peng S Z, Xing W Q, Yang T, Luo Y F, Yong B, Xu J Z. Reference evapotranspiration change and the causes across the Yellow River basin during 1957- 2008 and their spatial and seasonal differences. Water Resources Research, 2012, 48(5): 27- 35.

    [7] 王艷君, 劉波, 翟建青, 蘇布達, 羅勇, 張增信. 長江流域潛在蒸發(fā)量和實際蒸發(fā)量的關系. 氣候變化研究進展, 2011, 7(6): 393- 399.

    [8] Zuo D P, Xu Z X, Yang H. Spatiotemporal variations and abrupt changes of potential evapotranspiration and its sensitivity to key meteorological variables in the Wei River basin, China. Hydrological Processes, 2011, 26(8): 1149- 1160.

    [9] 張東, 張萬昌, 徐全芝. 漢江上游流域蒸散量計算方法的比較及改進. 資源科學, 2005, 27(1): 97- 103.

    [10] 蔣沖, 王飛, 穆興民, 李銳. 秦嶺南北潛在蒸散量時空變化及突變特征分析. 長江流域資源與環(huán)境, 2013, 22(5): 573- 580.

    [11] Robock A, Konstantin V Y, Srinivasan G, Entin J K, Hollinger S E, Speranskaya N A, Liu S, Namkhai A. The global soil moisture data bank. Bulletin of the American Meteorological Society, 2000, 81(6): 1281- 1299.

    [12] McVicar T R, Roderick M L, Donohue R J, Li L T, Van Niel T G, Thomas A, Grieser J, Jhajharia D, Himri Y, Mahowald N M, Mescherskaya A V, Kruger A C, Rehman S, Dinpashoh Y. Global review and synthesis of trends in observed terrestrial near-surface wind speeds: Implications for evaporation. Journal of Hydrology, 2012, 416- 417: 182- 205.

    [13] Roderick M L, Farquhar G D. The cause of decreased pan evaporation over the past 50 years. Science, 2002, 298(5597): 1410- 1411.

    [14] Yang H B, Yang D W. Climatic factors influencing changing pan evaporation across China from 1961 to 2001. Journal of Hydrology, 2012, 414- 415(2): 184- 193.

    [15] Han S J, Xu D, Wang S L. Decreasing potential evaporation trends in China from 1956 to 2005: Accelerated in regions with significant agricultural influence. Agricultural and Forest Meteorology, 2012, 154- 155: 44- 56.

    [16] 李鵬飛, 孫小明, 趙昕奕. 近50年中國干旱半干旱地區(qū)降水量與潛在蒸散量分析. 干旱區(qū)資源與環(huán)境, 2012, 26(7): 57- 63.

    [17] 徐彥偉, 康世昌, 張玉蘭, 張擁軍. 夏季納木錯湖水蒸發(fā)對當?shù)卮髿馑暙I的方法探討: 基于水體穩(wěn)定同位素的估算. 科學通報, 56(13): 1042- 1049.

    [18] 周旗, 卞娟娟, 鄭景云. 秦嶺南北1951- 2009年的氣溫與熱量資源變化. 地理學報, 2011, 66(9): 1211- 1218.

    [19] 蔣沖, 王飛, 穆興民, 李銳. 氣候變化對秦嶺南北植被凈初級生產(chǎn)力的影響(Ⅰ)——近52年秦嶺南北氣候時空變化特征分析. 中國水土保持科學, 2012, 10(5): 56- 63.

    [20] 王名才. 大氣科學常用公式. 北京: 氣象出版社, 1994: 2- 382.

    [21] 劉國緯. 水文循環(huán)的大氣過程. 北京: 科學出版社, 1997: 45- 50.

    [22] 錢璇, 姚永強, 李俊榮, 尹佳, 張永婧, 劉立勇, 王紅帥, 周云賀, 李林, 尤顯龍, 馬江龍. 全國大氣水汽含量分布. 天文學報, 2012, 53(4): 325- 340.

    [23] 張丹, 劉昌明, 付永鋒, 邱新法, 劉小莽. 基于MODIS數(shù)據(jù)的中國地面水汽壓模擬與分析. 資源科學, 2012, 34(1): 74- 80.

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

    [25] Golubev V S, Lawrimore J H, Groisman P Y, Speranskaya N A, Zhuravin S A, Menne M J, Peterson T C, Malone R W. Evaporation changes over the contiguous United States and the former USSR: a reassessment. Geophysical Research Letters, 2001, 28(13): 2665- 2668.

    [26] Cong Z T, Yang D W, Ni G H. Does evaporation paradox exist in China. Hydrology and Earth System Science, 2009, 13: 357- 366.

    [27] Teuling A J, Hirschi M, Ohmura, Wild M, Reichstein M, Ciais P, Buchmann N, Ammann, Montagnani L, Richardson D, Wohlfahrt G, Seneviratne S I. A regional perspective on trends in continental evaporation. Geophysical Research Letters, 2009, 36(2), L02404, doi: 10.1029/2008GL036584.

    [28] 楊漢波, 楊大文, 雷志棟, 孫福寶, 叢振濤. 蒸發(fā)互補關系在不同時間尺度上的變化規(guī)律及其機理. 中國科學E輯: 技術科學, 2009, 39(2): 333- 340.

    [29] 韓松俊, 胡和平, 楊大文, 劉群昌. 塔里木河流域山區(qū)和綠洲潛在蒸散發(fā)的不同變化及影響因素. 中國科學E輯: 技術科學, 2009, 39(8): 1375- 1383.

    [30] 楊漢波, 楊大文, 雷志棟, 孫福寶, 叢振濤. 蒸發(fā)互補關系的區(qū)域變異性. 清華大學學報: 自然科學版, 2008, 48(9): 1413- 1416.

    Spatial and temporal variation of absolute humidity and its relationship with potential evaporation in the northern and southern regions of Qinling Mountains

    JIANG Chong1,2, MU Xingmin1, MA Wenyong3, YU Xinyang3, LIU Xianfeng4, LI Jianguo5, LIU Sijie6, WANG Fei1,*

    1InstituteofSoilandWaterConservation,NorthwestA&FUniversity,Yangling712100,China2TheStateKeyLaboratoryofEarthSurfaceProcessesandResourceEcology,CollegeofGlobalChangeandEarthSystemScience,BeijingNormalUniversity,Beijing100875,China3InstituteofGeographicSciencesandNaturalResourcesResearch,ChineseAcademyofSciences,Beijing100101,China4CollegeofResourceScienceandTechnology,BeijingNormalUniversity,Beijing100875,China5ColdandAridRegionsEnvironmentalandEngineeringResearchInstitute,ChineseAcademyofSciences,Lanzhou730000,China6InstituteofRSandGIS,PekingUniversity,Beijing100871,China

    Based on the meteorological data from weather stations in the northern and southern regions of Qinling Mountains (NSQ), we used absolute humidity estimating formula and Penman-Monteith equation to estimate absolute humidity and potential evaporation, respectively. Furthermore, we analyzed the spatial and temporal distribution and variation of absolute humidity, and tried to give the possible explanations for the reciprocal feedback between absolute humidity and potential evaporation, the result are as following. (1) Absolute humidity presented a clear pattern that high in south and low in north, the order of it was Bawu Valley (BWV), Han River Basin (HRB), northern and southern regions of Qinling Mountains(NSQ), southern slope of Qinling Mountains (SSQ) and northern regions of Qinling Mountains (NRQ), the order of seasonal absolute humidity was summer, autumn, spring and winter. (2) In the past 52a, absolute humidity in most regions, except for BWV, had increasing trend. 1986 and 1998 are the change point, in the period of 1960—1986, it had insignificant decreasing trend, and then increased significantly until 1998, and decreased since 1998. (3) On the annual scale, spring and autumn, in most regions, except for BWV, absolute humidity negatively correlated with potential evaporation, meanwhile, in summer and winter, except for HRB and BWV, absolute humidity and potential evaporation also had the similar correlation relationship. Beside from this, the correlation coefficients on the annual scale, spring and autumn became smaller with the latitude decreased. In the period of 1960—2011, on the annual and seasonal scale, absolute humidity and potential evaporation had contrary change trend, however, in the period of 1960—1989, they showed same change trend. From 1990 to 2011, potential evaporation in year, spring and winter increased significantly, meanwhile, absolute humidity increased initially and then decreased. (4) The contrary change trend of absolute humidity and potential evapotranspiration can be explained by the evaporation complementary theory. The actual evaporation increasing will lead to the increase of water vapor content in the air, also for the absolute humidity, which will inhibit water evaporation. Absolute humidity and potential evaporation had complementary relationship in NRQ, SSQ, HRB, and part of BWV, and the correlation relationship decreased with water limitation reducing, ranging from negative correlation to insignificant positive correlation.

    the northern and southern regions of Qinling Mountains; absolute humidity; spatial and temporal variation; potential evaporation; complementary relationship

    中國科學院知識創(chuàng)新工程重點部屬項目(KZZD-EW- 04);中國-荷蘭政府間聯(lián)合主題研究項目(中國科學院對外合作重點項目,GJHZ1018和NOW,OND1339291); 國家自然科學基金面上項目(41171420);中國科學院西部之光人才培養(yǎng)計劃聯(lián)合學者項目([2013]165)

    2014- 01- 19;

    2014- 10- 27

    10.5846/stxb201401190145

    *通訊作者Corresponding author.E-mail: wafe@ms.iswc.ac.cn

    蔣沖,穆興民,馬文勇,于新洋,劉憲鋒,李建國,劉思潔,王飛.秦嶺南北地區(qū)絕對濕度的時空變化及其與潛在蒸發(fā)量的關系.生態(tài)學報,2015,35(2):378- 388.

    Jiang C, Mu X M, Ma W Y, Yu X Y, Liu X F, Li J G, Liu S J, Wang F.Spatial and temporal variation of absolute humidity and its relationship with potential evaporation in the northern and southern regions of Qinling Mountains.Acta Ecologica Sinica,2015,35(2):378- 388.

    猜你喜歡
    子區(qū)谷地蒸發(fā)量
    The Kathmandu Valley
    加德滿都谷地
    基于MFD的高鐵站周圍路網(wǎng)誘導-控制方法
    考慮超級街區(qū)的城市路網(wǎng)邊界控制策略研究
    基于狄利克雷問題的路網(wǎng)控制子區(qū)動態(tài)劃分
    基于網(wǎng)絡能耗與交通效率的多子區(qū)控制模型
    1958—2013年沽源縣蒸發(fā)量變化特征分析
    1981—2010年菏澤市定陶區(qū)蒸發(fā)量變化特征分析
    新疆民豐縣地表水面蒸發(fā)量分析
    湟水谷地世界屋脊上的隱秘詩意
    cao死你这个sao货| av片东京热男人的天堂| www国产在线视频色| 国产精品电影一区二区三区| 侵犯人妻中文字幕一二三四区| 免费少妇av软件| 国产成人精品在线电影| 国产极品粉嫩免费观看在线| 美女福利国产在线| 一级黄色大片毛片| 国产高清激情床上av| 午夜亚洲福利在线播放| 久久久精品国产亚洲av高清涩受| 精品国产乱子伦一区二区三区| 国产精品一区二区精品视频观看| 日本 av在线| 成人亚洲精品一区在线观看| 精品熟女少妇八av免费久了| 50天的宝宝边吃奶边哭怎么回事| 国产亚洲精品久久久久久毛片| 精品国产一区二区三区四区第35| 亚洲五月天丁香| 中文字幕色久视频| 无限看片的www在线观看| 97碰自拍视频| 嫩草影院精品99| 日本一区二区免费在线视频| 在线观看舔阴道视频| 精品电影一区二区在线| 99在线人妻在线中文字幕| 黄片大片在线免费观看| 国产av一区在线观看免费| 色老头精品视频在线观看| 国产精品久久久av美女十八| 久久久久久久精品吃奶| 高清av免费在线| 久久精品亚洲熟妇少妇任你| 丁香欧美五月| 99国产精品一区二区蜜桃av| 色尼玛亚洲综合影院| 亚洲av日韩精品久久久久久密| www.www免费av| 亚洲欧洲精品一区二区精品久久久| 日韩国内少妇激情av| 99国产极品粉嫩在线观看| 亚洲精品粉嫩美女一区| 国产麻豆69| 中文字幕av电影在线播放| 黄网站色视频无遮挡免费观看| 国产主播在线观看一区二区| 交换朋友夫妻互换小说| 久久久久精品国产欧美久久久| 精品福利观看| 亚洲,欧美精品.| 一边摸一边抽搐一进一小说| 久久人人爽av亚洲精品天堂| 纯流量卡能插随身wifi吗| 亚洲精品美女久久久久99蜜臀| 大香蕉久久成人网| 精品久久久久久久久久免费视频 | 另类亚洲欧美激情| 黄片大片在线免费观看| 亚洲 欧美一区二区三区| 国产精品免费视频内射| 宅男免费午夜| 99香蕉大伊视频| 看免费av毛片| 自线自在国产av| 午夜福利免费观看在线| 少妇的丰满在线观看| 国产成+人综合+亚洲专区| 久久久久久人人人人人| 999精品在线视频| 国产色视频综合| 久久中文字幕人妻熟女| 国产有黄有色有爽视频| 中文字幕人妻丝袜制服| 亚洲七黄色美女视频| 午夜福利在线观看吧| 国产熟女午夜一区二区三区| 中文字幕色久视频| 免费在线观看日本一区| 国产免费男女视频| 久久精品91蜜桃| 精品久久久久久,| 国产高清激情床上av| 欧美日韩中文字幕国产精品一区二区三区 | 一区在线观看完整版| 亚洲在线自拍视频| 老汉色∧v一级毛片| 十八禁网站免费在线| 91国产中文字幕| 国产在线精品亚洲第一网站| 欧美成人午夜精品| 精品电影一区二区在线| 国产深夜福利视频在线观看| 亚洲黑人精品在线| 老鸭窝网址在线观看| 日韩高清综合在线| 一进一出抽搐动态| 精品日产1卡2卡| 日日夜夜操网爽| 激情在线观看视频在线高清| 777久久人妻少妇嫩草av网站| 每晚都被弄得嗷嗷叫到高潮| 久久人人精品亚洲av| 丰满的人妻完整版| 91成年电影在线观看| 在线观看免费日韩欧美大片| 亚洲精品国产一区二区精华液| 无限看片的www在线观看| 久久狼人影院| 国产黄色免费在线视频| 日日摸夜夜添夜夜添小说| 亚洲精品中文字幕在线视频| 日韩视频一区二区在线观看| 亚洲免费av在线视频| 亚洲自偷自拍图片 自拍| 免费少妇av软件| 亚洲精华国产精华精| 欧美大码av| 欧美一区二区精品小视频在线| 久久精品人人爽人人爽视色| 免费久久久久久久精品成人欧美视频| 少妇 在线观看| 长腿黑丝高跟| 无人区码免费观看不卡| 手机成人av网站| 欧美午夜高清在线| av福利片在线| 99国产精品一区二区三区| 99久久人妻综合| 免费女性裸体啪啪无遮挡网站| 高清av免费在线| 一区二区三区激情视频| 深夜精品福利| 在线观看日韩欧美| 亚洲av片天天在线观看| 另类亚洲欧美激情| 女人精品久久久久毛片| 久久精品国产综合久久久| 亚洲国产精品合色在线| 99久久国产精品久久久| 久久久国产一区二区| 久久久久久久久免费视频了| 色播在线永久视频| 欧美精品亚洲一区二区| 欧美日韩黄片免| 黑人巨大精品欧美一区二区蜜桃| 日本精品一区二区三区蜜桃| 精品国产乱码久久久久久男人| 亚洲人成电影观看| 国产区一区二久久| 国产99久久九九免费精品| 亚洲国产欧美网| 欧美日韩亚洲国产一区二区在线观看| 亚洲在线自拍视频| 亚洲欧美精品综合一区二区三区| 欧美老熟妇乱子伦牲交| 午夜福利在线免费观看网站| 亚洲精品国产一区二区精华液| 女性被躁到高潮视频| av欧美777| 久久精品91无色码中文字幕| 一区二区三区国产精品乱码| 日韩精品免费视频一区二区三区| 后天国语完整版免费观看| 久久国产亚洲av麻豆专区| 一级a爱片免费观看的视频| 麻豆av在线久日| 精品第一国产精品| 国产深夜福利视频在线观看| 日韩高清综合在线| 久久精品亚洲熟妇少妇任你| 五月开心婷婷网| 欧美一级毛片孕妇| 久久人人97超碰香蕉20202| 欧美日韩乱码在线| 最近最新中文字幕大全电影3 | 免费在线观看黄色视频的| av电影中文网址| 久久天堂一区二区三区四区| 免费不卡黄色视频| 97人妻天天添夜夜摸| 18禁国产床啪视频网站| 日本一区二区免费在线视频| 国产亚洲av高清不卡| 久久久久久大精品| 国产人伦9x9x在线观看| 亚洲 国产 在线| 国产精品久久久久久人妻精品电影| 人人妻人人添人人爽欧美一区卜| 男人舔女人下体高潮全视频| 国产国语露脸激情在线看| 曰老女人黄片| 黄色视频,在线免费观看| 在线播放国产精品三级| 涩涩av久久男人的天堂| 满18在线观看网站| 久久久久久久久中文| 日本欧美视频一区| 亚洲一卡2卡3卡4卡5卡精品中文| 757午夜福利合集在线观看| 亚洲激情在线av| 国产成人av教育| 国产精品乱码一区二三区的特点 | 日韩视频一区二区在线观看| 午夜免费鲁丝| 国产免费av片在线观看野外av| 三级毛片av免费| 老熟妇乱子伦视频在线观看| 日韩欧美三级三区| 欧洲精品卡2卡3卡4卡5卡区| 久久久久亚洲av毛片大全| 一个人观看的视频www高清免费观看 | 久久久久久亚洲精品国产蜜桃av| 久久香蕉精品热| 午夜福利在线免费观看网站| 久久青草综合色| av在线播放免费不卡| 欧美日韩一级在线毛片| 亚洲三区欧美一区| 欧美在线一区亚洲| 黑人操中国人逼视频| 欧美老熟妇乱子伦牲交| 欧美日韩亚洲国产一区二区在线观看| 国产男靠女视频免费网站| 黄片大片在线免费观看| 欧美日韩精品网址| 性色av乱码一区二区三区2| 我的亚洲天堂| 国产欧美日韩一区二区三| 热99re8久久精品国产| 欧美精品亚洲一区二区| 日本a在线网址| 黄色a级毛片大全视频| 国产又爽黄色视频| 99在线视频只有这里精品首页| 两性午夜刺激爽爽歪歪视频在线观看 | 精品熟女少妇八av免费久了| 黄网站色视频无遮挡免费观看| 免费女性裸体啪啪无遮挡网站| 免费观看精品视频网站| 在线观看一区二区三区| 很黄的视频免费| 精品少妇一区二区三区视频日本电影| 亚洲av第一区精品v没综合| 欧美成人免费av一区二区三区| 亚洲一区二区三区色噜噜 | 黄色a级毛片大全视频| 啦啦啦免费观看视频1| 麻豆国产av国片精品| 女人爽到高潮嗷嗷叫在线视频| 欧美老熟妇乱子伦牲交| 18禁黄网站禁片午夜丰满| 成年版毛片免费区| 欧美日韩乱码在线| 韩国av一区二区三区四区| 日本精品一区二区三区蜜桃| 多毛熟女@视频| 国产三级在线视频| 女警被强在线播放| 成人av一区二区三区在线看| 又大又爽又粗| 一二三四在线观看免费中文在| 亚洲人成77777在线视频| 日本vs欧美在线观看视频| 亚洲熟女毛片儿| 午夜亚洲福利在线播放| 亚洲欧美一区二区三区久久| 国产精品二区激情视频| 熟女少妇亚洲综合色aaa.| 伊人久久大香线蕉亚洲五| 大型av网站在线播放| 可以在线观看毛片的网站| 搡老岳熟女国产| 日韩有码中文字幕| 黄色 视频免费看| 1024视频免费在线观看| 成人永久免费在线观看视频| videosex国产| 亚洲 欧美一区二区三区| 亚洲视频免费观看视频| 国产单亲对白刺激| 无遮挡黄片免费观看| 国产欧美日韩一区二区精品| 亚洲欧美日韩另类电影网站| 最近最新免费中文字幕在线| 亚洲国产精品合色在线| 日日干狠狠操夜夜爽| 在线国产一区二区在线| 国产免费现黄频在线看| 一级毛片高清免费大全| 91成人精品电影| cao死你这个sao货| 亚洲精品国产一区二区精华液| 亚洲熟女毛片儿| 老汉色av国产亚洲站长工具| 国产欧美日韩一区二区精品| 中文字幕av电影在线播放| 国产精品影院久久| 亚洲aⅴ乱码一区二区在线播放 | 老司机午夜十八禁免费视频| 国产伦一二天堂av在线观看| 国产主播在线观看一区二区| 精品乱码久久久久久99久播| av福利片在线| 欧美日韩国产mv在线观看视频| 精品人妻1区二区| 一夜夜www| 国产精品久久电影中文字幕| 精品国产超薄肉色丝袜足j| 18禁国产床啪视频网站| 成人手机av| 我的亚洲天堂| 欧美黑人欧美精品刺激| 亚洲在线自拍视频| 国产精品一区二区三区四区久久 | 亚洲aⅴ乱码一区二区在线播放 | 午夜免费激情av| 免费在线观看视频国产中文字幕亚洲| 欧美日韩乱码在线| 国产免费男女视频| 亚洲精品中文字幕一二三四区| 欧美乱码精品一区二区三区| 视频区欧美日本亚洲| 人人妻人人添人人爽欧美一区卜| 免费高清视频大片| av国产精品久久久久影院| 欧美黑人欧美精品刺激| 精品国产乱子伦一区二区三区| 国产精品免费一区二区三区在线| xxx96com| 男女做爰动态图高潮gif福利片 | 脱女人内裤的视频| 亚洲精品粉嫩美女一区| 日日夜夜操网爽| 十八禁网站免费在线| 亚洲情色 制服丝袜| 黄色女人牲交| www.熟女人妻精品国产| 啦啦啦免费观看视频1| 黄色片一级片一级黄色片| 精品福利观看| 久久人妻福利社区极品人妻图片| 黄色怎么调成土黄色| 天天影视国产精品| 成人三级做爰电影| 国产高清激情床上av| 校园春色视频在线观看| 久久狼人影院| 国产精品永久免费网站| 99re在线观看精品视频| 大码成人一级视频| 日本五十路高清| 嫁个100分男人电影在线观看| 国产精品98久久久久久宅男小说| 一边摸一边抽搐一进一小说| 天堂俺去俺来也www色官网| 日韩高清综合在线| 国产99白浆流出| 欧美精品一区二区免费开放| 每晚都被弄得嗷嗷叫到高潮| 人人妻人人澡人人看| 岛国在线观看网站| 人成视频在线观看免费观看| 老熟妇乱子伦视频在线观看| 免费看十八禁软件| 欧美精品亚洲一区二区| 亚洲人成电影观看| 亚洲精品中文字幕在线视频| 日韩有码中文字幕| 国产av在哪里看| 精品无人区乱码1区二区| 午夜精品国产一区二区电影| 免费在线观看亚洲国产| 精品久久久久久,| 成年女人毛片免费观看观看9| 欧美日本亚洲视频在线播放| 老司机午夜福利在线观看视频| 丁香欧美五月| 午夜成年电影在线免费观看| 丰满的人妻完整版| 我的亚洲天堂| 视频区欧美日本亚洲| 亚洲中文日韩欧美视频| 黄色视频,在线免费观看| 好看av亚洲va欧美ⅴa在| 中文字幕另类日韩欧美亚洲嫩草| 99国产精品一区二区三区| 日日干狠狠操夜夜爽| www日本在线高清视频| 51午夜福利影视在线观看| 欧美日韩瑟瑟在线播放| 夜夜爽天天搞| 一级a爱视频在线免费观看| 免费看十八禁软件| 欧美亚洲日本最大视频资源| 在线观看日韩欧美| 女人高潮潮喷娇喘18禁视频| 亚洲成国产人片在线观看| 国产又爽黄色视频| 国产av一区在线观看免费| 国产精品免费一区二区三区在线| 怎么达到女性高潮| 色综合站精品国产| 亚洲专区字幕在线| 一个人观看的视频www高清免费观看 | 国产精品野战在线观看 | 欧美日韩精品网址| 亚洲五月婷婷丁香| 天天躁夜夜躁狠狠躁躁| 国产区一区二久久| 在线观看免费视频网站a站| 欧美中文综合在线视频| 午夜福利,免费看| 淫妇啪啪啪对白视频| www.999成人在线观看| 看免费av毛片| 99国产综合亚洲精品| 不卡av一区二区三区| 91av网站免费观看| 黄色成人免费大全| 麻豆国产av国片精品| 女人被躁到高潮嗷嗷叫费观| 久久这里只有精品19| 91成人精品电影| 老司机福利观看| 在线永久观看黄色视频| 国产在线精品亚洲第一网站| videosex国产| 又黄又粗又硬又大视频| 91av网站免费观看| 精品福利观看| 男人操女人黄网站| 新久久久久国产一级毛片| 99久久综合精品五月天人人| av视频免费观看在线观看| 国产在线精品亚洲第一网站| 视频区欧美日本亚洲| 久久精品亚洲精品国产色婷小说| 欧美成人性av电影在线观看| 九色亚洲精品在线播放| 制服诱惑二区| 亚洲av成人不卡在线观看播放网| 久久久国产成人精品二区 | 亚洲情色 制服丝袜| 19禁男女啪啪无遮挡网站| 色综合婷婷激情| 免费观看人在逋| 欧美国产精品va在线观看不卡| 国产精品久久久久久人妻精品电影| 在线观看免费视频日本深夜| 如日韩欧美国产精品一区二区三区| 嫩草影视91久久| 亚洲一区高清亚洲精品| 少妇被粗大的猛进出69影院| 9191精品国产免费久久| 久久久精品欧美日韩精品| 色播在线永久视频| 高清黄色对白视频在线免费看| 黄色成人免费大全| 99在线视频只有这里精品首页| av天堂在线播放| 少妇的丰满在线观看| 午夜福利一区二区在线看| 天堂影院成人在线观看| 一a级毛片在线观看| av片东京热男人的天堂| 亚洲 欧美 日韩 在线 免费| 亚洲七黄色美女视频| 丝袜美足系列| 欧美激情极品国产一区二区三区| 欧美另类亚洲清纯唯美| 亚洲伊人色综图| 老司机午夜十八禁免费视频| 日韩精品青青久久久久久| 免费看十八禁软件| 欧美不卡视频在线免费观看 | 久久精品国产亚洲av香蕉五月| 十分钟在线观看高清视频www| 久久人人97超碰香蕉20202| 99riav亚洲国产免费| 亚洲av五月六月丁香网| 日韩欧美免费精品| 日日摸夜夜添夜夜添小说| 欧美黑人欧美精品刺激| 国产精品永久免费网站| 久久亚洲精品不卡| 欧美一级毛片孕妇| 午夜影院日韩av| 桃色一区二区三区在线观看| 大型黄色视频在线免费观看| 波多野结衣一区麻豆| 久久精品91蜜桃| 男人的好看免费观看在线视频 | 啦啦啦 在线观看视频| 国产熟女xx| 一级片免费观看大全| 国产无遮挡羞羞视频在线观看| 色婷婷av一区二区三区视频| 男人舔女人下体高潮全视频| 免费观看精品视频网站| 少妇裸体淫交视频免费看高清 | 午夜亚洲福利在线播放| 国产精品亚洲一级av第二区| 人人妻人人澡人人看| 国产不卡一卡二| 最近最新中文字幕大全电影3 | 乱人伦中国视频| 成人亚洲精品av一区二区 | 国产伦一二天堂av在线观看| 免费看十八禁软件| www.自偷自拍.com| 亚洲国产中文字幕在线视频| 欧美黑人欧美精品刺激| 老司机深夜福利视频在线观看| 色精品久久人妻99蜜桃| ponron亚洲| 成人亚洲精品av一区二区 | 亚洲一区二区三区不卡视频| 亚洲国产欧美一区二区综合| 自线自在国产av| 久久久久国产一级毛片高清牌| 深夜精品福利| 亚洲一区二区三区不卡视频| 国产成人影院久久av| 久久人人97超碰香蕉20202| 老司机福利观看| 国产一区二区三区综合在线观看| 免费在线观看视频国产中文字幕亚洲| 大香蕉久久成人网| 午夜影院日韩av| 亚洲欧洲精品一区二区精品久久久| 91成年电影在线观看| 精品福利观看| 三级毛片av免费| 在线观看免费视频网站a站| 亚洲精品av麻豆狂野| 日日干狠狠操夜夜爽| 老司机靠b影院| 最新在线观看一区二区三区| 亚洲av美国av| 中文字幕精品免费在线观看视频| tocl精华| 久久国产乱子伦精品免费另类| 国产高清国产精品国产三级| 亚洲av日韩精品久久久久久密| 成年版毛片免费区| 美女 人体艺术 gogo| 亚洲国产看品久久| 欧美最黄视频在线播放免费 | 中出人妻视频一区二区| 午夜精品国产一区二区电影| 国产在线观看jvid| 无人区码免费观看不卡| 亚洲 国产 在线| 国产欧美日韩综合在线一区二区| 999久久久国产精品视频| 很黄的视频免费| 国产精品久久久av美女十八| 色尼玛亚洲综合影院| 免费在线观看完整版高清| 国产高清国产精品国产三级| 婷婷六月久久综合丁香| e午夜精品久久久久久久| 丝袜在线中文字幕| 午夜福利在线观看吧| 国产亚洲精品久久久久久毛片| 久久久久久久久久久久大奶| 九色亚洲精品在线播放| 18美女黄网站色大片免费观看| 又紧又爽又黄一区二区| 很黄的视频免费| 欧美黄色片欧美黄色片| 香蕉久久夜色| 国产欧美日韩一区二区精品| av在线播放免费不卡| 丝袜美足系列| 国产成人影院久久av| 亚洲国产精品999在线| 黄色怎么调成土黄色| 国产成人精品久久二区二区91| 少妇粗大呻吟视频| 法律面前人人平等表现在哪些方面| 午夜福利免费观看在线| 欧美成人午夜精品| 中出人妻视频一区二区| 黄色视频不卡| 亚洲五月色婷婷综合| 国产麻豆69| 深夜精品福利| 色老头精品视频在线观看| 久久久精品国产亚洲av高清涩受| 亚洲视频免费观看视频| 欧美 亚洲 国产 日韩一| 波多野结衣一区麻豆| 精品人妻在线不人妻| 中文字幕高清在线视频| 又大又爽又粗| 亚洲欧美日韩无卡精品| 丁香六月欧美| 视频在线观看一区二区三区| 国产成人系列免费观看| 99精品在免费线老司机午夜| 18美女黄网站色大片免费观看| 亚洲成人免费av在线播放| 国产激情欧美一区二区| av网站免费在线观看视频| 新久久久久国产一级毛片| 我的亚洲天堂| a级毛片在线看网站| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲欧美日韩高清在线视频| 国产精品 国内视频|