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

    基于δD和δ18O的青海湖流域芨芨草水分利用來源變化研究

    2015-02-13 01:09:35吳華武李小雁蔣志云鄭肖然趙殿智
    生態(tài)學(xué)報(bào) 2015年24期
    關(guān)鍵詞:芨芨草土壤水表層

    吳華武, 李小雁, 蔣志云, 李 靜, 鄭肖然, 趙殿智

    1北京師范大學(xué)資源學(xué)院, 北京 100875 2中國科學(xué)院地理科學(xué)與資源研究所, 北京 100101 3中國科學(xué)院大學(xué), 北京 100049 4青海省三角城種羊場(chǎng), 剛察 812300

    基于δD和δ18O的青海湖流域芨芨草水分利用來源變化研究

    吳華武1,*, 李小雁1, 蔣志云1, 李 靜2,3, 鄭肖然1, 趙殿智4

    1北京師范大學(xué)資源學(xué)院, 北京 100875 2中國科學(xué)院地理科學(xué)與資源研究所, 北京 100101 3中國科學(xué)院大學(xué), 北京 100049 4青海省三角城種羊場(chǎng), 剛察 812300

    水分條件是限制干旱半干旱地區(qū)植物生長重要的生態(tài)因子,為了揭示青海湖流域典型生態(tài)系統(tǒng)下芨芨草植物的水分利用來源及其如何響應(yīng)水分條件的變化,選擇了自然和干旱控制條件下芨芨草植物,通過測(cè)定芨芨草植物莖水和各潛在水源(土壤水、地下水及降水)中δD、δ18O組成,并利用直接比較分析法和多源混合模型計(jì)算芨芨草植物對(duì)土壤水的利用比例。研究結(jié)果表明:表層土壤水分和土壤水中δD、δ18O值表征出較大波動(dòng)范圍,其直接受降水和蒸發(fā)作用影響,土壤蒸發(fā)線的斜率和截距明顯小于大氣水線斜率和截距,表明土壤水中同位素組成經(jīng)歷了強(qiáng)烈的蒸發(fā)分餾過程,而芨芨草莖水中δD、δ18O值都集中分布土壤水蒸發(fā)線附近,說明芨芨草根系主要利用不同深度的土壤水。自然條件下芨芨草在生長季初期(6月)利用表層土壤水(0—10cm,45.1%),8—9月份大降水事件影響土壤含水量和同位素組成,降水入滲深度較深且芨芨草根系對(duì)土壤水分吸收的比例相差不大,表明根系在土壤含水量較高時(shí)均能吸水不同深度土壤水。在干旱控制條件下芨芨草在7月初主要利用表層土壤水(0-30cm),隨著表層土壤水分的減少,根系吸收深度轉(zhuǎn)向較深土壤層,而灌溉后表層土壤水分明顯增加,其吸收深度又轉(zhuǎn)向表層,表明芨芨草根系吸收深度能敏感地響應(yīng)土壤水分的變化。另外還發(fā)現(xiàn)芨芨草在生長季內(nèi)并未直接利用地下水。

    氫氧穩(wěn)定同位素;降水;土壤水;芨芨草;水分利用來源;青海湖流域

    芨芨草系多年生密叢性禾本科植物,它既是優(yōu)良伺養(yǎng)植物、纖維植物,又具有良好的水土保持功能,廣泛地分布于我國干旱半干旱地區(qū)如內(nèi)蒙古高原、黃土高原、青藏高原東部地區(qū)(共和盆地、青海湖盆地及柴達(dá)木盆地)[1-2],芨芨草環(huán)青海湖呈環(huán)帶狀分布,占整個(gè)青海湖流域面積(29960km2,2005)的三分之一,其多生長于海拔高度在2900—3400m之間、開闊平坦的湖積平原、河谷階地和山前緩坡。目前,在氣候變化和人類活動(dòng)的影響下芨芨草草原生態(tài)系統(tǒng)遭到嚴(yán)重破壞,面積急劇減小、生態(tài)環(huán)境惡化,從而導(dǎo)致草地生產(chǎn)力下降,物種多樣性也發(fā)生變化[3-4]。目前,大量研究主要集中在芨芨草草原生態(tài)系統(tǒng)的物種多樣性、空間分布格局及生態(tài)效益功能[4-6],研究發(fā)現(xiàn)經(jīng)圍欄封育后芨芨草群落的物種多樣性及空間分布格局在圍欄內(nèi)外具有明顯的差別,且圍欄內(nèi)生態(tài)環(huán)境得到有效地改善[6]。然而,水分條件是限制芨芨草草原生態(tài)系統(tǒng)植物生長發(fā)育的關(guān)鍵因子,降水是該生態(tài)系統(tǒng)的主要水分補(bǔ)給來源,而降水、土壤水和地下水是芨芨草的主要可能利用水源,盡管芨芨草草原生態(tài)系統(tǒng)降水以小降水事件(<5mm)為主,這不穩(wěn)定的降水被其高效率地利用,如研究發(fā)現(xiàn)干旱區(qū)多年生深根系植物具有一定的可塑性功能,而一年生淺根系植物能高效地利用短期不穩(wěn)定的降水以避免干旱脅迫[7]。在水分相對(duì)短缺和降水時(shí)空分配不均的芨芨草草原生態(tài)系統(tǒng),確定芨芨草水分利用策略對(duì)全面認(rèn)識(shí)其生長發(fā)育與生態(tài)功能及評(píng)價(jià)其對(duì)整個(gè)青海湖流域的水資源尤其重要。到目前為止,國內(nèi)外對(duì)芨芨草水分利用的方式和機(jī)理尚不明確。

    穩(wěn)定同位素技術(shù)在植物與水分關(guān)系研究中已成為強(qiáng)有力的示蹤工具,且已在一定程度上替代了傳統(tǒng)方法(如挖掘法、植物生理生態(tài)等)[8-9]。穩(wěn)定氫氧同位素組成不僅可示蹤水循環(huán)過程,還可確定不同植物間的水分利用關(guān)系,這為認(rèn)識(shí)生態(tài)系統(tǒng)水文及物質(zhì)循環(huán)過程提供新的途徑。植物通過根系吸水及水分在植物體內(nèi)運(yùn)輸過程中同位素組成不會(huì)發(fā)生分餾現(xiàn)象,因此植物組織內(nèi)穩(wěn)定同位素組成是不同水源(地表水、地下水、土壤水)的綜合反映[8],從而有效地辨別和區(qū)分植物水分利用水源。但Asbjornsen等指出如果不同土壤水中同位素信息組成具有高度相似性,這增加辨別植物水分利用水源的難度[9],但可結(jié)合植物根系分布特征、土壤含水量和土壤水勢(shì)來準(zhǔn)確地識(shí)別植物的水分利用來源[10]。

    Walter′s 首次提出在干旱半干旱地區(qū)植物水分利用方式普遍符合兩層水分利用模式假設(shè)(two-layer model hypothesis)[11],即深根系木本植物(喬木或灌木)利用穩(wěn)定的深層水源(地下水、深層土壤水)而淺根系草本植物(多年生或一年生)依靠表層不穩(wěn)定的水源(表層土壤水、降水)[11-12]。大量研究發(fā)現(xiàn)干旱半干旱環(huán)境下的荒漠樹木(如杉木、槭樹林)的水分利用策略與兩層水分利用模式相似,深根性植物依靠穩(wěn)定的水源[11-12]。但Le Roux等和Liu等發(fā)現(xiàn)無論深根還是淺根系植物都主要依賴于表層土壤水[13-14]。芨芨草根系發(fā)達(dá)、入土較深且以須根為主,根系最大深度能達(dá)到1.5m[15],經(jīng)野外調(diào)查發(fā)現(xiàn)該研究區(qū)芨芨草深度最深深度可超過1.6m,這可能與地下水潛水位接近,因此芨芨草植物往往被作為當(dāng)?shù)啬撩竦叵滤闹甘局参?,芨芨草根系分布特征與不同水分條件下其水分利用方式關(guān)系及地下水是否真正能被芨芨草所利用目前尚不明確。鑒于此,本文選擇青海湖流域沙柳河下游的芨芨草植物為研究對(duì)象,利用穩(wěn)定同位素技術(shù)對(duì)生長季內(nèi)其水分利用來源進(jìn)行初步探討,從而揭示芨芨草水分利用來源的季節(jié)變化及其對(duì)土壤水分條件的響應(yīng)關(guān)系,為合理地解釋優(yōu)勢(shì)物種芨芨草的水分利用機(jī)理及適應(yīng)極端干旱環(huán)境的策略,將為退化的芨芨草草原生態(tài)系統(tǒng)的恢復(fù)與重建提供理論依據(jù)和技術(shù)指導(dǎo)。

    1 材料與方法

    1.1 研究區(qū)概況

    研究區(qū)位于青海湖流域沙柳河下游地區(qū)三角城種羊場(chǎng)內(nèi)的芨芨草草原(圖1),在20世紀(jì)80年代芨芨草草原開始被圍欄封育,主要作為牧民的冬季和春季牧場(chǎng)。在沙柳河下游地區(qū)植被類型以芨芨草草原、河谷灌木具鱗水柏枝及農(nóng)田生態(tài)系統(tǒng)(青稞、燕麥、油菜)分布為主,芨芨草草原植物群落以芨芨草斑塊為中心周圍主要分布著冷蒿、羊草、早熟禾、唐古韭、本氏針茅等植物。土壤為栗鈣土,土壤容重為1.21g/cm3,土壤粘粒含量34.8%—45.3%,淺層土壤田間持水量29.7% —32%[16]。其深居內(nèi)陸,屬于高原大陸性氣候,氣溫較低,太陽輻射強(qiáng)烈、氣溫日差較大,且大部分區(qū)域常年處于寒冷、干旱少雨、多風(fēng)的氣候條件下,寒冷季節(jié)漫長,而溫暖季節(jié)相對(duì)短暫,四季不分明。降水量稀少且主要集中在5—9月,蒸發(fā)量大,年平均降水量在291— 579mm之間,年蒸發(fā)量達(dá)到1200—2000mm,最大凍土深度為2.88m。

    圖1 青海湖流域芨芨草草原的空間分布(A)和搭棚干旱控制實(shí)驗(yàn)(B)Fig.1 Spatial distribution of typical Achnatherum splendens steppe in Lake Qinghai watershed (A) and The photo for drought-treatment site (B)

    1.2 樣品收集及實(shí)驗(yàn)設(shè)計(jì)

    在研究區(qū)內(nèi)選擇固定的芨芨草斑塊3叢,采集芨芨草植物地下根莖部分,將芨芨草根莖樣本切成小段(4—5cm)裝入30mL的玻璃瓶中并用Parafilm 封口膜密封并放入便攜式的冰箱,用于植物水穩(wěn)定同位素 (δD)組成分析芨芨草的水分利用來源。同時(shí),在其附近用土鉆收集不同深度的土壤樣品(0—10、10—20、20—30、30—40和40—60cm),每層3個(gè)重復(fù),將土壤樣品分成兩組,一組裝入30mL的玻璃瓶中并用Parafilm封口膜密封并放入便攜式的冰箱,另一組裝入自封袋中用烘干法測(cè)定土壤重量含水量。所有樣品帶入實(shí)驗(yàn)室內(nèi)在-4℃下冷藏直至室內(nèi)抽提分析。

    自2013年4月28日至9月28日,在實(shí)驗(yàn)站三角城種羊場(chǎng)利用自制裝置收集每次降水事件的樣品,共52個(gè)次降水樣,且每次降水事件收集時(shí)間都與安裝在芨芨草草原的自動(dòng)氣象站觀測(cè)(Campbell Scientific,Inc., USA)同步,記錄時(shí)間間隔為10min。其中,液態(tài)降水直接注入采樣瓶(30mL)中;固態(tài)降水(雪或冰雹)先裝入水樣袋內(nèi),扎緊口袋,待其在室溫下融化后再倒入30mL采樣瓶中密封保存。在每次收集植物及土壤樣時(shí),地下水樣取自離采樣點(diǎn)約0.5km處飲水井且井水水位低于2.7m。所有水樣均在-4℃下冷藏。

    實(shí)驗(yàn)時(shí)間集中植物的生長季(6—9月)進(jìn)行,為了更好地揭示優(yōu)勢(shì)植物芨芨草對(duì)土壤水分條件的響應(yīng)關(guān)系,進(jìn)行了為期近1個(gè)月(7月5日至8月1日:7月5日、7月8日、7月18日、7月23日、7月27日及8月1日)的干旱控制實(shí)驗(yàn)。在6月采樣點(diǎn)附近選擇芨芨草斑塊(6叢)較密集處,圍繞其四周挖溝至60cm (60cm以下為碎石和礫石),溝壕四周用透明塑料薄膜隔開并將原土填埋以防止土壤水的側(cè)向流動(dòng)。在所選擇的芨芨草斑塊正上方用透明塑料薄膜進(jìn)行遮雨處理,四周保持通風(fēng),遮雨面積大小為3.0m (長)×2.4m(寬)×0.9m(高)。其中在7月21日(250L)和7月31日(200L)用來自井水進(jìn)行兩次灌溉分析芨芨草的水分利用響應(yīng)土壤水分的變化。自然條件下采樣時(shí)間為6月10日、7月5日、7月18日、8月9日、8月23日及9月9日

    1.3 樣品測(cè)定

    利用傳統(tǒng)的低溫真空蒸餾法從植物和土壤樣中抽提土壤水和植物水[17-18]。所有植物樣、土壤樣和水樣中氫氧穩(wěn)定同位素均在北京師范大學(xué)地表過程與資源生態(tài)重點(diǎn)實(shí)驗(yàn)室完成,利用LGR液態(tài)同位素分析儀(DLT-100,型號(hào):908-0008) 進(jìn)行測(cè)式分析,所有測(cè)量結(jié)果均為與VSMOW同位素濃度的千分差來表示:

    δD(δ18O)=(R樣-RV-SMOW-1)×1000

    (1)

    式中,R樣和RV-SMOW分別是水樣和標(biāo)樣(V-SMOW)中的氫(或氧)穩(wěn)定同位素比率R(18O /16O或D/H),其中δD和δ18O測(cè)試的誤差分別為δD≤2‰(1s),δ18O≤0.3‰(1s)。

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

    利用直接對(duì)比法和多源混合模型(IsoSource)[19]分析芨芨草根系在生長季內(nèi)吸收深度,首先根據(jù)根據(jù)各水源間同位素變化特征初步判斷植物吸水的可能土壤層次;然后,利用多源混合模型(IsoSource)來定量地分析芨芨草對(duì)各水源所利用的相對(duì)貢獻(xiàn)率,模型運(yùn)行時(shí)增幅和質(zhì)量平衡公差分別設(shè)置為1%和0.1%,其中質(zhì)量平衡公差一般不小于來源增量與各可能水源同位素值之間最大差值的乘積的1/2。采用SPSS13.0中的Tukey法對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行分析,在0.05水平上檢驗(yàn)其顯著性。

    2 結(jié)果

    2.1 降水量、溫度及降水中δD和δ18O的變化特征

    對(duì)降水中δD、δ18O分析結(jié)果(圖2a)表明,降水中δD和δ18O存在較大的波動(dòng)范圍,δD的波動(dòng)范圍為:-96.9‰—45.2‰,加權(quán)平均值為-38.6‰;δ18O的波動(dòng)范圍為-14.1‰— 4.36‰,加權(quán)平均值為-6.71‰。其中在8、9月份降水中穩(wěn)定氫氧穩(wěn)定同位素組成較其他月份貧化,但小降水事件降水中δD、δ18O值出現(xiàn)正值的情況,這可能與小降水事件的雨滴在下降過程中受到強(qiáng)烈的蒸發(fā)過程[20-21]。在大降水事件(或雪樣)中穩(wěn)定同位素組成較貧化,這主要與降水的形成過程密切有關(guān)[20-21],如4月28日(雪:δ18O=-14.1‰,δD=-96.9‰)、5月17日(雨夾雪:δ18O=-12.8‰,δD=-90.3‰)和6月9日(雨夾雪:δ18O=-10.4‰,δD=-66.1‰)(圖2a)。

    青海湖流域?qū)儆谖覈湫偷母珊蛋敫珊?、高寒地區(qū),雨熱同期,降水主要集中5—9月且占全年的85% (圖2b)。2013年雨季(5—9月)降水量為350.6mm,大降水事件集中在8月份,8月20日—23日期間降水量為89.4mm,占雨季的25.6%,該月降水量在全年中最多,為129.1mm(圖2c)。在6月份有3次單場(chǎng)大降水事件(10mm)(如6月8、9、21日),而在7月份小降水事件(<10mm) 分布較多。溫度在2—17℃之間呈現(xiàn)波動(dòng)變化(圖2c)。

    圖2 芨芨草草原區(qū)溫度、降水量及降水中δD和δ18O的變化特征Fig.2 The variation characteristics of air temperature, precipitation amount, δD and δ18O of precipitation in the Achnatherum splendens steppe areaa: δD和δ18O;b: 多年平均降水量和溫度(1957—2013);c: 降水量和溫度;三角形表示雪樣;五角星表示雨夾雪樣;箭頭表示植物和土壤樣采集時(shí)間

    Craig[22]首次將大氣降水中δD與δ18O之間的關(guān)系稱為大氣水線,它對(duì)揭示水循環(huán)過程中穩(wěn)定同位素變化具有重要意義。在全球尺度下,全球大氣水線(GMWL)為:δD=8δ18O+10。不同區(qū)域大氣水線往往偏離全球大氣水線,此能較好地反映出不同地區(qū)降水變化規(guī)律。由于大氣降水、土壤水、植物水、地下水等在時(shí)空尺度上存在著一定的水力聯(lián)系,將各水體中氫氧穩(wěn)定同位素組成與當(dāng)?shù)卮髿馑€進(jìn)行對(duì)比分析具有重要意義。本文對(duì)該地區(qū)降水中δD與δ18O組成進(jìn)行線性回歸分析,得出青海湖流域大氣水線方程(LMWL)(圖3):

    δD=8.07δ18O+14.8

    (2)

    該大氣水線的斜率和截距與全球水線的相接近。根據(jù)在自然條件(式3, LASWL)和干旱控制條件(式4, LCSWL)下土壤水中δD與δ18O組成分析得出兩種不同條件下的土壤蒸發(fā)線:

    δD=4.88δ18O-21.6(n=30, r2=0.849)

    (3)

    δD=5.14δ18O-21.3(n=37, r2=0.923)

    (4)

    圖3 大氣降水、土壤水及植物水中δD-δ18O關(guān)系 Fig.3 The regression linear relationships between δD and δ18O in precipitation, soil and plant stem water

    從圖3可以看出,土壤水和植物水中δD、δ18O組成都位于大氣水線右下方,且在兩種生境條件下土壤蒸發(fā)線的斜率和截距相差不大,分別為4.88和5.14且介于4—6之間,它們明顯低于大氣水線的斜率和截距,表明土壤水中δD、δ18O經(jīng)歷了不同程度的二次蒸發(fā)分餾效應(yīng)使得δD、δ18O 發(fā)生富集作用[17,23-25]。植物莖水中δD、δ18O值都圍繞看土壤水分布(圖3),這表明芨芨草主要利用不同深度的土壤水[13]。

    2.2 土壤水分及土壤水中δD變化

    圖4 土壤水分、土壤水及植物水中δD的變化Fig.4 The variations of soil water content (SWC) and δD in soil and plant water

    土壤水分的時(shí)空異質(zhì)性不僅與降水和蒸發(fā)作用有關(guān),還與植被對(duì)水分吸收利用有關(guān)[26]。從圖4可以看出,兩種生境條件下土壤水分隨時(shí)間和深度的變化明顯,表層土壤含水量(0—20cm)變化明顯,而深層土壤含水量變化不明顯。在自然條件下,表層土壤水分(0—10cm)在13.2%— 34.5%之間波動(dòng)變化范圍明顯,在8月9日前有12d無降水事件,干旱時(shí)間長,從而土壤水分含量最低。在8月23日和9月9日受降水補(bǔ)給的影響,表層土壤含水量明顯增大,其中在8月23日前5d連續(xù)降水且降水量大(圖2a),降水量達(dá)86.3mm,且降水入滲深度較深,在40cm以上土壤水分受降水影響顯著。但在9月9日前有連續(xù)11d的降水過程,降水時(shí)間長且以小降水事件(<10mm)為主,土壤水分含量明顯比其他時(shí)間高。在干旱控制條件下,由于不受外界降水的影響,土壤水分含量減小明顯,在7月5日至18日期間,表層土壤含水量(0—10cm) 從21.5%減少到8.51%(圖4)。在7月21日和7月31日進(jìn)行了兩次灌溉后使得在7月23日和8月1日表層土壤含水量增大,但深層土壤水變化不明顯。

    圖4給出了兩種生境下土壤水中δD的變化特征,淺層土壤水中δD變化比深層土壤水中的δD變化明顯。在自然條件下,表層土壤水(0—10cm)中δD的差異性顯著(P<0.05),表現(xiàn)出季節(jié)性變化特征,在-14.3‰— -51.4‰中波動(dòng)變化,平均值-29.9‰。表層土壤水中穩(wěn)定同位素值明顯比深層富集,因?yàn)楸韺油寥浪艿綇?qiáng)烈的二次蒸發(fā)分餾作用從而使得其較富集于深層土壤水[23]。但在8月23日和9月9日土壤水中δD的變化呈S型變化(圖4),表明直接受降水中穩(wěn)定同位素組成影響,表層土壤水中δD與降水中同位素組成相接近,且土壤水中δD在垂直剖面上差異性不大。

    在干旱控制條件下,土壤水中穩(wěn)定同位素組成明顯富集于自然條件下,在7月5日至18日期間表層土壤水中δD從-23‰ 增大到-8.13‰(圖4),同樣表層土壤水中穩(wěn)定同位素組成受到強(qiáng)烈的蒸發(fā)作用,但在7月23日和8月1日土壤水中δD值較小,其直接受7月21日和7月31日灌溉水的影響,在7月27日表層土壤水中穩(wěn)定同位素組成經(jīng)歷了二次蒸發(fā)分餾作用使得δD比7月21日的δD更富集。然而,在20cm以下土壤水分同位素差異性較小,這是由于在此深度以下受蒸發(fā)作用逐漸減弱。

    2.3 芨芨草植物中δD變化特征及其對(duì)土壤水分利用比例

    在兩種生境條件下芨芨草植物中穩(wěn)定氫氧同位素(δD和δ18O)組成位于全球大氣水線(GMWL)和青海湖流域大氣水線(LMWL)右下側(cè)區(qū)域(圖3),表明植物所利用的水源受到一定程度蒸發(fā)作用影響,使得穩(wěn)定氫氧同位素發(fā)生分餾作用而重同位素富集。對(duì)兩種生境條件下芨芨草莖水中δD組成對(duì)比分析發(fā)現(xiàn):自然條件下芨芨草莖水在-18.5‰— -40.9‰ 之間波動(dòng)變化且其介于10—60cm土壤水之間,平均值為-31.4‰,表明芨芨草在整個(gè)生長季內(nèi)水分可能主要來源于不同深度的土壤水;在干旱控制條件下芨芨草莖水中δD的波動(dòng)變化范圍(-16.8‰ —-43.2‰,平均值:-31.6‰)比自然條件下的波動(dòng)范圍要大。芨芨草莖水中δD組成與地下水中δD組成存在一定的差異(圖4),表明芨芨草植物并沒有利用地下水。

    自然條件下芨芨草莖水中δD組成表現(xiàn)出季節(jié)變化特征(圖(4),芨芨草莖水中δD組成在生長初期(6月)與0—20cm的土壤水中δD組成相似,且土壤水分利用比例達(dá)65.1%(表 1)。隨著季節(jié)的變化,芨芨草對(duì)不同深度土壤水分利用比例也存在差異,在7月18日芨芨草對(duì)表層0—10cm的土壤水利用比例高達(dá)94.2%,且芨芨草莖水中δD值與降水中δD值比較接近,這表明芨芨草根系利用表層土壤水。在生長季末(8—9月)持續(xù)降水影響,表層土壤含水量明顯增大,芨芨草根系對(duì)淺層土壤水(0—20cm)利用比例分別為48.7%(8月24日)和36.9%(9月9日)。

    在干旱控制條件下芨芨草水分利用深度能較敏感地響應(yīng)土壤水分的變化(表1),在7月5日芨芨草根系吸水深度主要分布在0—30cm(77.8%),隨著干旱程度的加劇,芨芨草的吸水深度逐漸向深層土壤水轉(zhuǎn)移,在7月8日芨芨草利用表層土壤水,利用比例達(dá)67.9%;而在7月18日根系吸水深度在40—60cm且利用比例最高(57.4%)。7月23日和8月1日水源由深層土壤水(40—60cm)轉(zhuǎn)換到淺層土壤水(0—20cm,65.9%),這直接受灌溉水的影響使得表層根系水分利用增加,由此可以看出,芨芨草植物根系吸水功能具有可塑性,其能敏感地響應(yīng)土壤有效含水量的變化以改變水分利用方式,從而達(dá)到避免水分脅迫。

    3 討論

    3.1 植物各潛在水源的變化

    在干旱半干旱地區(qū)植物所利用的潛在水源主要有降水、土壤水、地下水和河水等,其中降水是土壤水和地下水主要補(bǔ)給水源,因此大氣降水中穩(wěn)定同位素信息能較好地反映土壤水中穩(wěn)定同位素信息[27]。青海湖流域大氣降水中δD、δ18O具有明顯季節(jié)變化,6—7月份降水中δD、δ18O值較高,而在其他月份降水中δD、δ18O值較低(圖2),這與不同水汽來源密切相關(guān),在6—7月降水中穩(wěn)定同位素偏高,這可能受到兩種因素的影響:(1)在此期間降水受到局地蒸發(fā)水汽補(bǔ)充,如青海湖大水體,該湖泊水穩(wěn)定同位素組成較富集,如青海湖湖水中的平均δD、δ18O值分別為10.0‰ 和1.97‰[28]。(2)受干熱的氣象條件的影響使得降落的雨滴產(chǎn)生強(qiáng)烈的蒸發(fā)作用,此變化特征與全球中高緯度地區(qū)降水中穩(wěn)定同位素變化特征相一致[29],而在4—5月和8—9月可能受來自北冰洋高緯度極地氣團(tuán)和季風(fēng)后期的水汽影響[30],其穩(wěn)定同位素含量較低。而4—6月固體降水(雪、雨夾雪)中穩(wěn)定同位素值含量較低,因?yàn)楣腆w降水在形成和下降過程中受二次蒸發(fā)分餾作用較弱[21]。

    表1 在不同生境條件下芨芨草對(duì)不同深度土壤水的利用比例(平均值(最小值—最大值)(月-日)

    Table 1 The plant-use water proportions from different depths of soil moisture in two habitats for Achnatherum splendens(mean (minimum—maximum)) (month-day)

    深度Depth/cm06-1007-0507-1808-0908-2309-09自然1045.1(21—63)22.4(0—77)94.2(93—96)19.6(0—59)23.2(0—68)20.3(0—55)Control2020(0—79)24.6(0—88)1.5(0—5)17.7(0—51)25.5(0—93)16.6(0—49)3013.8(0—55)25.3(0—99)1.9(0—7)22.5(0—74)12.6(0—31)15.5(0—45)4012.1(0—47)20.3(0—68)1.5(0—6)19.7(0—48)15.4(0—40)25.2(0—97)609.1(0—36)7.4(0—23)0.8(0—4)20.5(0—51)23.4(0—69)22.4(0—65)深度Depth/cm07-0507-0807-1807-2307-2708-01控制1025.8(0—78)67.9(49—77)5.5(0—22)33.3(0—67)3(0—12)8.1(0—33)Treatment2026.4(0—83)13(0—51)9.6(0—38)32.6(0—73)4(0—16)62.2(47—73)3025.6(0—98)6.8(0—27)10.5(0—42)13.9(0—45)12.9(0—50)13.2(0—53)4014.5(0—43)6.8(0—26)17(0—66)10.1(0—32)21.7(0—82)9.9(0—40)607.7(0—22)5.5(0—22)57.4(34—77)10.1(0—32)58.4(18—87)6.6(0—27)

    土壤水中δD組成主要受降水、蒸發(fā)、地下水等因素影響,在不同生境條件、時(shí)間、土層深度下土壤水中δD組成具有較大的差異性[25,31]。表層土壤水分和土壤水δD組成能敏感地響應(yīng)降水的變化,如在自然和干旱控制條件下表層土壤水含量和土壤水δD組成具有較大的差異性,7月份降水大部分以小降水事件為主且蒸發(fā)作用強(qiáng)烈,土壤水分低,δD發(fā)生強(qiáng)烈富集作用;另外,小降水事件受到二次蒸發(fā)作用影響使得其δD值偏大,因此表層土壤水中穩(wěn)定同位素值大,深層土壤水中δD值明顯小于表層的;而在8—9月受強(qiáng)降水的影響,表層土壤含水量高且δD值較低,這一研究結(jié)果與周海等[18]、田立德等[27]及Saxena[32]研究結(jié)果一致,研究發(fā)現(xiàn)表層土壤水中穩(wěn)定同位素組成直接受大氣降水的影響,且其與降水中穩(wěn)定同位素組成有相似的變化趨勢(shì),而深層土壤水受降水和蒸發(fā)作用影響小。土壤水的入滲深度與降水強(qiáng)度、歷時(shí)及土壤質(zhì)地、植物根系特征有關(guān),如在8月23日前期降水量大(t=5d,P=86.3mm)和9月9日前降水持續(xù)時(shí)間長(t=11d,P=30.7mm)(圖2a),這使得大量降水入滲到土壤層中,進(jìn)而影響著深層土壤水中δD值和土壤含水量。另外,土壤水穩(wěn)定同位素受蒸發(fā)分餾作用的影響,隨著土壤深度增加蒸發(fā)分餾作用逐漸減弱[23],因此土壤水中δD表征出隨著土壤深度增加而不斷減小(圖4)。在兩種實(shí)驗(yàn)條件下深層地下水和土壤水(60cm)的穩(wěn)定同位素組成的差異性不顯著,地下水中δD的波動(dòng)范圍小,這說明地下水受當(dāng)年降水影響較小。但8月23日深層(60cm)土壤水中δD值較小,這除了與降水特性有關(guān),還可能與芨芨草根系分布特征密切相關(guān),其使得降水快速通過優(yōu)勢(shì)流達(dá)到此土壤層[33]。

    3.2 兩種生境條件下芨芨草水分利用來源的比較分析

    兩種不同生境下芨芨草具有相同的水分利用方式,均主要依靠不同深度的土壤水,但它們對(duì)水分條件變化響應(yīng)程度具有顯著差異(圖4,表1)。芨芨草為深根性系禾本科草本植物,吸收根范圍主要分布在0—60cm,根系集中分布于表層土壤,分布在表層土壤中的根系能吸收降水補(bǔ)給土壤水,而深層根系主要利用冬春季降水或者地下水補(bǔ)給的土壤水[34]。在自然條件下芨芨草根系所吸收水分主要來自表層土壤水(0—30cm), 7月18日表層土壤含水量增加使得芨芨草根系吸水深度轉(zhuǎn)變?yōu)闇\層,進(jìn)而隨著土壤表層土壤水分的降低,其根系吸水深度逐漸轉(zhuǎn)向較深層土壤水(圖4),但這種水分利用方式在干旱控制條件下的芨芨草更顯著(圖4),土壤水分的變化會(huì)引起不同深度根系的活性發(fā)生變化,干旱控制條件下芨芨草灌溉后,表層根系的活性增強(qiáng)使得芨芨草表層根系吸收水分增加[35],從而表明芨芨草對(duì)干旱環(huán)境具有極強(qiáng)的適應(yīng)能力??梢钥闯觯杠覆莞祵?duì)土壤水分具有較強(qiáng)的可塑性,說明其根系吸水特征與干旱半干旱區(qū)植物(灌木、喬木)的二形態(tài)根系功能類似,在濕季側(cè)根利用淺層土壤水而干季主根吸收深層土壤水或地下水,這有利于植物避免受水分脅迫干擾[7,36]。在干旱稀樹草原地區(qū)樹木與草本的植物群落主要通過根系的可塑性來避免水分競爭,從而有利于植物在水分短缺的條件下共生生長[11]。深根芨芨草植物也能敏感地響應(yīng)水分條件變化來改變其水分利用方式以更好地適應(yīng)多變的水分環(huán)境,Asbjornsen等在愛荷華州中部地區(qū)研究發(fā)現(xiàn)須芒草(A.gerardii)系深根植物同樣在干濕季節(jié)通過根系的可塑性功能來改變根系吸水深度,使得深根植物在干旱環(huán)境下能高效率地利用有效水分和養(yǎng)分,從而其比淺根伴生植物具有較強(qiáng)的水分和養(yǎng)分競爭優(yōu)勢(shì)[9]。同樣,芨芨草根系發(fā)達(dá)且集中分布于表層土壤,這有利于其高效率地利用不穩(wěn)定的水分,這可能是芨芨草成為環(huán)青海湖生態(tài)系統(tǒng)優(yōu)勢(shì)植物之一。

    地下水往往是干旱半干旱區(qū)植物生長的穩(wěn)定水源,這為旱季植物能有效地避免因干旱脅迫而影響其正常生長[37],如在西澳大利亞的原生樹種(Banksiaprionotes)、青海湖流域河岸植物具鱗水柏枝(Myricariasquamosa)及荒漠生態(tài)系統(tǒng)中的白刺(Nitrariatangutorum)均在旱季以地下水作為主要的可利用水源[18,30,36-37]。而芨芨草植物無論在干季還是濕季均不以地下水作為水源(圖4),這與芨芨草作為當(dāng)?shù)啬撩駥ふ业叵滤吹闹甘局参锝Y(jié)論不一致。芨芨草主要生長在地下水位較淺的河岸沖積平原,土壤層較薄,在60cm以下多以碎石和礫石分布為主,這可能在一定程度上隔斷了以毛管上升水形式向上運(yùn)動(dòng)供植物利用,這可能影響分析了芨芨草根系吸收水分的范圍。鑒于此,需下一步控制實(shí)驗(yàn)向深層土壤中注入重氘或氧同位素來示蹤芨芨草根系是否利用地下水。

    4 結(jié)論

    基于氫氧穩(wěn)定同位素技術(shù)分析了青海湖流域典型芨芨草生態(tài)系統(tǒng)中優(yōu)勢(shì)植物芨芨草在不同生境條件下的水分利用來源變化特征,揭示了芨芨草植物的水分利用方式對(duì)水分條件響應(yīng)關(guān)系。研究結(jié)果表明:

    (1)土壤含水量及其穩(wěn)定同位素組成具有明顯的季節(jié)變化特征,在兩種生境條件下,表層土壤含水量和穩(wěn)定同位素組成波動(dòng)范圍比深層土壤顯著,表明表層土壤受蒸發(fā)作用和降水的影響比深層土壤大;

    (2)土壤水中δD、δ18O值明顯偏離青海湖流域大氣水線,表明了其經(jīng)歷了蒸發(fā)分餾作用,芨芨草植物水中δD、δ18O值都集中沿著土壤水蒸發(fā)線附近分布,這表征了芨芨草根系吸水主要來自不同深度的土壤水;

    (3)利用同位素多源混合模型(IsoSource)計(jì)算自然和干旱控制實(shí)驗(yàn)條件下芨芨草根系對(duì)土壤水分的利用情況,表明在生長季初期(6月)自然條件下芨芨草根系主要利用表層土壤水(0—10cm,45.1%),8—9月份由于受降水影響,土壤水分及其同位素組成差異性較小且根系對(duì)土壤水分吸收的比例相差不大,表明在土壤含水量較高時(shí)不同深度的根系均能吸水一定量的土壤水。而干旱控制實(shí)驗(yàn)條件下芨芨草根系吸水深度隨著土壤水分的變化而發(fā)生轉(zhuǎn)變,說明芨芨草根系吸水功能與其他灌木相似具有較強(qiáng)的可塑性,從而有利于其避免水分脅迫。

    [1] 吳珍蘭, 盧生蓮.論世界芨芨草屬(禾本科)的地理分布.植物分類學(xué)報(bào), 1995, 34(2): 152-161.

    [2] 陳桂琛, 彭敏.青海省芨芨草草原的群落特征及其分布規(guī)律.西北植物學(xué)報(bào), 1993, 13(2): 154-162.

    [3] Huai H Y, Wei W H, Zhang Y L.Characteristics of theAchnatherumsplendenscommunity along the Qinghai-Tibet Railway, China.Frontiers of Biology in China, 2008, 3(4): 477-483.

    [4] 周國英, 陳桂琛, 趙以蓮, 王順忠, 李 偉, 孫 菁, 彭 敏.青海湖地區(qū)芨芨草群落特征及其物種多樣性研究.西北植物學(xué)報(bào), 2003, 23(11): 1956-1962.

    [5] 周國英, 陳桂琛, 魏國良, 韓友吉, 祝存冠.青海湖地區(qū)芨芨草群落主要種群分布格局研究.西北植物學(xué)報(bào), 2006, 26(3): 579-584.

    [6] 周國英, 陳桂琛, 徐文華, 楊路存, 韓友吉, 李錦萍.圍欄封育對(duì)青海湖地區(qū)芨芨草草原生物量的影響.干旱區(qū)地理, 2010, 33(3): 434-441.

    [7] Ehleringer J R, Dawson T E.Water uptake by plants: perspectives from stable isotope composition.Plant, Cell &Environment, 1992, 15(9): 1073-1082.

    [8] Brunel J P, Walker G R, Kennett-Smith A K.Field validation of isotopic procedures for determining sources of water used by plants in a semi-arid environment.Journal of Hydrology, 1995, 167(1/4): 351-368.

    [9] Asbjornsen H, Shepherd G, Helmers M, Mora G.Seasonal patterns in depth of water uptake under contrasting annual and perennial systems in the Corn Belt Region of the Midwestern US.Plant and Soil, 2008, 308(1/2): 69-92.

    [10] Li S G, Tsujimura M, Sugimoto A, Sasaki L, Yamanaka T, Davaa G, Oyunbaatar D, Sugita M.Seasonal variation in oxygen isotope composition of waters for a montane larch forest in Mongolia.Trees, 2006, 20(1): 122-130.

    [11] Ward D, Wiegand K, Getzin S.Walter′s two-layer hypothesis revisited: back to the roots!.Oecologia, 2013, 172(3): 617-630.

    [12] McCole A A, Stern L A.Seasonal water use patterns ofJuniperusasheion the Edwards Plateau, Texas, based on stable isotopes in water.Journal of Hydrology, 2007, 342(3/4): 238-248.

    [13] Le Roux X, Bariac T, Mariotti A.Spatial partitioning of the soil water resource between grass and shrub components in a West African humid savanna.Oecologia, 1995, 104(2): 147-155.

    [14] Liu W J, Wang P Y, Li J T, Liu W Y, Li H M.Plasticity of source-water acquisition in epiphytic, transitional and terrestrial growth phases ofFicustinctoria.Ecohydrology, 2014, doi: 10.1002/eco.1475.

    [15] 王 珺, 劉茂松, 盛 晟, 徐 馳, 劉小愷, 王漢杰.干旱區(qū)植物群落土壤水鹽及根系生物量的空間分布格局.生態(tài)學(xué)報(bào), 2008, 28(9): 4120-4127.

    [16] 朱寶文, 陳曉光, 鄭有飛, 李英年.青海湖北岸天然草地小尺度地表徑流與降水關(guān)系.冰川凍土, 2009, 31(6): 1074-1079.

    [17] 王 鵬, 宋獻(xiàn)方, 袁瑞強(qiáng), 韓冬梅, 張應(yīng)華, 張兵.基于氫氧穩(wěn)定同位素的華北農(nóng)田夏玉米耗水規(guī)律研究.自然資源學(xué)報(bào), 2013, 28(3): 481-491.

    [18] 周 海, 鄭新軍, 唐立松, 李 彥.準(zhǔn)噶爾盆地東南緣多枝檉柳、白刺和紅砂水分來源的異同.植物生態(tài)學(xué)報(bào), 2013, 37(7): 665-673.

    [19] Phillips D L, Gregg J W.Source partitioning using stable isotopes: coping with too many sources.Oecologia, 2003, 136(2): 261-269.

    [20] Peng H, Mayer B, Harris S, Roy K H.A 10-yr record of stable isotope ratios of hydrogen and oxygen in precipitation at Calgary, Alberta, Canada.Tellus B, 2004, 56(2): 147-159.

    [21] Wu H W, Zhang X P, Li X Y, Li G, Huang Y M.Seasonal variations of deuterium and oxygen-18isotopes and their response to moisture source for precipitation events in the subtropical monsoon region.Hydrological Processes, 2015, 29(1): 90-102.

    [22] Craig H.Isotope variations in meteoric waters.Science, 1961, 133(3465): 1702-1703.

    [23] Barnes C J, Allison G B.Tracing of water movement in the unsaturated zone using stable isotopes of hydrogen and oxygen.Journal of Hydrology, 1988, 100(1/3): 143-176.

    [24] Gibson J J, Edwards T W D.Regional water balance trends and evaporation-transpiration partitioning from a stable isotope survey of lakes in northern Canada.Global Biogeochemical Cycles, 2002, 16(2): 10-1-10-14, doi: 10.1029/2001GB001839.

    [25] Midwood A J, Boutton T W, Archer S R, Watts S E.Water use by woody plants on contrasting soils in a savanna parkland: assessment withδ2H andδ18O.Plant and Soil, 1998, 205(1): 13-24.

    [26] 王 軍, 傅伯杰, 蔣小平.土壤水分異質(zhì)性的研究綜述.水土保持研究, 2002, 9(1): 1-5.

    [27] 田立德, 姚檀棟, Tsujimura M, 孫維貞.青藏高原中部土壤水中穩(wěn)定同位素變化.土壤學(xué)報(bào), 2002, 39(3): 289-295.

    [28] Zhang P, Zhang B, Yang W.Environmental evolution of the water body of Qinghai Lake since the postglacial age.Chinese Journal of Geochemistry, 1989, 8(2): 112-125.

    [29] 章新平, 劉晶淼, 孫維貞, 黃一民, 張劍明.中國西南地區(qū)降水中氧穩(wěn)定同位素比率與相關(guān)氣象要素之間關(guān)系的研究.中國科學(xué)(D輯), 2006, 36(9): 850-859.

    [30] An Z S, Colman S M, Zhou W, Li X Q, Brown E T, Jull A J T, Cai Y J, Huang, Y S, Lu X F, Chang H, Song Y, Sun Y, Xu H, Liu W, Jin Z, Liu X, Cheng P, Liu Y, Ai L, Li X, Liu X, Yan L, Shi Z, Wang X, Wu F, Qiang X, Dong J, Lu F, Xu X.Interplay between the Westerlies and Asian monsoon recorded in Lake Qinghai sediments since 32ka.Scientific Reports, 2012, 2: 619, doi: 10.1038/srep00619.

    [31] Rossatto D R, Silva L D R, Villalobos-Vega R, Sternberg L S L, Franco A C.Depth of water uptake in woody plants relates to groundwater level and vegetation structure along a topographic gradient in a neotropical savanna.Environmental and Experimental Botany, 2012, 77: 259-266.

    [32] Saxena R K.Seasonal variations of oxygen-18in soil moisture and estimation of recharge in esker and moraine formations.Nordic Hydrology, 1984, 15(4): 235-242.

    [33] Tang K L, Feng X H.The effect of soil hydrology on the oxygen and hydrogen isotopic compositions of plants′ source water.Earth and Planetary Science Letters, 2001, 185(3/4): 355-367.

    [34] Chimner R A, Cooper D J.Using stable oxygen isotopes to quantify the water source used for transpiration by native shrubs in the San Luis Valley, Colorado U.S.A.Plant and Soil, 2004, 260(1/2): 225-236.

    [35] Williams D G, Ehleringer J R.Intra-and interspecific variation for summer precipitation use in pinyon-juniper woodlands.Ecological Monographs, 2000, 70(4): 517-537.

    [36] Dawson T E, Pate J S.Seasonal water uptake and movement in root systems of Australian phraeatophytic plants of dimorphic root morphology: a stable isotope investigation.Oecologia, 1996, 107: 13-20.

    [37] 趙國琴, 李小雁, 吳華武, 張思毅, 李廣泳.青海湖流域具鱗水柏枝植物水分利用氫同位素示蹤研究.植物生態(tài)學(xué)報(bào), 2013, 37(12): 1091-1100.

    Variations in water use forAchnatherumsplendensin Lake Qinghai watershed, based onδD andδ18O

    WU Huawu1,*, LI Xiaoyan1, JIANG Zhiyun1, LI Jing2, 3, ZHENG Xiaoran1, ZHAO Dianzhi4

    1CollegeofResourcesScienceandTechnology,BeijingNormalUniversity,Beijing100875,China2InstituteofGeographicSciencesandNaturalResourcesResearch,CAS,Beijing100101,China3UniversityofChineseAcademyofSciences,Beijing100049,China4SanjiaochengSheepBreedingFarm,Gangcha812300,China

    Water availability is a key ecological factor for plant growth in arid and semi-arid regions.We identified the characteristics of water-uptake patterns byAchnatherumsplendensand how it responds to changes in water conditions in typical steppe ecosystem (Lake Qinghai watershed).We analyzed stable isotopes of plant stem water and potential water sources including soil water, groundwater, and precipitation via direct comparison with differences among potential water sources and a multi-source mixing model calculating the plant-use proportions from diverse soil depths under natural and drought treatment areas.The results showed that there were large fluctuations in the soil moisture and the isotopic contents (δD) of soil water in the surface soil layers, directly influenced by precipitation and evaporation.The slope and intercept of the evaporative soil water line were evidently smaller than those of the meteoric water line, indicating that the isotopic contents from soil water experienced stronger evaporation fractionation processes.The isotopic contents (δD andδ18O) of the stem were scattered along the soil evaporative water line indicating that the possible water sources were extracted from diverse soil water depths.A.splendensabsorbed mainly surface soil water (the proportion from 0—10cm was 45.1% in June) at the early growth stage.Large precipitation events significantly affected the soil water content and isotopic contents and infiltrated into deep soil layers.Differences in plant water use proportions were minimal during August and September, which revealed that the roots ofA.splendensextracted water from various depths when soil water availability was high.Under the drought treatment conditions,A.splendensused mostly surface soil water (0—30cm) in early July.The absorption water depth shifted to deep soil layers when the surface soil water became drier.However, the absorption depth again shifted from deep to surface soil layers when surface soil moisture sharply increased after irrigation.Groundwater was not a main water source forA.splendensduring the entire growing season.

    stable isotopes;precipitation;soil water;Achnatherumsplendens;water-use sources;Lake Qinghai watershed

    國家自然科學(xué)基金重點(diǎn)項(xiàng)目(41130640);長江學(xué)者和創(chuàng)新團(tuán)隊(duì)計(jì)劃( IRT1108) ;中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資助

    2014-06-23; < class="emphasis_bold">網(wǎng)絡(luò)出版日期:

    日期:2015-05-21

    10.5846/stxb201406231300

    *通訊作者Corresponding author.E-mail: wuhuawu416@163.com

    吳華武, 李小雁, 蔣志云, 李靜, 鄭肖然, 趙殿智.基于δD和δ18O的青海湖流域芨芨草水分利用來源變化研究.生態(tài)學(xué)報(bào),2015,35(24):8174-8183.

    Wu H W, Li X Y, Jiang Z Y, Li J, Zheng X R, Zhao D Z.Variations in water use forAchnatherumsplendensin Lake Qinghai watershed, based onδD andδ18O.Acta Ecologica Sinica,2015,35(24):8174-8183.

    猜你喜歡
    芨芨草土壤水表層
    竹子和芨芨草
    半潛式平臺(tái)表層卡套管處理與認(rèn)識(shí)
    海洋石油(2021年3期)2021-11-05 07:43:10
    水體表層沉積物對(duì)磷的吸收及釋放研究進(jìn)展
    改進(jìn)的PSO-RBF模型在土壤水入滲參數(shù)非線性預(yù)測(cè)中的應(yīng)用研究
    錦州市土壤水動(dòng)態(tài)過程及影響因素
    灌水定額對(duì)土壤水鹽分布及作物產(chǎn)量的影響
    甘肅省甘谷縣第三中學(xué)芨芨草文學(xué)社作品小輯
    芨芨草文學(xué)社
    大作文(2016年11期)2016-05-30 00:25:54
    氬弧熔覆原位合成Ti(C,N)-WC增強(qiáng)鎳基表層復(fù)合材料的研究
    焊接(2015年6期)2015-07-18 11:02:25
    超聲波光整強(qiáng)化40Cr表層顯微硬度研究
    满18在线观看网站| 非洲黑人性xxxx精品又粗又长| 黄色片一级片一级黄色片| 大香蕉久久成人网| 国产熟女午夜一区二区三区| 国产麻豆成人av免费视频| 特大巨黑吊av在线直播 | 老熟妇乱子伦视频在线观看| 欧美成人午夜精品| 欧美精品啪啪一区二区三区| 成年版毛片免费区| 日韩欧美 国产精品| 亚洲精品一区av在线观看| 丰满的人妻完整版| 极品教师在线免费播放| 夜夜爽天天搞| 午夜福利视频1000在线观看| 老司机午夜十八禁免费视频| 欧美成人免费av一区二区三区| 看免费av毛片| 国产亚洲精品综合一区在线观看 | 久久久久国产精品人妻aⅴ院| 老熟妇乱子伦视频在线观看| 国产精品一区二区精品视频观看| 欧美又色又爽又黄视频| 欧美激情 高清一区二区三区| 欧美日韩亚洲综合一区二区三区_| a在线观看视频网站| 国产精品 国内视频| 国产三级黄色录像| 男女下面进入的视频免费午夜 | 国产av不卡久久| 国产一区二区激情短视频| 宅男免费午夜| 国产主播在线观看一区二区| 日韩成人在线观看一区二区三区| 美女大奶头视频| 国产av在哪里看| 久久久久久免费高清国产稀缺| 亚洲中文字幕日韩| 久久久久久免费高清国产稀缺| 男人舔奶头视频| 亚洲片人在线观看| 亚洲av日韩精品久久久久久密| 免费人成视频x8x8入口观看| 国产极品粉嫩免费观看在线| 欧美一级毛片孕妇| 免费在线观看日本一区| 女性被躁到高潮视频| 久久婷婷成人综合色麻豆| 美女午夜性视频免费| 久久久水蜜桃国产精品网| 亚洲七黄色美女视频| 夜夜爽天天搞| 一级作爱视频免费观看| 久久久国产成人精品二区| 国产激情久久老熟女| 日本在线视频免费播放| 真人做人爱边吃奶动态| 一级黄色大片毛片| 村上凉子中文字幕在线| 国产欧美日韩精品亚洲av| 午夜福利在线观看吧| 国产成人av教育| 成人18禁在线播放| 99国产精品一区二区三区| 老汉色∧v一级毛片| 一a级毛片在线观看| 国产精品免费一区二区三区在线| 久久久久久久久久黄片| 亚洲一码二码三码区别大吗| 99精品久久久久人妻精品| 观看免费一级毛片| 国产伦一二天堂av在线观看| 亚洲熟女毛片儿| 巨乳人妻的诱惑在线观看| 天堂√8在线中文| 视频在线观看一区二区三区| 97超级碰碰碰精品色视频在线观看| 正在播放国产对白刺激| 欧美精品啪啪一区二区三区| 日日干狠狠操夜夜爽| 一a级毛片在线观看| 在线天堂中文资源库| 叶爱在线成人免费视频播放| 精品乱码久久久久久99久播| 午夜福利一区二区在线看| 午夜福利18| 琪琪午夜伦伦电影理论片6080| 一本久久中文字幕| 美女扒开内裤让男人捅视频| 国产熟女xx| 一区二区三区激情视频| 国产人伦9x9x在线观看| 久久热在线av| 天天一区二区日本电影三级| 欧美不卡视频在线免费观看 | 岛国在线观看网站| 亚洲精品在线观看二区| 深夜精品福利| 少妇裸体淫交视频免费看高清 | 精品国产乱码久久久久久男人| 欧美成人一区二区免费高清观看 | 亚洲国产精品999在线| 久久久久久久久久黄片| 男女之事视频高清在线观看| 视频在线观看一区二区三区| aaaaa片日本免费| 中文字幕久久专区| 在线国产一区二区在线| 亚洲精品国产区一区二| 美女大奶头视频| 国产精品爽爽va在线观看网站 | 亚洲国产精品久久男人天堂| 在线av久久热| 午夜激情福利司机影院| 18禁黄网站禁片午夜丰满| 最近最新中文字幕大全免费视频| 美女高潮喷水抽搐中文字幕| 亚洲精品在线观看二区| 超碰成人久久| 亚洲一区二区三区不卡视频| 一二三四在线观看免费中文在| 91字幕亚洲| 亚洲第一青青草原| 听说在线观看完整版免费高清| 午夜福利一区二区在线看| 好男人电影高清在线观看| 午夜两性在线视频| 国产熟女xx| 久久人妻av系列| 免费看日本二区| 亚洲人成伊人成综合网2020| 一本精品99久久精品77| 国产成人精品久久二区二区免费| 操出白浆在线播放| 一个人观看的视频www高清免费观看 | 激情在线观看视频在线高清| 久久青草综合色| 国产99久久九九免费精品| 不卡av一区二区三区| 午夜激情av网站| 亚洲专区中文字幕在线| 久久精品aⅴ一区二区三区四区| 国产精品一区二区三区四区久久 | 午夜精品在线福利| 国产熟女午夜一区二区三区| 免费看日本二区| 欧美丝袜亚洲另类 | 免费无遮挡裸体视频| 身体一侧抽搐| 成人三级黄色视频| 国产精品二区激情视频| 国内毛片毛片毛片毛片毛片| 久久久久九九精品影院| 两人在一起打扑克的视频| 精品久久久久久久久久久久久 | 美女国产高潮福利片在线看| 国产亚洲精品第一综合不卡| 老鸭窝网址在线观看| 中文在线观看免费www的网站 | 亚洲国产欧洲综合997久久, | 在线视频色国产色| 欧美性猛交╳xxx乱大交人| 一本久久中文字幕| 日日夜夜操网爽| 男女午夜视频在线观看| 母亲3免费完整高清在线观看| 黑人欧美特级aaaaaa片| 男人舔女人的私密视频| 午夜亚洲福利在线播放| 国产片内射在线| 日韩欧美 国产精品| 极品教师在线免费播放| 午夜久久久在线观看| 亚洲熟女毛片儿| 亚洲av日韩精品久久久久久密| 亚洲第一欧美日韩一区二区三区| 免费在线观看影片大全网站| 久久中文字幕人妻熟女| 国产成人啪精品午夜网站| 日韩中文字幕欧美一区二区| 午夜福利免费观看在线| 亚洲欧美精品综合久久99| 一级黄色大片毛片| 国产精品九九99| 国产成人欧美在线观看| 黑人巨大精品欧美一区二区mp4| 一二三四社区在线视频社区8| 97超级碰碰碰精品色视频在线观看| 免费无遮挡裸体视频| 亚洲国产欧洲综合997久久, | 男人舔女人下体高潮全视频| 一进一出抽搐动态| 丝袜在线中文字幕| 丝袜美腿诱惑在线| 久久伊人香网站| 每晚都被弄得嗷嗷叫到高潮| 日本 欧美在线| 久久青草综合色| 午夜影院日韩av| 18美女黄网站色大片免费观看| 亚洲 欧美一区二区三区| 日韩精品青青久久久久久| 男女之事视频高清在线观看| 可以在线观看毛片的网站| 国产精品影院久久| 亚洲九九香蕉| 国产单亲对白刺激| 亚洲国产日韩欧美精品在线观看 | 日韩免费av在线播放| 88av欧美| 一本久久中文字幕| 久久人妻av系列| 久久精品91蜜桃| 少妇被粗大的猛进出69影院| 欧美黑人巨大hd| 欧美另类亚洲清纯唯美| 久久精品国产综合久久久| 在线观看免费视频日本深夜| 欧美在线黄色| 国产一区二区激情短视频| 黄片大片在线免费观看| 精品日产1卡2卡| 99久久99久久久精品蜜桃| 在线观看舔阴道视频| 国产精品98久久久久久宅男小说| 久久婷婷成人综合色麻豆| 搡老岳熟女国产| 亚洲全国av大片| 韩国精品一区二区三区| 自线自在国产av| 日韩免费av在线播放| 女人高潮潮喷娇喘18禁视频| 亚洲av美国av| 给我免费播放毛片高清在线观看| 亚洲免费av在线视频| 国产亚洲欧美98| av片东京热男人的天堂| 亚洲免费av在线视频| 首页视频小说图片口味搜索| 久久久国产成人精品二区| 精品欧美一区二区三区在线| 91麻豆av在线| 精品福利观看| 亚洲五月婷婷丁香| 成人三级黄色视频| 精品国产一区二区三区四区第35| 啦啦啦免费观看视频1| 久久中文字幕一级| e午夜精品久久久久久久| 国产成人精品久久二区二区免费| 神马国产精品三级电影在线观看 | 少妇粗大呻吟视频| 1024香蕉在线观看| 一边摸一边抽搐一进一小说| 欧美日韩亚洲国产一区二区在线观看| 色播在线永久视频| 日韩免费av在线播放| 视频区欧美日本亚洲| 欧美日韩福利视频一区二区| 成人欧美大片| 久久久久久九九精品二区国产 | 久久中文字幕人妻熟女| 夜夜夜夜夜久久久久| 可以在线观看毛片的网站| 国产爱豆传媒在线观看 | 国产精品精品国产色婷婷| 在线看三级毛片| 亚洲国产日韩欧美精品在线观看 | 精品无人区乱码1区二区| 国产精品1区2区在线观看.| 亚洲狠狠婷婷综合久久图片| 亚洲欧美日韩高清在线视频| 一进一出好大好爽视频| av片东京热男人的天堂| 欧美亚洲日本最大视频资源| 亚洲黑人精品在线| 韩国av一区二区三区四区| 亚洲国产欧洲综合997久久, | 亚洲熟妇中文字幕五十中出| 久久久久精品国产欧美久久久| 色在线成人网| 麻豆国产av国片精品| 人人澡人人妻人| 欧美+亚洲+日韩+国产| 黄片大片在线免费观看| 又黄又爽又免费观看的视频| 叶爱在线成人免费视频播放| √禁漫天堂资源中文www| 精品国产超薄肉色丝袜足j| 国内毛片毛片毛片毛片毛片| 丰满人妻熟妇乱又伦精品不卡| 少妇粗大呻吟视频| a在线观看视频网站| 少妇熟女aⅴ在线视频| 国产午夜精品久久久久久| 正在播放国产对白刺激| 18禁黄网站禁片免费观看直播| 人人妻人人澡人人看| 亚洲电影在线观看av| 日日干狠狠操夜夜爽| 人人妻人人看人人澡| 日本精品一区二区三区蜜桃| 精品久久久久久久久久久久久 | 热99re8久久精品国产| 亚洲 欧美 日韩 在线 免费| 一夜夜www| 色播在线永久视频| 一二三四社区在线视频社区8| 日本在线视频免费播放| 一级a爱片免费观看的视频| 99精品久久久久人妻精品| 国产视频一区二区在线看| 亚洲av熟女| 香蕉久久夜色| 亚洲人成77777在线视频| 欧美色视频一区免费| 老司机在亚洲福利影院| 精品国产乱码久久久久久男人| 日日爽夜夜爽网站| 岛国视频午夜一区免费看| 亚洲国产精品久久男人天堂| 亚洲片人在线观看| 可以免费在线观看a视频的电影网站| 自线自在国产av| 国产av又大| 国产视频一区二区在线看| 免费观看人在逋| 看片在线看免费视频| 深夜精品福利| 黄色a级毛片大全视频| 国产高清激情床上av| 中出人妻视频一区二区| 亚洲精品一区av在线观看| 久99久视频精品免费| 亚洲熟女毛片儿| 亚洲五月天丁香| 欧美日韩瑟瑟在线播放| 制服丝袜大香蕉在线| 国产亚洲精品综合一区在线观看 | 久久性视频一级片| 变态另类成人亚洲欧美熟女| a级毛片在线看网站| 国产精品乱码一区二三区的特点| 精品免费久久久久久久清纯| 不卡av一区二区三区| 黄色视频,在线免费观看| 一级毛片高清免费大全| www.999成人在线观看| 国产又黄又爽又无遮挡在线| 宅男免费午夜| 国产精品一区二区免费欧美| 日韩欧美三级三区| 国产又黄又爽又无遮挡在线| 人人妻人人澡欧美一区二区| a级毛片在线看网站| 一级毛片女人18水好多| 韩国精品一区二区三区| 大型av网站在线播放| 欧美国产日韩亚洲一区| av天堂在线播放| 亚洲熟妇熟女久久| 一进一出抽搐动态| 99国产精品一区二区三区| 欧美在线黄色| 中亚洲国语对白在线视频| www.精华液| 久久国产精品影院| 两个人免费观看高清视频| 老鸭窝网址在线观看| 精品久久久久久久久久久久久 | 国产av一区二区精品久久| 麻豆国产av国片精品| 在线十欧美十亚洲十日本专区| 香蕉av资源在线| 国产成人av教育| 黄色成人免费大全| 精品久久久久久久毛片微露脸| 丰满人妻熟妇乱又伦精品不卡| 999精品在线视频| 免费观看人在逋| 日韩欧美国产在线观看| 亚洲真实伦在线观看| 嫁个100分男人电影在线观看| 亚洲av片天天在线观看| 国产免费av片在线观看野外av| 别揉我奶头~嗯~啊~动态视频| 中文字幕av电影在线播放| 免费在线观看影片大全网站| 国产av不卡久久| 成人三级做爰电影| 久久精品人妻少妇| 亚洲精华国产精华精| 黄色丝袜av网址大全| 侵犯人妻中文字幕一二三四区| 成人免费观看视频高清| 欧美久久黑人一区二区| av中文乱码字幕在线| 精品欧美国产一区二区三| 日日爽夜夜爽网站| 国产成人系列免费观看| av视频在线观看入口| 国产成人av激情在线播放| 波多野结衣高清无吗| 色精品久久人妻99蜜桃| 丁香六月欧美| 他把我摸到了高潮在线观看| aaaaa片日本免费| 国产伦在线观看视频一区| 亚洲精品色激情综合| 国产单亲对白刺激| 国产高清视频在线播放一区| 在线观看午夜福利视频| 日韩国内少妇激情av| 变态另类丝袜制服| 特大巨黑吊av在线直播 | 午夜福利在线观看吧| 看黄色毛片网站| 亚洲激情在线av| 国产单亲对白刺激| 成人国产一区最新在线观看| 黄片大片在线免费观看| 日韩国内少妇激情av| 90打野战视频偷拍视频| 国产成人啪精品午夜网站| 亚洲av成人不卡在线观看播放网| 这个男人来自地球电影免费观看| 国产亚洲精品av在线| 国产精品免费一区二区三区在线| 国产亚洲欧美精品永久| 日韩欧美 国产精品| 免费看十八禁软件| 亚洲专区字幕在线| 国产精品一区二区三区四区久久 | 黑人巨大精品欧美一区二区mp4| 日韩视频一区二区在线观看| 国产激情欧美一区二区| 国产成年人精品一区二区| 午夜成年电影在线免费观看| 国产av在哪里看| 久久午夜亚洲精品久久| 国产精品免费视频内射| 精品免费久久久久久久清纯| 精品国产国语对白av| 亚洲,欧美精品.| 国产亚洲欧美精品永久| 国产一区二区三区视频了| 久久久久久久久久黄片| 老司机福利观看| 午夜免费成人在线视频| tocl精华| 亚洲五月色婷婷综合| 性色av乱码一区二区三区2| 国产激情偷乱视频一区二区| 久久精品国产亚洲av香蕉五月| 精品欧美一区二区三区在线| 国产精品影院久久| 一级作爱视频免费观看| 黄网站色视频无遮挡免费观看| 免费人成视频x8x8入口观看| 欧美 亚洲 国产 日韩一| 午夜福利免费观看在线| 国产97色在线日韩免费| 麻豆久久精品国产亚洲av| 亚洲五月婷婷丁香| 亚洲全国av大片| 成人国语在线视频| 丰满的人妻完整版| 可以在线观看毛片的网站| 黄色视频,在线免费观看| 久久久久国产精品人妻aⅴ院| 18美女黄网站色大片免费观看| 亚洲精品中文字幕在线视频| 99久久精品国产亚洲精品| www日本黄色视频网| 亚洲欧洲精品一区二区精品久久久| 天堂动漫精品| 亚洲国产日韩欧美精品在线观看 | 国产亚洲精品一区二区www| 国产成人精品久久二区二区91| 在线观看免费视频日本深夜| 亚洲成人久久爱视频| 波多野结衣av一区二区av| 欧美国产日韩亚洲一区| 黄色视频不卡| 亚洲自偷自拍图片 自拍| 国产精品av久久久久免费| 久久国产亚洲av麻豆专区| 一二三四社区在线视频社区8| 欧美大码av| av在线播放免费不卡| 亚洲成a人片在线一区二区| 亚洲专区中文字幕在线| 午夜精品久久久久久毛片777| 女同久久另类99精品国产91| 久久精品成人免费网站| 久久国产精品影院| 母亲3免费完整高清在线观看| 少妇熟女aⅴ在线视频| 91av网站免费观看| 搡老熟女国产l中国老女人| 91成年电影在线观看| 一区二区三区国产精品乱码| 亚洲av日韩精品久久久久久密| 欧美最黄视频在线播放免费| 日韩欧美三级三区| 久久国产精品人妻蜜桃| 一区二区日韩欧美中文字幕| 日本在线视频免费播放| 精品国产一区二区三区四区第35| 久久精品国产99精品国产亚洲性色| 成人三级黄色视频| 久久精品国产99精品国产亚洲性色| 国产男靠女视频免费网站| 在线看三级毛片| 欧美+亚洲+日韩+国产| 99热6这里只有精品| 日本免费一区二区三区高清不卡| 欧美人与性动交α欧美精品济南到| 免费在线观看视频国产中文字幕亚洲| 色综合站精品国产| 欧美一级a爱片免费观看看 | av欧美777| av片东京热男人的天堂| 波多野结衣高清无吗| 非洲黑人性xxxx精品又粗又长| 日韩精品青青久久久久久| 三级毛片av免费| 久久精品人妻少妇| 在线观看免费午夜福利视频| 国产一级毛片七仙女欲春2 | 国产精品98久久久久久宅男小说| 国产伦人伦偷精品视频| 欧美绝顶高潮抽搐喷水| 欧美日韩乱码在线| 天天一区二区日本电影三级| 最近最新中文字幕大全免费视频| 日韩国内少妇激情av| 18禁美女被吸乳视频| 在线十欧美十亚洲十日本专区| 免费观看人在逋| 国产伦一二天堂av在线观看| 久久久久久久午夜电影| 黄片播放在线免费| 欧美黄色片欧美黄色片| 人人澡人人妻人| 美女高潮喷水抽搐中文字幕| 中文在线观看免费www的网站 | 午夜福利欧美成人| 欧美成人免费av一区二区三区| 巨乳人妻的诱惑在线观看| 成人亚洲精品av一区二区| 日本 av在线| 色av中文字幕| 美女午夜性视频免费| 成熟少妇高潮喷水视频| 亚洲片人在线观看| 日韩欧美国产在线观看| 美女高潮喷水抽搐中文字幕| 免费在线观看日本一区| 久久久久久国产a免费观看| 亚洲av成人一区二区三| 90打野战视频偷拍视频| 久久久久免费精品人妻一区二区 | 国产v大片淫在线免费观看| 一区福利在线观看| 欧美大码av| 长腿黑丝高跟| 中文字幕av电影在线播放| 女性生殖器流出的白浆| 变态另类成人亚洲欧美熟女| 国产日本99.免费观看| 国产免费av片在线观看野外av| 精品卡一卡二卡四卡免费| 日韩大码丰满熟妇| 人人妻人人澡欧美一区二区| 亚洲人成网站高清观看| 狠狠狠狠99中文字幕| 亚洲一区二区三区色噜噜| 午夜老司机福利片| 色在线成人网| 亚洲男人的天堂狠狠| 欧美激情 高清一区二区三区| 黄色a级毛片大全视频| 色综合站精品国产| 亚洲一区二区三区不卡视频| 久久久久久久精品吃奶| 侵犯人妻中文字幕一二三四区| 久久久久九九精品影院| 欧美日韩黄片免| 国产国语露脸激情在线看| 村上凉子中文字幕在线| 美女扒开内裤让男人捅视频| 亚洲激情在线av| 窝窝影院91人妻| 国产精品av久久久久免费| 激情在线观看视频在线高清| 一边摸一边抽搐一进一小说| 亚洲熟女毛片儿| 别揉我奶头~嗯~啊~动态视频| 神马国产精品三级电影在线观看 | 中文字幕精品亚洲无线码一区 | 成人亚洲精品一区在线观看| 欧美 亚洲 国产 日韩一| 精品久久久久久久末码| 国产99白浆流出| 中文字幕精品亚洲无线码一区 | 欧美在线一区亚洲| 亚洲中文字幕日韩| 日韩有码中文字幕| 亚洲欧美日韩无卡精品| 熟女少妇亚洲综合色aaa.| 久久久久久久久久黄片| 亚洲中文日韩欧美视频|