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

    基于三種空間預(yù)測模型的海南島土壤有機(jī)質(zhì)空間分布研究*

    2018-08-27 03:29:32姜賽平張懷志張認(rèn)連李兆君謝良商徐愛國
    土壤學(xué)報 2018年4期
    關(guān)鍵詞:變異土層變量

    姜賽平 張懷志 張認(rèn)連 李兆君 謝良商 徐愛國?

    (1 中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所,北京 100081)

    (2 海南省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與土壤研究所,???571100)

    海南島處于熱帶地區(qū),島內(nèi)地形復(fù)雜,土地利用多樣,具有豐富的熱帶作物資源,是我國重要的熱帶作物生產(chǎn)基地。研究該區(qū)土壤有機(jī)質(zhì)(SOM)含量的空間變異規(guī)律,對于了解熱帶氣候條件下、不同熱帶作物種植方式下土壤肥力的空間分布狀況具有重要意義。

    近些年來,隨著“3S”等信息技術(shù)的快速發(fā)展,數(shù)字土壤屬性制圖方法被廣泛應(yīng)用到土壤學(xué)領(lǐng)域中,用以描述土壤屬性的空間變異規(guī)律。常見的數(shù)字土壤屬性制圖方法有地統(tǒng)計學(xué)制圖方法、機(jī)器學(xué)習(xí)方法和混合模型方法等。地統(tǒng)計學(xué)制圖方法是應(yīng)用最廣泛也是最為成熟的方法。其中,普通克里格法(Ordinary kriging,OK)和回歸克里格法(Regression-kriging,RK)最具代表性。OK因方法簡單、容易操作而被普遍接受。但土壤屬性的空間變異同時受地形因子、土地利用類型、土壤類型等眾多因素的影響[1-3]。由于未考慮到輔助變量對土壤屬性的影響,因而成圖只能描述土壤屬性的整體空間分布規(guī)律,對局部信息描述不夠詳細(xì),且平滑效應(yīng)嚴(yán)重[4]。為了更好地揭示土壤屬性空間變異規(guī)律,結(jié)合輔助變量的RK被應(yīng)用于土壤屬性制圖當(dāng)中。趙永存等[5]采用多元線性回歸,泛克里格法和RK分別對河北省土壤有機(jī)碳的空間分布規(guī)律進(jìn)行預(yù)測,結(jié)果表明,結(jié)合地形屬性的RK預(yù)測精度最高,且對土壤有機(jī)碳局部變異信息描述地更加詳細(xì)。連綱等[3]采用該法結(jié)合地形因子與遙感指數(shù)對黃土高原丘陵溝壑區(qū)土壤屬性的空間分布規(guī)律進(jìn)行預(yù)測,結(jié)果表明,該法不但能夠提高制圖精度,且能消除部分平滑效應(yīng)。楊順華等[6]選取相對高程和匯流動力指數(shù)為輔助變量采用該法對湖北省枝城鎮(zhèn)SOM空間分布規(guī)律進(jìn)行預(yù)測,亦取得較高的精度。RK能夠較好地描述地形復(fù)雜區(qū)土壤屬性的空間分布規(guī)律,但該法要求數(shù)據(jù)滿足正態(tài)分布和內(nèi)蘊假設(shè)[7],給實際問題的處理帶來一定的困難。支持向量機(jī)是一種基于統(tǒng)計學(xué)習(xí)理論的機(jī)器學(xué)習(xí)方法,相對于地統(tǒng)計學(xué)方法,它不需要任何前提假設(shè),在解決小樣本、高維數(shù)、非線性以及局部極小點等問題方面有著較高的性能,但其要求輔助變量為連續(xù)型變量,且需對輔助變量進(jìn)行標(biāo)準(zhǔn)化處理,減小了輔助變量空間分布的異質(zhì)性,影響輔助變量與土壤屬性間的協(xié)同變化關(guān)系[8]。目前該法在數(shù)字土壤屬性制圖中應(yīng)用相對較少。神經(jīng)網(wǎng)絡(luò)模型、分類與回歸樹模型是近些年來發(fā)展較為迅速的機(jī)器學(xué)習(xí)方法,對于連續(xù)型變量和類別型變量均適用,且不需對數(shù)據(jù)進(jìn)行預(yù)處理,能夠很好地描述土壤屬性與環(huán)境變量間的非線性關(guān)系,雖有研究表明二者能夠提高制圖精度[9-10],但同時也存在著一些缺陷,例如:神經(jīng)網(wǎng)絡(luò)模型的參數(shù)較多,不易確定[11],模型的預(yù)測結(jié)果也不易解釋[12]。分類與回歸樹模型的結(jié)構(gòu)取決于樣本數(shù)據(jù),對于樣本量較小的數(shù)據(jù),訓(xùn)練樣本的微小變化可能導(dǎo)致分類節(jié)點的不同[13]。并且二者均易過度擬合[14-15]。因此并未廣泛應(yīng)用至土壤屬性制圖當(dāng)中。隨機(jī)森林 (Random Forest,RF)模型在分類與回歸樹模型基礎(chǔ)上發(fā)展而來,模型參數(shù)相對較少,計算簡單,且不易過度擬合[16],被引入數(shù)字土壤屬性制圖當(dāng)中,在地形相對復(fù)雜地區(qū),能夠取得較高的制圖精度[17-18]。

    目前涉及海南全島的SOM含量空間變異特征研究較少,結(jié)合土地利用類型、土壤類型等信息對不同模型的模擬結(jié)果進(jìn)行對比,從而選出適合該區(qū)SOM含量空間變異模型的研究更是鮮見[19-20]。本文根據(jù)地形復(fù)雜區(qū)模型適用情況,結(jié)合土地利用類型、土壤類型和地形因子等輔助變量,選用OK、RK和RF三種方法預(yù)測該區(qū)SOM含量的空間分布特征,通過精度驗證確定出最優(yōu)模型,從而為地處熱帶、較大區(qū)域尺度地形復(fù)雜區(qū)SOM含量空間變異制圖方法的選取提供依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    海南島位于1 8°1 0′~2 0°1 0′N,108°37′~111°3′E之間,是我國第二大島,面積32 900 km2。該島屬于熱帶季風(fēng)氣候,年均氣溫23~25℃,≥10℃年積溫8 200~9 200℃,年平均降水量為1 720mm。島內(nèi)地形復(fù)雜,由山地、丘陵、臺地、平原組成中間高聳、四周低平的環(huán)形層狀地貌。成土母質(zhì)有花崗巖、砂頁巖、淺海沉積物、玄武巖等10種。土壤類型包括磚紅壤、水稻土、赤紅壤、燥紅土、風(fēng)沙土、紫色土等15個土類。土地利用主要為水田、旱地、林地,島內(nèi)以種植熱帶作物為主,種類繁多。

    1.2 樣品采集與分析

    土壤樣品的采集時間為2012年11月—2012年12月,采用近似網(wǎng)格法,綜合考慮當(dāng)?shù)氐牡匦蔚孛?、土壤類型和土地利用類型等因素,以樣點具有代表性、均勻分布性為原則,在全島19個縣市布點,采樣間距平均為10 km,共采集樣點163個,每個樣點按0~5、0~20、20~40、40~60 cm分層采樣,分別研究表層、耕層、中層和中下層SOM含量空間分布狀況。其中,0~5和0~20 cm分別取樣,0~5 cm土鉆取樣點緊鄰0~20 cm土鉆點,記錄樣點經(jīng)緯度、地形和土地利用類型等信息。土壤類型包括磚紅壤、水稻土、赤紅壤、燥紅土、火山灰土、紫色土、風(fēng)沙土、石灰(巖)土等8個土類。土地利用類型包括水田、旱地、園地。

    土壤樣品SOM含量的測定采用重鉻酸鉀-硫酸消化法,具體步驟參見文獻(xiàn)[21]。

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

    輔助變量:RK和RF均在6個地形因子:高程(x1)、坡度(x2)、坡向(x3)、平面曲率(x4)、剖面曲率(x5)、地形濕度指數(shù)(x6),3個環(huán)境因子:歸一化植被指數(shù)(NDVI,x7),土地利用類型(x8)、土壤類型(x9)的基礎(chǔ)上進(jìn)行剔除。

    輔助變量數(shù)據(jù)來源:(1)海南島1︰50 000等高線矢量圖;(2)遙感影像:2013年10月26日和2013年12月6日的海南島Landsat8 OLI_TIRS四景衛(wèi)星影像,數(shù)據(jù)來源于中國科學(xué)院計算機(jī)網(wǎng)絡(luò)信息中心國際科學(xué)數(shù)據(jù)鏡像網(wǎng)站(http://www.gscloud.cn);(3)海南島1︰200 000土壤圖;(4)2010年1︰100 000遙感解譯土地利用圖,數(shù)據(jù)源自中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心。

    1.4 模型方法

    本研究采用OK、RK和RF對SOM含量的空間分布特征進(jìn)行預(yù)測。其中,RF的基本原理參見文獻(xiàn)[22]。樹節(jié)點預(yù)選的變量個數(shù)(mtry)和隨機(jī)森林中樹的個數(shù)(ntree)是RF的兩個重要參數(shù)。本文通過對不同個數(shù)的輔助變量逐次計算確定出最優(yōu)的mtry值[12]。設(shè)定ntree為500、1 000、1 500、2 000,結(jié)合mtry值,選擇使模型袋外誤差(OOB)最小的參數(shù)組合用于最終預(yù)測。

    1.5 數(shù)據(jù)處理與分析

    (1)數(shù)字高程模型的生成。采用1︰50 000等高線矢量圖在ArcGIS9.3中生成10 m分辨率的數(shù)字高程模型柵格影像。

    (2)啞變量賦值。由于土壤類型和土地利用類型為類別變量,不能直接用于回歸分析,因此采用啞變量賦值方法[23]對其賦值,賦值結(jié)果如下:(x91=1、x92=0、x93=0、x94=0、x95=0),(x91=0、x92=1、x93=0、x94=0、x95=0),(x91=0、x92=0、x93=1、x94=0、x95=0),(x91=0、x92=0、x93=0、x94=1、x95=0),(x91=0、x92=0、x93=0、x94=0、x95=1),(x91=0、x92=0、x93=0、x94=0、x95=0)分別代表磚紅壤、水稻土、赤紅壤、燥紅土、火山灰土和紫色土。(x81=1、x82=0),(x81=0、x82=1),(x81=0、x82=0)分別代表水田、旱地、園地。

    (3)異常值處理與驗證集的選取。采用閾值法將區(qū)間(平均值±3×標(biāo)準(zhǔn)差)外的數(shù)值視作異常值并剔除,并參見文獻(xiàn)[2],將風(fēng)沙土、石灰(巖)土樣點剔除,剩余有效樣點160個,從有效樣點中隨機(jī)選取128個樣點作為訓(xùn)練集,用于空間數(shù)據(jù)分析,其余32個樣點用于驗證。全部保留樣點土壤類型覆蓋海南島91.24%的陸域面積。樣點分布見圖1。

    (4)模型精度驗證。采用驗證集的平均預(yù)測誤差(ME)、均方根預(yù)測誤差(RMSE)和決定系數(shù)(R2)來評價預(yù)測的準(zhǔn)確性。

    (5)分析工具。地形因子的提取以及SOM含量空間制圖在ArcGIS9.3軟件中完成,遙感影像鑲嵌和NDVI的提取在ENVI5.1軟件中完成,地統(tǒng)計學(xué)分析和隨機(jī)森林規(guī)則分別在R3.3.1的gstat包和randomForest包中完成,統(tǒng)計分析使用SAS9.2軟件和Microsoft Excel 2010工具。

    圖1 研究區(qū)采樣點分布圖Fig. 1 Distribution of sampling sites in the study area

    2 結(jié)果與討論

    2.1 不同土層SOM含量的描述性統(tǒng)計特征

    由表1可知,隨著土層深度的增加,SOM含量的最小值、最大值逐漸減小,最小值由4.52減小至1.73 g kg-1,最大值由45.83減小至21.53 g kg-1。從表層(0~5 cm)到深層(40~60 cm),SOM含量的均值分別為19.67、15.89、10.30、8.07 g kg-1。按全國第二次土壤普查養(yǎng)分分級標(biāo)準(zhǔn),除40~60 cm土層SOM含量的均值處于五級水平外,其他三個土層均處于四級水平。四個土層SOM含量的標(biāo)準(zhǔn)誤差值均較小,表明各個土層數(shù)據(jù)波動不是很大,結(jié)果較可靠。從變異程度看,四個土層SOM含量的變異系數(shù)在47.16%~56.56%,均屬于中等變異。

    表1 SOM含量基本統(tǒng)計特征Table 1 Statistical characters of SOM contents

    2.2 基于OK的SOM含量空間分布

    對訓(xùn)練集數(shù)據(jù)進(jìn)行Kolmogorov-Smirnov檢驗發(fā)現(xiàn),四個土層原始SOM含量數(shù)據(jù)p均小于0.01,均不符合正態(tài)分布。對其取算數(shù)平方根后,0~5、0~20、20~40和40~60 cm,p>0.05,滿足正態(tài)分布。

    由圖2可知,0~5 cm土層SOM含量共分為6個等級,按占全島面積的比例進(jìn)行排序如下:20~25 (32.63%)>15~20 (25.83%)>10~15 (19.36%)>25~30 (12.99%)>6~10(8.74%)>30~40 g kg-1(0.45%),SOM含量>10 g kg-1的等級占全島面積的91.26%。其中高值(25~30、30~40 g kg-1)主要分布在島內(nèi)東北的瓊山南部、文昌西南部、定安東部、瓊海北部地區(qū)和島中部的瓊中地區(qū),低值(6~10 g kg-1)分布在島西部和西南部的昌江、東方沿海一帶。整體而言,該島SOM含量呈現(xiàn)出西部地區(qū)低于東部地區(qū)的空間分布趨勢,這與前人[19-20]的研究結(jié)果相一致。0~20 cm土層SOM含量的空間變異規(guī)律與表層相似,但SOM含量等級數(shù)減少,>10 g kg-1等級的土壤面積占全島的84.53%,其中10~15 g kg-1和15~20 g kg-1等級所占比例相當(dāng),二者共占全島面積的70.98%。20~40 cm土層>10 g kg-1等級的土壤面積占全島的44.54%。而40~60 cm土層SOM含量>10 g kg-1等級的土壤面積僅占全島的2.85%,該土層主要分布在6~10 g kg-1范圍內(nèi),占全島面積的72.12%。

    2.3 基于RK的SOM含量空間分布

    本文采用逐步線性回歸方程擬合SOM含量與輔助變量之間的趨勢,然后用簡單克里格對剔除趨勢后的殘差進(jìn)行估計,將二者結(jié)果疊加(即RK)作為最終預(yù)測結(jié)果。

    由表2可知,四個土層擬合方程對總方差的解釋率分別為:26.83%、17.12%、18.68%、18.04%,其回歸系數(shù)相對較小,這是因為本研究區(qū)面積較大,SOM與輔助變量之間的關(guān)系較復(fù)雜,而實際進(jìn)入方程的變量較少,故解釋率相對較低[17]。如果考慮氣溫、降水等因素,可能會獲得更高的解釋率。根據(jù)李燕麗等[24]的統(tǒng)計,逐步線性回歸模型對土壤屬性變異的解釋率在15%~82%,本研究結(jié)果在該范圍內(nèi)。且從概率來看,四個土層擬合模型均是極顯著的,表明方程能夠很好地描述SOM含量的變異規(guī)律。NDVI和土壤類型變量在四個土層中均進(jìn)入方程,表明這兩個變量是影響SOM含量的重要因素。有研究表明[1],NDVI和SOM含量呈極顯著正相關(guān)關(guān)系,即該指數(shù)越大,SOM含量越高。海南島西部NDVI值低于東部,與西部較東部干旱,植被生長茂密程度低于東部相一致。海南島土壤類型繁多,研究表明[25],在熱帶土壤中,不同的土壤類型,SOM含量差別較大。曾迪等[26]的研究也表明SOM含量的空間變異受土壤類型的影響。從表2中還可以看出,土地利用類型變量僅在0~5 cm和0~20 cm土層進(jìn)入方程,表明該變量對耕層的影響較大,而對深層影響較小。這是因為該區(qū)土地利用類型以水田、旱地和園地為主,受人為影響較大,隨著土層加深,施肥、耕作等人為擾動減小,所以對SOM含量變異的影響也減小。

    圖2 OK預(yù)測SOM含量空間分布圖Fig. 2 OK-based SOM content spatial distribution map

    表2 各土層逐步線性回歸方程擬合結(jié)果Table 2 Fitting of the soil layers with the stepwise linear regression equation

    由圖3可知,應(yīng)用RK預(yù)測SOM含量的空間變異規(guī)律與OK預(yù)測所得結(jié)果略有不同,呈現(xiàn)出西南、北部、東北高,西部、東南沿海地區(qū)低的空間分布趨勢。RK對細(xì)節(jié)刻畫地更加清晰,這是引入了輔助變量的緣故。該法對SOM含量劃分等級相對于OK增多,各個等級所占的面積有所變化,高低值分布略有不同。0~5 cm土層SOM含量各等級占全島面積的比例排序如下:15~20 (30.28%)>10~15 (22.40%)>20~25 (21.83%)>25~30(10.18%)>6~10 (6.75%)>30~40 (6.46%)>高于40 (1.18%)>低于6 g kg-1(0.84%),與OK相比,該法增加了高于40 g kg-1和低于6 g kg-1兩個等級,對SOM含量的預(yù)測范圍更廣。高值(25~30、30~40 g kg-1)較OK預(yù)測區(qū)域有所增加,主要分布在瓊山、儋州北部、定安南部、瓊中西南部、五指山東北部、保亭西部、三亞東北部、樂東西北部地區(qū),低值(6~10 g kg-1)與OK法所得結(jié)果相似,主要分布在島西部地區(qū)。這主要是引入輔助變量土壤類型所致。以下三個土層與0~5 cm分布規(guī)律相似,但等級數(shù)減少,根據(jù)表1統(tǒng)計,隨著土層深度的增加,SOM含量減少。0~20 cm土層SOM含量>10 g kg-1等級的土壤面積占全島的86.26%,與OK該等級占全島面積的比例相當(dāng),而10~15 g kg-1和15~20 g kg-1較OK有所減少,分別減少5.56%和3.02%,兩個等級共占全島面積的62.40%;20~40 cm土層>10 g kg-1等級的土壤面積占全島比例較OK增加10.55%,與OK相比增加了15~20、20~25、>25 g kg-1三個等級;40~60 cm土層主要分布在6~10 g kg-1(54.35%),與OK法所得結(jié)果相似。

    圖3 RK預(yù)測SOM含量空間分布圖Fig. 3 RK-based SOM content spatial distribution map

    2.4 基于 RF的SOM含量空間分布

    由圖4可知,RF預(yù)測的SOM含量空間分布圖與RK預(yù)測圖相似,呈現(xiàn)出西南、東北高,西部、東南沿海地區(qū)低的空間分布趨勢,因為二者均引入了輔助變量,能夠更好地描述SOM含量空間變異的細(xì)節(jié)信息。對SOM含量等級的劃分,該法在0~5 cm和0~20 cm土層與OK相同,40~60 cm土層與RK相同,20~40 cm土層介于兩者之間。RF四個土層>10 g kg-1等級的土壤面積占全島比例在三種方法中均最高,分別為99.79%、98.08%、60.32%、34.80%。0~5 cm土層SOM含量主要在15~25 g kg-1,占全島面積的74.6%;0~20 cm土層SOM含量主要在10~20 g kg-1(81.72%),與前兩種方法所得結(jié)果一致;20~40 cm土層SOM含量主要在6~15 g kg-1(84.2%),與RK法所得結(jié)果相同;40~60 cm土層SOM含量主要在6~15 g kg-1(89.15%),與前兩種方法相比,10~15 g kg-1等級所占比例明顯增加。

    2.5 不同預(yù)測模型的精度評價

    由表3可知,0~5 cm土層,三種預(yù)測方法中,從R2來看,RF > RK > OK,表明RF對SOM含量空間變異的解釋能力最強(qiáng),其次為RK法,OK法最弱;從RMSE來看,RF < RK < OK,表明RF法預(yù)測精度最高,RK法居中,而OK法預(yù)測精度最低;OK法的ME更接近于0,其次為RK法,而RF法的無偏估計相對較差;從以上三個參數(shù)綜合來看,該土層的最優(yōu)擬合模型為RF。

    圖4 RF預(yù)測SOM含量空間分布圖Fig. 4 RF-based SOM content spatial distribution map

    以此方法分析其他土層,確定0~20、20~40和40~60 cm三個土層的最優(yōu)擬合模型,分別為RF法、RF法和OK法。

    表3 各土層SOM含量預(yù)測模型精度Table 3 Accuracy of the prediction of SOM contents relative to soil layer

    雖然RF對該區(qū)SOM含量空間變異的描述在0~5、0~20和20~40 cm土層表現(xiàn)較好。但由表3可知,驗證集中決定系數(shù)(R2)在0.19~0.37之間,較前人的研究結(jié)果[17-18,27]偏低。這是因為以往的研究采樣密度在0.03~2.52個 km-2,是本研究區(qū)采樣密度(0.004個km-2)的8倍~610倍,采樣密度小可能是造成本研究區(qū)RF預(yù)測精度偏低的原因之一。此外,本研究選取的輔助因子相對較少,有研究表明[19,28],氣溫、降水和人為因素(如施肥等)是影響SOM含量空間變異的重要因素,如果增加上述因子可能會提高模型的解釋率。

    OK在本研究中僅適合40~60 cm土層,且從預(yù)測圖(圖2)來看,該法圖斑較大,僅能夠預(yù)測出SOM含量的整體空間分布規(guī)律,對SOM含量空間變異的細(xì)節(jié)信息描述地不夠詳細(xì)。這是因為本研究區(qū)地形復(fù)雜,而有研究表明,OK通常適用于土壤屬性變化較為均勻的區(qū)域,對地形復(fù)雜、土壤屬性變化較為強(qiáng)烈的區(qū)域,其制圖精度不太理想[29],40~60 cm土層受地形因子的影響相對表層要小,故OK法的制圖精度相對表層和其他方法均有所提高;此外,本研究采樣密度較小,對OK的制圖精度也會有一定的影響。

    RK在本研究區(qū)的應(yīng)用效果不是很理想,決定系數(shù)R2與前人[30]的研究相比較小,這可能是因為剔除趨勢的模型選擇不合適,在本研究中采用的是逐步線性回歸方程來剔除趨勢,而土壤屬性與輔助變量間的關(guān)系常常是復(fù)雜的非線性關(guān)系,且因子之間可能存在交互作用,逐步線性回歸模型難以描述上述信息。以后的研究中可以考慮采用其他回歸模型(如回歸樹或RF等)來剔除土壤屬性與輔助變量間的趨勢。

    從研究結(jié)果中還可以看出,不同土層所選擇的最優(yōu)模型有所不同,這可能是不同土層的影響因子略有不同,不同的模型根據(jù)自身機(jī)理,對于土壤屬性與輔助變量間的關(guān)系的解釋程度不同,因而導(dǎo)致最終的預(yù)測結(jié)果不同。

    3 結(jié) 論

    本研究在較大區(qū)域尺度地形復(fù)雜地區(qū),對海南島128個樣點數(shù)據(jù),采用OK、RK和RF結(jié)合高程、歸一化植被指數(shù)、土地利用類型和土壤類型等輔助變量對該區(qū)0~5、0~20、20~40、40~60 cm四個土層SOM含量的空間分布規(guī)律進(jìn)行預(yù)測,并以32個驗證點進(jìn)行驗證。主要結(jié)論如下:0~5、0~20和20~40 cm土層的最優(yōu)擬合模型均為RF,而40~60 cm土層的最優(yōu)擬合模型為OK,RK和RF相對于OK對SOM含量的局部變異信息描述地更加詳細(xì)。四個土層SOM含量的均值分別為19.67、15.89、10.30、8.07 g kg-1,隨著土層深度的增加,SOM含量逐漸減小。從空間預(yù)測圖來看,四個土層SOM含量均呈現(xiàn)出西南、東北高,西部、東南沿海地區(qū)低的空間分布趨勢。

    猜你喜歡
    變異土層變量
    土釘噴錨在不同土層的支護(hù)應(yīng)用及效果分析
    抓住不變量解題
    也談分離變量
    變異危機(jī)
    變異
    土層 村與人 下
    土層——伊當(dāng)灣志
    土層 沙與土 上
    SL(3,3n)和SU(3,3n)的第一Cartan不變量
    變異的蚊子
    百科知識(2015年18期)2015-09-10 07:22:44
    亚洲av日韩精品久久久久久密| 真实男女啪啪啪动态图| 欧美激情久久久久久爽电影| 久久国产乱子伦精品免费另类| 琪琪午夜伦伦电影理论片6080| 一本久久中文字幕| 一本精品99久久精品77| 久久国产精品影院| 亚洲,欧美精品.| 夜夜爽天天搞| 尤物成人国产欧美一区二区三区| 日本成人三级电影网站| 欧美bdsm另类| 欧美黑人巨大hd| 最新在线观看一区二区三区| 赤兔流量卡办理| 99久久久亚洲精品蜜臀av| 美女黄网站色视频| 真人一进一出gif抽搐免费| 色尼玛亚洲综合影院| 国内精品一区二区在线观看| 久久人人精品亚洲av| 久久久久久久午夜电影| 禁无遮挡网站| 国产亚洲欧美98| 国产aⅴ精品一区二区三区波| 欧美不卡视频在线免费观看| 九九久久精品国产亚洲av麻豆| 亚洲av熟女| 成年免费大片在线观看| 9191精品国产免费久久| 亚洲欧美日韩无卡精品| 狠狠狠狠99中文字幕| 亚洲一区高清亚洲精品| 国产精品国产高清国产av| 深夜a级毛片| 乱人视频在线观看| 女人被狂操c到高潮| 91在线观看av| 自拍偷自拍亚洲精品老妇| 综合色av麻豆| 男人狂女人下面高潮的视频| 欧美不卡视频在线免费观看| 国产精品一及| 国产人妻一区二区三区在| 国产色婷婷99| 成人毛片a级毛片在线播放| 五月玫瑰六月丁香| 亚洲人与动物交配视频| 日日摸夜夜添夜夜添小说| 麻豆久久精品国产亚洲av| 男女床上黄色一级片免费看| 国产成人aa在线观看| 欧美国产日韩亚洲一区| 国产毛片a区久久久久| 婷婷色综合大香蕉| 精品久久久久久久末码| 国产精品一区二区免费欧美| 一进一出抽搐gif免费好疼| 88av欧美| 少妇人妻精品综合一区二区 | av福利片在线观看| 美女高潮喷水抽搐中文字幕| 亚洲欧美日韩高清专用| 看片在线看免费视频| 国产乱人视频| 一个人免费在线观看电影| 免费在线观看成人毛片| 99国产精品一区二区三区| 91av网一区二区| 日日干狠狠操夜夜爽| 特大巨黑吊av在线直播| 久久国产精品影院| 国产高清激情床上av| 亚洲精品影视一区二区三区av| 日韩欧美在线二视频| 老熟妇乱子伦视频在线观看| 88av欧美| bbb黄色大片| 中文字幕高清在线视频| 村上凉子中文字幕在线| av视频在线观看入口| 精品99又大又爽又粗少妇毛片 | 久久国产乱子免费精品| a在线观看视频网站| 最近最新中文字幕大全电影3| 99视频精品全部免费 在线| 日韩欧美三级三区| 欧美色视频一区免费| 国产精品电影一区二区三区| 日本五十路高清| 日本精品一区二区三区蜜桃| 别揉我奶头~嗯~啊~动态视频| 可以在线观看毛片的网站| 久久精品人妻少妇| 嫩草影院精品99| 国产亚洲欧美98| 两性午夜刺激爽爽歪歪视频在线观看| eeuss影院久久| 欧美黑人巨大hd| 精品久久久久久久久久久久久| 九九在线视频观看精品| 国产精品一区二区三区四区久久| 欧美日韩福利视频一区二区| 国产免费男女视频| 国产私拍福利视频在线观看| 在现免费观看毛片| 小蜜桃在线观看免费完整版高清| 免费看a级黄色片| 亚洲av第一区精品v没综合| 色av中文字幕| av在线蜜桃| 琪琪午夜伦伦电影理论片6080| 老司机午夜十八禁免费视频| 色5月婷婷丁香| 亚洲精品成人久久久久久| 欧美在线一区亚洲| 一本精品99久久精品77| 日韩欧美国产一区二区入口| 黄色视频,在线免费观看| 欧美成人免费av一区二区三区| 国产精品国产高清国产av| 久久久成人免费电影| 久久九九热精品免费| 国产激情偷乱视频一区二区| 丁香欧美五月| 亚洲最大成人手机在线| 欧美日本视频| 久久久成人免费电影| 亚洲成av人片免费观看| 91九色精品人成在线观看| 五月伊人婷婷丁香| 午夜老司机福利剧场| 怎么达到女性高潮| 日韩成人在线观看一区二区三区| 高清毛片免费观看视频网站| 欧美高清成人免费视频www| 美女免费视频网站| 中亚洲国语对白在线视频| 真实男女啪啪啪动态图| 午夜影院日韩av| 亚洲黑人精品在线| 老女人水多毛片| 国产高清视频在线播放一区| 69av精品久久久久久| 国产色爽女视频免费观看| 日韩欧美 国产精品| 嫩草影院精品99| 一卡2卡三卡四卡精品乱码亚洲| 久久国产精品影院| 亚洲国产欧美人成| 国产精品久久电影中文字幕| 我的老师免费观看完整版| 在线免费观看的www视频| 禁无遮挡网站| 欧美潮喷喷水| 99国产极品粉嫩在线观看| 亚洲七黄色美女视频| 在线播放国产精品三级| 午夜福利高清视频| 亚洲国产精品成人综合色| 亚洲天堂国产精品一区在线| 精品一区二区三区视频在线观看免费| 天美传媒精品一区二区| 国产大屁股一区二区在线视频| 12—13女人毛片做爰片一| 深爱激情五月婷婷| 色尼玛亚洲综合影院| 每晚都被弄得嗷嗷叫到高潮| 成人午夜高清在线视频| 亚洲在线观看片| av天堂在线播放| 国产精品久久久久久人妻精品电影| 最近最新中文字幕大全电影3| 欧美+亚洲+日韩+国产| 岛国在线免费视频观看| 99精品在免费线老司机午夜| 热99在线观看视频| 老司机午夜十八禁免费视频| 成人特级黄色片久久久久久久| a级毛片免费高清观看在线播放| 极品教师在线视频| 亚洲中文字幕一区二区三区有码在线看| 噜噜噜噜噜久久久久久91| 蜜桃久久精品国产亚洲av| 1000部很黄的大片| 欧美黑人巨大hd| 国产欧美日韩精品一区二区| 久久婷婷人人爽人人干人人爱| 午夜免费激情av| 窝窝影院91人妻| 精华霜和精华液先用哪个| 国产亚洲av嫩草精品影院| 在线国产一区二区在线| 亚洲av美国av| 校园春色视频在线观看| 日韩精品青青久久久久久| 亚洲国产欧洲综合997久久,| 免费看日本二区| 日本黄色视频三级网站网址| 性色av乱码一区二区三区2| 午夜久久久久精精品| 日韩精品青青久久久久久| 天天躁日日操中文字幕| 成人欧美大片| 久久久久亚洲av毛片大全| 国产精品三级大全| 一个人观看的视频www高清免费观看| 99在线视频只有这里精品首页| 国产aⅴ精品一区二区三区波| 成年女人永久免费观看视频| 一区二区三区免费毛片| 成人鲁丝片一二三区免费| 成人精品一区二区免费| 欧美一区二区精品小视频在线| 五月玫瑰六月丁香| 最近在线观看免费完整版| 淫妇啪啪啪对白视频| 国产午夜福利久久久久久| 日本精品一区二区三区蜜桃| 青草久久国产| 久久久久久久亚洲中文字幕 | 亚洲aⅴ乱码一区二区在线播放| 人妻久久中文字幕网| 黄色日韩在线| 有码 亚洲区| 极品教师在线视频| av天堂在线播放| 日韩高清综合在线| 欧美激情久久久久久爽电影| 听说在线观看完整版免费高清| 高清在线国产一区| 3wmmmm亚洲av在线观看| 国产亚洲欧美98| 国产极品精品免费视频能看的| 免费观看人在逋| 免费观看精品视频网站| 宅男免费午夜| 非洲黑人性xxxx精品又粗又长| av专区在线播放| 成人国产一区最新在线观看| 日韩欧美免费精品| 久久精品国产99精品国产亚洲性色| 精品福利观看| 国产欧美日韩一区二区三| a在线观看视频网站| 精品久久久久久久久久免费视频| 757午夜福利合集在线观看| 99在线人妻在线中文字幕| 欧美中文日本在线观看视频| 国产精品久久久久久亚洲av鲁大| 麻豆成人午夜福利视频| 亚洲人成网站在线播放欧美日韩| 偷拍熟女少妇极品色| 国产精品一区二区免费欧美| 丰满的人妻完整版| 人妻久久中文字幕网| 91午夜精品亚洲一区二区三区 | www.熟女人妻精品国产| 精品人妻视频免费看| 亚洲精品久久国产高清桃花| 欧美激情久久久久久爽电影| 制服丝袜大香蕉在线| 免费在线观看日本一区| 亚洲av免费在线观看| 特大巨黑吊av在线直播| 中文字幕av成人在线电影| 午夜老司机福利剧场| 99国产精品一区二区三区| 国产三级黄色录像| 99热这里只有是精品50| 免费av毛片视频| 又爽又黄无遮挡网站| 成人毛片a级毛片在线播放| 欧美zozozo另类| 国产亚洲精品综合一区在线观看| 有码 亚洲区| 国产aⅴ精品一区二区三区波| 精品久久久久久久久亚洲 | 天堂av国产一区二区熟女人妻| 热99在线观看视频| 成人亚洲精品av一区二区| 婷婷精品国产亚洲av在线| 成年人黄色毛片网站| 国语自产精品视频在线第100页| 欧美性感艳星| 99久久99久久久精品蜜桃| 一二三四社区在线视频社区8| 精品久久久久久成人av| 亚洲国产精品999在线| 99精品久久久久人妻精品| 国产精品亚洲一级av第二区| 日本 欧美在线| 欧美激情国产日韩精品一区| 黄色视频,在线免费观看| 脱女人内裤的视频| 欧美午夜高清在线| 欧美国产日韩亚洲一区| 男女做爰动态图高潮gif福利片| 又爽又黄无遮挡网站| 亚洲精品亚洲一区二区| 国产亚洲精品久久久com| 婷婷亚洲欧美| 狠狠狠狠99中文字幕| 91在线观看av| 国产精品自产拍在线观看55亚洲| 国产精品人妻久久久久久| 一二三四社区在线视频社区8| 精品久久久久久久久久久久久| 日韩中字成人| 如何舔出高潮| 欧美在线黄色| 久久99热6这里只有精品| 国产高清视频在线播放一区| 国产淫片久久久久久久久 | 免费搜索国产男女视频| 白带黄色成豆腐渣| 色播亚洲综合网| 好看av亚洲va欧美ⅴa在| 国产精品日韩av在线免费观看| 男女那种视频在线观看| 亚洲av免费在线观看| 直男gayav资源| 婷婷色综合大香蕉| 色av中文字幕| 午夜a级毛片| 老司机午夜福利在线观看视频| 日韩欧美国产一区二区入口| x7x7x7水蜜桃| 国产乱人视频| 一级黄片播放器| ponron亚洲| 亚洲最大成人手机在线| 999久久久精品免费观看国产| 中国美女看黄片| 看十八女毛片水多多多| 一个人观看的视频www高清免费观看| 久久精品国产亚洲av涩爱 | 亚洲色图av天堂| 日本a在线网址| 国产伦一二天堂av在线观看| 一级毛片久久久久久久久女| 熟妇人妻久久中文字幕3abv| 性色av乱码一区二区三区2| 欧美日本视频| 我要搜黄色片| 国产探花极品一区二区| 精品免费久久久久久久清纯| 少妇熟女aⅴ在线视频| 18禁在线播放成人免费| 欧美日韩国产亚洲二区| 国产一区二区亚洲精品在线观看| a级一级毛片免费在线观看| 精品久久久久久久久av| 99热这里只有是精品50| 天堂√8在线中文| 欧美一级a爱片免费观看看| 国语自产精品视频在线第100页| 亚洲av一区综合| 国产午夜精品论理片| 90打野战视频偷拍视频| 美女被艹到高潮喷水动态| 天堂√8在线中文| 看免费av毛片| 精品久久久久久久久av| 真实男女啪啪啪动态图| www.999成人在线观看| 老女人水多毛片| 国产精品,欧美在线| 国产又黄又爽又无遮挡在线| 天堂动漫精品| 国产激情偷乱视频一区二区| 亚洲av电影在线进入| 色综合婷婷激情| 中文字幕熟女人妻在线| 熟妇人妻久久中文字幕3abv| 亚洲aⅴ乱码一区二区在线播放| 久久久久国产精品人妻aⅴ院| 精品免费久久久久久久清纯| 亚洲av五月六月丁香网| 亚洲精品久久国产高清桃花| 999久久久精品免费观看国产| 亚洲最大成人手机在线| 老熟妇仑乱视频hdxx| 在现免费观看毛片| 中文字幕av成人在线电影| 国产成人av教育| 97碰自拍视频| 国产亚洲欧美98| 少妇人妻一区二区三区视频| 老司机福利观看| avwww免费| 亚洲美女搞黄在线观看 | 亚洲欧美日韩无卡精品| 成人亚洲精品av一区二区| av欧美777| 国产久久久一区二区三区| 久久亚洲精品不卡| 久久国产精品人妻蜜桃| 久久精品影院6| 国产av麻豆久久久久久久| 久久久久久久午夜电影| 亚洲无线观看免费| av视频在线观看入口| 国产乱人伦免费视频| 亚洲七黄色美女视频| 婷婷亚洲欧美| 一本久久中文字幕| 老熟妇仑乱视频hdxx| 少妇人妻精品综合一区二区 | 国产欧美日韩一区二区精品| 久久精品人妻少妇| 国产成人欧美在线观看| 亚洲片人在线观看| 老女人水多毛片| 欧美精品国产亚洲| 色综合婷婷激情| 我要搜黄色片| 黄色日韩在线| 1024手机看黄色片| 精品人妻1区二区| 国产精品亚洲av一区麻豆| 又粗又爽又猛毛片免费看| 国产精品美女特级片免费视频播放器| 国产男靠女视频免费网站| 欧美3d第一页| 国产主播在线观看一区二区| 国产色爽女视频免费观看| 国产国拍精品亚洲av在线观看| 国产欧美日韩精品一区二区| 超碰av人人做人人爽久久| 成人高潮视频无遮挡免费网站| 国产精品久久久久久人妻精品电影| 又紧又爽又黄一区二区| 一级作爱视频免费观看| 99国产极品粉嫩在线观看| 老司机午夜福利在线观看视频| 日本 av在线| 中文字幕人妻熟人妻熟丝袜美| 他把我摸到了高潮在线观看| 日韩大尺度精品在线看网址| 床上黄色一级片| 小蜜桃在线观看免费完整版高清| 日本a在线网址| 欧美最新免费一区二区三区 | 亚洲经典国产精华液单 | 自拍偷自拍亚洲精品老妇| 欧美成人免费av一区二区三区| 国产免费av片在线观看野外av| 99热这里只有是精品50| 综合色av麻豆| 国产精品一区二区三区四区免费观看 | 国产极品精品免费视频能看的| 久99久视频精品免费| 国产av在哪里看| 校园春色视频在线观看| 欧美黄色片欧美黄色片| 中文字幕高清在线视频| 女人十人毛片免费观看3o分钟| 女同久久另类99精品国产91| 特级一级黄色大片| 啪啪无遮挡十八禁网站| 亚洲中文字幕一区二区三区有码在线看| 国产美女午夜福利| 婷婷色综合大香蕉| 日本一本二区三区精品| 亚洲狠狠婷婷综合久久图片| 99久久九九国产精品国产免费| 亚洲五月婷婷丁香| 人人妻,人人澡人人爽秒播| 亚洲av不卡在线观看| 免费电影在线观看免费观看| 99久久精品一区二区三区| 我要搜黄色片| 亚洲人成网站在线播放欧美日韩| 亚洲人与动物交配视频| 国产不卡一卡二| 国产 一区 欧美 日韩| 久久久久久久午夜电影| 黄色一级大片看看| 给我免费播放毛片高清在线观看| 欧美+亚洲+日韩+国产| 亚洲欧美精品综合久久99| 宅男免费午夜| netflix在线观看网站| 两个人视频免费观看高清| 99视频精品全部免费 在线| 亚洲成人久久爱视频| 日韩欧美精品v在线| 成熟少妇高潮喷水视频| 亚洲成人中文字幕在线播放| 国产精品亚洲av一区麻豆| 97碰自拍视频| 欧美色视频一区免费| h日本视频在线播放| 亚洲aⅴ乱码一区二区在线播放| 悠悠久久av| 亚洲无线在线观看| 午夜精品一区二区三区免费看| 久久精品国产清高在天天线| 日本三级黄在线观看| 3wmmmm亚洲av在线观看| 久久草成人影院| 亚洲欧美日韩无卡精品| 性欧美人与动物交配| 最好的美女福利视频网| 免费在线观看亚洲国产| av天堂在线播放| 99久久精品国产亚洲精品| 久久欧美精品欧美久久欧美| 无人区码免费观看不卡| 此物有八面人人有两片| 欧美xxxx性猛交bbbb| 大型黄色视频在线免费观看| 欧美日韩福利视频一区二区| 一进一出抽搐gif免费好疼| 亚洲国产色片| 中文字幕av成人在线电影| 精品一区二区三区视频在线| 日本五十路高清| 色尼玛亚洲综合影院| 琪琪午夜伦伦电影理论片6080| 国产毛片a区久久久久| 好男人在线观看高清免费视频| 精品人妻视频免费看| 国产高清三级在线| 欧美性感艳星| 欧美三级亚洲精品| 桃色一区二区三区在线观看| 色吧在线观看| 欧美+日韩+精品| 国产亚洲精品久久久com| 村上凉子中文字幕在线| 亚洲人成网站高清观看| 在线观看美女被高潮喷水网站 | 国产伦人伦偷精品视频| 国产一级毛片七仙女欲春2| 欧美3d第一页| 91麻豆精品激情在线观看国产| 免费av观看视频| 一区二区三区免费毛片| av中文乱码字幕在线| 亚洲天堂国产精品一区在线| 偷拍熟女少妇极品色| 亚洲电影在线观看av| 首页视频小说图片口味搜索| 深夜精品福利| 伊人久久精品亚洲午夜| av女优亚洲男人天堂| 亚洲精品成人久久久久久| 免费观看精品视频网站| 亚洲国产精品成人综合色| 精品久久国产蜜桃| 婷婷亚洲欧美| 午夜激情欧美在线| 色综合亚洲欧美另类图片| 少妇的逼水好多| 亚洲精品影视一区二区三区av| 此物有八面人人有两片| 一本精品99久久精品77| 一个人免费在线观看电影| 我的老师免费观看完整版| 校园春色视频在线观看| 午夜福利免费观看在线| 熟女电影av网| 蜜桃亚洲精品一区二区三区| 国产黄a三级三级三级人| 国产乱人伦免费视频| 精品无人区乱码1区二区| 国产私拍福利视频在线观看| 国产成+人综合+亚洲专区| 欧美乱妇无乱码| 国产一区二区在线观看日韩| 麻豆一二三区av精品| 国模一区二区三区四区视频| 久久人人精品亚洲av| 人妻丰满熟妇av一区二区三区| 69av精品久久久久久| 97碰自拍视频| av在线老鸭窝| 欧美激情在线99| 亚洲av成人不卡在线观看播放网| 国产综合懂色| 91在线精品国自产拍蜜月| 精品人妻偷拍中文字幕| 亚洲天堂国产精品一区在线| 性色av乱码一区二区三区2| 国产高清视频在线播放一区| 直男gayav资源| 欧美一区二区亚洲| 欧美日韩黄片免| 欧美另类亚洲清纯唯美| 亚洲av电影在线进入| 欧美三级亚洲精品| 欧美另类亚洲清纯唯美| 嫩草影院精品99| 亚洲一区二区三区色噜噜| 蜜桃久久精品国产亚洲av| 听说在线观看完整版免费高清| 国产一区二区三区视频了| 亚洲美女搞黄在线观看 | 亚洲内射少妇av| 亚洲精品一卡2卡三卡4卡5卡| 在线免费观看不下载黄p国产 | 亚洲av免费高清在线观看| av欧美777| 91狼人影院| 麻豆成人午夜福利视频| 他把我摸到了高潮在线观看| 国产探花极品一区二区| 无人区码免费观看不卡| 亚州av有码| 精品福利观看| 国产精品亚洲av一区麻豆| 国产一区二区亚洲精品在线观看|