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

    不同灌溉制度下河套灌區(qū)玉米膜下滴灌水熱鹽運(yùn)移規(guī)律

    2017-07-18 11:48:13孫貫芳屈忠義任中生劉安琪
    關(guān)鍵詞:河套土壤溫度鹽分

    孫貫芳,屈忠義,杜 斌,任中生,劉安琪

    ?

    不同灌溉制度下河套灌區(qū)玉米膜下滴灌水熱鹽運(yùn)移規(guī)律

    孫貫芳1,2,屈忠義1※,杜 斌3,任中生1,劉安琪1

    (1. 內(nèi)蒙古農(nóng)業(yè)大學(xué)水利與土木建筑工程學(xué)院,呼和浩特010018;2. 武漢大學(xué)水資源與水電工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,武漢 430072;3. 寧夏水利科學(xué)研究院,銀川750000)

    為探明干旱地區(qū)膜下滴灌玉米土壤水熱鹽效應(yīng)及秋澆洗鹽灌溉的效果,該文根據(jù)2014—2015年進(jìn)行的田間試驗(yàn),分析膜下滴灌不同灌溉制度下生育期土壤水分鹽分剖面分布特性、土壤溫度變化及對(duì)玉米產(chǎn)量品質(zhì)的影響和非生育期洗鹽灌溉(秋澆)效果。結(jié)果表明,不同滴灌制度下土壤剖面水分、鹽分剖面分布極不均勻,鹽分均由膜內(nèi)向膜外地表裸露區(qū)定向遷移,趨于膜外地表積累。膜下滴灌土壤溫度受氣溫、玉米葉面積指數(shù)、灌水及土壤含水率共同作用。灌水后各處理土壤溫度均劇烈下降,2~3 d后恢復(fù),玉米營(yíng)養(yǎng)生長(zhǎng)階段控制灌水下限為-30 kPa最有利于土壤溫度積累。玉米生育期各處理膜內(nèi)0~40 cm不積鹽,控制灌水下限為-10 kPa可有效淋濾0~100 cm土壤鹽分,而其他處理對(duì)0~100 cm土層鹽分的影響差異性短期內(nèi)不明顯,需對(duì)不同處理長(zhǎng)期膜下滴灌的鹽分進(jìn)一步觀測(cè)。非生育期洗鹽灌溉效果顯著,秋澆灌黃河水180 mm后,次年春播前0~100 cm土壤鹽分平均下降10.86%~26.14%,剖面分布較均勻。河套灌區(qū)膜下滴灌土壤鹽分調(diào)控建議為生育期滴灌灌溉制度和非生育期洗鹽灌溉雙重調(diào)控。玉米生育期灌水下限建議控制為-30 kPa,非生育期洗鹽灌溉由于河套灌區(qū)凍融影響及特殊的水文地質(zhì)條件,膜下滴灌鹽分累積到何種程度洗鹽灌溉及具體合理的洗鹽灌溉制度還需進(jìn)一步深入研究。

    土壤水分;溫度;含鹽量;河套灌區(qū); 膜下滴灌; 水熱鹽; 洗鹽; 秋澆

    0 引 言

    內(nèi)蒙古河套灌區(qū)地處干旱地區(qū),目前農(nóng)業(yè)灌溉引水量約43億m3~48億m3,農(nóng)田灌溉平均凈引水量45億m3[1],灌區(qū)內(nèi)灌水方式基本為地面灌溉,田間灌溉水利用效率低,灌溉水大量回補(bǔ)地下水[2],造成土壤次生鹽漬化。滴灌在提高灌溉水利用效率的同時(shí),可實(shí)現(xiàn)高產(chǎn)高效,且適宜的滴灌制度已被證明是在干旱區(qū)緩解土壤鹽漬化的有效辦法[3-5]。受不同地區(qū)氣候、土壤特性、作物品種、農(nóng)藝措施及作物栽培模式的影響[6-8],膜下滴灌技術(shù)具有很強(qiáng)的地域性,在河套灌區(qū)膜下滴灌關(guān)注最多的是經(jīng)濟(jì)效益較高的加工番茄[9-10]。玉米作為灌區(qū)內(nèi)最主要的糧食作物,近2 a剛引起人們的關(guān)注,且成果集中在膜下滴灌玉米耗水規(guī)律及水分生產(chǎn)率[11-12]、作物系數(shù)[13]、根系分布[14]等方面。河套灌區(qū)干旱少雨,土壤蒸發(fā)強(qiáng)烈,易造成土壤鹽分表聚和次生鹽漬化[15],在發(fā)展滴灌的過(guò)程中,土壤鹽分的問(wèn)題必須足夠關(guān)注。

    已有研究表明,土壤鹽分運(yùn)動(dòng)受制于土壤水分和溫度[16],干旱區(qū)膜下滴灌的微區(qū)環(huán)境是水鹽定向遷移的內(nèi)在機(jī)制,氣溫與蒸發(fā)的交互作用則進(jìn)一步增強(qiáng)了鹽分的定向遷移[17]。因此,周和平等[18]提出地表排鹽模式,經(jīng)研究發(fā)現(xiàn)采用田間兩膜之間的裸露地表排鹽溝會(huì)使鹽分下降36.7%~63.3%,效果良好,但其必須要求鹽分向農(nóng)田下游末端排泄有出路,在河套灌區(qū)真正實(shí)施有一定困難。Liu等[19]研究發(fā)現(xiàn),用電導(dǎo)率達(dá)到7.42 dS/m的微咸水進(jìn)行膜下滴灌,用150 mm秋澆定額淋洗后,土壤鹽分淋洗到60 cm以下土層,作物根區(qū)土壤電導(dǎo)率在下一年僅為0.2 dS/m,說(shuō)明秋澆能很好地控制滴灌土壤鹽分。同時(shí),Phocaides[20]也認(rèn)為在年降水量小于250 mm的區(qū)域采用微灌不推薦在每次灌水時(shí)加大灌水量進(jìn)行洗鹽,而是建議在生育期結(jié)束后對(duì)鹽分集中淋洗。非生育期洗鹽灌溉可能是控制膜下滴灌土壤鹽分的較佳途徑。課題組李金剛等[21]在河套灌區(qū)下游的長(zhǎng)勝試驗(yàn)站開(kāi)展了微咸水膜下滴灌鹽堿地玉米的大田試驗(yàn),詳細(xì)分析了生育期膜下滴灌膜內(nèi)膜外水鹽運(yùn)移過(guò)程,發(fā)現(xiàn)0~100 cm土層生育期鹽分不斷累積,春匯225 mm黃河水后各土層鹽分降低21.46%以上,對(duì)鹽堿地微咸水膜下滴灌的實(shí)踐具有一定的指導(dǎo)意義,但未深入探討膜下滴灌膜內(nèi)膜外土壤鹽分分布不均的機(jī)制及鹽分在各土層的平衡規(guī)律。另外,河套灌區(qū)屬于典型的寒旱灌區(qū),晝夜溫差大,土壤積溫對(duì)作物生長(zhǎng)至關(guān)重要,膜下滴灌不同灌水下限對(duì)土壤溫度的影響尚不清楚。因此,河套灌區(qū)膜下滴灌作物生育期灌水下限的選擇需綜合考慮土壤水熱鹽狀況及作物產(chǎn)量品質(zhì)等因素,所積累的土壤鹽分在非生育期通過(guò)地面灌溉集中淋洗可能是未來(lái)灌區(qū)控制膜下滴灌土壤鹽分的關(guān)鍵措施。本文通過(guò)田間試驗(yàn),研究河套灌區(qū)膜下滴灌生育期不同灌水下限(分別為-10、-20、-30、-40 kPa)下的土壤水熱鹽效應(yīng)及對(duì)玉米產(chǎn)量品質(zhì)的影響,探索不同灌水下限及非生育期秋澆洗鹽對(duì)土壤鹽分平衡的影響,以期確定河套灌區(qū)玉米膜下滴灌合理的生育期灌溉制度和適宜的鹽分調(diào)控模式,為膜下滴灌土壤水熱鹽的調(diào)控提供理論支持。

    1 材料與方法

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

    試驗(yàn)于2014年5月—2015年5月在內(nèi)蒙古河套灌區(qū)臨河區(qū)雙河鎮(zhèn)進(jìn)步村九莊農(nóng)業(yè)合作社試驗(yàn)基地(107°18′E,40°41¢N)進(jìn)行,試驗(yàn)區(qū)布置如圖1所示。

    該地深處內(nèi)陸,屬于中溫帶半干旱大陸性氣候,多年平均降水量140 mm,平均氣溫6.8 ℃,晝夜溫差大,日照時(shí)間長(zhǎng),多年日照時(shí)間平均值為3 229.9 h,是中國(guó)日照時(shí)數(shù)較長(zhǎng)的地區(qū)之一。光、熱、水同期,無(wú)霜期為130 d左右,適宜于農(nóng)作物生長(zhǎng)。該地一般每年11月中旬土壤開(kāi)始封凍,次年5月上旬融通。試驗(yàn)區(qū)以粉砂壤土為主,0~100 cm土壤容重為1.38 g/cm3(表1),平均田間持水率為28.5%,土壤全氮量、全磷量、全鉀量(質(zhì)量比)分別為0.093%、0.07%、1.60%、有機(jī)質(zhì)質(zhì)量比為1.2%,pH值為7.6。玉米生育期地下水埋深范圍為1.89~3.08 m,平均埋深2.84 m,非生育期受秋澆及凍融的影響,地下水埋深變化較為劇烈,平均為2.36 m(見(jiàn)圖2)。膜下滴灌所用的地下水含鹽量為1.07 g/L,pH值為8.0,Ca2+、Mg2+、Na++K+、Cl-、SO42-、CO32-、HCO4-質(zhì)量濃度分別為70.1、72.9、150、159.5、216.1、0、396.6 mg/L。玉米生育期(2014年5月1日—9月30日)降水量為104.6 mm,非生育期(2014年10月1日—2015年4月30日)降水量為62.8 mm(見(jiàn)圖3)。試驗(yàn)期間參考作物蒸散量(reference crop evapotranspiration,ET0)根據(jù)氣象數(shù)據(jù)用FAO-56 Penman- Monteith公式[22]計(jì)算,玉米生育期ET0為765.84 mm,非生育期為383.26 mm。

    表1 試驗(yàn)區(qū)土壤物理及水力參數(shù)

    圖3 試驗(yàn)期間降水量及參考作物蒸散量ET0

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

    膜下滴灌玉米品種為內(nèi)單314,采用1膜1管(滴灌帶)2行種植方式(見(jiàn)圖4),滴灌帶鋪設(shè)間距120 cm,地膜寬70 cm,玉米寬窄行種植,寬行70 cm,窄行50 cm,株距22 cm,密度75 000株/hm2。采用耐特菲姆內(nèi)鑲貼片式滴灌帶,管內(nèi)徑16 mm,流量1.60 L/h,滴頭間距30 cm。膜下滴灌玉米5月1日播種后統(tǒng)一灌出苗水30 mm,之后采用北京奧特思達(dá)科技有限公司生產(chǎn)的張力計(jì)控制灌溉,參考Kang等[23]研究成果,張力計(jì)埋設(shè)在膜內(nèi)滴頭下20 cm處(圖4),分別設(shè)置基質(zhì)勢(shì)為-10、-20、-30、-40 kPa 4個(gè)灌水下限處理,簡(jiǎn)記為D1、D2、D3、D4,每個(gè)處理為一個(gè)灌水單元(長(zhǎng)25 m,寬7.2 m,面積180 m3),當(dāng)張力計(jì)讀數(shù)達(dá)到相應(yīng)值時(shí)開(kāi)始灌水,灌溉水采用當(dāng)?shù)氐叵滤?,灌溉水的水溫一般?~12 ℃,設(shè)計(jì)灌水定額為30 mm,D1處理5、6、7、8、9月灌水次數(shù)分別為2、4、8、3、0次,D2處理6月、7月較D1各少灌1次水,D3處理7月、8月較D2少灌2次、1次水,D4處理5—7月較D3各少灌1次水,4個(gè)處理總灌水次數(shù)分別為17、15、12和9次,實(shí)際灌溉定額分別為518.0、444.7、368.3、268.3 mm。各處理的施肥制度為播種時(shí)統(tǒng)一施入底肥磷酸二銨600 kg/hm2,后期隨滴灌追尿素600 kg/hm2,鉀肥(硝酸鉀)90 kg/ hm2。其他的田間管理措施(除草,打藥等)各試驗(yàn)處理均相同,同當(dāng)?shù)亓?xí)慣。

    圖4 張力計(jì)、溫度傳感器埋設(shè)位置及土壤水鹽取樣位置示意圖

    參考李瑞平等[24]研究成果,2014年秋澆期不同滴灌制度試驗(yàn)小區(qū)于10月12日均灌水180 mm。具體過(guò)程為:各處理間用圍堰隔開(kāi),每個(gè)小區(qū)灌黃河水32.4 m3。于2015年5月6日播前取土測(cè)試土壤鹽分。

    1.3 觀測(cè)項(xiàng)目與方法

    1.3.1 土壤水分

    7月下旬玉米抽雄期灌水頻繁,有利于分析土壤剖面水分狀況。各處理于灌水前分別在垂直滴灌帶距滴頭水平距離0、17.5、35、60 cm位置取土,垂直取樣深度為0~60 cm每10 cm一層,>60~100 cm每20 cm一層,取樣點(diǎn)設(shè)置如圖4所示。各3次重復(fù),用烘箱105℃烘干測(cè)定含水率。

    1.3.2 土壤及大氣溫度

    土壤溫度采用邯鄲市益盟電子有限公司生產(chǎn)的地埋式8通路溫度記錄儀(YM-04)實(shí)時(shí)觀測(cè),溫度計(jì)埋設(shè)位置分別為膜內(nèi)滴頭正下方5、15、25、40 cm,膜外(即指兩地膜之間裸地正中間位置)5、15、25、40 cm處(見(jiàn)圖4)。土壤溫度記錄儀精度±0.02 ℃,溫度分辨率0.05 ℃。6月1日20:00點(diǎn)開(kāi)始自動(dòng)記錄,每1 h記錄1次溫度,9月24日結(jié)束。土壤溫度不設(shè)重復(fù)。大氣溫度采用試驗(yàn)區(qū)內(nèi)HOBO小型氣象站(美國(guó)Onset公司產(chǎn)品)溫度,每6 min記錄1次數(shù)據(jù)。

    1.3.3 土壤鹽分

    分別在播前、生育期末收獲后、次年播前取土測(cè)試土壤鹽分。播前鹽分取土方法為:土層深度0~60 cm每10 cm一層,>60~100 cm每20 cm一層,后2次取土方法同上述水分測(cè)試取土方法。風(fēng)干土經(jīng)碾磨、過(guò)2 mm篩后,用電導(dǎo)率儀(上海雷磁電導(dǎo)率儀DDSJ-308A)測(cè)試水土比為5∶1的土壤浸提液的電導(dǎo)率,后經(jīng)該試驗(yàn)區(qū)土壤含鹽率(%)與電導(dǎo)率(electrical conductivity,EC,dS/m)經(jīng)驗(yàn)公式=0.349EC計(jì)算土壤全鹽含量[25]。土壤鹽分取3次重復(fù)。

    1.3.4 玉米產(chǎn)量與品質(zhì)

    玉米成熟后,于2014年9月24日分別沿試驗(yàn)小區(qū)上段、中段、下段隨機(jī)選擇3個(gè)位置作為取樣點(diǎn),在中間3膜(試驗(yàn)小區(qū)共6膜)上連續(xù)每行取5株玉米共30株作為1個(gè)重復(fù),取回其玉米棒。先用游標(biāo)卡尺分別測(cè)量30個(gè)玉米棒的禿尖長(zhǎng)度,之后人工脫粒風(fēng)干稱量產(chǎn)量,并換算成單位面積的產(chǎn)量。從每個(gè)重復(fù)中隨機(jī)取100粒玉米籽粒稱量質(zhì)量即為百粒質(zhì)量。玉米產(chǎn)量、禿尖長(zhǎng)度及百粒質(zhì)量均是3個(gè)重復(fù)。玉米品質(zhì)是3個(gè)重復(fù)經(jīng)充分混合為一個(gè)樣后委托內(nèi)蒙古博測(cè)質(zhì)檢科技有限責(zé)任公司檢測(cè)粗脂肪、粗蛋白、粗淀粉,依據(jù)的標(biāo)準(zhǔn)分別是GB/T 6433-2006[26]、GB/T 6432-1994[27]、NY/T 11-1985[28]。

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

    1.4.1 土壤溫度的滑動(dòng)平均

    土壤溫度的滑動(dòng)平均有利于消除土壤溫度的不穩(wěn)定波動(dòng),把握溫度變化的趨勢(shì),是從一個(gè)有項(xiàng)的時(shí)間序列中來(lái)計(jì)算多個(gè)連續(xù)項(xiàng)序列的平均值。土壤溫度15 d滑動(dòng)平均是在長(zhǎng)序列的逐日土壤溫度資料中從第1天到第15天,第2天到第16天,第3天到第17天……,每相應(yīng)15 d的資料計(jì)算其平均值,具體如公式(1)所示:

    式中T為第1 d~第15 d共15 d的土壤溫度滑動(dòng)平均值,℃,文中的取值為0,1,2……,100;T為計(jì)算的15 d里每天的土壤溫度,℃。

    1.4.2 土壤鹽分平衡分析

    對(duì)1 m土層土壤鹽分進(jìn)行平衡分析是評(píng)價(jià)鹽分累積狀況的重要方法。土壤含鹽總量根據(jù)不同取樣點(diǎn)控制質(zhì)量加權(quán)計(jì)算[29],計(jì)算公式如(2):

    式中為鹽分總量,kg/ hm2;S為取樣點(diǎn)含鹽量,g/kg;L為取樣孔控制土層寬度,m;h為取樣點(diǎn)控制土層厚度,m;γ為相應(yīng)土層土壤容重,g/cm3;為計(jì)算深度內(nèi)取樣層數(shù),0~100、0~60、0~40 cm取樣層數(shù)分別為8、6、4;為整個(gè)剖面、膜內(nèi)、膜外土壤取樣孔數(shù),分別為4、2、2。

    根據(jù)質(zhì)量守恒定律,播前收后土層鹽分總量變化的方程為公式(3):

    式中D為播前收后土層鹽分的改變量,kg/hm2,D>0說(shuō)明土層鹽分增加;△≤0說(shuō)明生育期土層鹽分不增加。2為收獲后土層鹽分總量,kg/hm2;1為播種前土層鹽分總量,kg/hm2。

    秋澆前后土壤鹽分的變化直接用土壤電導(dǎo)率的均值來(lái)評(píng)價(jià)。

    1.5 統(tǒng)計(jì)方法

    利用Excel數(shù)據(jù)分析功能進(jìn)行成對(duì)檢驗(yàn),采用SPSS13.0軟件進(jìn)行方差分析。

    2 結(jié)果與分析

    2.1 滴灌制度對(duì)土壤水分鹽分分布的影響

    抽雄期玉米需水量大[13],根系活動(dòng)范圍最大[14],吸水旺盛,頻繁的灌水會(huì)使滴灌濕潤(rùn)體的影響半徑達(dá)到生育期最大值,對(duì)此時(shí)期灌前土壤剖面水分狀況的分析有利于進(jìn)一步明確土壤水分的散失過(guò)程及下一時(shí)刻灌水后土壤水分的運(yùn)動(dòng)趨勢(shì)。圖5為玉米抽雄期各處理達(dá)到灌水下限時(shí)的土壤含水率剖面分布圖,圖中“白色短劃線”是土壤水分等值線從低處向高處凸出部分的連線,即集水線,眾多的集水線構(gòu)成了如河網(wǎng)狀的水分路線,筆者定義為土壤水分集水線。在微觀尺度下土壤水分集水線是所有水分的低值點(diǎn),理論上也是土水勢(shì)最低的點(diǎn),周?chē)耐寥浪钊菀紫虼诉\(yùn)動(dòng),所以其可以表示土壤水分運(yùn)動(dòng)的趨勢(shì)。由圖5可知,D1、D2、D3、D4土壤剖面質(zhì)量含水率均值分別為25.16%、22.72%、22.36%、21.44%,越靠近地表遠(yuǎn)離滴頭,水分含量越低。D1、D2、D3、D4土壤水分垂直方向急劇變化的土層深度范圍(白色框線內(nèi)部分)分別是25~40、33~48、37~52、30~ 52 cm,D1、D2、D3深度隨基質(zhì)勢(shì)減小而增加,范圍大小基本一致,而D4范圍增大??拷乇硗翆蛹窗咨蛞陨喜糠植煌幚砟?nèi)(距滴頭距離0~30 cm)土壤水分以水平變化為主,膜邊緣及膜外(黑框內(nèi)部分)土壤水分由下層向上層、由膜內(nèi)向膜外依次遞減。由土壤水分集水線可知,D1、D2以剖面水分最大值為界上層土壤水分有由深層向地表,由膜內(nèi)向膜外運(yùn)動(dòng)的趨勢(shì),下層土壤水分將向下運(yùn)動(dòng)。D3、D4整個(gè)剖面土壤水分由下層向上層依次遞減,灌水前土壤水分在蒸發(fā)蒸騰作用下經(jīng)由地表散失。

    玉米生育期末即秋澆前對(duì)應(yīng)的土壤鹽分剖面分布狀況是整個(gè)生育期內(nèi)土壤鹽分隨水分遷移運(yùn)動(dòng)積累的結(jié)果如圖6。D1土壤鹽分在距滴頭0~20 cm正下方范圍內(nèi)鹽分基本低于0.3 dS/m,在距滴頭水平距離20~60 cm范圍內(nèi),鹽分自深層至表層依次增加,在地表0~20 cm從膜內(nèi)至膜外鹽分逐漸積累。D2土壤鹽分剖面分布與D1類似,但膜內(nèi)滴頭下方鹽分低值區(qū)和膜外鹽分高值區(qū)范圍相對(duì)較小。D3在距滴頭0~40 cm正下方0~40 cm范圍內(nèi)形成一個(gè)鹽分低于0.3 dS/m的環(huán)狀范圍,在深度40~80 cm及0~20 cm膜外位置處鹽分積累。D4在距滴頭水平距離0~30 cm地表下50~80 cm及距滴頭20~60 cm地表下0~30 cm位置處鹽分明顯高于其他位置。將圖5黑框平移至圖6,此范圍內(nèi)土壤鹽分變化與土壤水分變化規(guī)律近乎相反,即膜下滴灌土壤鹽分在地表淺層膜邊緣及膜外范圍內(nèi)由下層向上層、由膜內(nèi)向膜外依次遞增,這可能是滴灌膜內(nèi)及深層土壤水分經(jīng)膜邊緣及膜外散失,鹽分遷移運(yùn)動(dòng)趨于地表積累的結(jié)果。

    2.2 滴灌制度對(duì)土壤溫度的影響

    不同灌溉制度下土壤溫度剖面變化總體規(guī)律基本是一致的,以灌水下限為-30 kPa即D3處理的土壤溫度剖面分布為例(圖7),膜內(nèi)土壤溫度受灌溉的影響較大,灌水后越接近地表土壤溫度下降幅度越大,地表5 cm土壤溫度需2~3 d后恢復(fù)至灌前水平,而膜外基本不受灌溉的影響。受溫度較低的灌溉水的影響,膜內(nèi)土壤溫度總體較膜外低0.08~0.45 ℃,隨著時(shí)間的推移,膜內(nèi)外溫度差異越來(lái)越小。從時(shí)間上來(lái)看,各個(gè)剖面土壤溫度變化規(guī)律是相似的,呈持續(xù)下降的趨勢(shì),這主要是氣溫作用的結(jié)果[30]。

    注:D1~D4表示膜下滴灌灌水下限分別為10、-20、-30、-40 kPa,下同。白色短劃線為土壤集水線,黑色框是土壤水分由下層向上層、由膜內(nèi)向膜外依次遞減的區(qū)域;白色框是土壤水分垂直方向急劇變化的區(qū)域;圖e為各處理水平方向距滴頭不同距離位置處土壤水分含量的均值。

    注:黑色框是土壤鹽分由下層向上層、由膜內(nèi)向膜外依次遞減的區(qū)域。圖e為各處理水平方向距滴頭不同距離位置處土壤電導(dǎo)率的均值,S0是土壤初始電導(dǎo)率。

    圖7 膜下滴灌土壤溫度剖面分布特征(D3)

    不同滴灌制度下土壤溫度的差異以膜內(nèi)5 cm土壤溫度為例進(jìn)行分析。如圖8a所示,6月上旬至中旬玉米苗期地表裸露,膜內(nèi)土壤溫度在氣溫作用下持續(xù)升至最高。6月中旬以后,土壤溫度整體持續(xù)下降一方面是灌溉水水溫低(8~12 ℃)導(dǎo)致土壤溫度降低,另一方面是玉米葉面積指數(shù)(leaf area index,LAI)增大,逐漸覆蓋地面,地表陽(yáng)光輻射作用減弱[30]。土壤溫度15日滑動(dòng)平均結(jié)果表明:7月下旬前,D3處理表層5 cm土壤溫度最高,D4次之,D2土壤平均溫度在D3和D1之間波動(dòng),D1土壤溫度最低(見(jiàn)圖8b)。這主要是由于6月上旬—7月下旬間D1、D2較D3灌水次數(shù)分別多了4次、2次,頻繁的灌溉導(dǎo)致土壤溫度降低,雖D3、D4灌水次數(shù)相同,但D3土壤平均含水率高于D4,土壤比熱容大有利于熱量的吸收和保持。因?yàn)?月上旬—8月上旬氣溫達(dá)到同期最高值,土壤較高的含水率有利于土壤白天吸熱、晚上保持土壤溫度,所以7月下旬—8月上旬土壤平均溫度D2最高,同時(shí)D4>D3>D1主要是受灌水影響的結(jié)果。8月中旬至9月下旬氣溫持續(xù)降低,土壤熱量散失,各處理土壤溫度趨于一致。因此,滴灌不同灌溉制度下土壤溫度的變化是氣溫、玉米LAI、灌水及土壤含水率共同作用的結(jié)果。劉洋等[31]研究發(fā)現(xiàn)玉米田有效積溫與玉米LAI關(guān)系密切,呈開(kāi)口向下的2次拋物線關(guān)系,在玉米營(yíng)養(yǎng)生長(zhǎng)階段,有效積溫對(duì)促進(jìn)玉米生長(zhǎng)作用顯著。河套灌區(qū)屬于典型的寒旱灌區(qū),玉米生育期前期面臨土壤積溫不足的問(wèn)題,所以,單就土壤溫度而言,河套灌區(qū)膜下滴灌玉米營(yíng)養(yǎng)生長(zhǎng)階段(即7月下旬前)應(yīng)控制灌水下限為-30 kPa。

    2.3 滴灌制度-秋澆對(duì)土壤鹽分平衡的影響

    膜下滴灌不同灌溉制度下0~100、0~60、0~40 cm土層總鹽分、膜內(nèi)、膜外生育期鹽分改變量,如圖9所示。由圖9可知,土層總鹽分及膜外鹽分生育期內(nèi)0~100、0~60、0~40 cm均是增加的。D1、D2、D3、D4處理總鹽量分別增加2 011.57、3 751.42、5 377.69、4 427.89 kg/hm2,鹽分增加總量分別占到了灌溉水引入土壤鹽分總量(D1、D2、D3、D4引入的鹽量分別為5 542.6、4 758.3、3 940.8、2 870.8 kg/hm2)的36.29%、78.84%、136.46%、154.24%,說(shuō)明控制膜下滴灌灌水下限為-10和-20 kPa,可有效淋濾0~100 cm土層的鹽分,而D3、D4土層鹽分增加量高于灌溉水引入土壤鹽分的原因可能是河套灌區(qū)5月末6月初在強(qiáng)烈的蒸發(fā)作用下深層鹽分向上運(yùn)動(dòng)[32-33],而灌溉水量較少未能有效淋濾的結(jié)果。D1、D2、D3、D4各處理0~100 cm膜外土壤鹽分增加量占鹽分總增加量的百分比分別為97.25%、78.30%、70.69%、57.34%,說(shuō)明在滴灌局部濕潤(rùn)和強(qiáng)烈的蒸發(fā)作用下,膜外土壤成為鹽分的“承泄區(qū)”,土壤基質(zhì)勢(shì)越高(D1最高),灌水越頻繁,灌溉水向膜外濕潤(rùn)范圍越大,越多的鹽分遷移到膜外。D1、D2、D3、D4處理0~60 cm土層鹽分分別增加1 718.18、2 461.88、2 339.47、2 909.59 kg/hm2,分別占0~100 cm土層鹽分總增加量的85.42%、65.63%、43.50%、65.71%,0~60 cm積鹽率最低的為D3。0~60 cm土層膜內(nèi)和膜外土壤鹽分都較播前增加,但膜外增加更加明顯,其鹽量增加量占總增加量的106.38%、78.50%、71.85%、81.19%。0~40 cm各處理土壤總鹽分分別增加821.76、1 095.48、489.85、1 898.88 kg/hm2,分別占0~100 cm土層鹽分總增加量的47.83%、44.50%、20.94%、65.27%,D3處理0~40 cm積鹽率最低,膜外增加量占0~40 cm總增加量的145.84%、93.46%、144.12%、102.53%。各處理0~40 cm膜內(nèi)土壤沒(méi)有明顯積鹽。不同灌溉制度下各土層鹽分平衡及膜內(nèi)膜外鹽分分配結(jié)果表明,膜下滴灌土壤鹽分水平方向由膜內(nèi)向膜外地表裸露區(qū)遷移,趨于膜外積累。

    不同滴灌灌水下限及秋澆對(duì)整個(gè)土壤剖面及膜內(nèi)、膜外0~100 cm土層鹽分的影響見(jiàn)表2。秋澆前D1處理0~100 cm土壤鹽分不論整個(gè)剖面還是膜外均顯著小于其余3個(gè)處理(<0.05),說(shuō)明控制灌水下限為-10 kPa可有效淋濾0~100 cm 土層鹽分。秋澆前0~100 cm膜內(nèi)土壤鹽分D1、D2明顯低于D3、D4(<0.05),且D1、D2處理差異不顯著(>0.05),說(shuō)明生育期內(nèi)控制灌水下限為-10、-20 kPa可有效地淋濾0~100 cm膜內(nèi)土壤鹽分。但總體而言,生育期內(nèi)D1處理對(duì)鹽分淋洗效果最明顯,而D2、D3、D4處理對(duì)0~100 cm土層鹽分的影響差異性短期內(nèi)不明顯。秋澆后各處理間土壤電導(dǎo)率差異不顯著(>0.05),且土壤電導(dǎo)率變異系數(shù)較秋澆前均明顯降低,表明秋澆后鹽分分布更加均勻,這與土壤含鹽量越高淋濾效果越好有關(guān)[34],各土層鹽分均降至非鹽漬化水平[35]。秋澆前后D1、D2、D3、D4處理0~100 cm的土壤電導(dǎo)率均值經(jīng)成對(duì)檢驗(yàn),值均小于0.05,均值差異性顯著,秋澆后較秋澆前分別降低22.77%、10.86%、26.14%、12.59%。可見(jiàn),D3處理下滴灌-秋澆洗鹽效果最佳。然而,河套灌區(qū)地下水埋深淺,凍融影響使鹽分的問(wèn)題更加復(fù)雜,如何合理有效地淋洗鹽分使次年不影響苗期作物且使長(zhǎng)年滴灌條件下土壤不至發(fā)生鹽漬化是河套灌區(qū)發(fā)展膜下滴灌面臨的關(guān)鍵問(wèn)題。因此,河套灌區(qū)膜下滴灌鹽分累積到何種程度秋澆及具體合理的秋澆制度有重要的研究?jī)r(jià)值,還需進(jìn)一步深入研究。

    表2 滴灌灌水下限及秋澆對(duì)0~100 cm土層土壤鹽分的影響

    2.4 滴灌制度對(duì)玉米產(chǎn)量品質(zhì)的影響

    不同滴灌制度對(duì)玉米百粒質(zhì)量、禿尖、產(chǎn)量、品質(zhì)的影響如表3所示。D2處理禿尖最小,D2和D3處理百粒質(zhì)量最大,產(chǎn)量最高(分別為16 511.0和14 964.4 kg/hm2),與D1、D4產(chǎn)量差異性顯著,分別高16.03%、19.61%(<0.05)。由于品質(zhì)檢測(cè)使用各處理重復(fù)的混合樣,故未做顯著性分析??傮w來(lái)看,D1、D2、D3、D4處理粗脂肪、粗淀粉、粗蛋白的總量分別占總產(chǎn)量的81.79%、75.82%、80.51%、88.69%,為11 638.7、12 518.3、12 047.4、12 243.3 kg/hm2,差別不大。

    表3 不同滴灌制度下玉米產(chǎn)量品質(zhì)

    3 討 論

    秋澆洗鹽是為了滿足次年作物生長(zhǎng)的要求,因此,次年播種前后的土壤鹽分狀況是評(píng)價(jià)秋澆灌水制度優(yōu)劣的主要依據(jù)[24]。本研究以次年春播前土壤鹽分含量與秋澆前鹽分含量的差值來(lái)近似評(píng)估秋澆洗鹽的效果。秋澆鹽分的淋洗效果一是受秋澆定額的影響,秋澆定額越大,洗鹽效果越好[24];二是受秋澆前地下水位埋深的影響,埋深愈大,秋澆愈能將鹽分淋洗至土壤深層,淋洗效果也愈好[36]。另外,凍融作用對(duì)秋澆土壤鹽分淋洗過(guò)程也有復(fù)雜的影響。在凍結(jié)期,由于凍結(jié)是一個(gè)自上而下逐漸進(jìn)行的過(guò)程,其間非飽和土壤的內(nèi)排水過(guò)程,會(huì)使得鹽分隨水分向下遷移,這在一定程度上抑制了凍結(jié)期地下水及其下層土壤中鹽分的上移。彭振陽(yáng)等[37]的研究表明,在土壤凍結(jié)期間,在凍結(jié)和排水的共同作用下凍結(jié)層儲(chǔ)鹽總量變化不大,凍結(jié)層返鹽率與地下水的排水排鹽條件是密切相關(guān)的,排水排鹽條件越好,田間儲(chǔ)鹽量減少越多,會(huì)使得凍結(jié)層返鹽率越低,甚至?xí)霈F(xiàn)脫鹽現(xiàn)象。在消融期,一般認(rèn)為在凍結(jié)層以下部分土壤水轉(zhuǎn)化為地下水,土壤處于排水過(guò)程,鹽分也會(huì)隨之淋洗,凍結(jié)層以上由于蒸發(fā)作用,表層會(huì)積鹽。因此,秋澆前至次年春播前土壤鹽分的變化是秋澆和凍融共同作用的結(jié)果。秋澆前至次年春播前,土壤水鹽的變化及運(yùn)動(dòng)規(guī)律受到多方面因素的影響,是一個(gè)復(fù)雜的過(guò)程,還需進(jìn)一步深入研究。

    4 結(jié) 論

    通過(guò)對(duì)河套灌區(qū)膜下滴灌不同灌溉制度生育期土壤水熱鹽效應(yīng)及非生育期洗鹽灌溉研究,有如下結(jié)論:

    1)膜下滴灌鹽分隨水分定向遷移運(yùn)動(dòng),不同滴灌制度下土壤鹽分剖面分布極不均勻,各處理鹽分均由膜內(nèi)向膜外地表裸露區(qū)定向遷移,趨于膜外地表積累??刂乒嗨孪逓?10 kPa(D1處理)可有效淋濾0~100 cm土壤鹽分,而其他處理對(duì)0~100 cm土層鹽分的影響差異性短期內(nèi)不明顯,需要對(duì)不同處理長(zhǎng)期膜下滴灌的鹽分進(jìn)一步觀測(cè)研究。

    2)膜下滴灌土壤溫度受氣溫、玉米葉面積指數(shù)、灌水及土壤含水率共同作用,膜內(nèi)表層土壤溫度在玉米需水旺季劇烈波動(dòng),灌水后會(huì)劇烈下降,2~3 d后恢復(fù)。就土壤積溫而言,建議河套灌區(qū)膜下滴灌玉米抽雄期以前(即7月下旬前)控制灌水下限為-30 kPa。

    3)非生育期洗鹽灌溉(秋澆)效果顯著,秋澆灌黃河水180 mm后,次年春播前0~100 cm土壤鹽分下降10.86%~26.14%,剖面分布較均勻,是干旱半干旱地區(qū)控制膜下滴灌土壤鹽分的有效途徑。

    4)D2和D3處理百粒質(zhì)量最大,產(chǎn)量最高(分別為16 511.0和14 964.4 kg/hm2),與D1、D4產(chǎn)量差異性顯著(<0.05)。各處理粗脂肪、粗淀粉、粗蛋白含量差別不大。

    5)河套灌區(qū)膜下滴灌土壤鹽分調(diào)控建議為生育期滴灌灌溉制度和非生育期洗鹽灌溉雙重調(diào)控,在一年一秋澆的條件下,玉米生育期膜下滴灌灌水下限宜以作物高產(chǎn)為目標(biāo),玉米灌水下限可為-30 kPa。在多年秋澆的條件下,膜下滴灌鹽分累積到何種程度洗鹽灌溉及具體合理的膜下滴灌-秋澆洗鹽灌溉制度還需進(jìn)一步深入研究。

    致謝:感謝武漢大學(xué)楊金忠教授、查元源副研究員及內(nèi)蒙古農(nóng)業(yè)大學(xué)李仙岳教授對(duì)論文初稿所提的寶貴修改意見(jiàn)。

    [1] 屈忠義,楊曉,黃永江. 內(nèi)蒙古河套灌區(qū)節(jié)水工程改造效果分析評(píng)估[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2015,46(4):70-76. Qu Zhongyi, Yang Xiao, Huang Yongjiang. Analysis and assessment of water-saving project of Hetao Irrigation District in Inner Mongolia[J]. Transactions of the Chinese Society for Agricultural Machinery,2015,46(4): 70-76. (in Chinese with English abstract)

    [2] 張志杰,楊樹(shù)青,史海濱,等. 內(nèi)蒙古河套灌區(qū)灌溉入滲對(duì)地下水的補(bǔ)給規(guī)律及補(bǔ)給系數(shù)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2011,27(3):61-66. Zhang Zhijie, Yang Shuqing, Shi Haibin, et al. Irrigation infiltration and recharge coefficient in Hetao irrigation district in Inner Mongolia[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(3): 61-66. (in Chinese with English abstract)

    [3] Wang Ruoshui, Kang Yaohu, Wan Shuqin, et al. Salt distribution and the growth of cotton under different drip irrigation regimes in a saline area[J]. Agricultural Water Management, 2011, 100(1): 58-69.

    [4] Kang Yaohu, Wang Ruoshui, Wan Shuqin, et al. Effects of different water levels on cotton growth and water use through drip irrigation in an arid region with saline ground water of Northwest China[J]. Agricultural Water Management, 2012, 109(109): 117-126.

    [5] Wang Ruoshui, Kang Yaohu, Wan Shuqin, et al. Influence of different amounts of irrigation water on salt leaching and cotton growth under drip irrigation in an arid and saline area[J]. Agricultural Water Management, 2012, 110(7):109-117.

    [6] Mmolawa K, Or D. Root zone solute dynamics under drip irrigation: A review[J]. Plant and Soil, 2000, 222(1): 163-190.

    [7] Shannon M C, Grieve C M. Tolerance of vegetable crops to salinity[J]. Scientia Horticulturae, 1998, 78(1/2/3/4): 5-38.

    [8] Darwish T, Atallah T, Moujabber M E, et al. Salinity evolution and crop response to secondary soil salinity in two agro-climatic zones in Lebanon[J]. Agricultural Water Management, 2005, 78(1/2): 152-164.

    [9] 葉志勇,郭克貞,趙淑銀,等. 河套灌區(qū)加工型番茄膜下滴灌技術(shù)研究[J]. 灌溉排水學(xué)報(bào),2011,30(1):110-112. Ye Zhiyong, Guo Kezhen, Zhao Shuyin, et al. Mulched drip irrigation of processed tomato in Hetao Irrigation District[J]. Journal of Irrigation and Drainage, 2011, 30(1): 110-112. (in Chinese with English abstract)

    [10] 郭夢(mèng)吉,劉宇,任樹(shù)梅,等. 河套灌區(qū)微咸淡水交替灌溉對(duì)加工番茄根系生長(zhǎng)的影響[J]. 中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào),2016,21(2):65-72.Guo Mengji, Liu Yu, Ren Shumei, et al. Effect of alternative irrigation between brackish water and fresh water on root growth of processing tomoto in Hetao Irrigation District[J]. Journal of China Agricultural University, 2016, 21(2): 65-72. (in Chinese with English abstract)

    [11] 范雅君,呂志遠(yuǎn),田德龍,等. 河套灌區(qū)玉米膜下滴灌灌溉制度研究[J]. 干旱地區(qū)農(nóng)業(yè)研究,2015,33(1):123-129. Fan Yajun,Lü Zhiyuan,Tian Delong, et al. Irrigation regime for corn production with drip irrigation andplastic mulching in Hetao Irrigation District[J]. Agricultural Research in the Arid Areas, 2015, 33(1): 123-129. (in Chinese with English abstract)

    [12] Liu Haijun, Wang Xuming, Zhang Xian, et al. Evaluation on the responses of maize () growth, yield and water use efficiency to drip irrigation water under mulch condition in the Hetao irrigation District of China[J]. Agricultural Water Management, 2017, 179(1): 144-157.

    [13] 杜斌,屈忠義,于健,等. 河套灌區(qū)大田作物膜下滴灌作物系數(shù)試驗(yàn)研究[J]. 灌溉排水學(xué)報(bào),2014,33(4):16-20.Du Bin, Qu Zhongyi, Yu Jian, et al. Experimental study on crop coefficient under mulched drip irrigation in Hetao Irrigation District of Inner Mongolia[J]. Journal of Irrigation and Drainage, 2014, 33(4): 16-20. (in Chinese with English abstract)

    [14] 李仙岳,史海濱,龔雪文,等. 立體種植農(nóng)田不同生育期及土壤水分的根系分布特征[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2014,45(3):140-147. Li Xianyue, Shi Haibin, Gong Xuewen, et al. Root distribution in strip intercropping field under different growth period and different soil water[J]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(3): 140-147. (in Chinese with English abstract)

    [15] 史海濱,赤江剛夫,長(zhǎng)堀金造,等. 鹽分再分配條件下河套灌區(qū)沖洗水量的數(shù)值模擬[J]. 農(nóng)業(yè)工程學(xué)報(bào),2002,18(5):67-72. Shi Hai Bin, Akae Takeo, Nagahori Kinzo, et al. Simulation of leaching requirement for Hetao Irrigation District considering salt redistribution after irrigation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2002, 18(5): 67-72. (in English with Chinese abstract)

    [16] Noborio K, Mclnnes K J, Heilman J L. Two-dimensional model for water, heat, and solute transport in furrow-irrigated soil: I, II[J]. Soil Science Society of America Journal, 1996, 60(4): 1001-1021.

    [17] 周和平,王少麗,吳旭春. 膜下滴灌微區(qū)環(huán)境對(duì)土壤水鹽運(yùn)移的影響[J]. 水科學(xué)進(jìn)展,2014,25(6):816-824. Zhou Heping, Wang Shaoli, Wu Xuchun. Micro drip irrigation district environmental impact on soil water and salt transport[J]. Advances in Water Science, 2014, 25(6): 816-824. (in Chinese with English abstract)

    [18] 周和平,王少麗,姚新華,等. 膜下滴灌土壤水鹽定向遷移分布特征及排鹽效應(yīng)研究[J]. 水利學(xué)報(bào),2013,44(11):1380-1388. Zhou Heping, Wang Shaoli, Yao Xinhua, et al. Research on distribution characteristics and salt-discharging effect of directional migration of water and salt in soil through drip irrigation under plastic film[J]. Journal of Hydraulic Engineering, 2013, 44(11): 1380-1388. (in Chinese with English abstract)

    [19] Liu Meixian, Yang Jingsong, Li Xiaoming, et al. Distribution and dynamics of soil water and salt under different drip irrigation regimes in northwest China[J]. Irrigation Science, 2013, 31(4): 675-688.

    [20] Phocaides A. Water quality for irrigation[M]//Handbook on Pressurized Irrigation Techniques. 2nd ed. Rome, Italy: Food and Agriculture Organization of the United Nations, 2007.

    [21] 李金剛,屈忠義,黃永平,等. 微咸水膜下滴灌不同灌水下限對(duì)鹽堿地土壤水鹽運(yùn)移及玉米產(chǎn)量的影響[J]. 水土保持學(xué)報(bào),2017,31(1):217-223.Li Jingang, Qu Zhongyi, Huang Yongping, et al. Effects of control lower limit of saline water mulched drip irrigation on water salt movement and corn yield in saline soil[J]. Journal of Soil and Water Conservation, 2017, 31(1): 217-223.

    [22] Ren Xiaodong, Qu Zhongyi, Martins D S, et al. Daily reference evapotranspiration for hyper-arid to moist sub-humid climates in Inner Mongolia, China: I. Assessing temperature methods and spatial variability[J]. Water Resources Management, 2016, 30(11): 3769-3791.

    [23] Kang Yaohu, Chen Ming, Wan Shuqin. Effects of drip irrigation with saline water on waxy maize (var.ceratina Kulesh) in North China Plain[J]. Agricultural Water Management, 2010, 97(9):1303-1309

    [24] 李瑞平,史海濱,赤江剛夫,等. 基于SHAW 模型的內(nèi)蒙古河套灌區(qū)秋澆節(jié)水灌溉制度[J]. 農(nóng)業(yè)工程學(xué)報(bào),2010,26(2):31-36. Li Ruiping, Shi Haibin, Takeo Akae, et al. Scheme of water saving irrigation in autumn based on SHAW model in Inner Mongolia Hetao irrigation district[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(2): 31-36. (in Chinese with English abstract)

    [25] 孫貫芳,屈忠義,杜斌,等. 內(nèi)蒙古河套灌區(qū)不同灌溉模式對(duì)土壤溫度及鹽分的影響[J]. 節(jié)水灌溉,2016(2):28-31.Sun Guanfang, Qu Zhongyi, Du Bin, et al. Effects of different irrigation modes on soil temperature and salinity in Inner Mongolia Hetao Irrigation District[J]. Water Saving Irrigation, 2016(2): 28-31. (in Chinese with English abstract)

    [26] 中國(guó)國(guó)家標(biāo)準(zhǔn)化管理委員會(huì).飼料中粗脂肪的測(cè)定:GB/T 6433-2006[S]. 北京:中國(guó)標(biāo)準(zhǔn)出版社,2006.

    [27] 國(guó)家技術(shù)監(jiān)督局.飼料中粗蛋白測(cè)定方法:GB/T 6432-1994[S]. 北京:中國(guó)標(biāo)準(zhǔn)出版社,1994.

    [28] 國(guó)家標(biāo)準(zhǔn)局.谷物籽粒粗淀粉測(cè)定法:NY/T 11-1985[S]. 北京:中國(guó)標(biāo)準(zhǔn)出版社,1985.

    [29] 栗現(xiàn)文,靳孟貴,袁晶晶,等. 微咸水膜下滴灌棉田漫灌 洗鹽評(píng)價(jià)[J]. 水利學(xué)報(bào),2014,45(9):1091-1098. Li Xianwen, Jin Menggui, Yuan Jingjing, et al. Evaluation of soil salts leaching in cotton field after mulched drip irrigation with brackish water by freshwater flooding[J]. Journal of Hydraulic Engineering, 2014, 45(9): 1091-1098. (in Chinese with English abstract)

    [30] 孫貫芳,杜斌,屈忠義,等. 不同灌溉模式下土壤溫度的變化及對(duì)氣溫的響應(yīng)特征[J]. 土壤,2016,48(3):581-587. Sun Guanfang, Du Bin, Qu Zhongyi, et al. Changes of soil temperature and response to air temperature under different irrigation modes[J]. Soils, 2016, 48(3): 581-587. (in Chinese with English abstract)

    [31] 劉洋,栗巖峰,李久生,等. 東北半濕潤(rùn)區(qū)膜下滴灌對(duì)農(nóng)田水熱和玉米產(chǎn)量的影響[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2015,46(10):93-104. Liu Yang, Li Yanfeng, Li Jiusheng, et al. Effects of mulched drip irrigation on water and heat conditions in field and maize yield in sub-humid region of northeast China[J]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(10): 93-104. (in Chinese with English abstract)

    [32] 杜社妮,于健,耿桂俊,等. 定植孔密封方式對(duì)土壤水熱鹽及番茄苗存活率的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2012,28(12):110-116.Du Sheni, Yu Jian, Geng Guijun, et al. Effects of planting hole sealing methods on soil moisture, temperature, salinity and survival rate of tomato seedling[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(12): 110-116. (in Chinese with English abstract)

    [33] Ayars J E,Christen E W,Hornbuckle J W. Controlled drainage for improved water management in arid zone irrigated agriculture[J]. Agricultural Water Management, 2006, 86(1): 128-139.

    [34] 龐鴻斌. 賈大林論文選[M]. 北京:中國(guó)農(nóng)業(yè)科學(xué)技術(shù)出版社,2005:3-11.

    [35] 史海濱,李瑞平,楊樹(shù)青. 鹽漬化土壤水熱鹽遷移與節(jié)水灌溉理論研究[M]. 北京:中國(guó)水利水電出版社,2011:12-13.

    [36] 陳亞新,史海濱,田存旺.地下水與土壤鹽漬化關(guān)系的動(dòng)態(tài)模擬[J]. 水利學(xué)報(bào),1997(5):77-83. Chen Yaxin, Shi Haibin, Tian Cunwang. Dynamic simulation of the relationship between water table and salinization of soil[J]. Journal of Hydraulic Engineering, 1997(5): 77-83.

    [37] 彭振陽(yáng),黃介生,伍靖?jìng)?,? 秋澆條件下季節(jié)性凍融土壤鹽分運(yùn)動(dòng)規(guī)律[J]. 農(nóng)業(yè)工程學(xué)報(bào),2012,28(6):77-81. Peng Zhenyang, Huang Jiesheng, Wu Jingwei, et al. Salt movement of seasonal freezing-thawing soil under autumn irrigation condition[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(6): 77-81. (in Chinese with English abstract)

    Water-heat-salt effects of mulched drip irrigation maize with different irrigation scheduling in Hetao Irrigation District

    Sun Guanfang1,2, Qu Zhongyi1※, Du Bin3, Ren Zhongsheng1, Liu Anqi1

    (1.010018,; 2.430072,; 3.750000,)

    Along with the development of Northwest Water-saving and Benefit-increasing strategy in China, mulched drip irrigation is booming in Hetao Irrigation District because of its advantage in water saving and high production. However, soil salinization has been a potential risk factor in arid area, thus enough attention should be drawn in the management of soil salt. This research studied the influence of different mulched drip irrigation schedules (lower irrigation limits were-10,-20,-30 and-40 kPa, respectively) on water and salt profile distribution characteristics, soil temperature, maize production and quality. Moreover, the salt-leaching irrigation in non-crop growth period (i.e., autumn irrigation, which is a traditional irrigation custom in Hetao Irrigation District after crop was harvested in autumn to leach salt and store the water in soil for the next year) had also been investigated. The experiment was carried out in the year of 2014-2015. The irrigation quota for each treatment was 518.0, 444.7, 368.3 and 268.3 mm, respectively. The leaching irrigation amount was 180 mm. The soil samples were collected for moisture, salinity measurements. Soil temperature was also sequentially measured during the experiment. After harvest, the maize yield and quality were determined. The results showed that the distributions of soil water and salt profiles were uneven under different mulched drip irrigation schedules. Soil salt migrated from the central mulching soil to the bare field, tending to accumulate in surface soil layer outside the film during the mulched drip irrigation. Soil temperature was influenced by air temperature, leaf area index, irrigation practice and soil moisture comprehensively. Soil temperature dropped sharply after irrigation practice, and recovered 2-3 days later. The optimum accumulated soil temperature for the nutritive growth period of maize was-30 kPa. Salt accumulation at the depth of 0-40 cm inside the film in maize growth period was little for almost all the treatments. The treatment of-10 kPa could leach 0-100 cm soil salinity effectively, while influence of the other treatments on the 0-100 cm soil salt was not obvious in the short term, and long-term observation to study impact of different drip irrigation lower limit on soil salt is needed. The effect of salt-leaching irrigation with 180 mm Yellow River water in non-crop growth period was remarkable, and the soil salinity decreased by 10.86%-26.14% on average, resulting in a relatively uniform salt profile at 0-100 cm. A combination of irrigation schedule for water-saving and leaching irrigation in the non-crop growth period for salt-controlling was helpful to control the salt accumulation. The maize yield was highest in the treatment of-20 and-30 kPa. The yield was 16 511 and 14 964.4 kg/hm2, respectively. The quality of maize was similar for all the treatments. Considering soil moisture, temperature, salinity and maize yield, we suggested that the limit of soil suction to trigger irrigation was-30 kPa. However, due to the impacts of the frost-freeze phenomenon and particular hydrogeological condition in Hetao Irrigation District, the detailed leaching irrigation schedule to control soil salt induced by mulched drip irrigation need further investigation in the future.

    soil moisture; temperature; salinity; Hetao Irrigation District; mulched drip irrigation; water-heat-salt; salt leaching; autumn irrigation

    10.11975/j.issn.1002-6819.2017.12.019

    S278

    A

    1002-6819(2017)-12-0144-9

    2016-10-25

    2017-04-10

    國(guó)家自然科學(xué)基金(51469020);內(nèi)蒙古自治區(qū)水利科技重大專項(xiàng)([2014] 117)。

    孫貫芳,男,河南三門(mén)峽人,博士生,主要研究方向?yàn)楣?jié)水灌溉原理與新技術(shù)。呼和浩特 內(nèi)蒙古農(nóng)業(yè)大學(xué)水利與土木建筑工程學(xué)院,010018。Email:gfsun1990@126.com

    屈忠義,男,內(nèi)蒙古巴彥淖爾市人,教授,博士生導(dǎo)師,主要從事節(jié)水灌溉原理與技術(shù)方面的研究。呼和浩特 內(nèi)蒙古農(nóng)業(yè)大學(xué)水利與土木建筑工程學(xué)院,010018。Email:quzhongyi@imau.edu.cn

    孫貫芳,屈忠義,杜 斌,任中生,劉安琪. 不同灌溉制度下河套灌區(qū)玉米膜下滴灌水熱鹽運(yùn)移規(guī)律[J]. 農(nóng)業(yè)工程學(xué)報(bào),2017,33(12):144-152. doi:10.11975/j.issn.1002-6819.2017.12.019 http://www.tcsae.org

    Sun Guanfang, Qu Zhongyi, Du Bin, Ren Zhongsheng, Liu Anqi. Water-heat-salt effects of mulched drip irrigation maize with different irrigation scheduling in Hetao Irrigation District[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(12): 144-152. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2017.12.019 http://www.tcsae.org

    猜你喜歡
    河套土壤溫度鹽分
    河套地區(qū)南美白對(duì)蝦養(yǎng)殖試驗(yàn)
    近代晉陜移民與“河套”內(nèi)涵變遷
    遼東山區(qū)3種人工林土壤呼吸對(duì)土壤溫度和土壤水分的響應(yīng)
    管群間歇散熱的土壤溫度響應(yīng)與恢復(fù)特性
    長(zhǎng)期膜下滴灌棉田根系層鹽分累積效應(yīng)模擬
    攝影欣賞
    土壤濕度和土壤溫度模擬中的參數(shù)敏感性分析和優(yōu)化
    河套灌區(qū):興水利強(qiáng)基礎(chǔ) 惠民生促和諧
    河套灌區(qū)試用無(wú)人機(jī)獲取作物需水信息
    基于PLSR的陜北土壤鹽分高光譜反演
    人妻制服诱惑在线中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 国产91av在线免费观看| 久久久a久久爽久久v久久| 啦啦啦韩国在线观看视频| 久久久久国产网址| 久久精品夜夜夜夜夜久久蜜豆| 亚洲最大成人av| 久久久久网色| 成人二区视频| 女的被弄到高潮叫床怎么办| a级毛色黄片| 久久久久性生活片| 97热精品久久久久久| 一级毛片我不卡| 麻豆国产av国片精品| av女优亚洲男人天堂| 最近手机中文字幕大全| 亚洲国产色片| 插阴视频在线观看视频| 亚洲高清免费不卡视频| 有码 亚洲区| 久久久久久久久大av| 国内少妇人妻偷人精品xxx网站| 全区人妻精品视频| 少妇熟女aⅴ在线视频| 精品久久久久久久人妻蜜臀av| 在线国产一区二区在线| 中文字幕制服av| 国产精品久久久久久久久免| 久久午夜亚洲精品久久| 一个人免费在线观看电影| 一区二区三区四区激情视频 | 国产国拍精品亚洲av在线观看| 大型黄色视频在线免费观看| 日本免费一区二区三区高清不卡| 寂寞人妻少妇视频99o| 国产成人午夜福利电影在线观看| 极品教师在线视频| 男人的好看免费观看在线视频| 麻豆国产97在线/欧美| 少妇的逼好多水| 中国国产av一级| 女人十人毛片免费观看3o分钟| 亚洲最大成人手机在线| 美女 人体艺术 gogo| 夜夜夜夜夜久久久久| 亚洲图色成人| 亚洲欧美日韩高清在线视频| 成人午夜精彩视频在线观看| 天天一区二区日本电影三级| av天堂中文字幕网| 午夜福利在线在线| 成人国产麻豆网| 日韩成人av中文字幕在线观看| 好男人视频免费观看在线| 亚洲av电影不卡..在线观看| 三级男女做爰猛烈吃奶摸视频| 亚洲精品日韩av片在线观看| 日本爱情动作片www.在线观看| 欧美xxxx性猛交bbbb| 美女被艹到高潮喷水动态| 成人毛片60女人毛片免费| 亚州av有码| 亚洲精品久久国产高清桃花| av黄色大香蕉| 色哟哟哟哟哟哟| 黑人高潮一二区| 国产成人影院久久av| 又爽又黄a免费视频| 变态另类成人亚洲欧美熟女| 色吧在线观看| 人妻少妇偷人精品九色| 中国美白少妇内射xxxbb| 国产极品天堂在线| av视频在线观看入口| 人妻久久中文字幕网| 看免费成人av毛片| 少妇被粗大猛烈的视频| 一级毛片我不卡| 日韩欧美国产在线观看| 男的添女的下面高潮视频| 成人特级av手机在线观看| 一进一出抽搐动态| 欧美日本视频| 久久精品国产自在天天线| 欧美zozozo另类| 久久精品国产99精品国产亚洲性色| 国产成人一区二区在线| 热99在线观看视频| 一区二区三区免费毛片| 啦啦啦观看免费观看视频高清| 日本黄大片高清| 国产一区二区亚洲精品在线观看| 秋霞在线观看毛片| 国产白丝娇喘喷水9色精品| 国产高清激情床上av| 欧美成人a在线观看| 在线播放国产精品三级| 精品人妻偷拍中文字幕| 午夜亚洲福利在线播放| 99久国产av精品| 国产毛片a区久久久久| 久久婷婷人人爽人人干人人爱| 一区二区三区免费毛片| 国产日韩欧美在线精品| 日韩av不卡免费在线播放| 听说在线观看完整版免费高清| 国产视频内射| 看十八女毛片水多多多| 亚洲熟妇中文字幕五十中出| 午夜爱爱视频在线播放| 国产精品乱码一区二三区的特点| 国产精品免费一区二区三区在线| 在线免费观看不下载黄p国产| 国产亚洲精品av在线| 少妇人妻一区二区三区视频| 一本久久精品| 日韩欧美在线乱码| 日韩欧美精品v在线| 岛国在线免费视频观看| 亚洲人成网站高清观看| 只有这里有精品99| 两个人视频免费观看高清| 蜜桃久久精品国产亚洲av| 高清午夜精品一区二区三区 | 久久这里有精品视频免费| 亚洲av中文av极速乱| 美女xxoo啪啪120秒动态图| 免费看日本二区| 欧美最新免费一区二区三区| 日韩强制内射视频| 成年av动漫网址| av.在线天堂| АⅤ资源中文在线天堂| 成人二区视频| 国产免费一级a男人的天堂| 中国国产av一级| 中文字幕免费在线视频6| 夜夜看夜夜爽夜夜摸| 99在线人妻在线中文字幕| 久久久久性生活片| 在线a可以看的网站| 国产成年人精品一区二区| 国产精品一区www在线观看| 国产成人一区二区在线| 亚洲18禁久久av| 亚洲人与动物交配视频| 国产白丝娇喘喷水9色精品| 国产一区亚洲一区在线观看| 久久精品国产99精品国产亚洲性色| 网址你懂的国产日韩在线| 97超视频在线观看视频| 啦啦啦啦在线视频资源| 小说图片视频综合网站| 免费大片18禁| 亚洲av成人精品一区久久| 国产精品久久久久久亚洲av鲁大| 国产精品久久视频播放| 中文字幕av在线有码专区| 深夜a级毛片| 成人无遮挡网站| 国产伦精品一区二区三区四那| 国产精品福利在线免费观看| 免费观看的影片在线观看| 国产激情偷乱视频一区二区| 国产不卡一卡二| 午夜免费激情av| 国产伦精品一区二区三区四那| 91av网一区二区| 在线观看午夜福利视频| 日日摸夜夜添夜夜添av毛片| 成人国产麻豆网| 亚洲国产精品成人久久小说 | 97热精品久久久久久| 国产伦精品一区二区三区四那| 亚洲成人中文字幕在线播放| 国产美女午夜福利| av福利片在线观看| 最近视频中文字幕2019在线8| 欧美极品一区二区三区四区| 男女啪啪激烈高潮av片| videossex国产| 美女被艹到高潮喷水动态| 五月伊人婷婷丁香| 少妇人妻精品综合一区二区 | 国产精品爽爽va在线观看网站| 精品99又大又爽又粗少妇毛片| 精品一区二区三区人妻视频| 女人十人毛片免费观看3o分钟| 啦啦啦观看免费观看视频高清| 亚洲精品自拍成人| 又粗又硬又长又爽又黄的视频 | 亚洲精品456在线播放app| 丝袜喷水一区| 久久精品综合一区二区三区| 成熟少妇高潮喷水视频| 精品久久久久久久久久免费视频| 婷婷亚洲欧美| 最新中文字幕久久久久| 日本一二三区视频观看| 精品国内亚洲2022精品成人| 日本黄色视频三级网站网址| 国产蜜桃级精品一区二区三区| 两个人视频免费观看高清| 欧美极品一区二区三区四区| 亚洲成人av在线免费| 精品一区二区免费观看| 日本av手机在线免费观看| 日韩 亚洲 欧美在线| 国产午夜精品一二区理论片| 日韩一本色道免费dvd| 黄色视频,在线免费观看| 午夜爱爱视频在线播放| 国产又黄又爽又无遮挡在线| 白带黄色成豆腐渣| 亚洲第一区二区三区不卡| 精品国产三级普通话版| 91久久精品国产一区二区成人| 观看免费一级毛片| eeuss影院久久| 老司机影院成人| 日韩国内少妇激情av| 又粗又硬又长又爽又黄的视频 | av天堂中文字幕网| 久久久精品欧美日韩精品| 亚洲第一电影网av| 国产极品天堂在线| 人妻少妇偷人精品九色| 亚洲电影在线观看av| 国产精品女同一区二区软件| 欧美最黄视频在线播放免费| 99热网站在线观看| 成人亚洲精品av一区二区| 国产高潮美女av| 99久久无色码亚洲精品果冻| 日韩欧美在线乱码| 亚洲国产欧洲综合997久久,| 久久精品国产亚洲av天美| av天堂在线播放| 插阴视频在线观看视频| 亚洲av二区三区四区| a级毛片免费高清观看在线播放| 成年av动漫网址| 波多野结衣高清无吗| 中国美女看黄片| 国产伦在线观看视频一区| 亚洲av成人av| 天堂影院成人在线观看| 国产高清激情床上av| 青春草国产在线视频 | 18禁在线无遮挡免费观看视频| 欧美变态另类bdsm刘玥| 国产av不卡久久| 亚洲国产精品久久男人天堂| 国产精品福利在线免费观看| 在线免费观看不下载黄p国产| 亚洲国产高清在线一区二区三| 日本免费a在线| 精品国内亚洲2022精品成人| av又黄又爽大尺度在线免费看 | 亚洲丝袜综合中文字幕| 九九久久精品国产亚洲av麻豆| 国产中年淑女户外野战色| 小说图片视频综合网站| 天天躁日日操中文字幕| 日韩国内少妇激情av| 久久精品综合一区二区三区| 成人一区二区视频在线观看| 久久久久久伊人网av| 国产v大片淫在线免费观看| 亚洲精品乱码久久久久久按摩| 亚洲精品国产av成人精品| 不卡视频在线观看欧美| 中国国产av一级| 国产探花在线观看一区二区| 五月玫瑰六月丁香| 国内久久婷婷六月综合欲色啪| 国产一区二区三区av在线 | 久久精品国产自在天天线| 久久人人爽人人片av| 精品一区二区三区人妻视频| 国产亚洲精品久久久com| 国产精品一区二区三区四区免费观看| 深爱激情五月婷婷| 亚洲国产日韩欧美精品在线观看| 国产一区二区激情短视频| 亚洲国产精品国产精品| 亚洲精品日韩av片在线观看| 国产单亲对白刺激| 国内少妇人妻偷人精品xxx网站| 亚洲精品国产av成人精品| 婷婷亚洲欧美| 丝袜美腿在线中文| 国内精品一区二区在线观看| 狠狠狠狠99中文字幕| 国产精品不卡视频一区二区| 噜噜噜噜噜久久久久久91| 精品国内亚洲2022精品成人| АⅤ资源中文在线天堂| 99热精品在线国产| 哪里可以看免费的av片| 国产av一区在线观看免费| 3wmmmm亚洲av在线观看| 成年免费大片在线观看| 人人妻人人澡欧美一区二区| 狂野欧美白嫩少妇大欣赏| 亚洲一区二区三区色噜噜| 午夜老司机福利剧场| 亚洲欧美日韩东京热| 国产精品爽爽va在线观看网站| 免费人成在线观看视频色| 久久久久久久久大av| 国产精品蜜桃在线观看 | 日韩 亚洲 欧美在线| 永久网站在线| 久久久色成人| 亚洲成人久久性| 午夜福利视频1000在线观看| 亚洲第一区二区三区不卡| 成人无遮挡网站| 网址你懂的国产日韩在线| 亚洲七黄色美女视频| 亚洲性久久影院| 97在线视频观看| 亚洲av成人精品一区久久| 国产成人91sexporn| 免费黄网站久久成人精品| 欧美日韩精品成人综合77777| 国产亚洲欧美98| 久久久午夜欧美精品| av女优亚洲男人天堂| 亚洲欧美精品自产自拍| 亚洲四区av| 你懂的网址亚洲精品在线观看 | 亚洲天堂国产精品一区在线| 国产精品,欧美在线| 在线播放国产精品三级| 99久久中文字幕三级久久日本| 国产精品永久免费网站| 国产蜜桃级精品一区二区三区| 国产精品久久电影中文字幕| 精品一区二区免费观看| 国产大屁股一区二区在线视频| 丰满的人妻完整版| 99久久成人亚洲精品观看| 色噜噜av男人的天堂激情| 如何舔出高潮| eeuss影院久久| 99久国产av精品| 哪个播放器可以免费观看大片| 成人二区视频| 少妇丰满av| 欧美极品一区二区三区四区| 国产伦精品一区二区三区视频9| 精品久久久久久成人av| 久久精品国产自在天天线| 特大巨黑吊av在线直播| 欧美日韩精品成人综合77777| 99久国产av精品| 哪个播放器可以免费观看大片| 国产一级毛片在线| 乱人视频在线观看| 国产免费男女视频| 免费观看人在逋| 国产免费男女视频| 日韩av在线大香蕉| 小说图片视频综合网站| 国产亚洲精品久久久com| 亚洲第一区二区三区不卡| 日本一二三区视频观看| 免费观看在线日韩| 国产真实乱freesex| 国产高潮美女av| 在现免费观看毛片| 日韩av不卡免费在线播放| 三级毛片av免费| 波多野结衣高清无吗| 久久精品国产亚洲网站| 久久人人精品亚洲av| 青春草国产在线视频 | 国产精品美女特级片免费视频播放器| 天天躁夜夜躁狠狠久久av| 日韩 亚洲 欧美在线| 国产三级中文精品| 亚洲成a人片在线一区二区| 日本与韩国留学比较| a级毛色黄片| 中文欧美无线码| 色视频www国产| 亚洲国产精品sss在线观看| 久久久国产成人精品二区| 91aial.com中文字幕在线观看| 色视频www国产| 91av网一区二区| 中文字幕制服av| 亚洲欧洲国产日韩| 日韩欧美国产在线观看| 久久精品影院6| 成人一区二区视频在线观看| 国产一区二区三区在线臀色熟女| 伊人久久精品亚洲午夜| 少妇猛男粗大的猛烈进出视频 | 久久精品国产亚洲网站| 亚洲成人中文字幕在线播放| 亚洲国产精品国产精品| 我的女老师完整版在线观看| 一级毛片aaaaaa免费看小| 精品久久国产蜜桃| 亚洲成av人片在线播放无| 嫩草影院精品99| 特级一级黄色大片| 六月丁香七月| 久久精品91蜜桃| 简卡轻食公司| 亚洲自拍偷在线| 欧美另类亚洲清纯唯美| 在线免费观看的www视频| 国产色爽女视频免费观看| 亚洲国产欧洲综合997久久,| 在线国产一区二区在线| 69人妻影院| 黄色日韩在线| 毛片女人毛片| 在线观看av片永久免费下载| 99热这里只有是精品在线观看| 美女大奶头视频| 九九爱精品视频在线观看| 99久久中文字幕三级久久日本| 日韩大尺度精品在线看网址| 九九热线精品视视频播放| 久久亚洲国产成人精品v| 一夜夜www| 人人妻人人澡人人爽人人夜夜 | 九色成人免费人妻av| 18+在线观看网站| 校园人妻丝袜中文字幕| 欧美日韩综合久久久久久| 免费观看的影片在线观看| 午夜福利高清视频| 国产亚洲欧美98| 少妇丰满av| 观看免费一级毛片| 欧美性猛交黑人性爽| 欧美人与善性xxx| 毛片一级片免费看久久久久| 身体一侧抽搐| 久久人人爽人人片av| 亚洲精品成人久久久久久| 国产av一区在线观看免费| 国产亚洲91精品色在线| 97超碰精品成人国产| 午夜亚洲福利在线播放| 欧美区成人在线视频| 久久午夜亚洲精品久久| 自拍偷自拍亚洲精品老妇| 成年av动漫网址| 亚洲自偷自拍三级| 少妇的逼水好多| 亚洲七黄色美女视频| 精品不卡国产一区二区三区| 中文资源天堂在线| 欧美精品国产亚洲| 国产伦在线观看视频一区| 岛国在线免费视频观看| 日本-黄色视频高清免费观看| 午夜激情欧美在线| 欧美bdsm另类| 男女那种视频在线观看| 国产精品精品国产色婷婷| 色综合色国产| 亚洲最大成人手机在线| 自拍偷自拍亚洲精品老妇| 蜜桃久久精品国产亚洲av| av又黄又爽大尺度在线免费看 | 亚洲经典国产精华液单| 国产真实乱freesex| 91aial.com中文字幕在线观看| 99久久人妻综合| 国产国拍精品亚洲av在线观看| 国产伦理片在线播放av一区 | 伊人久久精品亚洲午夜| 午夜福利视频1000在线观看| 国产亚洲5aaaaa淫片| 中文字幕制服av| 五月玫瑰六月丁香| 国产伦在线观看视频一区| 久久午夜福利片| 国产精品综合久久久久久久免费| 如何舔出高潮| 国产精品一区二区三区四区久久| 国产亚洲精品av在线| 亚洲精品影视一区二区三区av| 亚洲欧美精品自产自拍| 免费看av在线观看网站| 在线免费观看的www视频| 18+在线观看网站| 国产精品久久久久久av不卡| 精品熟女少妇av免费看| 久久午夜亚洲精品久久| 少妇人妻一区二区三区视频| 又黄又爽又刺激的免费视频.| 哪里可以看免费的av片| 日韩,欧美,国产一区二区三区 | 欧美最黄视频在线播放免费| 91在线精品国自产拍蜜月| 天堂中文最新版在线下载 | 身体一侧抽搐| 中出人妻视频一区二区| 国产男人的电影天堂91| 免费看光身美女| 毛片女人毛片| 成人三级黄色视频| 久久久精品大字幕| 中文字幕熟女人妻在线| 欧美日韩综合久久久久久| 热99re8久久精品国产| 91aial.com中文字幕在线观看| 国产亚洲av片在线观看秒播厂 | 在线天堂最新版资源| 久久午夜亚洲精品久久| 国内少妇人妻偷人精品xxx网站| 少妇人妻精品综合一区二区 | 伊人久久精品亚洲午夜| 精品人妻视频免费看| 亚洲人成网站在线播| 草草在线视频免费看| 欧美极品一区二区三区四区| 亚洲欧美成人综合另类久久久 | 国产精品人妻久久久久久| 床上黄色一级片| 亚洲欧美日韩高清在线视频| 久久精品综合一区二区三区| 午夜精品一区二区三区免费看| 哪里可以看免费的av片| 高清在线视频一区二区三区 | 亚洲人成网站在线播放欧美日韩| 免费观看的影片在线观看| 1000部很黄的大片| 99riav亚洲国产免费| 亚洲经典国产精华液单| 国产视频首页在线观看| 老女人水多毛片| 成人国产麻豆网| 日韩欧美国产在线观看| 少妇熟女aⅴ在线视频| 国产黄色小视频在线观看| 欧美人与善性xxx| 成人午夜高清在线视频| 三级经典国产精品| 国产亚洲精品久久久com| 国产av在哪里看| 日本在线视频免费播放| 久久精品影院6| 尾随美女入室| 床上黄色一级片| 国产精品99久久久久久久久| 97热精品久久久久久| 国产免费一级a男人的天堂| 熟女人妻精品中文字幕| 只有这里有精品99| 午夜激情福利司机影院| 你懂的网址亚洲精品在线观看 | 亚洲五月天丁香| 亚洲,欧美,日韩| 91久久精品国产一区二区成人| 一级黄色大片毛片| 简卡轻食公司| av国产免费在线观看| 又黄又爽又刺激的免费视频.| or卡值多少钱| 欧美不卡视频在线免费观看| 午夜福利视频1000在线观看| 亚洲精品乱码久久久久久按摩| 午夜福利在线观看吧| 亚洲成人久久性| 99久国产av精品| 久久亚洲精品不卡| 天堂网av新在线| 久久久精品欧美日韩精品| 国产白丝娇喘喷水9色精品| 国产成人a∨麻豆精品| 一级黄色大片毛片| 国产淫片久久久久久久久| 国产av在哪里看| 亚洲图色成人| 欧美一区二区国产精品久久精品| 国产成人freesex在线| 日韩高清综合在线| 日韩人妻高清精品专区| 欧美极品一区二区三区四区| 中文精品一卡2卡3卡4更新| 亚洲国产日韩欧美精品在线观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲经典国产精华液单| 99在线视频只有这里精品首页| 久久99蜜桃精品久久| 中文字幕制服av| 天堂√8在线中文| 欧美人与善性xxx| 一本久久精品| 亚洲av中文av极速乱| 国产亚洲欧美98| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | kizo精华| 亚洲婷婷狠狠爱综合网| 三级国产精品欧美在线观看| 99热6这里只有精品| 亚洲国产精品sss在线观看| 真实男女啪啪啪动态图| 久久久久久大精品| 白带黄色成豆腐渣| 亚洲成人久久爱视频| 最好的美女福利视频网| 久久精品国产亚洲av天美|