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

    三江源區(qū)退化高寒草甸蒸散特征及凍融變化對其的影響

    2020-09-17 07:35:08田曉暉張立鋒陳之光唐艷鴻
    生態(tài)學(xué)報 2020年16期
    關(guān)鍵詞:源區(qū)草甸凍融

    田曉暉,張立鋒,張 翔,陳之光,趙 亮,李 奇,唐艷鴻,古 松,*

    1 南開大學(xué)生命科學(xué)學(xué)院,天津 300071 2 中南林業(yè)科技大學(xué)生命科學(xué)與技術(shù)學(xué)院,長沙 410004 3 中國科學(xué)院西北高原生物研究所,西寧 810008 4 北京大學(xué)城市與環(huán)境學(xué)院,北京 100871

    陸地生態(tài)系統(tǒng)蒸散(ET)是全球水量平衡的主要分量,也是生態(tài)系統(tǒng)水分消耗的主要方式,主要由土壤蒸發(fā)和植被蒸騰組成[1]。由于ET在土壤—植被—大氣系統(tǒng)中的重要紐帶作用,在全球變化科學(xué)研究中受到高度關(guān)注[2]。有研究表明,地表降水大約有70%通過蒸散返回大氣,在干旱區(qū)甚至達(dá)到95%以上[3]。ET變化受輻射、降水、溫度、植被等因子影響;反之,ET變化又影響生態(tài)系統(tǒng)植被生長狀況和水分收支等。水熱運動相互耦合,蒸散影響水分收支的同時,也會影響生態(tài)系統(tǒng)的能量(潛熱和顯熱能等)分配,從而對區(qū)域氣候產(chǎn)生一定影響[4]。草地是陸地生態(tài)系統(tǒng)中分布最廣泛的植被類型之一,約占陸地生態(tài)系統(tǒng)總面積的1/3[5]。然而,草地蒸散的年際之間存在較大的變化,近年來受到了科學(xué)界的廣泛關(guān)注[4,6]。因此準(zhǔn)確掌握草地生態(tài)系統(tǒng)的蒸散變化及環(huán)境因子對其的影響,對于更好地理解全球水分收支動態(tài)及氣候變化具有重要意義。

    青藏高原平均海拔超過4000 m,面積約250萬km2,被譽為“世界第三極”,由于其獨特的地理和氣候條件,成為全球氣候變化的敏感區(qū)和放大區(qū)[7]。草地是其最主要的植被,占青藏高原面積的70%,其中高寒草甸分布最廣,支撐著青藏高原地區(qū)的畜牧業(yè)發(fā)展[8]。由于青藏高原特殊的地理位置,高原上廣泛分布著多年凍土和季節(jié)凍土,近地表層土壤的季節(jié)凍融循環(huán)影響著高原地表能量和水分交換[9]。近幾十年來,在全球氣候變暖的背景下,高原的溫度呈明顯的上升趨勢,有研究表明,1960年以來,青藏高原氣溫以0.2—0.3℃/10 a的速度增加,而冬季增溫的幅度(0.3—0.5℃/10 a)遠(yuǎn)遠(yuǎn)高于全年平均增溫幅度[10]。不斷加劇的氣候變暖正改變凍土凍融發(fā)生過程和狀態(tài),總體上呈現(xiàn)最大凍土深度降低,凍結(jié)起始日期推遲,融化起始日期提前及凍結(jié)持續(xù)期縮短等的趨勢[11-12]。全球氣候變暖改變了青藏高原土壤凍融狀態(tài),這必將引起生態(tài)系統(tǒng)水分交換的變化。

    三江源區(qū)位于青藏高原腹地,是長江、黃河、瀾滄江的發(fā)源地,素有“中華水塔”之稱,在水源涵養(yǎng)、維持生物多樣性以及氣候調(diào)節(jié)等方面發(fā)揮著重要作用,成為我國乃至世界上影響力最大的生態(tài)調(diào)節(jié)區(qū)之一[13]。然而,在氣候變化和超載放牧的雙重壓力下,三江源區(qū)高寒草甸出現(xiàn)了嚴(yán)重的退化現(xiàn)象,部分地區(qū)形成了“黑土灘”,極大地降低了區(qū)域的水源涵養(yǎng)功能,進而嚴(yán)重影響到三江源區(qū)生態(tài)系統(tǒng)的結(jié)構(gòu)和功能,改變了區(qū)域的水分收支狀況,同時對三江源區(qū)的畜牧業(yè)乃至長江、黃河中下游地區(qū)的社會經(jīng)濟發(fā)展構(gòu)成巨大的威脅[13-14]。雖然關(guān)于青藏高原高寒草地的蒸散研究也有一些報道[6-7,15],但由于高原面積廣闊,地理環(huán)境差異較大,結(jié)果不盡相同。三江源區(qū)是我國重要的水源涵養(yǎng)地,也是典型的季節(jié)凍土區(qū),然而對該區(qū)域高寒草甸蒸散變化的詳盡研究仍然相對匱乏,更缺少凍融作用對其生態(tài)系統(tǒng)蒸散影響的相關(guān)報道。因此,定量研究三江源區(qū)退化高寒草甸蒸散變化及其凍融作用對其的影響對深入研究三江源區(qū)水分收支變化具有重要現(xiàn)實意義。

    渦度相關(guān)技術(shù)是通過測定和計算CO2和H2O等物理量脈動與垂直風(fēng)速脈動的協(xié)方差求算湍流通量的方法,目前已成為直接測定地表下墊面與大氣間氣體通量交換的理想方法,廣泛應(yīng)用于陸地生態(tài)系統(tǒng)的蒸散研究[1]。本研究利用渦度相關(guān)技術(shù)和微氣象系統(tǒng)對三江源區(qū)退化高寒草甸的水汽通量及相關(guān)環(huán)境因子進行了連續(xù)觀測,對2016年和2017年獲取的觀測數(shù)據(jù)進行統(tǒng)計分析,其主要目的是:(1)揭示三江源區(qū)退化高寒草甸生態(tài)系統(tǒng)蒸散的季節(jié)變化特征及年際差異;(2)探討分析土壤凍融交替變化對蒸散的影響;(3)闡明蒸散變化對環(huán)境因子的響應(yīng)。以期為全面探討三江源區(qū)退化草甸蒸散對氣候變化的響應(yīng)及其生態(tài)系統(tǒng)水分收支變化提供參考依據(jù)。

    1 研究地概況和研究方法

    1.1 研究地概況

    本試驗地位于青海省果洛州瑪沁縣大武鎮(zhèn)東南部的退化高寒草甸(34° 21′ N, 100° 29′ E),海拔3958 m。研究區(qū)地勢平坦,草地分布均勻,該研究站代表了三江源區(qū)高寒草地典型植被類型[16]。該區(qū)屬于典型的高原大陸性氣候,無四季之分,僅有冷暖季之別,冷季漫長,干燥而寒冷,暖季短暫,濕潤而涼爽;溫度年差較小而日差較大,研究地年平均氣溫為-0.2℃,最冷月1月平均氣溫為-12.3℃,最熱月7月的平均氣溫為10.1℃,全年無絕對無霜期。該區(qū)日照充足,太陽輻射強烈,年總輻射量為5500—6800 MJ/m2;年降水量為420—560 mm,85%集中在5—9月。該試驗地土壤類型以高山草甸土為主,建群種為矮蒿草(Kobresiahumilis),主要伴生種為小蒿草(Kobresiapygmaea)、垂穗披堿草(Elymusnutans)、早熟禾(Poaannua)、細(xì)葉亞菊(Ajaniatenuifolia)等[17]。根據(jù)之前對該研究地植被類型、蓋度等方面調(diào)查研究證實其處于退化狀態(tài)[18]。退化高寒草甸于4月底或5月初返青;地上生物量從5月開始增加,7月底至8月初達(dá)到最大,10月植物開始衰老。

    1.2 觀測方法

    渦度相關(guān)和微氣象觀測系統(tǒng)安裝在地勢平坦、視野開闊的退化高寒草甸上,“風(fēng)浪區(qū)”(半徑大于300 m)滿足渦度相關(guān)觀測要求。觀測要素主要包括輻射通量、顯熱和潛熱通量、土壤熱通量、空氣溫濕度、不同深度土壤溫度、不同深度土壤含水量、降水量等,觀測儀器及安裝高度見表1。渦度相關(guān)系統(tǒng)采樣頻率為10 Hz,所有觀測數(shù)據(jù)每15 min輸出一組平均值,并儲存在CR5000和CR23X (CSI, USA) 數(shù)據(jù)記錄儀。

    本研究于2016年6—8月和2017年6、7、9月每月中旬對植被地上部生物量(above-ground biomass, AGB)和葉面積指數(shù)(LAI)進行調(diào)查。生物量測定采取收割法,隨機選擇5個大小為50 cm×50 cm的樣方,將每一樣方內(nèi)綠色植物齊地面剪下后放進取樣袋內(nèi)并編號。之后將樣品帶回實驗室,放入65℃恒溫烘箱中72 h(烘干至恒重)稱重(g/m2)。同時選取5個大小為25 cm×25 cm的樣方,齊地面剪下地上部植物,用葉面積儀(Li- 3000, Li-Cor)測定植物葉面積。上述5個重復(fù)的平均值用于本研究的數(shù)據(jù)分析。

    1.3 數(shù)據(jù)分析

    首先對渦度相關(guān)系統(tǒng)觀測的原始數(shù)據(jù)進行坐標(biāo)軸旋轉(zhuǎn),溫度變換和WPL校正,此外,剔除由于降水、儀器故障等引起的錯誤數(shù)據(jù)以及摩擦風(fēng)較弱的數(shù)據(jù)(摩擦風(fēng)速u*<0.1 m/s)。由地表反照率的變化可知,該研究地偶有降雪發(fā)生(主要發(fā)生在11月至翌年3月),但降雪頻次和降雪量均較少且很快消融,由于蒸散量主要集中在生長季,因此對研究結(jié)果不會產(chǎn)生實質(zhì)性影響,本研究并未進行特別分析。對于缺失及剔除的異常數(shù)據(jù),按以下方法進行插補,當(dāng)缺失的通量數(shù)據(jù)時間小于2 h時,使用線性內(nèi)插法進行插補,而大于2 h的數(shù)據(jù)缺口則使用非線性回歸法進行插補[19-20]。最后采用Excel和SPSS等軟件進行統(tǒng)計分析。

    表1 觀測要素及其儀器

    圖1 退化高寒草甸渦度相關(guān)測定的顯熱與潛熱通量之和(H + LE)與能量平衡法測定的有效能(Rn - G)的關(guān)系 Fig.1 The relationship between the sum of sensible heat and latent heat fluxes (H + LE) by eddy covariance and the available energy (Rn - G) by energy balance method in degraded meadow

    生態(tài)系統(tǒng)能量平衡通常用下列公式來表示[20]:

    Rn-G=H+LE

    式中,Rn為凈輻射通量(W/m2);H為顯熱通量(W/m2);LE為潛熱通量(W/m2);G為土壤熱通量(W/m2)。

    能量閉合度是評價渦度相關(guān)技術(shù)的方法之一,通常用(H+LE)與(Rn-G)進行線性回歸,利用其斜率和截距來分析能量的閉合度[20]。由圖1可知,本研究能量閉合度約為0.71,說明存在能量不閉合現(xiàn)象,但該值在已報道的0.55-0.99范圍內(nèi)[21]。能量不閉合現(xiàn)象在渦度相關(guān)觀測中普遍存在,由于導(dǎo)致能量不閉合的原因很復(fù)雜,詳細(xì)原因有待于深入研究。

    波文比(β)定義為顯熱通量與潛熱通量的比值[22],即:

    β=H/LE

    冠層導(dǎo)度(gc)計算公式為[23]:

    1/gc=ρCPVPD/(γLE)+(βΔ/γ-1)/ga

    解耦系數(shù)(Ω)計算公式為[24]:

    Ω=(Δ+γ)/(Δ+γ(1+ga/gc))

    式中,ρ為空氣密度(kg/m3),Cp為空氣定壓比熱(MJ/kg℃),VPD為飽和水汽壓差(kPa),Δ飽和水汽壓—空氣溫度曲線的斜率,γ干濕表常數(shù),ga為空氣導(dǎo)度(mm/s),ga按下列公式計算[23]:

    1/ga=μ/μ*2+6.2μ*-0.67

    式中,u為風(fēng)速,u*為摩擦風(fēng)速(m/s)

    平衡蒸散(ETeq)由下式計算得出[25]:

    ETeq=Δ(Rn-G)/L(Δ+γ)

    式中,L為水的汽化潛熱系數(shù)(2.501 MJ/kg)。

    2 結(jié)果與分析

    2.1 環(huán)境因子季節(jié)變化

    2016與2017年的Rn季節(jié)變化趨勢相同(圖2),最高值分別為12.81(6月)和11.22 MJ m-2d-1(7月),最低值出現(xiàn)在12月左右。但各月的Rn值略有差異,生長季中,2016年生態(tài)系統(tǒng)接收的Rn略高于2017年,分別為1654.88 MJ/m2和1546.71 MJ/m2,均超過全年Rn的60%。

    兩年氣溫(Ta)和5 cm深度的土壤溫度(Ts)變化趨勢與Rn相似(圖2),最高值出現(xiàn)在7—8月,最低值均出現(xiàn)在1月。2016和2017年Ta的平均值分別為-0.35℃和-0.41℃,月均最高值分別出現(xiàn)在8月(11.0℃)和7月(9.7℃),而生長季中,2016年Ta的平均值為7.2℃,高于2017年平均值6.8℃。2016和2017年Ts的平均值分別為2.5℃和3.8℃,月均最高值分別出現(xiàn)在8月(12.9℃)和7月(13.4℃),而生長季中,2016年Ts的平均值為9.4℃,低于2017年的10.6℃??傮w而言,兩年的Ts高于Ta,且2016年的Ta略高于2017年,而2017年Ts高于2016年。

    圖2 2016和2017年退化高寒草甸的凈輻射(Rn)、降水量(P)、5 cm土壤含水量(SWC5)、空氣溫度(Ta)、5 cm土壤溫度(Ts)和飽和水汽壓差(VPD)年變化Fig.2 Annual variations of net radiation (Rn), precipitation (P), soil water content at 5 cm depth (SWC5), air temperature (Ta), soil temperature at 5 cm depth (Ts) and vapor pressure deficit (VPD) in degraded meadow for 2016 and 2017

    2016和2017年降水量相近(圖2),分別為451.8 mm和442.3 mm,且降水主要集中在生長季5—9月,前者為386.7 mm,略高于后者的368.4 mm,分別占全年降水量的85.6%和83.3%。然而,兩年生長季中月降水量分配存在較大年際差異,2016年降水量的最高值出現(xiàn)在5月(94.9 mm),其次為8月(94.7 mm),6、7、9月的降水量在60—70 mm之間;而2017年降水主要集中在8月(118.4 mm)和9月(87.7 mm),5—7月降水量相對較少,生長旺季的7月降水量僅為40.9 mm。

    5 cm深度土壤含水量(SWC5)的季節(jié)變化與降水量密切相關(guān)(圖2),最高值出現(xiàn)在降水較多的生長季,然而,由于該時期的ET較高,進而導(dǎo)致SWC5在降水相對較多的7—8月有所下降,類似結(jié)果在青藏高原其他草地研究中也有報道[5,28]。生長季中,2016和2017年SWC5的月最低值分別出現(xiàn)在生長旺季的8月(0.17 m3/m3)和7月(0.19 m3/m3),最高值分別出現(xiàn)在6月(0.28 m3/m3)和9月(0.31 m3/m3),且2016年生長季的SWC5平均值(0.24 m3/m3)略低于2017年(0.25 m3/m3)。總體而言,2016年5—7月的SWC5高于2017年同期,而8—9月卻低于2017年同期的SWC5,這與兩年的降水量變化基本一致。

    飽和水汽壓差(VPD)是表征空氣濕度的物理量,其變化受降水、輻射和氣溫等因子影響。由VPD的季節(jié)變化可知(圖2),總體上生長季的VPD高于非生長季,但各月的VPD值出現(xiàn)明顯的年際差異。2016年的月均最高值出現(xiàn)在8月(0.87 kPa),而2017年月均最高值出現(xiàn)在7月(0.86 kPa)。生長季中,2016年VPD的平均值為0.69 kPa,明顯高于2017的0.62 kPa,而非生長季中的平均VPD分別為0.48和0.42 kPa,總體上2016年VPD高于2017年。

    2.2 蒸散季節(jié)變化

    2016和2017年日蒸散量的年變化趨勢基本一致,最高值均出現(xiàn)在植物生長旺季的7—8月,最低值出現(xiàn)在土壤凍結(jié)的冬季(圖3)。土壤蒸發(fā)從3月末或4月初的凍融交替開始逐漸增加,之后隨著土壤的解凍、降水的增加和植物的生長,蒸散(ET)迅速增加。2016和2017年月蒸散總量最大值分別為71.8 mm(8月)和92.5 mm(7月),之后隨著輻射和溫度的降低、植物生長的結(jié)束,ET逐漸降低(圖3)。2016和2017年的ET分別為428.6 mm和485.6 mm,而生長季中的ET分別為314.4 mm和350.7 mm,分別占到全年總蒸散量的73%和72%,2017年蒸散量高于2016年,尤其是2017年7—10月ET明顯高于2016年同期。

    圖3 2016—2017年退化草甸日蒸散量和月蒸散量的年變化Fig.3 Annual variation of daily and monthly evapotranspiration (ET) in degraded meadow for 2016 and 2017

    2.3 土壤凍融與蒸散2.3.1 凍融時段的劃分

    參考Guo等[26]對凍融階段的劃分方法,本研究利用5 cm深土壤溫度日最高和最低值變化,將全年劃分為凍結(jié)期、凍融交替期、消融期3個時段:(1)土壤凍結(jié)期:土壤溫度日最高值小于0℃;(2)凍融交替期:土壤溫度日最高值大于0℃,而日最低值小于0℃;(3)土壤消融期:土壤溫度日最低值大于0℃。為避免隨機天氣過程對土壤凍融階段轉(zhuǎn)變的影響,只有當(dāng)連續(xù)三天滿足下一階段條件時,這三天中的第一天才作為下一階段起始日期。為了研究植被對蒸散的影響,進一步把5—9月作為植物生長季獨立劃分出來。由3個時段及生長季的劃分(表2)可知,2016年的凍結(jié)期和凍融交替期分別比2017年長8 d,而2017年的消融期比2016年多15 d。

    2.3.2不同凍融時段蒸散量和蒸散速率

    為闡明2016和2017年不同時段的ET差異,對上述劃分的3個時段及生長季(表2)的ET量和ET速率進行了統(tǒng)計(表3)。由表可知,2016年僅凍融交替期的ET量高于2017年,而其他3個時段均低于2017年,尤其在消融期和生長季,2017年比2016年分別增加了63.1 mm和36.3 mm。

    表2 2016和2017年土壤的凍結(jié)期、凍融交替期、消融期及生長季的劃分

    2016和2017年ET速率均在生長季達(dá)到最高值,分別為2.05 mm/d和2.29 mm/d;最低值出現(xiàn)在凍結(jié)期,分別為0.27 mm/d和0.33 mm/d。2016年凍融交替期的ET速率為0.97 mm/d,高于2017年的0.73 mm/d,且2016年凍融交替期的時間比2017年長8 d,因此該時段2016年ET量高于2017年,而在其他3個時段,2016年ET速率均低于2017年。

    表3 2016和2017年3個時段及生長季蒸散量和蒸散速率

    2.3.3不同凍融時段蒸散日變化

    2016和2017年各時段的蒸散均呈明顯的日變化(圖4),ET的日最高值出現(xiàn)在13:00 h左右,而夜間維持在零附近。生長季的蒸散速率明顯高于其他3個時段,而凍融交替期的蒸散速率介于消融期和凍結(jié)期之間。生長季期間,2017年蒸散速率的日最高值為0.30 mm/h,高于2016年同期的0.28 mm/h,而在土壤凍融交替時期,2017年蒸散速率的日最高值為0.11 mm/h,低于2016年同期的0.13 mm/h,在凍結(jié)期和消融期,2016年的蒸散速率均低于2017年。

    圖4 2016和2017年生長季、凍結(jié)期、凍融交替期、消融期的蒸散日變化Fig.4 Diurnal variations of evapotranspiration (ET) in growing season, frozen, freeze-thaw and thawed periods for 2016 and 2017

    2.4 平衡蒸散量

    平衡蒸散量(ETeq),是指在沒有水分限制的條件下,生態(tài)系統(tǒng)理論上ET的最大值[25]。為了進一步評價生態(tài)系統(tǒng)蒸散對水分收支的影響,我們對ETeq,ET和降水量(P)的年累計值年變化進行了比較(圖5)。由圖可知,2016與2017年3個變量的變化趨勢相同,均在生長季出現(xiàn)快速增加趨勢,而在其他時期相對緩慢,且ETeq的累計值高于ET和P。2016年ETeq的年累計值為626.9 mm,其中1—11月累計量為612.3 mm,高于2017年同期的573.7 mm。在植物生長季之前,2016年的ET累計值高于P,當(dāng)進入生長季后,P的累計值高于ET,而2017年的ET累計值始終高于P。

    圖5 2016與2017年蒸散量(ET)、降水量(P)和平衡蒸散量(ETeq)的累計值Fig.5 Cumulative evapotranspiration (ET), precipitation (P) and equilibrium evapotranspiration (ETeq) for 2016 and 2017

    3 討論

    3.1 植被對蒸散的影響

    蒸散包括土壤蒸發(fā)和植被蒸騰,生長季中的植被生長狀況直接影響植被蒸騰和土壤蒸發(fā)的變化[4]。由植物生長旺季(7月)的地上生物量(above-ground biomass,AGB)和葉面積指數(shù)(LAI)觀測數(shù)據(jù)可知,2016年的AGB為(192.7±33.1) g/m2,明顯高于2017年同期的(114.6±18.5) g/m2;且2016年LAI為(1.2±0.1) m2/m2,同樣高于2017年同期的(0.8±0.1) m2/m2。然而,2017年生長季中的ET卻明顯高于2016年,說明地上植被量的降低反而增加了該退化高寒草甸生態(tài)系統(tǒng)的蒸散,分析其可能的原因之一:Hu等[27]曾在研究中發(fā)現(xiàn),在植被LAI較低情況下,ET主要由土壤蒸發(fā)控制。因此在該退化草甸中,地上生物量減少使植被蓋度明顯降低,土壤裸露面積增加,提高了土壤蒸發(fā)量,從而使生態(tài)系統(tǒng)ET增加,這在一些研究中也得到證實[28-29]。Zhang等[28]利用模型對三江源區(qū)退化草甸和人工草地能量分配的研究結(jié)果發(fā)現(xiàn),退化草甸消耗的潛熱高于人工草地,主要原因是退化草甸的土壤蒸發(fā)遠(yuǎn)高于人工草地。

    ET/P是描述水分收支的重要參數(shù)之一[7]。通過與不同植被類型生態(tài)系統(tǒng)的蒸散結(jié)果比較可知(表4),本研究退化高寒草甸生態(tài)系統(tǒng)兩年ET/P分別為0.95和1.09,該值在已報道的草地生態(tài)系統(tǒng)ET/P的范圍(40%—150%)之間[34],但明顯高于Gu等報道的未退化海北高寒草甸的0.60[6]和Li等報道的日本溫帶草地的0.66[31],但低于內(nèi)蒙古退化草地的1.4[30]。本結(jié)果意味著植被退化可能會加劇該高寒草甸的水分散失,降低生態(tài)系統(tǒng)的水源涵養(yǎng)能力。

    3.2 環(huán)境因子對蒸散的影響3.2.1 輻射與溫度對蒸散的影響

    太陽輻射能(Rs)是驅(qū)動陸地生態(tài)系統(tǒng)蒸散的主要動力[6],由于青藏高原海拔高,因此接收的太陽總輻射遠(yuǎn)高于低海拔地區(qū)。凈輻射(Rn)是供給生態(tài)系統(tǒng)蒸散的有效能量,在不受水分限制條件下,ET主要受Rn控制[28]。2016和2017年ET隨Rn的變化趨勢相似(圖6),當(dāng)Rn<4 MJ m-2d-1時,由于土壤處于凍結(jié)狀態(tài),ET維持在相對較低水平,且基本不受Rn變化的影響;當(dāng)Rn≥ 4 MJ m-2d-1時,隨著土壤解凍的開始,ET隨Rn的增加而呈現(xiàn)直線上升的趨勢,該結(jié)果與張立鋒等[28]報道的相一致。但在相同Rn條件下,2017年ET值高于2016年,尤其在生長季Rn> 10 MJ m-2d-1后,兩年差異更明顯,且2017年擬合直線的斜率高于2016年,說明2017年ET隨Rn上升而增加的速率更快,即ET對Rn變化的響應(yīng)更敏感。

    表4 不同草地類型蒸散量和降水量的比較

    研究期間,該生態(tài)系統(tǒng)接收的太陽輻射年總量高達(dá)6594.32 MJ/m2,但由于高原上相對較低的大氣逆輻射和相對較高的反照率,反而導(dǎo)致接收的Rn較低[35]。本研究兩年Rn/Rs平均值0.43,接近于已報道的海北高寒草甸的0.44[6],明顯低于全球平均值0.61和一些已經(jīng)報道的草地系統(tǒng)[36]。因此該研究地驅(qū)動蒸散的有效能Rn反而比低海拔地區(qū)低[36],從而導(dǎo)致生態(tài)系統(tǒng)的蒸散對Rn的響應(yīng)異常敏感。另外,由于2017年高寒草甸的AGB和LAI低于2016年,使土壤蒸發(fā)在蒸散中所占的比例增加,Zhang等[29]指出三江源區(qū)高寒草甸的退化加劇了生態(tài)系統(tǒng)的蒸散量,這可能是導(dǎo)致2017年ET對Rn變化的響應(yīng)比2016年更為敏感的原因之一,此外相關(guān)環(huán)境因子的影響尚有待于進一步研究。

    圖6 2016和2017年退化草甸蒸散量(ET)與凈輻射(Rn)和空氣溫度(Ta)的關(guān)系Fig.6 Relationships between evapotranspiration (ET) and net radiation (Rn), air temperature (Ta) in degraded meadow for 2016 and 2017 ET與 Rn: 2016年y=0.213x-0.49 (R2=0.99) (P<0.01); 2017年y=0.255x-0.53 (R2=0.98)(P<0.01); ET與Ta: 2016年 y=0.84e0.098x (R2=0.94) (P<0.01); 2017年y=1.09e0.094x (R2=0.98) (P<0.01)

    溫度不僅控制植物生長,也是影響生態(tài)系統(tǒng)蒸散的重要因子,氣溫的升高可增加邊界層的水熱交換,從而提高蒸發(fā)的速率[37]。本研究ET均隨Ta升高而成指數(shù)上升趨勢(圖6),該結(jié)果與陳小平等[37]報道的草地生態(tài)系統(tǒng)相似。需要指出的是,許多平原地區(qū)的長期研究表明空氣溫度與蒸散相關(guān)度不高[38],而三江源區(qū)高寒草甸生態(tài)系統(tǒng)與之相反,主要是由于其長期處于低溫環(huán)境,因此ET對溫度的響應(yīng)更為敏感,ET隨溫度的變化速率更快。另外,由于2017年平均氣溫為-0.41℃低于2016年的-0.35℃,尤其在生長季2017年平均氣溫為6.8℃,明顯低于2016年的7.2℃(圖2),這可能是導(dǎo)致2017年ET隨溫度升高而增加的速率略高于2016年的原因之一。

    3.2.2飽和水氣壓差與土壤水分對蒸散的影響

    飽和水汽壓差(VPD)是反映空氣對水分的需求能力[39],VPD通過影響冠層導(dǎo)度gc,從而影響生態(tài)系統(tǒng)的蒸散[31]。本研究ET均隨VPD的增加而呈直線上升趨勢(圖7),這與已報道的一些研究結(jié)果類似[39-40]。然而,2017年擬合直線的斜率為2.64明顯高于2016年的1.46,說明2017年ET隨VPD的增加而上升的速率更快,進而導(dǎo)致在相同VPD條件下,2017年的蒸散量高于2016。低的VPD通??山档痛髿庹羯⑿枨?而高VPD表征空氣干燥,可促進土壤蒸發(fā)和植被蒸騰[5]。本研究VPD變化范圍在0.05-1.69 kPa之間,明顯低于已報道的其他草地生態(tài)系統(tǒng)2-5 kPa的變動范圍[4,41],這可能是導(dǎo)致本研究ET對VPD變化響應(yīng)敏感的原因之一。

    圖7 2016和2017年退化草甸蒸散量(ET)與飽和水汽壓差(VPD)和 5 cm土壤含水量(SWC5)的關(guān)系Fig.7 Relationships between evapotranspiration (ET) and vapor pressure deficit (VPD), soil water content at 5 cm depth (SWC5) in degraded meadow for 2016 and 2017 ET與VPD: 2016年: y=1.46x+0.36 (R2=0.82) (P<0.01); 2017年y=2.64x+0.09 (R2=0.92) (P<0.01); ET與SWC5: 2016年 y=31.73x2+18.06x-0.75 (R2=0.83) (P<0.01); 2017年y=-55.16x2+26.36x-1.03 (R2=0.84) (P<0.01)

    土壤水分是土壤蒸發(fā)和植物蒸騰的直接來源,是影響陸地生態(tài)系統(tǒng)ET的關(guān)鍵因子[28]。高寒草甸的根系主要分布在0—10 cm深度的表層土壤[5],因此生態(tài)系統(tǒng)的ET與5 cm的土壤含水量(SWC5)密切相關(guān)。當(dāng)SWC5低于0.25 m3/m3,兩年ET隨SWC5的升高而增加,當(dāng)高于0.25 m3/m3時,ET隨SWC5的升高呈平穩(wěn)(2016年)或下降(2017年)趨勢(圖7),這與已報道的一些草地研究結(jié)果類似[28,42]。通常,當(dāng)土壤水分受限時,ET隨土壤濕度的升高而增加,而當(dāng)土壤水分不受限時,ET隨土壤水分升高的響應(yīng)不敏感或呈下降趨勢[28]。由于2017年SWC5最高值出現(xiàn)在9月(圖2),而該時期驅(qū)動蒸散的凈輻射和溫度遠(yuǎn)低于7-8月的值(圖2),可能是導(dǎo)致2017年的蒸散隨SWC5的上升而降低的原因。

    3.2.3環(huán)境因子綜合作用對蒸散的影響

    為深入探討環(huán)境因子對ET的影響,對ET與主要環(huán)境因子(Rn、Ta、SWC5、VPD)進行逐步回歸分析。由相關(guān)系數(shù)(表5)可知,該生態(tài)系統(tǒng)ET與上述環(huán)境因子均呈極顯著相關(guān),其中Rn和Ta相關(guān)性最高。2016和2017年回歸方程分別為:ET=0.639Rn+0.291Ta+0.092SWC5-0.046 (R2=0.866,P=0.001);ET=0.532Rn+0.439Ta+0.353 (R2=0.787,P=0.000),表明2016年Rn對ET的影響最大,其次是Ta和SWC5;而2017年ET主要受Rn和Ta的影響。由此可見,該生態(tài)系統(tǒng)的ET變化主要受Rn和Ta的控制。

    3.3 冠層導(dǎo)度gc和解耦系數(shù)Ω對蒸散的影響

    冠層導(dǎo)度(gc)是影響蒸散的重要因子,其變化受LAI、SWC、VPD等環(huán)境因子的影響,通常ET隨gc的升高而增加[43]。為評價gc對ET的影響,本研究計算了4個時段晴天條件下[30](晴空指數(shù)大于0.7)的gc日變化(圖8),可以得出生長季中的gc最高,其次為消融期,凍融交替期,而凍結(jié)期的gc最低,該結(jié)果與已報道的加拿大的溫帶草地、海北高寒草甸等生態(tài)系統(tǒng)變化趨勢一致[4,6],且與不同時段的蒸散速率變化相一致(圖4),主要是由于高的ET一般與高的gc相關(guān),生長季植被生長茂盛以及較高的土壤濕度,進而導(dǎo)致gc升高,蒸散速率加快[44]。然而,2017年生長季gc的平均值為10.3 mm/s,明顯高于2016年的8.2 mm/s,這是導(dǎo)致2017年生長季ET高于2016年的原因之一。兩年凍融交替期和凍結(jié)期gc的平均值小于3 mm/s,且2016年與2017年相差不大。

    表5 2016和2017年退化草甸ET與主要環(huán)境因子的Pearson相關(guān)系數(shù)

    兩年的gc日變化趨勢相同,生長季和消融期gc最高值出現(xiàn)在正午前后,而另外兩個時段gc最高值出現(xiàn)在8:00—11:00之間,該結(jié)果與Wever等[4]和Zhang等[29]研究結(jié)果一致。在凍結(jié)期和凍融交替期,由于植物凋萎,地面有效輻射加強,溫度和VPD較低,這可能是導(dǎo)致gc在上午出現(xiàn)高值的原因[28];而在降水豐沛的生長季和消融期,由于生長季水分不受限制,溫度和VPD相對較高,日間gc隨VPD的增大而增強,但當(dāng)VPD過大時,導(dǎo)致氣孔關(guān)閉,gc降低,因此gc的最大值出現(xiàn)在正午前后[4,44],這也和ET的最大值出現(xiàn)的時間相吻合。

    圖8 2016年和2017年晴天條件下退化草甸凍結(jié)期、凍融交替期、消融期、生長季的冠層導(dǎo)度(gc)日變化Fig.8 Diurnal pattern of canopy conductance (gc) on clear days in growing season, frozen, freeze-thaw and thawed periods in degraded meadow for 2016 and 2017

    解耦系數(shù)Ω是評價地表與大氣之間水熱交換耦合狀況的參數(shù)[45],其值在0—1之間變化,當(dāng)Ω越接近0時,表明ET主要受VPD控制;Ω越接近1,表明ET主要受Rn主導(dǎo)[4,40]。兩年Ω的變化趨勢相同,最高值出現(xiàn)在生長季和消融期(圖9),2016和2017年生長季Ω的平均值分別為0.55和0.60,而凍結(jié)期和凍融交替期的Ω明顯低于生長季和消融期,這與已報道的諸多研究結(jié)果相同[6,17,28]。該結(jié)果說明生長季和消融期的ET主要受Rn影響,而凍結(jié)期和凍融交替期的Rn并不是驅(qū)動ET的主要因子。另外,2017年生長季和消融期的Ω均高于2016年,可見Rn對2017年ET的影響大于2016年。

    年Ω日變化均呈現(xiàn)上午高于下午的變化趨勢,該變化與Wever等[4]和Hao等[40]報道的草地結(jié)果相同,說明該生態(tài)系統(tǒng)Rn對午前ET的影響相對較大,而VPD對午后ET的影響相對較大。

    圖9 2016年和2017年晴天條件下退化草甸凍結(jié)期、凍融交替期、消融期、生長季的解耦系數(shù)(Ω)日變化Fig.9 Diurnal pattern of decoupling coefficient (Ω) on clear days in growing season, frozen, freeze-thaw and thawed periods in degraded meadow for 2016 and 2017

    圖10 2016年和2017年退化草甸晴天條件下日均Priestley-Taylor系數(shù)(ET/ETeq)與冠層導(dǎo)度(gc)之間的關(guān)系Fig.10 Relationship between the daily average Priestley-Taylor coefficient (α=ET / ETeq) and canopy conductance (gc) on clear days in degraded meadow for 2016 and 20172016年:y=0.2696ln(x)+0.2007 (R2=0.87) (P<0.01); 2017年:y=0.294ln(x)+0.2554 (R2=0.87) (P<0.01)

    Priestley-Taylor系數(shù)(α=ET/ETeq)為蒸散與平衡蒸散之比,其能夠反映土壤中可利用水分對蒸散的影響,通常水分供應(yīng)充足的生態(tài)系統(tǒng)ET/ETeq高于水分受限的生態(tài)系統(tǒng)[46]。晴天條件下,兩年ET/ETeq與gc均呈現(xiàn)對數(shù)增長趨勢(圖10),這與之前報道的很多結(jié)果相一致[5,29,32], 這是由于生長季的gc遠(yuǎn)高于凍結(jié)期和凍融交替期(圖8),而相對較高的gc導(dǎo)致生態(tài)系統(tǒng)更多的ET[6]。另外,2017年ET/ETeq值高于2016年,其中2016和2017年生長季的ET/ETeq分別為0.78和0.93,明顯高于受水分限制的半干旱草地生態(tài)系統(tǒng)[30,32],但與水分充足的森林生態(tài)系統(tǒng)相當(dāng)[47]。結(jié)果說明該生態(tài)系統(tǒng)ET受水分限制相對較小,且2017年生態(tài)系統(tǒng)受到的水分限制相對小于2016年,進而導(dǎo)致2016年ET低于2017年。

    3.4 凍融變化對蒸散的影響

    2016與2017年土壤凍融交替期、凍結(jié)期及消融期的長短不同(表2),凍融交替期分別為21 d和13 d,凍結(jié)期分別為140 d和132 d,而消融期分別為205 d和220 d。該結(jié)果說明2016年的凍結(jié)期和凍融交替期比2017年的時間長,而消融期則比2017年短15 d,進而導(dǎo)致2017年消融期(包括生長季)的ET比2016年高63.1 mm(表4);而2017年凍融交替期比2016年縮短8 d,其ET比 2016年低19.9 mm。

    由以上結(jié)果可知,凍融時間的變化可引起地氣之間能量和水分交換狀況的改變,從而導(dǎo)致生態(tài)系統(tǒng)蒸散的變化[26]。通常土壤凍結(jié)可降低土壤蒸發(fā),減少土壤水分的散失,當(dāng)土壤消融后,相對較高的土壤含水量及相對較低的地表反照率(使土壤吸收更多的太陽輻射),以及輻射和溫度的上升,進而加劇了土壤蒸發(fā)量[48]。有研究表明,溫度(特別是冬天)的升高,可能導(dǎo)致青藏高原土壤消融期提前、凍結(jié)推遲以及高原凍結(jié)天數(shù)持續(xù)下降[49]。因此在全球氣候變化背景下,三江源區(qū)土壤凍融時間將會發(fā)生改變,從而對蒸散變化產(chǎn)生深刻影響,而未來氣候變暖可能會加劇該退化高寒草甸生態(tài)系統(tǒng)系統(tǒng)的蒸散量,降低其水源涵養(yǎng)能力。

    4 結(jié)論

    三江源區(qū)退化高寒草甸生態(tài)系統(tǒng)的年蒸散量可達(dá)降水量的95%以上,且存在明顯的年際間差異,植被生物量相對較高的年份,其生態(tài)系統(tǒng)的蒸散量反而低于生物量相對較低的年份,由此可推測該區(qū)域高寒草甸的退化可能將導(dǎo)致其生態(tài)系統(tǒng)水源涵養(yǎng)功能的降低。同時,土壤凍融時間的改變對其生態(tài)系統(tǒng)的蒸散產(chǎn)生了重要影響,土壤凍結(jié)期的縮短將引起蒸散量的增加。由于影響該生態(tài)系統(tǒng)蒸散的最主要環(huán)境因子是凈輻射和溫度,因此在全球氣候變暖的背景下,三江源區(qū)退化高寒草甸生態(tài)系統(tǒng)的蒸散可能出現(xiàn)增加趨勢,進而對該區(qū)域的水分收支產(chǎn)生深刻的影響。

    猜你喜歡
    源區(qū)草甸凍融
    高山草甸
    The world's narrowest river
    山地草甸
    冬小麥蒸散源區(qū)代表性分析
    武功山山地草甸的成因調(diào)查
    活力(2019年21期)2019-04-01 12:17:10
    間苯三酚在凍融胚胎移植中的應(yīng)用
    渭河源區(qū)徑流量變化特征及趨勢分析
    反復(fù)凍融作用下巖橋破壞的試驗研究
    降調(diào)節(jié)方案在凍融胚胎移植周期中的應(yīng)用
    凍融后小鼠卵巢移植到雄鼠腎被膜下卵泡的生長發(fā)育
    成熟少妇高潮喷水视频| 中文字幕高清在线视频| 国产人伦9x9x在线观看| 欧美日韩亚洲国产一区二区在线观看| e午夜精品久久久久久久| 欧美人与性动交α欧美精品济南到| bbb黄色大片| 一卡2卡三卡四卡精品乱码亚洲| svipshipincom国产片| 久久伊人香网站| 久久人人爽av亚洲精品天堂| 色综合欧美亚洲国产小说| 久99久视频精品免费| 桃红色精品国产亚洲av| 日韩大码丰满熟妇| 欧美成人一区二区免费高清观看 | 久久久久久久久久久久大奶| 69av精品久久久久久| 久久中文字幕人妻熟女| 美女高潮到喷水免费观看| 亚洲精品美女久久av网站| 久久香蕉国产精品| 女性生殖器流出的白浆| 19禁男女啪啪无遮挡网站| 女生性感内裤真人,穿戴方法视频| 久久久国产精品麻豆| 一二三四社区在线视频社区8| 亚洲精品美女久久久久99蜜臀| 99在线人妻在线中文字幕| 日韩三级视频一区二区三区| 女性生殖器流出的白浆| 成在线人永久免费视频| 黄色视频,在线免费观看| 最近最新中文字幕大全免费视频| 中文字幕另类日韩欧美亚洲嫩草| 免费看美女性在线毛片视频| 免费人成视频x8x8入口观看| 久久天躁狠狠躁夜夜2o2o| 色精品久久人妻99蜜桃| 欧美在线一区亚洲| a在线观看视频网站| 后天国语完整版免费观看| 中文字幕人妻熟女乱码| 操出白浆在线播放| 很黄的视频免费| 女人被狂操c到高潮| 在线av久久热| 国产97色在线日韩免费| 婷婷丁香在线五月| 91麻豆av在线| 丰满人妻熟妇乱又伦精品不卡| 最新美女视频免费是黄的| 9色porny在线观看| 国内精品久久久久精免费| 亚洲在线自拍视频| av天堂久久9| 黑人欧美特级aaaaaa片| www.熟女人妻精品国产| 国产av一区在线观看免费| 午夜久久久在线观看| 色av中文字幕| 色综合欧美亚洲国产小说| 成人三级做爰电影| 老司机福利观看| 中出人妻视频一区二区| 女人被狂操c到高潮| 成人三级做爰电影| 在线av久久热| www.精华液| 一卡2卡三卡四卡精品乱码亚洲| 欧美另类亚洲清纯唯美| 午夜久久久在线观看| av天堂久久9| 免费观看精品视频网站| 中文字幕久久专区| 久久久久久久久免费视频了| 国产aⅴ精品一区二区三区波| 欧美大码av| 久久久久久久久免费视频了| 啦啦啦观看免费观看视频高清 | 日韩中文字幕欧美一区二区| 欧美日韩精品网址| 韩国av一区二区三区四区| 亚洲自偷自拍图片 自拍| 女人爽到高潮嗷嗷叫在线视频| 久久精品国产综合久久久| 国产日韩一区二区三区精品不卡| 久久香蕉激情| 1024视频免费在线观看| 一个人观看的视频www高清免费观看 | 亚洲五月天丁香| 国产色视频综合| 在线观看免费午夜福利视频| 欧美丝袜亚洲另类 | 黄色a级毛片大全视频| 我的亚洲天堂| 91九色精品人成在线观看| 亚洲午夜精品一区,二区,三区| 国产精品免费视频内射| 日本免费一区二区三区高清不卡 | 精品熟女少妇八av免费久了| 欧美黄色淫秽网站| 亚洲免费av在线视频| 婷婷丁香在线五月| 可以在线观看毛片的网站| 欧美绝顶高潮抽搐喷水| 成人国产综合亚洲| 黄色成人免费大全| 欧美一区二区精品小视频在线| 亚洲 国产 在线| 美女大奶头视频| 国产91精品成人一区二区三区| 久久久久久久午夜电影| 国产成人精品久久二区二区91| 午夜精品在线福利| 女人被狂操c到高潮| 午夜久久久久精精品| 窝窝影院91人妻| 熟妇人妻久久中文字幕3abv| 激情视频va一区二区三区| 国产男靠女视频免费网站| 久久天堂一区二区三区四区| 亚洲中文av在线| 大型av网站在线播放| 一区在线观看完整版| 国产精品,欧美在线| 亚洲精品粉嫩美女一区| 黄频高清免费视频| 又大又爽又粗| 久久这里只有精品19| 色av中文字幕| 国产av一区在线观看免费| 看免费av毛片| 精品国产国语对白av| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲天堂国产精品一区在线| 亚洲第一电影网av| 真人做人爱边吃奶动态| 少妇熟女aⅴ在线视频| 久久精品91蜜桃| 韩国精品一区二区三区| 婷婷六月久久综合丁香| 在线天堂中文资源库| 一夜夜www| 麻豆国产av国片精品| 老司机福利观看| 看免费av毛片| 日韩精品中文字幕看吧| 精品一区二区三区视频在线观看免费| 99国产极品粉嫩在线观看| 男女做爰动态图高潮gif福利片 | www.www免费av| 性少妇av在线| 伦理电影免费视频| 在线观看舔阴道视频| 国产精品亚洲一级av第二区| 一进一出抽搐gif免费好疼| 日韩成人在线观看一区二区三区| 午夜福利,免费看| 亚洲欧美一区二区三区黑人| 久久久久九九精品影院| 女人被狂操c到高潮| 人妻丰满熟妇av一区二区三区| 韩国av一区二区三区四区| 国产一区二区三区视频了| 国产又爽黄色视频| 女性生殖器流出的白浆| 麻豆av在线久日| 99香蕉大伊视频| svipshipincom国产片| 黑人操中国人逼视频| 亚洲专区国产一区二区| 91字幕亚洲| 在线永久观看黄色视频| 91av网站免费观看| 国产欧美日韩综合在线一区二区| 精品久久久久久成人av| 国产精品野战在线观看| 精品不卡国产一区二区三区| 欧美一区二区精品小视频在线| 日韩精品免费视频一区二区三区| 久久久久久久午夜电影| 视频区欧美日本亚洲| 黄频高清免费视频| av网站免费在线观看视频| 99riav亚洲国产免费| av片东京热男人的天堂| 怎么达到女性高潮| 久热爱精品视频在线9| 99国产极品粉嫩在线观看| 久久久水蜜桃国产精品网| 高清黄色对白视频在线免费看| 久99久视频精品免费| 国产片内射在线| 亚洲欧美日韩另类电影网站| 国产高清视频在线播放一区| 久久精品国产亚洲av香蕉五月| 国产精品av久久久久免费| 麻豆国产av国片精品| 日韩欧美在线二视频| 免费在线观看视频国产中文字幕亚洲| 麻豆久久精品国产亚洲av| 亚洲片人在线观看| 黄频高清免费视频| 精品一区二区三区四区五区乱码| 久久影院123| 黄片播放在线免费| 欧美绝顶高潮抽搐喷水| 韩国精品一区二区三区| 国产精品二区激情视频| 麻豆av在线久日| 一卡2卡三卡四卡精品乱码亚洲| 欧美日本视频| 不卡一级毛片| 老司机靠b影院| 男人舔女人下体高潮全视频| 久久中文字幕人妻熟女| 午夜精品久久久久久毛片777| xxx96com| 午夜两性在线视频| 亚洲精品一区av在线观看| 久久久精品欧美日韩精品| 日韩有码中文字幕| 久久国产精品人妻蜜桃| 国产亚洲精品久久久久久毛片| 搡老岳熟女国产| 少妇裸体淫交视频免费看高清 | 少妇熟女aⅴ在线视频| 99国产极品粉嫩在线观看| 久久香蕉国产精品| 欧美人与性动交α欧美精品济南到| 天堂影院成人在线观看| 一级a爱视频在线免费观看| 美女大奶头视频| 日韩欧美免费精品| 精品一区二区三区四区五区乱码| 精品高清国产在线一区| 99久久久亚洲精品蜜臀av| 中文字幕人妻丝袜一区二区| 久久久久国内视频| 国产单亲对白刺激| 无限看片的www在线观看| 亚洲精品一区av在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 两个人免费观看高清视频| 色播亚洲综合网| 校园春色视频在线观看| 久久久久国产精品人妻aⅴ院| 亚洲第一av免费看| 18禁美女被吸乳视频| 可以免费在线观看a视频的电影网站| 搡老熟女国产l中国老女人| 日韩欧美国产一区二区入口| 国产一区二区三区在线臀色熟女| 18禁黄网站禁片午夜丰满| 91精品国产国语对白视频| 一区二区三区高清视频在线| 欧美黄色片欧美黄色片| www.www免费av| 人人澡人人妻人| 亚洲国产高清在线一区二区三 | 成人av一区二区三区在线看| 亚洲精品美女久久av网站| 香蕉国产在线看| 精品国内亚洲2022精品成人| 国产欧美日韩一区二区三区在线| 一二三四社区在线视频社区8| 午夜免费激情av| 亚洲欧洲精品一区二区精品久久久| 久久国产精品影院| 日韩欧美三级三区| 欧美成狂野欧美在线观看| 制服人妻中文乱码| 国产亚洲精品久久久久久毛片| 久久中文看片网| 久99久视频精品免费| 男男h啪啪无遮挡| 此物有八面人人有两片| 国产精品二区激情视频| 久久午夜亚洲精品久久| 我的亚洲天堂| 免费在线观看影片大全网站| 欧美激情极品国产一区二区三区| 日本撒尿小便嘘嘘汇集6| 99久久99久久久精品蜜桃| 国产精品久久久av美女十八| 国产成人精品无人区| 成人18禁在线播放| 一边摸一边做爽爽视频免费| 欧美老熟妇乱子伦牲交| a在线观看视频网站| 啦啦啦观看免费观看视频高清 | 久99久视频精品免费| 精品卡一卡二卡四卡免费| 人人妻人人澡欧美一区二区 | 精品久久久久久久久久免费视频| 亚洲一区二区三区色噜噜| 亚洲精品国产精品久久久不卡| 制服诱惑二区| 色播在线永久视频| 亚洲全国av大片| 99在线人妻在线中文字幕| 午夜老司机福利片| 亚洲专区字幕在线| avwww免费| 一边摸一边抽搐一进一小说| 两个人看的免费小视频| 少妇粗大呻吟视频| av有码第一页| 非洲黑人性xxxx精品又粗又长| 又黄又粗又硬又大视频| 国产高清激情床上av| 成人手机av| 精品卡一卡二卡四卡免费| 在线视频色国产色| 亚洲精品久久国产高清桃花| 一进一出抽搐gif免费好疼| 亚洲熟妇中文字幕五十中出| 可以在线观看的亚洲视频| 女人爽到高潮嗷嗷叫在线视频| 国产成人系列免费观看| 久久精品成人免费网站| 亚洲在线自拍视频| 亚洲少妇的诱惑av| 亚洲一码二码三码区别大吗| 九色国产91popny在线| 日韩免费av在线播放| 制服丝袜大香蕉在线| 国产精品 欧美亚洲| 欧美一区二区精品小视频在线| 国产精品二区激情视频| av福利片在线| 精品卡一卡二卡四卡免费| 国产高清激情床上av| 一区在线观看完整版| av电影中文网址| av欧美777| 如日韩欧美国产精品一区二区三区| 丁香六月欧美| 亚洲精品在线观看二区| 99国产精品一区二区蜜桃av| 国产成人av激情在线播放| 欧美+亚洲+日韩+国产| 国内毛片毛片毛片毛片毛片| 校园春色视频在线观看| 妹子高潮喷水视频| 久久久久久久午夜电影| 亚洲avbb在线观看| 欧美亚洲日本最大视频资源| 国内毛片毛片毛片毛片毛片| 法律面前人人平等表现在哪些方面| 精品一区二区三区视频在线观看免费| 国产一区二区三区在线臀色熟女| 国产精品 欧美亚洲| 亚洲国产精品合色在线| 99国产极品粉嫩在线观看| 亚洲全国av大片| 国产精品亚洲一级av第二区| 亚洲精品国产一区二区精华液| 免费观看人在逋| av视频免费观看在线观看| 午夜福利一区二区在线看| 成人亚洲精品一区在线观看| 中出人妻视频一区二区| 中文亚洲av片在线观看爽| videosex国产| 亚洲av电影不卡..在线观看| 身体一侧抽搐| 欧美日韩乱码在线| 亚洲欧洲精品一区二区精品久久久| www.自偷自拍.com| 高清黄色对白视频在线免费看| 日韩大尺度精品在线看网址 | 国产欧美日韩一区二区精品| 国产伦一二天堂av在线观看| 一区二区日韩欧美中文字幕| 久久香蕉激情| 在线十欧美十亚洲十日本专区| 国产精品,欧美在线| av片东京热男人的天堂| 久热爱精品视频在线9| 亚洲av成人av| 国产97色在线日韩免费| 久久国产精品影院| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品自产拍在线观看55亚洲| 亚洲国产中文字幕在线视频| 成人特级黄色片久久久久久久| 免费看美女性在线毛片视频| 日韩欧美国产在线观看| 乱人伦中国视频| 这个男人来自地球电影免费观看| 搡老熟女国产l中国老女人| 1024香蕉在线观看| av片东京热男人的天堂| 精品国产亚洲在线| 久久中文字幕一级| 淫妇啪啪啪对白视频| 波多野结衣av一区二区av| 1024视频免费在线观看| 亚洲精品av麻豆狂野| 一区二区日韩欧美中文字幕| 久久香蕉激情| 一二三四在线观看免费中文在| 色综合站精品国产| 亚洲男人的天堂狠狠| 亚洲 国产 在线| 一级,二级,三级黄色视频| 亚洲国产欧美网| 欧美老熟妇乱子伦牲交| 女人精品久久久久毛片| 国产精品亚洲一级av第二区| 嫩草影院精品99| 成人18禁高潮啪啪吃奶动态图| 日韩成人在线观看一区二区三区| 制服诱惑二区| 欧美 亚洲 国产 日韩一| 中文字幕人妻丝袜一区二区| 中出人妻视频一区二区| 12—13女人毛片做爰片一| 久久久精品欧美日韩精品| 正在播放国产对白刺激| 午夜免费成人在线视频| 久久天躁狠狠躁夜夜2o2o| 久久久久久久久中文| 亚洲av美国av| 色综合欧美亚洲国产小说| 日本免费一区二区三区高清不卡 | av电影中文网址| 国产私拍福利视频在线观看| 女人被狂操c到高潮| 成人国产综合亚洲| 国产黄a三级三级三级人| 国产精品久久电影中文字幕| 亚洲精品中文字幕在线视频| 亚洲狠狠婷婷综合久久图片| 两个人看的免费小视频| 黑人巨大精品欧美一区二区蜜桃| 亚洲精华国产精华精| 国产成人免费无遮挡视频| 久久亚洲精品不卡| 国产成人av教育| 夜夜看夜夜爽夜夜摸| 搡老妇女老女人老熟妇| 亚洲欧美日韩高清在线视频| 国产主播在线观看一区二区| 热re99久久国产66热| 色播在线永久视频| 91麻豆精品激情在线观看国产| 日本撒尿小便嘘嘘汇集6| 妹子高潮喷水视频| 午夜福利免费观看在线| 久久热在线av| 午夜a级毛片| 日韩一卡2卡3卡4卡2021年| 999久久久精品免费观看国产| 69av精品久久久久久| 男男h啪啪无遮挡| 岛国视频午夜一区免费看| 国产成人欧美在线观看| 亚洲国产欧美日韩在线播放| 两个人视频免费观看高清| 国产1区2区3区精品| 精品一品国产午夜福利视频| 亚洲人成伊人成综合网2020| 国产视频一区二区在线看| 日韩欧美免费精品| 男男h啪啪无遮挡| 亚洲久久久国产精品| 精品日产1卡2卡| 精品久久久久久,| 久久香蕉激情| 一边摸一边抽搐一进一出视频| 十八禁网站免费在线| 国内久久婷婷六月综合欲色啪| 午夜视频精品福利| 18美女黄网站色大片免费观看| 国产亚洲精品综合一区在线观看 | 国产精品电影一区二区三区| 久久天堂一区二区三区四区| 十分钟在线观看高清视频www| 999久久久国产精品视频| 男人的好看免费观看在线视频 | 欧美成人性av电影在线观看| 国产精品野战在线观看| 男女下面进入的视频免费午夜 | 叶爱在线成人免费视频播放| 日韩大尺度精品在线看网址 | 首页视频小说图片口味搜索| 色哟哟哟哟哟哟| 免费看十八禁软件| 搞女人的毛片| 久久久精品欧美日韩精品| 亚洲欧洲精品一区二区精品久久久| 国产又爽黄色视频| 久久精品国产99精品国产亚洲性色 | 色播在线永久视频| 黄色视频不卡| 99精品在免费线老司机午夜| 亚洲,欧美精品.| 美女扒开内裤让男人捅视频| 12—13女人毛片做爰片一| 亚洲色图 男人天堂 中文字幕| 色老头精品视频在线观看| 久久中文字幕人妻熟女| 亚洲国产欧美网| 欧美日本中文国产一区发布| 91老司机精品| 在线免费观看的www视频| 级片在线观看| 久热这里只有精品99| 成熟少妇高潮喷水视频| 亚洲国产精品sss在线观看| 精品国产一区二区久久| 亚洲 国产 在线| 久久九九热精品免费| 亚洲三区欧美一区| 亚洲av熟女| 日韩精品中文字幕看吧| 亚洲狠狠婷婷综合久久图片| 免费观看精品视频网站| 黄片小视频在线播放| 欧美中文综合在线视频| 啦啦啦免费观看视频1| 成人国语在线视频| 久久久国产欧美日韩av| 精品熟女少妇八av免费久了| 麻豆成人av在线观看| 老鸭窝网址在线观看| 一本大道久久a久久精品| 亚洲精品中文字幕一二三四区| 欧美乱码精品一区二区三区| 亚洲狠狠婷婷综合久久图片| 国产激情欧美一区二区| 99在线人妻在线中文字幕| 窝窝影院91人妻| 亚洲无线在线观看| 亚洲中文日韩欧美视频| 国产精品爽爽va在线观看网站 | 十分钟在线观看高清视频www| 91成人精品电影| 久久久国产成人精品二区| 久久国产精品男人的天堂亚洲| 在线av久久热| 涩涩av久久男人的天堂| 国产成人欧美在线观看| 国产又色又爽无遮挡免费看| 国产成人啪精品午夜网站| 色综合亚洲欧美另类图片| 亚洲天堂国产精品一区在线| 国产精品免费一区二区三区在线| 欧美大码av| 国产精品一区二区三区四区久久 | 日本撒尿小便嘘嘘汇集6| 日本五十路高清| 亚洲免费av在线视频| 国产精品久久久av美女十八| 成熟少妇高潮喷水视频| 美女免费视频网站| 法律面前人人平等表现在哪些方面| 日韩欧美一区二区三区在线观看| 欧美大码av| 天堂动漫精品| 看片在线看免费视频| 久久精品人人爽人人爽视色| 欧美乱码精品一区二区三区| 极品教师在线免费播放| 不卡av一区二区三区| 亚洲黑人精品在线| 女同久久另类99精品国产91| 久久午夜亚洲精品久久| 男人舔女人的私密视频| 极品人妻少妇av视频| 一级a爱片免费观看的视频| 亚洲avbb在线观看| 久久精品国产清高在天天线| 淫妇啪啪啪对白视频| 女生性感内裤真人,穿戴方法视频| 亚洲在线自拍视频| 亚洲中文av在线| 亚洲中文日韩欧美视频| 亚洲狠狠婷婷综合久久图片| 超碰成人久久| 亚洲aⅴ乱码一区二区在线播放 | 99香蕉大伊视频| 99国产极品粉嫩在线观看| 国产单亲对白刺激| 国产蜜桃级精品一区二区三区| xxx96com| 国产一区二区激情短视频| 91九色精品人成在线观看| 两个人视频免费观看高清| 精品少妇一区二区三区视频日本电影| 国产精华一区二区三区| 亚洲专区中文字幕在线| 欧美丝袜亚洲另类 | 啪啪无遮挡十八禁网站| 国产乱人伦免费视频| 久久人妻福利社区极品人妻图片| 一级毛片女人18水好多| 欧美激情极品国产一区二区三区| 一级毛片精品| 午夜福利,免费看| 正在播放国产对白刺激| 久久香蕉国产精品| 成人永久免费在线观看视频| 欧美精品啪啪一区二区三区| 日本 欧美在线| 国产99白浆流出| 精品国产一区二区久久| 天天躁狠狠躁夜夜躁狠狠躁| 日本三级黄在线观看| av视频免费观看在线观看| 久久青草综合色| 欧美激情高清一区二区三区|