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

    柴達(dá)木盆地灌區(qū)枸杞根系水分吸收來(lái)源研究*

    2021-03-09 02:34:38周艷清高曉東王嘉昕趙西寧
    關(guān)鍵詞:利用

    周艷清, 高曉東, 王嘉昕, 趙西寧**

    柴達(dá)木盆地灌區(qū)枸杞根系水分吸收來(lái)源研究*

    周艷清1,2, 高曉東2,3, 王嘉昕1,2, 趙西寧2,3**

    (1. 西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院 楊凌 712100; 2. 西北農(nóng)林科技大學(xué)旱區(qū)農(nóng)業(yè)水土工程教育部重點(diǎn)實(shí)驗(yàn)室 楊凌 712100; 3. 中國(guó)科學(xué)院水利部水土保持研究所 楊凌 712100)

    探討柴達(dá)木盆地灌區(qū)作物水分利用特征, 可為灌溉系統(tǒng)設(shè)計(jì)和水資源高效利用提供理論依據(jù)。以柴達(dá)木盆地懷頭他拉灌區(qū)主要作物枸杞為對(duì)象, 將枸杞潛在水源劃分為淺層(0~20 cm)、中層(20~60 cm)和深層(60~100 cm)土壤水, 利用穩(wěn)定氧同位素示蹤技術(shù)以及MixSIAR模型定量研究不同田間管理模式下[平作裸地(CK)、平作覆膜(MF)和壟作覆膜(MR)]枸杞根系水分吸收特征。結(jié)果表明: 萌芽展葉期淺層土壤水分環(huán)境相對(duì)較好, CK和MF處理枸杞主要吸收淺層土壤水分, 利用比例分別為45.9%和37.7%, MR處理對(duì)淺層、中層和深層土壤水源的利用比例相當(dāng); 開(kāi)花坐果期各層土壤水分都有所提高, 相比于CK, MF處理淺層土壤水分利用比例增加13.5%, MR處理淺層土壤水分利用比例減少11.1%; 果熟期枸杞蒸騰耗水量增加, 相比于CK, MF和MR處理對(duì)淺層土壤水分的利用比例分別增加11.7%和24.0%; 落葉期淺層土壤含水量較低, 3個(gè)處理都主要吸收深層土壤水分。淺層土壤水源的利用比例與該層土壤水分呈正相關(guān), 枸杞對(duì)淺層土壤水分變化較敏感。3種管理模式下枸杞根系水分吸收來(lái)源差異明顯。覆膜和壟作措施均提高了土壤含水量。相比于CK, MF和MR處理枸杞都增加了對(duì)淺層水源的利用。MR處理下土壤水分狀況較好且枸杞根系水分吸收利用模式更加靈活, 對(duì)柴達(dá)木盆地灌區(qū)是一種較優(yōu)的田間水分管理模式。

    枸杞; 干旱灌區(qū); 根系水分吸收; 土壤水分; 管理模式; MixSIAR模型

    青藏高原自然資源豐富, 但是土壤沙漠化日趨嚴(yán)重, 對(duì)高原生態(tài)環(huán)境和社會(huì)經(jīng)濟(jì)產(chǎn)生了重大影響[1]。枸杞()不僅具有防風(fēng)固沙、改善土壤結(jié)構(gòu)的功能, 而且經(jīng)濟(jì)效益良好[2?4]。近年來(lái), 柴達(dá)木盆地枸杞產(chǎn)業(yè)發(fā)展迅猛, 已成為中國(guó)第二大枸杞種植區(qū)。但柴達(dá)木盆地屬于氣候干旱、沒(méi)有灌溉就沒(méi)有農(nóng)業(yè)的地區(qū), 降雨稀少且蒸發(fā)強(qiáng)烈, 水資源短缺成為制約作物生長(zhǎng)的主要瓶頸[5-7]。隨著枸杞種植規(guī)模擴(kuò)大, 該區(qū)水資源供需矛盾將進(jìn)一步加劇。因此探明枸杞生育期內(nèi)水分利用特征, 對(duì)灌區(qū)灌溉系統(tǒng)設(shè)計(jì)和水資源高效利用具有重要科學(xué)意義和應(yīng)用價(jià)值。

    穩(wěn)定氫氧同位素技術(shù)被廣泛應(yīng)用到植物水分來(lái)源研究[8-9]。一般植物可利用的水源(降水、徑流、不同深度的土壤水、地下水等)都具有不同的穩(wěn)定氫氧同位素特征[10-11], 因此對(duì)比植物木質(zhì)部水與各潛在水源的穩(wěn)定氫氧同位素組成可以判斷植物水分來(lái)源[12], 并且可采用多源混合模型(如IsoSource模型、貝葉斯混合模型)進(jìn)一步確定植物對(duì)各潛在水源的利用范圍或利用比例[8,13]。目前, 在干旱半干旱地區(qū)利用穩(wěn)定氫氧同位素技術(shù)研究植物土壤水分來(lái)源的報(bào)道主要集中在生態(tài)林、天然和人工灌木等[14-16], 但是對(duì)于生態(tài)經(jīng)濟(jì)型枸杞的研究較為薄弱。此外, 在計(jì)算方法上, IsoSource模型是目前估算多種水分來(lái)源對(duì)植物生長(zhǎng)貢獻(xiàn)的一種重要手段[8-9,17-18], 但該模型未考慮水源的不確定性。針對(duì)此問(wèn)題有學(xué)者提出了貝葉斯混合模型[17-20], 目前, 貝葉斯混合模型主要有3種, 分別為MixSIR、SIAR和MixSIAR, 其中MixSIAR模型融合了MixSIR和SIAR模型的特點(diǎn), 增加了固定和隨機(jī)效應(yīng)(作為協(xié)變量解釋混合比例變化)、過(guò)程誤差+殘差、整合不確定性源和明確比較競(jìng)爭(zhēng)模型或參數(shù)值的支持強(qiáng)度等功能, 使模型在分析植物水分來(lái)源及其貢獻(xiàn)比例上具有更高的準(zhǔn)確性[21]。

    雖然穩(wěn)定氫氧同位素技術(shù)在植物水分利用來(lái)源識(shí)別中取得巨大成功, 但是穩(wěn)定氫同位素容易在鹽堿和干旱環(huán)境下發(fā)生分餾[22-23], 相反, 目前鮮有穩(wěn)定氧同位素發(fā)生分餾的報(bào)道。因此, 本研究采用穩(wěn)定氧同位素技術(shù), 通過(guò)采樣和室內(nèi)測(cè)定分析獲得不同田間管理模式下枸杞木質(zhì)部水及其潛在水源的18O值, 使用貝葉斯混合模型(MixSIAR)揭示不同田間管理模式下枸杞根系水分吸收生育期內(nèi)變化特征, 以期為柴達(dá)木盆地灌區(qū)枸杞灌溉系統(tǒng)設(shè)計(jì)和水分高效利用提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 試驗(yàn)區(qū)概況

    試驗(yàn)區(qū)位于我國(guó)青藏高原東北部, 青海省海西蒙古族藏族自治州德令哈市懷頭他拉灌區(qū)(37°21′14.74″N, 96°41′41.06″E), 海拔2868 m。屬于高原大陸性氣候, 日照時(shí)長(zhǎng)3007 h, 年均氣溫6.5 ℃,多年平均降水量241.3 mm, 年內(nèi)分配不均, 72%的降水集中在6—9月。2018年6—9月降水量分別為50.2 mm、26.4 mm、49.2 mm和7.8 mm。地下水埋深大于50 m。試驗(yàn)區(qū)土壤為礫質(zhì)土, 入滲能力強(qiáng)。試驗(yàn)區(qū)土壤容重和養(yǎng)分含量如表1所示, 全氮、速效磷和速效鉀含量在0~20 cm土層最高, 有機(jī)質(zhì)含量在40~60 cm土層相對(duì)較高。植被組成簡(jiǎn)單, 以具有高抗旱能力的天然荒漠灌叢[駝絨藜(), 白刺()]和人工灌木林(枸杞)為主。

    表1 試驗(yàn)區(qū)0~100 cm土層土壤容重及養(yǎng)分含量

    1.2 試驗(yàn)設(shè)計(jì)

    以‘寧杞7號(hào)’枸杞作為供試樹(shù)種, 樹(shù)齡5年。地表滴灌為灌區(qū)主要的灌溉方式, 滴頭間距0.5 m, 滴頭流量3.75 L?h?1, 次灌水持續(xù)時(shí)間為4 h。2018年生育期中共灌5次(灌水時(shí)間分別為6月3日、7月4日、7月25日、8月5日和8月24日), 在生育期結(jié)束后進(jìn)行冬灌(2018年10月20日)。

    試驗(yàn)設(shè)置3個(gè)處理: 平作裸地(CK)、平作覆膜(MF)和壟作覆膜(MR)。處理間有1.5 m的保護(hù)行。其中平作裸地為灌區(qū)傳統(tǒng)的田間管理方式; 壟作覆膜處理中壟寬20 cm, 高30 cm, 溝寬20 cm, 枸杞種植在溝中, 所用地膜為黑色園藝地布, 經(jīng)久耐用, 可連續(xù)使用, 于2018年4月上旬整地覆蓋。試驗(yàn)采用隨機(jī)區(qū)組設(shè)計(jì), 每個(gè)處理3個(gè)重復(fù), 共9個(gè)小區(qū), 每個(gè)小區(qū)面積30 m2(1.5 m′20 m)。所有處理除壟作和覆蓋管理措施外, 其余耕作、果樹(shù)修剪、施肥、灌溉和病蟲(chóng)害防治處理均一致。試驗(yàn)枸杞基本信息如表2所示。

    表2 不同處理的枸杞基本信息

    每個(gè)處理選取9株測(cè)量樹(shù)高和地徑。Height and ground diameter are mean values of 9 trees in each treatment.

    1.3 樣品采集

    2018年5月1日至9月30日(植物生長(zhǎng)季), 監(jiān)測(cè)每次降雨歷時(shí)和降雨量, 并收集降雨樣品。分別在枸杞萌芽展葉期(5月25日)、開(kāi)花坐果期(7月9日)、果熟期(7月30日)和落葉期(9月22日)取樣, 每個(gè)小區(qū)分別選取1棵長(zhǎng)勢(shì)良好, 具有代表性的枸杞樹(shù)。

    采樣方法: 選取枝葉較多且已栓化的枝條, 除去枝條的外皮和韌皮部并剪去枝條兩端, 將木質(zhì)部剪成1~2 cm, 迅速裝入瓶?jī)?nèi), 裝滿后用封口膜密封放入冷凍室保存(?20 ℃), 每株植物各采集3瓶作為重復(fù)??紤]枸杞根系分布及冠幅的影響, 選擇在枸杞植株附近30 cm處用直徑40 mm土鉆取樣, 取樣深度為100 cm。0~20 cm土層每10 cm取樣1次, 20~100 cm土層每20 cm取樣1次, 每層取3個(gè)重復(fù), 將取得的每層土樣混合均勻。各土層樣品分為兩部分: 一部分迅速裝入樣品瓶中, 用封口膜密封保存, 取樣完畢后放入冷凍室保存(?20 ℃); 另一部分放入鋁盒, 帶回實(shí)驗(yàn)室用烘干稱質(zhì)量法測(cè)量土壤質(zhì)量含水量。試驗(yàn)區(qū)面積較小, 忽略降水、溫度等氣象要素在試驗(yàn)區(qū)內(nèi)的空間異質(zhì)性, 距試驗(yàn)區(qū)100 m左右的空曠處布設(shè)便攜式自動(dòng)氣象站(Automatic Weather Station AR5, Avolon Scientific, Jersey City, USA), 連續(xù)監(jiān)測(cè)降水量(mm)、大氣相對(duì)濕度(%)、空氣溫度(℃)和太陽(yáng)輻射(W?m?2)等氣象要素。

    1.4 水分提取與分析

    采用對(duì)樣品同位素組成影響不顯著的低溫真空蒸餾法(LI-2000土壤-植物水真空抽取系統(tǒng))提取土壤水和枸杞木質(zhì)部水[24-25]。在真空狀態(tài)下(或接近真空), 將萃取瓶中的樣品加熱到100 ℃, 水蒸汽經(jīng)過(guò)U型管進(jìn)入浸在液氮中的冷凝管中。West等[25]研究表明在真空蒸餾法中要獲取未分餾的水樣, 植物木質(zhì)部樣品需要加熱60~75 min、砂土需要30 min。因此為了獲取未分餾的水樣, 本研究中植物木質(zhì)部樣抽提時(shí)間大于150 min、土樣抽提時(shí)間大于90 min。提取的水樣用Parafilm膜密封在5 mL的玻璃瓶中, 放在冷藏箱(?2 ℃)中保存直至測(cè)樣, 時(shí)間均控制在1個(gè)月內(nèi)。測(cè)樣時(shí)需要將水樣經(jīng)過(guò)濾塞轉(zhuǎn)移到1 mL標(biāo)準(zhǔn)瓶中, 最后上機(jī)測(cè)定穩(wěn)定氧同位素值。

    土壤樣品中提取的水樣用液態(tài)水同位素分析儀(TIWA-45EP, Los Gatos Research, Mountain View, USA)測(cè)量, 由于低溫真空蒸餾法提取的木質(zhì)部水樣中可能存在有機(jī)物污染, 從而干擾液態(tài)水同位素分析儀的光譜信號(hào), 導(dǎo)致同位素測(cè)量的誤差超過(guò)儀器精度。因此植物木質(zhì)部水樣采用不受有機(jī)物污染影響的同位素比率質(zhì)譜儀(Finnigan MAT Delta V advantage)進(jìn)行測(cè)量[26]。樣品穩(wěn)定氧同位素值計(jì)算公式為:

    液態(tài)水同位素分析儀18O的測(cè)量精度為±0.2‰, 同位素比率質(zhì)譜儀18O的測(cè)量精度為±0.2‰。水分提取和樣品分析均在西北農(nóng)林科技大學(xué)中國(guó)旱區(qū)節(jié)水農(nóng)業(yè)研究院水同位素分析實(shí)驗(yàn)室完成。

    1.5 數(shù)據(jù)處理

    本研究采用貝葉斯混合模型(MixSIAR)量化不同土層的吸水比例。MixSIAR模型中木質(zhì)部水和各土層土壤水同位素值輸入全部使用原始數(shù)據(jù)(例如0~20 cm土層數(shù)據(jù)輸入為0~10 cm和10~20 cm的同位素?cái)?shù)據(jù)), TDF數(shù)據(jù)設(shè)為0, 默認(rèn)不發(fā)生同位素分餾, 馬爾可夫鏈蒙特卡羅(MCMC)運(yùn)行時(shí)長(zhǎng)設(shè)定為normal, 誤差結(jié)構(gòu)設(shè)定為Residual Only, 模型結(jié)果使用Gelman和Geweke診斷確定模型是否收斂, 模型不收斂則增加MCMC運(yùn)算時(shí)長(zhǎng)(long, very long), 模型輸出結(jié)果用平均值表示。不同土層的土壤水被認(rèn)為是枸杞的潛在水源, 將土壤含水量和土壤水的同位素組成作為潛在水源的劃分依據(jù), 淺層(0~20 cm)受降雨、灌溉水和蒸發(fā)影響較大, 土壤含水量和土壤水同位素組成變化較大; 中層(20~60 cm)土壤含水量和土壤水同位素組成變化相對(duì)較小; 深層(60~100 cm)土壤含水量變化最小, 土壤水的同位素組成較穩(wěn)定。將3個(gè)潛在水源的18O值代入貝葉斯混合模型(MixSIAR)確定各水源的貢獻(xiàn)比例。

    采用Microsoft Excel 2016軟件和SPSS 22.0軟件對(duì)數(shù)據(jù)進(jìn)行整理分析, 用Origin 2016軟件作圖。

    2 結(jié)果與分析

    2.1 降水量及降水、灌溉水同位素組成變化

    試驗(yàn)期間(2018年5月1日至9月30日), 試驗(yàn)區(qū)降水量、降水和灌溉水18O值與每日平均氣溫隨時(shí)間的變化如圖1所示。試驗(yàn)期間共有52次降水事件, 總降水量為146.8 mm。且76%以上的降水事件降水量小于5 mm, 對(duì)作物生長(zhǎng)而言屬于無(wú)效降雨。降水集中在6—8月, 降水總量占采樣期間降水總量的84.4%; 6月7日降水量為最大單次降水量, 達(dá)14.8 mm。降水中18O值變化幅度較大, 為?14.4‰~1.7‰, 加權(quán)平均值為?5.6‰。降水的18O值與降水量、降水日平均氣溫的相關(guān)分析表明, 無(wú)降水量效應(yīng), 有微弱的溫度效應(yīng):18O=0.3346+ 21.327 (為降水日平均氣溫,2=0.252,=0.07)。這可能是受大陸性氣團(tuán)的影響, 水汽主要來(lái)源于局地水汽再循環(huán)[27]。灌溉水中的18O值在試驗(yàn)研究期內(nèi)變幅較小, 試驗(yàn)區(qū)灌溉水取自深層地下水, 故灌溉水18O值較穩(wěn)定。

    2.2 不同管理模式土壤含水量與土壤水δ18O特征

    試驗(yàn)區(qū)不同管理模式枸杞樣地的土壤含水量隨土層深度和時(shí)間的變化如圖2所示。枸杞樣地土壤含水量主要受降水、灌溉水和田間管理模式等因素的共同影響。與CK相比, 整個(gè)生長(zhǎng)季MF和MR處理0~100 cm土層平均土壤含水量分別提高了17.7%和21.9%, 但3個(gè)處理的土壤含水量隨時(shí)間變化趨勢(shì)一致。各土層的土壤含水量情況大致可分為3個(gè)區(qū)間: 淺層(0~20 cm)土壤含水量受蒸發(fā)和降水、灌溉水補(bǔ)給的影響較大, 土壤含水量有較明顯的波動(dòng)。萌芽展葉期(5月25日)和落葉期(9月22日)水源補(bǔ)給較少, 3個(gè)處理的土壤含水量均維持在較低水平, 而開(kāi)花坐果期(7月9日)和果熟期(7月30日)在取樣前5 d都有灌溉水補(bǔ)給, 淺層土壤含水量相對(duì)較高。中層(20~ 60 cm)受蒸發(fā)的影響較小, 土壤含水量隨時(shí)間變化與淺層相比波動(dòng)較小, 6月份較豐富的降雨與采樣前5 d的灌溉水補(bǔ)給了開(kāi)花坐果期的土壤水分, 但CK中層仍沒(méi)有劇烈變化。3個(gè)處理在其他生育階段土壤含水量變化都較小。深層(60~100 cm)土壤含水量變化與淺層土壤含水量變化趨勢(shì)相似, 萌芽展葉期3個(gè)處理深層土壤含水量都處于較低水平, 開(kāi)花坐果期深層土壤含水量達(dá)最大值, 3個(gè)處理在果熟期和落葉期深層土壤含水量逐漸減小。

    CK: 平作裸地; MF: 平作覆膜; MR: 壟作覆膜。CK: conventional flat planting; MF: flat planting with full film mulching; MR: ridge-furrow planting with full film mulching.

    3個(gè)處理土壤水的18O值變化規(guī)律一致, 因此將3個(gè)處理土壤水18O值一同分析, 土壤水18O值時(shí)空變化規(guī)律如圖3所示。受降水、灌溉水和蒸發(fā)的影響, 淺層土壤水18O值呈現(xiàn)較為明顯的變化, 中層土壤水18O值波動(dòng)的幅度相對(duì)表層較小; 深層土壤水18O值各月間無(wú)顯著變化。土壤含水量和土壤水18O值顯著的季節(jié)性波動(dòng)特征意味著枸杞水分利用來(lái)源亦存在季節(jié)性變化。

    2.3 不同管理模式枸杞根系水分吸收來(lái)源的定量分析

    基于18O值的MixSIAR模型水源分析結(jié)果如圖4所示。CK處理在萌芽展葉期主要吸收淺層土壤水分, 利用比例達(dá)45.9%; 開(kāi)花坐果期主要吸收深層土壤水分, 利用比例為38.2%, 對(duì)淺層和中層的利用比例也接近30%; 果熟期主要吸收中層和深層的土壤水分, 利用比例分別為37.7%和37.3%; 落葉期又增加了對(duì)深層土壤水分的利用, 利用比例達(dá)43.4%。MF處理在萌芽展葉期、開(kāi)花坐果期和果熟期主要吸收淺層土壤水分, 利用比例分別為37.7%、46.1%和36.8%; 落葉期吸收的水分主要來(lái)源為深層, 占38.7%。MR處理枸杞在整個(gè)生育期的主要水分來(lái)源不斷轉(zhuǎn)換, 具有較強(qiáng)的靈活性, 萌芽展葉期對(duì)3個(gè)土壤水源的利用比例均大于30%, 且無(wú)顯著差異; 開(kāi)花坐果期主要吸收中層和深層土壤水分, 利用比例分別為38.4%和40.1%; 果熟期轉(zhuǎn)向利用淺層水源, 利用比例達(dá)49.1%; 落葉期主要利用中層和深層水源, 利用比例分別為35.9%和37.0%。為評(píng)估枸杞對(duì)各土壤水源利用比例與土壤水分之間的關(guān)系, 將每個(gè)采樣點(diǎn)下各層土壤水源利用比例與相應(yīng)土層的土壤含水量做相關(guān)性分析??傮w而言, 這些變量之間無(wú)顯著相關(guān)性, 但淺層水源利用比例與土壤含水量呈正相關(guān)關(guān)系(圖5)。

    CK: 平作裸地; MF: 平作覆膜; MR: 壟作覆膜。CK: conventional flat planting; MF: flat planting with full film mulching; MR: ridge-furrow planting with full film mulching.

    3 討論

    土壤含水量受灌溉水、降水補(bǔ)給與蒸發(fā)蒸騰作用的影響, 使土壤含水量在生育期內(nèi)不斷變化[9,28]。受蒸發(fā)過(guò)程影響, 土壤水18O值表現(xiàn)為隨土層深度的增加迅速下降并逐漸趨于穩(wěn)定的變化趨勢(shì), 這與前人研究結(jié)果相一致[29-31]。土壤含水量與土壤水18O值在不同生育期差異明顯, 根據(jù)穩(wěn)定同位素質(zhì)量守恒定律, 這意味著枸杞水分利用來(lái)源亦將呈現(xiàn)明顯季節(jié)性特征。

    萌芽展葉期(5月25日)氣溫低, 蒸發(fā)較弱, 淺層和中層土壤中貯存著冬灌和冬季融雪補(bǔ)充的水分, 使得土壤含水量較高[31], 又因枸杞剛剛度過(guò)休眠期, 葉芽、花芽還在發(fā)育, 僅長(zhǎng)出一些葉子, 枸杞進(jìn)行光合作用的能力較弱, 且枸杞利用淺層土壤水時(shí)消耗能量相對(duì)較少[32], 平作裸地和平作覆膜處理枸杞根系都優(yōu)先吸收淺層土壤水, 但壟作覆膜處理對(duì)3個(gè)土層水源利用比例差別不大, 可能是壟作和覆膜措施提高了土壤溫濕度[33], 從而提高了各土層根系活力[34]。開(kāi)花坐果期(7月9日), 枸杞根系生長(zhǎng)最快[35]且根系活性較強(qiáng)[36], 枸杞可以靈活地運(yùn)用各層土壤水源。平作裸地處理主要利用深層土壤水分, 其主要利用水分來(lái)源與土壤含水量的高值區(qū)吻合, 并且也可能是平作裸地處理枸杞適應(yīng)了長(zhǎng)期干旱少水的環(huán)境, 根系對(duì)灌溉水反應(yīng)敏感性降低, 只有降水或灌溉水增加到一定閾值時(shí), 淺層根系才有響應(yīng), 根系才能形成保持吸收淺層土壤水分的能力[31,37-38]。平作覆膜和壟作覆膜處理各土層土壤含水量都相對(duì)充足, 平作覆膜處理仍主要利用淺層土壤水分, 壟作覆膜處理則主要利用中層和深層的土壤水分, 田間土壤含水量的長(zhǎng)期變化特征可能影響著枸杞的水分利用來(lái)源, 一方面是由于在中層和深層, 壟作覆膜處理枸杞的根系生物量都大于平作覆膜處理的根系生物量[39], 在水分環(huán)境相對(duì)較好的情況下, 壟作覆膜處理枸杞會(huì)優(yōu)先選擇相對(duì)穩(wěn)定的中層和深層土壤水源[31]。果熟期(7月30日), 枸杞地上生物量達(dá)到最大值, 需要吸收更多的養(yǎng)分。取樣5 d前的灌溉水并沒(méi)有顯著增加土壤含水量, 是由于該時(shí)期枸杞蒸散量最大, 枸杞對(duì)土壤水分消耗速率加快[39]。平作裸地處理主要利用中層和深層土壤水分, 雖然表層養(yǎng)分含量最高, 由于淺層土壤水分流失較快, 枸杞仍需要從更深土層吸收土壤水分, 平作覆膜和壟作覆膜處理主要利用淺層土壤水分, 更有利于對(duì)養(yǎng)分的吸收[15]。落葉期(9月22日), 進(jìn)入秋季, 氣溫驟降, 降雨稀少, 枸杞蒸騰耗水量減少, 各處理都增大了對(duì)深層土壤水的利用比例, 可能是淺層土壤溫度和水分降低, 淺層土壤中根系活性也逐漸降低[39-40], 枸杞只能選擇相對(duì)充足和穩(wěn)定的中層和深層土壤水源。不同田間管理模式下枸杞生育期內(nèi)主要吸水土層不同, 證明其根系吸水模式具有較強(qiáng)的可塑性, 枸杞對(duì)土壤水分狀況的適應(yīng)能力較強(qiáng)。植物對(duì)各水源的利用比例與土壤水分有效性、功能根分布密切相關(guān)[41-42]。本研究表明枸杞根系對(duì)淺層水源較為敏感, 并能對(duì)淺層水源的變化做出快速反應(yīng), 這可能與水通道蛋白的活性和豐度有關(guān)[43]。也有研究表明土壤含水量與根系吸水量、木質(zhì)部水同位素值有很好的相關(guān)性[26]。當(dāng)然土壤水分的有效性也能誘導(dǎo)根系朝向土壤水分環(huán)境良好的區(qū)域生長(zhǎng)[44]。

    干旱地區(qū)水資源嚴(yán)重缺乏, 特別是在氣候變化、生態(tài)需水與生活需水增加的背景下, 水資源供需矛盾愈加突出, 灌溉水的高效利用愈加迫切[45], 干旱地區(qū)灌區(qū)的種植模式和管理方式亟需改變升級(jí)。與傳統(tǒng)的平作裸地模式相比, 壟作覆膜是一種簡(jiǎn)單有效提高土壤水分和溫度的方法[34,39,46], 在干旱半干旱地區(qū)具有很大的節(jié)水潛力, 與平作裸地處理相比, 整個(gè)生長(zhǎng)季中壟作覆膜處理顯著增加了土壤含水量, 提高了淺層和中層水分的利用比例, 有利于養(yǎng)分吸收[15]。因此推薦壟作覆膜種植模式作為柴達(dá)木灌區(qū)枸杞種植的田間管理模式。同時(shí), 為了降低土壤無(wú)效蒸發(fā), 灌區(qū)正在進(jìn)行灌溉系統(tǒng)升級(jí), 將目前的地表滴灌改為地下灌溉, 但如何確定地下灌溉深度尚無(wú)定論。本研究通過(guò)穩(wěn)定氧同位素技術(shù)確定了不同田間管理模式枸杞的主要水分利用土層, 據(jù)此可以確定地下灌水器的埋深, 為地下灌溉系統(tǒng)設(shè)計(jì)提供科學(xué)依據(jù)。

    4 結(jié)論

    本研究基于穩(wěn)定氧同位素技術(shù), 結(jié)合MixSIAR模型分析柴達(dá)木盆地灌區(qū)不同田間管理模式下枸杞根系水分吸收來(lái)源, 揭示其在不同土壤水分環(huán)境下的水分利用策略, 3種田間管理模式下枸杞生育期內(nèi)水分利用來(lái)源不同, 平作裸地處理根系水分吸收深層土壤水的平均利用比例為36.6%, 平作覆膜處理根系水分吸收淺層土壤水的平均利用比例為37.7%, 壟作覆膜處理對(duì)3層土壤水源的利用比例相當(dāng)。枸杞根系對(duì)淺層土壤含水量的變化較為敏感(<0.05), 根系吸水模式可塑性較強(qiáng), 淺層土壤含水量較低時(shí), 枸杞根系可以通過(guò)吸收更多的深層土壤水分來(lái)滿足植株蒸騰耗水的需求。壟作和覆膜措施提高了枸杞田間淺層土壤含水量, 同時(shí)枸杞也增加了對(duì)淺層土壤水源的利用比例。整個(gè)生育期中壟作覆膜處理土壤水分狀況較好, 枸杞水分利用模式更靈活。因此在枸杞田間水分管理時(shí), 壟作覆膜模式是柴達(dá)木盆地灌區(qū)較為適宜的種植模式。

    [1] ZHANG C L, LI Q, SHEN Y P, et al. Monitoring of aeolian desertification on the Qinghai-Tibet Plateau from the 1970s to 2015 using Landsat images[J]. Science of the Total Environment, 2018, 619/620: 1648–1659

    [2] 張?bào)w彬, 康躍虎, 萬(wàn)書(shū)勤, 等. 滴灌枸杞對(duì)龜裂堿土幾種酶活性的改良效應(yīng)[J]. 土壤學(xué)報(bào), 2015, 52(6): 1392–1400 ZHANG T B, KANG Y H, WAN S Q, et al. Ameliorative effect of croppingL. with drip irrigation on soil enzymes activities in takyric solonetz[J]. Acta Pedologica Sinica, 2015, 52(6): 1392–1400

    [3] 陳艷瑞, 尹林克. 人工防風(fēng)固沙林演替中群落組成和優(yōu)勢(shì)種群生態(tài)位變化特征[J]. 植物生態(tài)學(xué)報(bào), 2008, 32(5): 1126–1133 CHEN Y R, YIN L K. Community composition and niche change characteristics of dominant species in the wind-breaking and sand-fixing forest, Xinjiang, China[J]. Journal of Plant Ecology (Chinese Version), 2008, 32(5): 1126–1133

    [4] 胥生榮, 張恩和, 馬瑞麗, 等. 干旱脅迫及復(fù)水對(duì)耐旱枸杞水力學(xué)特性的影響[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 2017, 25(8): 1190–1197 XU S R, ZHANG E H, MA R L, et al. Hydraulic characteristics ofL. seedlings under drought stress and re-watering conditions[J]. Chinese Journal of Eco-Agriculture, 2017, 25(8): 1190–1197

    [5] GAO X D, WU P T, ZHAO X N, et al. Estimation of spatial soil moisture averages in a large gully of the Loess Plateau of China through statistical and modeling solutions[J]. Journal of Hydrology, 2013, 486: 466–478

    [6] LIANG H B, XUE Y Y, LI Z S, et al. Soil moisture decline following the plantation offorests: Evidence from the Loess Plateau[J]. Forest Ecology and Management, 2018, 412: 62–69

    [7] 楊榮金, 舒儉民, 李秀紅, 等. 柴達(dá)木盆地生態(tài)環(huán)境保護(hù)戰(zhàn)略與對(duì)策[J]. 科技導(dǎo)報(bào), 2017, 35(6): 115–119 YANG R J, SHU J M, LI X H, et al. Strategy and countermeasures of ecological environment protection in Qaidam Basin[J]. Science & Technology Review, 2017, 35(6): 115–119

    [8] HUO G P, ZHAO X N, GAO X D, et al. Seasonal water use patterns of rainfed jujube trees in stands of different ages under semiarid Plantations in China[J]. Agriculture, Ecosystems & Environment, 2018, 265: 392–401

    [9] 丁丹, 賈文雄, 馬興剛, 等. 祁連山亞高山灌叢優(yōu)勢(shì)植物水分來(lái)源[J]. 生態(tài)學(xué)報(bào), 2018, 38(4): 1348–1356 DING D, JIA W X, MA X G, et al. Water source of dominant plants of the subalpine shrubland in the Qilian Mountains, China[J]. Acta Ecologica Sinica, 2018, 38(4): 1348–1356

    [10] 鄧文平, 章潔, 張志堅(jiān), 等. 北京土石山區(qū)水分在土壤-植物-大氣連續(xù)體(SPAC)中的穩(wěn)定同位素特征[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2017, 28(7): 2171–2178 DENG W P, ZHANG J, ZHANG Z J, et al. Stable hydrogen and oxygen isotope compositions in soil-plant-atmosphere continuum (SPAC) in rocky mountain area of Beijing, China[J]. Chinese Journal of Applied Ecology, 2017, 28(7): 2171–2178

    [11] 石輝, 劉世榮, 趙曉廣. 穩(wěn)定性氫氧同位素在水分循環(huán)中的應(yīng)用[J]. 水土保持學(xué)報(bào), 2003, 17(2): 163–166 SHI H, LIU S R, ZHAO X G. Application of stable hydrogen and oxygen isotope in water circulation[J]. Journal of Soil Water Conservation, 2003, 17(2): 163–166

    [12] DAWSON T E, MAMBELLI S, PLAMBOECK A H, et al. Stable isotopes in plant ecology[J]. Annual review of Ecology and Systematics, 2002, 33(1): 507–559

    [13] 朱林, 許興, 毛桂蓮. 寧夏平原北部地下水埋深淺地區(qū)不同灌木的水分來(lái)源[J]. 植物生態(tài)學(xué)報(bào), 2012, 36(7): 618–628 ZHU L, XU X, MAO G L. Water sources of shrubs grown in the northern Ningxia Plain of China characterized by shallow groundwater table[J]. Chinese Journal of Plant Ecology, 2012, 36(7): 618–628

    [14] ZHAO L J, WANG L X, CERNUSAK L A, et al. Significant difference in hydrogen isotope composition between xylem and tissue water in[J]. Plant, Cell & Environment, 2016, 39(8): 1848–1857

    [15] 呂婷, 趙西寧, 高曉東, 等. 黃土丘陵區(qū)典型天然灌叢和人工灌叢優(yōu)勢(shì)植物土壤水分利用策略[J]. 植物生態(tài)學(xué)報(bào), 2017, 41(2): 175–185 LYU T, ZHAO X N, GAO X D, et al. Soil water use strategy of dominant species in typical natural and planted shrubs in loess hilly region[J]. Chinese Journal of Plant Ecology, 2017, 41(2): 175–185

    [16] 劉保清, 劉志民, 錢建強(qiáng), 等. 科爾沁沙地南緣主要固沙植物旱季水分來(lái)源[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2017, 28(7): 2093–2101 LIU B Q, LIU Z M, QIAN J Q, et al. Water sources of dominant sand-binding plants in dry season in southern Horqin Sandy Land, China[J]. Chinese Journal of Applied Ecology, 2017, 28(7): 2093–2101

    [17] 朱雅娟, 賈志清. 秋季巴丹吉林沙漠東南緣人工梭梭林水分來(lái)源[J]. 林業(yè)科學(xué), 2012, 48(8): 1–5 ZHU Y J, JIA Z Q. Water source ofplantations in autumn at the southeast edge of Badain Jaran Desert[J]. Scientia Silvae Sinicae, 2012, 48(8): 1–5

    [18] DAI Y, ZHENG X J, TANG L S, et al. Stable oxygen isotopes reveal distinct water use patterns of twospecies in the Gurbantonggut Desert[J]. Plant and Soil, 2015, 389(1/2): 73–87

    [19] MOORE J W, SEMMENS B X. Incorporating uncertainty and prior information into stable isotope mixing models[J]. Ecology Letters, 2008, 11(5): 470–480

    [20] PARNELL A C, INGER R, BEARHOP S, et al. Source partitioning using stable isotopes: Coping with too much variation[J]. PLoS One, 2010, 5(3): e9672

    [21] STOCK B C, JACKSON A L, WARD E J, et al. Analyzing mixing systems using a new generation of Bayesian tracer mixing models[J]. PeerJ, 2018, 6: e5096

    [22] ELLSWORTH P Z, WILLIAMS D G. Hydrogen isotope fractionation during water uptake by woody xerophytes[J]. Plant and Soil, 2007, 291(1/2): 93–107

    [23] BARBETA A, JONES S P, CLAVé L U, et al. Unexplained hydrogen isotope offsets complicate the identification and quantification of tree water sources in a riparian forest[J]. Hydrology and Earth System Sciences, 2019, 23(4): 2129–2146

    [24] 馬小軍, 靳靜靜, 司炳成, 等. 提取方法對(duì)土壤水同位素和植物水源分割的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2019, 30(6): 1840–1846 MA X J, JIN J J, SI B C, et al. Effects of extraction methods on soil water isotope and plant water source segmentation[J]. Chinese Journal of Applied Ecology, 2019, 30(6): 1840–1846

    [25] WEST A G, PATRICKSON S J, EHLERINGER J R. Water extraction times for plant and soil materials used in stable isotope analysis[J]. Rapid Communications in Mass Spectrometry, 2006, 20(8): 1317–1321

    [26] GAO X D, ZHAO X N, LI H C, et al. Exotic shrub species () is more resistant to extreme natural drought than native species () in a semiarid revegetated ecosystem[J]. Agricultural and Forest Meteorology, 2018, 263: 207–216

    [27] 吳華武, 李小雁, 趙國(guó)琴, 等. 青海湖流域降水和河水中18O和D變化特征[J]. 自然資源學(xué)報(bào), 2014, 29(9): 1552–1564 WU H W, LI X Y, ZHAO G Q, et al. The variation characteristics of δ18O and δD in precipitation and river water, Qinghai Lake basin[J]. Journal of Natural Resources, 2014, 29(9): 1552–1564

    [28] 霍高鵬, 趙西寧, 高曉東, 等. 黃土丘陵區(qū)棗農(nóng)復(fù)合系統(tǒng)土壤水分利用與競(jìng)爭(zhēng)[J]. 自然資源學(xué)報(bào), 2017, 32(12): 2043–2054 HUO G P, ZHAO X N, GAO X D, et al. Soil water use and competition in jujube-crop systems in loess hilly region[J]. Journal of Natural Resources, 2017, 32(12): 2043–2054

    [29] CAO X Q, YANG P L, ENGEL B A, et al. The effects of rainfall and irrigation on cherry root water uptake under drip irrigation[J]. Agricultural Water Management, 2018, 197: 9–18

    [30] 王紹飛, 趙西寧, 高曉東, 等. 黃土丘陵區(qū)盛果期蘋果樹(shù)土壤水分利用策略[J]. 林業(yè)科學(xué), 2018, 54(10): 31–38 WANG S F, ZHAO X N, GAO X D, et al. Water use source of apple trees with full productive Agein Loess Hilly Region[J]. Scientia Silvae Sinicae, 2018, 54(10): 31–38

    [31] 趙西寧, 李楠, 高曉東, 等. 基于18O示蹤的不同樹(shù)齡棗樹(shù)土壤水分利用特征分析[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2018, 34(3): 135–142 ZHAO X N, LI N, GAO X D, et al. Characteristics of soil water utilization for different stand ages of jujube trees based on18O tracking[J]. Transactions of the Chinese Society of Agricultural Engineering, 2018, 34(3): 135–142

    [32] HASSELQUIST N J, ALLEN M F. Increasing demands on limited water resources: Consequences for two endangered plants in Amargosa Valley, USA[J]. American Journal of Botany, 2009, 96(3): 620–626

    [33] GU X B, LI Y N, DU Y D. Continuous ridges with film mulching improve soil water content, root growth, seed yield and water use efficiency of winter oilseed rape[J]. Industrial Crops and Products, 2016, 85: 139–148

    [34] 胥生榮, 張恩和, 馬瑞麗, 等. 不同覆蓋措施對(duì)枸杞根系生長(zhǎng)和土壤環(huán)境的影響[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 2018, 26(12): 1802–1810 XU S R, ZHANG E H, MA R L, et al. Effects of mulching patterns on root growth and soil environment of[J]. Chinese Journal of Eco-Agriculture, 2018, 26(12): 1802–1810

    [35] 曹有龍, 何軍. 枸杞栽培學(xué)[M]. 銀川: 陽(yáng)光出版社, 2013: 24–26 CAO Y L, HE J.Cultivated Science[M]. Yinchuan: Sunshine Press, 2013: 24–26

    [36] 何昕孺, 戴國(guó)禮, 焦恩寧, 等. 限根栽培對(duì)枸杞根域溫度、生物量積累及營(yíng)養(yǎng)元素吸收的影響[J]. 西北農(nóng)業(yè)學(xué)報(bào), 2017, 26(2): 281–286 HE X R, DAI G L, JIAO E N, et al. Effects of root-limited cultivation on root zone temperature, biomass accumulation and nutrients absorption of wolfberry[J]. Acta Agriculturae Boreali-occidentalis Sinica, 2017, 26(2): 281–286

    [37] 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[J]. Oecologia, 1996, 107(1): 13–20

    [38] 孫守家, 孟平, 張勁松, 等. 利用氘同位素研究太行山南麓棗樹(shù)水分利用的季節(jié)性變化[J]. 林業(yè)科學(xué), 2011, 47(5): 46–53 SUN S J, MENG P, ZHANG J S, et al. Seasonal variation in water use ofin the south aspect of Taihang Mountains with deuterium isotope signature[J]. Scientia Silvae Sinicae, 2011, 47(5): 46–53

    [39] WANG J X, GAO X D, ZHOU Y Q, et al. Impact of conservation practices on soil hydrothermal properties and crop water use efficiency in a dry agricultural region of the Tibetan Plateau[J]. Soil and Tillage Research, 2020, 200: 104619

    [40] 胥生榮, 張恩和, 馬瑞麗, 等. 覆蓋對(duì)枸杞根系土壤環(huán)境和水分利用的影響[J]. 草業(yè)學(xué)報(bào), 2019, 28(2): 12–22 XU S R, ZHANG E H, MA R L, et al. Effects of mulching on soil environment and water utilization by roots of[J]. Acta Prataculturae Sinica, 2019, 28(2): 12–22

    [41] KULMATISKI A, ADLER P B, STARK J M, et al. Water and nitrogen uptake are better associated with resource availability than root biomass[J]. Ecosphere, 2017, 8(3): e01738

    [42] 郭京衡, 李嘗君, 曾凡江, 等. 2種荒漠植物根系生物量分布與土壤水分、養(yǎng)分的關(guān)系[J]. 干旱區(qū)研究, 2016, 33(1): 166–171 GUO J H, LI C J, ZENG F J, et al. Relationship between root biomass distribution and soil moisture, nutrient for two desert plant spices[J]. Arid Zone Research, 2016, 33(1): 166–171

    [43] VAN DER HEIJDEN G, DAMBRINE E, POLLIER B, et al. Mg and Ca uptake by roots in relation to depth and allocation to aboveground tissues: Results from an isotopic labeling study in a beech forest on base-poor soil[J]. Biogeochemistry, 2015, 122(2/3): 375–393

    [44] KULMATISKI A, BEARD K H. Root niche partitioning among grasses, saplings, and trees measured using a tracer technique[J]. Oecologia, 2013, 171(1): 25–37

    [45] HAMZEI J. Seed, oil, and protein yields of canola under combinations of irrigation and nitrogen application[J]. Agronomy Journal, 2011, 103(4): 1152–1158

    [46] 王自奎, 吳普特, 趙西寧, 等. 模擬壟溝灌溉土壤水分入滲特性試驗(yàn)研究[J]. 干旱地區(qū)農(nóng)業(yè)研究, 2011, 29(3): 24–28 WANG Z K, WU P T, ZHAO X N, et al. Simulation experiment on soil water infiltration characteristics under ridge-furrow irrigation[J]. Agricultural Research in the Arid Areas, 2011, 29(3): 24–28

    root water uptake characteristics in the Qaidam Basin irrigation*

    ZHOU Yanqing1,2, GAO Xiaodong2,3, WANG Jiaxin1,2, ZHAO Xining2,3**

    (1. College of Water Resources and Architectural Engineering, Northwest A & F University, Yangling 712100, China; 2. Key Laboratory of Agricultural Soil and Water Engineering in Arid and Semiarid Areas, Ministry of Education, Northwest A & F University, Yangling 712100, China; 3. Institute of Soil and Water Conservation, Chinese Academy of Sciences and Ministry of Water Resources, Yangling 712100, China)

    (Chinese wolfberry) helps ecosystems by providing storm protection and sand immobilization, and it is also a cash crop for Qaidam Basin farmers.is a common crop in arid regions, such as the Qaidam Basin, China, but drought and water scarcity have constrained industry development. Identifying crop water-use strategies is important for designing efficient irrigation systems. Seasonalwater-use was investigated in 2018 in the Qaidam Basin Huaiten Tula Irrigation Area. Natural oxygen stable isotope tracers were used to measure the oxygen stable isotope composition (18O) of water in the xylem and soil oforchardswith three field management practices (conventional flat planting, CK; flat-planting full-film mulching, MF; and ridge-furrow full-film mulching, MR). The soil water contribution to root water uptake was analyzed using the MixSIAR Bayesian mixing practices. The soil water and soil water18O profile distribution showed that all soil layers (shallow: 0–20 cm; middle: 20–60 cm; and deep: 60–100 cm) were potential water sources for. During the sprouting period (starting May 25, 2018), CK and MF plants used primarily shallow soil water, accounting for 45.9% and 37.7% of the total water use, respectively, due to the larger amounts of shallow-layer soil water. MR plants used the same amount of water from each layer. During the blossom and fruiting period (starting July 9, 2018), the soil water content increased in all soil layers. Compared with that in CK plants, the shallow-layer soil water use increased by 13.5% in MF plants and decreased by 11.1% in MR plants. During the fruit ripening period (starting July 30, 2018), water consumption increased. Compared with that in CK plants, shallow-layer soil water use increased in MF (by 11.7%) and MR (by 24.0%) plants. During the defoliation period (starting September 22, 2018), the water content in the shallow-layer soil was lower than that in the other soil layers, and all of the treatment groups used primarily deep-layer soil water. These results showed that shallow-layer soil contributions were positively correlated with water content, indicating thatis sensitive to changes in the shallow-layer soil water. There were significant differences in the root water uptake among the three field management practices. Mulching and ridge-furrow treatments increased the soil water content and, compared with that in CK plants, the MF and MR treatments increased the shallow-layer soil water use. The soil water content was higher and the water source was more flexible during the growth period with the MR treatment than with the CK and MF treatments. These results suggest that MR is a better field water management practice for promoting the sustainable and healthy development oforchards.

    ; Arid oasis irrigation area; Root water uptake; Soil water; Management practices; MixSIAR model

    10.13930/j.cnki.cjea.200435

    周艷清, 高曉東, 王嘉昕, 趙西寧. 柴達(dá)木盆地灌區(qū)枸杞根系水分吸收來(lái)源研究[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)(中英文), 2021, 29(2): 400-409

    ZHOU Y Q, GAO X D, WANG J X, ZHAO X N.root water uptake characteristics in the Qaidam Basin irrigation[J]. Chinese Journal of Eco-Agriculture, 2021, 29(2): 400-409

    Q948.1

    * 國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2017YFC0403605)資助

    趙西寧, 主要研究方向?yàn)檗r(nóng)業(yè)水土資源利用與調(diào)控。E-mail: xiningz@aliyun.com

    周艷清, 主要研究方向?yàn)檗r(nóng)業(yè)水土資源綜合利用。E-mail: zhouyanq021@163.com

    2020-06-08

    2020-08-03

    * This study was supported by the National Key Research and Development Project of China (2017YFC0403605).

    , E-mail: xiningz@aliyun.com

    Jun. 8, 2020;

    Aug. 3, 2020

    猜你喜歡
    利用
    利用min{a,b}的積分表示解決一類絕對(duì)值不等式
    利用倒推破難點(diǎn)
    如何利用基本不等式比較大小
    利用一半進(jìn)行移多補(bǔ)少
    利用口訣算除法
    利用數(shù)的分解來(lái)思考
    Roommate is necessary when far away from home
    利用
    回收木再利用——Piet Hein Eek
    低丘緩坡未利用地的開(kāi)發(fā)利用探討
    河北遙感(2015年4期)2015-07-18 11:05:06
    黄色视频不卡| 老司机深夜福利视频在线观看| 9色porny在线观看| 久久久久视频综合| 大型av网站在线播放| 99精国产麻豆久久婷婷| 久久精品aⅴ一区二区三区四区| 在线亚洲精品国产二区图片欧美| 精品高清国产在线一区| 黄色丝袜av网址大全| 亚洲成a人片在线一区二区| 51午夜福利影视在线观看| 国产午夜精品久久久久久| 捣出白浆h1v1| 亚洲美女黄片视频| 久久久久久免费高清国产稀缺| 久久精品亚洲精品国产色婷小说| 色综合欧美亚洲国产小说| 最近最新中文字幕大全免费视频| 色尼玛亚洲综合影院| 国产男女超爽视频在线观看| 精品高清国产在线一区| 99riav亚洲国产免费| 午夜福利视频在线观看免费| 中文字幕av电影在线播放| 国产av又大| 老司机影院毛片| 国产亚洲精品第一综合不卡| av电影中文网址| 国产一区二区三区在线臀色熟女 | 亚洲中文字幕日韩| 大型av网站在线播放| 精品视频人人做人人爽| 国产精品熟女久久久久浪| 成人国语在线视频| 免费女性裸体啪啪无遮挡网站| 国产精品免费一区二区三区在线 | 久久天堂一区二区三区四区| 欧美大码av| 亚洲黑人精品在线| 午夜福利欧美成人| 国产视频一区二区在线看| 午夜久久久在线观看| 母亲3免费完整高清在线观看| 人人妻人人澡人人爽人人夜夜| 国产真人三级小视频在线观看| 亚洲精品一卡2卡三卡4卡5卡| 午夜久久久在线观看| 国产区一区二久久| 99riav亚洲国产免费| 国产免费现黄频在线看| 欧美在线黄色| 亚洲国产成人一精品久久久| 在线观看免费高清a一片| 国产成人欧美在线观看 | 亚洲美女黄片视频| 免费在线观看视频国产中文字幕亚洲| 午夜老司机福利片| 中亚洲国语对白在线视频| 天堂动漫精品| 国产高清视频在线播放一区| 咕卡用的链子| 国产日韩一区二区三区精品不卡| 亚洲人成伊人成综合网2020| 五月开心婷婷网| 亚洲少妇的诱惑av| 动漫黄色视频在线观看| 亚洲精品国产色婷婷电影| 日本一区二区免费在线视频| 欧美日韩成人在线一区二区| 99精品欧美一区二区三区四区| 一个人免费看片子| 人人妻,人人澡人人爽秒播| 黑人巨大精品欧美一区二区mp4| 久久久久国内视频| 叶爱在线成人免费视频播放| 国产精品麻豆人妻色哟哟久久| 变态另类成人亚洲欧美熟女 | 亚洲成人免费av在线播放| 搡老乐熟女国产| 亚洲天堂av无毛| 露出奶头的视频| 日韩精品免费视频一区二区三区| 少妇 在线观看| 国产成人精品久久二区二区免费| 久久这里只有精品19| 久久久久久久久久久久大奶| 老汉色av国产亚洲站长工具| 91麻豆精品激情在线观看国产 | 新久久久久国产一级毛片| 中文字幕精品免费在线观看视频| 看免费av毛片| 精品福利永久在线观看| 九色亚洲精品在线播放| 亚洲成a人片在线一区二区| 肉色欧美久久久久久久蜜桃| 亚洲精品一二三| 亚洲美女黄片视频| 国产精品久久久久久精品电影小说| 久久久久网色| 岛国在线观看网站| 国产精品电影一区二区三区 | 久久精品亚洲精品国产色婷小说| 国产精品欧美亚洲77777| 久久精品91无色码中文字幕| 亚洲av电影在线进入| 国产不卡一卡二| av免费在线观看网站| 日日夜夜操网爽| 久久性视频一级片| 国产成人精品在线电影| 欧美老熟妇乱子伦牲交| 成人特级黄色片久久久久久久 | 色综合欧美亚洲国产小说| 久久中文看片网| 飞空精品影院首页| 亚洲三区欧美一区| 大片电影免费在线观看免费| 欧美日韩一级在线毛片| 免费观看av网站的网址| 色综合婷婷激情| 日韩大片免费观看网站| 久久久久久久久久久久大奶| 亚洲精品国产精品久久久不卡| 男女无遮挡免费网站观看| 成人永久免费在线观看视频 | 大片电影免费在线观看免费| 这个男人来自地球电影免费观看| 亚洲九九香蕉| 黄片小视频在线播放| 欧美日韩亚洲高清精品| 国产一区二区 视频在线| 国产精品美女特级片免费视频播放器 | 人人妻,人人澡人人爽秒播| 亚洲中文字幕日韩| 黄片播放在线免费| 在线观看免费日韩欧美大片| 99国产综合亚洲精品| 久久久久网色| 亚洲情色 制服丝袜| 欧美日韩福利视频一区二区| 国产精品av久久久久免费| 深夜精品福利| 巨乳人妻的诱惑在线观看| 亚洲精品在线观看二区| 免费在线观看影片大全网站| 成人国语在线视频| 日本撒尿小便嘘嘘汇集6| 午夜精品久久久久久毛片777| 嫁个100分男人电影在线观看| 交换朋友夫妻互换小说| 午夜福利,免费看| av视频免费观看在线观看| 丁香欧美五月| 九色亚洲精品在线播放| 免费日韩欧美在线观看| 黄色毛片三级朝国网站| 黄色毛片三级朝国网站| 欧美精品av麻豆av| 国产成人欧美在线观看 | 亚洲欧美激情在线| 啦啦啦视频在线资源免费观看| 亚洲专区字幕在线| 黑人操中国人逼视频| 日日夜夜操网爽| 大片电影免费在线观看免费| 午夜老司机福利片| 亚洲天堂av无毛| 电影成人av| 午夜福利免费观看在线| 日韩三级视频一区二区三区| 国产精品二区激情视频| 欧美变态另类bdsm刘玥| 久久久久久久国产电影| 欧美日韩av久久| 国产亚洲精品一区二区www | 国产av一区二区精品久久| 亚洲精品久久成人aⅴ小说| 午夜福利欧美成人| 欧美在线黄色| 国产亚洲一区二区精品| 日韩免费高清中文字幕av| 一本大道久久a久久精品| 女人爽到高潮嗷嗷叫在线视频| 午夜福利视频精品| 脱女人内裤的视频| svipshipincom国产片| 操出白浆在线播放| tocl精华| 欧美日韩成人在线一区二区| 久久久久久久精品吃奶| 日韩欧美国产一区二区入口| 国产成人欧美| 亚洲国产欧美日韩在线播放| 国产精品熟女久久久久浪| 欧美乱码精品一区二区三区| 久久ye,这里只有精品| 欧美成人免费av一区二区三区 | 国产片内射在线| 两人在一起打扑克的视频| 国产欧美日韩一区二区精品| 国产精品久久久久久人妻精品电影 | 欧美精品高潮呻吟av久久| 老司机深夜福利视频在线观看| 9191精品国产免费久久| 91麻豆精品激情在线观看国产 | 人人妻,人人澡人人爽秒播| 亚洲五月色婷婷综合| 亚洲中文字幕一区二区三区有码在线看 | 成人18禁在线播放| 久久精品91无色码中文字幕| 毛片女人毛片| 18禁裸乳无遮挡免费网站照片| 久久天躁狠狠躁夜夜2o2o| 精品免费久久久久久久清纯| 成人午夜高清在线视频| 国产黄色小视频在线观看| 久久中文字幕人妻熟女| 国产三级黄色录像| 97超视频在线观看视频| 精品不卡国产一区二区三区| 国产探花在线观看一区二区| 人人妻人人看人人澡| 国产熟女xx| 脱女人内裤的视频| 精品久久久久久久人妻蜜臀av| 亚洲av片天天在线观看| 搞女人的毛片| 久久伊人香网站| 国产久久久一区二区三区| 天堂影院成人在线观看| 久久人妻av系列| 国产精品 欧美亚洲| 久久久久精品国产欧美久久久| 黄片大片在线免费观看| 欧美黑人巨大hd| 十八禁人妻一区二区| 久久天躁狠狠躁夜夜2o2o| 亚洲欧美日韩卡通动漫| 在线免费观看的www视频| www日本黄色视频网| 欧美日韩精品网址| 国产乱人伦免费视频| 亚洲电影在线观看av| 手机成人av网站| 女同久久另类99精品国产91| 国产精品亚洲美女久久久| 亚洲成av人片在线播放无| 夜夜夜夜夜久久久久| 久久草成人影院| 中文字幕精品亚洲无线码一区| 十八禁人妻一区二区| 国产美女午夜福利| 久久精品国产99精品国产亚洲性色| 久99久视频精品免费| 少妇裸体淫交视频免费看高清| 波多野结衣高清作品| 在线视频色国产色| 亚洲真实伦在线观看| 成人国产一区最新在线观看| 久久人人精品亚洲av| 日韩中文字幕欧美一区二区| 999久久久精品免费观看国产| 亚洲成人免费电影在线观看| 亚洲人成网站在线播放欧美日韩| 女同久久另类99精品国产91| 久久久久免费精品人妻一区二区| 99久久99久久久精品蜜桃| 最新美女视频免费是黄的| 亚洲天堂国产精品一区在线| 九色国产91popny在线| 看免费av毛片| 天堂av国产一区二区熟女人妻| 窝窝影院91人妻| 特级一级黄色大片| 不卡av一区二区三区| 超碰成人久久| 高清在线国产一区| 欧美另类亚洲清纯唯美| 怎么达到女性高潮| 脱女人内裤的视频| 在线观看舔阴道视频| 国产亚洲精品综合一区在线观看| 亚洲一区高清亚洲精品| 最新在线观看一区二区三区| 亚洲av电影不卡..在线观看| www.熟女人妻精品国产| 欧美最黄视频在线播放免费| 99热精品在线国产| 久久久久国产精品人妻aⅴ院| 国语自产精品视频在线第100页| 观看免费一级毛片| 亚洲国产日韩欧美精品在线观看 | 国产午夜福利久久久久久| 欧美日本视频| 欧美不卡视频在线免费观看| 亚洲精品一卡2卡三卡4卡5卡| 变态另类成人亚洲欧美熟女| 欧美一区二区精品小视频在线| 成人三级黄色视频| 午夜免费观看网址| 久久久久久久午夜电影| 国产成人精品无人区| 精品一区二区三区视频在线 | 午夜a级毛片| av福利片在线观看| 一区二区三区高清视频在线| 中文字幕人成人乱码亚洲影| 一进一出好大好爽视频| 黄色女人牲交| 深夜精品福利| 亚洲精品中文字幕一二三四区| 别揉我奶头~嗯~啊~动态视频| 少妇的丰满在线观看| www.熟女人妻精品国产| 婷婷丁香在线五月| 亚洲人与动物交配视频| 久久久国产成人精品二区| 午夜福利成人在线免费观看| 亚洲在线观看片| 免费在线观看视频国产中文字幕亚洲| 又黄又粗又硬又大视频| 亚洲一区二区三区色噜噜| 人人妻人人澡欧美一区二区| 一区福利在线观看| 久久精品影院6| 两个人视频免费观看高清| 亚洲成人久久爱视频| 婷婷丁香在线五月| 国产aⅴ精品一区二区三区波| 久久久久免费精品人妻一区二区| 国产成人精品久久二区二区免费| 欧美日韩国产亚洲二区| 欧美三级亚洲精品| 搡老岳熟女国产| 亚洲美女黄片视频| 久久久水蜜桃国产精品网| 90打野战视频偷拍视频| 亚洲aⅴ乱码一区二区在线播放| 亚洲成人久久爱视频| 亚洲欧洲精品一区二区精品久久久| 国产在线精品亚洲第一网站| 亚洲一区二区三区不卡视频| 成年女人看的毛片在线观看| 两个人的视频大全免费| 国产主播在线观看一区二区| 久久久久久久午夜电影| 毛片女人毛片| 特大巨黑吊av在线直播| 天堂影院成人在线观看| 国产亚洲精品综合一区在线观看| 波多野结衣巨乳人妻| 搡老熟女国产l中国老女人| 久久精品91蜜桃| 欧美又色又爽又黄视频| 国产欧美日韩精品亚洲av| 俄罗斯特黄特色一大片| 国产亚洲av嫩草精品影院| 特大巨黑吊av在线直播| 级片在线观看| 国产精品99久久久久久久久| 熟妇人妻久久中文字幕3abv| 麻豆久久精品国产亚洲av| 丰满的人妻完整版| 久久久国产成人免费| 欧美日韩国产亚洲二区| 国产一区在线观看成人免费| a级毛片在线看网站| 中文字幕最新亚洲高清| 老汉色∧v一级毛片| 视频区欧美日本亚洲| 成年女人毛片免费观看观看9| 成人国产综合亚洲| 美女黄网站色视频| 精品欧美国产一区二区三| 免费一级毛片在线播放高清视频| 麻豆久久精品国产亚洲av| 九九热线精品视视频播放| 老司机午夜十八禁免费视频| 搞女人的毛片| 亚洲av成人不卡在线观看播放网| 亚洲国产精品sss在线观看| 国产又色又爽无遮挡免费看| 在线观看美女被高潮喷水网站 | 国产一级毛片七仙女欲春2| 男人舔女人的私密视频| 久久中文看片网| 亚洲欧美日韩卡通动漫| 女人被狂操c到高潮| 啦啦啦韩国在线观看视频| 国产精品亚洲一级av第二区| 村上凉子中文字幕在线| 国产欧美日韩一区二区三| 黄频高清免费视频| 日本黄色片子视频| 一级毛片高清免费大全| 国产视频一区二区在线看| 国产免费av片在线观看野外av| 成人特级av手机在线观看| 亚洲欧美日韩高清专用| 非洲黑人性xxxx精品又粗又长| 欧美日韩福利视频一区二区| 午夜精品一区二区三区免费看| 老熟妇乱子伦视频在线观看| 日韩精品青青久久久久久| 哪里可以看免费的av片| 给我免费播放毛片高清在线观看| 亚洲av电影不卡..在线观看| www.熟女人妻精品国产| 亚洲av五月六月丁香网| 国产野战对白在线观看| 18禁国产床啪视频网站| 啪啪无遮挡十八禁网站| 午夜福利成人在线免费观看| 99久久精品一区二区三区| 国内精品美女久久久久久| 亚洲激情在线av| 免费观看人在逋| 在线观看免费午夜福利视频| 国产精品电影一区二区三区| av天堂中文字幕网| 男女下面进入的视频免费午夜| 免费大片18禁| 国产真实乱freesex| 亚洲人成网站高清观看| 老汉色∧v一级毛片| 在线永久观看黄色视频| cao死你这个sao货| 久久久久国产一级毛片高清牌| netflix在线观看网站| 免费看光身美女| 亚洲av成人不卡在线观看播放网| 变态另类丝袜制服| 亚洲aⅴ乱码一区二区在线播放| 很黄的视频免费| 国产伦精品一区二区三区四那| 久久香蕉精品热| 色综合亚洲欧美另类图片| 91av网站免费观看| 亚洲在线观看片| 亚洲国产精品久久男人天堂| 又黄又爽又免费观看的视频| 成年版毛片免费区| 国产精品影院久久| 国产亚洲欧美在线一区二区| 99在线视频只有这里精品首页| 亚洲人成网站高清观看| 99国产精品一区二区三区| 国产又色又爽无遮挡免费看| 亚洲aⅴ乱码一区二区在线播放| 亚洲第一欧美日韩一区二区三区| 精品不卡国产一区二区三区| 美女高潮喷水抽搐中文字幕| 国产亚洲精品久久久com| 久久人妻av系列| 日本黄大片高清| 国产黄a三级三级三级人| 麻豆成人午夜福利视频| 免费观看精品视频网站| 老熟妇仑乱视频hdxx| 久久天堂一区二区三区四区| 亚洲九九香蕉| 桃色一区二区三区在线观看| 99国产精品99久久久久| 午夜免费成人在线视频| 国产精品九九99| 中文字幕精品亚洲无线码一区| 国产成人精品久久二区二区免费| 日韩欧美国产一区二区入口| 久久久久久大精品| 国产单亲对白刺激| 久久久久免费精品人妻一区二区| 亚洲精品一区av在线观看| 性色av乱码一区二区三区2| 十八禁人妻一区二区| 桃色一区二区三区在线观看| 欧美日本视频| 美女免费视频网站| 51午夜福利影视在线观看| 99热6这里只有精品| netflix在线观看网站| av欧美777| 久久久久国内视频| 亚洲五月婷婷丁香| 国产又黄又爽又无遮挡在线| 久久香蕉精品热| 久久久成人免费电影| 亚洲欧美日韩高清专用| 99精品在免费线老司机午夜| 色在线成人网| 18禁裸乳无遮挡免费网站照片| 中文字幕最新亚洲高清| 全区人妻精品视频| 国产黄a三级三级三级人| 一级作爱视频免费观看| 国产精品亚洲美女久久久| 久久精品aⅴ一区二区三区四区| 搡老妇女老女人老熟妇| 欧美中文综合在线视频| 一二三四在线观看免费中文在| 久久天堂一区二区三区四区| 好男人在线观看高清免费视频| av黄色大香蕉| 三级国产精品欧美在线观看 | 日本与韩国留学比较| 日韩av在线大香蕉| 亚洲中文日韩欧美视频| 国产精品一及| 亚洲av电影不卡..在线观看| 久久久国产欧美日韩av| 欧美成狂野欧美在线观看| 亚洲av五月六月丁香网| 免费大片18禁| 啦啦啦韩国在线观看视频| 国产免费av片在线观看野外av| 亚洲精品456在线播放app | 一进一出抽搐gif免费好疼| 中国美女看黄片| 国内精品一区二区在线观看| 亚洲性夜色夜夜综合| 免费一级毛片在线播放高清视频| 国产私拍福利视频在线观看| 人妻久久中文字幕网| 国产野战对白在线观看| 午夜精品在线福利| 女同久久另类99精品国产91| 岛国在线免费视频观看| 搞女人的毛片| 久久久国产精品麻豆| 中文字幕久久专区| 免费av毛片视频| 老司机午夜福利在线观看视频| 嫩草影视91久久| 日韩有码中文字幕| av在线蜜桃| 高清毛片免费观看视频网站| 婷婷亚洲欧美| 俄罗斯特黄特色一大片| 午夜精品一区二区三区免费看| 国产精品久久久久久精品电影| 久久久久国产一级毛片高清牌| 无遮挡黄片免费观看| 99久久国产精品久久久| 一进一出抽搐gif免费好疼| 男人和女人高潮做爰伦理| 日韩中文字幕欧美一区二区| av视频在线观看入口| 久久精品亚洲精品国产色婷小说| 免费高清视频大片| 在线播放国产精品三级| 色综合亚洲欧美另类图片| 级片在线观看| 国产精品 国内视频| 国产亚洲精品综合一区在线观看| 99精品欧美一区二区三区四区| 最近最新中文字幕大全电影3| 久久久成人免费电影| 亚洲精品456在线播放app | 久久久久久九九精品二区国产| 舔av片在线| 国产v大片淫在线免费观看| 精品福利观看| 久久中文看片网| 免费在线观看成人毛片| www日本在线高清视频| 男女下面进入的视频免费午夜| 禁无遮挡网站| 小蜜桃在线观看免费完整版高清| 免费电影在线观看免费观看| 欧美中文综合在线视频| 日韩成人在线观看一区二区三区| 首页视频小说图片口味搜索| 午夜免费成人在线视频| 97超级碰碰碰精品色视频在线观看| 亚洲 国产 在线| 亚洲欧美激情综合另类| av天堂在线播放| 在线看三级毛片| 村上凉子中文字幕在线| 国产成人欧美在线观看| 久久久国产欧美日韩av| 性欧美人与动物交配| 国产真人三级小视频在线观看| 91麻豆av在线| 巨乳人妻的诱惑在线观看| 人妻久久中文字幕网| 亚洲国产看品久久| 欧洲精品卡2卡3卡4卡5卡区| 91麻豆精品激情在线观看国产| www.精华液| 丰满人妻一区二区三区视频av | 久久精品91无色码中文字幕| av天堂中文字幕网| 18禁裸乳无遮挡免费网站照片| 国产精品久久久久久久电影 | 成年人黄色毛片网站| 久久久久久九九精品二区国产| 国产成人av教育| 巨乳人妻的诱惑在线观看| 成年女人毛片免费观看观看9| 啦啦啦免费观看视频1| 亚洲国产中文字幕在线视频| 香蕉av资源在线| 日韩大尺度精品在线看网址| 亚洲精品国产精品久久久不卡| 亚洲av中文字字幕乱码综合| a级毛片a级免费在线| 韩国av一区二区三区四区| 国产伦人伦偷精品视频| 手机成人av网站| 这个男人来自地球电影免费观看| 精品国产超薄肉色丝袜足j| 亚洲国产色片| 欧美又色又爽又黄视频| 女人被狂操c到高潮|