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

    湖南省柑橘園土壤pH 值空間變異特征、成因及養(yǎng)分有效性

    2024-07-09 00:00:00徐章倩周衛(wèi)軍崔浩杰周雨舟商貴鐸李敏岳晴晴

    摘要: 【目的】研究湖南省柑橘園土壤pH 值的空間變異特征、影響因素及主要養(yǎng)分的有效性,為減緩柑橘園土壤酸化,提升土壤肥力提供科學(xué)依據(jù)?!痉椒ā勘狙芯恳院鲜?8 個(gè)柑橘重點(diǎn)縣的15722 個(gè)柑橘園土壤為對(duì)象,利用地統(tǒng)計(jì)學(xué)、相關(guān)性分析、結(jié)構(gòu)方程模型等方法,探究了湖南省柑橘園土壤pH 值的空間變異特征,分析了土壤pH 值對(duì)地理、氣候、土壤性質(zhì)等因子的響應(yīng),并研究了土壤pH 值對(duì)土壤氮、磷、鉀有效性的影響?!窘Y(jié)果】湖南省柑橘園土壤pH 值范圍為3.00~8.50,均值為5.58,變異系數(shù)為18.06%,空間變異為中等程度。最適合柑橘生長(zhǎng)的弱酸性土壤面積占44.35%,而酸性、強(qiáng)酸性、中性和堿性土壤分別占48.92%、1.69%、4.95% 和0.08%。通過Pearson 相關(guān)性分析發(fā)現(xiàn),土壤pH 值與降水量呈顯著正相關(guān)(Plt;0.05),與速效鉀含量呈極顯著正相關(guān)(Plt;0.01),而與氣溫、堿解氮和有效磷含量呈極顯著負(fù)相關(guān)(Plt;0.01)。結(jié)構(gòu)方程模型模擬結(jié)果表明,土壤類型與成土母質(zhì)是解釋土壤pH 值空間變異的主導(dǎo)因子。除少量強(qiáng)堿性土壤中堿解氮和有效磷含量處于較低水平外,柑橘園土壤各酸度等級(jí)土壤中的堿解氮、有效磷和速效鉀含量皆處于適宜水平?!窘Y(jié)論】由于特有的土壤類型和成土母質(zhì),湖南省柑橘園超過一半(50.61%) 的土壤處于過酸狀態(tài)(酸性和強(qiáng)酸性)。湖南省柑橘園土壤速效氮、磷、鉀養(yǎng)分含量基本處于適宜水平,因此,應(yīng)重視改善土壤酸性。

    關(guān)鍵詞: 柑橘園土壤; pH 值; 空間特征; 成土母質(zhì); 地理因素; 速效養(yǎng)分

    湖南省是中國(guó)柑橘種植適宜區(qū),2018 年全省柑橘種植面積為38.4 萬hm2,產(chǎn)量達(dá)528.6 萬t,居全國(guó)第二,在全國(guó)乃至全世界都享有盛譽(yù)[1]。土壤酸堿度是土壤肥力和柑橘生長(zhǎng)的限制因子之一,在土壤質(zhì)量、土壤養(yǎng)分有效性和植物根系對(duì)土壤養(yǎng)分的吸收能力等方面起著至關(guān)重要的作用,從而影響柑橘品質(zhì)和產(chǎn)量[2?3]。研究發(fā)現(xiàn)中國(guó)柑橘園的土壤酸化普遍發(fā)生在湖南省[2]、廣東省[4]和福建省[5]等主要柑橘產(chǎn)區(qū)。其中,柑橘園土壤酸化導(dǎo)致土壤養(yǎng)分有效性降低[6],造成湖南省大多數(shù)果園土壤營(yíng)養(yǎng)不平衡[7]。大量研究表明,海拔、地形、氣候、成土母質(zhì)與土壤類型等因素均對(duì)土壤pH 值有顯著影響[8?10]。劉雪松等[11]對(duì)贛南地區(qū)土壤酸堿度空間變異的研究發(fā)現(xiàn),成土母質(zhì)是影響土壤pH 值的主要因子,表現(xiàn)為沉積巖發(fā)育的土壤pH 值顯著高于火成巖。Srivastava 等[12]對(duì)全球柑橘園土壤肥力退化與植物營(yíng)養(yǎng)狀況的統(tǒng)計(jì)分析表明,土壤類型對(duì)柑橘園土壤pH 值的影響較大。Wang 等[13]對(duì)三峽水庫區(qū)秭歸縣柑橘園的調(diào)查發(fā)現(xiàn),果園各季節(jié)土壤pH 值均隨海拔升高而降低。與此同時(shí),已有的研究多集中在我國(guó)村級(jí)或縣域等較小的空間區(qū)域,對(duì)較大規(guī)模柑橘園產(chǎn)區(qū)土壤pH 值空間變異和成因分析的研究很少。因此,探究湖南省土壤pH 值的空間變異及其影響因素對(duì)柑橘園土壤改良有重大意義。

    柑橘是一種常綠多年生植物,對(duì)氮(N)、磷(P)和鉀(K) 等大量營(yíng)養(yǎng)元素有很高的需求量,土壤缺氮抑制葉片和莖稈的生長(zhǎng),導(dǎo)致果實(shí)產(chǎn)量逐漸下降[14]。缺磷導(dǎo)致柑橘產(chǎn)量降低[15]。缺鉀造成柑橘果實(shí)葡萄糖、果糖和可溶性糖含量降低[16]。曹勝等[2]和Chen 等[7]研究表明,柑橘園土壤堿解氮、速效鉀與有效磷的含量受土壤pH 值的影響較大,導(dǎo)致多數(shù)果園土壤大量元素的有效態(tài)養(yǎng)分含量處于缺乏狀態(tài)[17]。唐志文等[18]在新寧縣采集并調(diào)查了402 個(gè)土壤樣品,通過相關(guān)分析發(fā)現(xiàn)土壤速效鉀含量隨土壤pH 值的升高而增加,但土壤有效磷含量隨土壤pH 值的升高而降低。所以,揭示湖南省柑橘園土壤速效養(yǎng)分對(duì)土壤pH 值的響應(yīng)機(jī)制,可以為提升柑橘園土壤肥力、增加速效養(yǎng)分有效性與化肥減施增效提供依據(jù)。

    因此,本研究以湖南省38 個(gè)重點(diǎn)柑橘縣(市) 的15722 個(gè)柑橘園土壤為研究對(duì)象,從地形因子、氣候條件和成土因子3 個(gè)方向,基于克里金插值、相關(guān)性分析和結(jié)構(gòu)方程模型分析,揭示湖南省柑橘主產(chǎn)區(qū)土壤pH 值的空間變異性和驅(qū)動(dòng)因子,探尋柑橘園速效養(yǎng)分與土壤pH 值的關(guān)系,為柑橘園土壤管理、土壤肥力提升、柑橘果品質(zhì)量提高提供科學(xué)指導(dǎo)與理論依據(jù),對(duì)農(nóng)業(yè)生態(tài)地理分區(qū)的發(fā)展和農(nóng)業(yè)宏觀決策的制定具有重要意義。

    1 材料與方法

    1.1 研究區(qū)域概況

    湖南省地處中國(guó)中部,東經(jīng)108°47′~114°15′,北緯24°38′~30°08′,擁有14 個(gè)地州市,122 個(gè)縣級(jí)區(qū)劃,總面積約為21.18 萬km2。境內(nèi)屬于亞熱帶季風(fēng)氣候,年平均氣溫為16℃~19℃,年平均降水量為1300~1600 mm,年日照為1300~1800 h,無霜期為253~311 天,冬季1 月平均溫度在4℃ 以上,日均溫度在0℃ 以下的天數(shù)不到10 天;海拔500~1000 m,地勢(shì)呈三面環(huán)山、朝北開口的馬蹄形地貌。

    湖南省柑橘產(chǎn)區(qū)集中分布在常德、湘西、懷化、邵陽、永州、郴州、張家界7 個(gè)市州38 個(gè)重點(diǎn)縣(國(guó)家柑橘優(yōu)勢(shì)區(qū)重點(diǎn)縣25 個(gè)),本研究依據(jù)2007年《湖南省各縣(市區(qū)) 耕地地力評(píng)價(jià)區(qū)域分布表》將38 個(gè)柑橘重點(diǎn)縣劃分為湘中湘東區(qū)、湘南區(qū)、湘西南區(qū)和湘北洞庭湖區(qū)4 個(gè)地區(qū)(圖1A)。

    1.2 樣品采集與分析

    本研究依托國(guó)家柑橘體系與測(cè)土配方施肥與耕地質(zhì)量評(píng)價(jià)項(xiàng)目,采集并收集了2015—2021 年的湖南省4 個(gè)區(qū)域38 個(gè)柑橘重點(diǎn)縣的柑橘園土壤樣點(diǎn)數(shù)據(jù)共計(jì)15722 個(gè)(表1、圖1B)。按照每個(gè)果園具體面積,隨機(jī)多點(diǎn)采集 0—40 cm 土壤樣品,混勻后采用四分法分取1 kg 左右土樣帶回實(shí)驗(yàn)室風(fēng)干,除去根、葉、石塊,磨細(xì),分別過2 mm 和0.149 mm篩,然后進(jìn)行土壤相關(guān)指標(biāo)的分析。

    土壤樣品測(cè)定參照鮑士旦[19]的方法。土壤pH 值采用電位法測(cè)定;土壤堿解氮采用擴(kuò)散法測(cè)定;土壤有效磷采用Olsen 法測(cè)定;土壤速效鉀采用乙酸銨浸提—火焰光度法測(cè)定。

    1.3 數(shù)據(jù)源

    選擇Landsat 8 遙感影像為數(shù)據(jù)源提取湖南省38 個(gè)柑橘重點(diǎn)縣柑橘果園的空間分布圖,引入數(shù)字高程模型(DEM) 作為地理輔助數(shù)據(jù),通過模型的高程、坡度、坡向等因子來描述研究區(qū)域的地形地貌特點(diǎn)。DEM 數(shù)據(jù)源于美國(guó)地質(zhì)調(diào)查局(https://glovis.usgs.gov/),空間分辨率為30 m。湖南省氣象數(shù)據(jù)(氣溫,降水量) 源于世界氣象數(shù)據(jù)庫(http://www.worldclim.org/),氣象數(shù)據(jù)空間分辨率1 km,用于提取湖南省近30 年的平均氣溫、年降水量數(shù)據(jù)。成土母質(zhì)與土壤類型數(shù)據(jù)源于第二次全國(guó)土壤普查1∶50000 的湖南省土壤圖(http://www.soilinfo.cn/map/index.aspx)。

    1.4 指標(biāo)分級(jí)標(biāo)準(zhǔn)

    柑橘園土壤pH 值分級(jí)標(biāo)準(zhǔn)根據(jù)2019 年農(nóng)業(yè)農(nóng)村部耕地質(zhì)量監(jiān)測(cè)保護(hù)中心印發(fā)的《湖南省耕地質(zhì)量監(jiān)測(cè)指標(biāo)分級(jí)標(biāo)準(zhǔn)》而定:pH≤4.5 (強(qiáng)酸性),4.5~5.5 (酸性),5.5~6.5 (弱酸性),6.5~7.5 (中性),7.5~8.5 (堿性),gt;8.5 (強(qiáng)堿性)。坡度分級(jí)標(biāo)準(zhǔn)依據(jù)1984 年中國(guó)農(nóng)業(yè)區(qū)劃委員會(huì)頒發(fā)的《土地利用現(xiàn)狀調(diào)查技術(shù)規(guī)程》而定:≤2° (一般無水土流失)、2°~6° (輕度土壤侵蝕)、6°~15° (易發(fā)生中度水土流失)、15°~25° (水土流失嚴(yán)重)、gt;25°(限制開荒)。坡向分級(jí)標(biāo)準(zhǔn)根據(jù)方位角劃分:135°~225° (陽坡)、225°~315° (半陽坡)、45°~135° (半陰坡)、0°~45°與315°~360° (陰坡)。將海拔高度按自然斷點(diǎn)法分為:0~85、85~170、170~255、255~340、340~425、425~490 m;年均氣溫劃分為:15.5℃~1 6℃、1 6℃~1 6 . 5℃、1 6 . 5℃~1 7℃、1 7℃~17.5℃、17.5℃~18℃、18℃~18.5℃;年均降水量劃分為:1180~1260、1260~1340、1340~1420、1420~1500、1500~1580、1580~1660 mm。養(yǎng)分分級(jí)標(biāo)準(zhǔn)參照魯劍巍[20]的研究方法(表2)。同時(shí)分別統(tǒng)計(jì)不同指標(biāo)對(duì)應(yīng)級(jí)別中土壤pH 平均值進(jìn)行數(shù)據(jù)分析與處理。

    1.5 結(jié)構(gòu)方程模型

    結(jié)構(gòu)方程模型(structural equation model,SEM),基于變量的協(xié)方差矩陣來分析變量之間關(guān)系的一種統(tǒng)計(jì)方法,是多元數(shù)據(jù)分析的重要工具。相關(guān)性分析得到的相關(guān)系數(shù)并不能表明因果關(guān)系,但是結(jié)構(gòu)方程模型卻能通過假設(shè)影響路徑,對(duì)復(fù)雜的多變量數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析,有效揭示變量之間的直接間接因果關(guān)系[21]。在SEM 模型中,總效應(yīng)包括直接效應(yīng)和間接效應(yīng),其中間接效應(yīng)的強(qiáng)度等于變量之間直接效應(yīng)的路徑系數(shù)的乘積,所有直接效應(yīng)與間接效應(yīng)的總和是總效應(yīng)強(qiáng)度。使用卡方檢驗(yàn)(χ2) 評(píng)估模型的擬合效果,當(dāng)Pgt;0.05,各參數(shù)中,χ2 /dflt;3,GFIgt;0.9,AGFIgt;0.9,RMSERlt;0.08 時(shí),SEM 擬合效果較為理想[22]。

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

    利用Excel 2016、SPSS 230 進(jìn)行數(shù)據(jù)統(tǒng)計(jì)與分析(ANOVA,Games-Howell 法),使用ArcGIS 10.6軟件進(jìn)行普通克里金插值法繪制土壤屬性空間分布圖,利用Amos 24 進(jìn)行結(jié)構(gòu)方程模型分析,采用Origin 2018 軟件作圖。在進(jìn)行數(shù)據(jù)分析之前,使用三倍標(biāo)準(zhǔn)差法(閾值法) 剔除異常值。

    2 結(jié)果與分析

    2.1 柑橘園土壤pH 值統(tǒng)計(jì)特征

    統(tǒng)計(jì)結(jié)果(圖2) 表明,湖南省柑橘園主產(chǎn)區(qū)土壤pH 值分布較為集中,服從偏正態(tài)分布,符合地統(tǒng)計(jì)學(xué)分析的要求。湖南省pH 值的變幅為3.00~8.50,均值為5.58,變異系數(shù)為18.06%,屬于中等程度的空間變異。各地區(qū)土壤pH 值除湘北洞庭湖區(qū)與湘西南區(qū)外差異顯著,平均值由大到小為:湘東湘中區(qū)(5.70)gt;湘西南區(qū)(5.61)gt;湘北洞庭湖區(qū)(5.52)gt;湘南區(qū)(5.50),均為弱酸性。其中,湘東湘中區(qū)與湘西南區(qū)土壤pH 均值顯著高于湖南省柑橘園土壤pH 均值,而湘北洞庭湖區(qū)與湘南區(qū)顯著低于全省pH 均值。從變異系數(shù)可知,各地區(qū)土壤pH 值具有中等變異水平(15.78%~19.22%),其中湘西南區(qū)變異系數(shù)高于其他地區(qū),說明該地區(qū)受到的干擾因素最為強(qiáng)烈。

    2.2 柑橘園土壤pH 值的空間分布

    由圖3 可知,湖南省柑橘園主產(chǎn)區(qū)大部分土壤呈酸性和弱酸性,面積共計(jì)189670 hm2,占比為93.3%,土壤普遍偏酸;其中,柑橘最適生長(zhǎng)在弱酸性土壤中,其分布面積為90191 hm2,面積占比為44.35%,主要分布于湘東中區(qū)、湘南區(qū)和湘西南區(qū),以湘西南分布面積最廣(34648 hm2),湘東中區(qū)面積占比最高(62.8%);酸性土壤分布面積為99479 hm2,面積占比為48.92%,重點(diǎn)分布于湘南區(qū)、湘西南區(qū)和湘北洞庭湖區(qū), 以湘西南的分布面積最大(33469 hm2 ),湘南與湘北洞庭湖區(qū)面積占比最高(分別占59.2% 和59.2%);不適宜柑橘生長(zhǎng)的強(qiáng)酸性土壤面積為3439 hm2,占總面積1.69%,主要分布在湘中東和湘西南區(qū), 湘西南區(qū)面積和占比最大(2952 hm2,3.73%);中性土壤的面積為10064 hm2,占總面積的4.95%,以湘西南區(qū)分布面積和占比最大(7990 hm2,10.08%);堿性土壤面積僅有167.8 hm2,占總面積的0.08%,全部分布于湘西南區(qū);各地區(qū)未有強(qiáng)堿性土壤分布。因此,柑橘園的堿性、中性土壤占比最少,有94.96% 的土壤為酸性土壤(土壤pH 值 lt;6.5),其中超過一半(50.6%) 的土壤處于過酸狀態(tài)(酸性和強(qiáng)酸性) 不利于柑橘生長(zhǎng)。同時(shí)Chen等[4]于2011—2019 采集361 個(gè)湖南省柑橘園土壤樣品,研究也發(fā)現(xiàn)79.0% 的臍橙園土壤pH 值(3.56~8.12) 和74.7% 的普通甜橙園土壤pH 值(2.75~8.15)發(fā)生了土壤酸化。由此說明研究區(qū)土壤酸化較為嚴(yán)重。

    2.3 湖南省柑橘園地形因子與氣候因子

    如圖4 所示,研究區(qū)地形以山地丘陵為主,海拔總體沿西北與東南向中部逐漸降低,海拔落差較大。而降水量由北到南逐步增加。

    由圖5 可見,湖南省柑橘園土壤pH 均值總體上隨海拔高度的上升呈先下降后升高的趨勢(shì)。而土壤pH 均值隨坡度的上升呈先下降后上升的勢(shì)態(tài),且不同坡度等級(jí)間柑橘園土壤pH 均值差異較大。不同坡向下土壤pH 均值無明顯變化規(guī)律,陰坡土壤pH 均值最高(5.37),半陽坡與半陰坡土壤pH 均值最低(5.34),二者相差僅0.03,說明坡向?qū)H 值的影響并不大。土壤pH 均值整體隨氣溫的升高呈降低趨勢(shì),隨降水量的增加呈升高趨勢(shì)。說明氣溫與降水量能一定程度影響土壤pH 均值。

    2.4 湖南省柑橘園成土母質(zhì)與土壤類型

    湖南省柑橘主產(chǎn)區(qū)土壤母質(zhì)共7 種類型,分別為第四紀(jì)紅色黏土、花崗巖風(fēng)化物、板頁巖風(fēng)化物、砂巖風(fēng)化物、石灰?guī)r風(fēng)化物、紫色砂頁巖風(fēng)化物、河湖沖(沉) 積物(圖6A);土壤類型主要為7 類,分別為紅壤、黃壤、紅色石灰土、黑色石灰土、紫色土、潮土和水稻土(圖6B)。不同成土母質(zhì)土壤pH 值范圍為5.11~5.60 (圖6C),不同土壤類型下土壤pH 值范圍為5.20~6.31 (圖6D)。其中,7 種成土母質(zhì)中pH 值差異顯著,依次為:石灰?guī)r風(fēng)化物(5.60)gt;河湖沖(沉) 積物(5.47)gt;紫色砂巖風(fēng)化物( 5 . 4 3 ) gt;板頁巖風(fēng)化物( 5 . 2 3 ) gt;第四紀(jì)紅色黏土(5.16)gt;砂巖風(fēng)化物(5.15)gt;花崗巖風(fēng)化物(5.11),僅第四紀(jì)紅色黏土的土壤pH 值與砂巖風(fēng)化物的土壤pH 值差異不顯著。其中,僅石灰?guī)r風(fēng)化物土壤屬于弱酸性土壤,其他母質(zhì)土壤均屬于酸性土壤。類似地,7 種土壤類型間pH 值差異顯著,依次為:黑色石灰土(6.31)gt;紅色石灰土(6.15)gt;潮土(5.49)gt;紫色土(5.43)gt;水稻土(5.27)gt;紅壤(5.25)gt;黃壤(5.20)。

    2.5 湖南省柑橘園土壤pH 值與影響因子相關(guān)性分析

    Pearson 相關(guān)性分析可以有效地反映土壤pH 值與影響因子之間的關(guān)系。如表3 所示,土壤pH 值與降水量呈顯著正相關(guān)關(guān)系(Plt;0.05),和速效鉀含量呈極顯著正相關(guān)關(guān)系(Plt;0.01);氣溫、堿解氮、成土母質(zhì)、土壤類型和有效磷與土壤pH 值呈極顯著負(fù)相關(guān)關(guān)系 (Plt;0.01)。而海拔高度、坡度和坡向與pH 值無顯著相關(guān)性。

    2.6 湖南省柑橘園土壤pH 值和影響因素的結(jié)構(gòu)方程模型

    進(jìn)一步揭示因子間相互作用,探索影響因子對(duì)土壤pH 值的直接和間接影響,我們使用了結(jié)構(gòu)方程模型分析(圖7)。模型的適配度檢驗(yàn)結(jié)果顯示:χ2 =0.324,df=2,P=0.850,GFI=1.000,AGFI=1.000,RMSEA=0.000,說明本研究中模型和數(shù)據(jù)的適配度較高、擬合情況理想,可以滿足研究和分析的要求。結(jié)果顯示,降水量、氣溫和海拔高度對(duì)土壤pH 值有直接和間接影響,土壤類型主要為直接影響,而成土母質(zhì)主要通過影響土壤類型進(jìn)而影響柑橘園土壤pH 值,其他因子對(duì)土壤pH 值的影響不顯著(圖7)。計(jì)算各因子對(duì)土壤pH 值的間接效應(yīng)和總效應(yīng),其中間接效應(yīng)的強(qiáng)度等于變量之間直接效應(yīng)的路徑系數(shù)的乘積,所有直接效應(yīng)與間接效應(yīng)的總和為總效應(yīng)強(qiáng)度。間接效應(yīng)結(jié)果顯示成土母質(zhì)、海拔高度、降水量與氣溫依次對(duì)土壤pH 值有極顯著效應(yīng),系數(shù)分別為?0.20、0.17、?0.10 和0.10 (Plt;0.01);總效應(yīng)結(jié)果表明土壤類型、成土母質(zhì)和氣溫分別對(duì)土壤pH 值有極顯著負(fù)效應(yīng),系數(shù)分別為?0.27、?0.20和?0.07 (Plt;0.01);降水量對(duì)土壤pH 值有極顯著正效應(yīng)(Plt;0.01)。而海拔高度對(duì)土壤pH 值的直接與間接效應(yīng)相互抵消,總效應(yīng)為0.00 (Plt;0.01),與相關(guān)性分析結(jié)果一致。因此土壤類型與成土母質(zhì)是本研究因子中解釋湖南省柑橘園土壤pH 值空間變異的主導(dǎo)因子。

    2.7 湖南省柑橘園土壤pH 值與速效養(yǎng)分含量變化趨勢(shì)分析

    柑橘園土壤養(yǎng)分狀況是制定果園土壤管理和施肥方案的重要依據(jù)之一,了解土壤pH 值對(duì)速效養(yǎng)分的影響對(duì)優(yōu)化施肥管理有重要作用。相關(guān)性分析結(jié)果表明,研究區(qū)堿解氮、有效磷含量與土壤pH 值存在極顯著負(fù)相關(guān)關(guān)系(?0.062**和?0.137**),速效鉀含量與土壤pH 值存在極顯著正相關(guān)關(guān)系(0.063**)。由圖8 可見,按柑橘土壤分級(jí)標(biāo)準(zhǔn)將土壤pH 值分為6 級(jí)(強(qiáng)酸性、酸性、弱酸性、中性、堿性和強(qiáng)堿性) 統(tǒng)計(jì)速效養(yǎng)分含量,并進(jìn)行不同土壤pH 值分級(jí)間速效養(yǎng)分含量的顯著性檢驗(yàn),顯示湖南省柑橘園土壤堿解氮與有效磷含量處于較低與適宜水平,并以適宜水平為主,而速效鉀含量皆處于適宜水平。其中,隨酸度等級(jí)降低,堿解氮與有效磷含量顯著降低,在強(qiáng)堿性時(shí)處于養(yǎng)分含量較低水平;而速效鉀含量整體呈現(xiàn)升高趨勢(shì)。

    3 討論

    本研究發(fā)現(xiàn)湖南省的柑橘主產(chǎn)區(qū)土壤pH 值范圍在3.00~8.50,均值為5.58,處于柑橘類作物生長(zhǎng)的最佳土壤pH 值范圍(5.5~6.5)[23],比較適宜柑橘生長(zhǎng)。具體表現(xiàn)為湘中湘東地區(qū)大約62.57% 的果園土壤為弱酸性,土壤酸堿度適宜柑橘生長(zhǎng)。其中,新寧縣的柑橘種植面積最大(52.27%),是我國(guó)四大出口臍橙生產(chǎn)基地之一[24],這充分說明了優(yōu)質(zhì)果品與土壤酸堿度是密不可分的。但超過一半(57.83%) 的采樣點(diǎn)以及50.61% 的面積土壤pH 值lt;5.5,處于酸性至強(qiáng)酸性環(huán)境中,主要分布在湘南區(qū)與湘西南區(qū)。盡管柑橘在酸性和弱酸性土壤中均能生長(zhǎng),但長(zhǎng)期在強(qiáng)酸性環(huán)境中,土壤容易發(fā)生板結(jié),透氣性變差,根系長(zhǎng)期處于缺氧狀態(tài),易引起樹體發(fā)育不良[25]。因此,湖南省柑橘主產(chǎn)區(qū)土壤pH 偏酸性的特點(diǎn),已成為制約柑橘產(chǎn)業(yè)發(fā)展的因素。探索湖南省柑橘園主產(chǎn)區(qū)土壤pH 值的影響因素,加強(qiáng)對(duì)果園土壤的酸化治理具有重要意義。量化區(qū)域土壤pH 值的空間異質(zhì)性并揭示異質(zhì)性形成的主要驅(qū)動(dòng)因子,可以有針對(duì)性地進(jìn)行土壤改良[26]。木合塔爾·艾買提等[27]基于山地丘陵交錯(cuò)地帶研究發(fā)現(xiàn),土壤pH 值在空間上存在一定的分布差異,受自然因素和人為因素的共同作用。本研究結(jié)果顯示各地區(qū)柑橘園土壤pH 值的變異系數(shù)為15.78%~19.22%,表明柑橘園土壤pH 值具有中等程度的空間變異性,因此湖南省柑橘園土壤pH 值空間變異是結(jié)構(gòu)性因素和隨機(jī)性因素綜合作用的結(jié)果。吳正祥等[28]研究發(fā)現(xiàn)成土母質(zhì)、海拔高度是十堰市耕層土壤pH 值空間變異的主控因素。在本研究結(jié)果中,不同土壤類型與成土母質(zhì)間土壤pH 值也具有顯著性差異;且結(jié)構(gòu)方程模型表明各影響因素中,土壤類型是土壤pH 值最重要的影響因素,其總效應(yīng)最強(qiáng)(0.27)。其次為成土母質(zhì),其總效應(yīng)為0.20。這主要是由于不同的土壤類型具有不同的成土條件和發(fā)育過程,從而影響土壤的物理結(jié)構(gòu)與鹽基物質(zhì)組成,導(dǎo)致土壤pH 值具有較大差異[29?30]。本研究中,黑色石灰土與紅色石灰土屬于適宜柑橘生長(zhǎng)的弱酸性土壤,是因?yàn)槭規(guī)r發(fā)育的土壤比較厚,土壤粘接性極強(qiáng),大部分為塊狀結(jié)構(gòu),土壤具有很強(qiáng)的緩沖酸化能力[8],主要集中在湘西南區(qū)的西部及中部地區(qū)。紫色土發(fā)育于紫色砂巖,在鹽基物質(zhì)的輕度淋失和母質(zhì)不斷風(fēng)化的快速補(bǔ)充作用下,土壤的鹽基飽和度較高,pH 值較高[28],在湘西南區(qū)中部和東部以及湘北洞庭湖區(qū)中部地區(qū)分布較多。而紅壤在高溫多雨條件下經(jīng)礦物加速風(fēng)化、鹽基大量淋失及脫硅富鋁化形成的發(fā)育程度較高的土壤,呈酸性[31],在湖南省柑橘園的面積高達(dá)61.68%,主要分布在湘東湘中區(qū)南部和西南部、湘北洞庭湖區(qū)中部和東北部、湘西南區(qū)西北部和東南部以及湘南區(qū)的大部分地區(qū);同時(shí),黃壤的發(fā)育程度較深,淋溶作用劇烈,土壤酸度高,呈酸性至強(qiáng)酸性反應(yīng)[29],但在湖南省柑橘園中分布面積較少(3.21%)。此外,水稻土的生境對(duì)土壤酸度有特殊的緩沖作用,具有一定的弱酸緩沖能力[32?33],其酸度顯著低于黃壤與紅壤。因此在柑橘園管理時(shí),針對(duì)土壤較酸的紅壤地區(qū)應(yīng)加強(qiáng)酸化治理。同樣地,土壤pH 值對(duì)母質(zhì)具有一定繼承性,不同母質(zhì)導(dǎo)致土壤pH 值呈現(xiàn)顯著性差異[34]。一般砂巖、板頁巖、花崗巖風(fēng)化過程比較完全,處于脫硅富鋁化階段時(shí)大量鹽基離子淋失,致使土壤呈酸性[35];而石灰?guī)r、紫色砂巖、河湖沖(沉) 積物多發(fā)育為堿性土壤,但卻在本研究中屬弱酸性或酸性,這與果園長(zhǎng)期單施化肥、柑橘自身根部分泌酸和酸沉降等因素有關(guān)[36?37]。2010—2018 年多年柚園調(diào)查結(jié)果發(fā)現(xiàn),過量施用化肥引起土壤鹽基離子的淋失,導(dǎo)致H+產(chǎn)生量急劇增加,造成土壤pH 值持續(xù)下降[38]。研究區(qū)屬亞熱帶季風(fēng)氣候,高溫高降雨,土壤侵蝕和養(yǎng)分流失嚴(yán)重[39?40]。相關(guān)結(jié)果也表明,湖南省柑橘園土壤pH 值與氣溫呈極顯著負(fù)相關(guān)關(guān)系。但也有研究發(fā)現(xiàn),坡度較小的地區(qū)土壤粘粒含量較高,膠體具有較大的比表面積,容易吸附陽離子,從而增加其酸緩沖能力[41];而坡度較大的地區(qū)土壤粘粒含量較低,土壤具有較弱的酸緩沖性質(zhì)[42]。而本研究發(fā)現(xiàn)坡向與坡度對(duì)柑橘園土壤pH 值影響都沒有顯著性差異,且與結(jié)構(gòu)方程模型結(jié)果一致。說明湖南省柑橘園土壤酸堿度主要受土壤自身因素的影響,而氣候和地形因素是次要因素。

    果園土壤養(yǎng)分對(duì)柑橘生長(zhǎng)發(fā)育、開花結(jié)果和果實(shí)品質(zhì)有著不可替代的作用[43]。曹勝等[2]采集84 個(gè)土壤樣品探究了湖南省柑橘園土壤pH 值和主要土壤養(yǎng)分特征的關(guān)系,發(fā)現(xiàn)土壤pH 值與土壤養(yǎng)分存在相關(guān)性,通過調(diào)節(jié)土壤酸堿度可以控制土壤養(yǎng)分的有效性。本研究發(fā)現(xiàn)湖南柑橘園土壤pH 值與有效態(tài)養(yǎng)分關(guān)系密切,土壤pH 值與堿解氮含量呈極顯著負(fù)相關(guān)關(guān)系,這可能是由于施入土壤中的氮素有14%~52% 以硝化反硝化、淋溶和氨揮發(fā)等形式損失,其中土壤為堿性時(shí)氮素以氨揮發(fā)形式損失,且總N 損失量占比高達(dá)9%~40%,堿解氮含量與土壤pH 值呈現(xiàn)負(fù)相關(guān)關(guān)系[44],土壤為酸性時(shí)果園長(zhǎng)期施用銨態(tài)氮肥,氮素化肥施用后在土壤中發(fā)生硝化作用產(chǎn)生NO3?,殘留H+引起土壤酸化[45];土壤pH 值與有效磷存在極顯著負(fù)相關(guān)關(guān)系,土壤磷在酸性條件下,容易與土壤中鐵錳氧化物形成磷酸鹽沉淀或閉蓄態(tài)的磷酸鹽導(dǎo)致有效性下降,在堿性條件下則容易形成磷酸鈣沉淀,導(dǎo)致磷的有效性下降。因此,一般認(rèn)為近中性條件下磷的有效性最高。但在柑橘的種植過程中,農(nóng)民為追求經(jīng)濟(jì)效益,在酸性土壤中為增加有效磷含量,過量施用磷肥的情況十分突出[46],造成柑橘園土壤磷素大量投入,進(jìn)而使得土壤中有效磷的含量提升,造成土壤酸性強(qiáng)時(shí)出現(xiàn)有效磷含量高這一結(jié)果,導(dǎo)致土壤pH 值與有效磷存在極顯著負(fù)相關(guān)關(guān)系,這一現(xiàn)象在果園土壤施肥中十分常見[22, 47]。而過量施用化肥雖然能增加土壤中的養(yǎng)分,但又會(huì)使土壤酸度增加,帶來的土壤酸化問題也會(huì)導(dǎo)致肥料利用率的下降。土壤pH 值與速效鉀含量呈極顯著正相關(guān)關(guān)系,這是因?yàn)殁涬x子在土壤中移動(dòng)性較強(qiáng),當(dāng)土壤過酸時(shí),鉀離子更容易淋溶損失;當(dāng)土壤pH 值增高時(shí),土壤含鹽基增加,鉀的有效性隨之增大[48]。說明湖南省柑橘園作為農(nóng)業(yè)土壤,受人為施肥因素干擾較大,土壤酸堿度限制著土壤養(yǎng)分的有效性,土壤肥料的施用又影響著土壤的酸堿性。本研究還發(fā)現(xiàn)各速效養(yǎng)分含量雖然在土壤pH 組間存在極顯著差異,但湖南省絕大部分柑橘園土壤,尤其是土壤pH 值為 5.6~6.3 弱酸性環(huán)境時(shí),各養(yǎng)分有效含量皆處于適宜水平,土壤堿解氮、有效磷、速效鉀含量適宜柑橘的生產(chǎn),但果樹對(duì)養(yǎng)分的積累不僅受土壤養(yǎng)分含量高低的影響,還與土壤對(duì)養(yǎng)分的涵養(yǎng)及樹體對(duì)養(yǎng)分的吸收能力有關(guān)[49?50]。研究發(fā)現(xiàn)施用有機(jī)肥一方面能調(diào)節(jié)土壤pH 值,緩解柑橘園土壤酸化[34];另一方面能提高土壤有機(jī)質(zhì)含量,增加土壤速效養(yǎng)分含量,涵養(yǎng)土壤養(yǎng)分[51]。因此,除不適宜柑橘種植的強(qiáng)堿性土壤外,對(duì)柑橘園調(diào)控土壤酸堿度進(jìn)行改良與培肥時(shí),不僅要考慮土壤的自身母質(zhì)基礎(chǔ)和土壤類型,還推薦適當(dāng)控施化肥,重視有機(jī)肥投入,并根據(jù)柑橘的需肥特性在特定的生長(zhǎng)階段補(bǔ)充部分養(yǎng)分,因地制宜地進(jìn)行果園土壤管理,以達(dá)到土壤養(yǎng)分涵養(yǎng)、土壤酸化改良、土壤質(zhì)量提升與柑橘健康生產(chǎn)的目的。

    4 結(jié)論

    湖南省柑橘園土壤pH 值普遍偏酸,主要受土壤類型與成土母質(zhì)的影響。除強(qiáng)堿性土壤外,整體上湖南省柑橘園土壤速效氮磷鉀養(yǎng)分處于適宜水平,適宜柑橘的生產(chǎn)。

    參 考 文 獻(xiàn):

    [ 1 ]喬進(jìn)超, 李旭瑋, 施蕊. 中國(guó)柑橘三大產(chǎn)區(qū)產(chǎn)業(yè)探究: 以廣西、湖南、湖北為例[J]. 云南科技管理, 2021, 34(1): 45?49.

    Qiao J C, Li X W, Shi R. Industrial exploration of the three majorcitrus producing areas in China: Taking Guangxi, Hunan and Hubeias an example[J]. Yunnan Science and Technology Management,2021, 34(1): 45?49.

    [ 2 ]曹勝, 歐陽夢(mèng)云, 周衛(wèi)軍, 等. 湖南省柑橘園土壤pH和主要養(yǎng)分特征及其相互關(guān)系[J]. 中國(guó)土壤與肥料, 2020, (1): 31?38.

    Cao S, Ouyang M Y, Zhou W J, et al. Soil pH and main nutrientcharacteristics of citrus orchards and their correlation in HunanProvince[J]. Soil and Fertilizer Sciences in China, 2020, (1): 31?38.

    [ 3 ]Long A, Zhang J, Yang L T, et al. Effects of low pH on photosynthesis,related physiological parameters, and nutrient profiles of citrus[J].Frontiers in Plant Science, 2017, 8: 185.

    [ 4 ]洪勵(lì)偉, 劉夢(mèng)夢(mèng), 李娟, 等. 紅壤砂糖橘果園土壤和葉片養(yǎng)分狀況研究[J]. 土壤科學(xué), 2017, 5(4): 44?52.

    Hong L W, Liu M M, Li J, et al. Study on nutrient status of soil andleaves in “Shatangju” orchards in red earth region[J]. Hans Journalof Soil Science, 2017, 5(4): 44?52.

    [ 5 ]Li Y, Han M Q, Lin F, et al. Soil chemical properties, 'Guanximiyou'pummelo leaf mineral nutrient status and fruit quality in the southernregion of Fujian Province, China[J]. Journal of Soil Science and PlantNutrition, 2015, 15(3): 615?628.

    [ 6 ]武松偉, 梁珊珊, 胡承孝, 等. 我國(guó)柑橘園“因土補(bǔ)肥”與化肥減施增效生態(tài)分區(qū)[J]. 華中農(nóng)業(yè)大學(xué)學(xué)報(bào), 2022, 41(2): 9?19.

    Wu S W, Liang S S, Hu C X, et al. Ecological region division of soil based supplementary fertilization and decrement fertilization in Chinacitrus orchards[J]. Journal of Huazhong Agricultural University,2022, 41(2): 9?19.

    [ 7 ]Chen Y W, Li F F, Wu Y C, et al. Profiles of citrus orchard nutritionand fruit quality in Hunan Province, China[J]. International Journal ofFruit Science, 2022, 22(1): 779?793.

    [ 8 ]張?jiān)啵?吳穎, 鄭雄偉. 湖北省土壤酸堿度趨勢(shì)分析及影響因素[J].資源環(huán)境與工程, 2018, 32(S): 30?34.

    Zhang Y P, Wu Y, Zheng X W. Trend analysis and influencingfactors of soil acidity and alkalinity in Hubei[J]. Resources Environmentamp; Engineering, 2018, 32(S): 30?34.

    [ 9 ]Habib-ur-rahman M, Raza A, Ahrends H E, et al. Impact of in-fieldsoil heterogeneity on biomass and yield of winter triticale in anintensively cropped hummocky landscape under temperate climateconditions[J]. Precision Agriculture, 2022, 23(3): 912?938.

    [10]Dindaroglu T, Tunguz V, Babur E, et al. Determination of therelationship among compound Topographic Index (CTI), soilproperties and land-use in Karst ecosystems[J]. Physical Geography,2023, 44(3): 307?329.

    [11]劉雪松, 張智印, 魏建朋. 贛南地區(qū)土壤pH值分布的影響因素及其強(qiáng)度變化[J]. 水土保持通報(bào), 2021, 41(4): 100?105.

    Liu X S, Zhang Z Y, Wei J P. Influencing factors and intensitychange of soil pH value in south Jiangxi Province[J]. Bulletin of Soiland Water Conservation, 2020, 41(4): 95?100.

    [12]Srivastava A K, Singh S. Citrus decline: Soil fertility and plantnutrition[J]. Journal of Plant Nutrition, 2009, 32(2): 197?245.

    [13]Wang P, Wang Y, Shu B, et al. Relationships between arbuscularmycorrhizal symbiosis and soil fertility factors in citrus orchardsalong an altitudinal gradient[J]. Pedosphere, 2015, 25(1): 160?168.

    [14]Mattos J D, Kadyampakeni D M, Oliver A Q, et al. Soil and nutritioninteractions[M]. The genus citrus. London : Taylor amp; Francis, 2020.

    [15]Zekri M, Obreza T A, Koo R. Irrigation, nutrition, and citrus fruitquality[J]. Institute of Food and Agricultural Sciences, 2003, 207:1?3.

    [16]Pettigrew W T. Potassium influences on yield and quality productionfor maize, wheat, soybean and cotton[J]. Physiologia Plantarum,2008, 133(4): 670?681.

    [17]劉東海, 鄭守會(huì), 陳云峰, 等. 丹江口柑橘園土壤和葉片營(yíng)養(yǎng)狀況調(diào)查分析[J]. 中國(guó)土壤與肥料, 2020, (4): 95?100.

    Liu D H, Zheng S H, Chen Y F, et al. Studies on nutrient contents ofcitrus orchard soils and leaves in Danjiangkou[J]. Soil and FertilizerSciences in China, 2020, (4): 95?100.

    [18]唐志文, 蔣雅琴, 譚善凱, 李桂華. 新寧縣臍橙園土壤養(yǎng)分特征及其培肥改良對(duì)策[J]. 湖南農(nóng)業(yè)科學(xué), 2020, (7): 41?46.

    Tang Z W, Jiang Y Q, Tan S K, Li G H. Soil nutrient characteristicsof navel orange orchards in Xinning County and countermeasures ofimproving fertility[J]. Hunan Agricultural Sciences, 2020, (7): 41?46.

    [19]鮑士旦. 土壤農(nóng)化分析(3版)[M]. 中國(guó)農(nóng)業(yè)出版社, 2000.

    Bao S D. Soil and agricultural chemistry analysis (3rd edn)[M].China Agricultural Press, 2000.

    [20]魯劍巍. 湖北省柑橘園土壤—植物養(yǎng)分狀況與柑橘平衡施肥技術(shù)研究[D]. 湖北武漢: 華中農(nóng)業(yè)大學(xué)博士學(xué)位論文, 2003.

    Lu J W. Study on soil and plant nutrition status and balancedfertilization techniques of the citrus orchards in Hubei[D]. Wuhan,Hubei: PhD Dissertation of Huazhong Agricultural University, 2003.

    [21]李慧, 汪景寬, 裴久渤, 李雙異. 基于結(jié)構(gòu)方程模型的東北地區(qū)主要旱田土壤有機(jī)碳平衡關(guān)系研究[J]. 生態(tài)學(xué)報(bào), 2015, 35(2):517?525.

    Li H, Wang J K, Pei J B, Li S Y. Equilibrium relationships of soilorganic carbon in the main croplands of Northeast China based onstructural equation modeling[J]. Acta Ecologica Sinica, 2015, 35(2):517?525.

    [22]彭健健, 徐堅(jiān), 王曉曉, 等. 楊梅主產(chǎn)區(qū)土壤肥力空間異質(zhì)性及其影響因素: 以浙江仙居和臨海為例[J]. 果樹學(xué)報(bào), 2023, 40(7):1421?1433.

    Peng J J, Xu J, Wang X X, et al. Spatial variation of soil fertility andits influencing factors in Myrica Rubra region: A case study in XianjuCounty and Linhai City[J]. Journal of Fruit Science, 2023, 40(7):1421?1433.

    [23]溫明霞, 石孝均, 聶振朋, 周鑫斌. 重慶市柑桔園土壤酸堿度及金屬元素含量的變化特征[J]. 水土保持學(xué)報(bào), 2011, 25(5): 191?194.

    Wen M X, Shi X J, Nie Z P, Zhou X B. Characteristics of soil pH andmetal element contents of citrus orchard in Chongqing[J]. Journal ofSoil and Water Conservation, 2011, 25(5): 191?194.

    [24]肖雨辰. 新寧縣石灣村臍橙特色產(chǎn)業(yè)發(fā)展報(bào)告[J]. 農(nóng)村經(jīng)濟(jì)與科技, 2019, 30(20): 167?168.

    Xiao Y C. Development report of navel orange characteristic industryin Shiwan Village, Xinning County[J]. Rural Economy and Science-Technology, 2019, 30(20): 167?168.

    [25]溫明霞, 鹿連明, 王鵬, 等. 浙江省柑橘園土壤養(yǎng)分調(diào)查[J]. 浙江農(nóng)業(yè)科學(xué), 2019, 60(2): 208?211.

    Wen M X, Lu L M, Wang P, et al. Status of soil nutrients in citrusorchards of Zhejiang Province[J]. Journal of Zhejiang AgriculturalSciences, 2019, 60(2): 208?211.

    [26]景海超, 劉穎慧, 賀佩, 等. 青藏高原典型區(qū)生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性及其影響因素: 以那曲市為例[J]. 生態(tài)學(xué)報(bào), 2022, 42(7): 2657?2673.

    Jing H C, Liu Y H, He P, et al. Spatial heterogeneity of ecosystemservices and it's influencing factors in typical areas of the Oinghai-Tibet Plateau: A case study of Nagqu City[J]. Acta Ecologica Sinica,2022, 42(7): 2657?2673.

    [27]木合塔爾·艾買提, 肖鵬南, 周勇, 等. 山地丘陵交錯(cuò)地帶耕層土壤pH和養(yǎng)分空間變異特征及影響因素分析[J]. 南方農(nóng)業(yè)學(xué)報(bào), 2019,50(7): 1432?1441.

    Muhtar A, Xiao P N, Zhou Y, et al. Spatial variability of croplandsoil pH and nutrients and their affecting factors in mountainous hillyinterlaced zone[J]. Journal of Southern Agriculture, 2019, 50(7):1432?1441.

    [28]吳正祥, 周勇, 木合塔爾·艾買提, 李晴. 鄂西北山區(qū)耕層土壤pH值空間變異特征及其影響因素研究[J]. 長(zhǎng)江流域資源與環(huán)境, 2020,29(2): 488?498.

    Wu Z X, Zhou Y, Muhtar A, Li Q. Spatial variability of soil pH valueand its influencing factors in the soil layer of northwestern HubeiProvince[J]. Resources and Environment in the Yangtze Basin, 2020,29(2): 488?498.

    [29]鄧茜, 方紅艷, 張?jiān)拢?等. 四川盆地耕地土壤pH空間變異特征及影響因素[J]. 土壤, 2022, 54(6): 1283?1290.

    Deng Q, Fang H Y, Zhang Y Y, et al. Spatial variability of croplandsoil pH and its influencing factors in Sichuan Basin[J]. Soils, 2022,54(6): 1283?1290.

    [30]熊杏, 熊清華, 郭熙, 等. 南方典型丘陵區(qū)耕地土壤全氮、有機(jī)碳和碳氮比空間變異特征及其影響因素[J]. 植物營(yíng)養(yǎng)與肥料學(xué)報(bào),2020, 26(9): 1656?1668.

    Xiong X, Xiong Q H, Guo X, et al. Spatial variation characteristics oftotal nitrogen, organic carbon and ratio of carbon to nitrogen ofcultivated land in typical hilly areas in south China and its influencingfactors[J]. Journal of Plant Nutrition and Fertilizers, 2020, 26(9):1656?1668.

    [31]焦天. 現(xiàn)代土壤空間變異性及成壤影響因素研究: 以桑干河陽原段流域?yàn)槔齕D]. 北京: 中國(guó)地質(zhì)科學(xué)院碩士學(xué)位論文, 2020.

    Jiao T. Research on the spatial variability of modern soil and theinfluencing factors of soil formation: Taking the Yangyuan Basin ofSanggan River as an example[D]. Beijing: MS Thesis of ChineseAcademy of Geological Sciences, 2020.

    [32]Wang Z N, Jia M Y, Li Z, et al. Acid buffering capacity of fourcontrasting metal-contaminated calcareous soil types: Changes in soilmetals and relevance to phytoextraction[J]. Chemosphere, 2020, 256:127045.

    [33]Wei H, Liu Y L, Xiang H M, et al. Soil pH responses to simulatedacid rain leaching in three agricultural soils[J]. Sustainability, 2020,12(1): 280.

    [34]王玨, 周衛(wèi)軍, 商貴鐸, 等. 不同母質(zhì)柑橘園土壤養(yǎng)分特征及肥力綜合評(píng)價(jià)[J]. 生態(tài)學(xué)雜志, 2022, 41(5): 933?940.

    Wang J, Zhou W J, Shang G D, et al. Comprehensive evaluation ofsoil nutrient characteristics and fertility of citrus orchard withdifferent parent materials[J]. Chinese Journal of Ecology, 2022,41(5): 933?940.

    [35]趙凱麗, 王伯仁, 徐明崗, 等. 我國(guó)南方不同母質(zhì)土壤pH剖面特征及酸化因素分析[J]. 植物營(yíng)養(yǎng)與肥料學(xué)報(bào), 2019, 25(8): 1308?1315.

    Zhao K L, Wang B R, Xu M G, et al. Changes in pH with depths ofsoils derived from different parent materials and analysis ofacidification in Southern China[J]. Journal of Plant Nutrition andFertilizers, 2019, 25(8): 1308?1315.

    [36]張超博, 邱潔雅, 王敏, 等. 桂北柑橘園土壤化學(xué)性狀研究[J]. 土壤,2020, 52(6): 1187?1195.

    Zhang C B, Qiu J Y, Wang M, et al. Study on soil chemicalproperties of citrus orchards in northern Guangxi[J]. Soils, 2020,52(6): 1187?1195.

    [37]Dong Y, Yang J L, Zhao X R, et al. Soil acidification and loss of basecations in a subtropical agricultural watershed[J]. Science of the TotalEnvironment, 2022, 827: 154338.

    [38]Chen X H, Yu W H, Cai Y Y, et al. How to identify and adoptcleaner strategies to improve the continuous acidification in orchardsoils?[J]. Journal of Cleaner Production, 2022, 330: 129826.

    [39]Fei K, Deng L Z, Sun T Y, et al. Runoff processes and lateraltransport of soil total carbon induced by water erosion in thehilly region of southern China under rainstorm conditions[J].Geomorphology, 2019, 340: 143?152.

    [40]Li Z W, Ning K, Chen J, et al. Soil and water conservation effectsdriven by the implementation of ecological restoration projects:Evidence from the red soil hilly region of China in the last threedecades[J]. Journal of Cleaner Production, 2020, 260(1/3): 121109.

    [41]林躍勝, 馬康, 周浩, 等. 基于地貌單元的安徽省耕地土壤pH空間變異及其驅(qū)動(dòng)因子分析[J]. 環(huán)境科學(xué)學(xué)報(bào), 2023, 43(7): 318?330.

    Lin Y S, Ma K, Zhou H, et al. Spatial variation and driving factors ofsoil pH in cultivated land of Anhui Province based on geomorphicunit[J]. Acta Scientiae Circumstantiae, 2023, 43(7): 318?330.

    [42]Yan P, Peng H, Yan L B, et al. Spatial variability in soil pH and landuse as the main influential factor in the red beds of the NanxiongBasin, China[J]. Peerj, 2019, 7: e6342.

    [43]鄭朝霞, 王穎, 石磊, 等. 陜西省蘋果主產(chǎn)區(qū)土壤有機(jī)質(zhì)、氮磷鉀養(yǎng)分含量與分布特征[J]. 植物營(yíng)養(yǎng)與肥料學(xué)報(bào), 2017, 23(5):1191?1198.

    Zheng Z X, Wang Y, Shi L, et al. Contents and distribution of soilorganic matter and nitrogen, phosphate, potassium in the main appleproduction regions of Shaanxi Province[J]. Journal of Plant Nutritionand Fertilizers, 2017, 23(5): 1191?1198.

    [44]康興蓉, 周旋, 彭建偉, 等. 包膜尿素與普通尿素配施對(duì)稻田氨揮發(fā)的影響[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2022, 41(11): 2539?2551.

    Kang X R, Zhou X, Peng J W, et al. Effects of combined applicationof coated ureas and ordinary ureas on ammonia volatilization inpaddy fields[J]. Journal of Agro-Environment Science, 2022, 41(11):2539?2551.

    [45]李娟, 張立成, 章明清, 等. 長(zhǎng)期施用尿素降低赤紅壤旱地耕層pH的特征與預(yù)測(cè)[J]. 植物營(yíng)養(yǎng)與肥料學(xué)報(bào), 2022, 28(12): 2161?2171.

    Li J, Zhang L C, Zhang M Q, et al. Characteristics and prediction oftopsoil acidification of latosolic red soil caused by long-term ureaapplication[J]. Journal of Plant Nutrition and Fertilizers, 2022, 28(12):2161?2171.

    [46]何悅, 漆雁斌. 農(nóng)戶過量施肥風(fēng)險(xiǎn)認(rèn)知及環(huán)境友好型技術(shù)采納行為的影響因素分析: 基于四川省380個(gè)柑橘種植戶的調(diào)查[J]. 中國(guó)農(nóng)業(yè)資源與區(qū)劃, 2020, 41(5): 8?15.

    He Y, Qi Y B. Analysis on the risk cognition of excess fertilizerapplication and the adoption behavior of environment-friendlytechnology and its reason: Base on the survey of 380 citrus grower inSichuan Province[J]. Chinese Journal of Agricultural Resources andRegional Planning, 2020, 41(5): 8?15.

    [47]王雪. 安徽省碭山縣黃桃根系土pH、有機(jī)質(zhì)及有效態(tài)元素相關(guān)性分析[J]. 四川地質(zhì)學(xué)報(bào), 2022, 42(2): 270?274.

    Wang X. Correlation analysis of pH, organic matter and effectivestate elements of yellow peach root system soil in Dangshan,Anhui[J]. Acta Geologica Sichuan, 2022, 42(2): 270?274.

    [48]李玲, 張少凱, 吳克寧, 等. 基于土壤系統(tǒng)分類的河南省土壤有機(jī)質(zhì)時(shí)空變異[J]. 土壤學(xué)報(bào), 2015, 52(5): 979?990.

    Li L, Zhang S K, Wu K N, et al. Analysis on spatio-temporalvariability of soil organic matter in Henan Province based on soiltaxonomy[J]. Acta Pedologica Sinica, 2015, 52(5): 979?990.

    [49]王慶萱, 阮云澤, 趙鵬飛. 基于因子?聚類分析的火龍果園土壤養(yǎng)分狀況評(píng)價(jià) [J]. 熱帶生物學(xué)報(bào), 2024, 15(2): 150?156.

    Wang Q X, Ruan Y Z, Zhao P F. Soil nutrient evaluation of pitayaorchards based on factor and cluster analysis[J]. Journal of TropicalBiolog, 2024, 15(2): 150?156.

    [50]張小卓, 秦太峰, 張乃明. 云南彌勒葡萄園土壤肥力特征與評(píng)價(jià)[J].西南農(nóng)業(yè)學(xué)報(bào), 2013, 26(6): 2380?2385.

    Zhang X Z, Qin T F, Zhang N M. Characteristics and evaluation ofvineyard soil fertility in Mile County, Yunnan[J]. Southwest ChinaJournal of Agricultural Sciences, 2013, 26(6): 2380?2385.

    [51]魏勇, 王帥, 趙敬坤, 等. 重慶花椒主產(chǎn)區(qū)土壤養(yǎng)分分布特征及綜合肥力評(píng)價(jià)[J]. 西南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2022, 44(7): 48?58.

    Wei Y, Wang S, Zhao J K, et al. Characteristics of nutrientdistribution and comprehensive evaluation of fertility of Zanthoxylumbungeanum orchards soil in Chongqing[J]. Journal of SouthwestUniversity (Natural Science Edition), 2022, 44(7): 48?58.

    基金項(xiàng)目:國(guó)家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項(xiàng) (CARS-26)。

    搡老乐熟女国产| 中文字幕人妻熟女乱码| 久久精品人人爽人人爽视色| 人人妻人人澡人人爽人人夜夜| 99热全是精品| 最新在线观看一区二区三区 | 亚洲三区欧美一区| 亚洲综合色网址| 亚洲在久久综合| 国产欧美日韩一区二区三区在线| 国产精品一区二区精品视频观看| 成人黄色视频免费在线看| 人体艺术视频欧美日本| 精品亚洲乱码少妇综合久久| 久久99热这里只频精品6学生| 久久这里只有精品19| 1024香蕉在线观看| 国产乱人偷精品视频| 国产日韩欧美视频二区| 亚洲国产最新在线播放| 两个人看的免费小视频| 黑人欧美特级aaaaaa片| 国产黄色免费在线视频| 黄色视频不卡| 精品国产一区二区三区四区第35| 亚洲国产精品成人久久小说| 91精品伊人久久大香线蕉| 哪个播放器可以免费观看大片| 精品国产露脸久久av麻豆| 精品一品国产午夜福利视频| 91成人精品电影| www.自偷自拍.com| 超碰成人久久| 狂野欧美激情性bbbbbb| 韩国精品一区二区三区| 天天躁日日躁夜夜躁夜夜| 国产麻豆69| 日韩视频在线欧美| 黄色怎么调成土黄色| av免费观看日本| 1024香蕉在线观看| 欧美日韩亚洲综合一区二区三区_| 电影成人av| 国产成人欧美在线观看 | 天天躁夜夜躁狠狠久久av| 成人手机av| 熟女少妇亚洲综合色aaa.| 精品少妇内射三级| 热re99久久精品国产66热6| 免费少妇av软件| 国产在线免费精品| 18在线观看网站| 日韩大片免费观看网站| 在线观看国产h片| 中文字幕av电影在线播放| 欧美精品亚洲一区二区| 久久久久久人人人人人| 一级黄片播放器| 一二三四中文在线观看免费高清| 咕卡用的链子| 嫩草影院入口| 最近的中文字幕免费完整| 中文精品一卡2卡3卡4更新| 人体艺术视频欧美日本| 亚洲国产看品久久| 精品少妇内射三级| 香蕉国产在线看| 最新在线观看一区二区三区 | 麻豆精品久久久久久蜜桃| 亚洲精品国产区一区二| 高清在线视频一区二区三区| 日韩欧美一区视频在线观看| 国产在线视频一区二区| 大话2 男鬼变身卡| bbb黄色大片| 免费av中文字幕在线| 国精品久久久久久国模美| 桃花免费在线播放| 欧美最新免费一区二区三区| 美女国产高潮福利片在线看| 日本午夜av视频| av有码第一页| 18禁动态无遮挡网站| 在现免费观看毛片| 9热在线视频观看99| 欧美成人精品欧美一级黄| 午夜福利视频在线观看免费| 日韩,欧美,国产一区二区三区| 中文字幕人妻丝袜制服| 国产伦人伦偷精品视频| 亚洲天堂av无毛| 亚洲av成人精品一二三区| 日韩av不卡免费在线播放| 久久久久精品久久久久真实原创| 亚洲av电影在线观看一区二区三区| 街头女战士在线观看网站| 国产精品无大码| 在线精品无人区一区二区三| 国产欧美亚洲国产| 不卡视频在线观看欧美| 最近的中文字幕免费完整| 看免费av毛片| 少妇人妻 视频| 肉色欧美久久久久久久蜜桃| 美女午夜性视频免费| 男女床上黄色一级片免费看| 18禁观看日本| 亚洲精品国产av成人精品| 亚洲精品成人av观看孕妇| 国产亚洲av片在线观看秒播厂| 亚洲,欧美,日韩| 国产成人精品久久二区二区91 | 国产成人av激情在线播放| av女优亚洲男人天堂| 老熟女久久久| a级毛片黄视频| 国产精品一区二区在线不卡| 在线看a的网站| 午夜福利视频在线观看免费| 亚洲精品成人av观看孕妇| 免费观看性生交大片5| 亚洲av福利一区| 亚洲一级一片aⅴ在线观看| 亚洲情色 制服丝袜| 91精品伊人久久大香线蕉| 最近手机中文字幕大全| 侵犯人妻中文字幕一二三四区| 一级黄片播放器| 韩国av在线不卡| 777久久人妻少妇嫩草av网站| 国产亚洲午夜精品一区二区久久| 久久精品国产亚洲av高清一级| 久久狼人影院| 免费黄色在线免费观看| 啦啦啦啦在线视频资源| 三上悠亚av全集在线观看| 国产日韩欧美亚洲二区| 欧美人与善性xxx| 成人毛片60女人毛片免费| 波多野结衣一区麻豆| 一级片'在线观看视频| 中国国产av一级| 在线观看免费高清a一片| 老司机亚洲免费影院| 亚洲欧美一区二区三区国产| 久久精品久久久久久久性| 99热全是精品| 国产亚洲精品第一综合不卡| 亚洲国产精品999| 国产精品香港三级国产av潘金莲 | 午夜福利,免费看| 欧美精品一区二区大全| 亚洲伊人久久精品综合| 黑人欧美特级aaaaaa片| 2018国产大陆天天弄谢| 成人国产av品久久久| 99热国产这里只有精品6| 又大又爽又粗| 满18在线观看网站| 亚洲人成网站在线观看播放| 日本欧美视频一区| 亚洲在久久综合| 涩涩av久久男人的天堂| 亚洲精华国产精华液的使用体验| 人人妻人人爽人人添夜夜欢视频| 国产免费福利视频在线观看| 国产精品成人在线| 日本午夜av视频| 精品久久蜜臀av无| 色视频在线一区二区三区| 无限看片的www在线观看| 国产亚洲一区二区精品| 午夜福利乱码中文字幕| 人人妻人人澡人人爽人人夜夜| 一级毛片黄色毛片免费观看视频| 成人亚洲欧美一区二区av| av在线观看视频网站免费| 欧美国产精品一级二级三级| 亚洲欧美一区二区三区久久| 亚洲专区中文字幕在线 | 国产精品成人在线| 亚洲综合色网址| 国产99久久九九免费精品| 国产免费一区二区三区四区乱码| 国产无遮挡羞羞视频在线观看| 欧美黑人精品巨大| 国产伦人伦偷精品视频| 看非洲黑人一级黄片| 精品第一国产精品| www日本在线高清视频| 国产爽快片一区二区三区| 中文字幕人妻丝袜一区二区 | 曰老女人黄片| 精品人妻一区二区三区麻豆| 亚洲成av片中文字幕在线观看| 婷婷色麻豆天堂久久| svipshipincom国产片| 99久国产av精品国产电影| 亚洲第一av免费看| 亚洲美女黄色视频免费看| 丁香六月欧美| 两个人免费观看高清视频| 国产精品无大码| 咕卡用的链子| 高清视频免费观看一区二区| 亚洲免费av在线视频| 老司机在亚洲福利影院| 波野结衣二区三区在线| 精品一区在线观看国产| videos熟女内射| 日韩一本色道免费dvd| 综合色丁香网| 精品视频人人做人人爽| 男人操女人黄网站| 下体分泌物呈黄色| 夫妻性生交免费视频一级片| 久久久久久久久免费视频了| 国产精品香港三级国产av潘金莲 | av国产精品久久久久影院| tube8黄色片| 亚洲精品,欧美精品| 毛片一级片免费看久久久久| 欧美黑人欧美精品刺激| av网站在线播放免费| 超碰97精品在线观看| 久久国产精品大桥未久av| 精品一区在线观看国产| 国产精品久久久av美女十八| av福利片在线| 日韩中文字幕欧美一区二区 | 爱豆传媒免费全集在线观看| 18禁动态无遮挡网站| 男女下面插进去视频免费观看| 亚洲精品久久午夜乱码| 国产精品国产av在线观看| 日本wwww免费看| 麻豆精品久久久久久蜜桃| 国产熟女欧美一区二区| 亚洲精品av麻豆狂野| 美女午夜性视频免费| 99香蕉大伊视频| 欧美日韩精品网址| 精品国产一区二区三区久久久樱花| 制服人妻中文乱码| 精品一品国产午夜福利视频| 国产97色在线日韩免费| 中文字幕亚洲精品专区| 亚洲一区二区三区欧美精品| 99久国产av精品国产电影| 捣出白浆h1v1| 黑人欧美特级aaaaaa片| 精品人妻在线不人妻| 成人国产麻豆网| 日韩欧美精品免费久久| 国产深夜福利视频在线观看| 热99久久久久精品小说推荐| 免费观看人在逋| 人体艺术视频欧美日本| 一级片免费观看大全| 性少妇av在线| 在线观看国产h片| 亚洲精品在线美女| 免费黄频网站在线观看国产| 成人免费观看视频高清| 亚洲色图 男人天堂 中文字幕| 欧美久久黑人一区二区| 少妇人妻精品综合一区二区| 久久亚洲国产成人精品v| 两性夫妻黄色片| 日韩免费高清中文字幕av| 久久精品熟女亚洲av麻豆精品| 一区二区三区精品91| 国产成人免费无遮挡视频| 国产一区有黄有色的免费视频| 美女大奶头黄色视频| 天堂8中文在线网| 免费在线观看完整版高清| 国产精品成人在线| 亚洲国产欧美在线一区| 天堂俺去俺来也www色官网| 欧美日本中文国产一区发布| 日本vs欧美在线观看视频| 成年人午夜在线观看视频| 久久久久久人人人人人| 久久精品国产亚洲av涩爱| 国产一区二区在线观看av| 午夜福利影视在线免费观看| 国产精品久久久人人做人人爽| 香蕉国产在线看| 一本久久精品| 宅男免费午夜| 99精品久久久久人妻精品| 9191精品国产免费久久| 丝袜喷水一区| 在线观看免费日韩欧美大片| 久久久久久久精品精品| 精品亚洲成国产av| 一级片免费观看大全| 久久97久久精品| 久久人人97超碰香蕉20202| www.av在线官网国产| 日本午夜av视频| 国产亚洲最大av| 999久久久国产精品视频| 欧美日韩亚洲综合一区二区三区_| 少妇的丰满在线观看| 少妇人妻精品综合一区二区| 欧美黑人精品巨大| 中文欧美无线码| 伦理电影免费视频| 亚洲欧美中文字幕日韩二区| 在线观看免费视频网站a站| 一级片免费观看大全| 成年女人毛片免费观看观看9 | 久久99精品国语久久久| 精品少妇久久久久久888优播| 亚洲国产成人一精品久久久| 亚洲欧美精品自产自拍| 国产高清国产精品国产三级| 中文字幕制服av| 在线天堂最新版资源| 高清av免费在线| 一边摸一边做爽爽视频免费| 国产毛片在线视频| 国产精品一国产av| 大片免费播放器 马上看| 亚洲美女搞黄在线观看| 欧美最新免费一区二区三区| 日日撸夜夜添| 五月天丁香电影| 国产av一区二区精品久久| 人人妻人人澡人人看| 水蜜桃什么品种好| 老司机影院毛片| 免费av中文字幕在线| 乱人伦中国视频| 亚洲av国产av综合av卡| 日日摸夜夜添夜夜爱| 一区二区日韩欧美中文字幕| 欧美日韩视频精品一区| 亚洲成人手机| 日韩精品免费视频一区二区三区| 欧美人与性动交α欧美精品济南到| 亚洲欧美清纯卡通| 亚洲精华国产精华液的使用体验| 日韩精品有码人妻一区| 亚洲精品日本国产第一区| 久久久欧美国产精品| 在线观看免费午夜福利视频| 亚洲四区av| 久久精品国产亚洲av涩爱| 少妇人妻 视频| 在现免费观看毛片| 午夜福利乱码中文字幕| av福利片在线| 婷婷色综合www| 亚洲在久久综合| av网站在线播放免费| 老汉色∧v一级毛片| videosex国产| 一级片'在线观看视频| 999精品在线视频| 叶爱在线成人免费视频播放| 少妇的丰满在线观看| 性色av一级| 免费观看人在逋| 久久久久国产一级毛片高清牌| 美女国产高潮福利片在线看| 亚洲av日韩精品久久久久久密 | 在线天堂最新版资源| 一级片'在线观看视频| 黄网站色视频无遮挡免费观看| 亚洲成色77777| 日本一区二区免费在线视频| 99久久99久久久精品蜜桃| videos熟女内射| 亚洲欧美精品自产自拍| 久久久久国产一级毛片高清牌| 国产片特级美女逼逼视频| 日韩 欧美 亚洲 中文字幕| 超色免费av| 色视频在线一区二区三区| 亚洲精品一区蜜桃| 日韩大码丰满熟妇| 丝袜美足系列| 国产精品麻豆人妻色哟哟久久| 午夜福利网站1000一区二区三区| 丰满乱子伦码专区| 麻豆精品久久久久久蜜桃| 美女脱内裤让男人舔精品视频| 天美传媒精品一区二区| 无遮挡黄片免费观看| 91aial.com中文字幕在线观看| 99久国产av精品国产电影| 国产欧美亚洲国产| 一边亲一边摸免费视频| 欧美成人精品欧美一级黄| 欧美在线一区亚洲| 丝袜美足系列| 宅男免费午夜| 亚洲精品乱久久久久久| 欧美成人午夜精品| 男女之事视频高清在线观看 | 国产精品无大码| 夜夜骑夜夜射夜夜干| 99国产精品免费福利视频| 日本猛色少妇xxxxx猛交久久| 人人妻,人人澡人人爽秒播 | 免费观看av网站的网址| 精品一区二区三区av网在线观看 | 国产无遮挡羞羞视频在线观看| 久久av网站| 国产精品人妻久久久影院| 麻豆乱淫一区二区| 久久女婷五月综合色啪小说| 国产成人啪精品午夜网站| 女人被躁到高潮嗷嗷叫费观| 51午夜福利影视在线观看| 成人漫画全彩无遮挡| 天天躁狠狠躁夜夜躁狠狠躁| 日本av免费视频播放| 免费在线观看视频国产中文字幕亚洲 | 精品国产露脸久久av麻豆| 国产精品一区二区在线不卡| av在线老鸭窝| 超色免费av| 中文字幕精品免费在线观看视频| 亚洲久久久国产精品| 午夜老司机福利片| 最近的中文字幕免费完整| 免费少妇av软件| 人人妻人人澡人人看| 久久久久久免费高清国产稀缺| 91aial.com中文字幕在线观看| 国产乱来视频区| 在线观看免费日韩欧美大片| 黄色 视频免费看| 欧美日韩一级在线毛片| 九九爱精品视频在线观看| 日本vs欧美在线观看视频| 国产爽快片一区二区三区| 国产麻豆69| 极品少妇高潮喷水抽搐| 菩萨蛮人人尽说江南好唐韦庄| 黑人巨大精品欧美一区二区蜜桃| 只有这里有精品99| 欧美日韩亚洲高清精品| 久久久精品国产亚洲av高清涩受| 91成人精品电影| 精品久久蜜臀av无| 精品人妻在线不人妻| 成人漫画全彩无遮挡| 99国产精品免费福利视频| 国精品久久久久久国模美| 免费高清在线观看日韩| 久久精品亚洲熟妇少妇任你| 亚洲精品日韩在线中文字幕| 丰满饥渴人妻一区二区三| 如日韩欧美国产精品一区二区三区| 国产成人欧美| 亚洲国产av影院在线观看| 亚洲一码二码三码区别大吗| 啦啦啦中文免费视频观看日本| 国产一区二区三区综合在线观看| 成人18禁高潮啪啪吃奶动态图| 国产成人一区二区在线| 国产免费视频播放在线视频| 欧美xxⅹ黑人| 国产乱来视频区| 又大又黄又爽视频免费| 国产野战对白在线观看| 一级毛片电影观看| 国产一区二区激情短视频 | 少妇猛男粗大的猛烈进出视频| 中文字幕人妻丝袜制服| 免费黄色在线免费观看| 国产精品国产三级专区第一集| 久久久久久人人人人人| 亚洲国产欧美在线一区| 欧美人与善性xxx| 丰满少妇做爰视频| 亚洲精品日韩在线中文字幕| 自线自在国产av| 精品午夜福利在线看| 欧美激情 高清一区二区三区| 91精品伊人久久大香线蕉| 国产福利在线免费观看视频| 婷婷成人精品国产| 一本—道久久a久久精品蜜桃钙片| 亚洲人成网站在线观看播放| 美女福利国产在线| 满18在线观看网站| 国产成人欧美| 精品午夜福利在线看| 欧美黄色片欧美黄色片| 精品久久久精品久久久| 最新的欧美精品一区二区| 伊人久久大香线蕉亚洲五| 亚洲国产欧美在线一区| 精品少妇黑人巨大在线播放| 亚洲欧美清纯卡通| 青草久久国产| 免费黄频网站在线观看国产| 精品第一国产精品| 国产有黄有色有爽视频| 亚洲精品第二区| 日韩伦理黄色片| 在线观看一区二区三区激情| 精品一区二区三区av网在线观看 | 国产成人a∨麻豆精品| 纵有疾风起免费观看全集完整版| 如何舔出高潮| av又黄又爽大尺度在线免费看| 人妻一区二区av| 一本久久精品| 久久性视频一级片| 久久久国产欧美日韩av| 国产一区二区 视频在线| 国产av精品麻豆| 亚洲精品成人av观看孕妇| 久久精品久久久久久噜噜老黄| 亚洲国产看品久久| 精品少妇久久久久久888优播| 免费观看av网站的网址| 超碰97精品在线观看| 国产成人欧美在线观看 | 99香蕉大伊视频| 国产色婷婷99| 啦啦啦中文免费视频观看日本| 国产成人精品久久久久久| 丰满饥渴人妻一区二区三| 桃花免费在线播放| 考比视频在线观看| 在线观看免费高清a一片| 伊人亚洲综合成人网| 晚上一个人看的免费电影| 成人三级做爰电影| 国产在线免费精品| 丝袜脚勾引网站| 在线观看免费午夜福利视频| 亚洲精品久久成人aⅴ小说| 久久97久久精品| 男女下面插进去视频免费观看| 性色av一级| 可以免费在线观看a视频的电影网站 | e午夜精品久久久久久久| 丝袜喷水一区| 人人澡人人妻人| 日韩熟女老妇一区二区性免费视频| 免费日韩欧美在线观看| 中文字幕人妻熟女乱码| 中文字幕亚洲精品专区| 日韩精品有码人妻一区| 婷婷色综合www| 国产精品免费视频内射| 欧美亚洲 丝袜 人妻 在线| 丰满少妇做爰视频| 欧美xxⅹ黑人| 亚洲精品成人av观看孕妇| 1024香蕉在线观看| 精品一区二区三区av网在线观看 | 亚洲第一青青草原| 国产一区二区 视频在线| 久久久欧美国产精品| 精品一区二区三区av网在线观看 | 国产日韩欧美在线精品| 久久久久国产一级毛片高清牌| 一级爰片在线观看| 成年女人毛片免费观看观看9 | 91精品国产国语对白视频| 午夜福利在线免费观看网站| 毛片一级片免费看久久久久| 精品国产一区二区三区久久久樱花| 久久国产亚洲av麻豆专区| 亚洲国产日韩一区二区| av卡一久久| 制服丝袜香蕉在线| 侵犯人妻中文字幕一二三四区| 国产熟女午夜一区二区三区| 高清在线视频一区二区三区| 又粗又硬又长又爽又黄的视频| 1024香蕉在线观看| 秋霞在线观看毛片| 各种免费的搞黄视频| 精品人妻熟女毛片av久久网站| 美女国产高潮福利片在线看| 日本wwww免费看| 高清视频免费观看一区二区| 老司机亚洲免费影院| 日韩不卡一区二区三区视频在线| 亚洲成人国产一区在线观看 | 精品亚洲成国产av| 最近中文字幕高清免费大全6| a级毛片在线看网站| 亚洲精品一区蜜桃| 国产日韩一区二区三区精品不卡| 亚洲av日韩在线播放| 亚洲图色成人| 我的亚洲天堂| 天天躁夜夜躁狠狠久久av| 亚洲精品一区蜜桃| 黄色怎么调成土黄色| 美女大奶头黄色视频| 精品亚洲成国产av| 国精品久久久久久国模美| 婷婷色麻豆天堂久久| 国产免费又黄又爽又色| 国产成人a∨麻豆精品| 成人手机av| 一区二区日韩欧美中文字幕| 午夜福利影视在线免费观看| 在线观看免费午夜福利视频| 宅男免费午夜| 亚洲国产精品成人久久小说| 少妇人妻久久综合中文|