• <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ā)育
    国产在线一区二区三区精| 久久久久网色| 国产精品欧美亚洲77777| 精品酒店卫生间| 国产日韩欧美视频二区| 综合色丁香网| 性少妇av在线| 交换朋友夫妻互换小说| tube8黄色片| 久久 成人 亚洲| 如日韩欧美国产精品一区二区三区| 高清视频免费观看一区二区| 成年人午夜在线观看视频| 午夜福利视频精品| 2018国产大陆天天弄谢| 国产精品偷伦视频观看了| 亚洲久久久国产精品| 晚上一个人看的免费电影| 日韩欧美精品免费久久| 又大又黄又爽视频免费| 国产av精品麻豆| 伊人久久大香线蕉亚洲五| 人妻人人澡人人爽人人| 精品午夜福利在线看| 日本av免费视频播放| 免费观看av网站的网址| 亚洲精品一二三| 天天躁日日躁夜夜躁夜夜| 国产精品国产三级国产专区5o| 日韩一卡2卡3卡4卡2021年| 欧美日韩精品成人综合77777| 亚洲精品美女久久久久99蜜臀 | 亚洲av综合色区一区| 国产成人av激情在线播放| 欧美日韩亚洲国产一区二区在线观看 | 久久精品国产亚洲av高清一级| 天天躁狠狠躁夜夜躁狠狠躁| 蜜桃国产av成人99| 亚洲精品成人av观看孕妇| 亚洲成人手机| 亚洲精品久久午夜乱码| 亚洲av男天堂| 1024视频免费在线观看| 母亲3免费完整高清在线观看 | 亚洲一码二码三码区别大吗| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久这里只有精品19| 电影成人av| 国产精品 国内视频| 边亲边吃奶的免费视频| 日本-黄色视频高清免费观看| 国产欧美亚洲国产| 国产av一区二区精品久久| 国产野战对白在线观看| 亚洲精品成人av观看孕妇| 女性被躁到高潮视频| 国产成人a∨麻豆精品| 黄色配什么色好看| 国产淫语在线视频| 天天躁夜夜躁狠狠躁躁| 性色avwww在线观看| 亚洲欧美精品自产自拍| 亚洲精华国产精华液的使用体验| 黄色 视频免费看| 国产在视频线精品| 国产精品熟女久久久久浪| 啦啦啦视频在线资源免费观看| 亚洲中文av在线| 91成人精品电影| 亚洲精品日本国产第一区| 卡戴珊不雅视频在线播放| 日韩中文字幕视频在线看片| 亚洲av欧美aⅴ国产| 十八禁网站网址无遮挡| 美女视频免费永久观看网站| 一级毛片我不卡| 青春草亚洲视频在线观看| 日本wwww免费看| 成人毛片60女人毛片免费| 亚洲人成网站在线观看播放| 在线亚洲精品国产二区图片欧美| 日本欧美国产在线视频| 久久久国产一区二区| 亚洲av成人精品一二三区| 精品酒店卫生间| 成年av动漫网址| 久久国产亚洲av麻豆专区| 一区二区av电影网| 看免费成人av毛片| 国产成人精品一,二区| 日本av手机在线免费观看| 最新的欧美精品一区二区| 欧美老熟妇乱子伦牲交| 老熟女久久久| 亚洲精品日本国产第一区| 又粗又硬又长又爽又黄的视频| 波多野结衣一区麻豆| 久久久国产精品麻豆| 在线亚洲精品国产二区图片欧美| 99久国产av精品国产电影| 国产精品亚洲av一区麻豆 | 成人影院久久| 香蕉国产在线看| 叶爱在线成人免费视频播放| 亚洲少妇的诱惑av| 少妇猛男粗大的猛烈进出视频| 日日撸夜夜添| 亚洲成人av在线免费| 考比视频在线观看| 18禁动态无遮挡网站| 2021少妇久久久久久久久久久| videossex国产| 9191精品国产免费久久| 99久久综合免费| 人人澡人人妻人| 国产精品 欧美亚洲| 超色免费av| 肉色欧美久久久久久久蜜桃| 一个人免费看片子| 下体分泌物呈黄色| 久久综合国产亚洲精品| 日韩大片免费观看网站| 国产麻豆69| 婷婷色麻豆天堂久久| 国产熟女欧美一区二区| 在线天堂中文资源库| 人妻 亚洲 视频| 久久这里有精品视频免费| 涩涩av久久男人的天堂| 久久久久网色| 亚洲av在线观看美女高潮| 精品少妇一区二区三区视频日本电影 | 欧美最新免费一区二区三区| 亚洲欧美成人精品一区二区| 国产精品久久久av美女十八| 久久人人97超碰香蕉20202| 曰老女人黄片| 欧美日韩亚洲国产一区二区在线观看 | 久久久精品区二区三区| 国产精品熟女久久久久浪| 国产精品一区二区在线不卡| 99热网站在线观看| 久久99蜜桃精品久久| 亚洲精品自拍成人| 国产精品偷伦视频观看了| 男女啪啪激烈高潮av片| 大陆偷拍与自拍| 午夜免费观看性视频| 久久久国产欧美日韩av| 国产爽快片一区二区三区| 国产精品久久久久久精品古装| 韩国av在线不卡| 王馨瑶露胸无遮挡在线观看| 青春草亚洲视频在线观看| 久久午夜综合久久蜜桃| 国产人伦9x9x在线观看 | 国产成人av激情在线播放| 91久久精品国产一区二区三区| 18禁裸乳无遮挡动漫免费视频| 精品国产一区二区久久| 黄片无遮挡物在线观看| 国产精品偷伦视频观看了| 天天操日日干夜夜撸| 成人毛片60女人毛片免费| 亚洲成人av在线免费| 国产精品一区二区在线观看99| 亚洲国产日韩一区二区| 国产伦理片在线播放av一区| 80岁老熟妇乱子伦牲交| 天天躁夜夜躁狠狠躁躁| 国产精品久久久久久精品古装| 黄网站色视频无遮挡免费观看| 丝袜脚勾引网站| 亚洲欧美日韩另类电影网站| 日韩 亚洲 欧美在线| 99久国产av精品国产电影| 夜夜骑夜夜射夜夜干| 色播在线永久视频| 国精品久久久久久国模美| 黄色怎么调成土黄色| 黄片播放在线免费| 久久影院123| 久久精品久久精品一区二区三区| 观看美女的网站| xxxhd国产人妻xxx| 中国三级夫妇交换| 精品久久久久久电影网| 日韩精品有码人妻一区| 日韩伦理黄色片| 婷婷色综合大香蕉| 精品国产乱码久久久久久小说| 中文字幕人妻熟女乱码| 亚洲国产av影院在线观看| 日产精品乱码卡一卡2卡三| 成年美女黄网站色视频大全免费| 美女脱内裤让男人舔精品视频| 中文字幕人妻丝袜制服| 久久人妻熟女aⅴ| 久久国产精品大桥未久av| 五月开心婷婷网| 亚洲中文av在线| 免费大片黄手机在线观看| 女人久久www免费人成看片| 精品人妻熟女毛片av久久网站| 国产午夜精品一二区理论片| av在线app专区| 伊人久久国产一区二区| 亚洲av综合色区一区| 中文乱码字字幕精品一区二区三区| 99久久综合免费| 99久国产av精品国产电影| 少妇 在线观看| 国产 精品1| 一本久久精品| 国产精品久久久久久av不卡| xxxhd国产人妻xxx| 啦啦啦视频在线资源免费观看| 黄色毛片三级朝国网站| 1024香蕉在线观看| 亚洲av在线观看美女高潮| 国产成人精品一,二区| 啦啦啦啦在线视频资源| 亚洲精品中文字幕在线视频| av国产久精品久网站免费入址| 美国免费a级毛片| 97在线视频观看| 久久热在线av| 黄片播放在线免费| 亚洲av福利一区| 搡老乐熟女国产| 国产成人精品婷婷| 亚洲色图综合在线观看| 亚洲精品av麻豆狂野| 国产成人a∨麻豆精品| 777久久人妻少妇嫩草av网站| 麻豆精品久久久久久蜜桃| 亚洲第一av免费看| 久久av网站| 边亲边吃奶的免费视频| 最新中文字幕久久久久| 久久综合国产亚洲精品| 国产精品一二三区在线看| 精品亚洲乱码少妇综合久久| 一区二区三区乱码不卡18| 午夜老司机福利剧场| 少妇 在线观看| 久久久久精品性色| 亚洲欧美精品自产自拍| 波野结衣二区三区在线| 人人妻人人澡人人爽人人夜夜| 午夜福利在线免费观看网站| 国产白丝娇喘喷水9色精品| 人人妻人人澡人人爽人人夜夜| 国产精品.久久久| 中文字幕另类日韩欧美亚洲嫩草| 国产免费福利视频在线观看| 亚洲欧美中文字幕日韩二区| 2018国产大陆天天弄谢| 少妇被粗大猛烈的视频| 久久精品国产自在天天线| 一级毛片 在线播放| 亚洲成国产人片在线观看| 亚洲成国产人片在线观看| 最新的欧美精品一区二区| 国产精品偷伦视频观看了| 日韩一本色道免费dvd| 人人妻人人澡人人看| 亚洲欧洲日产国产| 黑人欧美特级aaaaaa片| 精品第一国产精品| 免费观看无遮挡的男女| 国产精品女同一区二区软件| 91在线精品国自产拍蜜月| videossex国产| 性色avwww在线观看| 国产成人免费观看mmmm| 超色免费av| 国产野战对白在线观看| 国产色婷婷99| 最近2019中文字幕mv第一页| 亚洲欧美一区二区三区国产| 日韩一区二区三区影片| 中文字幕制服av| 亚洲欧美色中文字幕在线| 国产精品国产av在线观看| 18+在线观看网站| 日韩电影二区| 亚洲精华国产精华液的使用体验| 国产高清不卡午夜福利| 午夜日本视频在线| 成人亚洲欧美一区二区av| 你懂的网址亚洲精品在线观看| 欧美黄色片欧美黄色片| 久久久精品免费免费高清| 亚洲,欧美,日韩| 熟女电影av网| 天天操日日干夜夜撸| 91aial.com中文字幕在线观看| 国产一区二区三区av在线| 亚洲av男天堂| 久久午夜综合久久蜜桃| 日韩欧美精品免费久久| 久久久久网色| 男女国产视频网站| 成年女人毛片免费观看观看9 | 日本爱情动作片www.在线观看| 日日爽夜夜爽网站| 婷婷色综合www| 女的被弄到高潮叫床怎么办| 美女午夜性视频免费| 国产成人精品福利久久| 亚洲av日韩在线播放| 2021少妇久久久久久久久久久| 夜夜骑夜夜射夜夜干| 女性被躁到高潮视频| 晚上一个人看的免费电影| 免费黄频网站在线观看国产| 蜜桃在线观看..| 在线天堂最新版资源| 午夜av观看不卡| 国产精品人妻久久久影院| 亚洲 欧美一区二区三区| 欧美日韩av久久| 亚洲精品日本国产第一区| 久久久久久久国产电影| 波多野结衣av一区二区av| 性高湖久久久久久久久免费观看| 午夜福利在线观看免费完整高清在| 两性夫妻黄色片| av一本久久久久| 午夜久久久在线观看| videossex国产| 国产精品久久久久久久久免| 国产日韩欧美在线精品| 亚洲国产精品999| 亚洲av欧美aⅴ国产| 免费观看a级毛片全部| 欧美精品国产亚洲| 咕卡用的链子| 纯流量卡能插随身wifi吗| 少妇的逼水好多| 国产精品一区二区在线观看99| 亚洲av中文av极速乱| 七月丁香在线播放| 久久青草综合色| 一二三四在线观看免费中文在| 国产野战对白在线观看| 午夜福利,免费看| 美女大奶头黄色视频| av一本久久久久| 国产乱来视频区| 不卡视频在线观看欧美| 久久久久精品久久久久真实原创| 99九九在线精品视频| 男人爽女人下面视频在线观看| 精品人妻一区二区三区麻豆| 少妇的逼水好多| 中文天堂在线官网| 亚洲天堂av无毛| 多毛熟女@视频| av不卡在线播放| 国产野战对白在线观看| 18禁国产床啪视频网站| 亚洲国产精品一区二区三区在线| 日韩av在线免费看完整版不卡| 777米奇影视久久| 精品少妇内射三级| 亚洲成人一二三区av| 久热这里只有精品99| 亚洲视频免费观看视频| 边亲边吃奶的免费视频| 国产成人精品婷婷| 一级毛片我不卡| 欧美成人精品欧美一级黄| 国产不卡av网站在线观看| 亚洲av日韩在线播放| 国产一区亚洲一区在线观看| 亚洲av电影在线观看一区二区三区| 五月天丁香电影| 国产日韩欧美亚洲二区| 波野结衣二区三区在线| 哪个播放器可以免费观看大片| 久久综合国产亚洲精品| 女性被躁到高潮视频| tube8黄色片| 女性生殖器流出的白浆| 观看av在线不卡| 久久99热这里只频精品6学生| 一区二区三区激情视频| 成年动漫av网址| 午夜福利网站1000一区二区三区| 9色porny在线观看| 欧美成人午夜免费资源| 制服丝袜香蕉在线| 亚洲av福利一区| 激情视频va一区二区三区| 黑丝袜美女国产一区| 精品久久蜜臀av无| 亚洲五月色婷婷综合| 成年人免费黄色播放视频| 久久久久久伊人网av| 中文乱码字字幕精品一区二区三区| 久久精品国产亚洲av涩爱| 亚洲伊人久久精品综合| 最近的中文字幕免费完整| 在线天堂中文资源库| 亚洲成国产人片在线观看| 水蜜桃什么品种好| 精品一品国产午夜福利视频| 亚洲精品视频女| 国产人伦9x9x在线观看 | 性色av一级| 伊人久久国产一区二区| 老鸭窝网址在线观看| 亚洲四区av| 一个人免费看片子| 成年美女黄网站色视频大全免费| 777米奇影视久久| 久久久久精品性色| 久久久久久久大尺度免费视频| 久久久久久久大尺度免费视频| 成年人免费黄色播放视频| 欧美bdsm另类| 人人澡人人妻人| 久久久久久免费高清国产稀缺| 黄片小视频在线播放| 成人手机av| 中文字幕最新亚洲高清| 黄网站色视频无遮挡免费观看| 亚洲久久久国产精品| 69精品国产乱码久久久| 热99国产精品久久久久久7| 考比视频在线观看| 日韩av不卡免费在线播放| 欧美日韩视频高清一区二区三区二| 丝袜在线中文字幕| 午夜免费观看性视频| 伊人亚洲综合成人网| 少妇的丰满在线观看| 亚洲综合精品二区| 在线观看人妻少妇| 国产高清不卡午夜福利| 国产成人免费观看mmmm| 国产成人免费无遮挡视频| 亚洲一区中文字幕在线| 曰老女人黄片| 国产成人免费无遮挡视频| 亚洲国产色片| 人妻 亚洲 视频| 国产极品天堂在线| 久久热在线av| 精品国产超薄肉色丝袜足j| 国产精品二区激情视频| 如日韩欧美国产精品一区二区三区| 久久ye,这里只有精品| 最近2019中文字幕mv第一页| 另类亚洲欧美激情| 婷婷色麻豆天堂久久| 青春草亚洲视频在线观看| 另类精品久久| 成年女人在线观看亚洲视频| 日韩人妻精品一区2区三区| 男女国产视频网站| 韩国av在线不卡| 不卡视频在线观看欧美| 日韩不卡一区二区三区视频在线| 久久国产亚洲av麻豆专区| www.自偷自拍.com| 久久精品人人爽人人爽视色| 色吧在线观看| 我要看黄色一级片免费的| 久久狼人影院| 久久久久久久久久久久大奶| 日韩制服骚丝袜av| 性色avwww在线观看| 高清在线视频一区二区三区| 男女免费视频国产| 十分钟在线观看高清视频www| 婷婷色综合大香蕉| 亚洲 欧美一区二区三区| 18禁动态无遮挡网站| 看免费成人av毛片| 看非洲黑人一级黄片| 久久韩国三级中文字幕| 精品亚洲乱码少妇综合久久| 国产有黄有色有爽视频| 777米奇影视久久| 精品卡一卡二卡四卡免费| www.自偷自拍.com| 亚洲,欧美精品.| 亚洲色图 男人天堂 中文字幕| 久久 成人 亚洲| 97在线人人人人妻| 国产男女超爽视频在线观看| 久久久国产一区二区| 欧美成人午夜精品| 宅男免费午夜| 日本av免费视频播放| 青青草视频在线视频观看| 精品国产超薄肉色丝袜足j| 五月开心婷婷网| 国产成人免费观看mmmm| 最近2019中文字幕mv第一页| 午夜福利视频在线观看免费| 日韩人妻精品一区2区三区| 午夜免费鲁丝| 在线观看免费视频网站a站| av国产久精品久网站免费入址| 王馨瑶露胸无遮挡在线观看| 人人妻人人爽人人添夜夜欢视频| xxx大片免费视频| 久久精品国产亚洲av天美| av免费观看日本| 欧美成人午夜精品| 天天操日日干夜夜撸| 婷婷色综合www| 国产在线免费精品| 校园人妻丝袜中文字幕| 人妻 亚洲 视频| 精品久久蜜臀av无| 久久久久久人妻| 成人黄色视频免费在线看| 一区二区三区激情视频| av网站免费在线观看视频| 最近手机中文字幕大全| 大陆偷拍与自拍| 卡戴珊不雅视频在线播放| 久久鲁丝午夜福利片| 久久青草综合色| 深夜精品福利| 国产精品女同一区二区软件| 熟女av电影| 中文字幕最新亚洲高清| 赤兔流量卡办理| 精品一区二区免费观看| 女的被弄到高潮叫床怎么办| 欧美亚洲 丝袜 人妻 在线| 欧美精品亚洲一区二区| 欧美精品一区二区大全| 国产精品久久久久久精品电影小说| 肉色欧美久久久久久久蜜桃| av电影中文网址| 叶爱在线成人免费视频播放| av电影中文网址| 精品国产乱码久久久久久男人| 大片免费播放器 马上看| 伦理电影大哥的女人| 91在线精品国自产拍蜜月| 久久精品国产亚洲av涩爱| 亚洲 欧美一区二区三区| 日韩,欧美,国产一区二区三区| 校园人妻丝袜中文字幕| 99国产精品免费福利视频| 国产成人av激情在线播放| 亚洲精品一二三| 久久ye,这里只有精品| 成人午夜精彩视频在线观看| 黑人猛操日本美女一级片| 两个人免费观看高清视频| 亚洲精品在线美女| 国产av精品麻豆| 精品卡一卡二卡四卡免费| 久久久国产精品麻豆| 中文字幕最新亚洲高清| 一区二区日韩欧美中文字幕| 亚洲综合色网址| 免费少妇av软件| 亚洲成国产人片在线观看| 一个人免费看片子| 日韩一卡2卡3卡4卡2021年| 亚洲精品一区蜜桃| 国产成人免费观看mmmm| 天天躁狠狠躁夜夜躁狠狠躁| 日韩中文字幕视频在线看片| 久久狼人影院| 在线观看国产h片| 日本91视频免费播放| 大码成人一级视频| 丁香六月天网| 国产精品人妻久久久影院| 国产日韩欧美在线精品| 一边亲一边摸免费视频| 秋霞伦理黄片| 国产精品一区二区在线不卡| 亚洲,欧美,日韩| 观看美女的网站| 久久久久国产网址| 国产av国产精品国产| 亚洲中文av在线| 日日摸夜夜添夜夜爱| 亚洲久久久国产精品| 国产1区2区3区精品| 90打野战视频偷拍视频| 9191精品国产免费久久| 国产有黄有色有爽视频| 国产精品熟女久久久久浪| 老鸭窝网址在线观看| 日韩一区二区三区影片| 国产精品一国产av| 国产爽快片一区二区三区| 精品久久久久久电影网| 国产野战对白在线观看| av又黄又爽大尺度在线免费看| 少妇精品久久久久久久| 777米奇影视久久| 在线观看免费视频网站a站| 激情五月婷婷亚洲| 人人妻人人爽人人添夜夜欢视频| 啦啦啦中文免费视频观看日本| 国产精品.久久久| 亚洲欧美一区二区三区久久| 欧美xxⅹ黑人| 夜夜骑夜夜射夜夜干| 一级毛片黄色毛片免费观看视频| 亚洲av免费高清在线观看|