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

    西安市蔬菜種植區(qū)土壤屬性的空間變異與肥力適宜性

    2017-08-01 10:00:08楊聯(lián)安楊煜岑馮武煥于世鋒宋英強李新堯
    水土保持通報 2017年3期
    關鍵詞:土壤肥力西安市速效

    王 晶, 楊聯(lián)安, 楊煜岑, 馮武煥, 于世鋒, 宋英強, 張 彬, 李新堯

    (1.西北大學 城市與環(huán)境學院, 陜西 西安 710127;2.西安市農業(yè)技術推廣中心,陜西 西安 7100613; 3.西安市農產品質量安全檢驗監(jiān)測中心, 陜西 西安 710077)

    西安市蔬菜種植區(qū)土壤屬性的空間變異與肥力適宜性

    王 晶1, 楊聯(lián)安1, 楊煜岑1, 馮武煥2, 于世鋒3, 宋英強1, 張 彬1, 李新堯1

    (1.西北大學 城市與環(huán)境學院, 陜西 西安 710127;2.西安市農業(yè)技術推廣中心,陜西 西安 7100613; 3.西安市農產品質量安全檢驗監(jiān)測中心, 陜西 西安 710077)

    [目的] 對西安市蔬菜區(qū)土壤屬性空間變異與肥力適宜性進行評價,為相關部門土壤培肥和蔬菜種植提供決策依據(jù)。[方法] 選擇與蔬菜生長密切相關的土壤pH值、電導率、有機質、堿解氮、有效磷、速效鉀六個指標,利用半變異函數(shù)進行空間變異分析,應用模糊綜合法計算土壤肥力適宜性指數(shù)(SFI),對SFI進行Kriging插值,采用等間距法劃分等級,并對土壤肥力進行綜合評價。[結果] 半變異函數(shù)能較好的擬合蔬菜區(qū)土壤空間變異特征,人為施肥成為影響土壤電導率、堿解氮空間變異的主要因素,土壤pH值和其他養(yǎng)分則同時受到結構性因素和隨機性因素的作用。[結論] 西安市蔬菜區(qū)土壤肥力適宜性處于中等較好水平,限制蔬菜生長的主要因素是土壤有機質和堿解氮含量。臨潼區(qū)東南部,鄠邑區(qū)、周至縣北部地區(qū)以及長安區(qū)和鄠邑區(qū)北部交界處蔬菜區(qū)土壤肥力有待提高,菜農需要采取措施提高土壤有機質含量和補施氮肥。

    土壤肥力適宜性; 空間變異; 模糊數(shù)學; 蔬菜區(qū)

    文獻參數(shù): 王晶, 楊聯(lián)安, 楊煜岑, 等.西安市蔬菜種植區(qū)土壤屬性的空間變異與肥力適宜性[J].水土保持通報,2017,37(3):204-209.DOI:10.13961/j.cnki.stbctb.2017.03.034; Wang Jing, Yang Lianan, Yang Yucen, et al. Spatial variation and fertility suitability of vegetable field in Xi’an City[J]. Bulletin of Soil and Water Conservation, 2017,37(3):204-209.DOI:10.13961/j.cnki.stbctb.2017.03.034

    蔬菜產業(yè)是西安市種植業(yè)中最具競爭力的優(yōu)勢產業(yè)之一,對地方財政增長、農業(yè)增效、農民增收具有重要的現(xiàn)實意義。近年來,西安市大力推進蔬菜板塊建設和蔬菜產業(yè)帶建設,有力地推動了蔬菜產業(yè)快速發(fā)展。截止2014年底,西安市蔬菜播種面積達6.61×104hm2,年產蔬菜299×104t[1]。同時,蔬菜質量與安全越來越引起人們重視,為了獲得優(yōu)質、高產的蔬菜產品,就必須正確認識土壤肥力在蔬菜區(qū)環(huán)境中的重要意義。土壤肥力是蔬菜區(qū)土壤質量重要指標,是土壤的基本屬性和本質特征的反映,是土壤從營養(yǎng)條件與環(huán)境條件方面供應與協(xié)調作物生長能力。因此,通過開展蔬菜區(qū)土壤質量調查與綜合評價研究,對高效利用土地資源,提高土壤生產力,提升土壤養(yǎng)分合理管理水平具有重要意義。

    農業(yè)土壤屬性空間變異研究主要集中于土壤屬性空間分布特征與影響因素、土壤屬性空間變異特征及其綜合評價[2]、土壤重金屬污染分析[3]以及耕地土壤地力綜合評價[4]等,但對于蔬菜集中區(qū)(蔬菜區(qū))土壤研究較少,并且未與土壤肥力結合起來。前人[4-5]研究發(fā)現(xiàn)導致土壤屬性空間變異的因素包括地形、土壤母質、土壤類型、土壤質地等,種植類型、人為耕作管理措施等。土壤肥力定量綜合評價中常采用指數(shù)和法[6]、模糊綜合評價法[7]、層次分析法的物元模型[8]、修正的內梅羅法等[9]。本文選用模糊數(shù)學評價法,并與地統(tǒng)計學相結合,對西安市蔬菜區(qū)土壤pH值、電導率、有機質、堿解氮、有效磷、速效鉀6個指標進行土壤肥力適宜性綜合評價,旨在為西安市蔬菜區(qū)土壤養(yǎng)分管理和合理施肥決策提供科學依據(jù)。

    1 研究材料與數(shù)據(jù)采集

    1.1 研究區(qū)概況

    西安市位于黃河流域中部關中盆地,東經107°40′—109°49′和北緯33°42′—34°45′之間。轄境東西長約204 km,南北寬約116 km,面積10 108 km2,年平均氣溫13.0 ℃~ 13.7 ℃,年降水量522.4~719.5 mm,由北向南遞增。7月、9月為兩個明顯降水高峰月,年日照時數(shù)1 646.1~2 114.9 h。蔬菜區(qū)主要分布于該市東北部和中部,東北部蔬菜區(qū)包括灞橋區(qū)、高陵區(qū)、臨潼區(qū)、閻良區(qū)部分村鎮(zhèn),中部蔬菜區(qū)包括長安區(qū)、鄠邑區(qū)、周至縣部分村鎮(zhèn)。

    1.2 土壤樣品采集與分析

    以西安市蔬菜區(qū)作為試驗樣區(qū),為了提高估值精度,充分考慮不同土壤類型和地形特點,力求均勻布設采樣點。采樣點分布于灞橋區(qū)、高陵區(qū)、臨潼區(qū)、閻良區(qū)、長安區(qū)、鄠邑區(qū)、周至縣的菜田,并于2014年11月進行土壤樣品采集,所有樣點采集遵循“等量”、“隨機”、“多點混合”的原則,利用全球定位系統(tǒng)(GPS)對采樣點定位,記錄其經緯度位置,以“S”型或鋸齒形隨機采樣,以深度0—20 cm土層取10個樣點混合成1 kg土樣,共采取土樣369份。將土壤樣品自然風干、過篩后,參照《土壤分析技術規(guī)范》中相關方法進行了測定,土壤pH值采用玻璃電極法測定,電導率則通過電導法測定,采用重鉻酸鉀容量法(油浴加熱)測定有機質,采用堿解擴散法測定堿解氮,采用0.5 mol/L碳酸氫鈉浸提—鉬銻抗比色法測定有效磷,采用乙酸銨提取——火焰光度法測定速效鉀[10]。土壤質地、種植方式、種植年限等由西安市農業(yè)技術推廣中心提供。

    2 研究方法

    2.1 數(shù)據(jù)處理方法

    采用GS+和SPSS軟件對研究區(qū)土壤pH、電導率、有機質、堿解氮、有效磷、速效鉀空間變異進行研究。利用3倍標準差法篩選異常值,用最大、最小值替換異常值。由于空間變異分析的前提是:數(shù)據(jù)必須符合正態(tài)分布,否則會出現(xiàn)比例效應[11]。因此進行正態(tài)分布檢驗,并對非正態(tài)分布的變量進行轉換以滿足正態(tài)分布。利用SPSS18.0軟件對數(shù)據(jù)進行正態(tài)分布檢驗和描述性統(tǒng)計分析;GS+軟件分析蔬菜區(qū)各土壤養(yǎng)分空間變異函數(shù)、進行模型擬合;分析土壤養(yǎng)分空間變異特征。半變異函數(shù)模型在決定系數(shù)(R2)接近1,殘差(RSS)盡可能小的原則選取最優(yōu)模型。建立各指標隸屬度函數(shù),基于主成分分析計算不同指標權重值,采用模糊數(shù)學方法計算土壤肥力適宜性評價指數(shù)(soil fertility suitability index, SFI),利用ArcGIS10.3軟件和地統(tǒng)計模塊Kriging插值,得出土壤肥力適宜性空間分布。

    2.2 評價指標的隸屬度函數(shù)構建

    隸屬度函數(shù)F(x)是模糊數(shù)學中的一個概念,是一個值域在0~1之間的函數(shù),表示論域u中的元素x隸屬與u的模糊子集A的程度。F(x)越接近1,表示x屬于A的可能性越大。將隸屬度函數(shù)運用到土壤肥力的研究中,旨在判斷土壤養(yǎng)分所屬等級。土壤養(yǎng)分隸屬度函數(shù)越接近于1,表明土壤養(yǎng)分狀況越好[12]。

    2.3 蔬菜區(qū)土壤肥力適宜性分析

    蔬菜區(qū)土壤肥力適宜性采用模糊數(shù)學隸屬函數(shù)的數(shù)學模型進行評價。考慮到研究區(qū)土壤形成特點以及蔬菜生產性質,選擇土壤pH值、電導率、有機質、堿解氮、有效磷和速效鉀作為評價因子,并采用主成分分析法確定各因子權重。依據(jù)模糊數(shù)學原理構建蔬菜區(qū)土壤屬性和土壤養(yǎng)分指標的隸屬函數(shù),計算各土壤屬性和各土壤養(yǎng)分指標的隸屬度,再利用加成法計算研究區(qū)蔬菜區(qū)土壤肥力適宜性評價指數(shù)(ISF),計算公式為:

    (1)

    式中:Nij和Wij——第i個土壤樣品的第j個土壤屬性或土壤養(yǎng)分指標的隸屬度和相應的權重系數(shù)。ISF取值范圍為0~1,其值越高,表明土壤肥力對蔬菜的生長和品質的貢獻率越高。

    3 結果與分析

    3.1 蔬菜區(qū)土壤屬性基本特征

    利用SPSS軟件對經異常值處理的蔬菜區(qū)土壤屬性數(shù)據(jù)進行描述性統(tǒng)計和正態(tài)分布檢驗,統(tǒng)計結果見表1。表明蔬菜區(qū)土壤pH值、土壤有機質、堿解氮、速效鉀原數(shù)據(jù)符合正態(tài)分布。對電導率、有效磷進行了對數(shù)轉換,并對轉換后的數(shù)據(jù)進行正態(tài)檢驗,電導率、有效磷滿足對數(shù)分布。參考全國第二次土壤普查土壤養(yǎng)分分級標準,研究區(qū)蔬菜區(qū)土壤有機質含量屬于中等偏低水平,堿解氮含量處于中等水平,大部分采樣點有效磷、速效鉀含量豐富。

    表1 蔬菜區(qū)土壤指標含量統(tǒng)計性描述

    變異系數(shù)是指數(shù)據(jù)空間變異的大小,變異系數(shù)越大,說明它在空間上分布不均勻,受到外界影響程度越高,也就是說隨機性越強。一般認為,當變異系數(shù)≤10%時,為弱變異;10%<變異系數(shù)<100%時,表示中等變異;變異系數(shù)≥100%時,表示強變異。由表1可知電導率變異系數(shù)達到強變異,說明蔬菜區(qū)土壤電導率分布不均勻是由人為因素造成的。有機質、堿解氮、有效磷、速效鉀都屬于中等變異,且變異系數(shù)由低到高依次為堿解氮、有機質、速效鉀、有效磷,其中有效磷含量豐富且變異系數(shù)大,表明該養(yǎng)分含量與施用磷肥量有關。而蔬菜區(qū)土壤pH值變異系數(shù)只有4.76%,表明土壤酸堿性主要受到成土母質等自然因素的作用。

    3.2 蔬菜區(qū)土壤屬性空間變異分析

    采用GS+7.0軟件對滿足正態(tài)分布的采樣點數(shù)據(jù)進行半變異函數(shù)模型擬合,通過調節(jié)不同模型參數(shù)對比模型的擬合度和殘差[13],得到各土壤養(yǎng)分半變異函數(shù)理論模型及參數(shù)(如表2所示)。

    表2 蔬菜區(qū)土壤指標理論模型及參數(shù)

    塊金值表示土壤特性受人為干擾隨機性因素引起的變異,這四種土壤養(yǎng)分的塊金值都大于零說明在采樣的過程中存在人為因素引起的隨機變異,可能存在于采集或者土壤養(yǎng)分測定過程中。基臺值表示土壤特性的總變異,包括結構性變異和隨機性變異。塊金效應即塊金值與基臺值的比值,用它來反應土壤養(yǎng)分的空間依賴性,塊金效應值越高說明隨機因素是導致土壤養(yǎng)分空間變異的主要因素[5]。pH值、有機質、有效磷、速效鉀的塊金效應都介于25%~75%之間,屬于中等變異,說明引起土壤pH、有機質、有效磷、速效鉀空間變異的主導因素為結構性因素,比如土壤母質性質、氣候、地形等,而電導率、堿解氮的塊金效應分別達76.4%和88.4%,蔬菜區(qū)土壤電導率和堿解氮主要受隨機因素影響,表明該地區(qū)蔬菜養(yǎng)分的空間變異主要與菜農的種植制度和土壤培肥有關。變程表示的是土壤養(yǎng)分特性空間自相關范圍[14],表2反映電導率的變程距離較大,說明電導率空間變異以大塊變異為主,而其他屬性的變程較小,空間變異以小塊變異為主。

    3.3 蔬菜區(qū)土壤肥力適宜性評價

    3.3.1 評價指標的選取與權重確定 選取評價土壤肥力適宜性指標應遵循主導性、系統(tǒng)性、獨立性、空間變異性和指標體系區(qū)域性等原則。根據(jù)前人研究結果和蔬菜區(qū)土壤屬性、養(yǎng)分分布特點,選取對蔬菜生長發(fā)育影響較大的pH值、有機質、堿解氮、有效磷和速效鉀作為評價指標,采用模糊綜合評價法確定土壤肥力適宜性。主成分分析法實現(xiàn)了定性分析和定量分析相結合,解決了評價模型中確定權重的問題。因此,采用主成分分析法來確定權重[15],采用統(tǒng)計軟件SPSS進行分析,得出各個指標主成分的特征值和貢獻率,通過SPSS計算各個公因子方差占公因子方差總和的百分率作為單項評價指標的權重(表3)。其值的大小表示該項肥力指標對土壤肥力總體變異的貢獻。由表3可知:pH值的權重最低,這主要是由于研究區(qū)土壤受自然因素的作用較大,pH值屬于偏堿性土壤,變化較小,而且其本身對植物生長不起直接作用,主要通過酸堿性來調節(jié)植物對其他養(yǎng)分的吸收。除pH值之外,其他指標權重都介于0.1~0.3之間。

    表3 蔬菜區(qū)土壤評價指標的特征值、貢獻率和權重

    3.3.2 評價指標的隸屬函數(shù)構建及轉折點取值 計算不同評價指標的隸屬度值,首先要建立各評價指標之間相應的隸屬度函數(shù)來表示各評價指標的狀態(tài)值,同時實現(xiàn)評價指標量綱歸一化,使評價指標間具有可比性。由于土壤因子變化具有連續(xù)性,故各評價指標采用連續(xù)性質的隸屬度函數(shù)。結合王建國等[16]提出的模糊數(shù)學土壤肥力綜合評價模型,并根據(jù)研究區(qū)實際情況,采用以下3種隸屬度函數(shù)。pH值符合拋物線型,電導率為降半梯型,有機質、堿解氮、有效磷、速效鉀為S型,為了方便計算,將曲線函數(shù)轉換為折線函數(shù)。

    (2)

    (3)

    (4)

    結合蔬菜地土壤屬性含量實際情況和專家經驗,分別確定各指標隸屬度函數(shù)曲線轉折點取值(表4)[4]。

    表4 蔬菜區(qū)土壤評價指標的隸屬度函數(shù)及函數(shù)轉折點取值

    3.3.3 土壤屬性單因子適宜性評價 依據(jù)隸屬度函數(shù)計算蔬菜區(qū)各指標平均隸屬度,并分別與權重值相乘得單因子適宜性(表5)。平均隸屬度指某一指標對于蔬菜生長的平均優(yōu)劣度,平均隸屬度越高,說明該指標在土壤中的狀態(tài)越好,蔬菜越適宜該環(huán)境。由表5可知,平均隸屬度由高到底依次為:pH值>電導率>速效鉀>有效磷>堿解氮>有機質。有機質和堿解氮平均隸屬度都低于0.6,說明限制蔬菜生長的主要因素就是土壤有機質和堿解氮含量,這與土壤養(yǎng)分豐缺統(tǒng)計結果一致。隸屬度高的指標不一定對土壤肥力的貢獻最大,要同時考慮該因子在土壤肥力中的權重值。單因子肥力適宜性綜合隸屬度與權重,可分析各因子對綜合肥力的貢獻率。單因子適宜性得分越高,表明其對蔬菜生長貢獻率越高。由此可得,6個指標值貢獻率由高到低依次為:電導率>有效磷>堿解氮>有機質>速效鉀>pH值。pH值、有機質和速效鉀對土壤肥力的貢獻值都低于0.1,pH值對土壤肥力的貢獻值較低是受到其權重值得制約,pH值平均隸屬度指最好,表明菜田土壤pH值對蔬菜生長具有積極影響。速效鉀對土壤肥力的貢獻率較低的原因是速效鉀含量已經達到豐富水平,因此,建議菜農可以減少鉀肥的施用,甚至當季可以不施鉀肥。

    表5 蔬菜區(qū)評價指標隸屬度平均值及單因子適宜性

    3.3.4 土壤肥力適宜性綜合評價 土壤肥力適宜性受到各指標的共同影響,某采樣點土壤單因子適宜性值高并不能表示該地綜合肥力值也高,綜合肥力值可能受到其他因子影響而有所降低。利用公式(1)和計算得到的蔬菜區(qū)各指標權重、隸屬度計算得蔬菜區(qū)土壤肥力適宜性評價指數(shù)(ISF),ISF取值范圍為0~1,值越大,表明土壤肥力水平越高。結果表明蔬菜區(qū)SFI最小值為0.36,最大值為0.95,平均值為0.65,標準差為0.12,變異系數(shù)為18.43%。

    表6為西安市蔬菜區(qū)土壤養(yǎng)分適宜性指數(shù)的半變異函數(shù)擬合模型。依據(jù)半變異函數(shù)特征值表明蔬菜區(qū)土壤適宜性屬于中等變異程度,且自變異距離較小。

    表6 土壤肥力適宜性理論模型及參數(shù)

    利用ArcGIS 10.3軟件基于如上土壤肥力適宜性半變異模型,采用普通克里金插值方法[17],可以分析西安市蔬菜區(qū)土壤肥力適宜性空間分布格局,如圖1所示。經過插值的蔬菜區(qū)土壤肥力適宜性平滑掉小部分的樣點,依據(jù)等間距法將蔬菜區(qū)土壤肥力適宜性進行分級統(tǒng)計,結果如表7所示。

    圖1 蔬菜區(qū)土壤肥力適宜性分布表7 土壤肥力適宜性各等級面積比例統(tǒng)計

    等級SFI范圍面積比/%五級0.49~0.556.36四級0.56~0.6119.59三級0.62~0.6737.53二級0.68~0.7333.37一級0.74~0.803.14

    由表7可知,西安市蔬菜區(qū)SFI等級占面積比由大到小依次是:三級37.53%,二級33.37%,四級19.59%,五級6.36%,一級3.14%??傮w上說西安市蔬菜區(qū)肥力適宜性處于中等較好水平。圖1表明西安市蔬菜區(qū)ISF分布情況:長安區(qū)西北角、東南角,以及臨潼區(qū)中東部分布有少量ISF一級土壤;ISF二級土壤主要分布于閻良區(qū)和臨潼區(qū)交界處、高陵區(qū)中西部、長安區(qū)西部;ISF處于三級以下的土壤主要分布于臨潼區(qū)東南部、鄠邑區(qū)、周至縣北部地區(qū)以及長安區(qū)與鄠邑區(qū)北部交界處;其他廣闊地區(qū)土壤都處于ISF三級水平。ISF三級以下地區(qū)應重視蔬菜區(qū)土壤培肥。

    西安市蔬菜區(qū)ISF一、二級土壤各項指標狀況均處于良好水平,雖然閻良區(qū)、高陵區(qū)蔬菜種植年限較長,土壤自身肥力降低,但多采用溫室大棚的種植方式,水熱條件好,菜農施肥比較合理,因此土壤肥力適宜性好。而長安區(qū)蔬菜種植年限較短,且土壤質地為粉砂質黏壤土,有利于土壤有機質積累,土壤自身肥力條件較好,雖多為露地種植,但總體土壤肥力適應性較好。西安市蔬菜區(qū)ISF四、五級土壤與一、二級土壤相比,土壤有機質含量差異不大,相比差異較大的是有效磷和速效鉀含量,平均值含量相差100 mg/kg,因此,雖然整體上西安市蔬菜區(qū)有效磷、速效鉀含量豐富,但是各地分布不均勻,對于ISF四級、五級地區(qū)仍然要補施磷肥和鉀肥。臨潼區(qū)、鄠邑區(qū)北部ISF較低的土壤質地為壤土和中壤、沙壤,不利于土壤有機質積累,種植年限較長,土壤自身肥力較差,加上露地種植,所以土壤肥力適宜性較差。

    4 結 論

    (1) 土壤電導率、堿解氮受到隨機性因素影響較大,人為施肥成為影響其空間變異的主要因素。pH值、有機質、有效磷、速效鉀空間變異同時受到結構性因素(土壤母質、地形、氣候等自然因素)和隨機性因素共同作用。

    (2) 西安市蔬菜區(qū)土壤肥力適宜性處于中等較好水平,呈中等變異程度。臨潼區(qū)的東南部,鄠邑區(qū)、周至縣的北部地區(qū)以及長安區(qū)和鄠邑區(qū)北部交界處蔬菜區(qū)土壤肥力有待提高,應根據(jù)土壤肥力適宜性進行精準施肥決策,提高肥料利用率。

    (3) 西安市蔬菜區(qū)土壤肥力適宜性主要限制性因素是土壤有機質和堿解氮含量,西安市菜農存在氮肥表施、施肥盲目的問題。因此,菜農應重視土壤有機質含量的提升、積極制定蔬菜“補氮”措施,提升土壤肥力,保障蔬菜區(qū)蔬菜的品質,繼續(xù)發(fā)揮西安市蔬菜產業(yè)優(yōu)勢。

    [1] 陳琳.西安蔬菜產銷發(fā)展現(xiàn)狀的調查和思考[J].陜西農業(yè)科學,2015,61(11):105-107.

    [2] 張彬,楊聯(lián)安,楊粉莉,等.禮泉縣蘋果園土壤養(yǎng)分空間變異特征及綜合評價[J].土壤通報,2016,47(4):860-867.

    [3] 李雪芳,王文巖,上官宇先,等.西安市郊菜地土壤重金屬污染及其與蔬菜重金屬質量分數(shù)的相關性[J].西北農業(yè)學報,2014,23(8):173-181.

    [4] 易秀,谷曉靜,侯燕卿,等.陜西省涇惠渠灌區(qū)土壤肥力質量綜合評價[J].干旱區(qū)資源與環(huán)境,2011,25(2):132-137.

    [5] 胡艷霞,周連第,魏長山,等.水源保護地土壤養(yǎng)分空間變異特征及其影響因素分析[J].土壤通報,2013,44(5):1184-1191.

    [6] 趙明松,張甘霖,王德彩,等.徐淮黃泛平原土壤有機質空間變異特征及主控因素分析[J].土壤學報,2013,50(1):1-11.

    [7] 王軍艷,張鳳榮,王茹,等.應用指數(shù)和法對潮土農田土壤肥力變化的評價研究[J].生態(tài)與農村環(huán)境學報,2001,17(3):13-16.

    [8] 劉東海,梁國慶,周衛(wèi),等.基于神經網(wǎng)絡的土壤肥力綜合評價[J].中國土壤與肥料,2011,(5):12-19.

    [9] 趙串串,安若蘭,趙巧玉,等.物元模型在玉樹地區(qū)林地土壤養(yǎng)分評價中的應用[J].土壤,2014,46(6):984-988.

    [10] 趙映翠,耿增超,陳金海,等.宜君縣核桃經濟林地土壤養(yǎng)分調查與評價[J].干旱地區(qū)農業(yè)研究,2015,33(3):85-86.

    [11] 趙華富,周國蘭,劉曉霞,等.貴州茶區(qū)土壤養(yǎng)分狀況綜合評價[J].中國土壤與肥料,2012(3):30-34.

    [12] 史舟,李艷.地統(tǒng)計學在土壤學中的應用[M]:北京:中國農業(yè)出版社,2006.

    [13] 付同剛,張偉,陳洪松,等.喀斯特峰叢洼地土壤綜合肥力的空間分異特征及其影響因素[J].土壤通報,2013,44(2):385-391.

    [14] 王軍,傅伯杰,邱揚,等.黃土丘陵小流域土壤水分的時空變異特征:半變異函數(shù)[J].地理學報,2000,55(4):428-438.

    [15] 趙明松,張甘霖,吳運金,等.江蘇省土壤有機質含量時空變異特征及驅動力研究[J].土壤學報,2014,51(3):448-458.

    [16] 王建國,楊林章,單艷紅.模糊數(shù)學在土壤質量評價中的應用研究[J].土壤學報,2001,38(2):176-183.

    [17] 張彬,楊聯(lián)安,楊粉莉,等.蘋果主產區(qū)土壤養(yǎng)分空間分布特征及其影響因素[J].土壤,2016,48(4):777-784.

    Spatial Variation and Fertility Suitability of Vegetable Field in Xi’an City

    WANG Jing1, YANG Lianan1, YANG Yucen1, FENG Wuhuan2,YU Shifeng3, SONG Yingqiang1, ZHANG Bin1, LI Xinyao1

    (1.CollegeofUrbanandEnvironmentalSciences,NorthwestUniversity,Xi’an,Shaanxi710127,China; 2.Xi’anAgriculturalTechnologyExtensionCenter,Xian,Shaanxi710066,China;3.Xi'anAgriculturalProductQualitySafetyInspectionandMonitoringCenter,Xi’an,Shaanxi710077,China)

    [Objective] To evaluate the spatial variation and fertility suitability of vegetable field in Xi'an City, Shaanxi Province, in order to provide a basis for making decision of soil cultivation and agricultural production in relevant departments. [Methods] Six indexes that are closely related to vegetable growth including soil pH value, electrical conductivity, organic matter, alkali-hydrolysable nitrogen, available phosphorus and available potassium, were measured. Their spatial variability was analyzed using semi-variation function. And the soil fertility suitability was determined using fuzzy mathematics method. The soil fertility suitability was graded by equal spacing method after Kriging interpolation was conducted, and finally it was evaluated comprehensively. [Results] The semi-variogram could fit the spatial variability of soil in the vegetable field. The spatial distributions of the conductivity and alkali-hydrolysable nitrogen were mainly affected by the fertilization, and the soil pH value and other nutrient were affected by both anthropogenic and natural factors. [Conclusion] The soil fertility suitability of vegetable field in Xi’an City was at middle level, and the main limiting factors of vegetable growth were soil organic matter and alkali-hydrolysable nitrogen content. There are some areas where the soil fertility of vegetable field should be improved. These areas are the south-eastern of the Lintong District, the northern of the Huyi District and the Zhouzhi County, the northern part next to Chang’an District and the Huyi District. To the vegetable farmer, measures that can improve the contents of soil organic matter and nitrogen should be taken.

    soil fertility suitability; spatial variability; fuzzy mathematics; vegetable field

    2017-01-30

    2017-02-27

    教育部人文社會科學研究規(guī)劃基金項目“基于地統(tǒng)計學和 GIS的土壤網(wǎng)格采樣技術研究”(10YJA910010); 陜西省農業(yè)科技攻關項目(2011K02-11); 西安市科技計劃項目(NC1402,NC150201); 西北大學“211 工程”研究生自主創(chuàng)新項目(YZZ15013)

    王晶(1992—),女(漢族),內蒙古烏海市人,碩士研究生,研究方向為農業(yè)資源利用。E-mail:wangjingxdch@163.com。

    楊聯(lián)安(1968—),男(漢族),陜西省武功縣人,博士,副教授,主要從事地理信息系統(tǒng)開發(fā)與應用研究。E-mail:yanglianan@163.com。

    A

    1000-288X(2017)03-0204-06

    S158.2

    猜你喜歡
    土壤肥力西安市速效
    親子創(chuàng)意美工展
    西安市第四醫(yī)院
    123的幾種說法
    拉薩市土壤肥力變化趨勢與改良策略
    不同殺螨劑對柑桔紅蜘蛛速效與持效性防治田間試驗
    基于CA-Markov模型的西安市熱環(huán)境模擬研究
    土壤速效鉀測定方法的影響因素
    綠色科技(2016年16期)2016-10-11 06:53:51
    施用不同有機物料對山地果園紅壤速效氮磷鉀含量的影響
    豫中植煙區(qū)土壤速效養(yǎng)分的空間變異分析
    安吉白茶園土壤肥力現(xiàn)狀分析
    茶葉(2014年4期)2014-02-27 07:05:15
    久久久精品欧美日韩精品| 性色av乱码一区二区三区2| 国产一区二区三区在线臀色熟女| 女警被强在线播放| 成人18禁在线播放| 成人国产综合亚洲| 99国产极品粉嫩在线观看| 亚洲,欧美精品.| 免费搜索国产男女视频| 亚洲一区高清亚洲精品| 一区二区三区国产精品乱码| 伊人久久大香线蕉亚洲五| 欧美黑人欧美精品刺激| a在线观看视频网站| 美女扒开内裤让男人捅视频| 欧美国产精品va在线观看不卡| 日韩欧美 国产精品| 男女做爰动态图高潮gif福利片| 亚洲七黄色美女视频| 日本一区二区免费在线视频| 热re99久久国产66热| 丝袜在线中文字幕| 长腿黑丝高跟| 成人手机av| 日韩欧美在线二视频| 少妇熟女aⅴ在线视频| 久久人人精品亚洲av| 亚洲av电影在线进入| 欧美日韩福利视频一区二区| 精品一区二区三区av网在线观看| 老汉色∧v一级毛片| 亚洲,欧美精品.| 99riav亚洲国产免费| 日韩欧美在线二视频| 手机成人av网站| www日本黄色视频网| 香蕉丝袜av| 99久久精品国产亚洲精品| 亚洲中文日韩欧美视频| 97碰自拍视频| 成人免费观看视频高清| 人人妻人人澡人人看| 日日夜夜操网爽| 国产99白浆流出| 亚洲成人国产一区在线观看| 又大又爽又粗| www.熟女人妻精品国产| 亚洲熟女毛片儿| 午夜免费观看网址| 制服诱惑二区| 人妻丰满熟妇av一区二区三区| 成人特级黄色片久久久久久久| 久久亚洲精品不卡| 亚洲五月婷婷丁香| 精品久久久久久久久久久久久 | 国产黄a三级三级三级人| tocl精华| 高潮久久久久久久久久久不卡| 亚洲精品中文字幕在线视频| 国产精品日韩av在线免费观看| av福利片在线| 成人欧美大片| 亚洲 国产 在线| 亚洲无线在线观看| 首页视频小说图片口味搜索| 亚洲自偷自拍图片 自拍| 亚洲欧美日韩无卡精品| 老鸭窝网址在线观看| www.精华液| 国产乱人伦免费视频| 国内揄拍国产精品人妻在线 | 搡老熟女国产l中国老女人| 久久欧美精品欧美久久欧美| 免费看美女性在线毛片视频| www.精华液| 亚洲专区国产一区二区| 久久久水蜜桃国产精品网| 久9热在线精品视频| 精品国产乱子伦一区二区三区| 成年版毛片免费区| 一区二区三区高清视频在线| 久久久水蜜桃国产精品网| 免费在线观看黄色视频的| 日本三级黄在线观看| 看免费av毛片| 黑人巨大精品欧美一区二区mp4| 久久婷婷成人综合色麻豆| 天堂动漫精品| 别揉我奶头~嗯~啊~动态视频| 日日摸夜夜添夜夜添小说| 亚洲专区中文字幕在线| 欧美日韩亚洲综合一区二区三区_| 国产区一区二久久| 在线十欧美十亚洲十日本专区| 国产在线观看jvid| 日韩大码丰满熟妇| 国产精品野战在线观看| 免费看a级黄色片| 欧美激情高清一区二区三区| 伦理电影免费视频| 老熟妇仑乱视频hdxx| 久久久久久国产a免费观看| 国内精品久久久久久久电影| 757午夜福利合集在线观看| 一a级毛片在线观看| 亚洲黑人精品在线| 精品一区二区三区av网在线观看| 亚洲精品一区av在线观看| 免费观看人在逋| 亚洲av电影不卡..在线观看| 亚洲 欧美 日韩 在线 免费| 欧美日韩中文字幕国产精品一区二区三区| 少妇的丰满在线观看| 天天一区二区日本电影三级| 亚洲激情在线av| 禁无遮挡网站| 可以在线观看毛片的网站| 日韩欧美国产一区二区入口| 亚洲欧美一区二区三区黑人| 精品第一国产精品| 999久久久精品免费观看国产| 老汉色av国产亚洲站长工具| 日韩有码中文字幕| 美女国产高潮福利片在线看| 长腿黑丝高跟| 亚洲成a人片在线一区二区| 国内精品久久久久精免费| 侵犯人妻中文字幕一二三四区| 老司机在亚洲福利影院| 一边摸一边做爽爽视频免费| 两性午夜刺激爽爽歪歪视频在线观看 | 一个人观看的视频www高清免费观看 | 99在线人妻在线中文字幕| 天堂√8在线中文| 18禁国产床啪视频网站| 视频区欧美日本亚洲| 久久亚洲真实| 国产成人av教育| 久久精品国产99精品国产亚洲性色| 在线av久久热| 免费在线观看黄色视频的| 欧美国产日韩亚洲一区| 国产aⅴ精品一区二区三区波| 日日爽夜夜爽网站| 黄片播放在线免费| netflix在线观看网站| 又黄又粗又硬又大视频| 757午夜福利合集在线观看| 男人舔奶头视频| 久久狼人影院| 久久天堂一区二区三区四区| 亚洲成av片中文字幕在线观看| 亚洲av电影在线进入| 久久性视频一级片| 十分钟在线观看高清视频www| 久久精品国产99精品国产亚洲性色| 亚洲av电影不卡..在线观看| svipshipincom国产片| 精品久久久久久久久久免费视频| 男男h啪啪无遮挡| 又黄又爽又免费观看的视频| 国产视频一区二区在线看| 精品国产超薄肉色丝袜足j| 欧美+亚洲+日韩+国产| av在线播放免费不卡| 免费观看人在逋| 国产精品美女特级片免费视频播放器 | 欧美激情极品国产一区二区三区| 日韩欧美三级三区| 国产精品av久久久久免费| 亚洲国产欧美一区二区综合| 国产精品亚洲av一区麻豆| 婷婷精品国产亚洲av| 亚洲国产欧洲综合997久久, | 一a级毛片在线观看| 午夜福利视频1000在线观看| 欧美性长视频在线观看| 丰满的人妻完整版| 亚洲免费av在线视频| 国产黄片美女视频| 哪里可以看免费的av片| √禁漫天堂资源中文www| 999精品在线视频| 成在线人永久免费视频| 精品久久久久久久久久久久久 | 国产精华一区二区三区| 天天一区二区日本电影三级| 久久香蕉国产精品| 国产单亲对白刺激| 99久久久亚洲精品蜜臀av| 亚洲精品国产精品久久久不卡| 亚洲av成人av| 身体一侧抽搐| 黄片小视频在线播放| 成人手机av| 免费在线观看日本一区| 亚洲精品一卡2卡三卡4卡5卡| 精品久久久久久,| 欧美午夜高清在线| 制服诱惑二区| 一级毛片女人18水好多| 一进一出抽搐gif免费好疼| 长腿黑丝高跟| aaaaa片日本免费| 欧美 亚洲 国产 日韩一| 国产精品久久久av美女十八| 亚洲精华国产精华精| 欧洲精品卡2卡3卡4卡5卡区| 久久久久国产精品人妻aⅴ院| 一级a爱片免费观看的视频| 51午夜福利影视在线观看| 久久草成人影院| 99国产精品一区二区三区| 午夜福利18| 免费观看人在逋| 老司机午夜福利在线观看视频| 日本熟妇午夜| 成人三级黄色视频| 亚洲欧美一区二区三区黑人| 特大巨黑吊av在线直播 | 长腿黑丝高跟| 免费在线观看黄色视频的| 波多野结衣巨乳人妻| 神马国产精品三级电影在线观看 | 亚洲一码二码三码区别大吗| 色av中文字幕| 久久久久国产精品人妻aⅴ院| 国产又色又爽无遮挡免费看| 男女视频在线观看网站免费 | 国产精品野战在线观看| 亚洲 国产 在线| 制服丝袜大香蕉在线| 中国美女看黄片| 国产人伦9x9x在线观看| 中文亚洲av片在线观看爽| 免费在线观看完整版高清| 青草久久国产| 亚洲av中文字字幕乱码综合 | xxxwww97欧美| 91九色精品人成在线观看| 国产真人三级小视频在线观看| 国产在线观看jvid| 天天添夜夜摸| 亚洲精品国产一区二区精华液| 妹子高潮喷水视频| 国产三级黄色录像| xxxwww97欧美| 99国产精品一区二区三区| 久久国产精品男人的天堂亚洲| 91成年电影在线观看| 桃红色精品国产亚洲av| 99热这里只有精品一区 | 久久狼人影院| 欧美性猛交╳xxx乱大交人| 亚洲aⅴ乱码一区二区在线播放 | 伦理电影免费视频| 国产男靠女视频免费网站| 桃红色精品国产亚洲av| 变态另类成人亚洲欧美熟女| 久久中文字幕一级| 国产1区2区3区精品| 黄色女人牲交| 国产精品亚洲av一区麻豆| 亚洲国产欧美日韩在线播放| 视频在线观看一区二区三区| 亚洲国产精品合色在线| 欧美日韩亚洲综合一区二区三区_| 夜夜躁狠狠躁天天躁| 日韩成人在线观看一区二区三区| 无人区码免费观看不卡| 亚洲国产欧美日韩在线播放| 两性夫妻黄色片| 婷婷亚洲欧美| 久久人妻福利社区极品人妻图片| 看黄色毛片网站| 免费在线观看完整版高清| 色av中文字幕| 久久人人精品亚洲av| 午夜精品久久久久久毛片777| 国产精品99久久99久久久不卡| 久久中文字幕人妻熟女| 日日摸夜夜添夜夜添小说| 宅男免费午夜| 日韩欧美在线二视频| 人人妻人人澡欧美一区二区| 好男人在线观看高清免费视频 | 亚洲成人久久爱视频| 亚洲av五月六月丁香网| 一夜夜www| 久久国产亚洲av麻豆专区| 91麻豆精品激情在线观看国产| 欧美黄色淫秽网站| 国产精品久久久久久亚洲av鲁大| 国产成人av教育| 999久久久精品免费观看国产| 69av精品久久久久久| 成人精品一区二区免费| 99在线人妻在线中文字幕| 中文在线观看免费www的网站 | 一级片免费观看大全| 香蕉久久夜色| 国产一区二区在线av高清观看| 九色国产91popny在线| 日韩 欧美 亚洲 中文字幕| 亚洲成a人片在线一区二区| 一进一出抽搐gif免费好疼| 一本久久中文字幕| 侵犯人妻中文字幕一二三四区| av视频在线观看入口| 亚洲国产日韩欧美精品在线观看 | 免费在线观看视频国产中文字幕亚洲| 俄罗斯特黄特色一大片| 亚洲av成人一区二区三| 欧美性猛交黑人性爽| 欧美一级a爱片免费观看看 | 精品熟女少妇八av免费久了| 久久婷婷成人综合色麻豆| 俺也久久电影网| 成人国产综合亚洲| 色播亚洲综合网| 欧美绝顶高潮抽搐喷水| 日本免费a在线| 看黄色毛片网站| 久久久水蜜桃国产精品网| 成人18禁高潮啪啪吃奶动态图| 久久精品aⅴ一区二区三区四区| 女人被狂操c到高潮| 嫩草影视91久久| 在线观看免费午夜福利视频| 999久久久国产精品视频| 最近最新免费中文字幕在线| 日日夜夜操网爽| 国产av在哪里看| 色播亚洲综合网| 中文字幕久久专区| 亚洲 国产 在线| 在线天堂中文资源库| 操出白浆在线播放| 女人爽到高潮嗷嗷叫在线视频| 99re在线观看精品视频| 亚洲国产精品成人综合色| 国产成人啪精品午夜网站| 老熟妇乱子伦视频在线观看| 久久精品国产清高在天天线| 国产亚洲av嫩草精品影院| 男男h啪啪无遮挡| 久久人人精品亚洲av| 在线观看午夜福利视频| 可以在线观看毛片的网站| 亚洲av成人不卡在线观看播放网| 一夜夜www| 在线播放国产精品三级| 身体一侧抽搐| 久久中文字幕一级| 窝窝影院91人妻| 亚洲久久久国产精品| 两个人看的免费小视频| 日本精品一区二区三区蜜桃| 怎么达到女性高潮| 黄色视频不卡| 中文字幕高清在线视频| 黄色视频,在线免费观看| 亚洲av成人一区二区三| 中文字幕另类日韩欧美亚洲嫩草| 国产免费男女视频| 满18在线观看网站| 亚洲欧美精品综合一区二区三区| 一边摸一边抽搐一进一小说| 亚洲精品在线美女| 国产欧美日韩一区二区精品| 成人国产一区最新在线观看| 老司机午夜十八禁免费视频| 色在线成人网| 成熟少妇高潮喷水视频| 日韩欧美一区二区三区在线观看| 最近最新中文字幕大全电影3 | 可以免费在线观看a视频的电影网站| 亚洲国产毛片av蜜桃av| 成熟少妇高潮喷水视频| 国产蜜桃级精品一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 精品国产美女av久久久久小说| 国产人伦9x9x在线观看| 亚洲第一青青草原| 成人一区二区视频在线观看| 一个人免费在线观看的高清视频| 中文字幕精品免费在线观看视频| 色播亚洲综合网| 国产精品一区二区精品视频观看| 久久精品成人免费网站| 在线观看www视频免费| 制服诱惑二区| 伊人久久大香线蕉亚洲五| 亚洲人成伊人成综合网2020| 丁香欧美五月| 国产精品永久免费网站| 久久中文看片网| 精品一区二区三区av网在线观看| 哪里可以看免费的av片| 9191精品国产免费久久| 日日干狠狠操夜夜爽| 亚洲色图av天堂| 日本 欧美在线| 精品不卡国产一区二区三区| 欧美成人性av电影在线观看| 美国免费a级毛片| 亚洲欧美一区二区三区黑人| 久久国产精品影院| 久久天躁狠狠躁夜夜2o2o| 欧美成人午夜精品| 麻豆国产av国片精品| 午夜视频精品福利| 日韩一卡2卡3卡4卡2021年| 国产三级在线视频| 欧美日韩亚洲国产一区二区在线观看| 男女下面进入的视频免费午夜 | 国产亚洲精品第一综合不卡| 一卡2卡三卡四卡精品乱码亚洲| 久久婷婷人人爽人人干人人爱| 色av中文字幕| 在线永久观看黄色视频| 啦啦啦观看免费观看视频高清| 久久久久久久午夜电影| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲av熟女| 国产精品98久久久久久宅男小说| 欧美乱妇无乱码| 99久久久亚洲精品蜜臀av| 久久久精品国产亚洲av高清涩受| 99久久无色码亚洲精品果冻| 又紧又爽又黄一区二区| 99精品久久久久人妻精品| 亚洲精品在线观看二区| 精品久久蜜臀av无| 国产精品久久久久久精品电影 | 午夜老司机福利片| 国产三级在线视频| 欧美中文综合在线视频| 免费在线观看完整版高清| 可以在线观看毛片的网站| 听说在线观看完整版免费高清| 欧美性长视频在线观看| 又大又爽又粗| 999精品在线视频| 精品无人区乱码1区二区| 美女免费视频网站| 成人18禁高潮啪啪吃奶动态图| 变态另类丝袜制服| 国产伦在线观看视频一区| 18禁美女被吸乳视频| 黄色成人免费大全| 欧美人与性动交α欧美精品济南到| 亚洲 国产 在线| 少妇的丰满在线观看| 国产av一区二区精品久久| 久久人妻av系列| 国产三级黄色录像| 在线十欧美十亚洲十日本专区| 久久久久国内视频| 久久精品人妻少妇| 亚洲第一欧美日韩一区二区三区| 国产欧美日韩一区二区三| 成人一区二区视频在线观看| 午夜亚洲福利在线播放| 久久香蕉精品热| 日韩有码中文字幕| 国产又爽黄色视频| 精品久久久久久久久久免费视频| 天天添夜夜摸| 亚洲专区国产一区二区| 欧美日韩福利视频一区二区| 搡老岳熟女国产| 最新美女视频免费是黄的| 成人免费观看视频高清| 亚洲va日本ⅴa欧美va伊人久久| 亚洲五月天丁香| 亚洲精品中文字幕在线视频| 久久久久久亚洲精品国产蜜桃av| 免费av毛片视频| 日日摸夜夜添夜夜添小说| 日本免费a在线| 精品不卡国产一区二区三区| 无人区码免费观看不卡| 国内揄拍国产精品人妻在线 | 国产免费av片在线观看野外av| 看免费av毛片| 亚洲七黄色美女视频| 国产亚洲av高清不卡| 亚洲人成电影免费在线| 欧美在线黄色| 日韩大尺度精品在线看网址| 国产黄色小视频在线观看| 久久香蕉精品热| 久久久精品国产亚洲av高清涩受| 国产高清激情床上av| 啦啦啦观看免费观看视频高清| 国内毛片毛片毛片毛片毛片| 国产亚洲av高清不卡| 无人区码免费观看不卡| 国产主播在线观看一区二区| 黄色毛片三级朝国网站| 69av精品久久久久久| 麻豆久久精品国产亚洲av| 首页视频小说图片口味搜索| 免费在线观看完整版高清| 在线观看舔阴道视频| 18美女黄网站色大片免费观看| 999久久久国产精品视频| 成人亚洲精品av一区二区| 99热只有精品国产| 韩国av一区二区三区四区| 国产精品久久视频播放| 精品电影一区二区在线| 亚洲av电影不卡..在线观看| 国产av不卡久久| 在线十欧美十亚洲十日本专区| 国产色视频综合| 色在线成人网| 人人澡人人妻人| 丝袜美腿诱惑在线| 亚洲成人免费电影在线观看| 欧美日韩一级在线毛片| 999精品在线视频| 日韩有码中文字幕| 波多野结衣高清无吗| 老司机午夜十八禁免费视频| 韩国精品一区二区三区| 午夜激情福利司机影院| 久久欧美精品欧美久久欧美| 成人三级黄色视频| 一级作爱视频免费观看| 美国免费a级毛片| 免费高清视频大片| 国产高清videossex| 草草在线视频免费看| svipshipincom国产片| 亚洲自拍偷在线| 国产色视频综合| 午夜两性在线视频| www.www免费av| 老司机深夜福利视频在线观看| 老汉色av国产亚洲站长工具| 久久国产亚洲av麻豆专区| 嫩草影院精品99| 欧美乱色亚洲激情| 欧美zozozo另类| 侵犯人妻中文字幕一二三四区| 黄片大片在线免费观看| 久久婷婷成人综合色麻豆| 国产成人av教育| 国产野战对白在线观看| 特大巨黑吊av在线直播 | 在线观看免费午夜福利视频| 日韩精品青青久久久久久| 亚洲免费av在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产aⅴ精品一区二区三区波| 欧美日韩亚洲综合一区二区三区_| 亚洲国产高清在线一区二区三 | 男人舔女人的私密视频| 国内久久婷婷六月综合欲色啪| 视频区欧美日本亚洲| 亚洲五月色婷婷综合| 最新在线观看一区二区三区| 亚洲国产欧洲综合997久久, | 欧美色视频一区免费| 国产真人三级小视频在线观看| 成人欧美大片| 1024视频免费在线观看| 免费看a级黄色片| 亚洲在线自拍视频| 亚洲七黄色美女视频| 国产爱豆传媒在线观看 | 亚洲三区欧美一区| 老司机午夜福利在线观看视频| 啦啦啦 在线观看视频| 亚洲色图av天堂| 老司机午夜十八禁免费视频| 岛国视频午夜一区免费看| www日本在线高清视频| 日本一区二区免费在线视频| 久久天堂一区二区三区四区| 搞女人的毛片| 久99久视频精品免费| 丰满的人妻完整版| 又黄又粗又硬又大视频| 99精品久久久久人妻精品| 久久婷婷人人爽人人干人人爱| 久久精品国产综合久久久| 给我免费播放毛片高清在线观看| 欧美在线黄色| 色播在线永久视频| 成人手机av| 日韩精品免费视频一区二区三区| 精品国产国语对白av| 午夜精品在线福利| 亚洲国产看品久久| 9191精品国产免费久久| 一本一本综合久久| 日韩欧美一区视频在线观看| 搡老熟女国产l中国老女人| 好男人在线观看高清免费视频 | 久久国产亚洲av麻豆专区| 18禁裸乳无遮挡免费网站照片 | 在线看三级毛片| 国产精品一区二区三区四区久久 | 国产亚洲精品久久久久5区| 丝袜在线中文字幕| 国产片内射在线| 国产亚洲精品久久久久5区| 999久久久国产精品视频| 久久人妻av系列| 后天国语完整版免费观看| 午夜a级毛片| 国产亚洲精品一区二区www|