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

    九龍坡花椒種植區(qū)地形、土壤肥力與花椒產(chǎn)量的關(guān)系*

    2019-12-11 07:41:04楊仕曦呂廣斌向華輝王正銀
    關(guān)鍵詞:種植區(qū)土壤肥力產(chǎn)區(qū)

    楊仕曦, 呂廣斌, 黃 云**, 向華輝, 王正銀

    九龍坡花椒種植區(qū)地形、土壤肥力與花椒產(chǎn)量的關(guān)系*

    楊仕曦1, 呂廣斌1, 黃 云1**, 向華輝2, 王正銀1

    (1. 西南大學(xué)資源環(huán)境學(xué)院 重慶 400716; 2. 重慶市九龍坡區(qū)農(nóng)業(yè)技術(shù)推廣中心 重慶 400039)

    為了解九龍坡花椒種植區(qū)土壤養(yǎng)分狀況及該區(qū)地形因子、土壤肥力因子與花椒產(chǎn)量的關(guān)系, 為科學(xué)合理制定花椒高效施肥措施提供理論依據(jù), 本研究采用田間調(diào)查研究和室內(nèi)分析的方法, 研究了九龍坡花椒種植區(qū)低、中、高產(chǎn)區(qū)的海拔、坡度及土壤pH、有機(jī)質(zhì)、大量微量元素含量和交換性能的變化特征, 及其與花椒產(chǎn)量的關(guān)系。研究結(jié)果表明: 九龍坡花椒普遍種植于200~500 m海拔范圍, 高產(chǎn)區(qū)集中在300 m左右的海拔;從低產(chǎn)區(qū)到高產(chǎn)區(qū)坡度略有增加, 但未達(dá)顯著水平。土壤均屬酸性土, pH<6.5。土壤肥力總體屬高水平范圍, 但各養(yǎng)分因子差異很大, 其中土壤陽離子交換量(CEC)、有效磷、有效鈣、有效鎂、有效鐵、有效錳、有效銅、有效鋅含量豐富,分別為27.2 cmol(+)?kg-1、35.2 mg?kg-1、3 289.8 mg?kg-1、271.8 mg?kg-1、48.6 mg?kg-1、62.1 mg?kg-1、1.5 mg?kg-1、4.5 mg?kg-1; 有機(jī)質(zhì)、堿解氮、速效鉀、交換性酸屬適中水平,分別為19.1 mg?kg-1、114.9 mg?kg-1、107.0 mg?kg-1、8.1 cmol(+)?kg-1; 水溶性硼缺乏,為0.28 mg?kg-1。相關(guān)分析表明花椒產(chǎn)量與有效鈣、CEC、pH、有效錳、水溶性硼呈顯著正相關(guān); 通徑分析結(jié)果表明有效鈣、CEC、交換性酸、有效銅、有效鐵、有效鋅是影響花椒產(chǎn)量的主要因子,逐步回歸分析構(gòu)建了有效鈣(6)與花椒產(chǎn)量()的最優(yōu)回歸線性方程:=11.693+0.0036。綜上所述, 九龍坡花椒種植區(qū)土壤養(yǎng)分失衡較為嚴(yán)重,施肥應(yīng)注重養(yǎng)分的平衡,增施有機(jī)肥,改善土壤理化性狀,治理土壤酸化。

    花椒; 地形; 土壤肥力; 產(chǎn)量; 養(yǎng)分平衡

    花椒()為蕓香科(Rutaceae)花椒屬植物, 在我國已有兩千多年的栽培歷史; 花椒是我國特有的食用香辛料樹種, 廣泛種植于四川、重慶、甘肅、陜西等省(市), 其中以四川、重慶的花椒最為著名[1-2]。花椒根系發(fā)達(dá), 耐寒、耐干旱瘠薄[3], 是我國極力推廣的一種生態(tài)經(jīng)濟(jì)型樹種, 在促進(jìn)山區(qū)經(jīng)濟(jì)發(fā)展中具有重要作用[4]。近年來, 隨著花椒產(chǎn)業(yè)不斷發(fā)展, 農(nóng)戶為提高產(chǎn)量, 化學(xué)肥料用量不斷增加, 有機(jī)肥用量不足, 造成土壤養(yǎng)分失調(diào)、土壤酸化、土壤肥力下降的現(xiàn)象日益嚴(yán)重[5]。

    地形和土壤肥力是影響作物產(chǎn)量的兩個重要因素[6-7]。研究表明, 土壤養(yǎng)分含量的分布往往與地形的空間變異性有關(guān), 且獨(dú)特地形形成的小氣候也能對作物產(chǎn)量造成影響。土壤是最重要的農(nóng)業(yè)生產(chǎn)資料, 它不僅提供了作物吸收養(yǎng)分的環(huán)境條件, 同時(shí)也是產(chǎn)量的制約因素之一[8]。土壤的養(yǎng)分豐缺程度和分布情況直接影響到作物生長及營養(yǎng)特征, 了解土壤肥力狀況及其不同區(qū)域的分布差異, 對充分利用土地資源, 制定合理施肥措施, 實(shí)現(xiàn)優(yōu)質(zhì)生產(chǎn)和構(gòu)建生態(tài)經(jīng)濟(jì)型農(nóng)業(yè)有重要意義[9]。目前, 關(guān)于花椒的研究主要集中在花椒營養(yǎng)需求、品種選育、苗期管理、化感作用等方面, 但關(guān)于地形因子、土壤肥力與花椒產(chǎn)量的關(guān)系研究尚鮮見相關(guān)報(bào)導(dǎo)。本研究對九龍坡花椒種植區(qū)地形因子差異、土壤養(yǎng)分分布特征及其與花椒產(chǎn)量的關(guān)系進(jìn)行了研究分析, 確定土壤肥力是產(chǎn)量的主要限制因子, 旨在掌握九龍坡花椒種植區(qū)土壤養(yǎng)分豐缺狀況, 為科學(xué)、合理地制定花椒高效施肥措施提供理論依據(jù)。

    1 研究區(qū)域概況與研究方法

    1.1 研究區(qū)域概況

    本研究區(qū)選擇在重慶市九龍坡花椒生產(chǎn)區(qū)西彭鎮(zhèn)和陶家鎮(zhèn)(29°20′N, 106°19′E), 屬于中亞熱帶濕潤季風(fēng)氣候, 年平均氣溫19.2 ℃, 年平均降水量1 204 mm, 無霜期347 d。土壤類型以沙溪廟組發(fā)育的灰棕紫泥土和遂寧組發(fā)育的紅棕紫泥土為主。地貌類型以丘陵、山地為主。供試花椒品種為當(dāng)?shù)刂鲗?dǎo)品種‘九葉青’。

    1.2 樣品采集與測定方法

    通過大量實(shí)地調(diào)查, 根據(jù)研究區(qū)域地貌特點(diǎn)、土壤類型以及農(nóng)戶花椒種植情況, 結(jié)合農(nóng)戶多年種植經(jīng)驗(yàn), 將研究區(qū)劃分為高產(chǎn)區(qū)(產(chǎn)量>25 kg?株-1)、中產(chǎn)區(qū)(產(chǎn)量為15~25 kg?株-1)、低產(chǎn)區(qū)(產(chǎn)量<15 kg?株-1)。2017年3月, 經(jīng)過實(shí)地走訪調(diào)查, 隨機(jī)選擇花椒種植連片、種植規(guī)模相當(dāng)、花椒樹齡、種植模式、田間管理基本一致的30戶有代表性的農(nóng)戶作為采樣點(diǎn)。用GPS記錄該點(diǎn)的海拔, 同時(shí)測定坡度。采集的土壤為春季花期尚未施肥的基礎(chǔ)土壤, 選擇長勢基本一致的5棵花椒樹, 環(huán)繞每棵樹采集距離樹根15~25 cm處0~20 cm土層的土壤(不破壞花椒主根), 同時(shí)作采樣標(biāo)志。將采集的土壤樣品混合, 在自然條件下風(fēng)干、研磨, 過0.25 mm篩備用, 測定土壤理化性質(zhì)。在選擇這30個采樣點(diǎn)的同時(shí), 詳細(xì)調(diào)查每個農(nóng)戶近幾年花椒產(chǎn)量情況, 結(jié)合本次產(chǎn)量劃分出高、中、低產(chǎn)量水平各10戶。2017年8月, 選擇有采樣標(biāo)記的花椒樹, 收獲花椒, 除去莖、葉, 并測定花椒鮮重, 計(jì)算每棵花椒樹的平均產(chǎn)量。

    采集的土壤經(jīng)風(fēng)干、研磨、過篩后, 進(jìn)行養(yǎng)分含量測定。測定方法pH采用玻璃電極法[10], 有機(jī)質(zhì)采用重鉻酸鉀容量法[10], 堿解氮采用堿解擴(kuò)散法[10], 有效磷采用NaHCO3浸提-鉬銻抗比色法[10], 有效鉀采用乙酸銨浸提-火焰光度法[10], 陽離子交換量采用乙酸銨交換-原子吸收分光光度法[10], 土壤交換性酸采用氯化鋇-三乙醇胺法[10], 土壤有效鈣、有效鎂采用KCl浸提-原子吸收分光光度法[10], 有效鐵、有效錳、有效銅、有效鋅采用DTPA浸提-原子吸收分光光度法[10], 土壤水溶性硼采用姜黃素比色法[10]。

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

    所有數(shù)據(jù)采用Microsoft Excel 2010和SPSS 20.0統(tǒng)計(jì)軟件進(jìn)行分析, 采用描述性統(tǒng)計(jì)分析, ANOVA方差分析, LSD法進(jìn)行多重比較, Pearson相關(guān)分析法、通徑分析法和逐步回歸分析法。

    2 結(jié)果與分析

    2.1 不同產(chǎn)區(qū)的花椒產(chǎn)量水平狀況

    由表1可以看出, 不同產(chǎn)區(qū)的花椒產(chǎn)量差異達(dá)顯著水平(<0.05); 不同產(chǎn)區(qū)產(chǎn)量變幅為中產(chǎn)區(qū)≈高產(chǎn)區(qū)>低產(chǎn)區(qū)。從變異系數(shù)來看, 3個產(chǎn)區(qū)均屬中等變異, 變異系數(shù)最大的為中產(chǎn)區(qū), 其次為低產(chǎn)區(qū), 高產(chǎn)區(qū)最小, 表明花椒高產(chǎn)區(qū)產(chǎn)量整齊度優(yōu)于中低產(chǎn)區(qū)。

    表1 不同產(chǎn)區(qū)花椒產(chǎn)量水平

    產(chǎn)量表示為平均值±標(biāo)準(zhǔn)差, 同列不同小寫字母表示不同產(chǎn)區(qū)間差異顯著(<0.05)。Yield was expressed as mean ± S.D., and different lowercase letters in the same column indicate significant differences at 0.05 level among different yield areas.

    2.2 花椒不同產(chǎn)區(qū)地形因子狀況

    由圖1可以看出, 花椒不同產(chǎn)區(qū)海拔高度的變化幅度為低產(chǎn)區(qū)>中產(chǎn)區(qū)>高產(chǎn)區(qū), 低產(chǎn)區(qū)集中分布在海拔240~300 m, 中產(chǎn)區(qū)集中分布在海拔230~340 m, 高產(chǎn)區(qū)集中分布在海拔290~320 m。不同產(chǎn)區(qū)之間的海拔差異不顯著, 從3個產(chǎn)區(qū)海拔高度來看, 九龍坡花椒生產(chǎn)區(qū)的海拔普遍在200~500 m, 而高產(chǎn)區(qū)主要集中于300 m左右的海拔高度?;ń凡煌a(chǎn)區(qū)地面坡度差異不大, 在0°~15°均有分布, 且集中分布在3°~14°, 其中低產(chǎn)區(qū)平均坡度為7.7°, 中產(chǎn)區(qū)平均坡度為8.6°, 高產(chǎn)區(qū)平均坡度為10.4°(如箱圖中黑點(diǎn)所示), 隨著產(chǎn)量提高, 坡度略有增加, 但未達(dá)顯著差異。

    2.3 花椒不同產(chǎn)區(qū)土壤肥力特征

    2.3.1 pH

    由表2可見, 九龍坡花椒種植區(qū)土壤為酸性土, pH均在6.5以下。其中pH<5.5的強(qiáng)酸性土占66.7%, pH<4.5的極強(qiáng)酸性土占23.3%。從平均pH來看, 全產(chǎn)區(qū)的平均pH為5.1, 各產(chǎn)區(qū)平均pH表現(xiàn)為高產(chǎn)區(qū)>中產(chǎn)區(qū)>低產(chǎn)區(qū), 高產(chǎn)區(qū)與低產(chǎn)區(qū)之間差異顯著, 中產(chǎn)區(qū)與低、高產(chǎn)區(qū)的pH差異不顯著。3個產(chǎn)區(qū)pH變化幅度為低產(chǎn)區(qū)>中產(chǎn)區(qū)>高產(chǎn)區(qū)。從變異系數(shù)來看, 高產(chǎn)區(qū)變異系數(shù)最大, 低產(chǎn)區(qū)次之, 中產(chǎn)區(qū)最小, 且均屬中等變異。

    2.3.2 有機(jī)質(zhì)

    土壤有機(jī)質(zhì)在維持土壤肥力、土壤環(huán)境質(zhì)量、作物產(chǎn)量上有著重要的作用[11-12]。根據(jù)全國第2次土壤普查養(yǎng)分分級標(biāo)準(zhǔn)(下同), 九龍坡花椒種植區(qū)土壤有機(jī)質(zhì)含量處于中低等水平(表2), 其中<10 g?kg-1的土壤占10%, 10~20 g?kg-1的土壤占43.3%, 20~30 g?kg-1的土壤占36.7%, >30 g?kg-1的土壤占10%。全產(chǎn)區(qū)有機(jī)質(zhì)平均含量為19.1 g?kg-1, 各產(chǎn)區(qū)平均含量表現(xiàn)為中產(chǎn)區(qū)>高產(chǎn)區(qū)>低產(chǎn)區(qū), 3個產(chǎn)區(qū)之間差異不顯著, 低產(chǎn)區(qū)有機(jī)質(zhì)平均含量較低, 僅為16.6 g?kg-1,表明花椒低產(chǎn)地區(qū)應(yīng)適當(dāng)增加有機(jī)肥的施用。不同產(chǎn)區(qū)有機(jī)質(zhì)含量變化幅度為中產(chǎn)區(qū)>高產(chǎn)區(qū)>低產(chǎn)區(qū); 從變異系數(shù)來看, 低產(chǎn)區(qū)變異系數(shù)最大, 高產(chǎn)區(qū)次之, 中產(chǎn)區(qū)最小, 且均屬中等變異。

    圖1 不同花椒產(chǎn)區(qū)海拔與坡度因子分布狀況

    不同字母表示不同產(chǎn)區(qū)間差異顯著(<0.05)。Different lowercase letters indicate significant differences at 0.05 level among different yield areas.

    表2 花椒不同產(chǎn)區(qū)土壤肥力特征

    同行不同小寫字母表示不同產(chǎn)區(qū)間差異顯著(<0.05)。Different lowercase letters in the same line indicate significant differences at 0.05 level among different yield areas.

    2.3.3 大量元素

    九龍坡花椒產(chǎn)區(qū)大量元素含量水平變幅較大(表2)。其中堿解氮缺乏的土壤(<90 mg?kg-1)占36.7%, 堿解氮適中的土壤(90~120 mg?kg-1)占23.3%, 堿解氮豐富的土壤(>120 mg?kg-1)占40.0%。全產(chǎn)區(qū)堿解氮平均含量為114.9 mg?kg-1。從不同產(chǎn)區(qū)看, 堿解氮平均含量為低產(chǎn)區(qū)>中產(chǎn)區(qū)>高產(chǎn)區(qū), 但各產(chǎn)區(qū)之間堿解氮含量未達(dá)顯著差異。不同產(chǎn)區(qū)堿解氮變化幅度為高產(chǎn)區(qū)>中產(chǎn)區(qū)>低產(chǎn)區(qū), 變異系數(shù)為低產(chǎn)區(qū)>中產(chǎn)區(qū)>高產(chǎn)區(qū), 且均屬于中等變異。由此可見, 九龍坡花椒種植區(qū)土壤堿解氮缺乏的地區(qū)應(yīng)注意氮肥的施用, 而堿解氮含量豐富的區(qū)域應(yīng)該注意氮肥減量施用。

    九龍坡花椒種植區(qū)有效磷平均含量屬于豐富水平, 但變幅很大, 變異系數(shù)屬于強(qiáng)變異。有效磷極度缺乏的區(qū)域(<5 mg?kg-1)占10.0%, 有效磷缺乏的區(qū)域(5~10 mg?kg-1)占16.7%, 有效磷適中的區(qū)域(10~20 mg?kg-1)占20.0%, 有效磷含量豐富的區(qū)域(>20 mg?kg-1)占53.3%。從不同產(chǎn)區(qū)看, 有效磷平均含量為高產(chǎn)區(qū)最大, 中產(chǎn)區(qū)次之, 低產(chǎn)區(qū)最小, 但各產(chǎn)區(qū)之間差異不顯著。3個產(chǎn)區(qū)有效磷含量變化幅度為中產(chǎn)區(qū)>高產(chǎn)區(qū)>低產(chǎn)區(qū), 且變異系數(shù)均屬于強(qiáng)變異。

    花椒種植區(qū)速效鉀平均含量屬于適中水平, 但變幅較大, 其中約20%的土壤速效鉀含量屬于極缺乏水平(<50 mg?kg-1), 36.7%的土壤速效鉀屬缺乏水平, 23.3%的土壤速效鉀含量適中, 土壤速效鉀豐富的土壤約占20%。從不同產(chǎn)區(qū)來看, 速效鉀平均含量為中產(chǎn)區(qū)>高產(chǎn)區(qū)>低產(chǎn)區(qū), 中產(chǎn)區(qū)與低產(chǎn)區(qū)之間差異顯著, 而高產(chǎn)區(qū)與中、低產(chǎn)區(qū)之間差異均不顯著。不同產(chǎn)區(qū)變幅為高產(chǎn)區(qū)>中產(chǎn)區(qū)>低產(chǎn)區(qū), 從變異系數(shù)看, 3個產(chǎn)區(qū)均屬強(qiáng)變異。可見, 在速效鉀平均含量缺乏的低產(chǎn)區(qū), 應(yīng)注重鉀肥的合理施用。

    九龍坡花椒種植區(qū)土壤有效鈣平均含量屬于豐富水平, 其中, 有效鈣含量豐富的土壤(>1 200 mg?kg-1)約占90%, 變異系數(shù)屬強(qiáng)變異。從不同產(chǎn)區(qū)看, 土壤有效鈣平均含量為高產(chǎn)區(qū)>中產(chǎn)區(qū)>低產(chǎn)區(qū), 高產(chǎn)區(qū)與中、低產(chǎn)區(qū)之間差異顯著, 中、低產(chǎn)區(qū)之間差異不顯著。3個產(chǎn)區(qū)的平均有效鈣含量均屬豐富水平, 變幅為高產(chǎn)區(qū)>低產(chǎn)區(qū)>中產(chǎn)區(qū); 變異系數(shù)高產(chǎn)區(qū)>低產(chǎn)區(qū)>中產(chǎn)區(qū), 且均屬于中等變異。由此可見, 九龍坡花椒種植區(qū)土壤有效鈣含量豐富, 可以滿足作物對鈣的需求。

    花椒種植區(qū)土壤有效鎂平均含量屬豐富水平, 有效鎂含量豐富的土壤(>180 mg?kg-1)占90%, 變異系數(shù)為中等變異。從不同產(chǎn)區(qū)看, 從高產(chǎn)區(qū)到低產(chǎn)區(qū)土壤有效鎂平均含量呈逐步降低趨勢, 且高產(chǎn)區(qū)與低產(chǎn)區(qū)之間差異顯著, 而中產(chǎn)區(qū)與其他兩個產(chǎn)區(qū)差異均不顯著。不同產(chǎn)區(qū)的有效鎂含量變幅為中產(chǎn)區(qū)>低產(chǎn)區(qū)>高產(chǎn)區(qū); 3個產(chǎn)區(qū)的變異系數(shù)均屬中等變異。可見, 花椒種植區(qū)土壤有效鎂含量豐富, 可滿足作物對鎂的需求。

    2.3.4 土壤交換性能

    土壤陽離子交換量(CEC)是土壤的基本特性和重要肥力因素之一, 它反映了土壤保蓄、供應(yīng)和緩沖陽離子養(yǎng)分的能力[13]。九龍坡花椒種植區(qū)土壤平均CEC處于高水平范圍, CEC>20 cmol(+)?kg-1的土壤占80%, 變異系數(shù)為中等變異(表2)。從不同產(chǎn)區(qū)來看, 從高產(chǎn)區(qū)到低產(chǎn)區(qū)土壤的CEC呈逐漸降低趨勢, 其中高產(chǎn)區(qū)與中、低產(chǎn)區(qū)差異顯著, 中、低產(chǎn)區(qū)之間差異不顯著。土壤CEC變幅高產(chǎn)區(qū)>低產(chǎn)區(qū)>中產(chǎn)區(qū); 從變異系數(shù)看, 高產(chǎn)區(qū)最大, 其次為低產(chǎn)區(qū), 中產(chǎn)區(qū)最小。可見, 花椒種植區(qū)CEC處于高水平, 土壤的保肥能力較好, 但中低產(chǎn)區(qū)中有一部分土壤CEC較低, 這一部分土壤應(yīng)該注重優(yōu)化施肥, 防止土壤酸化現(xiàn)象加重。

    土壤交換性酸主要由土壤膠體吸附的H+和Al3+構(gòu)成[14], 花椒種植區(qū)平均土壤交換性酸總量為8.1 cmol(+)?kg-1, 從全產(chǎn)區(qū)來看, <5 cmol(+)?kg-1的土壤占16.7%, 5~10 cmol(+)?kg-1的土壤占56.6%, >10 cmol(+)?kg-1的土壤占26.7%; 交換性酸變化幅度為1.52~13.71, 變異系數(shù)為中等變異。從不同產(chǎn)區(qū)看, 交換性酸總量低產(chǎn)區(qū)>中產(chǎn)區(qū)>高產(chǎn)區(qū), 3個產(chǎn)區(qū)之間差異均不顯著; 變幅最大的是低產(chǎn)區(qū), 中產(chǎn)區(qū)次之, 最小的是高產(chǎn)區(qū)。變異系數(shù)低產(chǎn)區(qū)>高產(chǎn)區(qū)>中產(chǎn)區(qū), 且均為中等變異。由此可見, 花椒種植區(qū)土壤交換性酸總量為適中水平, 土壤交換性酸與pH之間存在極顯著負(fù)相關(guān)關(guān)系, 在pH較低的強(qiáng)酸性土壤, 交換性酸由交換鋁主導(dǎo), Al3+對作物根系有毒害作用, 因此, 在這一部分區(qū)域應(yīng)該注重土壤的酸性改良, 提高土壤pH。

    2.3.5 微量元素

    九龍坡花椒種植區(qū)有效鐵平均含量豐富, 變幅很大, 屬強(qiáng)變異(表2)。從有效鐵分布頻率來看, >10 mg?kg-1的土壤占46.7%, 4.5~10 mg?kg-1的土壤占16.6%, <4.5 mg?kg-1的土壤占36.6%。土壤有效鐵平均含量以低產(chǎn)區(qū)最高, 中產(chǎn)區(qū)次之, 高產(chǎn)區(qū)最低, 其中高產(chǎn)區(qū)與中、低產(chǎn)區(qū)之間差異均顯著, 中、低產(chǎn)區(qū)之間差異不顯著。3個產(chǎn)區(qū)變幅順序?yàn)榈彤a(chǎn)區(qū)>中產(chǎn)區(qū)>高產(chǎn)區(qū), 且變異系數(shù)均屬于強(qiáng)變異??梢? 中低產(chǎn)區(qū)的有效鐵含量極其豐富。

    花椒種植區(qū)有效錳平均含量豐富, 變異系數(shù)屬強(qiáng)變異, 從分布頻率看, >30 mg?kg-1的土壤占83.3%, 15~30 mg?kg-1的土壤占10.0%, <15 mg?kg-1的土壤占6.7%。土壤有效錳平均含量以高產(chǎn)區(qū)最高, 中產(chǎn)區(qū)次之, 低產(chǎn)區(qū)最低, 其中高、低產(chǎn)區(qū)之間差異顯著, 中產(chǎn)區(qū)與高、低產(chǎn)區(qū)差異均不顯著。3個產(chǎn)區(qū)有效錳含量變異幅度中產(chǎn)區(qū)>高產(chǎn)區(qū)>低產(chǎn)區(qū), 中產(chǎn)區(qū)變異系數(shù)屬中等變異, 高、低產(chǎn)區(qū)變異系數(shù)均屬強(qiáng)變異。

    九龍坡花椒種植區(qū)有效銅平均含量屬于較豐富水平, 變幅較大, 變異系數(shù)屬強(qiáng)變異。從分布頻率來看, <1.0 mg?kg-1的土壤占43.3%, 1.0~2.0 mg?kg-1的土壤占33.4%, >2.0 mg?kg-1的土壤占23.3%。3個產(chǎn)區(qū)有效銅平均含量以低產(chǎn)區(qū)最高, 中產(chǎn)區(qū)次之, 高產(chǎn)區(qū)最低, 3個產(chǎn)區(qū)之間差異均不顯著; 其變化幅度低產(chǎn)區(qū)>中產(chǎn)區(qū)>高產(chǎn)區(qū), 從變異系數(shù)看, 均屬于強(qiáng)變異。

    九龍坡花椒種植區(qū)有效鋅平均含量為4.5 mg?kg-1, 屬于含量豐富范圍, 變異系數(shù)屬強(qiáng)變異。從有效鋅分布頻率來看, >3.0 mg?kg-1的土壤占全區(qū)的66.7%, 1.0~3.0 mg?kg-1的土壤占26.6%, <1.0 mg?kg-1的土壤占6.7%。土壤有效鋅平均含量最高的是中產(chǎn)區(qū), 高產(chǎn)區(qū)次之, 低產(chǎn)區(qū)最低, 3個產(chǎn)區(qū)之間差異不顯著。土壤有效鋅變幅中產(chǎn)區(qū)>低產(chǎn)區(qū)>高產(chǎn)區(qū); 中產(chǎn)區(qū)的變異系數(shù)屬中等變異, 低、高產(chǎn)區(qū)為強(qiáng)變異。

    花椒種植區(qū)水溶性硼平均含量屬于缺乏水平, 變異系數(shù)為中等變異。從分布頻率來看, <0.3 mg?kg-1的土壤占56.7%, 0.3~0.5 mg?kg-1的土壤占40.0%, 0.5~1.0 mg?kg-1的土壤占3.3%。土壤水溶性硼平均含量從低產(chǎn)區(qū)到高產(chǎn)區(qū)呈逐漸升高趨勢, 高產(chǎn)區(qū)與低產(chǎn)區(qū)之間差異顯著, 而中產(chǎn)區(qū)與其他兩個產(chǎn)區(qū)均差異不顯著。土壤水溶性硼變幅高產(chǎn)區(qū)>中產(chǎn)區(qū)>低產(chǎn)區(qū), 3個產(chǎn)區(qū)的變異系數(shù)均屬強(qiáng)變異??梢? 九龍坡區(qū)土壤水溶性硼含量缺乏, 施肥時(shí)應(yīng)注意硼肥的施用。

    2.4 地形因子和土壤肥力對花椒產(chǎn)量的影響

    2.4.1 地形因子、土壤肥力因子間及其與產(chǎn)量的相關(guān)分析

    由表3可知, 花椒產(chǎn)量與土壤有效鈣和陽離子交換量(CEC)呈極顯著正相關(guān), 與pH、有效錳和水溶性硼呈顯著正相關(guān)。海拔、坡度因子與花椒產(chǎn)量及其他土壤肥力因子的相關(guān)性均不顯著。從土壤肥力因子間相互關(guān)系來看, 土壤pH與土壤有效鈣、有效鎂、有效錳呈顯著正相關(guān), 而與交換性酸、有效鐵、有效銅呈顯著負(fù)相關(guān); 土壤有機(jī)質(zhì)與堿解氮、速效鉀、有效鋅呈顯著正相關(guān); 土壤堿解氮與有效鋅呈顯著負(fù)相關(guān); 土壤有效磷與速效鉀、交換性酸呈顯著正相關(guān); 土壤速效鉀與有效鋅呈顯著正相關(guān), 而與有效鐵、有效銅呈顯著負(fù)相關(guān); 土壤有效鈣與有效鎂、陽離子交換量、有效錳呈顯著正相關(guān), 而與有效鐵、有效銅呈顯著負(fù)相關(guān); 土壤有效鎂與陽離子交換量、有效錳呈顯著正相關(guān); 土壤陽離子交換量與有效錳呈顯著正相關(guān), 而與有效鐵呈顯著負(fù)相關(guān); 土壤交換性酸與有效鋅呈顯著正相關(guān); 土壤有效鐵與有效銅呈顯著正相關(guān)??梢? 地形因子與產(chǎn)量及其他因子相關(guān)性不大, 可以先行排除; 而各土壤肥力因子之間存在著錯綜復(fù)雜的關(guān)系, 各因子相互影響、相互制約, 且通過彼此間的相互作用間接對花椒產(chǎn)量造成影響。因此, 不能簡單地看產(chǎn)量與各因子的相關(guān)性, 忽視各因子間的相互作用。

    2.4.2 土壤肥力與花椒產(chǎn)量的通徑分析

    為進(jìn)一步了解土壤肥力因子對花椒產(chǎn)量的影響, 以花椒產(chǎn)量()為因變量, 各土壤肥力因子(X)為自變量, 利用通徑分析[15], 將各因子(X)與產(chǎn)量()間的相關(guān)系數(shù)(r)分解為直接作用和間接作用, 計(jì)算每個因子對產(chǎn)量()的直接通徑系數(shù)(b)和間接通徑系數(shù)(rb), 評估每個因子對花椒產(chǎn)量的作用及因子間的相互影響(表4)。從直接通徑系數(shù)看, 各土壤肥力因子對花椒產(chǎn)量的直接作用有效鈣>CEC>交換性酸>有效銅>有效鐵>有效鋅>速效鉀>有效磷>有效鎂>有效錳>pH>水溶性硼>有機(jī)質(zhì)>堿解氮。表明在九龍坡花椒種植區(qū)各土壤肥力因子中, 有效鈣、CEC、交換性酸、有效銅、有效鐵、有效鋅是影響產(chǎn)量的主要因子。其中有效鈣、交換性酸、有效鐵、有效鋅對花椒產(chǎn)量的直接作用為正, CEC、有效銅對花椒產(chǎn)量直接作用為負(fù)。有效鈣的直接作用最大, 對產(chǎn)量的影響也最顯著。但從間接通徑系數(shù)來看, 有效鈣通過其他因子間接作用于花椒產(chǎn)量時(shí), 其間接通徑系數(shù)大多為負(fù), 表明在各因子的相互作用中, 有效鈣對產(chǎn)量的效應(yīng)有一部分被抑制了, 其中CEC、交換性酸、有效鐵、有效錳對有效鈣抑制作用最強(qiáng)烈, 這可能是因?yàn)樵撗芯繀^(qū)土壤有效鈣含量過于豐富, 從而導(dǎo)致部分有效鈣的效果沒有發(fā)揮作用, 造成鈣素的浪費(fèi); 同時(shí), 有效鈣的直接通徑系數(shù)大于其間接通徑系數(shù)之和, 表明有效鈣的作用是通過直接作用來體現(xiàn)。另一方面, pH、有效磷、速效鉀、有效鎂、CEC、有效錳、水溶性硼通過有效鈣間接作用于花椒產(chǎn)量時(shí), 使得這些因子的效應(yīng)得到增強(qiáng)。CEC對產(chǎn)量的直接通徑系數(shù)為負(fù), 但相關(guān)分析表明, CEC與產(chǎn)量之間有顯著的正相關(guān)關(guān)系, 原因是當(dāng)CEC通過pH、有效鈣、交換性酸、有效銅、有效鋅、水溶性硼間接作用于產(chǎn)量時(shí), CEC對產(chǎn)量的作用得到增強(qiáng), 且間接通徑系數(shù)之和大于直接通徑系數(shù), 說明CEC對產(chǎn)量的影響主要由間接作用體現(xiàn)。由此可見, 各土壤肥力因子之間的相互關(guān)系復(fù)雜, 施肥時(shí)應(yīng)注重優(yōu)化施肥, 改善土壤理化性狀, 調(diào)節(jié)各元素之間的平衡。

    表3 花椒產(chǎn)量與地形因子、土壤肥力因子的相關(guān)矩陣(r值)

    Y: 產(chǎn)量; Ele: 海拔; Slo: 坡度; OM: 有機(jī)質(zhì); AN: 堿解氮; AP: 有效磷; AK: 速效鉀; ACa: 有效鈣; AMg: 有效鎂; CEC: 陽離子交換量; EA: 交換性酸; AFe: 有效鐵; AMn: 有效錳; ACu: 有效銅; AZn: 有效鋅; WB: 水溶性硼。*表示在0.05水平上顯著相關(guān)(<0.05); **表示在0.01水平上顯著相關(guān)(<0.01)。Y: yield; Ele: elevation; Slo: slope; OM: organic matter; AN: alkali-hydrolysis nitrogen; AP: available phosphorus; AK: available potassium; ACa: available calcium; AMg: available magnesium; CEC: cation exchange capacity; EA: exchangeable acid; AFe: effective iron; AMn: effective manganese; ACu: effective copper; AZn: effective zinc; WB: water-soluble boron. * indicates a significant correlation at 0.05 level; ** indicates a significant correlation at 0.01 level.

    表4 土壤肥力因子對花椒產(chǎn)量(Y)的通徑系數(shù)

    1: pH;2: 有機(jī)質(zhì);3: 堿解氮;4: 有效磷;5: 速效鉀;6: 有效鈣;7: 有效鎂;8: 陽離子交換量;9: 交換性酸;10: 有效鐵;11: 有效錳;12: 有效銅;13: 有效鋅;14: 水溶性硼。1: pH;2: organic matter;3: alkali-hydrolysis nitrogen;4: available phosphorus;5: available potassium;6: available calcium;7: available magnesium;8: cation exchange capacity;9: exchangeable acid;10: available iron;11: available manganese;12: available copper;13: available zinc;14: water soluble boron.

    2.4.3 土壤肥力因子與花椒產(chǎn)量的逐步回歸分析

    為進(jìn)一步明確土壤肥力因子與花椒產(chǎn)量的數(shù)量關(guān)系, 在相關(guān)分析和通徑分析的基礎(chǔ)上, 利用逐步回歸分析方法, 剔除不顯著變量, 建立了有效鈣(6)與花椒產(chǎn)量()的最優(yōu)回歸方程:=11.693+0.0036(表5、表6); 經(jīng)檢驗(yàn),2=0.344,=14.712,<0.01, 回歸方程達(dá)極顯著水平; 回歸方程系數(shù)經(jīng)檢驗(yàn)也達(dá)到顯著水平, 說明回歸方程具有統(tǒng)計(jì)學(xué)意義。表明在九龍坡花椒種植區(qū), 土壤有效鈣含量與花椒產(chǎn)量有顯著正效應(yīng), 可以根據(jù)土壤有效鈣含量預(yù)測花椒產(chǎn)量水平。

    表5 花椒產(chǎn)量與土壤有效鈣的逐步回歸分析模型匯總

    表6 花椒產(chǎn)量與土壤有效鈣的逐步回歸方程系數(shù)及檢驗(yàn)

    a: 預(yù)測變量(常量);6: 有效鈣。a: predictor (constant);6: available calcium.

    3 討論

    地形的差異往往導(dǎo)致土壤養(yǎng)分的空間變異[16-17], 獨(dú)特地形條件產(chǎn)生的“小氣候”對作物生長也有影響[18]。本研究中, 花椒種植區(qū)普遍處于200~500 m海拔范圍內(nèi), 且高產(chǎn)區(qū)的海拔集中分布于300 m左右。相關(guān)分析表明, 海拔因子與土壤肥力、花椒產(chǎn)量的相關(guān)性不顯著, 這可能是九龍坡花椒種植區(qū)地貌類型為丘陵山地, 海拔差異不大。相關(guān)分析表明, 坡度因子與花椒產(chǎn)量關(guān)系更密切, 從低產(chǎn)區(qū)到高產(chǎn)區(qū), 坡度逐漸增大, 但未達(dá)顯著差異。地形因子與土壤肥力因子的相關(guān)關(guān)系不顯著, 這可能是該研究區(qū)范圍較小, 且地形地貌類型差異不大的原因。

    土壤養(yǎng)分主要由土壤pH、有機(jī)質(zhì)、有效養(yǎng)分和交換性能構(gòu)成, 與作物產(chǎn)量具有顯著的相關(guān)性[19-23]。本研究中, 花椒種植區(qū)土壤pH均<6.5, 從低產(chǎn)區(qū)到高產(chǎn)區(qū), pH逐漸增大; 相關(guān)分析表明, 土壤pH與花椒產(chǎn)量有顯著相關(guān)關(guān)系。土壤pH能影響作物營養(yǎng)成分的吸收利用, 適宜的pH能減少作物缺素癥和不良元素的毒害作用[24-25]。通徑分析表明, 適宜的土壤pH能增強(qiáng)有效鈣對花椒產(chǎn)量的促進(jìn)作用, 同時(shí)抑制交換性酸、有效銅的不良作用。花椒更合適在pH較高的土壤上生長, 在pH較低的花椒低產(chǎn)區(qū), 應(yīng)注重有機(jī)肥和堿性肥料的施用。本研究中通徑分析表明, 有機(jī)質(zhì)對產(chǎn)量的間接通徑系數(shù)之和大于直接通徑系數(shù), 表明有機(jī)質(zhì)主要通過間接作用影響產(chǎn)量; 但九龍坡花椒種植區(qū)土壤有機(jī)質(zhì)含量屬中低量水平, 特別是低產(chǎn)區(qū), 土壤有機(jī)質(zhì)水平低下, 因此應(yīng)注意有機(jī)肥的施用, 提高土壤理化性狀。

    土壤有效養(yǎng)分含量是土壤對作物養(yǎng)分供應(yīng)能力的體現(xiàn), 是評價(jià)土壤肥力的重要標(biāo)準(zhǔn)[26]。本研究中, 九龍坡花椒不同產(chǎn)區(qū)養(yǎng)分豐缺情況不同, 但總體來說高產(chǎn)區(qū)肥力水平最高, 中產(chǎn)區(qū)次之, 低產(chǎn)區(qū)最差。相關(guān)分析表明, 土壤有效鈣、有效錳、水溶性硼對花椒產(chǎn)量有極顯著的正效應(yīng)。鈣能增強(qiáng)植物運(yùn)輸能力和光合效率, 且能提高細(xì)胞壁的穩(wěn)定性[27-28]。通徑分析表明, 有效鈣對花椒產(chǎn)量的直接作用最大, 但豐富的鈣離子富集于花椒根系, 可能會使其他離子難以被根系吸收。錳是植物必需的微量營養(yǎng)元素, 對光合作用、維持細(xì)胞器正常結(jié)構(gòu)有重要作用[29]。九龍坡花椒種植區(qū)錳含量豐富, 能滿足作物對錳的需求, 從低產(chǎn)區(qū)到高產(chǎn)區(qū)錳含量逐漸升高。硼是作物必需的微量營養(yǎng)元素之一, 硼在穩(wěn)定細(xì)胞壁結(jié)構(gòu)、維持細(xì)胞生命活動中有重要意義[30]?;ń贩N植區(qū)土壤水溶性硼屬缺乏水平, 硼對花椒產(chǎn)量有顯著的正效應(yīng), 九龍坡花椒種植區(qū)硼含量缺乏, 施肥時(shí)應(yīng)該注意硼肥的施用。

    九龍坡花椒CEC屬于高水平范圍, 從低產(chǎn)區(qū)到高產(chǎn)區(qū)CEC逐漸增加。相關(guān)分析表明, CEC不僅與花椒產(chǎn)量存在顯著正相關(guān)性, 且與有效鈣、有效鎂、有效鐵、有效錳也存在顯著相關(guān)性, 說明CEC在土壤中與多個土壤養(yǎng)分因子有相互作用, 關(guān)系錯綜復(fù)雜。通徑分析表明, CEC對產(chǎn)量的直接作用為負(fù), 但CEC的直接通徑系數(shù)小于間接通徑系數(shù)之和, 因此, CEC對產(chǎn)量的作用主要由間接作用體現(xiàn)。交換性酸含量高的土壤中鋁、錳和重金屬的活性較強(qiáng), 容易導(dǎo)致作物缺素和元素毒害, 影響作物正常生長[31-33]。本研究中, 從低產(chǎn)區(qū)到高產(chǎn)區(qū)土壤交換性酸含量逐漸下降; 交換性酸與pH存在極顯著負(fù)相關(guān)關(guān)系, 酸度越高的土壤, 交換性酸含量越高; 且花椒產(chǎn)量與交換性酸也存在負(fù)相關(guān)效應(yīng), 通徑分析進(jìn)一步表明, 交換性酸通過其他因子對花椒產(chǎn)量產(chǎn)生作用時(shí), 其間接通徑系數(shù)為負(fù), 且間接通徑系數(shù)之和大于直接通徑系數(shù), 表明交換性酸以間接作用體現(xiàn), 且對花椒產(chǎn)量的影響為負(fù)效應(yīng)。說明在土壤酸化嚴(yán)重的地區(qū), 應(yīng)該注重施用堿性肥料, 如石灰、硅鈣鉀鎂肥等, 調(diào)節(jié)土壤pH, 防止土壤酸化程度進(jìn)一步加深, 同時(shí)配施有機(jī)肥, 調(diào)理土壤理化性狀, 改良土壤質(zhì)量。

    4 結(jié)論

    九龍坡花椒種植區(qū)高產(chǎn)區(qū)海拔集中在300 m左右, 坡度隨產(chǎn)量提高略有增加, 但未達(dá)顯著水平; 總的來說, 地形因子對花椒產(chǎn)量影響不顯著。九龍坡花椒種植區(qū)土壤均屬酸性土, 土壤肥力屬中高水平范圍, 土壤變異幅度較大, 高產(chǎn)區(qū)土壤肥力最高, 中產(chǎn)區(qū)次之, 低產(chǎn)區(qū)最差?;ń樊a(chǎn)量與土壤pH、有效鈣、CEC、有效錳、水溶性硼有顯著的正相關(guān)效應(yīng), 其中與有效鈣的顯著性最強(qiáng)。土壤肥力因子之間的相互作用復(fù)雜, 通徑分析表明, 有效鈣的直接作用最大, 但豐富的有效鈣含量會減弱其他土壤肥力因子對花椒產(chǎn)量的作用, 逐步回歸分析構(gòu)建了有效鈣與花椒產(chǎn)量的最優(yōu)回歸方程:=11.693+0.0036。綜上所述, 九龍坡花椒種植區(qū)施肥應(yīng)該注重土壤養(yǎng)分之間的平衡, 低產(chǎn)區(qū)應(yīng)該注重土壤酸化的改良, 合理使用堿性肥料(如石灰、硅鈣鉀鎂肥等), 配施有機(jī)肥, 調(diào)理土壤理化性狀。

    [1] 侯麗秀, 魏安智, 王麗華, 等. 花椒轉(zhuǎn)錄組SSR信息分析及其分子標(biāo)記開發(fā)[J]. 農(nóng)業(yè)生物技術(shù)學(xué)報(bào), 2018, 26(7): 1226–1236 HOU L X, WEI A Z, WANG L H, et al. Analysis of SSR loci and development of molecular markers intranscriptome[J]. Journal of Agricultural Biotechnology, 2018, 26(7): 1226–1236

    [2] 王景燕, 龔偉, 包秀蘭, 等. 水肥耦合對漢源花椒幼苗葉片光合作用的影響[J]. 生態(tài)學(xué)報(bào), 2016, 36(5): 1321–1330 WANG J Y, GONG W, BAO X L, et al. Coupling effects of water and fertilizer on diurnal variation of photosynthesis ofMaxim ‘hanyuan’ seedling leaf[J]. Acta Ecologica Sinica, 2016, 36(5): 1321–1330

    [3] 何騰兵, 劉元生, 李天智, 等. 貴州喀斯特峽谷水保經(jīng)濟(jì)植物花椒土壤特性研究[J]. 水土保持學(xué)報(bào), 2000, 14(2): 55–59 HE T B, LIU Y S, LI T Z, et al. Soil characteristics ofwith soil and water conservation benefits in Karst valley of Guizhou Province[J]. Journal of Soil and Water Conservation, 2000, 14(2): 55–59

    [4] 王進(jìn)闖, 潘開文, 吳寧, 等. 花椒農(nóng)林復(fù)合生態(tài)系統(tǒng)的簡化對某些相關(guān)因子的影響[J]. 應(yīng)用與環(huán)境生物學(xué)報(bào), 2005, 11(1): 36–39 WANG J C, PAN K W, WU N, et al. Effect of simplified agro-forest system of Sichuan pepper () on some ecological factors[J]. Chinese Journal of Applied and Environmental Biology, 2005, 11(1): 36–39

    [5] 張立新, 耿增超, 張朝陽, 等. 韓城市花椒園土壤養(yǎng)分狀況及施肥研究[J]. 干旱地區(qū)農(nóng)業(yè)研究, 2003, 21(4): 65–67 ZHANG L X, GENG Z C, ZHANG C Y, et al. Study on soil nutrient condition and fertilizer application in Chinese prickly ash orchard in Hancheng[J]. Agricultural Research in the Arid Areas, 2003, 21(4): 65–67

    [6] 鄭海峰, 趙新峰, 馬巖. 地形和土壤屬性對大豆產(chǎn)量的影響[J]. 土壤通報(bào), 2008, 39(6): 1348–1352 ZHENG H F, ZHAO X F, MA Y. Soybean yields related to topography and soil properties[J]. Chinese Journal of Soil Science, 2008, 39(6): 1348–1352

    [7] 周國朋, 謝志堅(jiān), 曹衛(wèi)東, 等. 稻草高茬-紫云英聯(lián)合還田改善土壤肥力提高作物產(chǎn)量[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2017, 33(23): 157–163 ZHOU G P, XIE Z J, CAO W D, et al. Co-incorporation of high rice stubble and Chinese milk vetch improving soil fertility and yield of rice[J]. Transactions of the CSAE, 2017, 33(23): 157–163

    [8] 于君寶, 劉景雙, 王金達(dá), 等. 典型黑土pH值變化對營養(yǎng)元素有效態(tài)含量的影響研究[J]. 土壤通報(bào), 2003, 34(5): 404–408 YU J B, LIU J S, WANG J D, et al. Effect of soil pH value on the effective content of nutrient elements in typical black soil[J]. Chinese Journal of Soil Science, 2003, 34(5): 404–408

    [9] LI J R, OKIN G S, ALVAREZ L, et al. Effects of wind erosion on the spatial heterogeneity of soil nutrients in two desert grassland communities[J]. Biogeochemistry, 2008, 88(1): 73–88

    [10] 鮑士旦. 土壤農(nóng)業(yè)化學(xué)分析[M]. 第3版. 北京: 中國農(nóng)業(yè)出版社, 2000 BAO S D. Soil Agrochemical Analysis[M]. Third Edition. Beijing: China Agriculture Press, 2000

    [11] GONG W, YAN X Y, WANG J Y, et al. Long-term manuring and fertilization effects on soil organic carbon pools under a wheat-maize cropping system in North China Plain[J]. Plant and Soil, 2009, 314(1/2): 67–76

    [12] 孟凡喬, 吳文良, 辛德惠. 高產(chǎn)農(nóng)田土壤有機(jī)質(zhì)、養(yǎng)分的變化規(guī)律與作物產(chǎn)量的關(guān)系[J]. 植物營養(yǎng)與肥料學(xué)報(bào), 2000, 6(4): 370–374 MENG F Q, WU W L, XIN D H. Changes of soil organic matter and nutrients and their relationship with crop yield in high yield farmland[J]. Plant Nutrition and Fertilizer Science, 2000, 6(4): 370–374

    [13] 范慶鋒, 虞娜, 張玉玲, 等. 設(shè)施蔬菜栽培對土壤陽離子交換性能的影響[J]. 土壤學(xué)報(bào), 2014, 51(5): 1132–1137 FAN Q F, YU N, ZHANG Y L, et al. Effects of vegetable cultivation on soil cation exchange capacity in greenhouse[J]. Acta Pedologica Sinica, 2014, 51(5): 1132–1137

    [14] 李忠意, 白穎艷, 程永毅, 等. 幾種土壤交換性酸測定方法的效果比較[J]. 土壤, 2017, 49(6): 1210–1215 LI Z Y, BAI Y Y, CHENG Y Y, et al. Comparison of soil exchangeable acidities by different measured methods[J]. Soils, 2017, 49(6): 1210–1215

    [15] Kozak M, Kang M S. Note on modern path analysis in application to crop science[J]. Communications in Biometry and Crop Science, 2006, 1(1): 32–34

    [16] 張曉光, 黃標(biāo), 張貝爾, 等. 鹽漬化平原區(qū)玉米產(chǎn)量空間變異與地形關(guān)系研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2012, 43(11): 51–57 ZHANG X G, HUANG B, ZHANG B E, et al. Spatial variability of maize yield and relations to terrain attributes in salinized plain region[J]. Transactions of the Chinese Society for Agricultural Machinery, 2012, 43(11): 51–57

    [17] 雷斯越, 趙文慧, 楊亞輝, 等. 不同坡位植被生長狀況與土壤養(yǎng)分空間分布特征[J]. 水土保持研究, 2019, 26(1): 86–91 LEI S Y, ZHAO W H, YANG Y H, et al. Spatial distribution characteristics of soil nutrients and vegetation growth status in different slopes[J]. Research of Soil and Water Conservation, 2019, 26(1): 86–91

    [18] 余玲, 王彥榮, 毛玉林, 等. 黃土高原地形對苜蓿種子產(chǎn)量和質(zhì)量的影響[J]. 草業(yè)學(xué)報(bào), 2002, 11(4): 62–67 YU L, WANG Y R, MAO Y L, et al. Effect of topography on seed yield and quality of alfalfa grown in the Loess Plateau[J]. Acta Prataculturae Sinica, 2002, 11(4): 62–67

    [19] 李春喜, 韓蕊, 邵云, 等. 黃淮海地區(qū)麥玉周年土壤肥力與玉米產(chǎn)量相關(guān)和通徑分析[J]. 玉米科學(xué), 2019, 27(3): 127–133 LI C X, HAN R, SHAO Y, et al. Correlation and path analysis between soil fertility and maize yield in wheat-maize rotation system of Huang-Huai-Hai region[J]. Journal of Maize Sciences, 2019, 27(3): 127–133

    [20] 彭衛(wèi)福, 呂偉生, 黃山, 等. 土壤肥力對紅壤性水稻土水稻產(chǎn)量和氮肥利用效率的影響[J]. 中國農(nóng)業(yè)科學(xué), 2018, 51(18): 3614–3624 PENG W F, LYU W S, HUANG S, et al. Effects of soil fertility on rice yield and nitrogen use efficiency in a red paddy soil[J]. Scientia Agricultura Sinica, 2018, 51(18): 3614–3624

    [21] 曾祥明, 韓寶吉, 徐芳森, 等. 不同基礎(chǔ)地力土壤優(yōu)化施肥對水稻產(chǎn)量和氮肥利用率的影響[J]. 中國農(nóng)業(yè)科學(xué), 2012, 45(14): 2886–2894 ZENG X M, HAN B J, XU F S, et al. Effect of optimized fertilization on grain yield of rice and nitrogen use efficiency in paddy fields with different basic soil fertilities[J]. Scientia Agricultura Sinica, 2012, 45(14): 2886–2894

    [22] 魯艷紅, 廖育林, 周興, 等. 長期不同施肥對紅壤性水稻土產(chǎn)量及基礎(chǔ)地力的影響[J]. 土壤學(xué)報(bào), 2015, 52(3): 597–606 LU Y H, LIAO Y L, ZHOU X, et al. Effect of long-term fertilization on rice yield and basic soil productivity in red paddy soil under double-rice system[J]. Acta Pedologica Sinica, 2015, 52(3): 597–606

    [23] 湯珧華, 潘建萍, 鄒福生, 等. 上海松柏古樹生長與土壤肥力因子的關(guān)系[J]. 植物營養(yǎng)與肥料學(xué)報(bào), 2017, 23(5): 1402–1408 TANG Y H, PAN J P, ZOU F S, et al. Relationship between soil fertility factors and growth of old pines and cypress in Shanghai[J]. Plant Nutrition and Fertilizer Science, 2017, 23(5): 1402–1408

    [24] 張強(qiáng), 魏欽平, 劉旭東, 等. 北京昌平蘋果園土壤養(yǎng)分、pH與果實(shí)礦質(zhì)營養(yǎng)的多元分析[J]. 果樹學(xué)報(bào), 2011, 28(3): 377–383 ZHANG Q, WEI Q P, LIU X D, et al. Multivariate analysis of relationship between soil nutrients, pH and fruit mineral nutrition in Fuji apple orchards of Changping, Beijing[J]. Journal of Fruit Science, 2011, 28(3): 377–383

    [25] MARTINSEN V, ALLING V, NURIDA N L, et al. pH effects of the addition of three biochars to acidic Indonesian mineral soils[J]. Soil Science and Plant Nutrition, 2015, 61(5): 821–834

    [26] 王倩, 李軍, 寧芳, 等. 渭北旱作麥田長期保護(hù)性耕作土壤肥力特征綜合評價(jià)[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2018, 29(9): 2925–2934 WANG Q, LI J, NING F, et al. Comprehensive assessment of soil fertility characteristics under different long-term conservation tillages of wheat field in Weibei Highland, China[J]. Chinese Journal of Applied Ecology, 2018, 29(9): 2925–2934

    [27] 王飛, 何春梅, 李清華, 等. 外源鈣水平與花生下針期不同土壤水分對植株生理特性的影響[J]. 植物營養(yǎng)與肥料學(xué)報(bào), 2013, 19(3): 623–631 WANG F, HE C M, LI Q H, et al. Effects of exogenous calcium and soil moisture at acicula forming stage of peanut on some physiological characteristics of plants[J]. Journal of Plant Nutrition and Fertilizer, 2013, 19(3): 623–631

    [28] 張大庚, 栗杰, 劉慧. 外源低分子量有機(jī)酸對土壤中鈣素遷移特征的影響[J]. 水土保持學(xué)報(bào), 2015, 29(5): 152–155 ZHANG D G, LI J, LIU H. Effect of exogenous low molecular weight organic acids on calcium migration in soil[J]. Journal of Soil and Water Conservation, 2015, 29(5): 152–155

    [29] 董珊珊, 李寧冉, 楊海燕, 等. 藍(lán)莓根系對土壤錳脅迫的生理影響[J]. 南京林業(yè)大學(xué)學(xué)報(bào): 自然科學(xué)版: 2019, 43(3): 169-173 DONG S S, LI N R, YANG H Y, et al. Physiological response of blueberry roots to Mn2+stress in soil[J]. Journal of Nanjing Forestry University:Natural Sciences Edition, 2019, 43(3): 169-173

    [30] 盧曉佩, 劉秀紅, 董肖昌, 等. 南豐和衢州柑橘產(chǎn)區(qū)土壤及葉片硼含量狀況分析[J]. 華中農(nóng)業(yè)大學(xué)學(xué)報(bào), 2015, 34(6): 61–65 LU X P, LIU X H, DONG X C, et al. Analyzing boron contents in soil and leaf oforchards in the Nanfeng and Quzhou[J]. Journal of Huazhong Agricultural University, 2015, 34(6): 61–65

    [31] 冀建華, 李絮花, 劉秀梅, 等. 硅鈣鉀鎂肥對南方稻田土壤酸度的改良作用[J]. 土壤學(xué)報(bào), 2019, 56(4): 895–906 JI J H, LI X H, LIU X M, et al. Effect of Si-Ca-K-Mg fertilizer remedying acid paddy soil in South China[J]. Acta Pedologica Sinica, 2019, 56(4): 895–906

    [32] Kirkham J M, Rowe B A, Doyle R B. Persistent improvements in the structure and hydraulic conductivity of a ferrosol due to liming[J]. Australian Journal of Soil Research, 2007, 45(3): 218–223

    [33] SARNO, IIJIMA M, LUMBANRAJA J, et al. Soil chemical properties of an Indonesian red acid soil as affected by land use and crop management[J]. Soil and Tillage Research, 2004, 76(2): 115–124

    Relationships ofyield with topography and soil fertility in Jiulongpo area*

    YANG Shixi1, LYU Guangbin1, HUANG Yun1**, XIANG Huahui2, WANG Zhengyin1

    (1. College of Resources and Environmental Sciences, Southwest University, Chongqing 400716, China; 2. Extension Centre of Agricultural Techniques of Jiulongpo District in Chongqing, Chongqing 400039, China)

    The relationships of soil nutrient conditions and topography with production in Jiulongpoplanting areas was studied through a field investigation-based research and lab analysis to provide a theoretical basis for scientific and effective fertilization practices. The characteristics of altitude and slope, and soil pH, organic matter content, macro- and micro-elements content, cation exchange capacity (CEC) and exchange acidity as well asyeild in low-, medium and high-yield areas ofwere investigated. The results showed thatwas generally planted at altitudes between 200 and 500 meters in Jiulongpo, and the altitude of the high-yield area was approximately 300 meters. The slope from low- to high-yield areas appeared to be slight increase but not obvious. The soil was acidic, with a pH < 6.5; and soil fertility was high, with large differences in nutrient factors. CEC and available phosphorus, calcium, magnesium, iron, manganese, copper, and zinc were abundant at 27.2 cmol(+)?kg-1, 35.2 mg?kg-1, 3 289.8 mg?kg-1, 271.8 mg?kg-1, 48.6 mg?kg-1, 62.1 mg?kg-1, 1.5 mg?kg-1, and 4.5 mg?kg-1, respectively. Organic matter and available nitrogen, potassium and exchange acid were at a moderate levels of 19.1 mg?kg-1, 114.9 mg?kg-1, 107.0 mg?kg-1and 8.1 cmol(+)?kg-1, respectively. The content of water-soluble boron was a relatively low 0.28 mg?kg-1.yield exhibited a positive correlation with available calcium, CEC, pH, available manganese, and water-soluble boron. Available calcium, CEC, exchange acidity, and available copper, iron, and zinc were the main influencing factors on yield. An optimal linear regression equation of available calcium (6) and yield () was established based on stepwise regression analysis (= 11.693 + 0.0036). In conclusion, soil nutrient in Jiulongpoarea was imbalanced seriously, more attention should be paid to balancing fertilization, increasing organic fertilizer, improving the physical and chemical properties of soil, and relieving soil acidification.

    ; Topography; Soil fertility; Yield; Nutrient balance

    S158.3

    2096-6237(2019)12-1823-10

    10.13930/j.cnki.cjea.190127

    * 國家自然科學(xué)基金重點(diǎn)項(xiàng)目(41530855)和國際植物營養(yǎng)研究所(IPNI)項(xiàng)目(2013-Chongqing-02)資助

    黃云, 主要研究方向?yàn)橹参餇I養(yǎng)與施肥。E-mail: huangyun1622@126.com

    楊仕曦, 研究方向?yàn)橹参餇I養(yǎng)與土壤酸化改良技術(shù)。E-mail: 751743470@qq.com

    2019-02-21

    2019-07-20

    * This study was supported by the National Natural Science Foundation of China (41530855) and the International Plant Nutrition Institute (IPNI) Funded Project (2013-Chongqing-02).

    , E-mail: huangyun1622@126.com

    Feb. 21, 2019;

    Jul. 20, 2019

    楊仕曦, 呂廣斌, 黃云, 向華輝, 王正銀. 九龍坡花椒種植區(qū)地形、土壤肥力與花椒產(chǎn)量的關(guān)系[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào)(中英文), 2019, 27(12): 1823-1832

    YANG S X, LYU G B, HUANG Y, XIANG H H, WANG Z Y. Relationships ofyield with topography and soil fertility in Jiulongpo area[J]. Chinese Journal of Eco-Agriculture, 2019, 27(12): 1823-1832

    猜你喜歡
    種植區(qū)土壤肥力產(chǎn)區(qū)
    C市主要草莓種植區(qū)土壤重金屬鎘、鉛現(xiàn)狀調(diào)查
    不同種植區(qū)隴東苜蓿營養(yǎng)價(jià)值的比較研究
    中國飼料(2021年17期)2021-11-02 08:15:24
    葡萄酒的產(chǎn)區(qū)品牌
    收藏界(2019年2期)2019-10-12 08:27:04
    葡萄酒的產(chǎn)區(qū)品牌(下)
    收藏界(2019年3期)2019-10-10 03:16:48
    草莓種植區(qū)土壤中典型有機(jī)氮化合物的分布及來源
    葡萄酒產(chǎn)區(qū)與『風(fēng)土』
    收藏界(2018年1期)2018-10-10 05:23:20
    拉薩市土壤肥力變化趨勢與改良策略
    新疆產(chǎn)區(qū)有機(jī)甜菜栽培技術(shù)探討
    中國糖料(2016年1期)2016-12-01 06:49:04
    河北昌黎縣葡萄種植區(qū)農(nóng)業(yè)地球化學(xué)特征
    安吉白茶園土壤肥力現(xiàn)狀分析
    茶葉(2014年4期)2014-02-27 07:05:15
    99热6这里只有精品| 乱人视频在线观看| 亚洲,欧美,日韩| 一边摸一边抽搐一进一小说| 欧美+日韩+精品| 日韩欧美在线乱码| 亚洲精华国产精华精| 特级一级黄色大片| 非洲黑人性xxxx精品又粗又长| 最近中文字幕高清免费大全6 | 在现免费观看毛片| 嫩草影视91久久| 欧美成人一区二区免费高清观看| 国产亚洲精品av在线| 国产在视频线在精品| av在线观看视频网站免费| or卡值多少钱| 久久天躁狠狠躁夜夜2o2o| 韩国av在线不卡| 久久99热这里只有精品18| 日本色播在线视频| 成年女人看的毛片在线观看| 亚洲国产欧美人成| 日韩欧美在线二视频| 人人妻,人人澡人人爽秒播| 亚洲国产色片| 成人性生交大片免费视频hd| 一区二区三区四区激情视频 | 免费在线观看影片大全网站| 性色avwww在线观看| 三级男女做爰猛烈吃奶摸视频| 国产精品日韩av在线免费观看| 中文字幕免费在线视频6| 国产精品久久久久久亚洲av鲁大| 男人和女人高潮做爰伦理| av.在线天堂| 久久久久九九精品影院| 国产乱人伦免费视频| 国产亚洲精品av在线| 日本三级黄在线观看| 亚洲精华国产精华精| 国产真实伦视频高清在线观看 | 午夜免费成人在线视频| 内射极品少妇av片p| 99久久无色码亚洲精品果冻| 白带黄色成豆腐渣| 99久久无色码亚洲精品果冻| 国产高清不卡午夜福利| 久久午夜福利片| av福利片在线观看| 成人无遮挡网站| 女同久久另类99精品国产91| videossex国产| 亚洲第一区二区三区不卡| 国产人妻一区二区三区在| 精品人妻1区二区| 一卡2卡三卡四卡精品乱码亚洲| 日韩,欧美,国产一区二区三区 | 久久欧美精品欧美久久欧美| 啦啦啦啦在线视频资源| 热99在线观看视频| 精品久久久久久久久久免费视频| 国产av不卡久久| netflix在线观看网站| 少妇的逼好多水| 精品日产1卡2卡| 欧美不卡视频在线免费观看| 午夜福利视频1000在线观看| 一进一出抽搐gif免费好疼| 免费看av在线观看网站| 亚洲av.av天堂| 在线国产一区二区在线| 色播亚洲综合网| 1000部很黄的大片| 久久草成人影院| 国产男人的电影天堂91| 高清毛片免费观看视频网站| 日本色播在线视频| 一个人看视频在线观看www免费| 免费av毛片视频| 国产精品人妻久久久久久| 校园人妻丝袜中文字幕| 国产精品一区二区三区四区免费观看 | 一进一出抽搐动态| 日韩精品有码人妻一区| 国产精品野战在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产熟女欧美一区二区| 在线看三级毛片| 国产精品久久久久久精品电影| 真人做人爱边吃奶动态| 少妇猛男粗大的猛烈进出视频 | 婷婷丁香在线五月| 精品久久久久久久久av| 欧美激情久久久久久爽电影| 精品一区二区三区视频在线| 午夜a级毛片| 国产亚洲精品av在线| 成年人黄色毛片网站| 美女免费视频网站| 美女被艹到高潮喷水动态| 欧美成人a在线观看| 日韩欧美国产一区二区入口| 久久久久精品国产欧美久久久| 日韩强制内射视频| 亚洲电影在线观看av| 啦啦啦啦在线视频资源| 床上黄色一级片| 婷婷精品国产亚洲av| 又爽又黄无遮挡网站| 少妇的逼好多水| 欧美区成人在线视频| 最新在线观看一区二区三区| 国产精品一区二区三区四区久久| 99久久九九国产精品国产免费| 午夜激情欧美在线| 国产精品一区二区三区四区久久| 免费大片18禁| 亚洲性夜色夜夜综合| 国产成人福利小说| 美女黄网站色视频| 男女做爰动态图高潮gif福利片| 亚洲av.av天堂| 日本黄色视频三级网站网址| 九九爱精品视频在线观看| 欧美zozozo另类| 无遮挡黄片免费观看| 国产精品人妻久久久久久| 熟女人妻精品中文字幕| 亚洲自偷自拍三级| 女同久久另类99精品国产91| 精品99又大又爽又粗少妇毛片 | 中文字幕人妻熟人妻熟丝袜美| 黄色欧美视频在线观看| 天堂影院成人在线观看| 亚洲久久久久久中文字幕| 欧美成人免费av一区二区三区| 国产aⅴ精品一区二区三区波| 国产精品久久电影中文字幕| 亚洲一区高清亚洲精品| 一个人免费在线观看电影| 少妇人妻一区二区三区视频| 两性午夜刺激爽爽歪歪视频在线观看| 国产亚洲欧美98| 深爱激情五月婷婷| 午夜福利在线在线| 成人精品一区二区免费| 简卡轻食公司| 国国产精品蜜臀av免费| 一区二区三区高清视频在线| 少妇被粗大猛烈的视频| 国产伦在线观看视频一区| 麻豆精品久久久久久蜜桃| 欧美一区二区国产精品久久精品| 悠悠久久av| 啦啦啦韩国在线观看视频| 麻豆成人午夜福利视频| 精品久久久久久久久av| 日韩人妻高清精品专区| 成人欧美大片| 能在线免费观看的黄片| 欧美精品国产亚洲| 国产高清视频在线播放一区| 精品国内亚洲2022精品成人| 精品欧美国产一区二区三| 日韩欧美三级三区| 午夜福利在线观看免费完整高清在 | 桃色一区二区三区在线观看| or卡值多少钱| 色在线成人网| 欧美激情国产日韩精品一区| 国产单亲对白刺激| 在线观看美女被高潮喷水网站| 午夜影院日韩av| 熟妇人妻久久中文字幕3abv| 特大巨黑吊av在线直播| 99国产精品一区二区蜜桃av| 成人国产麻豆网| 国产高清不卡午夜福利| 国产精品一区二区免费欧美| 亚洲av免费高清在线观看| 少妇的逼好多水| 欧美又色又爽又黄视频| 女生性感内裤真人,穿戴方法视频| 我的老师免费观看完整版| 亚洲av中文字字幕乱码综合| 很黄的视频免费| 精品免费久久久久久久清纯| 男女做爰动态图高潮gif福利片| 国产视频一区二区在线看| 观看免费一级毛片| 精品一区二区三区视频在线| 九九热线精品视视频播放| 日韩欧美精品v在线| 国产精品一区二区三区四区免费观看 | 精华霜和精华液先用哪个| 精品福利观看| av在线老鸭窝| 中国美白少妇内射xxxbb| 白带黄色成豆腐渣| a级一级毛片免费在线观看| 很黄的视频免费| 久久6这里有精品| 97超视频在线观看视频| 老熟妇仑乱视频hdxx| 搡老熟女国产l中国老女人| 91精品国产九色| 日韩高清综合在线| 少妇的逼好多水| 精品人妻视频免费看| 一个人看视频在线观看www免费| 国产精品永久免费网站| 久久久色成人| 免费大片18禁| 免费av毛片视频| 大型黄色视频在线免费观看| 中国美白少妇内射xxxbb| 日本三级黄在线观看| 在线免费十八禁| 日本免费a在线| 九九久久精品国产亚洲av麻豆| 免费观看在线日韩| 搡老熟女国产l中国老女人| 日韩欧美在线二视频| 成人av在线播放网站| 天美传媒精品一区二区| 久久久久久久久久黄片| 午夜免费成人在线视频| av在线观看视频网站免费| 国产精品爽爽va在线观看网站| 女人十人毛片免费观看3o分钟| 99riav亚洲国产免费| 我要看日韩黄色一级片| 久久久久久久久大av| 欧美又色又爽又黄视频| 熟女电影av网| 熟妇人妻久久中文字幕3abv| 午夜精品一区二区三区免费看| 又爽又黄无遮挡网站| 国产欧美日韩一区二区精品| 校园人妻丝袜中文字幕| 哪里可以看免费的av片| 久久精品久久久久久噜噜老黄 | 亚洲专区中文字幕在线| 亚洲av日韩精品久久久久久密| aaaaa片日本免费| 国产免费男女视频| 精品一区二区免费观看| 春色校园在线视频观看| 国产久久久一区二区三区| 国产精品一区二区性色av| 乱人视频在线观看| 精品日产1卡2卡| .国产精品久久| 国国产精品蜜臀av免费| 欧美色欧美亚洲另类二区| 午夜免费男女啪啪视频观看 | 亚洲人成网站在线播放欧美日韩| 国产91精品成人一区二区三区| 美女cb高潮喷水在线观看| 中文在线观看免费www的网站| av在线蜜桃| 99热只有精品国产| 18禁黄网站禁片午夜丰满| 国产精品国产高清国产av| 国产高清视频在线观看网站| 欧美3d第一页| 在线天堂最新版资源| 香蕉av资源在线| 18禁黄网站禁片午夜丰满| 老熟妇乱子伦视频在线观看| 九色成人免费人妻av| 成人鲁丝片一二三区免费| 看十八女毛片水多多多| 国产91精品成人一区二区三区| 韩国av在线不卡| 亚洲av二区三区四区| 亚洲内射少妇av| 婷婷丁香在线五月| 精品一区二区三区视频在线观看免费| 国产成人福利小说| 老司机福利观看| 亚洲第一区二区三区不卡| 国内毛片毛片毛片毛片毛片| 国产又黄又爽又无遮挡在线| av中文乱码字幕在线| 国产精品久久久久久av不卡| 级片在线观看| 亚洲av中文av极速乱 | 午夜视频国产福利| 最新中文字幕久久久久| bbb黄色大片| 成人国产综合亚洲| 精品久久久久久久久av| 精品无人区乱码1区二区| 一个人看视频在线观看www免费| 国产精品一及| 国产精品美女特级片免费视频播放器| 亚洲va日本ⅴa欧美va伊人久久| 国产精品免费一区二区三区在线| 国产欧美日韩一区二区精品| 国产av麻豆久久久久久久| 国产高清视频在线播放一区| 免费观看人在逋| 色尼玛亚洲综合影院| 国内少妇人妻偷人精品xxx网站| 嫩草影院精品99| 桃色一区二区三区在线观看| 欧美绝顶高潮抽搐喷水| avwww免费| 国产美女午夜福利| 日本熟妇午夜| 哪里可以看免费的av片| 特大巨黑吊av在线直播| 日韩欧美 国产精品| 噜噜噜噜噜久久久久久91| 午夜激情欧美在线| 亚洲精华国产精华液的使用体验 | 蜜桃久久精品国产亚洲av| 99视频精品全部免费 在线| 97超视频在线观看视频| 日本欧美国产在线视频| av在线亚洲专区| 美女黄网站色视频| 中文字幕精品亚洲无线码一区| 亚洲国产日韩欧美精品在线观看| 男人舔奶头视频| 免费人成视频x8x8入口观看| 亚洲三级黄色毛片| 又粗又爽又猛毛片免费看| 狠狠狠狠99中文字幕| 国产色爽女视频免费观看| 亚洲 国产 在线| 亚洲黑人精品在线| 精品国内亚洲2022精品成人| 美女大奶头视频| 日韩高清综合在线| 在线观看舔阴道视频| 久久精品国产清高在天天线| 国产中年淑女户外野战色| 婷婷精品国产亚洲av| 九九热线精品视视频播放| 国产精品久久电影中文字幕| 99久久久亚洲精品蜜臀av| 在线观看美女被高潮喷水网站| 男女那种视频在线观看| 亚洲精品一区av在线观看| 亚洲va在线va天堂va国产| 日韩在线高清观看一区二区三区 | 性色avwww在线观看| 婷婷精品国产亚洲av| 99久久无色码亚洲精品果冻| 少妇高潮的动态图| 夜夜爽天天搞| 精品不卡国产一区二区三区| 国产精华一区二区三区| av.在线天堂| 亚洲欧美清纯卡通| av在线蜜桃| 亚洲欧美日韩高清专用| 精品人妻偷拍中文字幕| 亚洲欧美激情综合另类| 日韩欧美精品v在线| 欧美3d第一页| 久久精品国产自在天天线| 亚洲欧美精品综合久久99| 国产 一区精品| 美女高潮的动态| 色在线成人网| 成人欧美大片| 狂野欧美激情性xxxx在线观看| 波多野结衣巨乳人妻| 丰满乱子伦码专区| 中国美白少妇内射xxxbb| 黄色视频,在线免费观看| 色综合婷婷激情| 日韩强制内射视频| 欧美激情国产日韩精品一区| 亚洲aⅴ乱码一区二区在线播放| 国产毛片a区久久久久| 18禁在线播放成人免费| 天天一区二区日本电影三级| 欧美zozozo另类| 精品福利观看| 一本久久中文字幕| 国产高潮美女av| 国产成人aa在线观看| 97碰自拍视频| 干丝袜人妻中文字幕| 日本色播在线视频| 日本在线视频免费播放| 亚洲av美国av| 天堂av国产一区二区熟女人妻| 国产成年人精品一区二区| 人妻丰满熟妇av一区二区三区| 国产男人的电影天堂91| 久久久久久久久中文| 最近视频中文字幕2019在线8| 免费av观看视频| 美女xxoo啪啪120秒动态图| 欧美+日韩+精品| 狠狠狠狠99中文字幕| 免费无遮挡裸体视频| 欧美xxxx黑人xx丫x性爽| 国产熟女欧美一区二区| 99精品久久久久人妻精品| 亚洲专区中文字幕在线| 啦啦啦啦在线视频资源| 3wmmmm亚洲av在线观看| 又爽又黄a免费视频| 婷婷六月久久综合丁香| 淫妇啪啪啪对白视频| 久久精品影院6| 波多野结衣巨乳人妻| 欧美性猛交╳xxx乱大交人| 国产午夜精品论理片| 亚洲七黄色美女视频| 网址你懂的国产日韩在线| 男女视频在线观看网站免费| 日本黄色视频三级网站网址| 亚洲成人久久爱视频| 精品人妻1区二区| 又爽又黄无遮挡网站| 男人舔女人下体高潮全视频| 乱码一卡2卡4卡精品| 天堂动漫精品| 午夜视频国产福利| 亚洲av第一区精品v没综合| 婷婷亚洲欧美| 国产亚洲精品久久久com| 老熟妇仑乱视频hdxx| 日本撒尿小便嘘嘘汇集6| 深夜a级毛片| 欧美高清性xxxxhd video| 99久久中文字幕三级久久日本| 亚洲三级黄色毛片| 欧美一区二区亚洲| 欧美不卡视频在线免费观看| 给我免费播放毛片高清在线观看| 此物有八面人人有两片| 又爽又黄无遮挡网站| 搡女人真爽免费视频火全软件 | 91在线观看av| 免费搜索国产男女视频| 国产极品精品免费视频能看的| 久9热在线精品视频| 婷婷精品国产亚洲av在线| 99久国产av精品| 免费av毛片视频| 免费在线观看日本一区| 乱人视频在线观看| h日本视频在线播放| 美女cb高潮喷水在线观看| 午夜福利高清视频| 成人鲁丝片一二三区免费| 国产高清激情床上av| 国产大屁股一区二区在线视频| 国产av不卡久久| 亚洲国产日韩欧美精品在线观看| 欧美激情在线99| 国产熟女欧美一区二区| 三级毛片av免费| 日韩欧美精品v在线| 日本在线视频免费播放| 有码 亚洲区| 久久久久久久久久成人| 最后的刺客免费高清国语| 成人二区视频| 老熟妇仑乱视频hdxx| 成人鲁丝片一二三区免费| 国产高清激情床上av| 午夜精品一区二区三区免费看| 日韩强制内射视频| 亚洲色图av天堂| 久久久色成人| 亚洲成人精品中文字幕电影| 久久久久精品国产欧美久久久| 九色成人免费人妻av| 亚洲,欧美,日韩| videossex国产| 老熟妇仑乱视频hdxx| 男人狂女人下面高潮的视频| 久久久久久久久大av| 国产伦精品一区二区三区四那| 国产黄色小视频在线观看| 欧美区成人在线视频| 一进一出抽搐动态| 男女下面进入的视频免费午夜| 精品一区二区三区av网在线观看| 深爱激情五月婷婷| 九九在线视频观看精品| 国产精品嫩草影院av在线观看 | a级毛片a级免费在线| 少妇熟女aⅴ在线视频| 狠狠狠狠99中文字幕| 亚洲欧美日韩东京热| 亚洲七黄色美女视频| 亚洲电影在线观看av| 免费在线观看成人毛片| www.色视频.com| 精品久久久久久久人妻蜜臀av| 成人午夜高清在线视频| 国产免费一级a男人的天堂| 亚洲一区高清亚洲精品| 美女黄网站色视频| 搡老熟女国产l中国老女人| videossex国产| 欧美+亚洲+日韩+国产| 中文字幕免费在线视频6| 亚洲熟妇中文字幕五十中出| 一个人看视频在线观看www免费| 男人的好看免费观看在线视频| 99国产极品粉嫩在线观看| 日本三级黄在线观看| 成年女人永久免费观看视频| 伦理电影大哥的女人| 99精品久久久久人妻精品| 精品久久久久久久久亚洲 | 超碰av人人做人人爽久久| 国产精品亚洲一级av第二区| 韩国av在线不卡| 国产免费av片在线观看野外av| 午夜视频国产福利| 亚洲五月天丁香| 老熟妇仑乱视频hdxx| 亚洲av电影不卡..在线观看| 久久人人爽人人爽人人片va| 久久精品国产亚洲av涩爱 | 亚洲精华国产精华精| 日韩人妻高清精品专区| 亚洲精品乱码久久久v下载方式| 少妇的逼好多水| 国产精品国产高清国产av| 国产淫片久久久久久久久| 亚洲成人中文字幕在线播放| 中文字幕久久专区| 国产高清有码在线观看视频| 91久久精品电影网| 欧美高清性xxxxhd video| 国产精品日韩av在线免费观看| 国产乱人伦免费视频| 51国产日韩欧美| 久久99热这里只有精品18| 午夜日韩欧美国产| 亚洲va在线va天堂va国产| 2021天堂中文幕一二区在线观| 观看免费一级毛片| 日本免费a在线| 色哟哟哟哟哟哟| 中出人妻视频一区二区| АⅤ资源中文在线天堂| 老熟妇乱子伦视频在线观看| 欧美xxxx性猛交bbbb| 99久久精品热视频| 日日摸夜夜添夜夜添av毛片 | 日韩人妻高清精品专区| 久久精品国产鲁丝片午夜精品 | 在线a可以看的网站| 国产淫片久久久久久久久| 午夜a级毛片| 国产淫片久久久久久久久| av国产免费在线观看| 窝窝影院91人妻| 性插视频无遮挡在线免费观看| 久久欧美精品欧美久久欧美| 中文在线观看免费www的网站| 国产精品女同一区二区软件 | 性欧美人与动物交配| 99热只有精品国产| 成人高潮视频无遮挡免费网站| 能在线免费观看的黄片| 日韩欧美一区二区三区在线观看| 韩国av一区二区三区四区| 国产亚洲精品久久久com| 婷婷色综合大香蕉| 真人做人爱边吃奶动态| 日韩高清综合在线| 美女大奶头视频| 天堂网av新在线| 校园人妻丝袜中文字幕| ponron亚洲| 又爽又黄a免费视频| 亚洲av一区综合| 999久久久精品免费观看国产| 特大巨黑吊av在线直播| 国产一区二区亚洲精品在线观看| 小说图片视频综合网站| 男人狂女人下面高潮的视频| 亚洲精品粉嫩美女一区| 麻豆成人av在线观看| 久久久国产成人精品二区| 欧美精品国产亚洲| 黄色日韩在线| 欧美3d第一页| 有码 亚洲区| 国产老妇女一区| 中出人妻视频一区二区| 欧美一区二区国产精品久久精品| 在线天堂最新版资源| 婷婷丁香在线五月| 九九热线精品视视频播放| 日本黄色视频三级网站网址| 成熟少妇高潮喷水视频| 亚洲狠狠婷婷综合久久图片| 国产一级毛片七仙女欲春2| 欧美日韩亚洲国产一区二区在线观看| or卡值多少钱| 身体一侧抽搐| 国产精品永久免费网站| 日韩欧美免费精品| 精品人妻偷拍中文字幕| 九色国产91popny在线| 狠狠狠狠99中文字幕| 久久久午夜欧美精品| 人妻丰满熟妇av一区二区三区|