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

    灌水量對(duì)冬小麥農(nóng)田土壤N2O 與CO2 排放的影響

    2020-09-28 02:03:42張保成李建柱石建華
    灌溉排水學(xué)報(bào) 2020年9期

    熊 浩,張保成,李建柱,張 旭,石建華

    (1.天津大學(xué) 水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室,天津 300350; 2.北京中水科苑工程咨詢(xún)有限公司,北京 100053)

    0 引 言

    【研究意義】全球變暖引發(fā)的極端氣候事件已經(jīng)成為人類(lèi)面臨的重要挑戰(zhàn)之一。人類(lèi)活動(dòng)引起的溫室氣體排放作為導(dǎo)致全球變暖的關(guān)鍵因素已成為相關(guān)領(lǐng)域普遍關(guān)注的焦點(diǎn)[1-2]。氧化亞氮(N2O)、二氧化碳(CO2)是2 種主要的溫室氣體,農(nóng)田土壤是其重要的排放源[3],據(jù)估計(jì),大氣中每年80%~90%的N2O和5%~20%的CO2來(lái)源于土壤[4],因此,研究農(nóng)田土壤N2O 和CO2的排放規(guī)律,在溫室氣體排放研究中具有十分重要的意義。

    【研究進(jìn)展】目前對(duì)農(nóng)田溫室氣體排放的研究主要集中在施肥、秸稈還田和耕作方式等方面[5-7]。李燕青等[8]發(fā)現(xiàn)施肥會(huì)顯著增加農(nóng)田土壤N2O 和CO2的排放總量。趙力瑩等[9]的研究發(fā)現(xiàn)冬小麥季農(nóng)田在長(zhǎng)期旋耕基礎(chǔ)上經(jīng)深松和翻耕后,N2O 和CO2排放有顯著變化。土壤水分是溫室氣體排放的關(guān)鍵驅(qū)動(dòng)因子[10],而人類(lèi)活動(dòng)中,灌溉是改變農(nóng)田土壤水分的主要措施之一。灌溉通過(guò)改變土壤水分狀況進(jìn)而影響農(nóng)田土壤N2O 和CO2排放,同時(shí)充足的灌水量也是作物保持高產(chǎn)的關(guān)鍵。灌溉量的多少、灌溉方式和灌溉的頻率等都會(huì)對(duì)土壤N2O 和CO2的排放產(chǎn)生明顯影響[11]?!厩腥朦c(diǎn)】目前研究灌溉對(duì)溫室氣體排放的影響主要在稻田、菜地和節(jié)水灌溉措施等方面[12-16]關(guān)于灌溉用水量對(duì)北方麥田溫室氣體排放的影響研究甚少。

    王玉英等[17]對(duì)太行山前平原冬小麥-夏玉米輪作農(nóng)田的試驗(yàn)研究表明施肥和灌溉會(huì)直接導(dǎo)致土壤N2O 和CO2排放通量增大,溫度升高和干濕交替都會(huì)增大土壤N2O 和CO2的排放速率。郭樹(shù)芳等[18]在華北平原西部的試驗(yàn)表明,溫度和土壤水分是影響土壤N2O 及CO2排放通量的主要因素,且微噴灌溉下土壤N2O 和CO2排放均高于傳統(tǒng)漫灌。劉晶晶等[19]的研究發(fā)現(xiàn)土壤呼吸與灌溉量正相關(guān),但過(guò)量灌溉會(huì)降低土壤呼吸,在一定范圍內(nèi)灌溉量的增加增大了土壤N2O 的年際排放總量。

    【擬解決的關(guān)鍵問(wèn)題】本文以華北地區(qū)冬小麥農(nóng)田土壤為研究對(duì)象,對(duì)不同灌溉水平下冬小麥全生育期土壤N2O 和CO2排放進(jìn)行田間監(jiān)測(cè),分析土壤N2O和CO2的排放特征和綜合溫室效應(yīng),結(jié)合不同灌溉水平下小麥產(chǎn)量,為制定節(jié)水、減排的灌水定額提供參考。

    1 材料與方法

    1.1 研究區(qū)概況

    選擇河北省石家莊市晉州市周家莊鄉(xiāng)科學(xué)灌溉項(xiàng)目農(nóng)田試驗(yàn)站作為試驗(yàn)區(qū),該試驗(yàn)區(qū)位于河北省中南部,滹沱河流域,屬東部季風(fēng)氣候區(qū),主要種植制度為冬小麥-夏玉米輪作,一年二熟。區(qū)域年平均氣溫6.8~10.5 ℃,多年平均降雨量506 mm,全年降水主要集中在7、8、9 月,雨熱同期。供試土壤為輕壤質(zhì)土,試驗(yàn)開(kāi)始前0~20 cm 土層土壤基本理化性質(zhì)為:土壤體積質(zhì)量1.26 g/cm3,有機(jī)質(zhì)量為14.02 g/kg,速效磷量為11 mg/kg,速效鉀量為85 mg/kg。

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

    以冬小麥農(nóng)田土壤為研究對(duì)象,按灌水單因素設(shè)計(jì)試驗(yàn),共設(shè)置6個(gè)灌水水平,分別為:CK(89.96 mm)、T1(80.96 mm,減少10%)、T2(71.97 mm,減少20%)、T3(62.97 mm,減少30%)、T4(53.98 mm,減少40%)、T5(44.98 mm,減少50%),CK 是按當(dāng)?shù)毓嗨~換算的灌水量,6 個(gè)灌水量處理,每個(gè)處理2 次重復(fù),共計(jì)12 個(gè)試驗(yàn)區(qū),小區(qū)面積120 m2(15 m×8 m),每個(gè)小區(qū)之間設(shè)1.5 m 寬隔離帶,試驗(yàn)小區(qū)周?chē)O(shè)置保護(hù)區(qū),保護(hù)區(qū)處理與相鄰試驗(yàn)小區(qū)一致。前茬作物玉米收獲后,旋耕機(jī)旋地2 次,深度20 cm,肥料為當(dāng)?shù)仄毡槭褂玫膹?fù)合肥(N-P2O5-K2O 18-17-10),播種同時(shí)按750 kg/hm2施基肥,拔節(jié)期按322.5 kg/hm2追施尿素。供試小麥品種為“冀麥418”,于2018 年10月17 日播種,2019 年6 月10 日收獲,大田生長(zhǎng)期為236 d。主要生育期劃分為苗期(10 月17 日—12月10 日)、越冬期(12 月11 日—3 月5 日)、返青—拔節(jié)期階段(3 月6 日—4 月30 日)、拔節(jié)—灌漿期段(5 月1 日—5 月26 日)、成熟期(5 月27 日—6月10 日)。為保證冬小麥順利出苗,試驗(yàn)各小區(qū)在播種前均灌水44.98 mm,整個(gè)生育期內(nèi)于11 月28 日、4 月16 日和5 月18 日共灌水3 次,4 月16 日灌水的同時(shí)追施尿素,其他田間管理措施與當(dāng)?shù)爻R?guī)管理保持一致。

    1.3 田間采樣與觀測(cè)

    試驗(yàn)用靜態(tài)箱法采集氣體(N2O、CO2)樣品。靜態(tài)箱底座于小麥播種后當(dāng)天埋設(shè)在各試驗(yàn)小區(qū)中部,不再移動(dòng)直至小麥?zhǔn)斋@,底座內(nèi)無(wú)小麥植株,底座上端有4 cm 深的凹槽用以放置靜態(tài)箱體,取樣時(shí)往凹槽內(nèi)注水密封防止周?chē)諝馀c箱內(nèi)氣體交換。箱體由板厚1 mm 的不銹鋼板制成,幾何尺寸為50 cm×50 cm×100 cm,箱體外表面包有隔熱泡沫板,防止取樣期間陽(yáng)光照射導(dǎo)致箱體內(nèi)溫度的劇烈升高,箱體內(nèi)壁裝有小風(fēng)扇,每次取樣前風(fēng)扇轉(zhuǎn)動(dòng)以使箱內(nèi)氣體混合均勻。箱體側(cè)部有溫度計(jì)插孔和取氣口(閥門(mén)),利用采氣泵連接取氣閥門(mén)抽出箱內(nèi)氣體并注入鋁箔真空采樣袋內(nèi),采集氣體的同時(shí)讀取安插在箱體側(cè)部的水銀溫度計(jì)讀數(shù),用于計(jì)算氣體排放通量。

    小麥播種后第7 天開(kāi)始田間原位氣體采樣,在冬小麥全生育期內(nèi)各個(gè)關(guān)鍵時(shí)期進(jìn)行氣體樣品的采集,平均7~10 d 取樣1 次。在有效降雨(大于5 mm)或施肥灌水后增加采樣頻率,分別在第1、3、5 天取樣,整個(gè)小麥季共進(jìn)行了28 次取樣。取樣時(shí)間在上午10:00—11:30,此時(shí)間段氣體排放通量最接近1 d 的平均排放通量[20],為保證充足的準(zhǔn)備時(shí)間,封箱(09:30)后30 min 開(kāi)始第1 次取樣,之后每隔20 min取樣1 次,每個(gè)采樣點(diǎn)共取樣4 次,已驗(yàn)證按此法采樣的氣體通量計(jì)算結(jié)果在合理誤差范圍內(nèi),該方法是可行的。采樣完成后將采樣袋依次分組并標(biāo)號(hào),帶回實(shí)驗(yàn)室進(jìn)行分析。

    1.4 指標(biāo)測(cè)定方法

    1.4.1 氣體排放通量測(cè)定

    氣體樣品采用Agilent hp5890 氣相色譜分析儀進(jìn)行測(cè)定分析,通過(guò)對(duì)每組4 個(gè)樣品的目標(biāo)氣體混合比與相對(duì)應(yīng)的采樣時(shí)間間隔0、20、40、60 min 進(jìn)行線性擬合,當(dāng)回歸系數(shù)R2>0.75 時(shí),視為有效數(shù)據(jù),求得目標(biāo)氣體的排放通量。

    CO2、N2O 排放通量(F)的計(jì)算式為:

    式中:F 為氣體的排放通量mg/(m2·h),負(fù)值表示土壤吸收該氣體,正值表示土壤排放該氣體。ρ 為標(biāo)準(zhǔn)狀態(tài)下的氣體密度(g/cm3);h 為采樣箱高度(1.0 m);T 為采樣時(shí)的箱內(nèi)溫度(℃);dc/dt 為箱內(nèi)氣體濃度的變化率。

    1.4.2 氣體排放總量和綜合溫室效應(yīng)計(jì)算

    CO2、N2O 累積排放量計(jì)算式為:

    式中:M 為全生育期內(nèi)氣體的累積排放量(kg/hm2);F 為氣體排放通量(mg/(m2·h));i 為采樣次數(shù);t為采樣時(shí)間(d)。

    對(duì)于100 年時(shí)間尺度的氣候變化,N2O 氣體相對(duì)于CO2的GWPs(全球增溫潛勢(shì))為298[2],將N2O 的排放量乘以298,可得出N2O 排放量的CO2當(dāng)量 (E-CO2),然后與CO2排放量相加,即得小麥季農(nóng)田CO2和N2O 排放的綜合溫室效應(yīng)(總CO2當(dāng)量)[21]。

    1.4.3 土壤溫度和水分測(cè)定

    土壤溫度測(cè)定:在采集氣體樣品的同時(shí),讀取埋設(shè)在各試驗(yàn)區(qū)采樣箱附近的地溫計(jì)讀數(shù),并記錄。

    土壤水分測(cè)定:每次采集氣體樣品后,在采樣箱底座周邊進(jìn)行土樣的采集。分別取0~10、10~20 cm土層土壤,用烘干法在(105±2)℃測(cè)定各土層土壤含水率,并計(jì)算土壤孔隙度含水率(WFPS)。

    1.4.4 小麥產(chǎn)量測(cè)定

    小麥?zhǔn)斋@時(shí),在每個(gè)試驗(yàn)區(qū)選取2 個(gè)1 m2的樣點(diǎn)進(jìn)行單獨(dú)收割,風(fēng)干脫粒后稱(chēng)取籽粒質(zhì)量并計(jì)算產(chǎn)量,最終結(jié)果以單位面積產(chǎn)量(kg/hm2)表示。

    2 結(jié)果與分析

    2.1 不同灌溉水平下土壤N2O 排放特征

    2.1.1 不同灌溉水平下N2O 排放通量變化

    圖1 為各灌水量處理下冬小麥土壤N2O 排放通量動(dòng)態(tài)變化圖。由圖1 可以看出,冬小麥整個(gè)生育期內(nèi)各處理土壤N2O 排放通量的變化規(guī)律基本一致。在苗期(2018 年10 月24 日左右)各處理N2O 排放出現(xiàn)峰值,這與播種前旋耕土地和施肥有關(guān),其中峰值最大的是T2 處理:872.29 μg/(m2·h),峰值最小的是T5 處理:530.43 μg/(m2·h)。隨后N2O 的排放通量快速降低,這可能與土壤肥力逐漸減小有關(guān),播種后38 d(2018 年11 月25 日)左右,各處理土壤N2O排放通量降至較低水平。11月26日—4月10日,各處理N2O 的排放通量都在較低水平范圍內(nèi)波動(dòng),波動(dòng)范圍為-55.07~109.05 μg/(m2·h)。4 月16 日施肥灌水后,各處理N2O 排放出現(xiàn)峰值,其中峰值最大的是T2 處理:120.68 μg/(m2·h),最小的是T5 處理:98.82 μg/(m2·h)。5 月18 日灌水后,各處理N2O 排放通量略有增大但變化不明顯。

    2.1.2 不同灌溉水平下冬小麥各生育期內(nèi)土壤N2O 排放量

    不同灌水量處理下冬小麥全生育期內(nèi)土壤N2O的累積排放量如圖2 所示。

    由圖2 可知,灌水量的減少能顯著降低小麥生育期土壤N2O 的排放強(qiáng)度,CK 全生育期土壤N2O 累積排放量最大,T1 和T2 處理次之,T5 處理N2O 累積排放量最小。T1 處理較CK 減少了24.99%(P=0.18);T2 處理較CK 減少了10.89%(P=0.53);T3 處理較CK 減少了39.14%(P=0.06);T4 處理較CK 減少了48.99%,差異顯著(P<0.05);T5 處理較CK 減少了56.11%,差異顯著(P<0.05)。

    圖1 不同灌水量下土壤N2O 排放通量的變化 Fig.1 Temporal variation of N2O emission flux at different irrigation schemes

    圖2 不同灌水量下土壤N2O 累積排放量 Fig.2 Cumulative N2O emission of soil under stageat different irrigation schemes

    冬小麥全生育期,各灌水量處理下土壤N2O 的平均排放通量表現(xiàn)為CK>T1、T2 處理>T3、T4 處理>T5 處理(表1)。T4 處理比CK 減小48.30%,差異顯著(P<0.05);T5 處理比CK 減小55.41%,差異顯著(P<0.05);T5 處理比T2 處理減小50.65%,差異顯著(P<0.05);其他處理之間無(wú)顯著性差異(P>0.05)。

    冬小麥不同生育階段內(nèi),各灌水量處理下土壤N2O 平均排放通量與整個(gè)生育期的變化規(guī)律基本一致。苗期,T1、T2 處理較CK 分別減小34.50%(P<0.05)和18.22%(P>0.05),T3、T4 處理較CK 分別減小58.37%和48.97%(P<0.05),T5 處理較CK 顯著減小64.43%(P<0.01)。越冬期,T1—T5 處理較CK 分別減小32.61%、39.87%、31.59%、59.44%(P<0.05)和36.79%,差異不顯著。返青至拔節(jié)期,T1—T3 處理土壤N2O 平均排放通量較CK 分別增大28.72 μg/(m2·h)和32.72 μg/(m2·h),T4 和T5 處理較CK 分別減小4.09 μg/(m2·h)和9.85 μg/(m2·h)。拔節(jié)至灌漿期和成熟期,各處理土壤N2O 排放通量都處于較低水平,CK、T1 和T3 處理中土壤N2O 的平均排放通量出現(xiàn)負(fù)值,即土壤為N2O 的匯。

    表1 不同灌水量下冬小麥各生育階段土壤N2O 平均排放通量 Table 1 Soil N2O fluxes from wheat fields under different irrigation schemes at different growth stages μg/(m2·h)

    2.2 不同灌溉水平下土壤CO2 排放特征

    2.2.1 不同灌溉水平下CO2排放通量變化

    各灌水量處理土壤CO2排放通量的變化規(guī)律基本一致(圖3),在冬小麥出苗期(2018 年10 月24日左右)出現(xiàn)排放峰,其中峰值最大的為T(mén)5 處理:566.71 mg/(m2·h),最小的為T(mén)2 處理:287.91 mg/(m2·h),隨后各處理CO2排放通量下降至較低水平,越冬期(2018 年12 月11 日—2019 年3 月5 日)各處理土壤CO2排放在較低水平范圍內(nèi)波動(dòng),波動(dòng)范圍為9.35~132.32 mg/(m2·h)。4 月16 日施肥和灌水后出現(xiàn)了CO2排放峰,峰值最大的是CK,268.45 mg/(m2·h)。5 月18 日灌水后出現(xiàn)CO2排放峰,峰值最大的是T2 處理:303.29 mg/(m2·h)。

    2.2.2 不同灌溉水平下冬小麥各生育期內(nèi)CO2排放量

    冬小麥季土壤CO2累積排放量的變化如圖4所示,各處理土壤CO2累積排放量的變化規(guī)律基本一致,土壤CO2累積排放量在冬小麥全生育期內(nèi)逐漸增加。苗期各處理土壤CO2累積排放量快速增加,進(jìn)入越冬期后累積增長(zhǎng)速率逐漸減慢,返青—拔節(jié)期,土壤CO2累積排放量增長(zhǎng)速率逐漸增大,拔節(jié)—灌漿階段累積增長(zhǎng)速率達(dá)到最大,直至冬小麥成熟。冬小麥全生育期內(nèi),土壤CO2累積排放量最大的是CK,7 923.83 kg/hm2,T1—T5處理較CK分別減少了5.47%、4.75%、9.65%、7.52%和12.33%(P>0.05),差異不顯著。各處理土壤CO2平均排放通量差異不顯著(表2)。

    圖3 不同灌水量下土壤CO2 排放通量的變化 Fig.3 Temporal variation of CO2 emission flux at different irrigation schemes

    圖4 不同灌水量下小麥季土壤CO2 累積排放量 Fig.4 Cumulative CO2 emission of soil in wheat growth stage at different irrigation schemes

    表2 不同灌水量下冬小麥各生育階段土壤CO2 平均排放通量 Table 2 Soil CO2 fluxes from wheat fields under different irrigation schemes at different growth stages mg/(m2·h)

    2.3 溫度和含水率對(duì)N2O 和CO2 排放影響分析

    2.3.1 不同灌溉水平下0~5 cm 土層溫度變化

    由圖5 可知,整個(gè)冬小麥生育期內(nèi),不同處理下0~5 cm 土壤溫度具有明顯的季節(jié)變化規(guī)律。冬小麥播種后,隨著氣溫的降低,土壤溫度不斷下降,最低溫度出現(xiàn)在越冬期。進(jìn)入返青期后,隨著氣溫的升高,土壤溫度也逐漸升高。隨后0~5 cm 土壤溫度一直呈上升趨勢(shì),直至冬小麥?zhǔn)斋@,灌水會(huì)使土壤溫度出現(xiàn)短暫的下降。

    2.3.2 不同灌溉水平下土層含水率變化

    圖6 為冬小麥季不同處理0~10 cm 和10~20 cm土壤含水率的動(dòng)態(tài)變化圖。由圖6 可知,不同灌水量處理下的土壤孔隙度含水率(WFPS)的變化規(guī)律基本一致,WFPS 的變化主要受灌水和降雨的影響,在灌水和降雨后,WFPS 顯著增大,隨后逐漸降低。各處理0~10 cm 和10~20 cm 土壤WFPS 基本表現(xiàn)為CK、T1 處理最大,T4、T5 處理最小,且在灌水后這種現(xiàn)象較為明顯。

    圖5 不同處理0~5 cm 土壤溫度 Fig.5 Soil temperature of 0~5 cm under different treatments

    圖6 不同處理土層孔隙度含水率的變化 Fig.6 Soil WFPS changes under different treatments

    表3 不同處理土壤CO2 排放與土壤溫濕度的相關(guān)系數(shù) Table 3 Correlation coefficients between CO2 emission and soil temperature and humidity under different treatments

    2.3.3 土壤N2O、CO2排放通量與土壤溫度、含水率相關(guān)性分析

    考慮播種后土壤肥力對(duì)氣體排放的影響占的權(quán)重較大,取氣體排放通量穩(wěn)定后的23 組數(shù)據(jù)進(jìn)行相關(guān)性分析。表3 和表4 分別為不同灌水量處理土壤CO2排放與土壤溫濕度的相關(guān)系數(shù)表和土壤N2O 排放與土壤溫濕度的相關(guān)系數(shù)表。由表3 可知,土壤CO2排放通量與0~10 cm 土壤WFPS 存在一定正相關(guān)關(guān)系,相關(guān)系數(shù)在0.002~0.334;與10~20 cm 土壤WFPS 存在一定負(fù)相關(guān)關(guān)系, 相關(guān)系數(shù)在-0.410~-0.113;與土壤0~5 cm 地溫間存在一定正相關(guān)關(guān)系,相關(guān)系數(shù)0.247~0.657。由表4 可知,土壤N2O排放通量與0~10 cm 土壤WFPS 存在一定正相關(guān)關(guān)系,相關(guān)系數(shù)0.077~0.377;與10~20 cm 土壤WFPS 存在一定正相關(guān)關(guān)系,相關(guān)系數(shù)0.095~0.433;與土壤0~5 cm 地溫間存在一定正相關(guān)關(guān)系,相關(guān)系數(shù)0.080~0.428。

    表4 不同處理土壤N2O 排放與土壤溫濕度的相關(guān)系數(shù) Table 4 Correlation coefficients between N2O emission and soil temperature and humidity under different treatments

    2.4 不同灌溉水平對(duì)土壤N2O、CO2 排放綜合溫室效應(yīng)和冬小麥產(chǎn)量的影響

    由綜合排放溫室效應(yīng)(表6)可知,與CK 相比,T1—T5 處理N2O 和CO2的綜合排放溫室效應(yīng)分別減少了8.48%、5.70%、13.72%、13.91%和19.08%。從冬小麥產(chǎn)量情況來(lái)看(表5),T1、T3、T4 處理較CK 相比增產(chǎn)了3.4%、6.8%和3.4%,T2、T5 處理較CK 相比減產(chǎn)了3.4%和5.1%,T5 處理小麥產(chǎn)量最小(7 000 kg/hm2),隨著灌水量的增加,產(chǎn)量有增加的趨勢(shì),T3 處理下產(chǎn)量達(dá)到最大(7 875 kg/hm2),此后隨著灌水量的增加產(chǎn)量減小。綜合考慮經(jīng)濟(jì)效益和環(huán)境效益,T3 處理的灌溉用水量最佳。

    表5 冬小麥全生育期N2O 和CO2 累積排放量及產(chǎn)量 Table 5 Cumulative emission of N2O and CO2 and wheat yield during wheat entire growth period

    表6 冬小麥全生育期麥田N2O 和CO2 綜合排放溫室效應(yīng) Table 6 GWPs of N2O and CO2 emission during wheat entire growth period kg/hm2

    3 討 論

    3.1 不同灌溉水平對(duì)冬小麥土壤N2O 排放的影響

    農(nóng)田土壤N2O 排放主要來(lái)自氮素的硝化和反硝化作用,土壤水分、溫度、質(zhì)地等都會(huì)對(duì)N2O 的產(chǎn)生和擴(kuò)散產(chǎn)生影響[22],農(nóng)業(yè)生產(chǎn)中施肥和灌水是影響土壤N2O 排放的重要因素[23-24]。本試驗(yàn)中,各處理N2O 排放主要集中在施肥播種后,播種后7 d 第1 次取樣N2O 排放處于較高水平,到播種后26 d,各處理土壤N2O 排放通量均驟減到較低水平,這可能是因?yàn)椴デ肮嗟讐勊褪┘踊?,使土壤N2O 排放通量較大。隨著苗期小麥生長(zhǎng),土壤肥力減弱,N2O 排放通量便隨之迅速降低并維持在較低水平。各處理N2O 排放通量的峰值均出現(xiàn)在降雨或灌水以后,且峰值遠(yuǎn)小于播種前期排放水平,可見(jiàn)在不施肥的情況下,灌水量對(duì)土壤N2O 排放峰值的影響較施肥情況下小。4 月16 日灌水后各處理N2O 出現(xiàn)的峰值要高于5 月18 日灌水后N2O 的峰值,是因?yàn)? 月16 日灌水的同時(shí)追施尿素,而5 月18 日僅灌水沒(méi)有施肥。同一處理,不同生育階段土壤N2O 平均排放通量表現(xiàn)為苗期>越冬期,返青至拔節(jié)期>拔節(jié)至灌漿期>成熟期,苗期平均排放通量最大與初期土壤肥力充足和土壤含水率高有關(guān)。

    土壤濕度通過(guò)影響硝化及反硝化作用影響土壤N2O 的排放[25]。封克等[26]研究表明,在土壤WFPS為45%~75%時(shí),硝化和反硝化共同作用產(chǎn)生較多的N2O。本試驗(yàn)小麥全生育期內(nèi),以CK 為例,高于整個(gè)觀測(cè)期平均值97.1 μg/(m2·h)的N2O 排放通量,對(duì)應(yīng)的 0~10 cm WFPS 在 37.39%~47.0%之間,10~20cm WFPS 在37.0%~55.03%之間。而低于32 μg/(m2·h)的N2O 排放通量中,72.7%對(duì)應(yīng)的0~10 cm WFPS 在10.99%~31.37%之間,10~20 cm WFPS 在13.78%~31.65%之間,可以看出,WFPS>37%的濕度條件利于N2O 的排放。

    在旱地階段,溫度是影響N2O 排放季節(jié)變化的關(guān)鍵因子[27]。Bremner 等[28]研究發(fā)現(xiàn)N2O 排放速率的變化幾乎與表層土壤溫度同步。本試驗(yàn)相關(guān)性分析表明,各處理N2O 排放與土壤0~5 cm 地溫間存在一定正相關(guān)關(guān)系,但未達(dá)到顯著水平。同時(shí),本試驗(yàn)中監(jiān)測(cè)到N2O 排放通量出現(xiàn)負(fù)值的現(xiàn)象,關(guān)于這種現(xiàn)象出現(xiàn)的原因目前尚無(wú)定論。劉曄等[29]對(duì)北京森林生態(tài)系統(tǒng)的研究表明,當(dāng)?shù)乇頊囟雀哂?5 ℃時(shí),N2O排放通量通常為正值,當(dāng)溫度在7~15 ℃時(shí)N2O 出現(xiàn)了負(fù)排放。Barton 等[30]發(fā)現(xiàn)在含氮量較低和土壤含水率較高條件下會(huì)出現(xiàn)N2O 的吸收。對(duì)此還有待進(jìn)一步研究。

    3.2 不同灌溉水平對(duì)冬小麥土壤CO2 排放的影響

    小麥季土壤CO2排放通量具有明顯的季節(jié)變化,播種、施肥和灌水后出現(xiàn)排放高峰,這與土壤水分、土壤溫度、土壤微生物呼吸和作物根呼吸變化密切相關(guān)[31]。播種初期由于土壤肥力充足,加之旋耕使土壤疏松,各處理CO2排放均處于較高水平,隨著土壤肥力減小和越冬期溫度降低,CO2排放通量一直在較低水平范圍內(nèi)波動(dòng)(圖3)。進(jìn)入返青期后,隨著氣溫的升高,土壤溫度不斷增加,促進(jìn)了土壤微生物活性和有機(jī)質(zhì)降解等。土壤微生物呼吸逐漸增強(qiáng),作物的根系呼吸也逐漸增強(qiáng)[32],使得CO2排放通量逐漸增大并維持在較高水平。施肥和灌水后都出現(xiàn)明顯排放高峰。

    土壤水分和溫度是影響土壤呼吸排放CO2的二大主要因素[33]。土壤水分通過(guò)影響作物根系呼吸、土壤微生物呼吸等多個(gè)方面影響土壤CO2的產(chǎn)生與排放[34],其與土壤CO2排放之間的關(guān)系較為復(fù)雜,目前不同研究者對(duì)二者之間關(guān)系的研究結(jié)果有所不同。楊凡等[35]在關(guān)中地區(qū)夏玉米試驗(yàn)中發(fā)現(xiàn)土壤CO2排放通量對(duì)0~10 cm 土層WFPS 變化較10~20 cm 土層WFPS 的變化響應(yīng)敏感。楊碩歡等[32]在玉米-小麥輪作試驗(yàn)中觀測(cè)到土壤CO2排放通量與土壤水分(0~10、10~20 cm 土層)之間呈負(fù)相關(guān)關(guān)系,且相關(guān)性不顯著。土壤CO2排放峰值出現(xiàn)時(shí)間同當(dāng)次灌溉時(shí)間有滯后性,且高水分處理CO2排放通量并非全高于中、低水分處理[36],這與本試驗(yàn)結(jié)果一致。過(guò)高的土壤含水率會(huì)阻礙土壤中O2擴(kuò)散,抑制植物根系和好氧微生物的活動(dòng)[37],因此土壤CO2排放并非隨土壤濕度增大而一直增大。有研究指出,土壤水分對(duì)土壤CO2排放量的影響取決于與土壤溫度的相互協(xié)調(diào)作用[38]。本試驗(yàn)結(jié)果與之相符,從不同生育階段各處理CO2平均排放通量的方差分析來(lái)看,除越冬期外,其他階段各處理CO2排放沒(méi)有顯著性差異(表2)??梢?jiàn)在越冬期溫度較低時(shí),CO2排放處于較低水平,此時(shí)氣體排放強(qiáng)度對(duì)土壤水分的改變有明顯響應(yīng)。Mcculley 等[39]研究表明,灌溉會(huì)顯著增加土壤呼吸的溫度敏感性,但當(dāng)土壤含水率超過(guò)某個(gè)閾值,土壤呼吸的溫度敏感性反而會(huì)降低;相關(guān)分析表明,本試驗(yàn)各處理CO2排放通量與土壤WFPS 間相關(guān)關(guān)系不明顯,而與0~5 cm 土壤溫度間存在一定正相關(guān)關(guān)系,且在T5 處理下相關(guān)性極顯著(P<0.01)。

    4 結(jié) 論

    1)小麥季各處理土壤N2O 排放季節(jié)性波動(dòng)較小,變化趨勢(shì)比較一致,均由苗期的峰值下降到較低水平范圍內(nèi)波動(dòng);不同灌水水平對(duì)小麥季土壤N2O 平均排放通量有一定影響,且灌水量對(duì)土壤N2O 排放的影響在施肥之后表現(xiàn)得更加顯著。T4、T5 處理較CK 相比N2O累積排放通量顯著減少48.99%和56.11%(P<0.05)。

    2)小麥季農(nóng)田土壤CO2排放具有明顯季節(jié)變化規(guī)律,各處理均呈先減小至較低水平后波動(dòng)性升高的“U”形曲線;不同灌水水平土壤CO2平均排放通量存在差異,但均未達(dá)到顯著水平。土壤CO2累積排放量最大的是CK:7 923.83 kg/hm2,T1—T5 處理較CK分別減少了5.47%、4.75%、9.65%、7.52%和12.33%(P>0.05)。

    3)土壤N2O 排放與0~10、10~20 cm 土壤WFPS和0~5 cm 土壤溫度之間均存在一定正相關(guān)關(guān)系,土壤CO2排放與0~5 cm 土壤溫度存在正相關(guān)關(guān)系,且在T5處理下相關(guān)性極顯著(P<0.01)。

    4)與CK 相比T1—T5 處理,土壤N2O 和CO2的綜合排放溫室效應(yīng)減少了5.70%~19.08%,T3 處理小麥產(chǎn)量最大(7 875 kg/hm2),T3 處理對(duì)應(yīng)的灌溉水量能夠在穩(wěn)產(chǎn)的前提下達(dá)到節(jié)水、減排效果。

    老鸭窝网址在线观看| 男男h啪啪无遮挡| 日本av免费视频播放| 国产欧美日韩综合在线一区二区| 亚洲一级一片aⅴ在线观看| 久久ye,这里只有精品| 国产综合精华液| 性色avwww在线观看| 一边摸一边做爽爽视频免费| 波多野结衣av一区二区av| 日韩一区二区视频免费看| 欧美黄色片欧美黄色片| 国产av码专区亚洲av| 日韩人妻精品一区2区三区| 制服诱惑二区| 久久久精品94久久精品| 美女国产视频在线观看| 9色porny在线观看| 成年人午夜在线观看视频| 亚洲精品中文字幕在线视频| 天天影视国产精品| 久久久久久久久久久免费av| 午夜日本视频在线| 日韩制服骚丝袜av| 色哟哟·www| 一区二区三区四区激情视频| 侵犯人妻中文字幕一二三四区| 欧美 日韩 精品 国产| 欧美av亚洲av综合av国产av | 免费观看性生交大片5| 中文字幕精品免费在线观看视频| 成人二区视频| 久久久欧美国产精品| 欧美激情极品国产一区二区三区| 日韩av在线免费看完整版不卡| 少妇人妻 视频| 国产精品三级大全| 国语对白做爰xxxⅹ性视频网站| 亚洲欧美日韩另类电影网站| 国产日韩欧美在线精品| 精品视频人人做人人爽| 大香蕉久久网| av在线观看视频网站免费| 丰满乱子伦码专区| 老汉色∧v一级毛片| 久久av网站| 1024视频免费在线观看| 欧美激情高清一区二区三区 | 日韩 亚洲 欧美在线| 老司机影院毛片| 久久精品国产鲁丝片午夜精品| 精品第一国产精品| 两个人免费观看高清视频| 久久久久视频综合| 国产一区二区激情短视频 | 热re99久久精品国产66热6| 国产高清国产精品国产三级| 日韩一卡2卡3卡4卡2021年| 亚洲精品国产一区二区精华液| 777米奇影视久久| 成年女人毛片免费观看观看9 | 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品av麻豆狂野| 伊人亚洲综合成人网| 国产又爽黄色视频| 亚洲四区av| 永久网站在线| 性色av一级| 精品少妇一区二区三区视频日本电影 | 亚洲欧美中文字幕日韩二区| 18禁裸乳无遮挡动漫免费视频| 捣出白浆h1v1| 十八禁高潮呻吟视频| 国产激情久久老熟女| 国产精品人妻久久久影院| 久久精品人人爽人人爽视色| 国产精品久久久久久精品电影小说| 搡女人真爽免费视频火全软件| 免费黄网站久久成人精品| 一本色道久久久久久精品综合| 久久久久久人人人人人| 亚洲成国产人片在线观看| 国产国语露脸激情在线看| 极品人妻少妇av视频| 久久久久国产一级毛片高清牌| 精品久久蜜臀av无| 国产一区二区三区av在线| 狠狠婷婷综合久久久久久88av| 欧美 日韩 精品 国产| 久久午夜福利片| 亚洲国产av影院在线观看| 亚洲精品久久成人aⅴ小说| 一本大道久久a久久精品| 各种免费的搞黄视频| 91aial.com中文字幕在线观看| 久久久久久久久久人人人人人人| 永久网站在线| 久久av网站| 国产一级毛片在线| a级毛片黄视频| 91在线精品国自产拍蜜月| 国产成人精品婷婷| 一级毛片电影观看| 欧美bdsm另类| 最新的欧美精品一区二区| 99国产精品免费福利视频| 国产色婷婷99| 国产乱来视频区| 亚洲欧美日韩另类电影网站| 国产精品三级大全| xxx大片免费视频| 亚洲综合色网址| 精品午夜福利在线看| 精品卡一卡二卡四卡免费| 777久久人妻少妇嫩草av网站| 乱人伦中国视频| 成年人免费黄色播放视频| a级毛片在线看网站| 中国三级夫妇交换| 国产av国产精品国产| 天堂俺去俺来也www色官网| 国产精品久久久av美女十八| 国产熟女欧美一区二区| 国产成人91sexporn| 亚洲情色 制服丝袜| 亚洲欧美中文字幕日韩二区| 在线观看免费视频网站a站| 人体艺术视频欧美日本| 免费在线观看黄色视频的| 久久人人爽av亚洲精品天堂| 如何舔出高潮| 亚洲精品美女久久av网站| 成人毛片a级毛片在线播放| 男女边摸边吃奶| 亚洲综合色网址| 好男人视频免费观看在线| 欧美日韩亚洲高清精品| 美女高潮到喷水免费观看| 中文欧美无线码| 亚洲欧洲日产国产| 丝袜脚勾引网站| 亚洲精品一区蜜桃| 色94色欧美一区二区| 亚洲色图 男人天堂 中文字幕| 国产高清国产精品国产三级| 欧美成人午夜精品| 亚洲第一青青草原| 精品人妻偷拍中文字幕| 搡老乐熟女国产| 欧美精品亚洲一区二区| 国产精品国产av在线观看| 久久av网站| 精品人妻偷拍中文字幕| 色哟哟·www| 欧美国产精品va在线观看不卡| 久久久久久人妻| 尾随美女入室| 香蕉国产在线看| 亚洲精品国产色婷婷电影| 啦啦啦啦在线视频资源| 一个人免费看片子| 亚洲一码二码三码区别大吗| 亚洲经典国产精华液单| 国产极品粉嫩免费观看在线| 亚洲欧美中文字幕日韩二区| 久久精品国产亚洲av涩爱| 国产97色在线日韩免费| 亚洲成人一二三区av| 黑丝袜美女国产一区| √禁漫天堂资源中文www| 婷婷成人精品国产| 国产成人精品久久二区二区91 | 激情五月婷婷亚洲| 亚洲美女搞黄在线观看| 丝袜人妻中文字幕| 啦啦啦视频在线资源免费观看| 精品亚洲成国产av| 麻豆av在线久日| 啦啦啦视频在线资源免费观看| 国产精品久久久av美女十八| 精品国产乱码久久久久久小说| 亚洲av免费高清在线观看| 久久av网站| av片东京热男人的天堂| 一二三四在线观看免费中文在| 国产成人精品福利久久| 欧美精品av麻豆av| 在线观看人妻少妇| 亚洲天堂av无毛| 欧美 亚洲 国产 日韩一| 国产精品不卡视频一区二区| 爱豆传媒免费全集在线观看| 三级国产精品片| 欧美变态另类bdsm刘玥| 老司机影院毛片| 久久鲁丝午夜福利片| 天天躁狠狠躁夜夜躁狠狠躁| 一区二区三区激情视频| 啦啦啦在线免费观看视频4| 亚洲精品久久午夜乱码| 国产片内射在线| 国产黄频视频在线观看| 色94色欧美一区二区| 日韩精品免费视频一区二区三区| 日韩制服丝袜自拍偷拍| 热99久久久久精品小说推荐| 18在线观看网站| 国产综合精华液| 欧美成人精品欧美一级黄| 国产白丝娇喘喷水9色精品| 久久精品国产自在天天线| 黑人猛操日本美女一级片| 国产成人a∨麻豆精品| 老司机亚洲免费影院| 天天躁日日躁夜夜躁夜夜| 欧美日韩亚洲国产一区二区在线观看 | 男人舔女人的私密视频| 亚洲精品成人av观看孕妇| 精品人妻一区二区三区麻豆| 人体艺术视频欧美日本| 亚洲,欧美精品.| 久久精品人人爽人人爽视色| 成人手机av| 精品国产一区二区久久| 国产免费一区二区三区四区乱码| 一区福利在线观看| 性少妇av在线| 香蕉精品网在线| av.在线天堂| 日韩成人av中文字幕在线观看| 国产亚洲av片在线观看秒播厂| 青草久久国产| 亚洲精品第二区| 精品国产一区二区三区久久久樱花| 亚洲国产精品一区三区| av在线app专区| 欧美精品一区二区免费开放| 国产成人aa在线观看| 亚洲人成77777在线视频| 亚洲精品国产色婷婷电影| 久久精品国产亚洲av天美| 男的添女的下面高潮视频| 毛片一级片免费看久久久久| 捣出白浆h1v1| 黑人欧美特级aaaaaa片| 午夜免费男女啪啪视频观看| 高清av免费在线| av线在线观看网站| 一区在线观看完整版| 久久久久网色| 欧美日韩一区二区视频在线观看视频在线| 日本av手机在线免费观看| 久久国产精品男人的天堂亚洲| 国产欧美日韩综合在线一区二区| 成人二区视频| 免费人妻精品一区二区三区视频| 97人妻天天添夜夜摸| 色网站视频免费| 美女xxoo啪啪120秒动态图| 欧美在线黄色| 亚洲欧美成人综合另类久久久| 国产日韩欧美视频二区| 91午夜精品亚洲一区二区三区| 一本久久精品| 90打野战视频偷拍视频| 最近中文字幕高清免费大全6| 精品人妻在线不人妻| 99久久中文字幕三级久久日本| 亚洲三区欧美一区| 亚洲精品日韩在线中文字幕| videossex国产| 国产成人午夜福利电影在线观看| 一级片'在线观看视频| 国产乱人偷精品视频| 久久精品亚洲av国产电影网| 少妇熟女欧美另类| 男女边摸边吃奶| av电影中文网址| 欧美+日韩+精品| 免费av中文字幕在线| 欧美av亚洲av综合av国产av | 国产男女内射视频| 久久久精品区二区三区| 久久免费观看电影| 国产免费视频播放在线视频| 国产日韩欧美亚洲二区| 久久久精品国产亚洲av高清涩受| 热re99久久精品国产66热6| 精品福利永久在线观看| 欧美国产精品va在线观看不卡| 亚洲av电影在线进入| 精品一区在线观看国产| 欧美变态另类bdsm刘玥| 97在线视频观看| 久久久国产欧美日韩av| 成人国产av品久久久| 精品一品国产午夜福利视频| 色吧在线观看| 91成人精品电影| 免费播放大片免费观看视频在线观看| 中文字幕制服av| 欧美97在线视频| 欧美人与性动交α欧美精品济南到 | 日本欧美视频一区| 国产精品人妻久久久影院| 乱人伦中国视频| 中文字幕人妻熟女乱码| 国产精品国产三级国产专区5o| 纵有疾风起免费观看全集完整版| 五月天丁香电影| 美女脱内裤让男人舔精品视频| av在线播放精品| 午夜福利,免费看| 亚洲欧洲日产国产| 欧美精品高潮呻吟av久久| 精品国产一区二区久久| 国产深夜福利视频在线观看| 日本wwww免费看| 女的被弄到高潮叫床怎么办| 91国产中文字幕| 97精品久久久久久久久久精品| 免费女性裸体啪啪无遮挡网站| 伦理电影免费视频| 欧美变态另类bdsm刘玥| 伦理电影免费视频| 午夜免费男女啪啪视频观看| 亚洲欧洲日产国产| 少妇的逼水好多| 有码 亚洲区| 国产免费又黄又爽又色| 少妇的丰满在线观看| 久久精品国产a三级三级三级| 亚洲天堂av无毛| 国产精品一二三区在线看| 久久免费观看电影| 精品少妇一区二区三区视频日本电影 | 国产白丝娇喘喷水9色精品| 精品少妇久久久久久888优播| 尾随美女入室| 在线观看国产h片| 另类亚洲欧美激情| 中国国产av一级| 热99久久久久精品小说推荐| 国产精品一区二区在线不卡| 亚洲,欧美精品.| 久久久久人妻精品一区果冻| 1024香蕉在线观看| 久久人人爽av亚洲精品天堂| 人妻一区二区av| 国产日韩一区二区三区精品不卡| 国产精品二区激情视频| 看免费av毛片| 五月开心婷婷网| 国产精品女同一区二区软件| 国产精品久久久av美女十八| 久久精品国产自在天天线| av电影中文网址| 亚洲欧美色中文字幕在线| 亚洲视频免费观看视频| 亚洲精品自拍成人| 啦啦啦在线观看免费高清www| 黑人欧美特级aaaaaa片| 欧美 日韩 精品 国产| 日韩欧美一区视频在线观看| 香蕉丝袜av| 婷婷色综合大香蕉| 在线观看三级黄色| 亚洲欧美一区二区三区久久| 性色av一级| 日韩,欧美,国产一区二区三区| 国产黄频视频在线观看| 我要看黄色一级片免费的| 久久久国产精品麻豆| 成人18禁高潮啪啪吃奶动态图| 国产成人精品一,二区| 五月伊人婷婷丁香| 国产精品久久久久成人av| 考比视频在线观看| 国产成人91sexporn| 97精品久久久久久久久久精品| 久久女婷五月综合色啪小说| 少妇的丰满在线观看| 99热全是精品| 国产av精品麻豆| 2018国产大陆天天弄谢| 九色亚洲精品在线播放| 亚洲久久久国产精品| 久久久亚洲精品成人影院| 少妇的丰满在线观看| 高清欧美精品videossex| 精品亚洲成国产av| 国产有黄有色有爽视频| av片东京热男人的天堂| 一级毛片我不卡| 中文字幕人妻丝袜一区二区 | 成人黄色视频免费在线看| 不卡av一区二区三区| 热99久久久久精品小说推荐| 国产成人免费观看mmmm| 国产在视频线精品| 涩涩av久久男人的天堂| 午夜福利在线免费观看网站| 国产国语露脸激情在线看| 99久久综合免费| 亚洲美女搞黄在线观看| 在线精品无人区一区二区三| 最新的欧美精品一区二区| 国产成人91sexporn| 色视频在线一区二区三区| 精品少妇一区二区三区视频日本电影 | 黄色怎么调成土黄色| 一级爰片在线观看| 国产精品秋霞免费鲁丝片| 国产 精品1| 亚洲视频免费观看视频| 五月伊人婷婷丁香| 亚洲五月色婷婷综合| 亚洲第一区二区三区不卡| 夫妻午夜视频| 久热久热在线精品观看| 中文字幕人妻丝袜制服| 97在线视频观看| 免费黄网站久久成人精品| 国产成人精品福利久久| 日韩中字成人| 99九九在线精品视频| 新久久久久国产一级毛片| 国产成人精品久久久久久| www.熟女人妻精品国产| 女人被躁到高潮嗷嗷叫费观| 亚洲人成电影观看| 亚洲,一卡二卡三卡| 久久精品aⅴ一区二区三区四区 | 2021少妇久久久久久久久久久| 国产精品一区二区在线观看99| 看免费av毛片| 精品亚洲成国产av| 大香蕉久久成人网| 国精品久久久久久国模美| 天天影视国产精品| 日日撸夜夜添| 国产精品久久久久久精品古装| 久久久久久久久久人人人人人人| 亚洲精品久久午夜乱码| 久久久久精品人妻al黑| 一区二区三区乱码不卡18| 边亲边吃奶的免费视频| 另类精品久久| 亚洲综合色惰| 久久久精品免费免费高清| 久久久久久久久免费视频了| 免费高清在线观看日韩| 少妇被粗大的猛进出69影院| 国产高清国产精品国产三级| 亚洲中文av在线| 好男人视频免费观看在线| 大话2 男鬼变身卡| 一级黄片播放器| 午夜av观看不卡| 亚洲伊人色综图| 久久久国产一区二区| 18+在线观看网站| 啦啦啦中文免费视频观看日本| 大陆偷拍与自拍| 下体分泌物呈黄色| 免费观看a级毛片全部| 女人高潮潮喷娇喘18禁视频| 97在线视频观看| 老汉色av国产亚洲站长工具| 久久精品国产综合久久久| 国产成人精品无人区| 亚洲三区欧美一区| 精品一区二区免费观看| 两个人免费观看高清视频| 最近中文字幕高清免费大全6| 亚洲av电影在线进入| 亚洲国产色片| av免费观看日本| 在线观看人妻少妇| 久久久久人妻精品一区果冻| 性高湖久久久久久久久免费观看| 国产片特级美女逼逼视频| 国产极品天堂在线| 亚洲伊人色综图| 午夜久久久在线观看| 人人妻人人爽人人添夜夜欢视频| 亚洲欧美一区二区三区久久| 亚洲激情五月婷婷啪啪| 伊人久久大香线蕉亚洲五| 国产成人精品婷婷| 久久久久精品人妻al黑| 亚洲国产欧美网| 满18在线观看网站| 亚洲精品在线美女| 亚洲一级一片aⅴ在线观看| 国产成人午夜福利电影在线观看| 嫩草影院入口| 91午夜精品亚洲一区二区三区| 国产一级毛片在线| av又黄又爽大尺度在线免费看| 欧美激情 高清一区二区三区| 国产成人免费无遮挡视频| 亚洲av日韩在线播放| 亚洲三区欧美一区| 搡女人真爽免费视频火全软件| 欧美av亚洲av综合av国产av | 天堂8中文在线网| 久热久热在线精品观看| 久久久欧美国产精品| 亚洲国产欧美在线一区| 欧美av亚洲av综合av国产av | 国产黄频视频在线观看| 麻豆乱淫一区二区| 欧美日韩成人在线一区二区| 三级国产精品片| 91国产中文字幕| 欧美亚洲日本最大视频资源| 只有这里有精品99| 人妻系列 视频| 18在线观看网站| 另类精品久久| 亚洲天堂av无毛| 久久久久久人人人人人| 成人国产麻豆网| 伊人久久大香线蕉亚洲五| 日韩不卡一区二区三区视频在线| av线在线观看网站| 日日撸夜夜添| 免费少妇av软件| 午夜久久久在线观看| 午夜福利一区二区在线看| 国产日韩一区二区三区精品不卡| 人体艺术视频欧美日本| 韩国高清视频一区二区三区| 亚洲成国产人片在线观看| 伦精品一区二区三区| 一区二区三区乱码不卡18| 欧美日韩成人在线一区二区| 一个人免费看片子| 校园人妻丝袜中文字幕| 永久免费av网站大全| 亚洲国产成人一精品久久久| 在线观看三级黄色| 黄色配什么色好看| 水蜜桃什么品种好| 少妇人妻 视频| 国产又色又爽无遮挡免| 一本久久精品| 国产精品一国产av| 男人添女人高潮全过程视频| 国产1区2区3区精品| 精品人妻偷拍中文字幕| 久久精品国产鲁丝片午夜精品| av天堂久久9| 一级毛片黄色毛片免费观看视频| 亚洲成人一二三区av| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久国产网址| 日韩 亚洲 欧美在线| 菩萨蛮人人尽说江南好唐韦庄| av视频免费观看在线观看| av国产久精品久网站免费入址| 只有这里有精品99| 天天影视国产精品| 国产成人免费无遮挡视频| 男女边摸边吃奶| 涩涩av久久男人的天堂| 人人妻人人添人人爽欧美一区卜| 成人毛片60女人毛片免费| 中文字幕制服av| 一级毛片电影观看| 香蕉丝袜av| 国产精品久久久久久久久免| 男女边吃奶边做爰视频| 99精国产麻豆久久婷婷| 成人二区视频| 有码 亚洲区| 免费高清在线观看日韩| 久久久久久免费高清国产稀缺| 午夜激情久久久久久久| av有码第一页| 亚洲三区欧美一区| 亚洲综合色网址| 成人亚洲精品一区在线观看| 亚洲五月色婷婷综合| 日韩三级伦理在线观看| av有码第一页| 国产精品国产三级国产专区5o| 国产一区亚洲一区在线观看| 精品酒店卫生间| 麻豆乱淫一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 日韩熟女老妇一区二区性免费视频| 韩国av在线不卡| 亚洲国产欧美在线一区| videosex国产| 久久久精品区二区三区| 亚洲精品国产一区二区精华液| 久久国产精品大桥未久av| 成人二区视频| 欧美另类一区| 国产熟女欧美一区二区| 18禁动态无遮挡网站| 久久久久人妻精品一区果冻| 久久精品国产亚洲av高清一级| www.自偷自拍.com| 制服人妻中文乱码| 男女下面插进去视频免费观看| 伊人久久国产一区二区| 久久人人爽av亚洲精品天堂| 男的添女的下面高潮视频| 国产野战对白在线观看| 免费观看无遮挡的男女| 一区二区三区精品91| 成年动漫av网址|