• <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表層顯微硬度研究
    精品99又大又爽又粗少妇毛片| 哪个播放器可以免费观看大片| 精品国产一区二区久久| 成人18禁高潮啪啪吃奶动态图| 国产免费视频播放在线视频| 另类精品久久| 国产精品人妻久久久久久| 99视频精品全部免费 在线| 蜜臀久久99精品久久宅男| 极品人妻少妇av视频| 国产高清国产精品国产三级| 99久国产av精品国产电影| 18禁动态无遮挡网站| 久久国产精品男人的天堂亚洲 | 男女高潮啪啪啪动态图| 午夜激情av网站| 久久人人97超碰香蕉20202| 午夜免费观看性视频| 欧美最新免费一区二区三区| 一区在线观看完整版| 精品久久国产蜜桃| 女性被躁到高潮视频| 精品亚洲乱码少妇综合久久| 国产一区二区三区综合在线观看 | 国产无遮挡羞羞视频在线观看| 亚洲三级黄色毛片| 免费高清在线观看日韩| 亚洲国产av新网站| 久久韩国三级中文字幕| 国产精品久久久久久精品古装| 久久狼人影院| 日本黄大片高清| 亚洲国产av新网站| 啦啦啦视频在线资源免费观看| a级毛色黄片| 青春草视频在线免费观看| 欧美最新免费一区二区三区| 国产一区二区在线观看日韩| 久久影院123| 久久久久网色| 欧美97在线视频| 久久女婷五月综合色啪小说| 日本wwww免费看| 久久午夜综合久久蜜桃| 欧美xxxx性猛交bbbb| 日日啪夜夜爽| 精品酒店卫生间| 久久久亚洲精品成人影院| 十分钟在线观看高清视频www| 校园人妻丝袜中文字幕| 成人二区视频| 在线天堂中文资源库| 亚洲精品,欧美精品| 搡老乐熟女国产| 少妇 在线观看| 日韩一区二区视频免费看| 永久网站在线| 国产精品免费大片| 国产片特级美女逼逼视频| 极品少妇高潮喷水抽搐| 天堂中文最新版在线下载| 亚洲av成人精品一二三区| 亚洲av电影在线观看一区二区三区| 亚洲国产色片| h视频一区二区三区| 大陆偷拍与自拍| 亚洲欧美成人综合另类久久久| 桃花免费在线播放| 久久99热这里只频精品6学生| 国产日韩欧美视频二区| 国产伦理片在线播放av一区| 两性夫妻黄色片 | 国产精品久久久av美女十八| 久久99一区二区三区| 七月丁香在线播放| 我要看黄色一级片免费的| av线在线观看网站| 免费av不卡在线播放| 人人妻人人澡人人看| 啦啦啦中文免费视频观看日本| 一级黄片播放器| 日韩精品免费视频一区二区三区 | 一二三四中文在线观看免费高清| 国产69精品久久久久777片| 另类亚洲欧美激情| 1024视频免费在线观看| 青春草亚洲视频在线观看| 黄色视频在线播放观看不卡| 日韩在线高清观看一区二区三区| 超色免费av| 国产片内射在线| 51国产日韩欧美| 日韩精品免费视频一区二区三区 | 高清不卡的av网站| 成人二区视频| 国产精品一区二区在线观看99| 日韩伦理黄色片| 啦啦啦视频在线资源免费观看| 18禁观看日本| 肉色欧美久久久久久久蜜桃| 大码成人一级视频| 一二三四在线观看免费中文在 | 国产精品国产三级国产av玫瑰| 久久午夜福利片| 久久人人97超碰香蕉20202| 黄色配什么色好看| 亚洲国产av影院在线观看| 精品国产国语对白av| 免费黄网站久久成人精品| 国产欧美亚洲国产| 桃花免费在线播放| 国产精品久久久久久精品电影小说| 欧美97在线视频| 母亲3免费完整高清在线观看 | 成人影院久久| 亚洲精品成人av观看孕妇| av在线老鸭窝| 全区人妻精品视频| 日韩大片免费观看网站| 最近2019中文字幕mv第一页| 亚洲精品一二三| 亚洲激情五月婷婷啪啪| 最近中文字幕高清免费大全6| 婷婷色av中文字幕| 久久女婷五月综合色啪小说| 成人国产av品久久久| 国产极品天堂在线| 校园人妻丝袜中文字幕| av黄色大香蕉| 久久99一区二区三区| 国产成人精品婷婷| 午夜激情av网站| 秋霞伦理黄片| 欧美日韩成人在线一区二区| 青春草亚洲视频在线观看| 黄色 视频免费看| 精品酒店卫生间| 少妇人妻 视频| 美女内射精品一级片tv| 美女中出高潮动态图| 男女啪啪激烈高潮av片| 两个人看的免费小视频| 自线自在国产av| 99国产综合亚洲精品| 成人国产av品久久久| 中文字幕最新亚洲高清| 欧美日韩亚洲高清精品| 寂寞人妻少妇视频99o| 在线观看三级黄色| 99热这里只有是精品在线观看| 在线精品无人区一区二区三| 亚洲伊人久久精品综合| 9热在线视频观看99| 国产又色又爽无遮挡免| 少妇被粗大猛烈的视频| 在线 av 中文字幕| 久久毛片免费看一区二区三区| 蜜臀久久99精品久久宅男| 在线观看美女被高潮喷水网站| 中文字幕亚洲精品专区| 国产高清不卡午夜福利| 精品久久国产蜜桃| 三上悠亚av全集在线观看| 90打野战视频偷拍视频| 爱豆传媒免费全集在线观看| 日韩一区二区三区影片| 免费av不卡在线播放| 日韩中文字幕视频在线看片| 日韩一区二区三区影片| 国产精品久久久久久久电影| 久久午夜综合久久蜜桃| 国产乱人偷精品视频| 肉色欧美久久久久久久蜜桃| 9色porny在线观看| 国产色爽女视频免费观看| 欧美精品一区二区大全| 久久久久久久精品精品| 亚洲国产毛片av蜜桃av| 国产黄色免费在线视频| av国产久精品久网站免费入址| 免费黄频网站在线观看国产| 波多野结衣一区麻豆| 国产免费视频播放在线视频| 男人添女人高潮全过程视频| 亚洲天堂av无毛| 午夜福利视频精品| 丁香六月天网| 日产精品乱码卡一卡2卡三| 色视频在线一区二区三区| 国产成人欧美| 熟女电影av网| 欧美日韩视频高清一区二区三区二| 九九爱精品视频在线观看| 成人亚洲欧美一区二区av| 精品久久蜜臀av无| 亚洲精品中文字幕在线视频| videosex国产| 亚洲国产精品一区三区| 免费少妇av软件| 桃花免费在线播放| 欧美日韩亚洲高清精品| 久久99一区二区三区| 制服人妻中文乱码| 亚洲少妇的诱惑av| 波多野结衣一区麻豆| 又黄又粗又硬又大视频| 在线精品无人区一区二区三| 国产激情久久老熟女| 美女视频免费永久观看网站| 亚洲精品aⅴ在线观看| 黑人欧美特级aaaaaa片| 中文字幕亚洲精品专区| 1024视频免费在线观看| 国产精品久久久久成人av| 国产精品人妻久久久久久| 少妇人妻 视频| 黄色视频在线播放观看不卡| 99国产综合亚洲精品| 亚洲欧洲精品一区二区精品久久久 | 欧美人与性动交α欧美精品济南到 | 免费大片18禁| 亚洲国产看品久久| 午夜av观看不卡| 免费观看性生交大片5| 亚洲国产精品成人久久小说| 午夜福利视频精品| 97精品久久久久久久久久精品| 精品一品国产午夜福利视频| 欧美精品人与动牲交sv欧美| 日韩伦理黄色片| 黑人欧美特级aaaaaa片| 51国产日韩欧美| 18+在线观看网站| 这个男人来自地球电影免费观看 | 日韩精品免费视频一区二区三区 | 亚洲欧美日韩卡通动漫| 国产成人欧美| 亚洲,欧美,日韩| 哪个播放器可以免费观看大片| 国产成人av激情在线播放| 国产无遮挡羞羞视频在线观看| 中文字幕人妻丝袜制服| 亚洲,一卡二卡三卡| 国产xxxxx性猛交| 亚洲色图 男人天堂 中文字幕 | 久久久久久久久久人人人人人人| 青春草亚洲视频在线观看| 乱码一卡2卡4卡精品| 国产极品天堂在线| 夜夜骑夜夜射夜夜干| 国产在线视频一区二区| 97在线视频观看| 国产深夜福利视频在线观看| 国产一区二区三区av在线| www日本在线高清视频| 欧美激情极品国产一区二区三区 | 国产白丝娇喘喷水9色精品| 国产亚洲午夜精品一区二区久久| 天天躁夜夜躁狠狠躁躁| 亚洲av中文av极速乱| 国产精品久久久久久av不卡| 亚洲欧美成人精品一区二区| 一本色道久久久久久精品综合| 亚洲精品第二区| 久久久久久久久久久免费av| 久久久精品94久久精品| tube8黄色片| 亚洲色图 男人天堂 中文字幕 | 中文字幕人妻丝袜制服| 丝袜美足系列| 在线精品无人区一区二区三| 激情五月婷婷亚洲| 丝袜喷水一区| 国产国语露脸激情在线看| 少妇 在线观看| 99精国产麻豆久久婷婷| av在线播放精品| 亚洲丝袜综合中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 欧美3d第一页| 欧美精品一区二区大全| 99久久综合免费| 午夜福利,免费看| 大码成人一级视频| 精品亚洲成a人片在线观看| 欧美日韩视频精品一区| 亚洲人与动物交配视频| 在线看a的网站| 亚洲欧美日韩卡通动漫| 亚洲人成网站在线观看播放| 亚洲欧洲国产日韩| 亚洲四区av| 久久久久久久久久成人| 亚洲av.av天堂| 亚洲国产最新在线播放| 免费看光身美女| 夫妻午夜视频| 青春草视频在线免费观看| www.色视频.com| 少妇熟女欧美另类| 欧美精品人与动牲交sv欧美| 2021少妇久久久久久久久久久| 久久精品国产亚洲av涩爱| 中文精品一卡2卡3卡4更新| 丝袜在线中文字幕| 街头女战士在线观看网站| 中文乱码字字幕精品一区二区三区| 女性生殖器流出的白浆| 日日爽夜夜爽网站| 久久久久精品人妻al黑| 9色porny在线观看| 国产不卡av网站在线观看| 欧美精品国产亚洲| 久久久a久久爽久久v久久| 五月伊人婷婷丁香| 一边摸一边做爽爽视频免费| 午夜日本视频在线| 人妻少妇偷人精品九色| 成人毛片60女人毛片免费| 9191精品国产免费久久| 国产欧美亚洲国产| 一区二区三区四区激情视频| 高清视频免费观看一区二区| 久久久久人妻精品一区果冻| 国产av一区二区精品久久| 99热网站在线观看| 18禁动态无遮挡网站| 亚洲精品久久久久久婷婷小说| 99久国产av精品国产电影| 久久久久久久久久人人人人人人| 国产成人a∨麻豆精品| 亚洲欧美成人精品一区二区| 亚洲精品国产色婷婷电影| 巨乳人妻的诱惑在线观看| 日日撸夜夜添| 伊人久久国产一区二区| 精品国产乱码久久久久久小说| a级毛色黄片| 成人毛片a级毛片在线播放| 国语对白做爰xxxⅹ性视频网站| 在线天堂中文资源库| 久久午夜综合久久蜜桃| 欧美老熟妇乱子伦牲交| 免费av中文字幕在线| 午夜影院在线不卡| 成人国语在线视频| 日日啪夜夜爽| www.色视频.com| 18禁动态无遮挡网站| 国产精品国产三级专区第一集| 日韩 亚洲 欧美在线| 你懂的网址亚洲精品在线观看| 精品少妇内射三级| 波野结衣二区三区在线| 免费观看无遮挡的男女| 欧美激情极品国产一区二区三区 | 国产av码专区亚洲av| 国产av国产精品国产| 免费在线观看完整版高清| 有码 亚洲区| 国产一区有黄有色的免费视频| 国产极品粉嫩免费观看在线| 久久久久久久国产电影| 国产熟女午夜一区二区三区| 婷婷成人精品国产| 亚洲熟女精品中文字幕| 晚上一个人看的免费电影| 2022亚洲国产成人精品| 亚洲精品日本国产第一区| 色网站视频免费| 亚洲婷婷狠狠爱综合网| 久久国产精品大桥未久av| 毛片一级片免费看久久久久| 色94色欧美一区二区| 日本av手机在线免费观看| 尾随美女入室| 久久综合国产亚洲精品| 国产不卡av网站在线观看| √禁漫天堂资源中文www| videos熟女内射| 久久久久久伊人网av| 九九爱精品视频在线观看| 国产一区二区在线观看av| 丝瓜视频免费看黄片| 99热国产这里只有精品6| 久久久久久久大尺度免费视频| 午夜av观看不卡| 国产精品蜜桃在线观看| 久久久久精品久久久久真实原创| 国产乱人偷精品视频| 成人国产av品久久久| 水蜜桃什么品种好| 2022亚洲国产成人精品| 国产精品久久久久久精品古装| av免费观看日本| 中文字幕制服av| 日韩不卡一区二区三区视频在线| 男女午夜视频在线观看 | 久久毛片免费看一区二区三区| av一本久久久久| 国产不卡av网站在线观看| 在线观看国产h片| 久久久欧美国产精品| 制服人妻中文乱码| 看免费成人av毛片| 大片免费播放器 马上看| 99视频精品全部免费 在线| 婷婷色综合www| 国产福利在线免费观看视频| 黄色视频在线播放观看不卡| 亚洲精品,欧美精品| 亚洲精品中文字幕在线视频| 日本wwww免费看| 国产日韩欧美亚洲二区| 少妇的逼好多水| 午夜免费观看性视频| 伦理电影免费视频| 亚洲成av片中文字幕在线观看 | 国产精品久久久久成人av| 亚洲国产精品一区二区三区在线| 丝袜脚勾引网站| 一级,二级,三级黄色视频| 日韩一区二区三区影片| 欧美人与性动交α欧美精品济南到 | 精品少妇黑人巨大在线播放| 亚洲色图 男人天堂 中文字幕 | 女性被躁到高潮视频| 伊人亚洲综合成人网| 精品少妇内射三级| 在线天堂最新版资源| 精品少妇久久久久久888优播| 国产成人精品福利久久| 搡老乐熟女国产| 乱人伦中国视频| 美女国产视频在线观看| 高清毛片免费看| 成人国语在线视频| 亚洲美女搞黄在线观看| 日韩大片免费观看网站| 久久精品国产鲁丝片午夜精品| 高清黄色对白视频在线免费看| 最近2019中文字幕mv第一页| 国产精品久久久久久久电影| 久久久久精品性色| 国产av精品麻豆| 成人无遮挡网站| videosex国产| 成人影院久久| 国产白丝娇喘喷水9色精品| 久久影院123| 国产成人精品福利久久| 亚洲av综合色区一区| 午夜福利在线观看免费完整高清在| 亚洲 欧美一区二区三区| 日韩电影二区| 日韩大片免费观看网站| 欧美xxxx性猛交bbbb| 制服诱惑二区| 性色av一级| 亚洲美女视频黄频| 飞空精品影院首页| 一本—道久久a久久精品蜜桃钙片| 晚上一个人看的免费电影| 亚洲av.av天堂| 欧美3d第一页| 亚洲av.av天堂| 婷婷成人精品国产| av播播在线观看一区| 熟妇人妻不卡中文字幕| 成人黄色视频免费在线看| 最近2019中文字幕mv第一页| 中文字幕人妻熟女乱码| 麻豆乱淫一区二区| 国产麻豆69| 欧美成人精品欧美一级黄| www.色视频.com| 国产精品蜜桃在线观看| 国产精品成人在线| 久久精品国产鲁丝片午夜精品| 欧美人与性动交α欧美精品济南到 | 精品国产一区二区久久| 美女中出高潮动态图| 日日撸夜夜添| 午夜免费男女啪啪视频观看| 日本猛色少妇xxxxx猛交久久| 国产伦理片在线播放av一区| 成人漫画全彩无遮挡| www.av在线官网国产| 久久99热这里只频精品6学生| 99精国产麻豆久久婷婷| 国产成人免费无遮挡视频| 天堂中文最新版在线下载| 欧美日韩国产mv在线观看视频| 大陆偷拍与自拍| 国产精品久久久久久av不卡| 久久精品熟女亚洲av麻豆精品| 久久人人爽av亚洲精品天堂| 亚洲成人手机| 香蕉丝袜av| 亚洲精品第二区| 久久久久久久久久久久大奶| 亚洲精品456在线播放app| 超色免费av| 久久精品国产亚洲av天美| 国产一区二区在线观看日韩| 一级爰片在线观看| 18禁裸乳无遮挡动漫免费视频| 另类亚洲欧美激情| 中文字幕最新亚洲高清| 欧美最新免费一区二区三区| 中文字幕av电影在线播放| 性高湖久久久久久久久免费观看| 国产1区2区3区精品| 日韩成人伦理影院| av国产久精品久网站免费入址| 伦理电影大哥的女人| 亚洲国产精品成人久久小说| 男女免费视频国产| 久热久热在线精品观看| 男人操女人黄网站| 免费少妇av软件| 午夜激情久久久久久久| www.av在线官网国产| 久久韩国三级中文字幕| h视频一区二区三区| www.熟女人妻精品国产 | 国产亚洲最大av| 91成人精品电影| av天堂久久9| 久久毛片免费看一区二区三区| 久久韩国三级中文字幕| 国产精品一区二区在线观看99| 亚洲欧美一区二区三区国产| 一级片免费观看大全| 飞空精品影院首页| 亚洲天堂av无毛| 亚洲第一av免费看| 精品久久蜜臀av无| 精品人妻熟女毛片av久久网站| 亚洲av.av天堂| 国产av一区二区精品久久| 国产片内射在线| 两个人免费观看高清视频| 97在线视频观看| 日本91视频免费播放| 国产男女内射视频| 精品一区二区三卡| 中文字幕av电影在线播放| 看免费成人av毛片| 水蜜桃什么品种好| 男女啪啪激烈高潮av片| 黄色 视频免费看| 又粗又硬又长又爽又黄的视频| 久久久a久久爽久久v久久| 久久综合国产亚洲精品| 欧美人与善性xxx| 99久国产av精品国产电影| 制服丝袜香蕉在线| 女人被躁到高潮嗷嗷叫费观| 亚洲欧美成人综合另类久久久| 人妻 亚洲 视频| 一级a做视频免费观看| 岛国毛片在线播放| 高清不卡的av网站| 精品人妻偷拍中文字幕| 侵犯人妻中文字幕一二三四区| 亚洲综合色惰| av黄色大香蕉| 综合色丁香网| 国产精品.久久久| 99热全是精品| 69精品国产乱码久久久| 青春草视频在线免费观看| 国产国拍精品亚洲av在线观看| 九草在线视频观看| 久久韩国三级中文字幕| 国产黄色免费在线视频| 国产乱来视频区| 色吧在线观看| 赤兔流量卡办理| 黑人欧美特级aaaaaa片| 在线看a的网站| 亚洲国产精品成人久久小说| 国产探花极品一区二区| www日本在线高清视频| 边亲边吃奶的免费视频| 国产深夜福利视频在线观看| 亚洲av男天堂| 精品少妇黑人巨大在线播放| 日本av手机在线免费观看| 男的添女的下面高潮视频| av黄色大香蕉| 久久久久视频综合| 久久国产精品男人的天堂亚洲 | 国产在视频线精品| 色网站视频免费| 欧美亚洲日本最大视频资源| 国产精品成人在线| 热99国产精品久久久久久7| 夫妻性生交免费视频一级片| 韩国av在线不卡| 欧美日韩视频精品一区| 一区二区日韩欧美中文字幕 | 色婷婷久久久亚洲欧美| 99热国产这里只有精品6| 内地一区二区视频在线| 久久久久人妻精品一区果冻| 亚洲久久久国产精品| 日韩 亚洲 欧美在线| 免费人妻精品一区二区三区视频| 在线观看一区二区三区激情| 午夜激情久久久久久久| 免费高清在线观看日韩| 久久99热这里只频精品6学生| 一级片免费观看大全|