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

    華北典型區(qū)域土壤耕作方式對(duì)土壤特性和作物產(chǎn)量的影響*

    2019-11-11 09:08:16關(guān)劼兮陳素英邵立威張玉銘張喜英閆宗正
    關(guān)鍵詞:欒城實(shí)度夏玉米

    關(guān)劼兮, 陳素英, 邵立威, 張玉銘, 張喜英, 路 楊, 閆宗正

    華北典型區(qū)域土壤耕作方式對(duì)土壤特性和作物產(chǎn)量的影響*

    關(guān)劼兮1,2, 陳素英1**, 邵立威1, 張玉銘1, 張喜英1, 路 楊1,2, 閆宗正1,2

    (1. 中國(guó)科學(xué)院遺傳與發(fā)育生物學(xué)研究所農(nóng)業(yè)資源研究中心/中國(guó)科學(xué)院農(nóng)業(yè)水資源重點(diǎn)實(shí)驗(yàn)室/河北省節(jié)水農(nóng)業(yè)重點(diǎn)實(shí)驗(yàn)室 石家莊 050022; 2. 中國(guó)科學(xué)院大學(xué) 北京 100049)

    華北平原是我國(guó)重要的小麥玉米種植區(qū), 長(zhǎng)期土壤旋耕免耕和秸稈全量還田帶來(lái)耕層變淺、犁底層變厚和上移、土壤養(yǎng)分表聚等現(xiàn)象, 通過(guò)耕作方式改變, 解決上述問(wèn)題對(duì)維持區(qū)域糧食生產(chǎn)有重要意義。試驗(yàn)以冬小麥-夏玉米輪作系統(tǒng)為研究對(duì)象, 分別在代表華北平原高產(chǎn)區(qū)的欒城試驗(yàn)區(qū)和代表中低產(chǎn)區(qū)的南皮試驗(yàn)區(qū)進(jìn)行, 設(shè)置冬小麥播種前進(jìn)行土壤深耕、深松、窄深松3種處理, 以生產(chǎn)上常用的旋耕為對(duì)照。所有處理夏玉米季均采用土壤免耕播種, 測(cè)定項(xiàng)目包括土壤容重、作物根系、作物產(chǎn)量和水分利用效率。結(jié)果表明, 不同耕作方式對(duì)土壤特性和作物產(chǎn)量的影響具有區(qū)域差異。南皮試驗(yàn)區(qū)土壤深耕(松)顯著地(<0.05)提高了作物產(chǎn)量, 深耕、深松和窄深松處理的冬小麥產(chǎn)量比旋耕分別增加16.5%、19.3%和13.1%, 夏玉米產(chǎn)量分別增加17.3%、16.2%和21.9%, 周年產(chǎn)量分別增加16.9%、17.6%和17.8% ; 深耕、深松和窄深松處理間作物產(chǎn)量差異不顯著。欒城試驗(yàn)區(qū)冬小麥、夏玉米產(chǎn)量和周年產(chǎn)量各處理之間差異不顯著。土壤深耕、深松、窄深松和旋耕均能降低0~20 cm土層土壤緊實(shí)度和土壤容重。冬小麥播種后, 與土壤耕作前比較, 土壤深耕、深松和旋耕處理土壤緊實(shí)度南皮試驗(yàn)區(qū)分別平均降低71.6%和68.2%, 欒城試驗(yàn)區(qū)分別降低88.8%和-7.7%, 常用的旋耕模式在欒城試區(qū)沒(méi)有降低土壤緊實(shí)度。小麥?zhǔn)斋@時(shí)不同耕作方式0~40 cm土層的土壤容重均低于土壤耕作前的土壤容重, 至夏玉米收獲時(shí)不同耕作處理的土壤容重與耕作前基本一致, 不同耕作處理對(duì)土壤容重的影響差異不顯著。在南皮試驗(yàn)區(qū), 3種耕作方式與旋耕相比, 均顯著提高了冬小麥和夏玉米水分利用效率; 在欒城試驗(yàn)區(qū), 各處理冬小麥和夏玉米水分利用效率差異不顯著。本研究結(jié)果顯示在華北平原高產(chǎn)區(qū)連續(xù)實(shí)施土壤旋耕模式?jīng)]有影響作物產(chǎn)量, 而在中低產(chǎn)區(qū)實(shí)施土壤深耕或者深松模式更利于作物產(chǎn)量提高。

    小麥玉米一年兩熟; 高產(chǎn)區(qū); 中低產(chǎn)區(qū); 耕作方式; 土壤容重; 產(chǎn)量構(gòu)成; 水分利用效率

    華北平原是我國(guó)重要的糧食生產(chǎn)基地之一, 70%以上的耕地實(shí)施冬小麥()和夏玉米()一年兩熟輪作制。20世紀(jì)80年代開始, 隨著秸稈焚燒的禁止, 秸稈機(jī)械化還田和土壤旋耕作業(yè)迅速發(fā)展, 土壤耕作采用冬小麥旋耕后播種, 夏玉米免耕播種, 秸稈全量還田。由于長(zhǎng)期實(shí)施單一土壤免耕、旋耕和秸稈還田, 加上農(nóng)業(yè)全程機(jī)械化生產(chǎn)對(duì)土壤的壓實(shí)作用, 導(dǎo)致耕層變淺、犁底層加厚、土壤容重增大、耕層有效土量減少、土壤養(yǎng)分表聚等土壤退化特征突出[1-2], 嚴(yán)重阻礙了作物根系的深層分布和水肥吸收功能[3-5], 降低了作物的水肥資源利用率, 限制了作物產(chǎn)量的提高[6-9]。因此, 改進(jìn)土壤耕作方式對(duì)緩解土壤退化, 提升糧食持續(xù)生產(chǎn)具有重要意義。

    傳統(tǒng)的深耕(松)可以拓展土壤耕層深度, 打破犁底層, 改善耕層土壤物理性狀, 從而促進(jìn)作物生長(zhǎng), 提高作物的水分利用效率和產(chǎn)量[10-13]。通過(guò)深耕可以將聚集在上層土壤中多年未分解的秸稈翻壓到耕層以下, 加速其分解, 解決養(yǎng)分表聚問(wèn)題[14-15]。同時(shí), 深耕可減少作物土傳病蟲害的發(fā)生[16]; 深松可有效打破犁底層的限制, 促進(jìn)作物根系深扎[17], 降低土壤緊實(shí)度和容重[18], 構(gòu)建疏松的土壤耕層, 提高土壤的蓄水和保水能力, 提高作物產(chǎn)量[19]。通過(guò)深耕(松)改善土壤理化性狀、提高土壤蓄水保水保肥能力已成為該領(lǐng)域研究熱點(diǎn)。

    土壤深松技術(shù)在華北平原實(shí)施多年, 關(guān)于深松對(duì)土壤理化性狀、作物產(chǎn)量、水分利用效率、作物根系和深松深施肥等方面的影響前人已經(jīng)進(jìn)行了大量研究, 但大部分研究是同一地點(diǎn)進(jìn)行, 對(duì)不同區(qū)域不同土壤類型、氣候條件和灌溉方式下的比較研究較少, 不同區(qū)域土壤最適宜的土壤耕作模式也不明確。本文通過(guò)研究了不同區(qū)域土壤深耕(松)對(duì)作物產(chǎn)量、產(chǎn)量構(gòu)成、土壤容重和作物根系的影響, 提出適合不同區(qū)域的土壤耕作模式。

    1 材料與方法

    1.1 試驗(yàn)地點(diǎn)描述

    試驗(yàn)于2016年10月—2017年9月分別在華北平原的中國(guó)科學(xué)院欒城農(nóng)業(yè)生態(tài)系統(tǒng)試驗(yàn)站(以下簡(jiǎn)稱欒城試驗(yàn)區(qū), LC)和中國(guó)科學(xué)院南皮生態(tài)農(nóng)業(yè)試驗(yàn)站(以下簡(jiǎn)稱南皮試驗(yàn)區(qū), NP)進(jìn)行。根據(jù)河北省2015年農(nóng)村統(tǒng)計(jì)年鑒, 南皮縣的冬小麥和夏玉米產(chǎn)量分別為6 049.5 kg?hm-2和5 367.0 kg?hm-2, 欒城縣分別為7 606.8 kg?hm-2和7 876.8 kg?hm-2。南皮縣產(chǎn)量較低的主要原因一是該區(qū)歷史上土壤鹽漬化較嚴(yán)重, 20世紀(jì)90年代以后隨著區(qū)域地下水位下降, 土壤逐漸脫鹽, 但大部分土壤仍為輕度鹽漬化; 二是該區(qū)域農(nóng)田化肥投入一直處于較低水平, 2015年統(tǒng)計(jì)數(shù)據(jù)表明, 南皮縣的化肥投入平均為306.6 kg?hm-2, 欒城縣為476.8 kg?hm-2。根據(jù)作物產(chǎn)量水平, 欒城試驗(yàn)區(qū)代表華北平原的高產(chǎn)區(qū), 南皮試驗(yàn)區(qū)代表中低產(chǎn)區(qū)。

    欒城試驗(yàn)區(qū)位于37°53¢N、114°41¢E, 海拔50.1 m。屬于暖溫帶半濕潤(rùn)、半干旱季風(fēng)氣候, 雨熱同期, 年平均氣溫12.0 ℃, 年日照時(shí)數(shù)2 608 h, 多年平均降雨量480 mm, 土壤為褐土類灰黃土種, 土壤質(zhì)地以壤土為主, 地下水埋深42 m[20]。南皮試驗(yàn)區(qū)位于38°06′N、116°40′E, 海拔20 m。屬于暖溫帶半濕潤(rùn)季風(fēng)氣候區(qū), 年均氣溫12.3 ℃, 年日照時(shí)數(shù)2 938.6 h, 年均降水量524.5 mm, 土壤類型為潮土, 土壤質(zhì)地以壤土為主[21-22], 淺層地下水位埋深2~10 m, 為微咸水(礦化度≤3 g?L-1)[23]。兩個(gè)試驗(yàn)區(qū)的種植制度均為冬小麥-夏玉米一年兩熟制, 兩季秸稈長(zhǎng)期連續(xù)全量還田。冬小麥秸稈覆蓋還田, 即冬小麥機(jī)械收獲時(shí)同時(shí)將秸稈拋撒地表, 夏玉米免耕播種, 夏玉米生育期粉碎的小麥秸稈一直在地表覆蓋。夏玉米秸稈粉碎還田, 即夏玉米收獲后利用秸稈粉碎機(jī)機(jī)械粉碎后, 土壤旋耕(7~10 cm), 粉碎的玉米秸稈與表層0~10 cm土壤混合。兩個(gè)站點(diǎn)的試驗(yàn)地土壤養(yǎng)分背景值詳見表1。0~20 cm土壤有機(jī)質(zhì)、全氮、速效氮、速效鉀和速效磷含量欒城試驗(yàn)區(qū)比南皮試驗(yàn)區(qū)分別高43.1%、34.3%、57.1%、3.4%和-27.5%, 南皮試驗(yàn)區(qū)土壤速效磷含量高于欒城試驗(yàn)區(qū), 其他土壤養(yǎng)分含量均低于欒城試驗(yàn)區(qū)。

    表1 南皮試驗(yàn)區(qū)和欒城試驗(yàn)區(qū)土壤養(yǎng)分含量(2016年10月)

    1.2 試驗(yàn)處理

    兩個(gè)試驗(yàn)區(qū)的試驗(yàn)處理相同, 都是在連續(xù)多年秸稈全量還田和土壤旋耕的基礎(chǔ)上進(jìn)行。試驗(yàn)共設(shè)置4個(gè)處理: 1)土壤常規(guī)旋耕(CK), 9月底至10月初夏玉米收獲后, 用玉米秸稈粉碎機(jī)(165型)將玉米秸稈粉碎2遍, 用旋耕犁(1GKNB-420型變速)旋耕2遍, 深度為7~10 cm, 粉碎的玉米秸稈與表層0~10 cm土壤混合, 立即播種冬小麥; 2)土壤深耕(deep tillage, DT), 夏玉米秸稈粉碎后用1L-435深耕犁深耕1遍,深度為20~25 cm, 隨著表層土壤翻轉(zhuǎn), 粉碎的玉米秸稈分布于0~25 cm土壤, 旋耕1遍后播種小麥; 3)土壤深松(subsoiling, SS), 夏玉米秸稈粉碎后用深松犁(1S-200L)深松, 深度為30~35 cm, 2個(gè)深松犁之間的寬度為60 cm, 深松后旋耕1遍播種小麥, 深松處理的玉米秸稈分布于0~10 cm土壤; 4)土壤窄深松(narrow subsoiling, NSS): 作業(yè)方式和機(jī)具與SS處理相同, 由于2個(gè)深松犁之間的寬度為30 cm, 相對(duì)于SS處理較窄。4個(gè)處理的冬小麥采用2BXF-18小麥播種機(jī)播種, 播種后立即擦地。玉米采用2BYFSF玉米播種機(jī)播種。試驗(yàn)小區(qū)面積7 m × 8 m, 4次重復(fù)。

    1.3 試驗(yàn)管理

    欒城試驗(yàn)區(qū)和南皮試驗(yàn)區(qū)試驗(yàn)的冬小麥品種分別為‘科農(nóng)199’和‘衡5229’, 夏玉米品種分別為‘鄭單958’和‘華農(nóng)866’, 均為當(dāng)?shù)刂髟云贩N。灌溉方式欒城站為噴灌, 南皮站為畦灌, 小麥生育期灌溉均為拔節(jié)期灌溉1水, 灌溉時(shí)用水表計(jì)量, 灌溉量為 750 m3?hm-2。欒城試驗(yàn)區(qū)冬小麥底肥施用磷酸二胺復(fù)合肥和尿素, N為125 kg?hm-2, P2O5為158 kg?hm-2,拔節(jié)期追尿素N 105?hm-2; 玉米一次底施玉米專用肥(28-5-7) 715 kg?hm-2。南皮試驗(yàn)區(qū)冬小麥底肥施用磷酸二胺復(fù)合肥和尿素, N為140 kg?hm-2, P2O5為150 kg?hm-2, 拔節(jié)期追尿素N 140 kg?hm-2; 玉米一次底施玉米專用肥(28-5-7)900 kg?hm-2。

    欒城試驗(yàn)區(qū)小麥播種和收獲時(shí)間分別為2016年10月8日和2017年6月8日, 播種量為300 kg?hm-2; 夏玉米播種和收獲時(shí)間分別為2017年6月10日和2017年9月29日, 播種量為45 kg?hm-2。南皮試驗(yàn)區(qū)小麥播種和收獲時(shí)間分別為2016年10月21日和2017年6月10日, 播種量為450 kg?hm-2; 夏玉米的播種和收獲時(shí)間分別為2017年6月15日和2017年10月5日, 播種量為45 kg?hm-2。

    1.4 參數(shù)測(cè)定及計(jì)算方法

    土壤容重: 分別在試驗(yàn)開始前、冬小麥?zhǔn)斋@期和夏玉米收獲期, 用環(huán)刀法分層測(cè)定0~40 cm土層土壤容重, 每10 cm一層, 每個(gè)處理測(cè)定3次重復(fù)。將不銹鋼環(huán)刀中的原狀土樣轉(zhuǎn)移到鋁盒中, 在 105 ℃下烘 24 h 稱重, 計(jì)算單位體積中烘干土的質(zhì)量。

    土壤緊實(shí)度: 在試驗(yàn)開始前和冬小麥播種后, 采用SC-900土壤堅(jiān)實(shí)度儀測(cè)定(SPECTRUM, USA), 測(cè)定總深度為0~45 cm, 每2.5 cm測(cè)定1個(gè)數(shù)據(jù), 每個(gè)處理測(cè)定6次重復(fù)。

    小麥根長(zhǎng)密度: 在冬小麥成熟時(shí),用直徑10 cm的根鉆取0~10 cm土壤樣品, 取樣時(shí)根鉆中心置于作物行上和行間交界處, 10 cm為1個(gè)層次, 共10個(gè)深度。每個(gè)處理層4個(gè)重復(fù)。取樣后立即帶回實(shí)驗(yàn)室, 用0.25 mm土壤篩洗凈根系, 去除雜質(zhì), 然后用1.27 cm刻度根盤測(cè)量根長(zhǎng)。用測(cè)定的不同深度根長(zhǎng)值, 計(jì)算根長(zhǎng)密度(RLD)。根長(zhǎng)密度定義為單位土壤體積中的根長(zhǎng), 單位為cm?cm-3。不同取樣地點(diǎn)某一特定深度的根長(zhǎng)密度為4個(gè)取樣重復(fù)的平均值。

    土壤含水量: 采用土鉆法取土烘干, 每20 cm為1層, 每個(gè)處理測(cè)定3次重復(fù)。

    農(nóng)田耗水(ET)根據(jù)水量平衡公式計(jì)算:

    ET=SWD+++--R (1)

    式中: ET為作物實(shí)際耗水量(mm), SWD為生育期土壤水分消耗量(mm),為降雨量(mm),為灌溉量(mm),為土壤毛細(xì)管提升水(mm),為水分下滲量(mm),為地表徑流量(mm)。試驗(yàn)期間沒(méi)有發(fā)現(xiàn)地表徑流,忽略不計(jì); 欒城試驗(yàn)區(qū)水位埋深40 m, 南皮試驗(yàn)區(qū)的地下水位埋深雖然較淺, 本文中只計(jì)算了0~1 m的土壤消耗, 土壤毛細(xì)管提升值可以忽略。在本試驗(yàn)中和值設(shè)定為0。因此, 公式(1)可以簡(jiǎn)化為ET=SWD++–。公式中的根區(qū)土壤水分下滲量()根據(jù)文獻(xiàn)[24]計(jì)算。

    水分利用效率(WUE, kg?cm-3)計(jì)算公式為:

    WUE=/ET (2)

    式中:為經(jīng)濟(jì)產(chǎn)量(kg?hm-2), ET為農(nóng)田耗水量(mm)。

    冬小麥和夏玉米產(chǎn)量、生物量、產(chǎn)量構(gòu)成測(cè)定: 冬小麥和夏玉米收獲時(shí)各小區(qū)單獨(dú)收獲, 脫粒曬干后計(jì)算小麥和玉米籽粒產(chǎn)量, 折算成單產(chǎn)。并隨機(jī)在每個(gè)小區(qū)內(nèi)選擇60穗小麥和3株玉米考種, 考種項(xiàng)目包括生物量、穗數(shù)、穗粒數(shù)、千粒重。

    1.5 氣象數(shù)據(jù)

    欒城試驗(yàn)區(qū)的氣象數(shù)據(jù)來(lái)自本試驗(yàn)站建有的標(biāo)準(zhǔn)氣象站, 氣象站距離試驗(yàn)地50 m。南皮試驗(yàn)區(qū)的氣象數(shù)據(jù)來(lái)自距離該試驗(yàn)站最近的國(guó)家標(biāo)準(zhǔn)氣象站——泊頭氣象站(38°08¢N, 116°55¢E, 海拔13.2 m)。

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

    各處理取多次重復(fù)測(cè)定的數(shù)據(jù)計(jì)算平均值, 并計(jì)算標(biāo)準(zhǔn)偏差; 方差和相關(guān)性分析用SPSS 11.0進(jìn)行, 用SPSS軟件分析不同處理間顯著性。用Excel軟件做圖, 處理之間的顯著性分析均在<0.05水平。

    2 結(jié)果與分析

    2.1 試驗(yàn)期間氣象條件

    南皮和欒城試驗(yàn)區(qū)氣象條件存在較大差異。表2為試驗(yàn)期間冬小麥和夏玉米生育期主要?dú)庀笠蛩乇容^。可以看出, 兩試驗(yàn)區(qū)冬小麥生育期≥10 ℃積溫、日均溫和降雨量基本相同; 日照時(shí)數(shù)相差較大, 南皮試驗(yàn)區(qū)生育期總?cè)照諘r(shí)數(shù)比欒城試驗(yàn)區(qū)高20.0%。夏玉米生育期≥10 ℃積溫基本相同, 南皮試驗(yàn)區(qū)比欒城試驗(yàn)區(qū)降雨量高116.2%, 日均溫高20.2%, 日照時(shí)數(shù)高8.3%。冬小麥和夏玉米生育期每月日照時(shí)數(shù)和每天的溫度日較差見圖1, 南皮試驗(yàn)區(qū)的月日照時(shí)數(shù)和溫度日較差均高于欒城試驗(yàn)區(qū)。

    表2 南皮試驗(yàn)區(qū)和欒城試驗(yàn)區(qū)冬小麥和夏玉米生育期氣象因素

    圖1 南皮試驗(yàn)區(qū)(南皮)和欒城試驗(yàn)區(qū)(欒城)逐月日照時(shí)數(shù)和溫度日較差(2016年10月至2017年9月)

    2.2 耕作方式對(duì)土壤緊實(shí)度和容重的影響

    土壤緊實(shí)度和土壤容重是土壤耕性的重要指標(biāo), 緊實(shí)度和容重大的土壤影響水分入滲、根系下扎和肥料利用效率。2016年11月小麥播種后進(jìn)行了土壤緊實(shí)度的測(cè)定, 結(jié)果表明, 土壤深耕(松)可以降低土壤耕層的緊實(shí)度(圖2)。在南皮試驗(yàn)區(qū), 與土壤耕作前相比, 不同土壤耕作均降低0~20 cm土壤緊實(shí)度, 深耕、深松、窄深松和旋耕分別降低69.7%、72.7%、72.5%和68.2%, 處理間差異不顯著。20~45 cm分別降低6.7%、14.0%、28.7%和-7.0%, 土壤旋耕處理的緊實(shí)度大于耕作前, 窄深松處理與其他3個(gè)處理間差異顯著(<0.05), 其他3個(gè)處理間差異不顯著。0~20 cm土壤緊實(shí)度降低程度大于20~45 cm。在欒城試驗(yàn)區(qū), 與土壤耕作前相比, 不同土壤耕作均降低0~20 cm和20~45 cm的土壤緊實(shí)度。0~20 cm深耕、深松、窄深松和旋耕分別降低52.0%、53.3%、161.0%和-7.7%, 土壤旋耕緊實(shí)度大于土壤耕作前, 窄深松的緊實(shí)度下降最明顯。20~45 cm深耕、深松、窄深松和旋耕分別降低31.0%、35.1%、65.0%和41.7%, 土壤緊實(shí)度降低程度大于南皮試驗(yàn)區(qū)。兩個(gè)區(qū)域土壤緊實(shí)度比較, 南皮的土壤緊實(shí)度最大值位于2 000~2 500 kPa, 欒城的土壤緊實(shí)度最大值位于2 500~3 000 kPa, 欒城的土壤緊實(shí)度最大值大于南皮。

    土壤深耕(松)可以降低土壤耕層的容重(圖3)。南皮試驗(yàn)區(qū)冬小麥播種前(2016年10月)進(jìn)行土壤深耕、深松、窄深松和旋耕后, 冬小麥?zhǔn)斋@時(shí)(2017年6月)的土壤容重與播種前比較, 0~40 cm土壤容重均出現(xiàn)了顯著降低(<0.05); 夏玉米收獲時(shí)(2017年10月)土壤容重與播種前相差不大。土壤深耕和深松處理在夏玉米收獲時(shí), 20 cm土層土壤容重仍低于播種前, 而窄深松和旋耕處理則與播種前的土壤容重基本相同。欒城試驗(yàn)區(qū)的土壤容重變化與南皮站基本一致, 小麥?zhǔn)斋@時(shí)與播種前相比土壤容重都出現(xiàn)了降低, 但降低不顯著(<0.05)。夏玉米收獲時(shí)土壤容重與播種前的容重相差不大, 20 cm處土壤深耕(松)處理的土壤容重低于播種前處理。

    2.3 耕作方式對(duì)冬小麥根系的影響

    根長(zhǎng)密度是反映作物根系生長(zhǎng)發(fā)育狀況的指標(biāo)之一, 不同耕作方式下冬小麥?zhǔn)斋@時(shí)的根長(zhǎng)密度詳見圖4??梢钥闯? 兩個(gè)試驗(yàn)區(qū)冬小麥根長(zhǎng)密度分布趨勢(shì)一致, 均表現(xiàn)為表層大于下層。0~100 cm平均根長(zhǎng)密度南皮試驗(yàn)區(qū)為1.26 cm?cm-3, 欒城試驗(yàn)區(qū)為1.42 cm?cm-3, 欒城試驗(yàn)區(qū)比南皮高12.7%。分層來(lái)看0~20 cm土層平均根長(zhǎng)密度兩試驗(yàn)區(qū)基本相同, 但20~40 cm和40~100 cm平均根長(zhǎng)密度欒城比南皮分別高20.0%和21.4%。從耕作方式對(duì)冬小麥根長(zhǎng)密度的影響來(lái)看, 欒城試驗(yàn)區(qū)和南皮試驗(yàn)區(qū)耕作方式對(duì)冬小麥根長(zhǎng)密度的影響均不顯著。

    圖2 不同耕作措施對(duì)南皮試驗(yàn)區(qū)(a)和欒城試驗(yàn)區(qū)(b)小麥苗期土壤緊實(shí)度的影響(2016年11月)

    CK: rotary tillage; DT: deep tillage; SS: subsoiling; NSS: narrow subsoiling.

    圖3 不同耕作措施對(duì)南皮試驗(yàn)區(qū)和欒城試驗(yàn)區(qū)冬小麥和夏玉米收獲時(shí)土壤容重的影響

    DT: deep tillage; SS: subsoiling; NSS: narrow subsoiling; CK: rotary tillage.

    2.4 耕作方式對(duì)土壤水分和作物水分利用效率的影響

    土壤水分是供給作物用水的重要來(lái)源。不同耕作方式對(duì)土壤水分的影響具有相同的趨勢(shì), 圖5為代表性的深耕和旋耕處理對(duì)土壤水分的影響??梢钥闯? 無(wú)論是旋耕還是深耕處理, 小麥從播種到收獲, 南皮試驗(yàn)區(qū)0~2 m土壤水分高于欒城試驗(yàn)區(qū)。小麥?zhǔn)斋@時(shí), 欒城0~2 m土壤含水量從上至下為10%左右, 整體處于較低狀態(tài); 南皮試驗(yàn)區(qū)0~1 m土層含水量很低, 表層基本達(dá)到凋萎濕度, 但1~2 m土壤含水量較高。至夏玉米收獲時(shí), 隨著夏玉米生育期灌溉和降雨的補(bǔ)充, 土壤含水量整體比冬小麥?zhǔn)斋@時(shí)提高, 南皮試驗(yàn)區(qū)由于降雨遠(yuǎn)高于欒城試驗(yàn)區(qū), 其1~2 m土層的土壤含水量較高。

    圖4 耕作方式對(duì)欒城試驗(yàn)區(qū)和南皮試驗(yàn)區(qū)冬小麥?zhǔn)斋@期根長(zhǎng)密度的影響

    CK: rotary tillage; DT: deep tillage; SS: subsoiling; NSS: narrow subsoiling.

    圖5 旋耕和深耕處理下南皮試驗(yàn)區(qū)和欒城試驗(yàn)區(qū)不同時(shí)間土壤含水量的變化

    在南皮試驗(yàn)區(qū), 深耕(松)處理顯著增加了冬小麥和夏玉米的水分利用效率(表3)。與CK相比, 深耕、深松和窄深松處理使冬小麥的水分利用效率分別提高27.8%、23.8%和25.8%, 夏玉米的水分利用效率分別提高14.1%、11.6%和18.7%, 周年水分利用效率分別提高20.1%、17.2%和22.4%, 與CK處理均達(dá)到顯著水平(<0.05)。在欒城試驗(yàn)區(qū), 深耕(松)處理對(duì)冬小麥和夏玉米水分利用效率的影響不顯著。

    表3 不同耕作方式對(duì)南皮試驗(yàn)區(qū)和欒城試驗(yàn)區(qū)冬小麥和夏玉米水分利用效率的影響

    同列不同小寫字母表示不同處理間在<0.05水平差異顯著。Different lowercase letters in the same column mean significant differences among treatments at 0.05 level.

    2.5 耕作方式對(duì)冬小麥和夏玉米產(chǎn)量及產(chǎn)量構(gòu)成因素的影響

    對(duì)長(zhǎng)期秸稈還田和土壤旋耕的農(nóng)田, 冬小麥播種前進(jìn)行土壤深耕和深松處理, 兩個(gè)試驗(yàn)區(qū)域?qū)Χ←満拖挠衩桩a(chǎn)量的影響不同(表4)。

    南皮試驗(yàn)區(qū)不同耕作方式對(duì)冬小麥和夏玉米均具有增產(chǎn)效應(yīng), 夏玉米的增產(chǎn)效應(yīng)稍高于冬小麥。深耕、深松和窄深松處理冬小麥產(chǎn)量分別比CK顯著增加16.5%、19.3%和13.1%(<0.05), 但深耕、深松和窄深松處理間差異不顯著。土壤耕作方式對(duì)夏玉米產(chǎn)量也具有顯著(<0.05)增產(chǎn)作用, 深耕、深松和窄深松處理分別比CK增產(chǎn)17.3%、16.2%和21.9%。深松和深耕均提高了夏玉米的穗數(shù), 深耕、深松和窄深松處理穗數(shù)分別比CK提高8.0%、6.2%和21.8%, 窄深松處理與CK差異顯著。與CK相比不同耕作處理夏玉米千粒重差異不顯著。

    欒城試驗(yàn)區(qū)不同土壤耕作方式對(duì)冬小麥和夏玉米的增產(chǎn)效應(yīng)不顯著。窄深松處理的冬小麥產(chǎn)量最高, 比CK增產(chǎn)7.6%, 深耕、深松和CK之間產(chǎn)量差異不顯著, 與窄深松產(chǎn)量差異顯著(<0.05)。土壤深耕松對(duì)玉米產(chǎn)量的增產(chǎn)效果不顯著。深耕、深松和窄深松與CK處理的夏玉米產(chǎn)量差異不顯著。深耕和窄深松均提高了夏玉米的穗數(shù)(<0.05), 分別比CK提高17.9%和14.5%, 深松處理的穗數(shù)高于CK, 差異不顯著。

    南皮試驗(yàn)區(qū)土壤深耕、深松和窄深松的周年產(chǎn)量分別比CK顯著(<0.05)提高16.9%、17.6%和17.8%; 欒城試驗(yàn)區(qū)增產(chǎn)不顯著。

    表4 不同耕作處理對(duì)冬小麥和夏玉米產(chǎn)量及產(chǎn)量構(gòu)成的影響

    CK: 常規(guī)旋耕; DT: 深耕; SS: 深松; NSS: 窄深松。同列同一試驗(yàn)區(qū)不同小寫字母表示<0.05水平差異顯著。CK: rotary tillage; DT: deep tillage; SS: subsoiling; NSS: narrow subsoiling. Different lowercase letters in the same column for the same experiment site mean significant differences among treatments at 0.05 level.

    3 結(jié)論與討論

    3.1 討論

    研究結(jié)果表明, 針對(duì)華北平原冬小麥-夏玉米一年兩熟制長(zhǎng)期單一旋耕和免耕的農(nóng)田, 在冬小麥播種前進(jìn)行深耕(松)處理不僅提高了當(dāng)季冬小麥的產(chǎn)量, 對(duì)下茬夏玉米也具有增產(chǎn)作用, 這與其他同類型地區(qū)的研究結(jié)果一致[9,25]。但是耕作方式有區(qū)域適應(yīng)性, 土壤耕作對(duì)產(chǎn)量的效應(yīng)受土壤、氣候等環(huán)境因子等影響, 在不同區(qū)域的增產(chǎn)效應(yīng)不同[26]。李永平等[27]研究發(fā)現(xiàn), 深耕和深松處理比常規(guī)旋耕分別提高玉米產(chǎn)量20.4%和3.7%; 梁金鳳等[28]研究不同耕作方式對(duì)土壤理化性狀及玉米根系生長(zhǎng)和產(chǎn)量的影響, 發(fā)現(xiàn)耕作方式對(duì)玉米的產(chǎn)量無(wú)顯著影響。本研究2個(gè)試驗(yàn)區(qū)的不同耕作方式對(duì)產(chǎn)量的影響也不相同。南皮試驗(yàn)區(qū), 深耕(松)對(duì)冬小麥和夏玉米的增產(chǎn)效應(yīng)均達(dá)顯著水平(<0.05), 欒城試驗(yàn)區(qū)深耕(松)處理的冬小麥和夏玉米產(chǎn)量與常規(guī)旋耕(CK)的差異不顯著, 增產(chǎn)效應(yīng)不顯著。南皮試區(qū)和欒城試區(qū)的差異表現(xiàn)在以下幾個(gè)方面:

    1)耕作方式對(duì)土壤容重和作物根系的影響。與欒城試驗(yàn)區(qū)相比, 南皮試驗(yàn)區(qū)的土壤肥力較低。0~20 cm和20~40 cm土壤有機(jī)質(zhì)比欒城試驗(yàn)區(qū)分別低47.2%和44.8%, 全氮低34.1%和58.3%, 速效氮低63.6%和70.9%, 速效鉀低5.1%和25.3%; 僅速效磷含量南皮試驗(yàn)區(qū)高于欒城試驗(yàn)區(qū), 0~20 cm和20~40 cm土層分別高34.7%和34.4%。在土壤肥力較低的情況下, 土壤深松和深耕顯著降低了土壤緊實(shí)度和土壤容重, 促進(jìn)了作物根系的生長(zhǎng)和作物產(chǎn)量, 深耕、深松和窄深松的作物產(chǎn)量異不顯著, 但與常規(guī)旋耕之間差異顯著。欒城試驗(yàn)區(qū)由于長(zhǎng)期采用噴灌灌溉, 1~2 m土壤含水量難以得到補(bǔ)充, 深層土壤含水量較低, 雖然促進(jìn)了冬小麥下層根系的生長(zhǎng), 但是生育后期土壤水分供應(yīng)不足, 影響根系吸水和產(chǎn)量提升。

    2)土壤水分條件和水分利用效率。冬小麥生育期降雨量不能滿足其生長(zhǎng)發(fā)育的需求, 并且隨著水資源日益短缺, 華北平原區(qū)的欒城和南皮試驗(yàn)區(qū)冬小麥生產(chǎn)均采取了調(diào)虧灌溉, 因此, 生育期的土壤水分供應(yīng)能力對(duì)冬小麥穩(wěn)產(chǎn)和高產(chǎn)起決定作用。南皮試驗(yàn)區(qū)冬小麥生育期土壤水分含量高于欒城試驗(yàn)區(qū), 尤其是冬小麥生育后期, 深層土壤水分高低對(duì)冬小麥的灌漿具有重要作用。南皮試驗(yàn)區(qū)冬小麥生育期深耕、深松、窄深松和旋耕(CK)的土壤耗水分別為169.1 mm、177.3 mm、150.4 mm和194.7 mm, CK的土壤耗水量最高, 但處理間差異不顯著。在深(松)耕處理增產(chǎn)的前提下, 水分利用效率也相應(yīng)提高, 與CK間達(dá)到顯著水平。欒城試驗(yàn)區(qū)處理間的土壤耗水和產(chǎn)量均沒(méi)有差異, 最終水分利用效率差異也不顯著。

    從不同土壤耕作方式的作業(yè)過(guò)程和動(dòng)力消耗來(lái)看, 深耕、深松的作業(yè)程序和動(dòng)力消耗基本相同, 深耕或者深松后進(jìn)行旋耕播種, 窄深松則相當(dāng)于深松2次進(jìn)行旋耕后播種, 動(dòng)力消耗較高, 深耕、深松和窄深松對(duì)土壤特性和作物產(chǎn)量的影響之間差異不顯著。因此, 生產(chǎn)上不推薦使用窄深松。

    3.2 結(jié)論

    本研究結(jié)果表明, 針對(duì)華北平原冬小麥-夏玉米一年兩熟種植中長(zhǎng)期采用單一土壤旋耕和免耕農(nóng)田, 實(shí)施土壤深耕、深松和窄深松處理對(duì)作物產(chǎn)量的影響不同。在代表華北平原高產(chǎn)區(qū)土壤深耕(松)對(duì)冬小麥、夏玉米和全年產(chǎn)量的增產(chǎn)效應(yīng)不顯著, 可以繼續(xù)實(shí)施土壤旋耕作業(yè)。但中低產(chǎn)區(qū)域土壤深耕(松)多少年后對(duì)降低作物產(chǎn)量和高產(chǎn)區(qū)需要不需要進(jìn)行土壤深耕(松)需要做進(jìn)一步研究。

    [1] 杜章留, 高偉達(dá), 陳素英, 等. 保護(hù)性耕作對(duì)太行山前平原土壤質(zhì)量的影響[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 2011, 19(5): 1134–1142 DU Z L, GAO W D, CHEN S Y, et al. Effect of conservation tillage on soil quality in the piedmont plain of Mount Taihang[J]. Chinese Journal of Eco-Agriculture, 2011, 19(5): 1134–1142

    [2] BECERRA A T, BOTTA G F, BRAVO X L, et al. Soil compaction distribution under tractor traffic in almond (L.) orchard in Almería Espa?a[J]. Soil and Tillage Research, 2010, 107(1): 49–56

    [3] BENGOUGH A G, MCKENZIE B M, HALLETT P D, et al. Root elongation, water stress, and mechanical impedance: a review of limiting stresses and beneficial root tip traits[J]. Journal of Experimental Botany, 2011, 62(1): 59–68

    [4] LIPIEC J, HATANO R. Quantification of compaction effects on soil physical properties and crop growth[J]. Geoderma, 2003, 116(1/2): 107–136

    [5] 李潮海, 李勝利, 王群, 等. 下層土壤容重對(duì)玉米根系生長(zhǎng)及吸收活力的影響[J]. 中國(guó)農(nóng)業(yè)科學(xué), 2005, 38(8): 1706–1711 LI C H, LI S L, WANG Q, et al. A study on corn root growth and activities at different soil layers with special bulk density[J]. Scientia Agricultura Sinica, 2005, 38(8): 1706–1711

    [6] VAN DEN AKKER J J H, ARVIDSSON J, HORN R. Introduction to the special issue on experiences with the impact and prevention of subsoil compaction in the European Union[J]. Soil and Tillage Research, 2003, 73(1/2): 1–8

    [7] POWERS R F, ANDREW SCOTT D, SANCHEZ F G, et al. The North American long-term soil productivity experiment: Findings from the first decade of research[J]. Forest Ecology and Management, 2005, 220(1/2/3): 31–50

    [8] 謝迎新, 靳海洋, 孟慶陽(yáng), 等. 深耕改善砂姜黑土理化性狀提高小麥產(chǎn)量[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2015, 31(10): 167–173 XIE Y X, JIN H Y, MENG Q Y, et al. Deep tillage improving physical and chemical properties of soil and increasing grain yield of winter wheat in lime concretion black soil farmland[J]. Transactions of the CSAE, 2015, 31(10): 167–173

    [9] 趙亞麗, 劉衛(wèi)玲, 程思賢, 等. 深松(耕)方式對(duì)砂姜黑土耕層特性、作物產(chǎn)量和水分利用效率的影響[J]. 中國(guó)農(nóng)業(yè)科學(xué), 2018, 51(13): 2489–2503 ZHAO Y L, LIU W L, CHENG S X, et al. Effects of pattern of deep tillage on topsoil features, yield and water use efficiency in lime concretion black soil[J]. Scientia Agricultura Sinica, 2018, 51(13): 2489–2503

    [10] 李霞. 小麥-玉米周年生產(chǎn)體系中播前耕作對(duì)夏玉米產(chǎn)量及生理特性的影響[D]. 泰安: 山東農(nóng)業(yè)大學(xué), 2014 LI X. Effect of pre-planting tillage of winter wheat and summer maize on yield and physiological characteristics of summer maize[D]. Tai’an: Shandong Agricultural University, 2014

    [11] 孔曉民, 韓成衛(wèi), 曾蘇明, 等. 不同耕作方式對(duì)土壤物理性狀及玉米產(chǎn)量的影響[J]. 玉米科學(xué), 2014, 22(1): 108–113 KONG X M, HAN C W, ZENG S M, et al. Effects of different tillage managements on soil physical properties and maize yield[J]. Journal of Maize Sciences, 2014, 22(1): 108–113

    [12] 楊永輝, 武繼承, 張潔梅, 等. 耕作方式對(duì)土壤水分入滲、有機(jī)碳含量及土壤結(jié)構(gòu)的影響[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 2017, 25(2): 258–266 YANG Y H, WU J C, ZHANG J M, et al. Effect of tillage method on soil water infiltration, organic carbon content and structure[J]. Chinese Journal of Eco-Agriculture, 2017, 25(2): 258–266

    [13] 劉衛(wèi)玲, 程思賢, 周金龍, 等. 深松(耕)時(shí)機(jī)與方式對(duì)土壤物理性狀和玉米產(chǎn)量的影響[J]. 河南農(nóng)業(yè)科學(xué), 2018, 47(3): 7–13 LIU W L, CHENG S X, ZHOU J L, et al. Effects of time and pattern of deep tillage (subsoiling) on soil physical properties and maize yield[J]. Journal of Henan Agricultural Sciences, 2018, 47(3): 7–13

    [14] 于淑婷, 趙亞麗, 王育紅, 等. 輪耕模式對(duì)黃淮海冬小麥-夏玉米兩熟區(qū)農(nóng)田土壤改良效應(yīng)[J]. 中國(guó)農(nóng)業(yè)科學(xué), 2017, 50(11): 2150–2165YU S T, ZHAO Y L, WANG Y H, et al. Improvement effects of rotational tillage patterns on soil in the winter wheat-summer maize double cropping area of Huang-Huai-Hai region[J]. Scientia Agricultura Sinica, 2017, 50(11): 2150–2165

    [15] 謝迎新, 靳海洋, 孟慶陽(yáng), 等. 深耕改善砂姜黑土理化性狀提高小麥產(chǎn)量[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2015, 31(10): 167–173 XIE Y X, JIN H Y, MENG Q Y, et al. Deep tillage improving physical and chemical properties of soil and increasing grain yield of winter wheat in lime concretion black soil farmland[J]. Transactions of the CSAE, 2015, 31(10): 167–173

    [16] 趙俊靖, 李洪奎, 李華春, 等. 深耕對(duì)小麥病蟲害的防治效果[J]. 植物醫(yī)生, 2017, 30(11): 51–53ZHAO J J, LI H K, LI H C, et al. Effects of deep tillage on wheat diseases and insect pests[J]. Plant Doctor, 2017, 30(11): 51–53

    [17] 雷友, 曹國(guó)鑫, 牛新勝, 等. 土壤深耕對(duì)冬小麥根系在土壤剖面分布的影響[J]. 現(xiàn)代農(nóng)業(yè)科技, 2011(8): 272–273 LEI Y, CAO G X, NIU X S, et al. Influence of deep tillage on distribution of wheat root system in soil profile[J]. Modern Agricultural Sciences and Technology, 2011(8): 272–273

    [18] 宮秀杰, 錢春榮, 于洋, 等. 深松免耕技術(shù)對(duì)土壤物理性狀及玉米產(chǎn)量的影響[J]. 玉米科學(xué), 2009, 17(5): 134–137 GONG X J, QIAN C R, YU Y, et al. Effects of subsoiling and No-tillage on soil physical characters and corn yield[J]. Journal of Maize Sciences, 2009, 17(5): 134–137

    [19] CAI H G, MA W, ZHANG X Z, et al. Effect of subsoil tillage depth on nutrient accumulation, root distribution, and grain yield in spring maize[J]. The Crop Journal, 2014, 2(5): 297–307

    [20] ZHANG Y C, SHEN Y J, SUN H Y, et al. Evapotranspiration and its partitioning in an irrigated winter wheat field: A combined isotopic and micrometeorologic approach[J]. Journal of Hydrology, 2011, 408(3/4): 203–211

    [21] LIU X W, FEIKE T, CHEN S Y, et al. Effects of saline irrigation on soil salt accumulation and grain yield in the winter wheat-summer maize double cropping system in the low plain of North China[J]. Journal of Integrative Agriculture, 2016, 15(12): 2886–2898

    [22] 呂夢(mèng)宇, 王仕琴, 齊永青, 等. 華北低平原區(qū)降水與坑塘蓄水響應(yīng)關(guān)系研究——以河北省南皮縣為例[J]. 自然資源學(xué)報(bào), 2018, 33(10): 1796–1805 LYU M Y, WANG S Q, QI Y Q, et al. Study on the relationship between precipitation and pond water storage in lowland area of North China Plain — A case study in Nanpi County, Hebei Province[J]. Journal of Natural Resources, 2018, 33(10): 1796–1805

    [23] 孫宏勇, 劉小京, 張喜英. 鹽堿地水鹽調(diào)控研究[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 2018, 26(10): 1528–1536 SUN H Y, LIU X J, ZHANG X Y. Regulations of salt and water of saline-alkali soil: A review[J]. Chinese Journal of Eco-Agriculture, 2018, 26(10): 1528–1536

    [24] ZHANG X Y, CHEN S Y, SUN H Y, et al. Changes in evapotranspiration over irrigated winter wheat and maize in North China Plain over three decades[J]. Agricultural Water Management, 2011, 98(6): 1097–1104

    [25] 王群, 王建, 張學(xué)林, 等. 不同耕作模式下小麥玉米周年生產(chǎn)及土壤養(yǎng)分變化特征[J]. 河南農(nóng)業(yè)大學(xué)學(xué)報(bào), 2015, 49(4): 429–437 WANG Q, WANG J, ZHANG X L, et al. Change characteristics of wheat and maize anniversary production and soil nutrient content under different rotation tillage patterns[J]. Journal of Henan Agricultural University, 2015, 49(4): 429–437

    [26] TOLIVER D K, LARSON J A, ROBERTS R K, et al. Effects of no-till on yields as influenced by crop and environmental factors[J]. Agronomy Journal, 2012, 104(2): 530

    [27] 李永平, 王孟本, 史向遠(yuǎn), 等. 不同耕作方式對(duì)土壤理化性狀及玉米產(chǎn)量的影響[J]. 山西農(nóng)業(yè)科學(xué), 2012, 40(7): 723–727 LI Y P, WANG M B, SHI X Y, et al. Influence of different tillage methods on soil physical and chemical properties and maize yield[J]. Journal of Shanxi Agricultural Sciences, 2012, 40(7): 723–727

    [28] 梁金鳳, 齊慶振, 賈小紅, 等. 不同耕作方式對(duì)土壤性質(zhì)與玉米生長(zhǎng)的影響研究[J]. 生態(tài)環(huán)境學(xué)報(bào), 2010, 19(4): 945–950 LIANG J F, QI Q Z, JIA X H, et al. Effects of different tillage managements on soil properties and corn growth[J]. Ecology and Environmental Sciences, 2010, 19(4): 945–950

    Soil tillage practices affecting the soil characteristics and yield of winter wheat and summer maize in North China*

    GUAN Jiexi1,2, CHEN Suying1**, SHAO Liwei1, ZHANG Yuming1, ZHANG Xiying1, LU Yang1,2, YAN Zongzheng1,2

    (1. Center for Agricultural Resources Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences / Key Laboratory of Agricultural Water Resources, Chinese Academy of Sciences / Hebei Key Laboratory of Agricultural Water-saving, Shijiazhuang 050022, China; 2. University of Chinese Academy of Sciences, Beijing 100049, China)

    The North China Plain (NCP) is one of the most intensively farmed agricultural regions in China, with approximately 70% of the total cultivated land being used for an annual double-cropping system of winter wheat and summer maize. Owing to the long-term rotary and no tillage practices accompanying with the whole straw of winter wheat and summer maize return to field for several years, soil physical characteristics are gradually changing in terms of the increased soil pan depth, bulk density and content of soil nutrients in the surface soil layer. Improving soil quality by changing the tillage practices might help to maintain crop productivity in this region. An experiment was conducted for the winter wheat-summer maize rotation system in Luancheng County, which represented a high yield region, and in Nanpi County, which represented a medium and low yield region, in the NCP. Four treatments — soil deep tillage (DT), subsoiling (SS), narrow subsoiling (NSS), and rotary tillage (control, CK) — before winter wheat sowing and no tillage before summer maize sowing to all treatments were simultaneously conducted at the two areas. Soil bulk density, crop root growth, soil water use, yield and water use efficiency (WUE) were monitored throughout. Results showed that the effects of different tillage practices on soil and crop were different in the two regions. At Nanpi, deep tillage and subsoiling significantly increased crop yield. Compared with traditional rotary tillage, winter wheat yield was improved by 16.5% under DT, 19.3% under SS, and 13.1% under NSS. Yield of summer maize was increased by 17.3%, 16.2%, and 21.9%, respectively, with annual yield increases of 16.9%, 17.6% and 17.8%, respectively. Yield differences were not observed among the DT, SS, and NSS treatments. However, no significant difference in crop yield among the four treatments was found at Luancheng. Furthermore, four tillage practices reduced soil penetration resistance and bulk density for the 0-20-cm soil layer in both Luancheng and Nanpi. At Nanpi, after sowing winter wheat, the soil penetration resistance of the 0-20-cm soil layer under DT, SS, NSS and CK decreased by 69.7%, 72.7%, 72.5% and 68.2%, respectively. At Luancheng, soil penetration resistance of the 0-20-cm soil layer was reduced by 88.8% averagely under treatments of deep tillage and subsoiling, and slightly increased by 7.7% under CK. Soil bulk density of the 0-40-cm soil layer under the four tillage treatments were all lower at wheat harvest compared with that before tillage. Until the summer maize harvest, soil bulk density under different tillage treatments was essentially similar to that before tillage, and there was no significant difference among the four tillage treatments. At Nanpi, WUE of winter wheat and summer maize was significantly increased under DT, SS, and NSS compared with that under CK. At Luancheng, the WUE of winter wheat and summer maize was similar among the four treatments. These results indicated that different tillage practices in the low yield regions benefited crop production and water productivity. However, in the high yield regions, the three tillage practices did not enhance crop performance compared with traditional tillage practice. Therefore, it is suggested that the model of soil rotary tillage can be continuously implemented in the high yield regions of the NCP, whereas DT can be beneficially implemented in the medium and low yield regions.

    Winter wheat-summer maize rotation system; High-yield area; Medium and low-yield area; Tillage pattern; Soil bulk density; Yield composition; Water use efficiency

    關(guān)劼兮, 陳素英,邵立威, 張玉銘, 張喜英, 路楊, 閆宗正. 華北典型區(qū)域土壤耕作方式對(duì)土壤特性和作物產(chǎn)量的影響[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)(中英文), 2019, 27(11): 1663-1672

    GUAN J X, CHEN S Y, SHAO L W, ZHANG Y M, ZHANG X Y, LU Y, YAN Z Z. Soil tillage practices affecting the soil characteristics and yield of winter wheat and summer maize in North China[J]. Chinese Journal of Eco-Agriculture, 2019, 27(11): 1663-1672

    * 國(guó)家重點(diǎn)研發(fā)計(jì)劃專項(xiàng)(2016YFD0300808, 2016YFC0401403)和國(guó)家自然科學(xué)基金項(xiàng)目(31371578)資助

    陳素英, 主要從事農(nóng)田節(jié)水和保護(hù)性耕作技術(shù)研究。E-mail: csy@sjziam.ac.cn

    關(guān)劼兮, 主要從事農(nóng)田節(jié)水機(jī)理與技術(shù)研究。E-mail: gjx391214674@qq.com

    2019-04-02

    2019-06-04

    * This study was supported by the National Key R&D Program of China (2016YFD0300808, 2016YFC0401403) and the National Natural Science Foundation of China (31371578).

    , E-mail: csy@sjziam.ac.cn

    Jun. 4, 2019

    Apr. 2, 2019;

    S342

    2096-6237(2019)11-1663-10

    10.13930/j.cnki.cjea.190246

    猜你喜歡
    欒城實(shí)度夏玉米
    “小家站、大民生”繪出欒城高質(zhì)量發(fā)展新畫卷——河北省石家莊市欒城區(qū)人大常委會(huì)以家站為依托搞活代表工作
    深入學(xué)習(xí)貫徹習(xí)近平法治思想 奮力打造公共法律服務(wù)欒城品牌
    石家莊市欒城區(qū)人大常委會(huì)組織“百名代表”視察民生實(shí)事辦理工作
    公民與法治(2022年1期)2022-07-26 05:58:22
    石家莊市欒城區(qū)人大常委會(huì)高度重視少先隊(duì)員提案工作
    公民與法治(2022年5期)2022-02-05 02:43:40
    土壤緊實(shí)度對(duì)毛竹冬筍和春筍產(chǎn)量的影響
    小麥?zhǔn)崭钪?如何種植夏玉米才能高產(chǎn)
    土壤緊實(shí)度對(duì)裸土侵蝕強(qiáng)度影響的實(shí)驗(yàn)與分析
    夏玉米高產(chǎn)的關(guān)鍵栽培技術(shù)措施
    基于加速度補(bǔ)償?shù)耐寥谰o實(shí)度測(cè)量方法與傳感器設(shè)計(jì)
    內(nèi)陸沙區(qū)丘間地土壤緊實(shí)度特征及其與植被的關(guān)系
    乱码一卡2卡4卡精品| 看免费av毛片| 久久草成人影院| 亚洲精品456在线播放app | 午夜激情福利司机影院| 成人午夜高清在线视频| 能在线免费观看的黄片| 亚洲成a人片在线一区二区| 人妻久久中文字幕网| 极品教师在线视频| 亚洲欧美日韩东京热| 欧美色欧美亚洲另类二区| aaaaa片日本免费| 直男gayav资源| 久久精品国产亚洲av天美| 亚洲av电影不卡..在线观看| 国产精品,欧美在线| 亚洲成人久久性| 国产精品一区二区三区四区免费观看 | 亚洲熟妇中文字幕五十中出| 在线播放无遮挡| 久久午夜福利片| 波多野结衣高清作品| 成年版毛片免费区| 99国产综合亚洲精品| 性色av乱码一区二区三区2| 国产av不卡久久| 一区二区三区四区激情视频 | 性欧美人与动物交配| 日本在线视频免费播放| 国产 一区 欧美 日韩| 色5月婷婷丁香| 国产精品综合久久久久久久免费| 久久亚洲精品不卡| 国产三级黄色录像| 久久久久九九精品影院| 欧美日韩亚洲国产一区二区在线观看| 啦啦啦观看免费观看视频高清| 伦理电影大哥的女人| 国产精品影院久久| 波多野结衣高清作品| 天天躁日日操中文字幕| 色视频www国产| 国产精品野战在线观看| 精品国产亚洲在线| 99久久精品热视频| 最近中文字幕高清免费大全6 | 岛国在线免费视频观看| 午夜激情福利司机影院| 亚洲av成人精品一区久久| 国产探花极品一区二区| 99久久久亚洲精品蜜臀av| 国内揄拍国产精品人妻在线| av在线天堂中文字幕| 国产欧美日韩一区二区精品| 一卡2卡三卡四卡精品乱码亚洲| 可以在线观看毛片的网站| 亚洲经典国产精华液单 | 蜜桃亚洲精品一区二区三区| 亚洲av中文字字幕乱码综合| 日本成人三级电影网站| 午夜福利高清视频| 久久精品国产自在天天线| 国产一区二区三区在线臀色熟女| 国产欧美日韩一区二区精品| 又爽又黄a免费视频| 亚洲av成人精品一区久久| 日日摸夜夜添夜夜添小说| .国产精品久久| 国产午夜精品论理片| 狠狠狠狠99中文字幕| 男女之事视频高清在线观看| 十八禁人妻一区二区| 精品午夜福利视频在线观看一区| 窝窝影院91人妻| 99热这里只有是精品50| 小蜜桃在线观看免费完整版高清| 国产乱人伦免费视频| 人人妻,人人澡人人爽秒播| 亚洲精品亚洲一区二区| 波多野结衣高清无吗| 国产真实伦视频高清在线观看 | 国产三级黄色录像| 国产午夜精品久久久久久一区二区三区 | 国产精品国产高清国产av| 在线十欧美十亚洲十日本专区| 国产白丝娇喘喷水9色精品| 亚洲中文字幕一区二区三区有码在线看| 最好的美女福利视频网| 欧美午夜高清在线| 国产一区二区三区在线臀色熟女| 男女视频在线观看网站免费| 美女 人体艺术 gogo| 国产单亲对白刺激| 国产aⅴ精品一区二区三区波| 国产在线男女| 在线天堂最新版资源| 久久亚洲精品不卡| 亚洲精品粉嫩美女一区| 欧美最新免费一区二区三区 | 久久国产乱子免费精品| 成人av一区二区三区在线看| ponron亚洲| 18禁黄网站禁片免费观看直播| 国产亚洲欧美在线一区二区| 国产精品三级大全| 有码 亚洲区| 俄罗斯特黄特色一大片| av视频在线观看入口| 在现免费观看毛片| 男人的好看免费观看在线视频| 欧美成狂野欧美在线观看| xxxwww97欧美| 尤物成人国产欧美一区二区三区| 少妇的逼水好多| 亚洲av二区三区四区| 国产乱人视频| 久久这里只有精品中国| 国产一区二区亚洲精品在线观看| 精品久久久久久久末码| 很黄的视频免费| 午夜久久久久精精品| 特级一级黄色大片| 久久久久久九九精品二区国产| 婷婷六月久久综合丁香| 欧美性猛交╳xxx乱大交人| 欧美成狂野欧美在线观看| 级片在线观看| 亚洲av成人不卡在线观看播放网| 欧美黑人巨大hd| 亚洲人成网站在线播| 精品不卡国产一区二区三区| 亚洲国产欧美人成| 校园春色视频在线观看| 国内少妇人妻偷人精品xxx网站| 黄色视频,在线免费观看| 欧美极品一区二区三区四区| 国产精品综合久久久久久久免费| 成人午夜高清在线视频| av专区在线播放| 2021天堂中文幕一二区在线观| 亚洲七黄色美女视频| 色尼玛亚洲综合影院| 国产精品久久久久久久电影| 国产精品爽爽va在线观看网站| 久久精品国产99精品国产亚洲性色| 99久久九九国产精品国产免费| 国产精品一及| 午夜激情欧美在线| 九九在线视频观看精品| a级毛片a级免费在线| 成年女人毛片免费观看观看9| h日本视频在线播放| 亚洲精品色激情综合| 熟女人妻精品中文字幕| 少妇高潮的动态图| 亚洲自偷自拍三级| 欧美黄色淫秽网站| 亚洲国产高清在线一区二区三| 色哟哟·www| 国产精品伦人一区二区| 三级国产精品欧美在线观看| а√天堂www在线а√下载| 色尼玛亚洲综合影院| 国产精品一区二区性色av| 亚洲熟妇中文字幕五十中出| 成熟少妇高潮喷水视频| 国产极品精品免费视频能看的| 床上黄色一级片| 国产单亲对白刺激| 欧美一区二区精品小视频在线| 制服丝袜大香蕉在线| 国产久久久一区二区三区| 国产一区二区三区视频了| 欧美性猛交╳xxx乱大交人| 国产久久久一区二区三区| 99久久久亚洲精品蜜臀av| 一个人免费在线观看电影| av黄色大香蕉| 欧美一区二区精品小视频在线| 亚洲av一区综合| 别揉我奶头 嗯啊视频| 免费看光身美女| 欧美在线黄色| 三级毛片av免费| 成人国产综合亚洲| 丰满乱子伦码专区| 中文在线观看免费www的网站| 91字幕亚洲| 久久亚洲真实| 亚洲无线观看免费| 久久精品国产清高在天天线| 国产精品亚洲av一区麻豆| 国产老妇女一区| 露出奶头的视频| 日本三级黄在线观看| 国产高清视频在线播放一区| 丝袜美腿在线中文| 中文资源天堂在线| 高清日韩中文字幕在线| 狂野欧美白嫩少妇大欣赏| 欧美不卡视频在线免费观看| 亚洲熟妇熟女久久| 亚洲av电影不卡..在线观看| 99国产极品粉嫩在线观看| 亚洲精品影视一区二区三区av| 成人欧美大片| 国产91精品成人一区二区三区| x7x7x7水蜜桃| 国产精品一区二区性色av| 搡老妇女老女人老熟妇| 天堂av国产一区二区熟女人妻| www日本黄色视频网| 精品欧美国产一区二区三| 欧美+日韩+精品| 在线免费观看的www视频| 两人在一起打扑克的视频| 在线a可以看的网站| 欧美日韩中文字幕国产精品一区二区三区| 国产又黄又爽又无遮挡在线| 美女大奶头视频| 麻豆国产av国片精品| 在线观看舔阴道视频| 久久精品国产清高在天天线| 久久草成人影院| 身体一侧抽搐| 午夜福利在线观看吧| 午夜视频国产福利| 免费av毛片视频| 老司机午夜十八禁免费视频| 两个人的视频大全免费| 欧美zozozo另类| 99久久精品一区二区三区| 精品久久久久久久久亚洲 | 精品99又大又爽又粗少妇毛片 | 亚洲国产欧美人成| 国产精品一区二区免费欧美| 亚洲美女黄片视频| 亚洲国产色片| 亚洲色图av天堂| 国产亚洲欧美在线一区二区| 99热这里只有是精品50| 热99在线观看视频| 国产精品亚洲av一区麻豆| 人妻夜夜爽99麻豆av| 色播亚洲综合网| 日韩欧美在线乱码| 99热只有精品国产| 色尼玛亚洲综合影院| 最好的美女福利视频网| 精品熟女少妇八av免费久了| 国产午夜福利久久久久久| 一个人免费在线观看电影| 国产探花在线观看一区二区| 在线观看一区二区三区| 女人被狂操c到高潮| 麻豆久久精品国产亚洲av| 毛片一级片免费看久久久久 | 亚洲精品一卡2卡三卡4卡5卡| 久久久久久大精品| 国产色爽女视频免费观看| 亚洲欧美日韩无卡精品| 免费在线观看成人毛片| 成人av一区二区三区在线看| 午夜福利视频1000在线观看| 欧美日韩瑟瑟在线播放| 一个人免费在线观看电影| 国产精品野战在线观看| 国产av不卡久久| 在线播放国产精品三级| avwww免费| 免费在线观看亚洲国产| 国产极品精品免费视频能看的| 久久久久久久精品吃奶| 一个人看视频在线观看www免费| 亚洲熟妇中文字幕五十中出| 最近最新中文字幕大全电影3| 精华霜和精华液先用哪个| 亚洲欧美日韩高清专用| 国产亚洲精品久久久久久毛片| 国产精品女同一区二区软件 | 12—13女人毛片做爰片一| 亚洲五月天丁香| 91狼人影院| 床上黄色一级片| 国产又黄又爽又无遮挡在线| 国产精品久久久久久精品电影| 性欧美人与动物交配| 一进一出抽搐gif免费好疼| 亚洲精华国产精华精| 亚洲,欧美,日韩| 亚洲av日韩精品久久久久久密| 国产精品一区二区免费欧美| 亚洲精品日韩av片在线观看| 亚洲一区高清亚洲精品| 国产野战对白在线观看| 非洲黑人性xxxx精品又粗又长| 全区人妻精品视频| 亚洲五月天丁香| 国产精品乱码一区二三区的特点| 天美传媒精品一区二区| 久久久久久久久久成人| 琪琪午夜伦伦电影理论片6080| 国产在线精品亚洲第一网站| 亚洲美女视频黄频| 日本三级黄在线观看| 日本五十路高清| 能在线免费观看的黄片| 国产爱豆传媒在线观看| av女优亚洲男人天堂| 中文亚洲av片在线观看爽| 一个人观看的视频www高清免费观看| 久久久久久久久久成人| 欧美性猛交╳xxx乱大交人| 亚洲三级黄色毛片| 国产成年人精品一区二区| 尤物成人国产欧美一区二区三区| 久久国产乱子免费精品| 成人鲁丝片一二三区免费| 久久精品人妻少妇| 欧美成狂野欧美在线观看| 欧美日韩综合久久久久久 | 国产精品嫩草影院av在线观看 | 我的老师免费观看完整版| 18美女黄网站色大片免费观看| 亚洲va日本ⅴa欧美va伊人久久| 淫秽高清视频在线观看| 国产激情偷乱视频一区二区| 高清日韩中文字幕在线| 真实男女啪啪啪动态图| 欧美黄色片欧美黄色片| av在线老鸭窝| 国产成年人精品一区二区| 天堂影院成人在线观看| 九九久久精品国产亚洲av麻豆| 精品午夜福利视频在线观看一区| 国产野战对白在线观看| 欧美黄色淫秽网站| 美女免费视频网站| 成年女人毛片免费观看观看9| 人人妻人人澡欧美一区二区| 日本黄色片子视频| 午夜福利18| 在线免费观看不下载黄p国产 | 亚洲精品粉嫩美女一区| 脱女人内裤的视频| 天堂√8在线中文| 精品久久久久久久久久免费视频| 免费高清视频大片| 老女人水多毛片| 亚洲av.av天堂| 久久伊人香网站| 久久精品国产自在天天线| 亚洲第一欧美日韩一区二区三区| 日韩欧美精品v在线| 亚洲欧美激情综合另类| 亚洲av成人不卡在线观看播放网| 国产视频内射| 色在线成人网| 成年免费大片在线观看| 免费看a级黄色片| xxxwww97欧美| 精品久久久久久,| 成年免费大片在线观看| 一级a爱片免费观看的视频| 国产亚洲欧美98| 村上凉子中文字幕在线| 午夜精品久久久久久毛片777| 永久网站在线| 最近在线观看免费完整版| 亚洲第一区二区三区不卡| 久久九九热精品免费| 欧美日韩黄片免| 欧美乱妇无乱码| 熟女电影av网| 欧美成狂野欧美在线观看| 国内精品美女久久久久久| 禁无遮挡网站| 国产一级毛片七仙女欲春2| 给我免费播放毛片高清在线观看| 999久久久精品免费观看国产| 国产亚洲av嫩草精品影院| 日韩中文字幕欧美一区二区| 久久久色成人| 日韩精品青青久久久久久| 非洲黑人性xxxx精品又粗又长| 久久久久亚洲av毛片大全| 久久久久性生活片| 国产野战对白在线观看| 精品午夜福利视频在线观看一区| 亚洲片人在线观看| 乱码一卡2卡4卡精品| 在线免费观看的www视频| 国产成人a区在线观看| 国产高清视频在线观看网站| 亚洲av中文字字幕乱码综合| 免费大片18禁| 少妇高潮的动态图| 一本久久中文字幕| 成人毛片a级毛片在线播放| 日韩亚洲欧美综合| 久久热精品热| 欧美黄色淫秽网站| 国产主播在线观看一区二区| 老女人水多毛片| 日本精品一区二区三区蜜桃| 国产av麻豆久久久久久久| 亚洲无线在线观看| 久久精品国产亚洲av涩爱 | h日本视频在线播放| 国模一区二区三区四区视频| 国产极品精品免费视频能看的| 亚洲成av人片免费观看| 日本黄色视频三级网站网址| 国产精品嫩草影院av在线观看 | 精华霜和精华液先用哪个| 精品久久久久久成人av| av在线蜜桃| 日韩 亚洲 欧美在线| 婷婷精品国产亚洲av在线| 女人被狂操c到高潮| 久久精品久久久久久噜噜老黄 | 国产视频一区二区在线看| 男人的好看免费观看在线视频| 欧美日韩黄片免| 可以在线观看的亚洲视频| 久久久久亚洲av毛片大全| 又爽又黄无遮挡网站| 国产av在哪里看| 中出人妻视频一区二区| 丰满的人妻完整版| 怎么达到女性高潮| 噜噜噜噜噜久久久久久91| 亚洲午夜理论影院| 亚洲七黄色美女视频| 白带黄色成豆腐渣| 99riav亚洲国产免费| av在线蜜桃| 日韩 亚洲 欧美在线| 午夜免费成人在线视频| 国产麻豆成人av免费视频| 亚洲精品亚洲一区二区| 最新在线观看一区二区三区| 91av网一区二区| 国产av一区在线观看免费| 五月玫瑰六月丁香| 亚洲国产日韩欧美精品在线观看| 性色av乱码一区二区三区2| 免费无遮挡裸体视频| 亚洲五月婷婷丁香| 波多野结衣高清无吗| 国产69精品久久久久777片| 亚洲精品成人久久久久久| 男人和女人高潮做爰伦理| 亚洲无线观看免费| 最新中文字幕久久久久| 日本一本二区三区精品| 久久久久久久午夜电影| 日本熟妇午夜| 国产精华一区二区三区| 男人舔奶头视频| 婷婷丁香在线五月| 亚洲av免费在线观看| 麻豆一二三区av精品| 听说在线观看完整版免费高清| 搞女人的毛片| 一本久久中文字幕| 国产精品久久久久久亚洲av鲁大| 色av中文字幕| 波多野结衣高清作品| 9191精品国产免费久久| 天天躁日日操中文字幕| 日韩欧美 国产精品| 国产69精品久久久久777片| 亚洲美女视频黄频| 欧美潮喷喷水| 欧美成人性av电影在线观看| 国产精品久久久久久久电影| 免费在线观看成人毛片| 久久精品国产亚洲av香蕉五月| 国内精品久久久久精免费| 欧美性猛交╳xxx乱大交人| 亚洲色图av天堂| 国产精品嫩草影院av在线观看 | 色噜噜av男人的天堂激情| 久久国产精品人妻蜜桃| 色哟哟·www| 亚洲va日本ⅴa欧美va伊人久久| 99riav亚洲国产免费| 又粗又爽又猛毛片免费看| 国产成人欧美在线观看| 宅男免费午夜| 国产亚洲精品av在线| 亚洲精华国产精华精| 国产高清视频在线播放一区| 91麻豆av在线| 日本精品一区二区三区蜜桃| 国产成+人综合+亚洲专区| 欧美xxxx黑人xx丫x性爽| 午夜福利高清视频| 99久国产av精品| 在线免费观看的www视频| 国产亚洲精品综合一区在线观看| 亚洲五月天丁香| 欧美一区二区精品小视频在线| 蜜桃亚洲精品一区二区三区| 日本精品一区二区三区蜜桃| 久久久久久久久久成人| 国产亚洲精品久久久com| 免费看a级黄色片| 99久久成人亚洲精品观看| 亚洲精品粉嫩美女一区| 成人国产综合亚洲| 能在线免费观看的黄片| av欧美777| 欧美潮喷喷水| 色播亚洲综合网| 老熟妇仑乱视频hdxx| 中出人妻视频一区二区| 一二三四社区在线视频社区8| 欧美绝顶高潮抽搐喷水| 亚洲av免费高清在线观看| 色在线成人网| 精品久久久久久久久久免费视频| 亚洲av一区综合| 日本黄色视频三级网站网址| 欧美3d第一页| 高清日韩中文字幕在线| 久久久久免费精品人妻一区二区| 久久久久国内视频| 亚洲真实伦在线观看| 精品久久久久久久久久久久久| 综合色av麻豆| 在线a可以看的网站| 免费在线观看亚洲国产| 国产91精品成人一区二区三区| 韩国av一区二区三区四区| 精品久久久久久久末码| 婷婷精品国产亚洲av在线| av黄色大香蕉| 欧美性感艳星| 午夜福利在线观看免费完整高清在 | 精品福利观看| 午夜激情欧美在线| 久久久久久久久大av| 91九色精品人成在线观看| 97热精品久久久久久| av中文乱码字幕在线| 九色国产91popny在线| 超碰av人人做人人爽久久| 国产极品精品免费视频能看的| 国产精品伦人一区二区| 99在线人妻在线中文字幕| 少妇高潮的动态图| 热99re8久久精品国产| 九九久久精品国产亚洲av麻豆| 日韩欧美在线乱码| 嫁个100分男人电影在线观看| 色在线成人网| 啦啦啦观看免费观看视频高清| 女人十人毛片免费观看3o分钟| ponron亚洲| 日本在线视频免费播放| 亚洲性夜色夜夜综合| 全区人妻精品视频| 老司机午夜福利在线观看视频| АⅤ资源中文在线天堂| 国语自产精品视频在线第100页| 非洲黑人性xxxx精品又粗又长| 亚洲激情在线av| 久久久久久久久久成人| 国产淫片久久久久久久久 | 在线看三级毛片| 欧美黄色片欧美黄色片| 少妇熟女aⅴ在线视频| 精品人妻熟女av久视频| 黄色女人牲交| 一本综合久久免费| 亚洲欧美日韩高清专用| 国产不卡一卡二| 小蜜桃在线观看免费完整版高清| 中国美女看黄片| 日韩欧美精品v在线| 啦啦啦韩国在线观看视频| 亚洲精品影视一区二区三区av| 亚洲精品亚洲一区二区| 真实男女啪啪啪动态图| 欧美成人性av电影在线观看| 久久精品人妻少妇| 久久人妻av系列| 国产高清有码在线观看视频| 91字幕亚洲| 久久人妻av系列| 两个人视频免费观看高清| 又粗又爽又猛毛片免费看| 观看免费一级毛片| 女人十人毛片免费观看3o分钟| 男人舔奶头视频| 亚洲人成网站高清观看| 亚洲国产精品成人综合色| 国产高潮美女av| 天堂影院成人在线观看| 日本免费a在线| 日韩人妻高清精品专区| 中文字幕人妻熟人妻熟丝袜美| 最近中文字幕高清免费大全6 | 在线观看一区二区三区| 91在线观看av| 99在线视频只有这里精品首页| 日本 av在线| 国产亚洲精品av在线| 1024手机看黄色片| 婷婷六月久久综合丁香| 特级一级黄色大片| 亚洲人与动物交配视频|