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

    黑河流域高寒草甸生態(tài)系統(tǒng)水分收支及蒸散發(fā)拆分研究

    2018-11-15 05:14:30童雅琴佩1李小雁1張賜成
    生態(tài)學(xué)報(bào) 2018年20期
    關(guān)鍵詞:貯水量黑河草甸

    童雅琴,王 佩1,,李小雁1,,*,張賜成,白 巖

    1 北京師范大學(xué)地表過程與資源生態(tài)國家重點(diǎn)實(shí)驗(yàn)室,北京 100875 2 北京師范大學(xué)地理科學(xué)學(xué)部自然資源學(xué)院,北京 100875

    根據(jù)IPCC第五次評估報(bào)告,自工業(yè)革命以來,全球氣候變化是當(dāng)前人類面臨的重要議題,其主要表現(xiàn)為氣溫升高和降水格局的變化[1],是改變水循環(huán)的重要驅(qū)動因素。水分收支是對水循環(huán)要素降水、蒸發(fā)蒸騰、徑流以及土壤儲水量變化等的定量刻畫[2]。在水分收支過程中,降水是陸地生態(tài)系統(tǒng)水分的主要來源,直接影響植被生產(chǎn)力;蒸散發(fā)(ET)包括土壤蒸發(fā)(E)和植被蒸騰(T)[3],是陸地生態(tài)系統(tǒng)水分的主要耗散形式,土壤蒸發(fā)和植被蒸騰受不同的生物物理過程控制[4],對于其在生態(tài)系統(tǒng)中的相對貢獻(xiàn)的定量研究對準(zhǔn)確預(yù)測生態(tài)系統(tǒng)功能對氣候變化的響應(yīng)至關(guān)重要,當(dāng)前關(guān)于生態(tài)系統(tǒng)蒸散發(fā)拆分的方法包括實(shí)地觀測法[5- 6],同位素法[7- 10]以及模型估算法[11- 13];土壤水分是土壤-植被-大氣連續(xù)體的關(guān)鍵因子[14- 15],土壤貯水量對維持“土壤水庫”功能,調(diào)節(jié)降水分布不均等具有重要作用[16]。高寒草甸/草原占中國西部高寒區(qū)面積的64%[17],被認(rèn)為是對氣候變化最為敏感的區(qū)域[18- 19]。在氣候變化背景下,準(zhǔn)確評估其水分收支過程對理解氣候、植被、水文交互作用至關(guān)重要。

    黑河流域是典型的干旱區(qū)內(nèi)陸河流域,是以水循環(huán)過程為紐帶的冰川/凍土、綠洲和荒漠為主要特征的多元生態(tài)系統(tǒng)[20],上游高寒區(qū)被認(rèn)為是干旱區(qū)水塔,水循環(huán)過程復(fù)雜。已有的研究多是從流域尺度運(yùn)用模型進(jìn)行水分收支估算,或者側(cè)重水分收支的個別要素[21- 24],如Jia等[25]用分布式水文模型分析了黑河中上游水分收支和能量收支平衡,表明在上游地區(qū)年均蒸散發(fā)和年均徑流分別占年均降水的63%和37%;Qin等[26]表明在黑河上游52.46%的出山徑流量產(chǎn)生于裸地景觀,上游森林景觀幾乎不產(chǎn)流;Yang等[27]運(yùn)用分布式生態(tài)水文模型研究表明,黑河上游地下水徑流和壤中流占主導(dǎo),地表徑流極少,且水分收支特征受植被類型的影響顯著;以上模型雖然能很好的分析流域尺度水分收支,但是,一般流域尺度上空間異質(zhì)性較強(qiáng),缺少在生態(tài)系統(tǒng)尺度上對水分收支過程的驗(yàn)證,且缺少直接觀測的生態(tài)系統(tǒng)實(shí)際蒸散發(fā),可能導(dǎo)致對生態(tài)系統(tǒng)水分以及植被需水量的估算存在偏差。渦動相關(guān)技術(shù)實(shí)地觀測的水汽通量數(shù)據(jù),為準(zhǔn)確的分析生態(tài)系統(tǒng)水分收支提供了可能,且當(dāng)前運(yùn)用雙源模型對蒸散發(fā)進(jìn)行拆分來評價(jià)生態(tài)系統(tǒng)水分收支狀況的研究較少。因此,本文以黑河流域高寒草甸生態(tài)系統(tǒng)為研究對象,運(yùn)用蒸散發(fā)拆分模型和渦動相關(guān)技術(shù)直接觀測的水汽通量數(shù)據(jù),估算2014和2015年黑河上游高寒草甸生態(tài)系統(tǒng)水分收支及其組分的季節(jié)變化和年際變化旨在定量評估黑河流域高寒草甸生態(tài)系統(tǒng)的水分收支狀況,為流域水資源承載力評估及水資源管理提供科學(xué)依據(jù)。

    1 研究區(qū)概況與數(shù)據(jù)來源

    1.1 研究區(qū)概況

    黑河是我國第二大內(nèi)陸河,發(fā)源于青海省祁連山區(qū)的冰川積雪帶,流域總面積14.3萬km2,是集合高山冰川、森林草原、平原綠洲及戈壁荒漠的復(fù)合生態(tài)系統(tǒng)[28]。本研究站點(diǎn)位于黑河上游青海省祁連縣中東地區(qū)的阿柔鄉(xiāng),黑河上游支流八寶河南側(cè)的河谷高地上,經(jīng)緯度為100.46°E,38.05°N,海拔3044 m,周圍地勢相對平坦開闊,自東南向西北略有傾斜下降,南北兩側(cè)約3km外是連綿的山丘和高山。研究區(qū)多年平均降水為400.3 mm[21],降水量年內(nèi)分配不均,主要集中在生長季(5—9月)。研究區(qū)屬于高寒草甸區(qū),植被類型以小嵩草草甸為主。

    1.2 數(shù)據(jù)來源

    研究數(shù)據(jù)包括2014和2015年自動氣象站數(shù)據(jù),渦動相關(guān)數(shù)據(jù)及相關(guān)的遙感數(shù)據(jù),本研究的氣象數(shù)據(jù)和渦度相關(guān)數(shù)據(jù)來源于阿柔凍融觀測站,黑河綜合遙感聯(lián)合試驗(yàn):水文氣象觀測網(wǎng)數(shù)據(jù)集(http://hiwater.westgis.ac.cn/)[29]。其中,渦動相關(guān)數(shù)據(jù)用來計(jì)算實(shí)際蒸散發(fā)和檢驗(yàn)?zāi)P偷男Ч?該數(shù)據(jù)是經(jīng)過Edire軟件后處理的30 min通量數(shù)據(jù),處理過程包括野點(diǎn)剔除、坐標(biāo)旋轉(zhuǎn)、頻率響應(yīng)修正、WPL修正以及初步質(zhì)量控制等,本文對缺失的渦動數(shù)據(jù)運(yùn)用查表法和平均晝夜法進(jìn)行插補(bǔ)[30]。氣象數(shù)據(jù)包括氣溫、降水、相對濕度、太陽輻射、風(fēng)速、土壤溫濕度(0—320 cm)等[31- 32]。其中,自動氣象站的降水采用翻斗式雨量計(jì)觀測,翻斗式雨量計(jì)安裝在阿柔超級站40 m觀測塔上。土壤溫濕度探頭埋設(shè)在地表0 cm和地下2、4、6、10、15、20、30、40、60、80、120、160、200、240、280、320 cm處,用于土壤水分和溫度的觀測。葉面積指數(shù)(LAI)數(shù)據(jù)來自2014和2015黑河生態(tài)水文遙感試驗(yàn):黑河流域1km分辨率5d合成葉面積指數(shù)(LAI)數(shù)據(jù)集(http://www.heihedata.org/data)[33- 35],模型中需要的日葉面積指數(shù)(LAI)采用三次樣條函數(shù)插值所得,植被高度數(shù)據(jù)參考黑河流域生態(tài)水文樣帶調(diào)查:2013年上游植被數(shù)據(jù)[36]。

    2 研究方法

    2.1 水分收支計(jì)算

    陸地生態(tài)系統(tǒng)水分收支的基本方程[37]為:

    P+Sd+I+R上入+R下入=ET+R上出+R下出+D+ΔW

    (1)

    式中,P為降水量,Sd為地下深層補(bǔ)給量,I為灌溉量,R上入、R下出分別為計(jì)算時(shí)段內(nèi)經(jīng)地面和地下流入生態(tài)系統(tǒng)內(nèi)的徑流量;ET(ET=E+T)為蒸散發(fā),包含有土壤蒸發(fā)E和植物蒸騰T,R上出、R下入分別為計(jì)算時(shí)段內(nèi)經(jīng)地面、地下流出的徑流量,ΔW為土壤貯水量的變化量,D為土壤水分的下滲。在祁連山高寒草甸沒有灌溉且地下水埋藏較深,可將I和Sd忽略,本文用于計(jì)算土壤貯水量的是0—320 cm處的土壤水分,所以可將D忽略,R上入、R下入和R上出、R下出可分別表示為R入和R出,所以水分收支公式可以簡化為:

    P+R入=ET+R出+ΔW

    (2)

    簡化后公式的物理意義可以理解為到達(dá)生態(tài)系統(tǒng)的水分消耗于土壤和植被的蒸散發(fā)、徑流以及土壤貯水量的變化,若P-ET+ΔW>R出-R入>0表示生態(tài)系統(tǒng)水分盈余;反之,表明生態(tài)系統(tǒng)水分虧缺。

    2.2 土壤貯水量及其變化的估算

    土壤貯水量(W)的計(jì)算采用各層土壤體積含水量的平均值與其代表的深度乘積累加求得[38],具體計(jì)算公式如下:

    (3)

    2.3 蒸散發(fā)組分估算

    本文使用Wang和Yamanaka雙源模型[39],該模型以植被冠層和土壤表面能量平衡為基本原理,包括太陽輻射傳輸過程,植被氣孔導(dǎo)度的模擬,而且其在不同的生態(tài)系統(tǒng)對蒸散發(fā)的模擬與分割都取得了很好的結(jié)果[39- 41]。模型的主要公式表征如下:

    (4)

    (5)

    式中,RnV是植被冠層的凈輻射(W/m2),HV是植被冠層的顯熱通量(W/m2),T是蒸騰量 (kg m-2s-1),fV為植被冠層的電介常數(shù),Sd是向下的短波輻射(W/m2),Ld是向下的長波輻射(W/m2),σ為 Stefan-Boltzmann常數(shù)(=5.67×10-8W m-2K-4),TG和TL分別為地表溫度和葉片表面溫度(℃),RnG為地表凈輻射(W/m2),G是土壤熱通量(W/m2),HG是地表顯熱通量(W/m2),E是蒸發(fā)通量(kg/ m2s),αG為地表反照率,αV是植冠層反照率,αG,αV可近似為常數(shù)0.1和0.2。

    能量平衡中的顯熱通量可以通過以下等式表示:

    HV=ρa(bǔ)(TL-Ta)/raV

    (6)

    HG=ρa(bǔ)(TG-Ta)/raG

    (7)

    式中,cp為干燥空氣的熱容量,raV為植被冠層的空氣動力阻抗,raG是地表的空氣動力阻抗(s/m),它們主要與風(fēng)速,植被高度和風(fēng)速測量儀的高度有關(guān),詳細(xì)描述可見王佩等[39]研究。

    模型結(jié)果中的土壤蒸發(fā)(E)和植被蒸騰計(jì)算公式如下:

    T=ρa(bǔ)([q(TL)-qa])/(raV+rc)

    (8)

    E=ρa(bǔ)([q(TG)-qa])/(raG+rSS)

    (9)

    式中,qsat(TL) 為葉片溫度的飽和比濕 (kg/kg),qsat(TG)為地表溫度的飽和比濕,qa為空氣比濕,raV是植被冠層的空氣動力阻抗,rc為氣孔導(dǎo)度(s/m),以及rSS是土壤表面阻抗(s/m)。

    rc=rst/LAI

    (10)

    式中,rst為葉片氣孔阻抗。

    (11)

    式中,rst_min和rst_max分別為土壤完全濕潤時(shí)葉片氣孔阻力的最大值和最小值,公式中的常數(shù)CSW可表示為:

    (12)

    式中,θ為土壤體積含水量(m3/m3),θmax是土壤最大體積含水量。土壤的表面阻力與土壤含水量滿足指數(shù)關(guān)系式:

    rsoil=a(θs/θ)b+c

    (13)

    式中,a,b,c為經(jīng)驗(yàn)系數(shù),a=3.5,b=2.3 和c=433.5。

    2.4 模型的敏感性分析

    為了評價(jià)模型指定參數(shù)和輸入變量可能存在的誤差,對模型進(jìn)行敏感性分析,敏感性系數(shù)si[39,42- 43]定義如下:

    (14)

    式中,Pi是指影響模型輸出結(jié)果(如ET或者T/ET)的i個驅(qū)動參數(shù)或變量,其偏導(dǎo)?O/(?Pi)計(jì)算公式如下:

    (15)

    3 結(jié)果分析

    3.1 模型的有效性和敏感性分析

    圖1為高寒草甸渦動相關(guān)系統(tǒng)觀測的和模型模擬的蒸散發(fā)以及地表溫度值的1∶1線性擬合關(guān)系,結(jié)果發(fā)現(xiàn)模型模擬值與實(shí)測值之間擬合度高,相關(guān)系數(shù)為分別0.90和0.89,表明該雙源模型對高寒草甸生態(tài)系統(tǒng)蒸散發(fā)的模擬是有效的。

    圖1 實(shí)測值與模型模擬蒸散發(fā)(潛熱)值與溫度對比2014—2015Fig.1 Comparison of evapotranspiration (latent heat flux) and temperature between measured and the predicted by model during 2014—2015

    表1總結(jié)了高寒草甸潛熱通量(LE)和蒸散比(T/ET)對模型輸入變量的敏感性。結(jié)果表明潛熱通量(LE)對相對濕度(RH)最為敏感,為-0.74 ± 0.32,即相對濕度5%的誤差會導(dǎo)致潛熱通量平均3.7%的誤差,對于蒸散比T/ET,葉面積指數(shù)(LAI)和土壤水分(θ)對其較為敏感,其5%的誤差分別可導(dǎo)致蒸散比T/ET2.2%和1.7%的誤差。因此可以看出模型估算的蒸散發(fā)(LE)和蒸散比(T/ET)存在的誤差都在可以忽略的范圍之內(nèi),表明模型模擬的結(jié)果較為可靠。

    表1模型輸入?yún)?shù)對蒸散發(fā)(ET)和蒸散比(T/ET)敏感性系數(shù)(Si)的平均值±標(biāo)準(zhǔn)偏差(SD)

    Table1Meanandstandarddeviation(SD)ofthesensitivitycoefficients(Si)ofevapotranspiration(ET)andtranspirationfraction(T/ET)fortheinputvariablesinmode

    變量名稱Variables 平均值±標(biāo)準(zhǔn)差Mean ± Standard deviation潛熱通量 Latent heat Fluex LE蒸散發(fā)比 Transpiration fraction T/ET氣溫 Temperature Ta/(℃)0.52±0.380.19±0.27相對濕度 Relative humid RH/%-0.74±0.320.00±0.22葉面積指數(shù) Leaf area index LAI/(m3/m3)0.16±0.090.44±0.28根系層土壤水分 Soil moisture θ/%0.22±0.060.34±0.28

    3.2 高寒草甸生態(tài)系統(tǒng)水分收支要素特征

    3.2.1 降水特征分析

    高寒草甸生態(tài)系統(tǒng)2014和2015年降水量為分別為520 mm和401 mm,降水量年內(nèi)分配不均,降水主要集中在5—9月(圖2),2014年和2015年生長季降水分別占全年總降水量的88.75%和90.21%,2014和2015年均溫度分別-0.43℃和0.28℃,氣溫、飽和氣壓差(VPD)以及潛在蒸散發(fā)(PET)的季節(jié)變化趨勢一致,呈現(xiàn)單峰變化規(guī)律,都是從5月開始上升,7月達(dá)到最大值,隨后開始緩慢下降。

    圖2 2014—2015年研究區(qū)氣象要素的季節(jié)變化Fig.2 Monthly variation of meteorological variables in study area from 2014 to 2015

    圖3 2014和2105年模擬蒸散發(fā)與觀測蒸散發(fā)對比Fig.3 Comparison of Predicted evapotranspiration and simulated evapotranspiration in 2014 and 2015

    3.2.2 蒸散發(fā)及其組分變化特征分析

    通過渦動相關(guān)觀測的高寒草甸生態(tài)系統(tǒng)2014年和2015年的實(shí)際蒸散發(fā)量分別為520 mm和551 mm,生長季(5—9)蒸散發(fā)分別為408 mm和422 mm,分別占全年蒸散發(fā)總量的78%和76%,圖3顯示為高寒草甸蒸散發(fā)的時(shí)間序列,從中可以看出觀測蒸散發(fā)、模型模擬蒸散發(fā)在時(shí)間變化上具有較好的一致性。日蒸散發(fā)量在全年尺度上呈現(xiàn)明顯的單峰變化規(guī)律,4月份蒸散發(fā)量開始增加,7月達(dá)到全年最大值,2014年和2015年月最大蒸散發(fā)量ET分別為85 mm和89 mm,日最大蒸散發(fā)量分別為5.30 mm/d和5.24 mm/d,隨后開始下降,這與生長季植被物候變化特征相一致。

    已有研究表明植被蒸騰(T)是陸地生態(tài)系統(tǒng)水分收支的重要組分,蒸散比(T/ET)可以很好的指示生態(tài)系統(tǒng)的生產(chǎn)力。通過模型對蒸散發(fā)的拆分發(fā)現(xiàn),在高寒草甸生態(tài)系統(tǒng)2014年和2015年蒸散比(T/ET)的日變化也呈現(xiàn)明顯的單峰變化規(guī)律(圖4),在非生長季,缺少植被覆蓋,蒸散發(fā)(ET)主要是土壤蒸發(fā),蒸散比(T/ET)很小0.3左右;在生長季(5—9月),隨著植被的生長,植被蓋度增加,植被蒸騰在蒸散發(fā)中占主導(dǎo),T/ET迅速增加,2014和2015年生長季平均蒸散比(T/ET)分別為0.74和0.79,日最大值可達(dá)0.91。

    3.2.3 土壤水分變化特征

    2014年和2015年高寒草甸生態(tài)系統(tǒng)0—320 cm土壤貯水量季節(jié)及年際變化特征如圖5所示,從圖中可以看出該高寒草甸生態(tài)系統(tǒng)土壤貯水量存在年內(nèi)波動,呈現(xiàn)雙峰變化規(guī)律。每年3月開始,土層解凍再加上積雪融化,土壤貯水量快速增加,到5—6月份達(dá)到每年的最大值,7—8月草甸生態(tài)系統(tǒng)具有較強(qiáng)的植被蒸騰,植被的生長消耗了土壤中儲存的水分,使得土壤貯水量略有下降,9—10月蒸騰減弱,土壤貯水量又有所上升,每年11月至次年3月份,受土層凍結(jié)的影響,實(shí)測土壤含水量僅為土壤中的液態(tài)含水量,固態(tài)水無法探測,該時(shí)段內(nèi)的土壤貯水量很小,即土壤貯水量受凍融作用影響顯著。圖6為不同深度凍融前后土壤體積含水量的變化,表層至10 cm處土壤體積含水量凍融前后變化顯著,消融后土壤體積含水量明顯高于凍結(jié)前,20 cm以下至120 cm處則相反,凍融前后土壤體積含水量變化不大,至160 cm處土壤無凍融,表明冬季凍結(jié)對高寒草甸生態(tài)系統(tǒng)表層土壤水分保持具有重要作用。

    圖4 2014和2015年高寒草甸蒸散比(T/ET)的日變化 Fig.4 The diurnal variation of transpiration fraction T/ET of Alpine meadow ecosystem in 2014 and 2015

    圖5 2014和2015年高寒草甸0—320 cm土壤貯水量的日變化 Fig.5 Diurnal variation of soil water storage (0—320 cm) of Alpine meadow in 2014 and 2015

    圖6 高寒草甸生態(tài)系統(tǒng)不同土層深度凍結(jié)前消融后土壤體積含水量Fig.6 Compare to the soil water content between the before freezing and after thawing at different soil depth

    在垂直變化上,高寒草甸生態(tài)系統(tǒng)土壤體積含水量在0—40 cm處波動較大(圖7),越往表層,土壤含水量越高,40 cm以下土壤體積含水量小(10%以下)且變化不大,說明0—40 cm為高寒草甸生態(tài)系統(tǒng)主要的貯水層。從年際變化上看,2014年該高寒草甸生態(tài)系0—320 cm處的土壤貯水量大于2015年,表明土壤貯水量的變化深受當(dāng)年降水量的影響。

    圖7 不同深度處的土壤體積含水量的變化Fig.7 Variation of volumetric soil water content at different depths

    3.3 高寒草甸生態(tài)系統(tǒng)水分收支特征

    圖8為2014和2015年高寒草甸生態(tài)系統(tǒng)水分收支要素降水(P),蒸散發(fā)(ET)和土壤貯水量變化(ΔW)的月變化及水分收支狀況,每年10月至次年2月,觀測到的降水量與蒸散發(fā)量很小,土壤貯水量的變化由于冬季土層凍結(jié),使(P-ET-ΔW)為正值;在每年3—5月,春季氣溫回升,再加上降水少,該時(shí)段凍結(jié)的土層及冬季降雪相繼融化,補(bǔ)充土壤水分,土壤儲水量急劇增加,導(dǎo)致(P-ET-ΔW)為負(fù)值,可見,水分收支的季節(jié)變化受凍融過程影響顯著;在完全不受凍融影響的月份6—9月(表2),2014年降水較多,生態(tài)系統(tǒng)有地表徑流形成,6—9月地表徑流量分別為6、15、0.45和22 mm;2015年降水量少,6—8月份,蒸散發(fā)不僅消耗全部的降水,還消耗儲存在土壤中的水分,生態(tài)系統(tǒng)呈現(xiàn)水分虧缺,6—8月虧缺量分別為22、12、13 mm,9月份到植被生長的末期,植被蒸騰減弱,該時(shí)段生態(tài)系統(tǒng)產(chǎn)流量為7 mm。

    圖8 2014和2015高寒草甸生態(tài)系統(tǒng)水分收支月變化Fig.8 Monthly variation of water budget in Alpine meadow ecosystem in 2014 and 2015

    從表3可以看到,在全年尺度上,2014年和2015年高寒草甸呈現(xiàn)不同程度的水分虧缺,2014年降水多,虧缺量為18 mm,2015年虧缺量為134 mm;在完全不受凍融影響的月份,2014年該高寒草甸生態(tài)系統(tǒng)產(chǎn)生地表徑流量為42 mm;2015年降水量少,生態(tài)系統(tǒng)呈現(xiàn)為水分虧缺,總的虧缺量為26 mm。

    表2 阿柔高寒草甸2014至2015年凍融起止時(shí)間

    根據(jù)不同深度土壤最高溫度(Tmax)和最低溫度(Tmin),秋季始凍期(Tmax>0℃且Tmin<0℃)、完全凍結(jié)期(Tmax<0 ℃且Tmin<0℃)和春季解凍期(Tmax>0℃且Tmin<0℃),完全解凍(Tmax>0℃且Tmin>0℃)。

    表3 2014和2015高寒草甸生態(tài)系統(tǒng)水分收支要素年變化

    4 討論

    4.1 蒸散比T/ET及影響因素分析

    圖9 葉面積指數(shù)(LAI)與蒸散比(T/ET)的關(guān)系 Fig.9 Relationship of LAI with daily mean transpiration fraction (T/ET)

    通過對蒸散發(fā)ET的拆分發(fā)現(xiàn)在高寒草甸生態(tài)系統(tǒng)植被蒸騰,尤其生長季植被蒸騰是生態(tài)系統(tǒng)水分消耗的主要形式,2014和2015年高寒草甸生態(tài)系統(tǒng)年均蒸散比(T/ET)分別為0.54和0.53。Wei等[44]研究表明在全球尺度上T/ET為0.57 ± 0.06。Hu等[4]基于S-W雙源模型研究的2004年高寒草甸的T/ET為0.44,王海波等[45]基于Penman-Monteith模型研究發(fā)現(xiàn)2008—2009年高寒草甸T/ET約為0.6。本研究估算T/ET與其他相關(guān)研究基本接近,但也存在一定的差異,可能與LAI數(shù)據(jù)集有關(guān),已有研究表明T/ET在季節(jié)尺度上主要受植被葉面積指數(shù)(LAI)的影響[4,39,44,46- 48],這與我們的研究結(jié)果相一致(圖9),葉面積指數(shù)LAI與T/ET呈現(xiàn)明顯的對數(shù)關(guān)系。Wei等[44]表明,不同的來源的LAI數(shù)據(jù)集,可能會導(dǎo)致T/ET存在一定的誤差,因此,高精度的LAI數(shù)據(jù)對蒸散發(fā)ET拆分影響重大。

    4.2 水分收支估算

    本研究利用水量平衡方程對高寒草甸生態(tài)系統(tǒng)2014和2015年水分收支進(jìn)行估算。在每年的非生長季,各月蒸散發(fā)量略大于降水量,生態(tài)系統(tǒng)水分虧缺,這與Yao等[49]運(yùn)用2004—2011實(shí)測數(shù)據(jù)對海北高寒草甸生態(tài)系統(tǒng)水分收支研究結(jié)果相一致。在2014年不受凍融影響的生長季,生態(tài)系統(tǒng)水分盈余。在正常年份的生長季,生態(tài)系統(tǒng)水分呈現(xiàn)虧缺,蒸發(fā)和蒸騰不僅消耗全部的降水,還消耗儲存在土壤中的水分,儲存在土壤中的水分對生態(tài)系統(tǒng)水分收支起重要調(diào)節(jié)作用。Yan等[50]通過對森林生態(tài)系統(tǒng)的研究發(fā)現(xiàn),土壤水分是旱季蒸散發(fā)的重要水源,這與本研究結(jié)果較為一致。在年尺度上,阿柔站高寒草甸生態(tài)系統(tǒng)并無地表徑流產(chǎn)生,只有在降水量大的濕潤年份的生長季,有地表徑流形成。Yang等[51]運(yùn)用分布式水文模型對黑河上游水量平衡研究表明2001—2012年,地下水徑流占總徑流量的80%左右,壤中流占20%左右,地表徑流量極小,占徑流量的2%—3%。陳仁升等[17]在黑上游高寒地區(qū)通過長時(shí)間小流域觀測實(shí)驗(yàn)表明,高山寒漠帶為山區(qū)流域的主要產(chǎn)流區(qū),而高寒草甸/草原區(qū)徑流貢獻(xiàn)很少,其水源涵養(yǎng)功能大于水文功能;通常平緩地形下的高寒草甸/草原具有較多的水源涵養(yǎng)功能,對徑流的貢獻(xiàn)少[17],這與本研究結(jié)果相當(dāng)。

    在高寒草甸生態(tài)系統(tǒng),土壤水分的劇烈變化主要發(fā)生在0—40 cm處,與郭淑海[45]在天山高寒草甸的試驗(yàn)結(jié)果一致,0—40 cm土壤貯水,是蒸散發(fā)消耗的重要水源,對維持植被生長、發(fā)育具有重要作用。此外,在生態(tài)系統(tǒng)水分收支中,外界徑流補(bǔ)給是生態(tài)系統(tǒng)水分消耗的重要來源,而且已有研究[52- 53]表明高寒草甸生態(tài)系統(tǒng)地下水埋藏較深,且壤中流較為發(fā)育,已有相關(guān)研究[26,54-55]表明黑河上游,受凍土的影響,流域產(chǎn)流以壤中流為主,因此,高寒草甸生態(tài)系統(tǒng)虧缺的水分可能是壤中流補(bǔ)給。

    5 結(jié)論

    本文運(yùn)用雙源模型對黑河上游高寒草甸生態(tài)系統(tǒng)蒸散發(fā)組分進(jìn)行拆分及以生態(tài)系統(tǒng)水量平衡方程為原理,估算了黑河流域高寒草甸生態(tài)系統(tǒng)2014年和2015年水分收支狀況,得到以下結(jié)論:(1)生長季植被蒸騰是生態(tài)系統(tǒng)水分主要的消耗形式,在季節(jié)尺度上,蒸散比(T/ET)主要受葉面積指數(shù)(LAI)的控制;(2)0—40 cm土壤儲水對維持生長季旺期植被蒸騰具有重要的意義;(3)在濕潤年(2014年)高寒草甸生態(tài)系統(tǒng)水分收支基本平衡,且生長季有產(chǎn)流;在正常年份(2015年)生態(tài)系統(tǒng)呈現(xiàn)水分虧缺,虧缺的水分主要來自外界壤中流的補(bǔ)給??傮w而言,實(shí)測數(shù)據(jù)可以準(zhǔn)確的刻畫生態(tài)系統(tǒng)水分收支過程,但也存在一定的不確定性,如儀器的觀測精度和觀測年限等,因此,在未來的研究中應(yīng)加強(qiáng)對生態(tài)系統(tǒng)的連續(xù)觀測,為遙感數(shù)據(jù)和模型模擬研究提供更精確地驗(yàn)證數(shù)據(jù)。

    猜你喜歡
    貯水量黑河草甸
    高山草甸
    The world's narrowest river
    山地草甸
    不同生育期調(diào)虧灌溉對景電灌區(qū)春玉米農(nóng)田土壤水分含量的影響
    綠肥壓青茬口麥田貯水量與紅蜘蛛、莖基腐病發(fā)生的相關(guān)性分析
    黑河的“護(hù)衛(wèi)隊(duì)”
    等高反坡階措施對滇中紅壤坡耕地土壤貯水量的影響
    到張掖看黑河
    文學(xué)港(2019年5期)2019-05-24 14:19:42
    武功山山地草甸的成因調(diào)查
    活力(2019年21期)2019-04-01 12:17:10
    張掖黑河濕地國家級自然保護(hù)區(qū)
    国产一区二区在线观看日韩| 在线免费观看不下载黄p国产| 国产单亲对白刺激| 久久久久免费精品人妻一区二区| 久久欧美精品欧美久久欧美| 一级毛片电影观看 | 波多野结衣巨乳人妻| av福利片在线观看| 久久亚洲精品不卡| videossex国产| 国产精品国产三级国产专区5o | 中国国产av一级| 精品欧美国产一区二区三| 久久久久久久国产电影| 日韩国内少妇激情av| 亚洲精品亚洲一区二区| 国产视频首页在线观看| 1024手机看黄色片| 可以在线观看毛片的网站| 在线天堂最新版资源| 在线观看一区二区三区| 综合色丁香网| 中文字幕精品亚洲无线码一区| 国产免费又黄又爽又色| 你懂的网址亚洲精品在线观看 | 国产探花在线观看一区二区| 亚洲精品aⅴ在线观看| 久久久久久伊人网av| 成人午夜高清在线视频| 一级av片app| 97超碰精品成人国产| 看十八女毛片水多多多| 看非洲黑人一级黄片| 国产精品女同一区二区软件| 男人狂女人下面高潮的视频| 成人高潮视频无遮挡免费网站| 亚洲精品,欧美精品| 国产精品日韩av在线免费观看| 午夜精品在线福利| 日韩一区二区视频免费看| 日本免费在线观看一区| 美女国产视频在线观看| 日本一二三区视频观看| 国产免费又黄又爽又色| 伦理电影大哥的女人| 国产精品不卡视频一区二区| 毛片一级片免费看久久久久| 亚洲欧美日韩卡通动漫| 99九九线精品视频在线观看视频| 乱人视频在线观看| 欧美zozozo另类| 国产亚洲精品久久久com| 秋霞在线观看毛片| 亚洲人成网站在线播| 国产成人精品一,二区| 欧美高清成人免费视频www| 99久久九九国产精品国产免费| 日本免费在线观看一区| 汤姆久久久久久久影院中文字幕 | 国产成人a区在线观看| av天堂中文字幕网| 国产欧美另类精品又又久久亚洲欧美| 人妻少妇偷人精品九色| 久久亚洲国产成人精品v| 午夜日本视频在线| 中国国产av一级| 美女高潮的动态| 色5月婷婷丁香| 亚洲自拍偷在线| 校园人妻丝袜中文字幕| 别揉我奶头 嗯啊视频| .国产精品久久| 99久久精品热视频| 午夜福利视频1000在线观看| 久久精品国产99精品国产亚洲性色| 少妇猛男粗大的猛烈进出视频 | 最近2019中文字幕mv第一页| 午夜视频国产福利| 日本色播在线视频| 亚洲欧美清纯卡通| 麻豆一二三区av精品| 欧美又色又爽又黄视频| 99国产精品一区二区蜜桃av| 最近最新中文字幕大全电影3| 成人亚洲欧美一区二区av| 在线免费观看的www视频| 大香蕉97超碰在线| 国产激情偷乱视频一区二区| 成人三级黄色视频| 最后的刺客免费高清国语| 男女那种视频在线观看| ponron亚洲| 国产欧美日韩精品一区二区| 舔av片在线| 一个人免费在线观看电影| 综合色丁香网| 日本爱情动作片www.在线观看| 少妇猛男粗大的猛烈进出视频 | 国产精品综合久久久久久久免费| 男女下面进入的视频免费午夜| 国产成人福利小说| 欧美激情在线99| 简卡轻食公司| 九九爱精品视频在线观看| 99久久精品一区二区三区| 你懂的网址亚洲精品在线观看 | 日本三级黄在线观看| 高清午夜精品一区二区三区| 国产成年人精品一区二区| 午夜爱爱视频在线播放| 看黄色毛片网站| 久久99精品国语久久久| 日日摸夜夜添夜夜爱| 国产在线一区二区三区精 | 99久久无色码亚洲精品果冻| av在线老鸭窝| 免费看a级黄色片| 一夜夜www| 蜜桃久久精品国产亚洲av| 最近中文字幕2019免费版| 精品午夜福利在线看| 水蜜桃什么品种好| 精品酒店卫生间| 99久久精品热视频| 乱系列少妇在线播放| 久久久久免费精品人妻一区二区| 久久久久久伊人网av| 国产老妇伦熟女老妇高清| 人人妻人人看人人澡| 亚洲欧美精品专区久久| 性色avwww在线观看| 免费观看精品视频网站| 欧美日韩精品成人综合77777| 丝袜美腿在线中文| 精品国内亚洲2022精品成人| 亚洲精品久久久久久婷婷小说 | 97超碰精品成人国产| 亚洲精品日韩av片在线观看| 亚洲四区av| 少妇的逼水好多| 精品人妻视频免费看| av卡一久久| 国产成人精品一,二区| 婷婷色麻豆天堂久久 | 村上凉子中文字幕在线| 男女边吃奶边做爰视频| 久久精品国产99精品国产亚洲性色| 亚洲av二区三区四区| 国产欧美另类精品又又久久亚洲欧美| 久久99热6这里只有精品| 国内精品一区二区在线观看| 免费电影在线观看免费观看| 麻豆精品久久久久久蜜桃| 久久久久国产网址| 亚洲av福利一区| 狂野欧美白嫩少妇大欣赏| 日韩成人伦理影院| 久久99热这里只有精品18| 99在线视频只有这里精品首页| 免费一级毛片在线播放高清视频| 日韩精品有码人妻一区| 人妻少妇偷人精品九色| 老司机影院毛片| 亚洲精品乱久久久久久| 久久精品国产99精品国产亚洲性色| h日本视频在线播放| 国产亚洲91精品色在线| 午夜精品国产一区二区电影 | 日本一本二区三区精品| 日韩欧美精品v在线| 美女内射精品一级片tv| 久久久久免费精品人妻一区二区| 久久精品国产自在天天线| 一级毛片电影观看 | 超碰97精品在线观看| 午夜精品一区二区三区免费看| 国产伦一二天堂av在线观看| 九九在线视频观看精品| 色哟哟·www| 老司机福利观看| 全区人妻精品视频| 国产精品女同一区二区软件| 国产色婷婷99| 看非洲黑人一级黄片| 国产精品嫩草影院av在线观看| 亚洲成av人片在线播放无| 免费播放大片免费观看视频在线观看 | 干丝袜人妻中文字幕| 黄色日韩在线| 国产精品精品国产色婷婷| 亚洲精品一区蜜桃| 亚洲欧美日韩无卡精品| 国产美女午夜福利| 最近中文字幕2019免费版| 天堂√8在线中文| 国产精品熟女久久久久浪| 国产亚洲最大av| 高清日韩中文字幕在线| 嫩草影院精品99| 一本久久精品| 日本三级黄在线观看| 网址你懂的国产日韩在线| 26uuu在线亚洲综合色| 亚洲av熟女| 亚洲精品456在线播放app| 久久久久免费精品人妻一区二区| 日产精品乱码卡一卡2卡三| 久久久精品欧美日韩精品| 免费黄色在线免费观看| 亚洲av中文字字幕乱码综合| ponron亚洲| or卡值多少钱| 免费看av在线观看网站| 免费看光身美女| 亚洲精品日韩av片在线观看| 国产黄色视频一区二区在线观看 | 嫩草影院新地址| 亚洲人成网站高清观看| 久久久久精品久久久久真实原创| 能在线免费观看的黄片| 亚洲国产高清在线一区二区三| 免费观看人在逋| 欧美一级a爱片免费观看看| 只有这里有精品99| 能在线免费看毛片的网站| 久久亚洲国产成人精品v| 日韩大片免费观看网站 | 国产精品综合久久久久久久免费| 国语对白做爰xxxⅹ性视频网站| 日韩在线高清观看一区二区三区| 又粗又硬又长又爽又黄的视频| 成人亚洲精品av一区二区| ponron亚洲| 一区二区三区免费毛片| 亚洲精品456在线播放app| 亚洲成av人片在线播放无| 亚洲av.av天堂| 久99久视频精品免费| 亚洲av二区三区四区| 我要看日韩黄色一级片| 99热精品在线国产| 国产精品国产三级国产av玫瑰| 国产精品伦人一区二区| 2022亚洲国产成人精品| 久久久色成人| 欧美日韩综合久久久久久| 噜噜噜噜噜久久久久久91| 乱码一卡2卡4卡精品| 亚洲国产精品sss在线观看| 亚洲欧美中文字幕日韩二区| 插逼视频在线观看| 好男人视频免费观看在线| 中文资源天堂在线| 国产一区二区亚洲精品在线观看| 国产精品麻豆人妻色哟哟久久 | 一个人观看的视频www高清免费观看| 久久久久久久久中文| 七月丁香在线播放| 精品久久久久久成人av| 2022亚洲国产成人精品| 日韩 亚洲 欧美在线| 在线天堂最新版资源| 女人十人毛片免费观看3o分钟| 国产激情偷乱视频一区二区| 国产亚洲91精品色在线| 国产片特级美女逼逼视频| h日本视频在线播放| 亚洲av电影在线观看一区二区三区 | 毛片一级片免费看久久久久| 久久久久国产网址| 欧美性感艳星| 18+在线观看网站| 国产精品一区二区三区四区久久| 亚洲不卡免费看| 日韩欧美精品v在线| 能在线免费观看的黄片| 国产女主播在线喷水免费视频网站 | 99久国产av精品国产电影| 色综合色国产| 亚洲av成人精品一二三区| 18禁裸乳无遮挡免费网站照片| 麻豆av噜噜一区二区三区| 国产又色又爽无遮挡免| 人妻少妇偷人精品九色| 欧美激情国产日韩精品一区| 国产不卡一卡二| 亚洲综合精品二区| 国产精品电影一区二区三区| 只有这里有精品99| 亚洲欧美日韩卡通动漫| 国产精品麻豆人妻色哟哟久久 | 中文欧美无线码| 18禁在线无遮挡免费观看视频| av在线老鸭窝| 天天躁夜夜躁狠狠久久av| av在线天堂中文字幕| 国产精品久久电影中文字幕| 高清午夜精品一区二区三区| 97在线视频观看| 特大巨黑吊av在线直播| 九九在线视频观看精品| 亚洲国产欧美人成| 欧美色视频一区免费| 青春草亚洲视频在线观看| 丰满少妇做爰视频| 国产色爽女视频免费观看| 精品国产一区二区三区久久久樱花 | 亚洲精华国产精华液的使用体验| 欧美日韩精品成人综合77777| 日韩高清综合在线| 少妇的逼水好多| 男人舔女人下体高潮全视频| 午夜精品在线福利| 久久亚洲国产成人精品v| 亚洲欧美成人综合另类久久久 | 1000部很黄的大片| 亚洲av电影在线观看一区二区三区 | 99热网站在线观看| 亚洲美女搞黄在线观看| 婷婷六月久久综合丁香| 亚洲精品一区蜜桃| 午夜a级毛片| 亚洲精品456在线播放app| 国产亚洲5aaaaa淫片| 午夜激情欧美在线| 日本五十路高清| 亚洲成人久久爱视频| 黄色日韩在线| 一级毛片我不卡| 欧美激情国产日韩精品一区| 国产伦精品一区二区三区视频9| www日本黄色视频网| 亚洲乱码一区二区免费版| 中文字幕熟女人妻在线| 51国产日韩欧美| 日韩国内少妇激情av| 国产精品,欧美在线| 在线免费观看不下载黄p国产| 大香蕉久久网| 天堂√8在线中文| 美女高潮的动态| 老司机影院毛片| 欧美精品国产亚洲| 色视频www国产| 国产高清不卡午夜福利| 成年女人永久免费观看视频| 麻豆精品久久久久久蜜桃| 国语自产精品视频在线第100页| 天堂√8在线中文| 国产91av在线免费观看| 国产不卡一卡二| 国产精品久久久久久久久免| 亚洲国产精品sss在线观看| 亚洲一区高清亚洲精品| 婷婷色综合大香蕉| 亚洲精品乱码久久久v下载方式| www日本黄色视频网| 如何舔出高潮| 草草在线视频免费看| 欧美另类亚洲清纯唯美| 国产精品av视频在线免费观看| 成人亚洲精品av一区二区| eeuss影院久久| 晚上一个人看的免费电影| 哪个播放器可以免费观看大片| 干丝袜人妻中文字幕| 亚洲欧美日韩无卡精品| 久久久久久久国产电影| 色综合亚洲欧美另类图片| 亚洲图色成人| 国产私拍福利视频在线观看| 亚洲在线观看片| 欧美精品国产亚洲| 波多野结衣巨乳人妻| 综合色av麻豆| 国产成人一区二区在线| 久久6这里有精品| 国产精品乱码一区二三区的特点| 卡戴珊不雅视频在线播放| 网址你懂的国产日韩在线| 亚洲高清免费不卡视频| 日本黄色视频三级网站网址| 精品一区二区三区视频在线| 岛国毛片在线播放| 国产免费视频播放在线视频 | 91久久精品电影网| 欧美xxxx性猛交bbbb| 久久热精品热| 国产69精品久久久久777片| av播播在线观看一区| 桃色一区二区三区在线观看| 五月伊人婷婷丁香| 在线播放国产精品三级| 亚州av有码| 免费观看精品视频网站| 狂野欧美激情性xxxx在线观看| 国产又色又爽无遮挡免| 久久99热这里只频精品6学生 | 青春草亚洲视频在线观看| 伦精品一区二区三区| 日本免费a在线| 免费看光身美女| 丝袜美腿在线中文| 欧美区成人在线视频| 禁无遮挡网站| 国内少妇人妻偷人精品xxx网站| 在线免费观看不下载黄p国产| 久久99精品国语久久久| 欧美另类亚洲清纯唯美| 亚洲一级一片aⅴ在线观看| 国产高清三级在线| 日本猛色少妇xxxxx猛交久久| 天堂av国产一区二区熟女人妻| 日韩强制内射视频| 精品少妇黑人巨大在线播放 | 亚洲第一区二区三区不卡| 99在线视频只有这里精品首页| 纵有疾风起免费观看全集完整版 | 国产精品久久久久久久久免| 亚洲欧美成人综合另类久久久 | 亚洲最大成人av| 中文资源天堂在线| 少妇猛男粗大的猛烈进出视频 | 日韩,欧美,国产一区二区三区 | 国产精品综合久久久久久久免费| 久久久久性生活片| 成人三级黄色视频| 亚洲国产精品国产精品| 国产精品一区二区三区四区久久| 赤兔流量卡办理| 深爱激情五月婷婷| .国产精品久久| 毛片女人毛片| 午夜福利在线在线| 特大巨黑吊av在线直播| 美女黄网站色视频| 看黄色毛片网站| 日韩成人av中文字幕在线观看| 亚洲色图av天堂| 丰满少妇做爰视频| 国产精品爽爽va在线观看网站| 中文资源天堂在线| 亚洲国产精品成人综合色| 亚洲av二区三区四区| 男插女下体视频免费在线播放| 身体一侧抽搐| 中文字幕精品亚洲无线码一区| 欧美日韩在线观看h| 亚洲美女搞黄在线观看| 长腿黑丝高跟| 欧美日韩国产亚洲二区| 亚洲高清免费不卡视频| 国产免费一级a男人的天堂| 99久久九九国产精品国产免费| 桃色一区二区三区在线观看| 蜜臀久久99精品久久宅男| 特大巨黑吊av在线直播| 免费不卡的大黄色大毛片视频在线观看 | 又黄又爽又刺激的免费视频.| 一边亲一边摸免费视频| 亚洲,欧美,日韩| 一个人观看的视频www高清免费观看| 日韩,欧美,国产一区二区三区 | 日日摸夜夜添夜夜添av毛片| 亚洲熟妇中文字幕五十中出| 国产av一区在线观看免费| 午夜日本视频在线| 亚洲四区av| 高清午夜精品一区二区三区| 亚洲五月天丁香| 在线a可以看的网站| 美女cb高潮喷水在线观看| 中文天堂在线官网| 国产v大片淫在线免费观看| 亚洲伊人久久精品综合 | 搡女人真爽免费视频火全软件| 久久久久久久亚洲中文字幕| 国产精品爽爽va在线观看网站| 黄色配什么色好看| 2021少妇久久久久久久久久久| 波多野结衣高清无吗| 欧美极品一区二区三区四区| 97超视频在线观看视频| 97超碰精品成人国产| 国产免费视频播放在线视频 | 有码 亚洲区| 亚洲,欧美,日韩| 能在线免费看毛片的网站| 中文字幕av成人在线电影| 国产精品蜜桃在线观看| 久久久久久久久久黄片| 欧美日本视频| 日韩欧美精品免费久久| 男女啪啪激烈高潮av片| 精品99又大又爽又粗少妇毛片| 亚洲18禁久久av| 国产精品不卡视频一区二区| av国产久精品久网站免费入址| 纵有疾风起免费观看全集完整版 | 最近中文字幕高清免费大全6| 成人欧美大片| 国产精品,欧美在线| 中文字幕人妻熟人妻熟丝袜美| 精品久久久久久久久av| 日韩国内少妇激情av| 国产亚洲av嫩草精品影院| 神马国产精品三级电影在线观看| 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产欧美在线一区| 日韩国内少妇激情av| 国产亚洲av嫩草精品影院| 亚洲真实伦在线观看| 亚洲美女搞黄在线观看| 色综合站精品国产| 亚洲成色77777| 国产美女午夜福利| 亚洲美女视频黄频| 午夜视频国产福利| 亚洲婷婷狠狠爱综合网| 亚洲精品日韩在线中文字幕| 亚洲不卡免费看| 欧美一区二区国产精品久久精品| 国产精品一区二区三区四区免费观看| 男人和女人高潮做爰伦理| 亚洲精品,欧美精品| videos熟女内射| 激情 狠狠 欧美| 色网站视频免费| 最近的中文字幕免费完整| 免费播放大片免费观看视频在线观看 | 亚洲中文字幕日韩| 青春草国产在线视频| 我的女老师完整版在线观看| 91狼人影院| 免费观看的影片在线观看| 国产午夜精品久久久久久一区二区三区| 老女人水多毛片| 亚洲av一区综合| 在线天堂最新版资源| 久久精品久久久久久久性| 我的女老师完整版在线观看| 麻豆成人av视频| 亚洲av免费高清在线观看| 哪个播放器可以免费观看大片| 夜夜爽夜夜爽视频| 日日摸夜夜添夜夜爱| 国内少妇人妻偷人精品xxx网站| 精品久久久久久成人av| 男女国产视频网站| 美女脱内裤让男人舔精品视频| 久久99热6这里只有精品| 在线免费观看的www视频| 高清视频免费观看一区二区 | 免费黄色在线免费观看| 汤姆久久久久久久影院中文字幕 | 国产成人a∨麻豆精品| 成年女人看的毛片在线观看| 韩国高清视频一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 国产精品乱码一区二三区的特点| 简卡轻食公司| 精品一区二区免费观看| 色5月婷婷丁香| 一本一本综合久久| 久久精品久久久久久噜噜老黄 | 日韩av在线大香蕉| 久久久亚洲精品成人影院| 在线播放国产精品三级| 欧美一区二区亚洲| 三级国产精品欧美在线观看| 99热这里只有精品一区| 日韩强制内射视频| 免费观看在线日韩| 亚洲av电影不卡..在线观看| 97在线视频观看| 欧美区成人在线视频| 久久久国产成人精品二区| 国产私拍福利视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 成人av在线播放网站| 久久99热这里只有精品18| 久久久久免费精品人妻一区二区| 午夜精品在线福利| 午夜a级毛片| 亚洲精华国产精华液的使用体验| 国产精华一区二区三区| 日本午夜av视频| 精品久久久久久久人妻蜜臀av| videossex国产| 国产真实伦视频高清在线观看| 成人一区二区视频在线观看| 久久精品久久精品一区二区三区| 久久精品人妻少妇| 乱码一卡2卡4卡精品| 69人妻影院| 国产三级中文精品| 国产一区有黄有色的免费视频 | 99在线视频只有这里精品首页| 搡女人真爽免费视频火全软件| 久久久久免费精品人妻一区二区| 欧美色视频一区免费| 亚洲怡红院男人天堂| 欧美人与善性xxx| 水蜜桃什么品种好| 日韩大片免费观看网站 | 18禁在线播放成人免费| 1000部很黄的大片| 国产乱人偷精品视频| 校园人妻丝袜中文字幕| 亚洲激情五月婷婷啪啪| 亚洲人成网站高清观看| 男女视频在线观看网站免费| 国产精品国产三级专区第一集| 久久久精品大字幕| 久久午夜福利片|