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

    綠洲鹽化潮土有效鋅含量與鹽分離子的相關(guān)性及通徑分析①

    2017-07-21 05:05:48楊思存王成寶姜萬禮
    土壤 2017年3期
    關(guān)鍵詞:潮土負(fù)效應(yīng)缺鋅

    楊思存,霍 琳,王成寶,姜萬禮

    ?

    綠洲鹽化潮土有效鋅含量與鹽分離子的相關(guān)性及通徑分析①

    楊思存,霍 琳,王成寶,姜萬禮

    (甘肅省農(nóng)業(yè)科學(xué)院土壤肥料與節(jié)水農(nóng)業(yè)研究所,農(nóng)業(yè)部甘肅耕地保育與農(nóng)業(yè)環(huán)境科學(xué)觀測實驗站,蘭州 730070)

    以張掖市甘州區(qū)和臨澤縣綠洲鹽化潮土為研究對象,通過采集具有代表性、正常耕作種植、前茬為玉米的典型地塊0 ~ 20 cm耕層土樣分析化驗,研究了區(qū)域土壤有效鋅和鹽分離子含量特征,運(yùn)用通徑分析方法研究了它們之間的相關(guān)性及不同鹽分離子對土壤有效鋅含量的影響程度。結(jié)果表明,河西綠洲灌區(qū)鹽化潮土有效鋅含量介于0.289 ~ 0.736 mg/kg之間,平均為0.473 mg/kg,變異系數(shù)為18.68%,超過70% 的土壤低于缺鋅臨界值(0.5 mg/kg),且土壤有效鋅含量有隨含鹽量升高而降低的趨勢。土壤有效鋅含量與HCO– 3、Mg2+、Na+含量之間呈極顯著負(fù)相關(guān),與Cl–、Ca2+、K+含量之間呈顯著負(fù)相關(guān),與SO2– 4含量之間的相關(guān)性不顯著。Mg2+對土壤有效鋅含量表現(xiàn)出強(qiáng)烈的直接和間接負(fù)效應(yīng),HCO– 3表現(xiàn)出很大的直接負(fù)效應(yīng)和強(qiáng)烈的間接正效應(yīng),Cl–、K+和Na+均表現(xiàn)出較強(qiáng)的直接和間接負(fù)效應(yīng),SO2– 4表現(xiàn)出較弱的直接和間接正效應(yīng),Ca2+表現(xiàn)出較弱的直接和間接負(fù)效應(yīng)。Mg2+、HCO– 3和Na+是影響土壤有效鋅含量變化的3個主導(dǎo)因素,對決定系數(shù)2的貢獻(xiàn)超過了68%。

    綠洲鹽化潮土;有效鋅含量;鹽分離子;相關(guān)分析;通徑分析

    鋅是植物的必需營養(yǎng)元素,植物所需的鋅主要來自土壤[1]。全世界有大范圍的缺鋅土壤,亞洲有200萬hm2缺鋅稻田[2-5],由此所引發(fā)的作物籽粒和人體缺鋅問題日益引起全球的廣泛關(guān)注[6-9]。我國有40% 的缺鋅土壤,主要分布在北方石灰性土壤中,河西走廊的綠洲鹽化潮土也是典型缺鋅區(qū)域之一[10-11]。蘆滿濟(jì)等[12-13]對該地區(qū)非鹽化土壤的調(diào)查結(jié)果表明,耕層土壤有效鋅含量為0.74 mg/kg,有90.9% 的樣本高于土壤缺鋅臨界值(0.5 mg/kg),在田間也很少發(fā)現(xiàn)有玉米缺鋅植株。但在該地區(qū)鹽化土壤上的研究結(jié)果表明,土壤有效鋅平均含量只有0.34 mg/kg,玉米缺鋅癥狀非常普遍,而且不同類型、不同程度鹽漬化土壤中有效鋅含量和植株缺鋅率的差異也較大。在中度氯化物-硫酸鹽鹽漬土上,土壤有效鋅含量達(dá)到0.57 mg/kg,高于土壤缺鋅臨界值,但仍然有18.3% 的玉米植株表現(xiàn)出缺鋅癥狀;而在重度鎂質(zhì)鹽漬土上,土壤有效鋅含量只有0.22 mg/kg,玉米植株缺鋅率達(dá)到了37.4%,即便是施入150 kg/hm2硫酸鋅,仍不能完全解決苗期生長緩慢、全生育期花葉病普遍發(fā)生、后期果穗發(fā)育不良和產(chǎn)量低的問題[14]。同時還發(fā)現(xiàn),在施用了鈣鎂磷肥、磷酸鎂銨、硫酸鎂等含鎂肥料的土壤上,玉米的“花白苗”現(xiàn)象越嚴(yán)重,干枯死亡率在10% ~ 20%,表現(xiàn)出典型缺鋅癥狀。由此我們推測,河西綠洲灌區(qū)鹽漬化環(huán)境(鹽分濃度和鹽分離子類型)及不合理的施肥措施可能是造成土壤有效鋅含量降低和作物缺鋅的重要因素。

    目前關(guān)于鹽分離子對土壤-作物系統(tǒng)鋅營養(yǎng)影響的研究報道,主要考慮的是單鹽或單一鹽分離子的影響[15-20],在研究方法上也大多局限于簡單相關(guān)、多元回歸分析等。然而,簡單相關(guān)并不能全面考察變量間的相互關(guān)系,其結(jié)果往往帶有一定的片面性;多元回歸分析雖然在一定程度上能夠消除變量之間的多重共線性,能夠真實地表現(xiàn)出各個自變量和因變量的關(guān)系,但由于偏回歸系數(shù)帶有單位,使各自變量對因變量的效應(yīng)不能直接進(jìn)行比較,限制了人們對鋅鹽關(guān)系更深入和準(zhǔn)確的認(rèn)識。而通徑分析是對多元回歸分析的拓展,由美國數(shù)量遺傳學(xué)家Sewall Wright 于1921年提出,它在多元回歸的基礎(chǔ)上將自變量與因變量的單相關(guān)系數(shù)加以分解,分解出自變量對因變量的直接影響力和間接影響力。因而,通徑系數(shù)是介于相關(guān)系數(shù)與回歸系數(shù)之間的統(tǒng)計量,其經(jīng)過標(biāo)準(zhǔn)化后去掉了單位,彼此間可以相互比較,從而可以反映各自變量對因變量的影響程度和相對重要性[21-24]。因此,將相關(guān)分析、回歸分析與通徑分析相結(jié)合,可以提供更多、更精確的因果信息,對多變量資料的統(tǒng)計分析更為合理[25]。本文以河西綠洲灌區(qū)鹽化潮土為研究對象,在典型地區(qū)采樣,在對土壤有效鋅含量與鹽分離子含量相關(guān)分析的基礎(chǔ)上,進(jìn)一步運(yùn)用通徑分析方法,探討了各鹽分離子對土壤有效鋅含量的影響程度和作用機(jī)理,旨在找出影響土壤鋅有效性含量的主要因素,為綠洲鹽化潮土的改良利用提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    研究區(qū)位于甘肅省張掖市甘州區(qū)和臨澤縣(38°54′ ~ 39°18′N,99°57′ ~ 100°42′E),河西走廊的中部,黑河流經(jīng)其中,是典型的綠洲農(nóng)業(yè)和大型灌溉農(nóng)業(yè)區(qū)。該區(qū)處在溫帶氣候區(qū)和暖溫帶偏干旱荒漠氣候交匯地帶,甘州區(qū)海拔1 370 ~ 2 200 m,年均降水量129.0 mm、蒸發(fā)量2047.9 mm,年均氣溫7.3℃,日照時數(shù)3 085.1 h,無霜期153 d;臨澤縣海拔1 380 ~ 2 278 m,年均降水量113.4 mm,蒸發(fā)量2 341.0 mm,年均氣溫7.6℃,日照時數(shù)2 965 h,無霜期178 d。甘州區(qū)總灌溉面積約5.09×104hm2,其中耕地鹽漬化面積1.41×104hm2,占27.74%;臨澤縣總灌溉面積約1.97×104hm2,其中耕地鹽漬化面積1.49×104hm2,占75.47%。研究區(qū)大部分耕地土壤是由鹽化草甸土經(jīng)耕作熟化演變而來,在土壤潮化、耕作熟化及土壤鹽漬化3種過程共同作用下形成的鹽化潮土[26]。研究區(qū)的作物類型以小麥、玉米、葵花、馬鈴薯等為主,在河西綠洲灌區(qū)有典型代表性。

    1.2 土樣采集

    土壤樣品采集于2012年秋季作物收獲后,冬灌開始前進(jìn)行(10月28日—11月1日),此時的土壤鹽分變化受灌水、凍融影響較小,相對比較穩(wěn)定,基本上能反映河西綠洲灌區(qū)耕地土壤鹽漬化狀況。樣點(diǎn)的布局涵蓋了鹽化潮土的主要類型(黑潮土、灰潮土、青潮土)和程度(輕度、中度、重度)的差異。在各采樣區(qū)選擇具有代表性、無污染、正常耕作種植、前茬為玉米的典型地塊作為采樣點(diǎn),每個樣點(diǎn)采用S形采樣法采集5個耕層(0 ~ 20 cm)樣品,以充分混合后的土壤代表該樣點(diǎn)土樣。共采集的樣點(diǎn)有41個,帶回實驗室后自然風(fēng)干、磨碎,過2 mm 篩后備用。

    1.3 測定項目與方法

    土壤鹽分離子包含Ca2+、Mg2+、K+、Na+、SO2– 4、Cl–、CO32–和HCO– 3,Ca2+、Mg2+用EDTA滴定法測定;K+、Na+用火焰光度法測定;Cl–用AgNO3滴定法測定;SO2– 4用EDTA間接滴定法測定;CO2– 3和HCO– 3用雙指示劑鹽酸滴定法測定;DTPA-Zn以水土比為2∶1浸提(pH 7.3),用原子吸收分光光度計(AA320 CRT)測定[27]。

    1.4 數(shù)據(jù)處理及分析

    數(shù)據(jù)處理和統(tǒng)計分析在 Microsoft Office Excel 2003 和 SPSS 18.0 軟件中完成。先分別計算出土壤鹽分離子和有效鋅含量的平均數(shù)()、標(biāo)準(zhǔn)差()、變異系數(shù)()及相關(guān)系數(shù),再采用逐步回歸分析方法建立土壤鹽分離子對有效鋅含量的多元回歸方程,對顯著性影響因子進(jìn)行識別,然后根據(jù)回歸分析結(jié)果對識別出的顯著性影響因子進(jìn)行通徑分析,計算變量標(biāo)準(zhǔn)化后的偏回歸系數(shù)(即通徑系數(shù)),再將通徑系數(shù)分解為直接作用和間接作用,比較各鹽分離子對土壤有效鋅含量的影響程度和相對重要性,最后確定影響土壤有效鋅含量的主要因素。

    2 結(jié)果與分析

    2.1 土壤有效鋅含量特征

    西綠洲灌區(qū)鹽化潮土的有效鋅含量總體比較低,介于0.289 ~ 0.736 mg/kg,平均為0.473 mg/kg,大部分樣點(diǎn)在土壤缺鋅臨界值附近(0.5 mg/kg),占到了樣本總量的83%(圖1)。按照王遵親等[28]的分級標(biāo)準(zhǔn)對41個土樣的鹽漬化程度進(jìn)行分類(表1),可以看出土壤含鹽量及其離子組成對土壤有效鋅含量的影響較大,輕度鹽漬化土壤介于0.477 ~ 0.736 mg/kg,平均為0.546 mg/kg;中度鹽漬化土壤介于0.364 ~ 0.553 mg/kg,平均為0.463 mg/kg;重度鹽漬化土壤介于0.289 ~ 0.454 mg/kg,平均為0.379 mg/kg,說明隨著土壤鹽漬化程度的逐步加重,土壤有效鋅含量在逐漸降低。

    表1 耕層土壤鹽分離子及有效鋅含量特征

    2.2 土壤鹽分離子分布特征

    2.2.1 陽離子分布特征 從耕層土壤陽離子含量和組成比例來看(圖2A),以Ca2+含量最高,在0.152 ~ 4.588 g/kg之間,平均為2.007 g/kg,占陽離子總量的64.8%;其次是Na+,含量在0.128 ~ 2.112 g/kg之間,平均為0.489 g/kg,占15.8%;再次是Mg2+,含量在0.122 ~ 0.955 g/kg之間,平均為0.376 g/kg,占12.2%;K+含量最低,在0.122 ~ 0.375 g/kg之間,平均只有0.223 g/kg,只占7.2%。

    變異系數(shù)是反映變量離散程度的重要指標(biāo),在一定程度上揭示了變量的空間分布特性。從耕層土壤陽離子含量的變異系數(shù)來看(圖2B),Na+的變異系數(shù)最高,達(dá)到了108.6%,表明土壤中Na+的變異性很強(qiáng),說明河西綠洲灌區(qū)鹽化潮土的Na+含量分布不均勻,空間異質(zhì)性較強(qiáng)。其次是Ca2+和Mg2+,變異系數(shù)分別為76.8% 和55.0%。K+的變異系數(shù)最小,只有30.5%,說明K+受環(huán)境因素的影響較小,在土壤中的分布比較均勻。

    2.2.2 陰離子分布特征 從耕層土壤陰離子含量和組成比例來看(圖3A),河西綠洲灌區(qū)鹽化潮土中沒有CO2– 3,HCO– 3含量最低,在0.195 ~ 0.711 g/kg之間,平均為0.413 g/kg,約占陰離子總量的10% 左右;Cl–含量較高,在0.884 ~ 4.320 g/kg之間,平均為0.820 g/kg,約占陰離子總量的20%;SO2– 4含量最高,在0.263 ~ 9.696 g/kg之間,平均為2.894 g/kg,占到了陰離子總量的70%。

    從圖3B可以看出,Cl–的變異系數(shù)最高,達(dá)到了137.4%,表明土壤中Cl–的變異性很強(qiáng),這主要是因為Cl–容易隨水移動,受土體構(gòu)型和灌溉水平影響較大。其次是SO2– 4,變異系數(shù)為109.2%,變異程度也較大。HCO– 3的變異系數(shù)最小,只有32.6%,說明HCO– 3在土壤中的分布比較均勻。

    2.3 土壤有效鋅含量與鹽分離子的通徑分析

    2.3.1 土壤有效鋅含量與鹽分離子之間的相關(guān)性 相關(guān)分析結(jié)果表明(表2),土壤有效鋅含量與HCO– 3、Mg2+、Na+含量呈極顯著負(fù)相關(guān),與Cl–、Ca2+、K+含量呈顯著負(fù)相關(guān),與SO2– 4含量呈正相關(guān),但達(dá)不到顯著水平,相關(guān)系數(shù)大小順序為:Mg2+>HCO– 3>Na+>K+>Cl–>Ca2+>SO2– 4。從鹽分離子之間的相關(guān)性來看,HCO– 3與SO2– 4、Ca2+、Mg2+、K+之間存在著極顯著的負(fù)相關(guān),與Cl–之間存在著顯著的負(fù)相關(guān);Cl–與Mg2+、K+、Na+之間,SO2– 4與Ca2+、Mg2+、Na+之間,Ca2+與Mg2+之間,Mg2+與K+、Na+之間,K+與Na+之間都存在著極顯著的正相關(guān)。這些相關(guān)性的信息與人們的普遍認(rèn)識基本一致,體現(xiàn)了鹽漬化環(huán)境中鹽分離子的分化關(guān)系,以及不同鹽分離子對土壤有效鋅含量的直接和間接影響。

    2.3.2 土壤有效鋅含量與鹽分離子之間的逐步多元回歸分析 以鹽分離子HCO– 3(1),Cl–(2),SO2– 4 (3),Ca2+(4),Mg2+(5),K+(6),Na+(7)為自變量,以土壤有效鋅含量(A-Zn,)為因變量,進(jìn)行多元逐步回歸分析。偏回歸系數(shù)顯著性檢驗結(jié)果:= 0.922 3,= 47.389,<0.01,表明土壤有效鋅與鹽分離子間的復(fù)相關(guān)關(guān)系和多元回歸關(guān)系真實存在,可以進(jìn)行通徑分析,進(jìn)而得到如下標(biāo)準(zhǔn)多元回歸方程:

    = 0.720 6 – 0.309 41– 0.278 52+ 0.086 63– 0.117 54– 0.432 75– 0.234 96– 0.215 67

    由方程可以看出,所選的7種鹽分離子共同解釋了綠洲鹽化潮土有效鋅含量變異的92.23%,其中HCO– 3、Mg2+、Na+含量的影響極顯著(<0.01),Cl–、Ca2+、K+含量的影響顯著(<0.05),這與簡單相關(guān)分析的結(jié)果是一致的。

    表2 土壤有效鋅含量與各鹽分離子間的相關(guān)系數(shù)

    注:* 表示相關(guān)性達(dá)到<0.05顯著水平,** 表示相關(guān)性達(dá)到<0.01顯著水平。

    2.3.3 土壤有效鋅含量與鹽分離子之間的通徑分析 由于鹽分離子之間存在較強(qiáng)的相關(guān)性(表3),在進(jìn)行回歸分析時可能出現(xiàn)共線性現(xiàn)象,并且由于各鹽分離子含量的變動范圍不同,因此采用逐步回歸分析還不能直觀地體現(xiàn)各鹽分離子對有效鋅含量的貢獻(xiàn)大小,而通過標(biāo)準(zhǔn)化回歸系數(shù)的方法計算通徑系數(shù),并將相關(guān)系數(shù)分解為直接通徑系數(shù)與間接通徑系數(shù)的代數(shù)和,就能較直觀地反映各鹽分離子對土壤有效鋅含量的影響作用,其結(jié)果見表3。由表3可以看出,Mg2+含量對土壤有效鋅含量的直接效應(yīng)最大,表現(xiàn)為負(fù)效應(yīng),直接通徑系數(shù)為–0.432 7;Cl–、SO2– 4、Ca2+、K+、Na+等其他鹽分離子對有效鋅含量也產(chǎn)生了強(qiáng)烈的負(fù)效應(yīng)的間接作用(–1.074 2),使得Mg2+含量與有效鋅含量表觀上顯示出一種顯著的負(fù)相關(guān)關(guān)系;但通過HCO– 3對有效鋅含量產(chǎn)生的是正效應(yīng)(0.340 1),這種正效應(yīng)部分抵消了Mg2+對土壤有效鋅的強(qiáng)烈負(fù)效應(yīng)。HCO– 3表現(xiàn)出很大的直接負(fù)效應(yīng)(–0.309 4),但其通過Cl–、SO2– 4等鹽分離子所表現(xiàn)出的是強(qiáng)烈的間接正效應(yīng)(1.053 4),而且這種正效應(yīng)遠(yuǎn)遠(yuǎn)超出了直接負(fù)效應(yīng)。因此,若由相關(guān)系數(shù)–0.747 4、0.701 6(表2)就簡單地認(rèn)為Mg2+、HCO– 3含量與土壤有效鋅含量僅具有顯著的直接負(fù)影響是不恰當(dāng)?shù)?,必須綜合考慮其直接影響和間接影響。此外,Cl–、K+和Na+含量對土壤有效鋅含量均表現(xiàn)出較強(qiáng)的直接和間接負(fù)效應(yīng),SO2– 4表現(xiàn)出較弱的直接和間接正效應(yīng),Ca2+表現(xiàn)出較弱的直接和間接負(fù)效應(yīng)。

    表3 土壤鹽分離子對有效鋅含量影響的通徑系數(shù)

    注:間接作用表示鹽分離子兩兩之間對土壤有效鋅含量的作用,最左邊一列的鹽分離子在間接作用中起主要作用。

    2.3.4 鹽分離子對土壤有效鋅含量的決定程度分析 從決定系數(shù)計算結(jié)果可以看出(表4),按絕對值大小,在鹽化潮土上對土壤有效鋅含量的影響順序為:D5(0.187 2)>D5x6(0.140 9)>D5x7(0.118 6)>D1x5(0.118 3)>D2x7(0.115 3)>D1x4(0.104 7)>D1(0.095 7)>D1x3(0.093 6)>D2x6(0.087 6)>其他。表明Mg2+、HCO– 3和Cl–是影響鹽化潮土中有效鋅含量變化的主導(dǎo)因素,由它們間的相互作用共同控制著有效鋅轉(zhuǎn)化的方向與程度。其他因子的決定程度相對較小,不再作進(jìn)一步分析和討論。

    2.3.5 鹽分離子對回歸方程估測可靠程度的貢獻(xiàn)分析 從鹽分離子對有效鋅含量的決定系數(shù)R的貢獻(xiàn)來看(表5),Mg2+、HCO– 3和Na+通過直接作用和間接作用成為影響土壤有效鋅含量變化的3個主導(dǎo)因素,對決定系數(shù)R的貢獻(xiàn)超過了68%;其次是Cl–和K+,Ca2+和SO2– 4對土壤有效鋅含量的影響相對較弱。

    表4 各鹽分離子通徑分析的決定系數(shù)

    注:表中對角線數(shù)值為單個鹽分離子對有效鋅含量的決定系數(shù),對角線以上數(shù)值為兩種鹽分離子共同對有效鋅含量的決定系數(shù)。

    表5 土壤鹽分離子對有效鋅含量的決定系數(shù)R2的貢獻(xiàn)

    3 討論

    關(guān)于石灰性土壤鋅的有效性問題,國內(nèi)外學(xué)者開展了大量研究,普遍認(rèn)為土壤pH、有機(jī)質(zhì)、黏粒含量、碳酸鈣、養(yǎng)分元素、共存的金屬離子等都會對土壤-作物系統(tǒng)鋅營養(yǎng)產(chǎn)生影響[29–34],但在鹽漬化環(huán)境對土壤有效鋅含量的影響方面,開展的相關(guān)研究并不多。Hajiboland等[6]認(rèn)為較高的pH和大量CaCO3對鋅的吸附固定是造成鹽漬化土壤鋅活性較低的主要原因;LU等[11]卻認(rèn)為土壤中添加大量CaCO3并不一定會產(chǎn)生大量的HCO– 3,進(jìn)而未能對土壤中的有效鋅含量產(chǎn)生明顯影響;孫桂芳等[35]更是認(rèn)為,在含鈉較高的堿性土壤上,隨著pH的升高,更容易形成溶解性較好的鋅酸鈉,從而增加土壤鋅的有效性。本研究表明,河西綠洲灌區(qū)鹽化潮土的有效鋅含量介于0.289 ~ 0.736 mg/kg之間,平均為0.473 mg/kg,低于土壤缺鋅臨界值(0.5 mg/kg),隨著含鹽量的升高,土壤有效鋅有降低的趨勢,這是因為河西綠洲灌區(qū)的鹽化潮土是典型的石灰性土壤,土壤pH大多在8.5以上,有時甚至高達(dá)10.0,土壤有效鋅含量較低也是必然。土壤有效鋅含量與鹽分離子間的相關(guān)性分析表明,與HCO– 3、Mg2+、Na+含量呈極顯著負(fù)相關(guān),與Cl–、Ca2+、K+含量呈顯著負(fù)相關(guān),這與蘆滿濟(jì)等[14]、緱倩倩等[36]的研究結(jié)果也是一致的,但土壤有效鋅含量與SO2– 4含量之間呈不顯著的正相關(guān),這一研究結(jié)果與人們對鹽堿地的普遍認(rèn)識不同,還有待更進(jìn)一步深入研究。

    本論文通徑分析結(jié)果表明,鹽分離子通過直接和間接作用共同影響著土壤有效鋅含量,但不同鹽分離子的作用機(jī)理和影響強(qiáng)度不同。Mg2+在土壤中的含量雖然不高,但在鹽化潮土中所起的作用卻不小,不僅直接影響著土壤有效鋅含量,而且通過Cl–、SO2– 4、Ca2+、K+、Na+等其他鹽分離子對土壤有效鋅產(chǎn)生間接的負(fù)效應(yīng),這也進(jìn)一步印證了Mg2+和Zn2+之間拮抗作用的存在[37-39]。HCO– 3在土壤中的含量也不高,對土壤有效鋅含量產(chǎn)生了較大的直接負(fù)效應(yīng),這可以理解為一部分Zn2+通過化學(xué)沉淀或吸附的方式降低了土壤鋅的有效性[6, 11, 18, 40-41],但其通過Cl–、SO2– 4等鹽分離子所表現(xiàn)出的是強(qiáng)烈的間接正效應(yīng),而且這種正效應(yīng)遠(yuǎn)遠(yuǎn)超出了直接負(fù)效應(yīng),說明HCO– 3在降低土壤有效鋅的同時,Cl–、SO2– 4等鹽分離子卻對HCO– 3產(chǎn)生了拮抗作用。SO2– 4和Ca2+是土壤中含量最多的鹽分離子,但在對土壤有效鋅含量的影響方面卻不是最主要的,SO2– 4表現(xiàn)出較弱的直接和間接正效應(yīng),這有可能是ZnSO4是土壤鋅存在的較穩(wěn)定的有效狀態(tài),但這種狀態(tài)的鋅含量太少;Ca2+表現(xiàn)出較弱的直接和間接負(fù)效應(yīng),這有可能是因為Zn2+與Ca2+的離子半徑比較接近,可以與土壤晶格中的Ca2+發(fā)生同晶置換作用,從而降低土壤有效Zn含量,但相比較而言,土壤晶格中Mg2+和OH–間的結(jié)合力更弱,因此Mg2+對有效鋅含量的影響比Ca2+更顯著[42]。雖然目前關(guān)于鹽分離子影響有效鋅含量的研究并不多,其影響機(jī)理還需要今后進(jìn)一步深入研究,但總體而言,通徑分析比簡單相關(guān)和多元回歸分析更能客觀地反映各鹽分離子與土壤有效鋅含量之間的關(guān)系,能夠為綠洲鹽化潮土缺鋅問題的解決提供決策依據(jù)。在改良利用措施方面,一方面要通過灌溉洗鹽等措施來減少土壤鹽分離子含量,特別是HCO– 3的含量;另一方面,要盡量避免施用鈣鎂磷肥、磷酸鎂銨、硫酸鎂等含鎂的肥料,以減少土壤中鋅的吸附固定。

    4 結(jié)論

    1) 河西綠洲灌區(qū)鹽化潮土的有效鋅含量介于0.289 ~ 0.736 mg/kg之間,平均為0.473 mg/kg,低于土壤缺鋅臨界值(0.5 mg/kg),隨著含鹽量的升高,土壤有效鋅含量有降低的趨勢。

    2) 土壤有效鋅含量與鹽分離子間的簡單相關(guān)性分析表明,與HCO– 3、Mg2+、Na+含量呈極顯著負(fù)相關(guān),與Cl–、Ca2+、K+含量呈顯著負(fù)相關(guān),與SO2– 4含量之間的相關(guān)性不顯著。

    3) 通徑分析結(jié)果表明,Mg2+對土壤有效鋅含量表現(xiàn)出強(qiáng)烈的直接和間接負(fù)效應(yīng),HCO– 3表現(xiàn)出很大的直接負(fù)效應(yīng)和強(qiáng)烈的間接正效應(yīng),Cl–、K+和Na+對土壤有效鋅含量均表現(xiàn)出較強(qiáng)的直接和間接負(fù)效應(yīng),SO2– 4表現(xiàn)出較弱的直接和間接正效應(yīng),Ca2+表現(xiàn)出較弱的直接和間接負(fù)效應(yīng)。

    4) Mg2+、HCO– 3和Na+通過直接作用和間接作用成為影響土壤有效鋅含量變化的3個主導(dǎo)因素,對決定系數(shù)R的貢獻(xiàn)超過了68%。

    [1] Sommer A L, Lipman C B. Evidence on the indispensible nature of zinc and boron for higher green plant[J]. Plant Physiology, 1926, 1: 231–249

    [2] Takkar P N, Walker C D. The Distribution and Correction of Zinc Deficiency[M]. Netherlands: Springer, 1993: 151–165

    [3] Alloway B J. Zinc in Soils and Crop Nutrition[M]. Paris: International Fertilizer Industry Association, 2008

    [4] Shivay Y S, Kumar D, Prasad R, et al. Relative yield and zinc uptake by rice from zinc sulphate and zinc oxide coatings onto urea[J]. Nutrient Cycling in Agroecosystems, 2008, 80(2) : 181–188

    [5] Muhanned A, Farah A N, Tayel E H, et al. Zinc adsor-ption–desorption isotherms: Possible effects on the calcare-ous vertisol soils from Jordan[J]. Environmental Earth Sciences, 2012, 65(7): 2079–2085

    [6] Hajiboland R, Yang X E, Romheld V. Effects of bicarbonate and high pH on growth of Zn-efficient and Zn-inefficient genotypes of rice, wheat and rye[J]. Plant and Soil, 2003, 250: 349–357

    [7] Zou C Q, Zhang Y Q, Rashid A,et al. Biofortification of wheat with zinc through zinc fertilization in seven countries[J]. Plant Soil, 2012, 361: 119–130

    [8] Ma G S, Jin Y, Li Y P, et al. Iron and zinc deficiencies in China: What is a feasible and cost-effective strategy?[J]. Public Health Nutrition, 2008, 11(6): 632–638

    [9] Cakmak I. Enrichment of cereal grains with zinc: Agronomic or genetic biofortification?[J]. Plant and Soil, 2008, 302(1): 1–17

    [10] 劉錚. 我國土壤中鋅含量的分布規(guī)律[J]. 中國農(nóng)業(yè)科學(xué), 1994, 27(1): 30–37

    [11] Lu X C, Cui J, Tian X H, et al. Effects of zinc fertilization on zinc dynamics in potentially zinc-deficient calcareous soil[J]. Agronomy Journal, 2012, 104: 963–969

    [12] 蘆滿濟(jì), 祁國元, 楊思存, 等. 綠洲鹽化潮土的鋅肥效應(yīng)及適宜用量[J]. 西北農(nóng)業(yè)學(xué)報, 1999, 8(1): 69–73

    [13] 蘆滿濟(jì), 吳惠蘭, 胡梅. 河西綠洲潮土有效鋅含量特征及施鋅條件[J]. 甘肅農(nóng)業(yè)科技, 1999(2): 39–40

    [14] 蘆滿濟(jì), 祁國元, 楊思存. 土壤鹽分對土壤植物鋅營養(yǎng)和施鋅效果的影響// 李生秀. 土壤—植物營養(yǎng)研究文集[C]. 西安: 陜西科學(xué)技術(shù)出版社, 1999: 746–752

    [15] Yang X E, Romheld V, Marschner H. Effect of bicarbonate on root growth and accumulation of organic acids in Zn-inefficient and Zn-efficient rice cultivars(L.)[J]. Plant and Soil, 1994, 164: 1–7

    [16] Kipcak A S, Senberber F T, Derun E M, et al. Characterization and thermal dehydration kinetics of zinc borates synthesized from zinc sulfate and zinc chloride[J]. Research on Chemical Intermediates, 2015, 41: 9129–9143

    [17] Imran M, Rehim A, Sarwar N, et al. Zinc bioavailability in maize grains in response of phosphorous–zinc interaction[J]. Journal of Plant Nutrition and Soil Science, 2016, 179(1): 60–66

    [18] 錢金紅, 謝振翅. 碳酸鹽對土壤鋅解吸影響的研究[J]. 土壤學(xué)報, 1994, 31(1): 105–108

    [19] 田霄鴻, 買文選, 陸欣春, 等. 重碳酸根對不同小麥基因型生長及鋅營養(yǎng)的影響[J]. 植物營養(yǎng)與肥料學(xué)報, 2008, 14(1): 9–16

    [20] Shahab Q, Afzal M, Sarfaraz Q. Response of maize to different rates and methods of zinc application[J]. Environment and Plant Systems, 2015, 1: 43–47

    [21] Sarawgi A K, Rastogi N K, Soni D K. Correlation and path analysis in rice accessions from Madhya Pradesh[J]. Field Crops Research, 1997, 52: 161–167

    [22] Zhao H, Zhang Z B, Shao H B, et al. Genetic correlation and path analysis of transpiration efficiency for wheat flag leaves[J]. Environmental and Experimental Botany, 2008, 64: 128–134

    [23] Lenzen M. Structural path analysis of ecosystem net- works[J]. Ecological Modelling, 2006, 200: 334–342

    [24] 徐強(qiáng), 遲鳳琴, 匡恩俊, 等. 基于通徑分析的土壤性質(zhì)與硒形態(tài)的關(guān)系—— 以黑龍江省主要類型土壤為例[J]. 土壤, 2016, 48(5): 992- 999

    [25] 蓋鈞鎰. 試驗統(tǒng)計方法 [M]. 4版. 北京: 中國農(nóng)業(yè)出版社, 2013

    [26] 甘肅省土壤普查辦公室編著. 甘肅土壤[M]. 北京: 中國農(nóng)業(yè)出版社, 1991

    [27] 魯如坤. 土壤農(nóng)業(yè)化學(xué)分析方法[M]. 北京: 中國農(nóng)業(yè)科技出版社, 1999

    [28] 王遵親, 祝壽泉, 俞仁培, 等. 中國鹽漬土[M]. 北京: 科學(xué)出版社, 1993

    [29] Gupta R K, Elshout S V D, Abbol I P. Effect of pH on zinc adsorption-precipitation reactions in an alkali soil[J]. Soil Science, 1987, 143(3): 198–204

    [30] Obrador A, Alvarez J M, Lopez V L M, et al. Relationships of soil properties with Mn and Zn distribution in acidic soils and their uptake by a barley crop[J]. Geoderma, 2007, 137(3): 432–443

    [31] Zahedifar M, Karimian N, Yasrebi J. Influence of applied zinc and organic matter on zinc desorption kinetics in calcareous soils[J]. Archives of Agronomy and Soil Science, 2012, 58(2): 169–178

    [32] Kamali S, Ronaghi A, Karimian N. Soil zinc transfor-mations as affected by applied zinc and organic mate-rials[J]. Communications in Soil Science and Plant Anal-ysis, 2011, 42(9): 1038–1049

    [33] Sun Q, Chen J, Ding S M, et al. Comparison of diffusive gradients in thin film technique with traditional methods for evaluation of zinc bioavailability in soils[J]. Environmental Monitoring and Assessment, 2014, 186(10): 6553–6564

    [34] Hussain S, Maqsood M A, Aziz T, et al. Zinc bioa-vailability response curvature in wheat grains under incre-mental zinc applications[J]. Archives of Agronomy and Soil Science, 2012, 59(7):1–16

    [35] 孫桂芳, 楊光穗. 土壤-植物系統(tǒng)中鋅的研究進(jìn)展[J]. 華南熱帶農(nóng)業(yè)大學(xué)學(xué)報, 2002, 8(2): 22–30

    [36] 緱倩倩, 韓致文, 屈建軍, 等. 秦王川灌區(qū)農(nóng)田土壤鹽分離子特征分析[J]. 土壤, 2014, 46(1): 100–106

    [37] 楊思存, 霍琳, 王成寶, 等. 綠洲鹽化潮土施鎂對玉米幼苗生長、活性氧自由基代謝和鋅營養(yǎng)的影響[J]. 植物營養(yǎng)與肥料學(xué)報, 2015, 21(2): 354–361

    [38] Shula L M, Rathan P K. Zinc magnesium interactions in soil and plants as measured by crop response[J]. Journal of Tropical Agriculture, 1983(1): 309–316

    [39] Jansen S, Blust R, van Leeuwen H P. Metal speciation dynamics and bioavailability: Zn (II) and Cd (Ⅱ) uptake by mussel () and carp () [J]. Environmental Science & Technology, 2002, 36: 2164– 2170

    [40] Behera S K, Singh M V, Singh K N, et al. Distribution variability of total and extractable zinc in cultivated acid soils of India and their relationship with some selected soil properties[J]. Geoderma, 2011, 162(3-4): 242–250

    [41] Marschner H. Zinc uptake from soil // Robson, A D. Zinc in Soils and Plants[C]. Dordrecht: Kluwer Academic Pub-l-is-hers, 1993, 59–77

    [42] 劉錚, 朱其清. 微量元素的農(nóng)業(yè)化學(xué)[M]. 北京: 農(nóng)業(yè)出版社, 1991

    Correlation and Path Analyses of Available Zinc Contents and Salt Ions in Saline Fluvo-aquic Soil of Hexi Oasis Area

    YANG Sicun, HUO Lin, WANG Chengbao, JIANG Wanli

    (Institute of Soil, Fertilizer and Water-saving Agriculture, Gansu Academy of Agricultural Sciences, Lanzhou, Gansu Scientific Observing and Experiment Station of Agro-Environment and Arable Land Conservation, Ministry of Agriculture, Lanzhou 730070, China)

    Typical saline fluvo-aquic topsoil samples (0-20 cm) were collected under the conventional cultivation of maize in Ganzhou district and Linze county of Zhangye City in Hexi oasis irrigation area, the characteristics and correlation of the contents of available zinc and salt ions as well as the effects of various ions on the content of available zinc were studied with the method of path analysis. The results showed that available zinc content ranged from 0.289 to 0.736 mg/kg with an average of 0.473 mg/kg and a variation coefficient of 18.68%. Available zinc contents of more than 70% of the studied soil samples were lower than the deficient level (0.5 mg/kg). Available zinc content had negative correlation with salt content, it was extremely significantly correlated with the contents of HCO– 3, Mg2+and Na+, and significantly correlated with the contents of Cl–, Ca2+and K+, but not significantly correlated with the content of SO2– 4. Mg2+showed strong direct and indirect negative effects on available zinc content, HCO– 3 showed more strong direct negative effects and strong indirect positive effects; Cl–, K+and Na+showed more strong direct and indirect negative effects, SO2– 4 showed mild direct and indirect positive effects, Ca2+showed mild direct and indirect negative effects. Mg2+, HCO– 3and Na+were the three dominant affecting factors on available zinc contents, and they in total decided 68% more of the determinative coefficient2.

    Saline fluvo-aquic soil; Available zinc contents; Salt ions; Correlation analysis; Path analysis

    10.13758/j.cnki.tr.2017.03.018

    S156.4

    A

    國家自然科學(xué)基金項目 (41261072)、農(nóng)業(yè)部公益性行業(yè)科研專項(200903001) 和甘肅省農(nóng)業(yè)科學(xué)院科技創(chuàng)新專項(2015GAAS03)資助。

    楊思存(1971—),男,甘肅靖遠(yuǎn)人,副研究員,主要從事土壤養(yǎng)分管理與鹽堿地改良利用研究。E-mail: yangsicun@sina.com

    猜你喜歡
    潮土負(fù)效應(yīng)缺鋅
    分析微信公眾平臺新聞傳播正負(fù)效應(yīng)
    新聞傳播(2018年7期)2018-12-06 08:57:42
    微信公眾平臺對新聞傳播產(chǎn)生的正負(fù)效應(yīng)
    新聞傳播(2018年12期)2018-09-19 06:26:52
    不同土壤類型小麥測土配方施肥對當(dāng)季土壤供氮的影響
    糖醇鋅葉面肥
    THE SYMMETRIC POSITIVE SOLUTIONS OF 2n-ORDER BOUNDARY VALUE PROBLEMS ON TIME SCALES??
    法制報道“負(fù)效應(yīng)”的規(guī)避與防范
    新聞傳播(2016年4期)2016-07-18 10:59:22
    親子(2016年5期)2016-05-14 15:15:52
    缺鋅脅迫對水稻葉綠體抗氧化系統(tǒng)的影響*
    長期施鉀和秸稈還田對河北潮土區(qū)作物產(chǎn)量和土壤鉀素狀況的影響
    典型潮土N2O排放的DNDC模型田間驗證研究
    男人和女人高潮做爰伦理| 久久人人爽人人片av| 亚洲精品成人av观看孕妇| 国产日韩欧美亚洲二区| 国产一区二区三区综合在线观看 | 色播亚洲综合网| 亚洲成人久久爱视频| 超碰97精品在线观看| 欧美97在线视频| 十八禁网站网址无遮挡 | 在线播放无遮挡| 亚洲成人av在线免费| 亚洲成人中文字幕在线播放| 听说在线观看完整版免费高清| 国产精品熟女久久久久浪| 99久久人妻综合| 一个人看视频在线观看www免费| 久久久亚洲精品成人影院| av免费在线看不卡| 午夜福利视频1000在线观看| 免费黄频网站在线观看国产| 国产午夜福利久久久久久| 亚洲成人中文字幕在线播放| 69av精品久久久久久| 日本黄色片子视频| 人人妻人人爽人人添夜夜欢视频 | 亚洲av中文字字幕乱码综合| 日韩人妻高清精品专区| 欧美国产精品一级二级三级 | 91在线精品国自产拍蜜月| 久久久久性生活片| 日韩成人伦理影院| 国产精品无大码| 久久久久久久午夜电影| 最近中文字幕高清免费大全6| 免费观看的影片在线观看| 久久精品国产鲁丝片午夜精品| 色综合色国产| 又大又黄又爽视频免费| 亚洲久久久久久中文字幕| 王馨瑶露胸无遮挡在线观看| 亚洲精品乱码久久久久久按摩| 国产免费福利视频在线观看| 亚洲精品久久久久久婷婷小说| 高清日韩中文字幕在线| 久久97久久精品| 韩国av在线不卡| 久久久a久久爽久久v久久| 中国美白少妇内射xxxbb| 亚洲av二区三区四区| 精品人妻视频免费看| 久久精品综合一区二区三区| 亚洲综合精品二区| 亚洲精品456在线播放app| 免费av观看视频| 免费观看在线日韩| 看非洲黑人一级黄片| 国产一区二区亚洲精品在线观看| 亚洲高清免费不卡视频| 97精品久久久久久久久久精品| 乱码一卡2卡4卡精品| 免费看不卡的av| 五月开心婷婷网| 国产成人精品福利久久| 久久久久国产网址| 久久久久久久久久久免费av| 日韩欧美 国产精品| 日本与韩国留学比较| 国产成人免费观看mmmm| 精品久久久久久电影网| 精品久久久久久久久av| 色视频在线一区二区三区| 国产黄片视频在线免费观看| 禁无遮挡网站| 午夜激情福利司机影院| 亚洲不卡免费看| av专区在线播放| 亚洲av日韩在线播放| 国产亚洲91精品色在线| 少妇熟女欧美另类| 99re6热这里在线精品视频| 国产有黄有色有爽视频| 2021少妇久久久久久久久久久| 自拍欧美九色日韩亚洲蝌蚪91 | 国产伦理片在线播放av一区| 大陆偷拍与自拍| 亚洲激情五月婷婷啪啪| 18禁裸乳无遮挡动漫免费视频 | 国产大屁股一区二区在线视频| 欧美变态另类bdsm刘玥| 久久鲁丝午夜福利片| 成人亚洲欧美一区二区av| 日韩av不卡免费在线播放| 美女主播在线视频| 免费电影在线观看免费观看| 在线观看三级黄色| 亚洲最大成人av| 久久精品国产亚洲网站| 日本熟妇午夜| 国产成人91sexporn| 国产在视频线精品| 国产成人午夜福利电影在线观看| 99热国产这里只有精品6| 国产精品久久久久久av不卡| 国产在线一区二区三区精| 精品少妇黑人巨大在线播放| 大话2 男鬼变身卡| 免费看a级黄色片| 日韩三级伦理在线观看| 午夜激情久久久久久久| 久久影院123| 人妻系列 视频| 久久亚洲国产成人精品v| 国产真实伦视频高清在线观看| 国产黄片视频在线免费观看| 亚洲人成网站高清观看| 欧美亚洲 丝袜 人妻 在线| 日韩成人av中文字幕在线观看| 亚洲精品国产av蜜桃| 亚洲经典国产精华液单| 国产成人精品婷婷| 青春草国产在线视频| 国产成人freesex在线| 草草在线视频免费看| 日韩免费高清中文字幕av| 国产色婷婷99| 日日啪夜夜撸| 王馨瑶露胸无遮挡在线观看| 联通29元200g的流量卡| 日本色播在线视频| 午夜老司机福利剧场| 国产成人a区在线观看| 男的添女的下面高潮视频| 国产视频首页在线观看| 啦啦啦啦在线视频资源| 久久久久久久久久久免费av| 欧美国产精品一级二级三级 | 丰满少妇做爰视频| 亚洲av中文字字幕乱码综合| 精品国产三级普通话版| 乱系列少妇在线播放| av在线亚洲专区| 精华霜和精华液先用哪个| 嫩草影院入口| 国产免费又黄又爽又色| 中文欧美无线码| 亚洲国产精品成人综合色| 国产综合懂色| 国产真实伦视频高清在线观看| 天天躁日日操中文字幕| 国产成人午夜福利电影在线观看| 国产精品久久久久久精品电影| 国产精品熟女久久久久浪| 久久久欧美国产精品| 日本欧美国产在线视频| 联通29元200g的流量卡| 少妇人妻 视频| 亚洲四区av| 国产精品国产av在线观看| 国产真实伦视频高清在线观看| 精品久久久久久久末码| 成人欧美大片| freevideosex欧美| 日韩免费高清中文字幕av| 波多野结衣巨乳人妻| av免费观看日本| videossex国产| 青春草视频在线免费观看| 自拍偷自拍亚洲精品老妇| 嫩草影院入口| av在线亚洲专区| 九九久久精品国产亚洲av麻豆| 国产老妇伦熟女老妇高清| 色5月婷婷丁香| 精品久久久久久电影网| 国产美女午夜福利| 在线a可以看的网站| 国产一区二区亚洲精品在线观看| 直男gayav资源| 国产免费一级a男人的天堂| 欧美一区二区亚洲| 老师上课跳d突然被开到最大视频| 日韩一本色道免费dvd| 国产av国产精品国产| 真实男女啪啪啪动态图| 精品人妻视频免费看| 日本一二三区视频观看| 免费大片黄手机在线观看| 国产成年人精品一区二区| 久久久久性生活片| 身体一侧抽搐| 嫩草影院新地址| 丝袜喷水一区| 亚洲熟女精品中文字幕| 亚洲欧美一区二区三区国产| 精品酒店卫生间| 午夜激情久久久久久久| 直男gayav资源| 男人舔奶头视频| 大片免费播放器 马上看| av专区在线播放| 少妇丰满av| 男人爽女人下面视频在线观看| 高清视频免费观看一区二区| 日本猛色少妇xxxxx猛交久久| 18禁动态无遮挡网站| 我的老师免费观看完整版| 嫩草影院新地址| 精品久久久噜噜| 亚洲欧美日韩无卡精品| 人妻少妇偷人精品九色| 全区人妻精品视频| 亚洲aⅴ乱码一区二区在线播放| 欧美97在线视频| 成人黄色视频免费在线看| 亚洲最大成人手机在线| 国产成人一区二区在线| 欧美国产精品一级二级三级 | 亚洲精品中文字幕在线视频 | 18禁动态无遮挡网站| 免费人成在线观看视频色| 在线亚洲精品国产二区图片欧美 | 久久久久国产网址| 97超视频在线观看视频| 性色avwww在线观看| 国产成人a区在线观看| 制服丝袜香蕉在线| 日韩大片免费观看网站| 高清欧美精品videossex| 久久精品国产鲁丝片午夜精品| 日韩成人伦理影院| 91午夜精品亚洲一区二区三区| 国产亚洲91精品色在线| 一本久久精品| 女人被狂操c到高潮| 亚洲,欧美,日韩| 亚洲精品亚洲一区二区| 亚洲欧洲日产国产| 五月天丁香电影| 成年女人在线观看亚洲视频 | 亚洲国产精品成人久久小说| 亚洲人成网站高清观看| 亚洲国产高清在线一区二区三| 国产亚洲91精品色在线| 丝袜脚勾引网站| 国产伦精品一区二区三区四那| 日本午夜av视频| 熟女人妻精品中文字幕| 亚洲电影在线观看av| 天美传媒精品一区二区| 国产淫语在线视频| 97超视频在线观看视频| 尤物成人国产欧美一区二区三区| 婷婷色综合大香蕉| 超碰97精品在线观看| 国产高潮美女av| 一个人观看的视频www高清免费观看| 午夜视频国产福利| 五月玫瑰六月丁香| 中文在线观看免费www的网站| 亚洲天堂av无毛| 国产成人aa在线观看| 精品国产露脸久久av麻豆| 偷拍熟女少妇极品色| 国产成人a区在线观看| 最新中文字幕久久久久| 特大巨黑吊av在线直播| 精品人妻视频免费看| 熟女人妻精品中文字幕| 纵有疾风起免费观看全集完整版| 午夜福利视频精品| 亚洲人成网站高清观看| 免费av不卡在线播放| 少妇的逼水好多| 我的老师免费观看完整版| 一级a做视频免费观看| 午夜爱爱视频在线播放| 亚洲国产高清在线一区二区三| 69人妻影院| 小蜜桃在线观看免费完整版高清| 高清av免费在线| 亚洲天堂国产精品一区在线| 搡女人真爽免费视频火全软件| 青春草国产在线视频| 日韩欧美 国产精品| 亚洲va在线va天堂va国产| 日韩,欧美,国产一区二区三区| 夫妻午夜视频| 青春草视频在线免费观看| 亚洲av.av天堂| 晚上一个人看的免费电影| 国产黄色免费在线视频| 国产中年淑女户外野战色| 免费观看av网站的网址| 插阴视频在线观看视频| 欧美日韩综合久久久久久| 久久99热这里只频精品6学生| 午夜视频国产福利| 成人特级av手机在线观看| 国产一区有黄有色的免费视频| 日韩伦理黄色片| 国产69精品久久久久777片| 亚洲不卡免费看| 日韩一本色道免费dvd| 亚洲四区av| 网址你懂的国产日韩在线| 亚洲av一区综合| 亚洲av中文字字幕乱码综合| 夫妻午夜视频| 欧美变态另类bdsm刘玥| 激情五月婷婷亚洲| 赤兔流量卡办理| 中文精品一卡2卡3卡4更新| 欧美 日韩 精品 国产| 国产黄色免费在线视频| 麻豆乱淫一区二区| 汤姆久久久久久久影院中文字幕| 午夜激情福利司机影院| 神马国产精品三级电影在线观看| 狂野欧美白嫩少妇大欣赏| 纵有疾风起免费观看全集完整版| 麻豆成人午夜福利视频| 永久免费av网站大全| 国产免费一区二区三区四区乱码| 亚洲av二区三区四区| 国产av国产精品国产| 国产精品国产av在线观看| 亚洲内射少妇av| 99热这里只有是精品50| 99久久中文字幕三级久久日本| 亚洲不卡免费看| 26uuu在线亚洲综合色| 又爽又黄无遮挡网站| 国产亚洲5aaaaa淫片| 欧美xxxx性猛交bbbb| 成人毛片a级毛片在线播放| 午夜福利视频1000在线观看| 天堂中文最新版在线下载 | tube8黄色片| 国产精品av视频在线免费观看| 国产一级毛片在线| 欧美日韩精品成人综合77777| 美女高潮的动态| 久久久午夜欧美精品| 97精品久久久久久久久久精品| 久久久久久九九精品二区国产| 国产片特级美女逼逼视频| 国产精品久久久久久久电影| 国产精品熟女久久久久浪| 欧美三级亚洲精品| 成年人午夜在线观看视频| av.在线天堂| 午夜老司机福利剧场| av天堂中文字幕网| 插阴视频在线观看视频| 亚洲av日韩在线播放| 99久国产av精品国产电影| 高清午夜精品一区二区三区| 国产亚洲av嫩草精品影院| 最近中文字幕2019免费版| 嫩草影院入口| 夜夜看夜夜爽夜夜摸| 精品久久久精品久久久| 亚洲第一区二区三区不卡| 毛片女人毛片| 高清午夜精品一区二区三区| 老司机影院成人| 国内精品美女久久久久久| 日本欧美国产在线视频| 欧美精品人与动牲交sv欧美| 99热全是精品| 亚洲熟女精品中文字幕| 熟妇人妻不卡中文字幕| 高清日韩中文字幕在线| 欧美一级a爱片免费观看看| 国产一区亚洲一区在线观看| 免费人成在线观看视频色| 国产男女内射视频| 国产又色又爽无遮挡免| 在线观看美女被高潮喷水网站| 日日摸夜夜添夜夜添av毛片| 校园人妻丝袜中文字幕| 国产伦理片在线播放av一区| 亚洲自偷自拍三级| 亚洲欧美中文字幕日韩二区| av女优亚洲男人天堂| 亚洲在久久综合| 国产男人的电影天堂91| 99热这里只有是精品在线观看| 伦理电影大哥的女人| 97在线视频观看| 在线亚洲精品国产二区图片欧美 | 日本与韩国留学比较| 美女内射精品一级片tv| 人妻 亚洲 视频| 亚洲欧洲国产日韩| 水蜜桃什么品种好| 黄色一级大片看看| 国产精品人妻久久久影院| 欧美成人精品欧美一级黄| 国产精品一区二区性色av| 午夜免费鲁丝| 国产淫语在线视频| 亚洲精品成人久久久久久| 日韩欧美精品免费久久| 看免费成人av毛片| 青春草亚洲视频在线观看| 男人添女人高潮全过程视频| 国产成人91sexporn| 国产亚洲av片在线观看秒播厂| 一级毛片 在线播放| 免费少妇av软件| 夫妻午夜视频| 在线精品无人区一区二区三 | 亚洲不卡免费看| 国产精品一及| 一级爰片在线观看| 亚洲精品国产成人久久av| 尾随美女入室| 毛片一级片免费看久久久久| 在线天堂最新版资源| 久久久色成人| 亚洲精品日韩在线中文字幕| 亚洲成人一二三区av| eeuss影院久久| 久久精品国产鲁丝片午夜精品| 一本色道久久久久久精品综合| 国产免费又黄又爽又色| 搡老乐熟女国产| 免费看光身美女| 国产精品久久久久久av不卡| 国产成人精品久久久久久| 亚洲欧美一区二区三区国产| 亚洲熟女精品中文字幕| 精品人妻一区二区三区麻豆| 少妇的逼好多水| 欧美老熟妇乱子伦牲交| 在线观看免费高清a一片| 99精国产麻豆久久婷婷| 国产精品一区www在线观看| 2022亚洲国产成人精品| 久久久久久久国产电影| 一级毛片 在线播放| 国产淫语在线视频| 综合色丁香网| 美女脱内裤让男人舔精品视频| 成年人午夜在线观看视频| 2022亚洲国产成人精品| 成年av动漫网址| 爱豆传媒免费全集在线观看| kizo精华| 性插视频无遮挡在线免费观看| 插逼视频在线观看| 午夜福利视频1000在线观看| 亚洲成人av在线免费| 五月开心婷婷网| 18禁裸乳无遮挡免费网站照片| 亚洲精品aⅴ在线观看| 99热6这里只有精品| av在线观看视频网站免费| 午夜精品国产一区二区电影 | 在线免费十八禁| 国产黄频视频在线观看| 亚洲精品一区蜜桃| 精品久久久久久久久av| 免费看光身美女| 18禁裸乳无遮挡免费网站照片| 久久久久久九九精品二区国产| 中文欧美无线码| 直男gayav资源| 国产av不卡久久| 日韩在线高清观看一区二区三区| 国产精品久久久久久av不卡| 成人毛片a级毛片在线播放| 欧美丝袜亚洲另类| 日韩伦理黄色片| 国产精品麻豆人妻色哟哟久久| 色哟哟·www| 欧美+日韩+精品| 日韩欧美精品免费久久| 我的老师免费观看完整版| 建设人人有责人人尽责人人享有的 | 91精品国产九色| 七月丁香在线播放| 波多野结衣巨乳人妻| 女人十人毛片免费观看3o分钟| 一级毛片黄色毛片免费观看视频| 免费在线观看成人毛片| 精品国产三级普通话版| 成人欧美大片| 日日摸夜夜添夜夜爱| 寂寞人妻少妇视频99o| 青春草国产在线视频| 亚洲精品日韩av片在线观看| tube8黄色片| 舔av片在线| 久久鲁丝午夜福利片| 国内精品宾馆在线| av播播在线观看一区| 中文字幕人妻熟人妻熟丝袜美| 九色成人免费人妻av| 内射极品少妇av片p| 婷婷色综合www| av在线观看视频网站免费| 亚洲精品一区蜜桃| 综合色av麻豆| 国产高清不卡午夜福利| 秋霞在线观看毛片| 久久久久久伊人网av| 国产又色又爽无遮挡免| 亚洲婷婷狠狠爱综合网| 男人舔奶头视频| 亚洲精品乱码久久久v下载方式| 久久人人爽人人爽人人片va| 久久鲁丝午夜福利片| 在线 av 中文字幕| 欧美日韩精品成人综合77777| 亚洲在线观看片| 黄色怎么调成土黄色| 夜夜看夜夜爽夜夜摸| 性色av一级| 七月丁香在线播放| 天堂网av新在线| 日韩av免费高清视频| 久久这里有精品视频免费| 欧美亚洲 丝袜 人妻 在线| 舔av片在线| 国产综合懂色| 一本色道久久久久久精品综合| 亚洲人成网站高清观看| 国产黄a三级三级三级人| 18+在线观看网站| 视频中文字幕在线观看| 国产欧美亚洲国产| 亚洲婷婷狠狠爱综合网| 亚洲欧美精品自产自拍| 亚洲欧美中文字幕日韩二区| 久久久午夜欧美精品| 亚洲精品中文字幕在线视频 | 日韩一区二区三区影片| 欧美日韩精品成人综合77777| 久久精品国产a三级三级三级| 国产欧美另类精品又又久久亚洲欧美| 午夜免费观看性视频| 一级毛片久久久久久久久女| 看十八女毛片水多多多| 美女cb高潮喷水在线观看| 男的添女的下面高潮视频| 日本黄色片子视频| 日韩不卡一区二区三区视频在线| 久久久久久久久久久丰满| 成人国产av品久久久| 亚洲av福利一区| 青春草视频在线免费观看| 国产视频首页在线观看| 国产精品不卡视频一区二区| 国产综合精华液| 精品酒店卫生间| 18禁裸乳无遮挡免费网站照片| 国产午夜福利久久久久久| 国产精品一区www在线观看| 免费看a级黄色片| 亚洲av不卡在线观看| 搡女人真爽免费视频火全软件| 在线观看人妻少妇| 我的女老师完整版在线观看| 精品视频人人做人人爽| 日韩强制内射视频| 十八禁网站网址无遮挡 | 亚洲精品久久久久久婷婷小说| 免费看av在线观看网站| 亚洲av成人精品一二三区| 99热这里只有是精品50| 欧美少妇被猛烈插入视频| 精品少妇黑人巨大在线播放| 欧美成人a在线观看| 美女内射精品一级片tv| 91aial.com中文字幕在线观看| 国产精品一二三区在线看| 亚洲精品,欧美精品| 久久精品久久久久久久性| 日韩中字成人| 极品少妇高潮喷水抽搐| 97在线视频观看| 国产成人精品福利久久| 亚洲四区av| 2021天堂中文幕一二区在线观| 内地一区二区视频在线| 少妇被粗大猛烈的视频| 真实男女啪啪啪动态图| 激情 狠狠 欧美| 欧美人与善性xxx| 91精品伊人久久大香线蕉| 99精国产麻豆久久婷婷| 久久97久久精品| 肉色欧美久久久久久久蜜桃 | 国产淫语在线视频| 校园人妻丝袜中文字幕| 波野结衣二区三区在线| 成人无遮挡网站| 22中文网久久字幕| 最近的中文字幕免费完整| 国产精品国产三级专区第一集| 三级国产精品片| 国产一区二区在线观看日韩| 观看免费一级毛片| 色哟哟·www| 欧美人与善性xxx| 亚洲精品影视一区二区三区av| 免费看光身美女| 五月天丁香电影| 成人黄色视频免费在线看| 极品少妇高潮喷水抽搐| 国产黄a三级三级三级人| 日本一本二区三区精品|