• <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网站在线观看| 国产av一区二区精品久久| 五月开心婷婷网| 高清黄色对白视频在线免费看| 黄片无遮挡物在线观看| 日韩制服骚丝袜av| 欧美国产精品一级二级三级| 国产片内射在线| 国产av国产精品国产| 国产精品国产三级国产专区5o| tube8黄色片| 日本黄色日本黄色录像| 哪个播放器可以免费观看大片| 激情五月婷婷亚洲| 水蜜桃什么品种好| 国产午夜精品一二区理论片| freevideosex欧美| 国产日韩一区二区三区精品不卡| av女优亚洲男人天堂| tube8黄色片| 精品国产乱码久久久久久男人| 99热国产这里只有精品6| 最近中文字幕高清免费大全6| 免费黄色在线免费观看| 尾随美女入室| 一区在线观看完整版| 青春草亚洲视频在线观看| 午夜激情av网站| 欧美在线黄色| 欧美精品一区二区大全| 激情五月婷婷亚洲| 97在线视频观看| 中文乱码字字幕精品一区二区三区| 亚洲图色成人| 夫妻性生交免费视频一级片| 成人国产麻豆网| 久久99热这里只频精品6学生| 人妻系列 视频| 国产精品国产av在线观看| 精品福利永久在线观看| 国产亚洲最大av| 午夜激情av网站| 一区二区日韩欧美中文字幕| 日产精品乱码卡一卡2卡三| 97在线视频观看| 九色亚洲精品在线播放| 亚洲欧美成人精品一区二区| 一级,二级,三级黄色视频| 国产精品欧美亚洲77777| 男人爽女人下面视频在线观看| 国产精品久久久av美女十八| 老司机影院毛片| 午夜av观看不卡| 久久人人爽av亚洲精品天堂| 18禁观看日本| 亚洲成色77777| 日韩不卡一区二区三区视频在线| 黄色怎么调成土黄色| 2021少妇久久久久久久久久久| 国产片特级美女逼逼视频| 九色亚洲精品在线播放| 亚洲久久久国产精品| 亚洲人成电影观看| 亚洲 欧美一区二区三区| 成人二区视频| 久久久国产精品麻豆| 亚洲天堂av无毛| xxx大片免费视频| 久久久久久久久免费视频了| 三上悠亚av全集在线观看| 日本黄色日本黄色录像| 日韩中字成人| 亚洲综合精品二区| 一区福利在线观看| 汤姆久久久久久久影院中文字幕| 黄片播放在线免费| 亚洲国产欧美日韩在线播放| 波野结衣二区三区在线| 成人国产av品久久久| 人人澡人人妻人| 妹子高潮喷水视频| 高清视频免费观看一区二区| 久久99蜜桃精品久久| 成人漫画全彩无遮挡| 国产精品国产三级国产专区5o| 国产成人一区二区在线| 99香蕉大伊视频| 一级,二级,三级黄色视频| 久久99蜜桃精品久久| 欧美精品国产亚洲| 观看美女的网站| 少妇猛男粗大的猛烈进出视频| 亚洲精品国产一区二区精华液| 男女国产视频网站| 国产男女内射视频| 国产成人免费无遮挡视频| a级毛片黄视频| 日韩欧美精品免费久久| 男人爽女人下面视频在线观看| 韩国av在线不卡| 欧美成人精品欧美一级黄| 亚洲国产成人一精品久久久| 亚洲美女搞黄在线观看| 一级爰片在线观看| 伦理电影大哥的女人| 欧美精品一区二区免费开放| 亚洲婷婷狠狠爱综合网| 精品一区二区三卡| 亚洲精品视频女| 成人漫画全彩无遮挡| 国产精品久久久久久久久免| 大陆偷拍与自拍| av又黄又爽大尺度在线免费看| av女优亚洲男人天堂| tube8黄色片| 99国产综合亚洲精品| 菩萨蛮人人尽说江南好唐韦庄| 午夜老司机福利剧场| 一级爰片在线观看| 国产免费福利视频在线观看| 男女高潮啪啪啪动态图| 婷婷成人精品国产| 99国产综合亚洲精品| 爱豆传媒免费全集在线观看| 一级a爱视频在线免费观看| 热99久久久久精品小说推荐| 欧美精品一区二区免费开放| 精品卡一卡二卡四卡免费| 国产日韩欧美在线精品| 免费人妻精品一区二区三区视频| 大片免费播放器 马上看| 午夜福利视频在线观看免费| 国产黄频视频在线观看| 亚洲国产最新在线播放| 久久毛片免费看一区二区三区| 久久精品亚洲av国产电影网| 日韩免费高清中文字幕av| 午夜福利乱码中文字幕| 两性夫妻黄色片| 最近中文字幕高清免费大全6| 亚洲av在线观看美女高潮| 免费观看av网站的网址| 欧美精品av麻豆av| 久久精品国产a三级三级三级| 肉色欧美久久久久久久蜜桃| 精品第一国产精品| 久久精品国产亚洲av天美| 亚洲国产精品一区二区三区在线| 国产精品女同一区二区软件| 伊人久久大香线蕉亚洲五| 精品国产乱码久久久久久男人| 欧美少妇被猛烈插入视频| 热re99久久精品国产66热6| 五月开心婷婷网| 成年女人毛片免费观看观看9 | 国产亚洲午夜精品一区二区久久| 中文字幕亚洲精品专区| av.在线天堂| 激情五月婷婷亚洲| 国产野战对白在线观看| 1024视频免费在线观看| 国产成人免费无遮挡视频| 日本91视频免费播放| 人成视频在线观看免费观看| 精品第一国产精品| 美女午夜性视频免费| 亚洲人成电影观看| 亚洲欧美一区二区三区久久| 宅男免费午夜| 亚洲国产av影院在线观看| 亚洲av免费高清在线观看| 热99久久久久精品小说推荐| 一本色道久久久久久精品综合| 九色亚洲精品在线播放| videosex国产| 亚洲精品国产色婷婷电影| 在线观看www视频免费| 国产av精品麻豆| 多毛熟女@视频| 丝袜在线中文字幕| 大话2 男鬼变身卡| 黄色毛片三级朝国网站| 美女中出高潮动态图| 在线观看免费日韩欧美大片| 91精品国产国语对白视频| 男女无遮挡免费网站观看| 免费人妻精品一区二区三区视频| 国产精品嫩草影院av在线观看| 宅男免费午夜| 中文字幕色久视频| 欧美亚洲日本最大视频资源| 国产女主播在线喷水免费视频网站| 人人妻人人添人人爽欧美一区卜| 美女中出高潮动态图| 亚洲熟女精品中文字幕| 国产精品熟女久久久久浪| 欧美老熟妇乱子伦牲交| 日韩中文字幕视频在线看片| 免费黄网站久久成人精品| 永久免费av网站大全| 综合色丁香网| 精品一品国产午夜福利视频| 十分钟在线观看高清视频www| 三级国产精品片| 国产熟女欧美一区二区| 国产一区二区 视频在线| 狠狠精品人妻久久久久久综合| 中文字幕另类日韩欧美亚洲嫩草| 国产男女内射视频| 亚洲国产精品国产精品| 日产精品乱码卡一卡2卡三| 99九九在线精品视频| 成年女人在线观看亚洲视频| 亚洲国产av影院在线观看| 久久精品夜色国产| 久久99一区二区三区| 国产av国产精品国产| 美女午夜性视频免费| 欧美日韩精品成人综合77777| 亚洲国产色片| 国产白丝娇喘喷水9色精品| 国产精品免费大片| 男女午夜视频在线观看| 国产一区二区三区av在线| 最近手机中文字幕大全| 免费大片黄手机在线观看| 午夜久久久在线观看| 自线自在国产av| 国产精品 国内视频| 午夜免费鲁丝| 国产色婷婷99| 色吧在线观看| 国产深夜福利视频在线观看| 亚洲精品久久久久久婷婷小说| 好男人视频免费观看在线| 少妇的逼水好多| 久久精品久久久久久噜噜老黄| 曰老女人黄片| 熟女电影av网| 精品少妇内射三级| 亚洲,欧美,日韩| 国产一区亚洲一区在线观看| 亚洲美女搞黄在线观看| 婷婷色综合大香蕉| 秋霞在线观看毛片| 亚洲一码二码三码区别大吗| 亚洲欧美一区二区三区久久| 人人妻人人爽人人添夜夜欢视频| 91久久精品国产一区二区三区| 久久久a久久爽久久v久久| 亚洲美女黄色视频免费看| 高清欧美精品videossex| 亚洲在久久综合| 一级片'在线观看视频| 狠狠精品人妻久久久久久综合| 伊人亚洲综合成人网| 精品少妇内射三级| 国产探花极品一区二区| 国产人伦9x9x在线观看 | 国产精品三级大全| 老女人水多毛片| 国产精品偷伦视频观看了| 国产老妇伦熟女老妇高清| av有码第一页| 人人妻人人添人人爽欧美一区卜| 欧美老熟妇乱子伦牲交| 国产亚洲最大av| 在线观看一区二区三区激情| 日韩一卡2卡3卡4卡2021年| 亚洲av综合色区一区| 亚洲精品aⅴ在线观看| 色哟哟·www| 久久人人97超碰香蕉20202| 国产精品偷伦视频观看了| 国产精品久久久久久精品电影小说| 免费高清在线观看日韩| 热re99久久国产66热| 熟女av电影| 在线观看免费高清a一片| 老女人水多毛片| 国产欧美日韩综合在线一区二区| 国产精品国产三级国产专区5o| 亚洲国产精品一区二区三区在线| 大片免费播放器 马上看| 国产精品二区激情视频| 国产综合精华液| 成年av动漫网址| 黄片无遮挡物在线观看| 亚洲精品第二区| 精品国产露脸久久av麻豆| 午夜免费鲁丝| 十八禁网站网址无遮挡| 人人妻人人添人人爽欧美一区卜| 人人妻人人爽人人添夜夜欢视频| 久久热在线av| 亚洲精品一二三| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲美女搞黄在线观看| 人人妻人人澡人人看| 麻豆精品久久久久久蜜桃| 高清黄色对白视频在线免费看| 午夜91福利影院| 免费大片黄手机在线观看| 激情五月婷婷亚洲| 热99久久久久精品小说推荐| 伦理电影免费视频| 亚洲视频免费观看视频| 男女高潮啪啪啪动态图| 久久久欧美国产精品| 男人操女人黄网站| 少妇 在线观看| 亚洲国产精品一区二区三区在线| 亚洲国产欧美网| 国产一区二区三区综合在线观看| 国产精品女同一区二区软件| 久久亚洲国产成人精品v| 久久久国产一区二区| 精品亚洲乱码少妇综合久久| 日韩精品有码人妻一区| 国产无遮挡羞羞视频在线观看| 边亲边吃奶的免费视频| 在线 av 中文字幕| 中文天堂在线官网| 久久久久精品性色| 久久久亚洲精品成人影院| 日本欧美国产在线视频| 久久久国产精品麻豆| 一区二区三区四区激情视频| 亚洲精华国产精华液的使用体验| 91久久精品国产一区二区三区| 久久婷婷青草| 97在线视频观看| 精品视频人人做人人爽| 男女边吃奶边做爰视频| 久久这里有精品视频免费| 日韩三级伦理在线观看| 欧美97在线视频| 成人国语在线视频| 丝袜喷水一区| 一本色道久久久久久精品综合| 十八禁高潮呻吟视频| 日韩精品有码人妻一区| 日日啪夜夜爽| 亚洲精品乱久久久久久| 天天躁狠狠躁夜夜躁狠狠躁| 一本色道久久久久久精品综合| 精品亚洲成a人片在线观看| 精品久久久久久电影网| 老熟女久久久| 日韩大片免费观看网站| 欧美在线黄色| 夫妻午夜视频| 另类亚洲欧美激情| 久久午夜综合久久蜜桃| 亚洲久久久国产精品| 人人澡人人妻人| 亚洲国产精品999| 成人漫画全彩无遮挡| a级片在线免费高清观看视频| 成人影院久久| 少妇被粗大的猛进出69影院| 欧美精品国产亚洲| 久久婷婷青草| 久久久久人妻精品一区果冻| 99国产精品免费福利视频| 丰满迷人的少妇在线观看| 精品一品国产午夜福利视频| 黑人猛操日本美女一级片| 精品视频人人做人人爽| 91久久精品国产一区二区三区| 国产午夜精品一二区理论片| 9色porny在线观看| 国产av一区二区精品久久| 成人国产麻豆网| 亚洲欧美日韩另类电影网站| 中文字幕另类日韩欧美亚洲嫩草| 日本欧美国产在线视频| av网站免费在线观看视频| 欧美日韩av久久| 大话2 男鬼变身卡| 欧美精品亚洲一区二区| 男女边摸边吃奶| 人妻少妇偷人精品九色| 久久免费观看电影| 国产精品国产三级国产专区5o| 免费人妻精品一区二区三区视频| 香蕉精品网在线| 日本色播在线视频| 乱人伦中国视频| 日日撸夜夜添| 国产成人精品一,二区| 一级片'在线观看视频| 成年女人在线观看亚洲视频| 中国国产av一级| 宅男免费午夜| 久久精品久久精品一区二区三区| 女人久久www免费人成看片| 毛片一级片免费看久久久久| 午夜福利一区二区在线看| 一级毛片 在线播放| 午夜福利视频在线观看免费| 有码 亚洲区| 精品人妻在线不人妻| 免费av中文字幕在线| 80岁老熟妇乱子伦牲交| 欧美日韩视频高清一区二区三区二| 伦精品一区二区三区| 啦啦啦中文免费视频观看日本| 成人免费观看视频高清| 精品久久久精品久久久| 人成视频在线观看免费观看| 日韩中文字幕欧美一区二区 | 欧美激情高清一区二区三区 | 考比视频在线观看| 欧美日韩视频高清一区二区三区二| 欧美激情 高清一区二区三区| 一本色道久久久久久精品综合| 性色av一级| 欧美xxⅹ黑人| 国产成人a∨麻豆精品| 热99久久久久精品小说推荐| 亚洲一级一片aⅴ在线观看| 香蕉精品网在线| 最新的欧美精品一区二区| 国产在视频线精品| 国产精品.久久久| 欧美少妇被猛烈插入视频| 精品人妻一区二区三区麻豆| 麻豆av在线久日| 一本大道久久a久久精品| 日韩,欧美,国产一区二区三区| 夫妻午夜视频| 久久久久久久亚洲中文字幕| 日本wwww免费看| 久久久久久免费高清国产稀缺| 国产成人精品无人区| 咕卡用的链子| 精品人妻熟女毛片av久久网站| 国产亚洲午夜精品一区二区久久| 久久久国产欧美日韩av| 亚洲成人手机| 一区二区日韩欧美中文字幕| 一本色道久久久久久精品综合| 九九爱精品视频在线观看| 久久99精品国语久久久| 人妻少妇偷人精品九色| 免费观看在线日韩| 国产免费一区二区三区四区乱码| 久久精品国产亚洲av高清一级| 久久 成人 亚洲| 国产探花极品一区二区| 制服人妻中文乱码| 亚洲三级黄色毛片| 免费不卡的大黄色大毛片视频在线观看| 亚洲美女搞黄在线观看| av免费观看日本| 99热全是精品| 一级,二级,三级黄色视频| 亚洲一区中文字幕在线| 精品国产一区二区三区久久久樱花| 午夜福利视频精品| 国产白丝娇喘喷水9色精品| 一区二区三区四区激情视频| 久久精品国产亚洲av涩爱| 丰满乱子伦码专区| 久久久精品免费免费高清| 97在线视频观看| 精品国产一区二区三区四区第35| 日韩制服骚丝袜av| 三上悠亚av全集在线观看| 国产一区二区三区av在线| 国产精品国产三级国产专区5o| 日韩人妻精品一区2区三区| 免费在线观看视频国产中文字幕亚洲 | 少妇精品久久久久久久| 精品国产乱码久久久久久小说| 亚洲精品在线美女| 大香蕉久久成人网| 麻豆精品久久久久久蜜桃| 国产精品国产三级专区第一集| 亚洲精品第二区| 26uuu在线亚洲综合色| 老女人水多毛片| 国产亚洲最大av| 国产精品.久久久| 在现免费观看毛片| 制服诱惑二区| 国产一区有黄有色的免费视频| 免费av中文字幕在线| 久久久久国产网址| 国产片特级美女逼逼视频| 国产一区二区三区av在线| 91成人精品电影| 国产人伦9x9x在线观看 | 国产精品久久久久久精品电影小说| 久久99蜜桃精品久久| 伊人久久国产一区二区| 精品一区在线观看国产| 人人妻人人澡人人爽人人夜夜| 曰老女人黄片| 最黄视频免费看| 亚洲国产毛片av蜜桃av| 久久久久人妻精品一区果冻| 亚洲成人av在线免费| 欧美日韩精品网址| 日韩精品有码人妻一区| 一区二区日韩欧美中文字幕| 少妇猛男粗大的猛烈进出视频| 纵有疾风起免费观看全集完整版| 女人被躁到高潮嗷嗷叫费观| 大片电影免费在线观看免费| 国产欧美亚洲国产| 日本欧美国产在线视频| 一区在线观看完整版| 午夜影院在线不卡| 亚洲精品久久午夜乱码| av.在线天堂| 久久久精品区二区三区| 各种免费的搞黄视频| 宅男免费午夜| 日韩免费高清中文字幕av| 一区二区三区精品91| 乱人伦中国视频| 亚洲精品国产色婷婷电影| 午夜福利在线免费观看网站| 欧美日韩精品网址| 亚洲美女视频黄频| 午夜免费男女啪啪视频观看| 亚洲欧美一区二区三区久久| av在线观看视频网站免费| 亚洲精品一区蜜桃| www.熟女人妻精品国产| 成人免费观看视频高清| 少妇 在线观看| 日韩一卡2卡3卡4卡2021年| 亚洲视频免费观看视频| 天天躁日日躁夜夜躁夜夜| 超色免费av| 一个人免费看片子| xxxhd国产人妻xxx| 中文天堂在线官网| 久久精品国产亚洲av天美| av一本久久久久| 色94色欧美一区二区| 亚洲精品第二区| 桃花免费在线播放| 99re6热这里在线精品视频| 国产精品熟女久久久久浪| 五月伊人婷婷丁香| 午夜免费观看性视频| 亚洲成人手机| av在线观看视频网站免费| 波多野结衣av一区二区av| 午夜免费鲁丝| 交换朋友夫妻互换小说| 中文字幕制服av| 欧美国产精品一级二级三级| 久久久久久久久久人人人人人人| 建设人人有责人人尽责人人享有的| 一本色道久久久久久精品综合| 久久精品国产亚洲av涩爱| 亚洲一区中文字幕在线| 亚洲美女视频黄频| 亚洲,欧美精品.| 婷婷色综合www| 哪个播放器可以免费观看大片| 成年女人在线观看亚洲视频| 毛片一级片免费看久久久久| 国产极品天堂在线| 久久久国产欧美日韩av| 久久青草综合色| 日韩人妻精品一区2区三区| 日日啪夜夜爽| 色视频在线一区二区三区| 国产成人精品久久久久久| 精品99又大又爽又粗少妇毛片| 伦理电影免费视频| av.在线天堂| 在线亚洲精品国产二区图片欧美| 男女高潮啪啪啪动态图| 欧美激情高清一区二区三区 | 成人免费观看视频高清| 欧美成人午夜精品| 日韩大片免费观看网站| 日韩,欧美,国产一区二区三区| 国产精品99久久99久久久不卡 | 一区二区三区乱码不卡18| av在线观看视频网站免费| 欧美在线黄色| 精品卡一卡二卡四卡免费| 亚洲美女黄色视频免费看| 成年动漫av网址| 日韩欧美精品免费久久| av免费观看日本| 午夜免费男女啪啪视频观看| 久久国内精品自在自线图片| 欧美人与性动交α欧美精品济南到 | 精品福利永久在线观看| 精品亚洲乱码少妇综合久久| 天堂8中文在线网| 老鸭窝网址在线观看| 久久精品国产鲁丝片午夜精品| 色视频在线一区二区三区| 日韩人妻精品一区2区三区| 国产亚洲av片在线观看秒播厂| 一区二区三区四区激情视频| 这个男人来自地球电影免费观看 | 精品福利永久在线观看| 极品少妇高潮喷水抽搐| 亚洲色图综合在线观看| 亚洲精品久久久久久婷婷小说|