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

    江西省耕地土壤氮磷生態(tài)化學(xué)計量空間變異特征及其影響因素*

    2017-12-11 05:32:58江葉楓葉英聰李偉峰
    土壤學(xué)報 2017年6期
    關(guān)鍵詞:施用量氮磷氮肥

    江葉楓 葉英聰 郭 熙? 饒 磊 孫 凱 李偉峰

    (1 江西農(nóng)業(yè)大學(xué)國土資源與環(huán)境學(xué)院,南昌 330045)

    (2 江西省鄱陽湖流域農(nóng)業(yè)資源與生態(tài)重點實驗室,南昌 330045)

    江西省耕地土壤氮磷生態(tài)化學(xué)計量空間變異特征及其影響因素*

    江葉楓1,2葉英聰2郭 熙1,2?饒 磊1,2孫 凱1,2李偉峰1,2

    (1 江西農(nóng)業(yè)大學(xué)國土資源與環(huán)境學(xué)院,南昌 330045)

    (2 江西省鄱陽湖流域農(nóng)業(yè)資源與生態(tài)重點實驗室,南昌 330045)

    基于江西省2012年測土配方項目采集的16 582個耕地表層(0~20 cm)土壤樣點數(shù)據(jù),探討省域尺度下耕地土壤氮磷生態(tài)化學(xué)計量空間變異特征及其影響因素。運用相關(guān)性分析、回歸分析、方差分析和普通克里格法分析了成土母質(zhì)、土壤類型、耕地利用方式、秸稈還田方式和氮肥施用量對耕地土壤氮磷生態(tài)化學(xué)計量空間變異的影響。結(jié)果表明,土壤氮磷比(N/P)在0.28~13.63之間,均值為3.38,變異系數(shù)為45.56%,呈中等程度的變異??臻g分布上,高值區(qū)主要分布在海拔相對較低或鄱陽湖平原地區(qū),低值區(qū)主要分布在九江市。成土母質(zhì)、土壤類型、耕地利用方式、秸稈還田方式和氮肥施用量對江西省耕地土壤N/P空間變異影響顯著(p<0.01),但影響程度有所不同。氮肥施用量能影響27.2%的土壤N/P空間變異,是引起江西省耕地土壤氮磷生態(tài)化學(xué)計量空間變異的主要因素。

    江西省;土壤氮磷比;生態(tài)化學(xué)計量;空間變異;氮肥施用量

    生態(tài)化學(xué)計量學(xué)(Ecological stoichiometry)是一門主要研究化學(xué)元素平衡和生態(tài)系統(tǒng)能量平衡的新興生態(tài)學(xué)科[1],為研究土壤—植物生態(tài)交互作用以及土壤碳(C)、氮(N)、磷(P)循環(huán)方面提供了新思路,廣泛應(yīng)用于植物動態(tài)變化、種群演替、物種競爭、化學(xué)元素循環(huán)乃至全球氣候變化,并取得了顯著成就[2-3]。其中,土壤系統(tǒng)中C、N、P元素循環(huán)已成為全球氣候變化背景下土壤學(xué)、環(huán)境學(xué)和生態(tài)計量學(xué)研究學(xué)者關(guān)注的焦點之一[4]。土壤N和P是植物礦物質(zhì)養(yǎng)分、有機養(yǎng)分和土壤肥力的重要組成部分,作為陸地生態(tài)系統(tǒng)中的2個限制性元素,在植物的生長發(fā)育、群落組成以及生態(tài)系統(tǒng)的結(jié)構(gòu)和功能等方面發(fā)揮著重要作用[5]。受隨機性因素與結(jié)構(gòu)性因素的共同作用,土壤N、P在空間分布上往往表現(xiàn)出一定的異質(zhì)性,會產(chǎn)生異質(zhì)的土壤環(huán)境,進而影響植物動態(tài)變化、種群演替和物種競爭[6]。目前,關(guān)于土壤N、P的研究主要集中在土壤N、P自身變異特征,而忽略了土壤N、P耦合平衡的關(guān)系。研究表明,與僅考慮土壤N、P自身變異特征相比,氮磷比(N/P)能夠更全面地闡述土壤N、P變異[7]。

    由于不同生態(tài)系統(tǒng)的土壤微生物活性差異、土壤N與P的循環(huán)與轉(zhuǎn)換以及有機質(zhì)分解與積累速度的不同,導(dǎo)致土壤N/P變化明顯。與海洋生態(tài)系統(tǒng)相比較,陸地生態(tài)系統(tǒng)具有更為多變的環(huán)境因子、種群特征以及人為干擾,所以,土壤N/P變化更為復(fù)雜[8]。近年來,有關(guān)陸地生態(tài)系統(tǒng)土壤N、P化學(xué)計量特征空間變異及其影響因素的報道越來越多。盧同平等[9]對我國典型自然帶土壤N/P空間變異特征及其影響因素進行研究,發(fā)現(xiàn)溫帶荒漠帶森林土壤N/P主要受海拔影響,而青藏高原高寒植被帶主要受氣候和地形的影響。陳思宇等[10]對甘南高原草地土壤N/P空間異質(zhì)性研究結(jié)果表明:土壤N/P具有較強的空間變異性,在一定程度上受海拔和植被的共同影響。鄔畏等[6]研究表明,群落演替明顯影響土壤N/P的空間異質(zhì)性,且土壤N/P也是群落演替的一個重要影響因素。上述研究為認識我國陸地生態(tài)系統(tǒng)土壤N、P化學(xué)計量學(xué)特征提供了寶貴的資料。目前,對陸地生態(tài)系統(tǒng)土壤N/P的研究大多集中在濕地、草地與森林等,而對于耕地生態(tài)系統(tǒng)的報道并不多見。受成土母質(zhì)、土壤類型以及耕地利用方式、秸稈還田方式、氮肥施用量等人為活動因素影響,導(dǎo)致耕地生態(tài)系統(tǒng)的土壤N、P較其他陸地生態(tài)系統(tǒng)空間異質(zhì)性更為明顯[4]。研究陸地生態(tài)系統(tǒng)中耕地土壤N/P空間變異特征及其影響因素是研究陸地生態(tài)系統(tǒng)中營養(yǎng)結(jié)構(gòu)變化、種群演替和化學(xué)元素循環(huán)的基本依據(jù)。

    江西省地處我國東南部,是長江中下游地區(qū)農(nóng)業(yè)人口較為密集的區(qū)域,也是我國從未間斷向外輸出商品糧的兩個省份之一。耕地土壤N、P化學(xué)計量特征是表征耕地質(zhì)量的重要因素,掌握耕地土壤N/P空間變異特征及其影響因素對農(nóng)業(yè)可持續(xù)發(fā)展、現(xiàn)代農(nóng)業(yè)管理與生態(tài)環(huán)境保護具有十分迫切的現(xiàn)實意義。本文基于江西省2012年測土配方施肥項目采集的16 582個耕地表層(0~20 cm)土壤樣點數(shù)據(jù),以耕地土壤N、P生態(tài)化學(xué)計量特征為研究對象,結(jié)合成土母質(zhì)、土壤類型、耕地利用方式、秸稈還田方式和氮肥施用量,運用單因素方差分析、回歸分析和普通克里格法,探討江西省耕地表層土壤N/P空間變異特征及其影響因素。以期為江西省耕地表層土壤氮磷調(diào)控、耕地土壤化學(xué)元素分異與循環(huán)乃至全球氣候變化響應(yīng)模式提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    研究區(qū)位于江西省,年平均氣溫20~30℃,雨量充沛,年均降水量1 340~2 000 mm,為亞熱帶濕潤氣候;地理坐標介于24°29′14″N~30°04′41″N,113°34′36″E~118°28′58″E之間,總面積為1.69×105km2,人口4 566萬,轄11個地級行政區(qū)、100個縣級行政區(qū)、1個國家級新區(qū)。省內(nèi)東、西、南三面環(huán)山,中部丘陵和河谷平原交錯分布,北部則為鄱陽湖平原;土地利用類型以耕地和林地為主(圖1(a)),其中,耕地面積為3.09×106hm2,約占全省土地總面積的23.21%。耕地利用方式有一季水田、兩季水田、水旱輪作、一季旱地和兩季旱地等。土壤類型主要有水稻土、紅壤、黃褐土、石灰土和潮土等。秸稈還田存在不還田、翻壓還田和覆蓋還田3種還田方式。糧食作物以水稻為主,水稻種植面積約占全省耕地總面積的80%。植被以常綠闊葉林為主,具有典型的亞熱帶森林植物群落。

    1.2 土壤樣品采集與測定

    2012年11月至2013年2月按照《全國耕地地力調(diào)查與質(zhì)量評價技術(shù)規(guī)程》采集土壤樣品。在考慮空間分布均勻性和樣點代表性的前提下,每個采樣點采集4個點并進行充分混合,用四分法留取1 kg土樣,并運用GPS記錄該樣點經(jīng)緯度,同時,詳細記錄該點成土母質(zhì)、土壤類型、土屬類型、耕地利用方式、秸稈還田方式和氮肥施用量等信息,共采集耕地表層(0~20 cm)土壤樣品16 582個(圖1(b))。土壤樣品經(jīng)過自然風干后,在實驗室磨碎過篩,采用全自動凱氏定氮儀(Kjeltec8400,F(xiàn)oss,瑞典)測定全氮(TN)含量[11],全磷(TP)含量采用鉬銻抗比色法[12]進行測定。

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

    圖1 江西省土地利用類型圖和樣點分布圖Fig. 1 Land-use type map andsoil sampling site distribution map of Jiangxi Province

    受采樣及實驗室分析誤差的影響,土壤氮磷含量測量結(jié)果存在異常值,從而對土壤N/P空間總體分布與統(tǒng)計分析產(chǎn)生影響。本文采用域值法(3倍的標準差)對采樣點數(shù)據(jù)的土壤TN與TP含量異常值進行剔除[13],得到剔除后總樣點為16 109個,本文相關(guān)研究均采用剔除異常值后的數(shù)據(jù)?;谌?∶10萬的地形圖生成數(shù)字高程模型(Digital elevation model,DEM)。高程、坡度、坡度變率、坡向、曲率、地形起伏度、河流動能指數(shù)均由DEM在ArcGIS10.2軟件中處理提取,各地形因子的計算公式見參考文獻[14]。本研究中各縣市氮肥施用量空間分布來源于1985—2015年《江西省統(tǒng)計年鑒》中各縣市平均氮肥施用折純量。成土母質(zhì)、土壤類型、耕地利用方式、秸稈還田方式和氮肥施用量為定性變量,本文采用虛擬變量[15]進行賦值。

    2 結(jié) 果

    2.1 江西省土壤氮、磷化學(xué)計量描述性統(tǒng)計

    土壤氮磷描述性統(tǒng)計分析結(jié)果見表1。江西省耕地表層土壤TN、TP和N/P(質(zhì)量比)均值分別為1.58 g kg-1、0.52 g kg-1和3.38。TN較我國平均水平(1.54 g kg-1)[16]高0.04 g kg-1,TP處于我國中等水平(0.2~1.1 g kg-1)[11],N/P略低于我國平均水平(4.65)[9];從變異系數(shù)可以看出,N/P變異系數(shù)最大,其次為TP,變異系數(shù)最小的為TN,三者均呈中等程度的變異性。

    2.2 江西省土壤氮、磷空間結(jié)構(gòu)特征

    在GS+軟件中對TN、TP和N/P進行半方差函數(shù)的擬合,用半方差函數(shù)描述其空間結(jié)構(gòu)特征。從半方差擬合結(jié)果可以看出(圖2和表2),江西省土壤TN、TP和N/P的最優(yōu)模型分別為球狀模型、高斯模型和球狀模型,各模型的擬合效果較好,決定系數(shù)均在0.85以上。從模型的參數(shù)來看(表2),TN、TP和N/P的塊金效應(yīng)(隨機性因素引起的空間變異與系統(tǒng)總空間變異的比值)[15,17]分別為62.55%、64.52%和60.54%,表明TN、TP和N/P均呈弱空間相關(guān)性,隨機性因素引起的空間變異程度大于結(jié)構(gòu)性因素。TN、TP和N/P的變程分別為93.78、33.41和18.36km,表明TN的空間自相關(guān)范圍較大,而TP和N/P的空間自相關(guān)范圍較小。

    表1 江西省土壤全氮、全磷和氮磷比的描述性統(tǒng)計特征Table 1 Descriptive statistics of total nitrogen(TN),total phosphorus(TP)and N/P ratio in the soil of Jiangxi Province

    圖2 江西省土壤TN、TP和N/P半方差函數(shù)圖Fig. 2 Isotropic semivariogram of TN,TP and N/P ratio in the soil of Jiangxi Province

    表2 土壤TN、TP和N/P的半方差函數(shù)參數(shù)Table 2 Semivariance parameters of TN,TP and N/P ratio in the soil of Jiangxi

    2.3 江西省土壤氮、磷空間變異特征

    為直觀反映江西省土壤N/P的空間變異特征,在半方差模型擬合的基礎(chǔ)上運用普通克里格法對研究區(qū)土壤TN、TP和N/P進行空間插值,進而得到三者的空間分布。如圖3所示,TN與N/P的空間分布圖較為相似(圖3(a)和圖3(c)),高低分布趨勢相當,空間分布為突變而非漸變,能在一定程度上體現(xiàn)二者的空間異質(zhì)性;而TP分布較平滑,其中,含量在0.4~0.6 g kg-1之間所占面積比例最大。TN與N/P的高值區(qū)域主要分布在上饒市婺源縣、南昌市進賢縣、萍鄉(xiāng)市湘東區(qū)和蓮花縣、新余市分宜縣、撫州市黎川縣以及贛州市定南縣,即主要分布在海拔相對較低或鄱陽湖平原地區(qū)(圖1(b));低值區(qū)主要分布在贛州市信豐縣和九江市九江縣、德安縣、湖口縣、永修縣與彭澤縣(圖3(c))。大部分區(qū)域N/P位于2.0~5.0之間。

    2.4 土壤N/P空間變異的影響因素

    成土母質(zhì)通過影響土壤團聚體數(shù)量及其穩(wěn)定性、物理化學(xué)組成和風化淋溶進程等引起土壤氮磷含量的差異。成土母質(zhì)對江西省土壤N/P空間變異影響極顯著(p<0.01,表3)。不同成土母質(zhì)類型間N/P均值差異明顯(表4)。以碳酸巖類風化物最高(3.66),下蜀系黃土狀物最低(2.33)。從變異系數(shù)可以看出(表4),不同母質(zhì)間N/P均呈中等變異性,變異系數(shù)最高為67.38%(下蜀系黃土狀物),最低為45.14%(紅砂巖類風化物)。

    不同土壤類型具有不同的成土過程、礦物組成和發(fā)育程度,從而影響土壤氮磷含量。土壤類型對江西省土壤N/P空間變異影響極顯著(p<0.01,表3)。5種土壤類型中N/P均值大小依次為石灰土(3.84)>水稻土(3.40)>紅壤(3.30)>潮土(2.97)>黃褐土(2.02)。變異系數(shù)處于為33.44%~73.27%之間,以黃褐土變異系數(shù)最大(73.27%),麻沙泥田變異系數(shù)最?。?3.44%),均屬于中等變異性。在水稻土的3種亞類中,以淹育型水稻土均值最低(3.14),而潛育型水稻土均值最高(3.59),這反映了水稻土不同亞類之間由于附加的成土過程導(dǎo)致其N/P均值已出現(xiàn)了一定的差異。從土屬來看,土屬值域范圍為1.94~4.04,差異較大。

    圖3 江西省土壤TN、TP、N/P和氮肥施用量空間分布Fig. 3 Spatial distributions of TN,TP,N/P ratio and nitrogen fertilization rate in the soil of Jiangxi Province

    表3 不同因素對土壤N/P的回歸分析Table 3 Regression analysis of N/P ratio with affecting factors

    表4 不同成土母質(zhì)類型土壤N/P描述性統(tǒng)計特征Table 4 Descriptive statistic characteristics of N/P ratio relative to parent material

    耕地利用方式影響土壤肥料和植物殘體的數(shù)量和性質(zhì)、土壤水分含量和耕作管理措施,因而,不同利用方式下土壤氮磷含量存在差異性。耕地利用方式對江西省土壤N/P空間變異影響極顯著(p<0.01,表3)。從表5可以看出,均值大小依次為:水旱輪作(3.47)=一季旱地(3.47)>兩季水田(3.39)>一季水田(3.37)>兩季旱地(2.48)。其原因一方面與耕地休耕有關(guān);另一方面,與水田和旱地的土壤微生物活性及土壤固氮和固磷的能力有關(guān)。兩季旱地由于土地開發(fā)利用強度大、施肥量少于水田且地表作物大多被人為收取,歸還量小,使得土壤N/P最低。水旱輪作條件下雖然土壤開發(fā)利用強度大,但水田耕作條件下氮肥施用量也大,而且長期的干濕交替會降低土壤磷的有效性[18],因此,水旱輪作與一季旱地均值相近。兩季水田與一季水田均值含量亦相近。變異系數(shù)處于44.84%~58.06%之間,均表現(xiàn)為中等程度的變異性。

    秸稈還田通過改良土壤結(jié)構(gòu)與物理化學(xué)性質(zhì),促進微生物活力與作物根系發(fā)育并提高氮肥利用率,從而提高土壤氮含量,改善耕地土壤肥力。從表6可以看出,秸稈不還田與還田土壤N/P差異明顯,秸稈還田較不還田土壤N/P增加52.28%~53.11%,這是因為,作物秸稈N、P含量不一且土壤氮含量要高于磷含量,導(dǎo)致土壤氮的積累速度要大于磷,因此,秸稈還田較不還田土壤N/P呈增加趨勢。翻壓和覆蓋兩種還田方式下土壤N/P均值無顯著差異(p=0.700)。從變異系數(shù)的范圍(40.33%~53.11%)可以得出,不同秸稈還田方式均呈中等程度的變異性。

    表5 不同耕地利用方式土壤N/P描述性統(tǒng)計特征Table 5 Descriptive statistics of N/P ratio relative to farmland-use type

    表6 不同秸稈還田方式下土壤N/P描述性統(tǒng)計特征值Table 6 Descriptive statistics of N/P ratio relative to straw returningpattern

    氮肥施用量通過控制土壤氮積累速度來影響土壤N/P空間變異。氮肥施用量對江西省土壤N/P空間變異影響極顯著(p<0.001,表3)。具體表現(xiàn)為:隨著氮肥施用量增大,土壤N/P均值呈上升趨勢(表7)。從變異系數(shù)看,氮肥施用量大于300 kg hm-2的變異系數(shù)最小,為42.44%;氮肥施用量小于100 kg hm-2的變異系數(shù)最大,為59.30%;不同氮肥施用量土壤N/P均呈中等變異性。

    表7 不同氮肥施用量土壤N/P描述性統(tǒng)計特征值Table 7 Descriptive statistics of N/P ratio relative to nitrogen fertilization rate

    為定量揭示成土母質(zhì)、土壤類型、耕地利用方式、秸稈還田方式和氮肥施用量等因素對研究區(qū)耕地土壤N/P空間變異的影響水平,對影響因素進行回歸分析(表3)。結(jié)果表明:成土母質(zhì)、土壤類型、耕地利用方式、秸稈還田方式和氮肥施用量對江西省土壤N/P空間變異影響程度不一。成土母質(zhì)對土壤N/P空間變異的影響程度較低,僅為2.4%,這與相關(guān)研究結(jié)果[19]一致。這是因為,江西省成土母質(zhì)主要為第四紀紅色黏土、第三紀紅砂巖和花崗巖的紅壤巖層,該巖層多為花崗巖、石灰?guī)r和紅色砂巖[18],不同巖組主要因沉淀環(huán)境和各巖組所占比例有所分異,但巖性總體上較為相似,導(dǎo)致土壤機械組成相近[19],成土母質(zhì)主要通過影響土壤機械組成來影響土壤N、P含量分異,因此,成土母質(zhì)對研究區(qū)土壤N/P影響有限。土壤類型的影響程度隨土壤分類級別降低而升高,土類、亞類和土屬的影響程度分別為1.8%、2.2%和3.0%,這與其他區(qū)域研究結(jié)果[7-8]一致。因為級別越低,反映的成土過程、母巖特性、耕作管理措施和環(huán)境背景等信息越多,其影響程度越高[20],但具體如何影響尚有待進一步研究。耕地利用方式的影響程度為2.2%,遠遠低于其他區(qū)域的研究結(jié)果[21]。究其原因:①土地利用方式的改變勢必會打破陸地生態(tài)系統(tǒng)原有的氮磷平衡[22-23],但本文只研究了耕地,無法驗證和比較與其他土地利用方式下土壤N/P的差異性。②研究表明,地上植被變化往往會導(dǎo)致土壤N、P含量及土壤N/P發(fā)生分異[24],其他土地利用方式如園地、林地和草地等土壤N/P變異主要源于母質(zhì)與植物凋落物,而耕地除了種植作物的殘骸與母質(zhì)外,氮肥與磷肥的施用對江西省耕地土壤N/P嚴重不協(xié)調(diào)的地區(qū)進行了一定程度上的調(diào)節(jié),導(dǎo)致僅耕地一種土地利用方式獨立解釋土壤N/P空間異質(zhì)性的能力較低。秸稈還田方式能影響土壤N/P空間變異的13.0%,影響程度較高。氮肥施用量影響土壤N/P空間的程度最高,達到27.2%,遠遠高于其他影響因素,是影響江西省耕地土壤N/P空間變異的主要因素。

    3 討 論

    3.1 地形因子、TN、TP與土壤N/P的關(guān)系

    對研究區(qū)土壤N/P與地形因子、TN及TP進行相關(guān)性分析,從表8可以看出,高程與土壤N/P相關(guān)性顯著,相關(guān)系數(shù)為-0.38,表明土壤N/P隨高程上升逐漸減少,這與其他研究區(qū)域結(jié)果[9-10]一致。不同高程下,水熱條件和人為活動強度的差異引起土壤N/P空間變異。土壤N/P與其他地形因子相關(guān)性不顯著,其原因可能在省域尺度下地形因子本身的變異較小,同時1∶10萬的地形圖生成的DEM數(shù)據(jù)精度較低,其派生的地形變量精度也相對較低,這也在一定程度影響地形因子與土壤N/P的相關(guān)性。TN和TP與土壤N/P的相關(guān)系數(shù)分別為0.681和-0.665(p<0.01),表明土壤N/P主要取決于TN與TP含量。

    3.2 秸稈還田方式和氮肥施用量對土壤N/P的影響

    20世紀70年代,江西省作物生產(chǎn)就開始由傳統(tǒng)的僅僅依賴土壤自身肥力與有限的有機肥轉(zhuǎn)變?yōu)榻斩掃€田與化肥施用(尤其是氮肥)中來。20世紀80年代以來,土壤N、P含量具有不同程度的提升,其原因一方面與江西省大力推廣秸稈還田有關(guān);另一更重要的原因是江西省“高氮高產(chǎn)”的種植觀念[25]。秸稈還田對土壤N/P空間變異影響相對較大,不同還田方式之間也有一定的差異。造成這種現(xiàn)象主要是由于耕地土壤N/P變化取決于土壤氮和磷的輸入與輸出的相對關(guān)系,但磷一般相對穩(wěn)定,有效性高,圖3(b)也說明了這一點,因此,耕地土壤N/P變化在很大程度上由土壤氮元素含量決定。土壤氮大多數(shù)以有機氮的形式存在于有機質(zhì)中,目前的研究均表明,秸稈還田能為耕地提供豐富的碳源,刺激微生物活性,降低作物吸收氮素的能力,減少氮素淋洗損失,但不同秸稈還田方式對微生物活性、秸稈腐解與礦化速率影響程度不一,導(dǎo)致對秸稈有機物質(zhì)釋放和土壤氮積累效果形成一定的差異。從空間分布圖可以看出(圖3(c)和圖3(d)),土壤N/P高值區(qū)域主要分布在氮肥施用量大于200 kg hm-2的縣市(贛州市全南縣、南昌市進賢縣、九江市瑞昌市、上饒市鄱陽縣),而低值區(qū)域主要分布在氮肥施用量小于100 kg hm-2的縣市(如九江市湖口縣、彭澤縣和德安縣、萍鄉(xiāng)市袁州區(qū)、撫州市樂安縣),表現(xiàn)為隨著氮肥施用量增加,土壤N/P升高的總體趨勢,與前人研究結(jié)果[26-27]一致。究其原因:①氮肥的施用提高了土壤TN含量,使得N/P整體上呈明顯上升的趨勢,江西省“高氮高產(chǎn)”的種植觀念導(dǎo)致農(nóng)戶傾向于氮肥的施用而忽略磷肥的施用,陳琦[26]研究表明,2011年江西省晚稻平均施氮量和施磷量分別為226.5和84.96 kg hm-2,2012年江西省早稻平均施氮量和施磷量分別為197.7和93.34 kg hm-2。②土壤酸堿度是影響土壤固磷的重要影響因子之一[13],氮肥的施用降低了土壤酸堿度[27],使磷更容易發(fā)生沉淀反應(yīng)和吸附反應(yīng),形成羥基磷灰石或氟磷灰石,增加了磷的固定作用,減少了磷的有效性。氮肥施用量在本研究所有影響因素中對土壤N/P空間變異的影響程度最高,達到27.2%,遠高于其他因素。

    表8 土壤N/P與地形因子、全氮、全磷之間的相關(guān)系數(shù)Table 8 Correlation coefficients of N/P ratio with terrain factor,TN and TP

    3.3 土壤N/P變化的重要意義

    生物固氮是生態(tài)系統(tǒng)中主要的氮源。研究表明,無論是在陸地生態(tài)系統(tǒng)還是海洋生態(tài)系統(tǒng)中,生物固氮量隨土壤N/P增大而減少,減少而增大[28]。從表9可以看出,江西省耕地土壤生態(tài)系統(tǒng)的N/P均值要高于青海封育芨芨草原帶、高原過渡帶、溫帶草原帶、青藏高原高寒草原帶和暖溫帶闊葉落葉帶,但低于東北農(nóng)田、閩江河口濕地、熱帶季雨林、亞熱帶常綠闊葉林、青藏高原高寒植被帶和溫帶荒漠帶。這表明,江西省耕地生態(tài)系統(tǒng)土壤N/P處于中等偏下水平,同時也說明耕地生態(tài)系統(tǒng)生物固氮量處于較高水平?!督魇〗y(tǒng)計年鑒》數(shù)據(jù)顯示,近10年來,江西省氮肥施用量從2005年的47.75萬t下降至2015年的42.36萬t,氮肥施用量呈下降趨勢,而糧食產(chǎn)量卻在逐年上升,反映出在江西省目前N/P的情況下,若繼續(xù)加大氮肥的施用對土壤N、P耦合平衡及提高耕地生態(tài)系統(tǒng)生物固氮能力方面意義甚微,且過量地施用氮肥不但對生物固氮與固碳能力是一種削弱,而且氮肥殘留量或損失至環(huán)境中的量會顯著增加,造成土壤酸化和水體污染等農(nóng)業(yè)面源污染問題[29]。因此,保持穩(wěn)定的N/P是實現(xiàn)糧食增產(chǎn)與耕地生態(tài)系統(tǒng)生物固氮與固碳能力提升的主要途徑,一方面,在高氮高產(chǎn)的種植觀念下應(yīng)加大宣傳和相關(guān)技術(shù)培訓(xùn)力度,提高農(nóng)戶對化肥(尤其是氮肥)的認知水平,同時大力推廣測土配方施肥技術(shù),促進農(nóng)戶合理施肥;另一方面,增加對有機肥施用的補貼和政策優(yōu)惠,加快綠肥產(chǎn)業(yè)發(fā)展,同時配合秸稈還田。

    表9 耕地生態(tài)系統(tǒng)與其他地區(qū)的土壤N/PTable 9 Farmland ecosystem and soil N/P in other areas

    4 結(jié) 論

    江西省耕地表層土壤N/P表現(xiàn)為中等程度的變異性,空間分布格局受隨機性因素和結(jié)構(gòu)性因素共同影響,但主要受隨機性因素影響??臻g分布方面,高值區(qū)域主要分布在上饒市婺源縣、南昌市進賢縣、萍鄉(xiāng)市湘東區(qū)和蓮花縣、新余市分宜縣、撫州市黎川縣以及贛州市定南縣,低值區(qū)域主要分布在贛州市信豐縣和九江市九江縣、德安縣、湖口縣、永修縣與彭澤縣。成土母質(zhì)、土壤類型、耕地利用方式、秸稈還田方式和氮肥施用量對江西省耕地土壤N/P空間變異均有顯著影響(p<0.01)。在所有因素中,氮肥施用量對土壤N/P空間變異的影響程度最大,遠遠大于其他因素,是影響江西省耕地土壤N/P空間變異的主要因素。

    [1] Sterner R W,Elser J J. Ecological stoichiometry:The biology of elements from molecules to the biosphere.Princeton,New Jersey:Princeton University Press,2002:225—226

    [2] Tessier J T,Raynal D J. Use of nitrogen to phosphorus ratios in plant tissue as an indicator of nutrient limitation and nitrogen saturation. Journal of Applied Ecology,2003,40(3):523—534

    [3] Verhoeven J T,Koerselman W,Meuleman A F.Nitrogen-or phosphorus-limited growth in herbaceous,wet vegetation:Relations with atmospheric inputs and management regimes. Trends in Ecology amp; Evolution,1996,11(12):494—497

    [4] 王紹強,于貴瑞. 生態(tài)系統(tǒng)碳氮磷元素的生態(tài)化學(xué)計量學(xué)特征. 生態(tài)學(xué)報,2008,28(8):3937—3947 Wang S Q,Yu G R. Ecological stoichiometry characteristics of ecosystem carbon,nitrogen and phosphorus elements(In Chinese). Acta Ecologica Sinica,2008,28(8):3937—3947

    [5] Han W,F(xiàn)ang J,Guo D,et al. Leaf nitrogen and phosphorus stoichiometry across 753 terrestrial plant species in China. New Phytologist,2005,168(2):377—385

    [6] 鄔畏,何興東,周啟星. 生態(tài)系統(tǒng)氮磷比化學(xué)計量特征研究進展. 中國沙漠,2010,30(2):296—302 Wu W,He X D,Zhou Q X. Review on N:P stoichiometry in eco-system(In Chinese). Journal of Desert Research,2010,30(2):296—302

    [7] 王建林,鐘志明,王忠紅,等. 青藏高原高寒草原生態(tài)系統(tǒng)土壤氮磷比的分布特征. 應(yīng)用生態(tài)學(xué)報,2013,24(12):3399—3406 Wang J L,Zhong Z M,Wang Z H,et al. Soil N/P ratio distribution characteristics of alpine grassland ecosystem in Qinghai-Tibet Plateau(In Chinese). Chinese Journal of Applied Ecology,2013,24(12):3399—3406

    [8] Jenny H J. The factors of soil formation. New York:McGraw-Hill,1941

    [9] 盧同平,張文翔,牛潔,等. 典型自然帶土壤氮磷化學(xué)計量空間分異特征及其驅(qū)動因素研究. 土壤學(xué)報,2017,54(3):681—691 Lu T P,Zhang W X,Niu J,et al. Study on spatial variability and driving factors of stoichiometry of nitrogen and phosphorus in soils of typical natural zones of China(In Chinese). Acta Pedologica Sinica,2017,54(3):681—691

    [10] 陳思宇,于惠,馮琦勝,等. 甘南高原土壤氮磷比空間異質(zhì)性研究. 草地學(xué)報,2013,21(1):30—36 Chen S Y,Yu H,F(xiàn)eng Q S,et al. Spatial heterogeneity of soil N/P ratio on Gannan Plateau(In Chinese). Acta Agrestia Sinica,2013,21(1):30—36

    [11] 黃昌勇. 土壤學(xué). 北京:中國農(nóng)業(yè)出版社,2000:199—200 Huang C Y. Soil science(In Chinese). Beijing:China Agriculture Press,2000:199—200

    [12] 魯如坤. 土壤農(nóng)業(yè)化學(xué)分析方法. 北京:中國農(nóng)業(yè)科學(xué)技術(shù)出版社,2000 Lu R K. Analyticalmethodsforsoil and agro-chemistry(In Chinese). Beijing:China Agricultural Science and Technology Press,2000

    [13] 王丹丹,史學(xué)正,于東升,等. 東北地區(qū)旱地土壤有機碳密度的主控自然因素研究. 生態(tài)環(huán)境學(xué)報,2009,18(3):1049—1053 Wang D D,Shi X Z,Yu D S,et al. Main natural factors controlling soil organic carbon density in the uplands of Northeast China(In Chinese). Ecology and Environmental Sciences,2009,18(3):1049—1053

    [14] 張素梅,王宗明,張柏,等. 利用地形和遙感數(shù)據(jù)預(yù)測土壤養(yǎng)分空間分布. 農(nóng)業(yè)工程學(xué)報,2010,26(5):188—194 Zhang S M,Wang Z M,Zhang B,et al. Prediction of spatial distribution of soil nutrients using terrain attributes and remote sensing data(In Chinese).Transactions of the Chinese Society of Agricultural Engineering,2010,26(5):189—194

    [15] 羅由林,李啟權(quán),王昌全,等. 四川省仁壽縣土壤有機碳空間分布特征及其主控因素. 中國生態(tài)農(nóng)業(yè)學(xué)報,2015,23(1):34—42 Luo Y L,Li Q Q,Wang C Q,et al. Spatial variability of soil organic carbon and related controlling factors in Renshou County,Sichuan Province(In Chinese).Chinese Journal of Eco-Agriculture,2015,23(1):34—42

    [16] 李啟權(quán),岳天祥,范澤孟,等. 中國表層土壤全氮的空間模擬分析. 地理研究,2010,29(11):1981—1992 Li Q Q,Yue T X,F(xiàn)an Z M,et al. Spatial simulation of topsoil TN at the national scale in China(In Chinese).Geographical Research,2010,29(11):1981—1992

    [17] Zhu Q,Lin H S. Comparing ordinary kriging and regression kriging for soil properties in contrasting landscapes. Pedosphere,2010,20(5):594—606

    [18] 江西省土地利用管理局. 江西土壤. 北京:中國農(nóng)業(yè)科技出版社,1991 Administration Bureau of Land,Soil Survey Staff Office of Jiangxi Province. Soils of Jiangxi Province(In Chinese). Beijing:China Agricultural Science and Technology Press,1991

    [19] 胡玉福,鄧良基,張世熔,等. 川中丘陵區(qū)城墻巖群和蓬萊鎮(zhèn)組紫色巖上土壤性質(zhì)研究. 土壤通報,2007,38(6):1076—1080 Hu Y F,Deng L J,Zhang S R,et al. The soil properties of soils developed from K1cgand J3pin central hill region of Sichuan basin(In Chinese). Chinese Journal of Soil Science,2007,38(6):1076—1080

    [20] Bell M J,Worrall F. Estimating a region’s soil organic carbon baseline:The undervalued role of land management. Geoderma,2009,152(1/2):74—84

    [21] 夏品華,孔祥量,喻理飛. 草海濕地小流域土地利用與景觀格局對氮、磷輸出的影響. 環(huán)境科學(xué)學(xué)報,2016,36(8):2983—2989 Xia P H,Kong X L,Yu L F. Effects of land-use and landscape pattern on nitrogen and phosphorus exports in Caohai wetland watershed(In Chinese). Acta Scientiae Circumstantiae,2016,36(8):2983—2989

    [22] Guo L B,Gifford R M. Soil carbon stocks and land use change:Ameta analysis. Global Change Biology,2002,8(4):345—360

    [23] 周正虎,王傳寬,張全智. 土地利用變化對東北溫帶幼齡林土壤碳氮磷含量及其化學(xué)計量特征的影響. 生態(tài)學(xué)報,2015,35(20):6694—6702 Zhou Z H,Wang C K,Zhang Q Z. The effect of land use change on soil carbon,nitrogen,and phosphorus contents and their stoichiometry in temperate sapling stands in northeastern China(In Chinese). Acta Ecologica Sinica,2015,35(20):6694—6702

    [24] 張文彥,樊江文,鐘華平,等. 中國典型草原優(yōu)勢植物功能群氮磷化學(xué)計量學(xué)特征研究. 草地學(xué)報,2010,18(4):503—509 Zhang W Y,F(xiàn)an J W,Zhong H P,et al. The nitrogen:phosphorus stoichiometry of different plant functional groups for dominant species of typical steppes in China(In Chinese). Acta Agrectir Sinica,2010,18(4):503—509

    [25] 江葉楓,郭熙,孫凱,等. 江西省耕地土壤碳氮比空間變異特征及其影響因素. 環(huán)境科學(xué),2017,38(9):3840—3850 Jiang Y F,Guo X,Sun K,et al. Spatial variability of farmland soil C/N ratio of Jiangxi Province(In Chinese). Environmental Science,2017,38(9):3840—3850

    [26] 陳琦. 農(nóng)戶水稻施肥現(xiàn)狀調(diào)查與分析. 南京:南京農(nóng)業(yè)大學(xué),2013 Chen Q. Investigation and analysis of farmers’fertilization status on rice-taking Jiangxi Province and Jiangsu Province as examples(In Chinese). Nanjing:Nanjing Agricultural University,2013

    [27] Blake L,Goulding K W T,Mott C J B,et al. Changes in soil chemistry accompanying acidification over more than 100 years under woodland and grass at Rothamsted Experimental Station,UK. European Journal of Soil Science,1999,50(3):401—412

    [28] Smith V H. Effects of nitrogen:Phosphorus supply ratios on nitrogen fixation in agricultural and pastoral ecosystems. Biogeochemistry,1992,18(1):19—35

    [29] 巨曉棠. 氮肥有效率的概念及意義——兼論對傳統(tǒng)氮肥利用率的理解誤區(qū). 土壤學(xué)報,2014,51(5):921—933 Ju X T. The concept and meanings of nitrogen fertilizer availability ratio:Discussing misunderstanding of traditional nitrogen use efficiency(In Chinese). Acta Pedologica Sinica,2014,51(5):921—933

    [30] 張興義,隋躍宇,張少良,等. 薄層農(nóng)田黑土全量碳及氮磷鉀含量的空間異質(zhì)性. 水土保持通報,2008,28(2):1—5 Zhang X Y,Sui Y Y,Zhang S L,et al. Spatial heterogeneities of total carbon,nitrogen,phosphorus and potassium content in black thin-layer soil(In Chinese). Bulletin of Soil and Water Conservation,2008,28(2):1—5

    [31] 王維奇,仝川,賈瑞霞,等. 不同淹水頻率下濕地土壤碳氮磷生態(tài)化學(xué)計量學(xué)特征. 水土保持學(xué)報,2010,24(3):238—242 Wang W Q,Tong C,Jia R X,et al. Ecological stoichiometry characteristics of wetland soil carbon,nitrogenand phosphorus in different water–flooded frequency(In Chinese). Journal of Soil and Water Conservation,2010,24(3):238—242

    [32] 李亞娟,曹廣民,龍瑞軍. 不同草地利用方式對土壤有機碳、全氮和全磷的影響. 草原與草坪,2012,32(5):26—35 Li Y J,Cao G M,Long R J. Effects of grassland using models on soil carbon,nitrogen and phosphorus content(In Chinese). Grassland and Turf,2012,32(5):26—35

    [33] 王晶苑,王紹強,李紉蘭,等. 中國四種森林類型主要優(yōu)勢植物的C:N:P化學(xué)計量學(xué)特征. 植物生態(tài)學(xué)報,2011,35(6):587—595 Wang J Y,Wang S Q,Li R L,et al. C:N:P stoichiometric characteristics of four forest types’dominant tree species in China(In Chinese). Chinese Journal of Plant Ecology,2011,35(6):587—595

    (責任編輯:陳榮府)

    Spatial Variability of Ecological Stoichiometry of Soil Nitrogen and Phosphorus in Farmlands of Jiangxi Province and Its Inf l uencing Factors

    JIANG Yefeng1,2YE Yingcong2GUO Xi1,2?RAO Lei1,2SUN Kai1,2LI Weifeng1,2
    (1 Academy of Land Resource and Environment,Jiangxi Agricultural University,Nanchang 330045,China)
    (2 Key Laboratory of Poyang Lake Watershed Agricultural Resources and Ecology of Jiangxi Province,Nanchang 330045,China)

    【Objective】Nitrogen(N)and Phosphorus(P)are two essential mineral nutrients for plant growth,and important parts of organic nutrient and soil fertility. As two restrictive elements in the terrestrial ecosystem,they play important roles in growth and development of plants,community composition of the vegetation and structure and function of the ecosystem as well. As ecosystems vary in soil microbial activity,they do in recycling and transformation of soil N and P and in decomposition and accumulation rateof organic matter,thus leading to significant variation of the ecological stoichiometry of soil N/P.Therefore,the study on spatial variability of the ecological stoichiometry of soil N/P and its influencing factors in the terrestrial ecosystemat a provincial scale is the fundamental basis of the study on changes in trophic structure,population succession and chemical elements recycling in the terrestrial ecosystem.【Method】Based on the data of the 16 582 soil samples collected in the topsoil layer(0~20 cm)of the farmlands during implementation of the project of soil-test-based formulated fertilization in Jiangxi Province in 2012,values were assigned to the qualitative factors,i.e. parent material,soil type,farmland-use type,straw incorporation pattern and nitrogen fertilization(divided into 4 levels),and one-way variance analysis,ordinary kriging analysis and regression analysis were done of the data to quantify the ability of each factortoexplainspatial variability of the ecological stoichiometry of soil N/P ratio independently.【Result】Results show as follows:①The soil N/P ratio in the topsoil layer ranged from 0.28 to 13.63,with an average of 3.38,which indicates that the soil N/P in the cultivated land ecosystem of Jiangxi Province is at the level of medium on the lower side,and that the biological nitrogen fixation rate in the cultivated land ecosystem is at a quite high level. The coefficient of variation was 45.56%,indicating that soil N/P is moderate in variability. The nuggest/sill ratio was 60.54%,indicating that soil N and P is weak in spatial dependence,varying between structural factor and stochastic factor,and the latter played a more important role in spatial variability of soil N/P ratio. ②Soil N/P ratiowas significantly related to elevation(p<0.01,r=-0.38),TN(p<0.01,r=0.681)and TP(p<0.01,r=-0.665),suggesting that elevation is a major factor affecting spatial distribution of soil N/P and that soil N/P depends mainly on TN and TP contents.③The spatial distribution of soil N/P ratioappeared to be quite smooth relatively,with high values over in Wuyuanof Shangrao,Jinxianof Nanchang,Xiangdong and Lianhuaof Pingxiang,F(xiàn)enyiof Xinyu,Lichuanof Fuzhou and Dingnanof Ganzhou. ④The five factors,parent material,soil type,farmland-use type,straw incorporation pattern and nitrogen fertilization rate,all had significant impacts on spatial variability of soil N/P ratio(p<0.01),but varied in degree. Parent materialexplained 2.4% of the variability;type,sub-type and genus of the soil did 1.8%,2.2% and 3.0%,respectively;farmland-use type did 2.2%;straw incorporation pattern did 13.0%;and N fertilization rate did 27.2%.【Conclusion】The factor of N fertilization rate is far more capable of explaining independently spatial variability ofthe soil N/P ratio than the other factors,indicating that it is the main factor determining soil N/P ratio in the topsoil layer of Jiangxi Province. Therefore,in order to improve the ability of the farmland ecosystem to fix nitrogen biologically and sequestrate nitrogen,here are two suggestions. One is to devote more efforts to improving farmers’ knowledge about the use of chemical fertilizers,especially Nitrogen fertilizer,through propaganda and related technical training,vigorously extrapolating the soil-test-based formulated fertilization technology and stimulating farmers to rationalize fertilizer usage. And the other is to raise subsidies and implement preferential policies for use of organic manure and accelerate development of the green manure industry coupled with the extension of straw incorporation.

    Jiangxi Province;Soil N/P ratio;Ecological stoichiometry;Spatial variability;Nitrogen fertilization rate

    S154.1;Q148

    A

    10.11766/trxb201705310212

    * 國家自然科學(xué)基金項目(41361049)、江西省自然科學(xué)基金項目(20122BAB204012)和江西省贛鄱英才“555”領(lǐng)軍人才項目(201295) 資助 Supported by the National Natural Science Foundation of China(No. 41361049),the Natural Science Foundation of Jiangxi Province(No. 20122BAB204012)and the GanPo“555”Talent Research Funds of Jiangxi Province(No.201295)

    ? 通訊作者 Corresponding author,E-mail:xig435@163.com

    江葉楓(1994—),男,江西余干人,碩士研究生,主要從事土壤環(huán)境與系統(tǒng)模擬方面研究。E-mail:jiangyf0308@163.com

    2017-05-31;

    2017-07-10;優(yōu)先數(shù)字出版日期(www.cnki.net):2017-08-14

    猜你喜歡
    施用量氮磷氮肥
    玻利維亞擬建新的氮磷鉀肥料工廠
    春大豆施鉬條件下最適氮磷肥施用量研究
    氮肥供應(yīng)充足 春耕生產(chǎn)有保障
    江淮小氮肥 耕耘六十年——紀念安徽小氮肥誕生六十周年
    抓住機遇 主動作為 努力推進我國氮肥市場穩(wěn)步前行
    2017春季各地氮肥市場掠影
    復(fù)合微生物菌劑施用量對烤煙產(chǎn)量和質(zhì)量的影響
    植物需要氮磷鉀
    檳榔芋不同氮肥施用量試驗
    三料磷不同施用量對棉花生長和產(chǎn)量的影響
    国产又色又爽无遮挡免| 久久久久久久国产电影| 久久毛片免费看一区二区三区| 久久国产精品大桥未久av| 亚洲精品视频女| 一区二区av电影网| 色视频在线一区二区三区| 亚洲国产色片| 亚洲三级黄色毛片| 欧美激情极品国产一区二区三区 | 精品国产露脸久久av麻豆| 日本与韩国留学比较| 在线观看三级黄色| 大香蕉97超碰在线| 晚上一个人看的免费电影| 精品人妻在线不人妻| 尾随美女入室| 男女啪啪激烈高潮av片| 国产精品三级大全| 下体分泌物呈黄色| 男女无遮挡免费网站观看| 亚洲av二区三区四区| 最新的欧美精品一区二区| 黑人高潮一二区| 91久久精品国产一区二区成人| 在线免费观看不下载黄p国产| 极品少妇高潮喷水抽搐| 久久久久国产网址| 最黄视频免费看| 日韩强制内射视频| 少妇熟女欧美另类| 中文字幕久久专区| 我的老师免费观看完整版| 国产一区二区在线观看av| 久久精品久久精品一区二区三区| 嘟嘟电影网在线观看| 国产精品人妻久久久久久| 国产成人精品无人区| 国产成人91sexporn| av国产精品久久久久影院| 欧美97在线视频| 亚洲无线观看免费| 午夜日本视频在线| 搡老乐熟女国产| 下体分泌物呈黄色| 亚洲人成网站在线观看播放| 亚洲av国产av综合av卡| 中文字幕最新亚洲高清| 建设人人有责人人尽责人人享有的| 欧美一级a爱片免费观看看| 黑丝袜美女国产一区| 国产午夜精品久久久久久一区二区三区| 亚州av有码| 久久午夜福利片| 国产精品蜜桃在线观看| 韩国av在线不卡| 久久韩国三级中文字幕| 日本猛色少妇xxxxx猛交久久| 777米奇影视久久| 大片电影免费在线观看免费| 国产综合精华液| 人人妻人人添人人爽欧美一区卜| 性色avwww在线观看| 最后的刺客免费高清国语| 国语对白做爰xxxⅹ性视频网站| 一区二区三区乱码不卡18| 欧美日韩国产mv在线观看视频| 精品久久久精品久久久| 久久人人爽人人片av| 国产69精品久久久久777片| 香蕉精品网在线| 黑人巨大精品欧美一区二区蜜桃 | 黄色欧美视频在线观看| 欧美激情国产日韩精品一区| 免费观看在线日韩| 日本av免费视频播放| 日本午夜av视频| 国产精品一区二区在线不卡| 久久久久久久大尺度免费视频| 在线播放无遮挡| 超色免费av| 国产日韩一区二区三区精品不卡 | 国产有黄有色有爽视频| 自拍欧美九色日韩亚洲蝌蚪91| 我要看黄色一级片免费的| 国产亚洲精品第一综合不卡 | 国产av一区二区精品久久| 欧美日韩一区二区视频在线观看视频在线| 97超视频在线观看视频| 天天操日日干夜夜撸| 一本—道久久a久久精品蜜桃钙片| av国产精品久久久久影院| 91精品国产九色| 国产欧美亚洲国产| 国产精品国产三级国产av玫瑰| 日韩免费高清中文字幕av| 国产日韩欧美亚洲二区| 亚洲人成77777在线视频| 国产一区二区在线观看日韩| 中文字幕精品免费在线观看视频 | 人体艺术视频欧美日本| 一级片'在线观看视频| 丝袜脚勾引网站| 精品人妻在线不人妻| 午夜日本视频在线| 久久影院123| 男人爽女人下面视频在线观看| 成年女人在线观看亚洲视频| 91久久精品电影网| 亚洲欧洲日产国产| 美女视频免费永久观看网站| 成人免费观看视频高清| 精品国产乱码久久久久久小说| 国产高清有码在线观看视频| 中文字幕免费在线视频6| 黄色配什么色好看| 国产永久视频网站| 国产精品国产三级国产av玫瑰| 国产亚洲精品久久久com| 三上悠亚av全集在线观看| 精品少妇黑人巨大在线播放| 中文字幕人妻丝袜制服| 男女啪啪激烈高潮av片| 观看av在线不卡| 亚洲精品乱久久久久久| 男女高潮啪啪啪动态图| 啦啦啦中文免费视频观看日本| 午夜老司机福利剧场| 国产精品久久久久久精品电影小说| 亚洲五月色婷婷综合| 久热久热在线精品观看| 国产乱人偷精品视频| 欧美97在线视频| 高清视频免费观看一区二区| 日韩熟女老妇一区二区性免费视频| 在线精品无人区一区二区三| 亚洲精华国产精华液的使用体验| 青春草视频在线免费观看| 考比视频在线观看| www.av在线官网国产| 免费观看av网站的网址| 亚洲激情五月婷婷啪啪| 国产极品粉嫩免费观看在线 | 一本大道久久a久久精品| 97精品久久久久久久久久精品| 久热久热在线精品观看| tube8黄色片| 日本黄色片子视频| 欧美 日韩 精品 国产| 日日摸夜夜添夜夜添av毛片| 日本av免费视频播放| 国产高清有码在线观看视频| 内地一区二区视频在线| 国产白丝娇喘喷水9色精品| 精品国产一区二区久久| 嫩草影院入口| 新久久久久国产一级毛片| av.在线天堂| 国产亚洲欧美精品永久| 熟女人妻精品中文字幕| 午夜免费男女啪啪视频观看| 久久99蜜桃精品久久| 人人妻人人添人人爽欧美一区卜| 国产69精品久久久久777片| 人妻制服诱惑在线中文字幕| 欧美bdsm另类| 精品少妇内射三级| 中文字幕免费在线视频6| 人人妻人人澡人人爽人人夜夜| 久久久久精品久久久久真实原创| 好男人视频免费观看在线| 老熟女久久久| 99视频精品全部免费 在线| 亚洲美女视频黄频| 最近中文字幕高清免费大全6| 国产精品一区www在线观看| 亚洲无线观看免费| 国产极品粉嫩免费观看在线 | 日韩成人伦理影院| av播播在线观看一区| 日本爱情动作片www.在线观看| 久久久久久久国产电影| 久久久久久久久久人人人人人人| 草草在线视频免费看| 一级毛片我不卡| 久久久国产欧美日韩av| 久久99一区二区三区| 国产亚洲午夜精品一区二区久久| 99热6这里只有精品| 午夜福利,免费看| 国产国语露脸激情在线看| 欧美日韩国产mv在线观看视频| 一级爰片在线观看| 边亲边吃奶的免费视频| 伦理电影免费视频| 香蕉精品网在线| 九九在线视频观看精品| 久久久国产一区二区| 久久99热6这里只有精品| 不卡视频在线观看欧美| 搡女人真爽免费视频火全软件| 午夜老司机福利剧场| 亚洲三级黄色毛片| 久久久午夜欧美精品| 亚洲情色 制服丝袜| 欧美 日韩 精品 国产| 黑丝袜美女国产一区| 亚洲天堂av无毛| 高清不卡的av网站| 少妇精品久久久久久久| 欧美精品一区二区大全| av播播在线观看一区| 日日爽夜夜爽网站| 国产精品麻豆人妻色哟哟久久| 永久免费av网站大全| 丰满迷人的少妇在线观看| 亚洲人成网站在线播| 寂寞人妻少妇视频99o| 黑人猛操日本美女一级片| 少妇熟女欧美另类| 国产白丝娇喘喷水9色精品| 国产国拍精品亚洲av在线观看| 国产在视频线精品| 如何舔出高潮| 一区二区日韩欧美中文字幕 | 国产成人精品一,二区| 国产成人freesex在线| 纵有疾风起免费观看全集完整版| 久久久久精品久久久久真实原创| 中文精品一卡2卡3卡4更新| 丰满乱子伦码专区| 久久久精品94久久精品| 777米奇影视久久| 久久 成人 亚洲| 午夜激情福利司机影院| 日韩视频在线欧美| 久久国产精品男人的天堂亚洲 | 国产精品99久久久久久久久| 天天躁夜夜躁狠狠久久av| 日韩中文字幕视频在线看片| 免费观看的影片在线观看| 亚洲精品一二三| 菩萨蛮人人尽说江南好唐韦庄| 国产 精品1| 国产精品一区二区三区四区免费观看| 啦啦啦在线观看免费高清www| 蜜桃国产av成人99| 亚洲av日韩在线播放| 观看美女的网站| 老女人水多毛片| 亚洲精品久久久久久婷婷小说| 丝袜美足系列| 考比视频在线观看| 亚洲美女视频黄频| 最近的中文字幕免费完整| 欧美最新免费一区二区三区| 男女高潮啪啪啪动态图| 亚洲av在线观看美女高潮| 精品人妻偷拍中文字幕| 好男人视频免费观看在线| 日日摸夜夜添夜夜添av毛片| 成人国语在线视频| 国产69精品久久久久777片| 最近中文字幕2019免费版| 亚洲国产精品一区三区| 乱码一卡2卡4卡精品| 街头女战士在线观看网站| 午夜老司机福利剧场| 人妻一区二区av| 亚洲av福利一区| 午夜精品国产一区二区电影| av网站免费在线观看视频| 卡戴珊不雅视频在线播放| 一区二区日韩欧美中文字幕 | 亚洲国产欧美在线一区| 18禁动态无遮挡网站| 日韩人妻高清精品专区| 午夜91福利影院| 尾随美女入室| 国产国拍精品亚洲av在线观看| 香蕉精品网在线| 久久久久久久久久久久大奶| 91在线精品国自产拍蜜月| 天天影视国产精品| 五月伊人婷婷丁香| 亚洲欧美一区二区三区国产| 国产高清不卡午夜福利| 亚洲情色 制服丝袜| 一级毛片aaaaaa免费看小| 精品国产一区二区久久| 涩涩av久久男人的天堂| 91成人精品电影| 免费播放大片免费观看视频在线观看| 成人无遮挡网站| 久久精品熟女亚洲av麻豆精品| 午夜福利在线观看免费完整高清在| 亚洲精品色激情综合| 人妻人人澡人人爽人人| 一级爰片在线观看| 成年美女黄网站色视频大全免费 | 国产精品99久久99久久久不卡 | 少妇被粗大猛烈的视频| 简卡轻食公司| 久久久久久久久久成人| 精品久久久久久电影网| 亚洲欧美一区二区三区国产| 少妇高潮的动态图| 免费不卡的大黄色大毛片视频在线观看| 五月天丁香电影| 国产在线视频一区二区| 亚洲一级一片aⅴ在线观看| 国产视频内射| 夫妻性生交免费视频一级片| 日本猛色少妇xxxxx猛交久久| 亚洲精品456在线播放app| 久久人妻熟女aⅴ| 有码 亚洲区| 亚洲av不卡在线观看| 夫妻性生交免费视频一级片| xxx大片免费视频| 国产成人aa在线观看| 国产欧美亚洲国产| 在线观看三级黄色| 99久国产av精品国产电影| 一级爰片在线观看| 永久免费av网站大全| 2018国产大陆天天弄谢| 黄色视频在线播放观看不卡| 18+在线观看网站| av免费观看日本| 性色avwww在线观看| 久久精品熟女亚洲av麻豆精品| 久久99蜜桃精品久久| 久久国产亚洲av麻豆专区| 国产乱人偷精品视频| 边亲边吃奶的免费视频| 国产精品一区二区在线不卡| 欧美日韩亚洲高清精品| 制服丝袜香蕉在线| 免费av中文字幕在线| 国产精品一二三区在线看| 人妻制服诱惑在线中文字幕| 久久久久久伊人网av| 久久精品夜色国产| 看十八女毛片水多多多| 少妇人妻久久综合中文| 国产一级毛片在线| 久久久午夜欧美精品| 丝袜脚勾引网站| 丝袜在线中文字幕| 18在线观看网站| 亚洲欧美一区二区三区国产| 亚洲av欧美aⅴ国产| 中文字幕精品免费在线观看视频 | 99热国产这里只有精品6| 少妇熟女欧美另类| 婷婷色麻豆天堂久久| 亚洲国产精品成人久久小说| 中文字幕人妻丝袜制服| 婷婷色综合www| 夫妻性生交免费视频一级片| 啦啦啦视频在线资源免费观看| 国产成人freesex在线| 国产精品女同一区二区软件| 考比视频在线观看| 精品久久国产蜜桃| 免费人成在线观看视频色| 大码成人一级视频| 国产高清三级在线| av免费在线看不卡| 十八禁高潮呻吟视频| 日韩不卡一区二区三区视频在线| 中文字幕制服av| 免费观看a级毛片全部| 99热这里只有是精品在线观看| 日韩一本色道免费dvd| 亚洲国产精品成人久久小说| 天天操日日干夜夜撸| 国产老妇伦熟女老妇高清| 寂寞人妻少妇视频99o| 欧美 亚洲 国产 日韩一| 色网站视频免费| 好男人视频免费观看在线| 久久精品人人爽人人爽视色| 中文字幕av电影在线播放| 国产一区二区三区综合在线观看 | 国产精品蜜桃在线观看| 日本午夜av视频| 久热久热在线精品观看| 国产 一区精品| 超碰97精品在线观看| 久久99精品国语久久久| 久久久精品区二区三区| 午夜av观看不卡| 性色av一级| 亚洲美女黄色视频免费看| 国产精品一二三区在线看| 2022亚洲国产成人精品| 三级国产精品片| 日产精品乱码卡一卡2卡三| 一边摸一边做爽爽视频免费| 日产精品乱码卡一卡2卡三| 99热全是精品| 欧美日韩在线观看h| 欧美3d第一页| 国产 一区精品| 纯流量卡能插随身wifi吗| 3wmmmm亚洲av在线观看| 熟妇人妻不卡中文字幕| 美女国产视频在线观看| 亚洲少妇的诱惑av| 一级毛片 在线播放| 亚洲av成人精品一区久久| 黑人高潮一二区| 国产免费又黄又爽又色| 国产精品三级大全| 久久毛片免费看一区二区三区| 亚洲国产精品一区三区| 国产精品99久久久久久久久| av国产精品久久久久影院| 久久久国产精品麻豆| 毛片一级片免费看久久久久| 久久久久久久久久人人人人人人| 亚洲高清免费不卡视频| 国产有黄有色有爽视频| 热99国产精品久久久久久7| 纵有疾风起免费观看全集完整版| 熟女人妻精品中文字幕| 这个男人来自地球电影免费观看 | 免费播放大片免费观看视频在线观看| 久久精品国产亚洲av天美| 日韩制服骚丝袜av| 午夜影院在线不卡| 国产男女内射视频| 考比视频在线观看| 成人漫画全彩无遮挡| 一个人看视频在线观看www免费| 一级a做视频免费观看| 日韩欧美一区视频在线观看| 亚洲精品国产色婷婷电影| 欧美激情国产日韩精品一区| 一区二区三区四区激情视频| 国产精品一区二区在线观看99| 夜夜骑夜夜射夜夜干| 国产免费福利视频在线观看| 亚洲欧美清纯卡通| 交换朋友夫妻互换小说| 99热全是精品| 国产 一区精品| 久久精品国产亚洲av天美| 我的女老师完整版在线观看| 夜夜爽夜夜爽视频| av女优亚洲男人天堂| 你懂的网址亚洲精品在线观看| videossex国产| 九草在线视频观看| 大片免费播放器 马上看| 国产有黄有色有爽视频| 国产精品99久久99久久久不卡 | 在线观看人妻少妇| 有码 亚洲区| 日韩成人伦理影院| 大片免费播放器 马上看| 你懂的网址亚洲精品在线观看| 中文乱码字字幕精品一区二区三区| freevideosex欧美| 国产欧美日韩一区二区三区在线 | 久热这里只有精品99| 欧美性感艳星| 女性生殖器流出的白浆| 国产老妇伦熟女老妇高清| 国产黄色免费在线视频| 999精品在线视频| 人人妻人人澡人人看| 少妇精品久久久久久久| 日韩强制内射视频| 九草在线视频观看| 色婷婷久久久亚洲欧美| 这个男人来自地球电影免费观看 | 亚洲av福利一区| 久久毛片免费看一区二区三区| 天美传媒精品一区二区| 国产男女超爽视频在线观看| 水蜜桃什么品种好| 欧美人与性动交α欧美精品济南到 | 日本免费在线观看一区| 欧美精品亚洲一区二区| 免费观看a级毛片全部| 日日啪夜夜爽| 国国产精品蜜臀av免费| 国产成人freesex在线| 丰满少妇做爰视频| 亚洲av国产av综合av卡| 国产爽快片一区二区三区| 69精品国产乱码久久久| 嘟嘟电影网在线观看| 七月丁香在线播放| 人妻 亚洲 视频| 久久综合国产亚洲精品| 亚洲经典国产精华液单| 国产 精品1| 亚洲熟女精品中文字幕| 久久久a久久爽久久v久久| 欧美另类一区| av有码第一页| 亚洲欧美中文字幕日韩二区| 国产午夜精品久久久久久一区二区三区| 欧美成人午夜免费资源| 夫妻午夜视频| 人妻制服诱惑在线中文字幕| av不卡在线播放| 亚洲精品日韩在线中文字幕| 在线亚洲精品国产二区图片欧美 | a级片在线免费高清观看视频| 亚洲欧美日韩卡通动漫| 久久人妻熟女aⅴ| 国产亚洲一区二区精品| 精品久久久久久电影网| 热99久久久久精品小说推荐| 国产成人a∨麻豆精品| 亚洲av日韩在线播放| 超碰97精品在线观看| 久久久久久久久大av| 午夜视频国产福利| 亚洲高清免费不卡视频| 精品人妻偷拍中文字幕| 狠狠精品人妻久久久久久综合| 国产欧美另类精品又又久久亚洲欧美| 欧美激情 高清一区二区三区| 国产精品久久久久久av不卡| 男的添女的下面高潮视频| 婷婷成人精品国产| 色吧在线观看| 国产毛片在线视频| 午夜福利影视在线免费观看| 免费人妻精品一区二区三区视频| 十八禁高潮呻吟视频| 男女国产视频网站| 日韩一本色道免费dvd| 99热这里只有精品一区| √禁漫天堂资源中文www| 日韩熟女老妇一区二区性免费视频| 女性被躁到高潮视频| 国产精品无大码| 少妇精品久久久久久久| 久久人人爽av亚洲精品天堂| 免费观看在线日韩| 国产69精品久久久久777片| 七月丁香在线播放| 国产毛片在线视频| 免费观看的影片在线观看| 国产一区二区三区av在线| 91成人精品电影| 成人18禁高潮啪啪吃奶动态图 | 狂野欧美激情性xxxx在线观看| 看十八女毛片水多多多| 最近2019中文字幕mv第一页| 国产一区二区在线观看av| 赤兔流量卡办理| 国精品久久久久久国模美| 22中文网久久字幕| 亚洲精品国产av蜜桃| 午夜激情久久久久久久| 国产成人精品福利久久| 国产精品一区www在线观看| 亚洲美女视频黄频| 国产精品国产av在线观看| 精品卡一卡二卡四卡免费| 日韩av不卡免费在线播放| 九九爱精品视频在线观看| 男女高潮啪啪啪动态图| 老司机影院毛片| 午夜激情福利司机影院| 国产成人91sexporn| 国产精品不卡视频一区二区| 婷婷色综合www| 午夜av观看不卡| 一级毛片aaaaaa免费看小| 少妇猛男粗大的猛烈进出视频| av网站免费在线观看视频| 亚洲av在线观看美女高潮| 国产亚洲最大av| 久久久久久久久久人人人人人人| 日本黄色片子视频| xxx大片免费视频| 性高湖久久久久久久久免费观看| 欧美日韩av久久| 久久久久久久久大av| 色5月婷婷丁香| 精品一区二区三区视频在线| 男人操女人黄网站| 成人午夜精彩视频在线观看| 18禁在线无遮挡免费观看视频| 2022亚洲国产成人精品| 午夜影院在线不卡| 校园人妻丝袜中文字幕| 国产免费一级a男人的天堂| 欧美 日韩 精品 国产| 国产精品欧美亚洲77777| 欧美最新免费一区二区三区| 日本黄色片子视频| 欧美精品国产亚洲| .国产精品久久| 天堂俺去俺来也www色官网| 久久人人爽人人爽人人片va| 女性生殖器流出的白浆| 国产精品一二三区在线看| 亚洲国产欧美在线一区| 色94色欧美一区二区| 建设人人有责人人尽责人人享有的| 啦啦啦在线观看免费高清www| 看十八女毛片水多多多| 99热全是精品| 色吧在线观看| 五月开心婷婷网|