王曉軍, 武江濤, 王 兵, 溫奮翔

(山西大學(xué) 環(huán)境與資源學(xué)院， 山西 太原 030006)

基于影像分析的黃土丘陵溝壑區(qū)土壤水分采樣研究

王曉軍, 武江濤, 王 兵, 溫奮翔

(山西大學(xué) 環(huán)境與資源學(xué)院， 山西 太原 030006)

摘要：[目的] 以山西省壽陽縣三眼井小流域土壤水分采樣為例，對(duì)采樣點(diǎn)進(jìn)行了均勻性和代表性分析，以此驗(yàn)證土壤采樣設(shè)計(jì)的合理性。[方法] 利用統(tǒng)計(jì)分析SPSS 13.0軟件和ArcGIS軟件，以規(guī)則網(wǎng)格采樣點(diǎn)作為初步布設(shè)圖，然后結(jié)合遙感影像，把自動(dòng)生成的規(guī)則網(wǎng)格采樣點(diǎn)在兼顧土地利用類型和可操作性的前提下進(jìn)行微調(diào)，即采用大均勻小隨機(jī)布點(diǎn)法進(jìn)行土壤采樣點(diǎn)布設(shè)。[結(jié)果] 最近鄰點(diǎn)指數(shù)和變異系數(shù)均顯示采樣點(diǎn)分布較均勻；從采樣點(diǎn)的坡度代表性來看，由于大均勻小隨機(jī)采樣考慮到實(shí)際可操作性，采樣點(diǎn)布設(shè)多集中在0°～15°，坡度的代表性較差；從采樣點(diǎn)的坡向代表性來看，采樣點(diǎn)代表性較好；從土地利用類型代表性來看，大均勻小隨機(jī)法得到的采樣點(diǎn)對(duì)草地與耕地的代表性較好，對(duì)林地的代表性較差，主要因?yàn)榱值囟喾植荚诙钙律想y以采樣。[結(jié)論] 總體來看，大均勻小隨機(jī)采樣法得到的樣點(diǎn)均勻性和代表性都比較好，且由于黃土丘陵溝壑區(qū)溝壑縱橫，地形復(fù)雜，利用此方法可以兼顧均勻性、代表性和可操作性。

關(guān)鍵詞：GIS; 大均勻小隨機(jī)布點(diǎn); 規(guī)則網(wǎng)格布點(diǎn); 均勻性; 代表性

土壤水分作為植物生存的一個(gè)基本因子，在水土流失嚴(yán)重的黃土高原丘陵溝壑區(qū)卻是制約黃土高原植被建設(shè)的關(guān)鍵限制因子，因此土壤水分的研究一直是黃土高原水土資源合理利用和生態(tài)恢復(fù)的主要內(nèi)容[1-2]。要了解研究區(qū)真實(shí)的土壤水分狀況，依賴于所獲得的土壤水分?jǐn)?shù)據(jù)的可靠性和有效性，因此要求采樣點(diǎn)盡可能多，但是過多的采樣不僅費(fèi)時(shí)費(fèi)力，其時(shí)效性也受到影響，而釆樣點(diǎn)過少，則可能丟失重要的土壤性質(zhì)空間信息。因此，建立科學(xué)的土壤采樣布點(diǎn)模式，對(duì)于降低土壤采樣成本、提高預(yù)測(cè)制圖質(zhì)量、客觀真實(shí)地反映研究區(qū)現(xiàn)狀等具有重要的意義[3-6]，特別是針對(duì)黃土高原地區(qū)小流域這種復(fù)雜的地形條件，如何合理制定采樣方案顯得尤為重要。

目前國內(nèi)外常用的土壤采樣布點(diǎn)方法主要有主觀判斷采樣(judgmental sampling)、簡(jiǎn)單隨機(jī)采樣(simple random sampling)、規(guī)則網(wǎng)格采樣(grid sampling)、分區(qū)采樣(stratified sampling)以及混合采樣(composite sampling)等[7-8]。而在土壤采樣模式中，Wang等[9]認(rèn)為主要有簡(jiǎn)單隨機(jī)采樣模式、規(guī)則網(wǎng)格采樣模式。簡(jiǎn)單隨機(jī)采樣雖然保證每個(gè)被抽取的樣本概率相同，且相互獨(dú)立，保證了抽取的無偏差性，但缺點(diǎn)之一是常出現(xiàn)樣本不均勻成堆分布的情況，以及在相同花費(fèi)下采樣方差比其他采樣模式大。而規(guī)則網(wǎng)格采樣則可確保樣點(diǎn)能覆蓋到整個(gè)區(qū)域，能更好的獲取整體信息。但是網(wǎng)格法比較適用于某個(gè)坡面或某種單一的土地利用方式下環(huán)境因素相對(duì)均一、簡(jiǎn)單的土壤空間異質(zhì)性的研究，可以較好地反映研究區(qū)土壤性質(zhì)的空間相關(guān)性[10-13]。因此，如果土壤采樣點(diǎn)同時(shí)結(jié)合這兩種布點(diǎn)方法，在規(guī)則網(wǎng)格布點(diǎn)的基礎(chǔ)上兼顧采樣實(shí)際可操作性對(duì)采樣點(diǎn)進(jìn)行調(diào)整，即大均勻小隨機(jī)方法，則可更加合理地布置采樣點(diǎn)。目前，考慮到黃土高原丘陵溝壑區(qū)地形破碎，景觀高度異質(zhì)性，土壤采樣多根據(jù)不同的土地利用方式、地形狀況以及典型坡面法來布設(shè)樣點(diǎn)[14-15]，相比較而言，利用規(guī)則網(wǎng)格和隨機(jī)布點(diǎn)相結(jié)合的土壤采樣在黃土丘陵溝壑區(qū)的研究中并不多見。

本研究在遙感影像的基礎(chǔ)上，利用統(tǒng)計(jì)分析方法和GIS空間分析技術(shù)，在土壤樣品采集前初步對(duì)大均勻小隨機(jī)樣點(diǎn)分布的均勻性和代表性進(jìn)行分析，以驗(yàn)證此次土壤采樣設(shè)計(jì)的合理性，為黃土高原丘陵溝壑區(qū)更合理地進(jìn)行土壤采樣提供方法和依據(jù)。

1研究區(qū)域概況

本文研究區(qū)域位于山西省壽陽縣的三眼井流域(113°1′43″—113°2′43″E，37°45′45″—37°45′53″N)，流域長(zhǎng)約2 km，寬約1 km，面積約2 km2；主溝道長(zhǎng)約2 km，寬約30 m，主溝方向?yàn)楸蹦舷?，北低南高，次?jí)溝谷多呈東西向分布。該區(qū)主要為山區(qū)與平原區(qū)的轉(zhuǎn)折部分，地表黃土覆蓋，第四系的馬蘭黃土和離石黃土廣泛發(fā)育，是典型的黃土丘陵區(qū)，主要由梁和深度切溝組成，海拔在996～1 160 m。該地區(qū)年均氣溫7.3 ℃，年均降雨量為581.3 mm，降水主要集中在7—9月，全年蒸發(fā)量介于1 800～2 000 mm。該區(qū)的土地利用類型以旱作梯田為主，玉米為主要的種植作物，耕作方式以牲畜耕作為主。植被主要有松樹、楊樹、刺槐以及蒿類、禾本科草類等灌木和草本植物。

2材料與方法

2.1　采樣點(diǎn)的布設(shè)

采用高清晰正射影像和實(shí)測(cè)1∶2 000地形圖，利用漁網(wǎng)工具生成150 m×150 m規(guī)則網(wǎng)格，以自動(dòng)生成的規(guī)則網(wǎng)格采樣點(diǎn)為初步布設(shè)圖，結(jié)合遙感影像，兼顧實(shí)際可操作性，對(duì)規(guī)則網(wǎng)格中心樣點(diǎn)進(jìn)行適當(dāng)?shù)恼{(diào)整，盡量使調(diào)整后的樣點(diǎn)在其原來對(duì)應(yīng)的網(wǎng)格內(nèi)，相對(duì)均勻的布設(shè)101個(gè)采樣點(diǎn)，其中在水井周圍布設(shè)3個(gè)采樣點(diǎn)，打印并塑封采樣點(diǎn)圖，便于在野外實(shí)際采樣過程中及時(shí)調(diào)整。

本文的土壤剖面水分樣品采集于2013年春季，采用傳統(tǒng)土鉆法鉆取0—300 cm剖面土壤，測(cè)定土壤水分。在采樣過程中，首先觀察草圖的采樣點(diǎn)是否適宜采樣，然后在采樣的同時(shí)記錄《采樣點(diǎn)調(diào)查表》和《壽陽三眼井流域土壤水分測(cè)定表》。如果草圖中的采樣點(diǎn)難以實(shí)行，例如位于極陡的坡上等情況，則在附近重新選取采樣點(diǎn)。通過室內(nèi)布設(shè)的采樣點(diǎn)草圖和實(shí)際采樣點(diǎn)圖比較，有8個(gè)采樣點(diǎn)進(jìn)行了重新調(diào)整，采樣點(diǎn)共計(jì)101個(gè)。實(shí)際采樣點(diǎn)分布情況如圖1所示。

圖1　土壤剖面水分樣品實(shí)際采樣點(diǎn)分布

2.2　數(shù)據(jù)處理與分析方法

數(shù)據(jù)的處理與分析采用Excel制作表格，SPSS 13.0進(jìn)行經(jīng)典統(tǒng)計(jì)分析，運(yùn)用GIS的空間分析技術(shù)提取采樣點(diǎn)的土地利用類型與地形信息。

2.2.1采樣點(diǎn)均勻性分析方法一方面將采用大均勻小隨機(jī)法、規(guī)則網(wǎng)格法和集中布點(diǎn)法3種方法得到采樣點(diǎn)分布，以規(guī)則網(wǎng)格為單元統(tǒng)計(jì)每個(gè)網(wǎng)格中的采樣點(diǎn)個(gè)數(shù)，計(jì)算出平均值、標(biāo)準(zhǔn)差和變異系數(shù)，通過變異系數(shù)來分析此次采樣點(diǎn)空間分布的均勻性；另一方面還用平均最鄰近距離工具[16]對(duì)采樣點(diǎn)分布均勻性進(jìn)行分析。

2.2.2采樣點(diǎn)代表性分析方法利用ArcGIS軟件結(jié)合實(shí)地調(diào)查繪制土地利用圖，根據(jù)實(shí)測(cè)1∶2 000地形圖生成DEM數(shù)字高程模型，并利用空間分析方法繪制坡度和坡向圖，將大均勻小隨機(jī)法和規(guī)則網(wǎng)格法得到的采樣點(diǎn)分布圖與土地利用圖、坡度圖、坡向圖疊加提取相關(guān)信息[17]，對(duì)采樣點(diǎn)的分布狀況與流域中各種土地利用類型、坡度類型、坡向類型分布狀況進(jìn)行比較分析其代表性。

本文采用黃河水利委員會(huì)綏德水保試驗(yàn)站在黃土丘陵溝壑區(qū)韭園溝多年試驗(yàn)結(jié)果制定的臨界坡度分級(jí)系統(tǒng)，按照3°，8°，15°，25°，35°及45°的臨界坡度值將流域內(nèi)的坡度分成7個(gè)等級(jí)[18-19]。坡向分為平坦、東、西、南、北、東北、東南、西南和西北9個(gè)等級(jí)。

3結(jié)果與分析

3.1　采樣點(diǎn)分布均勻性分析

李梅等[7]在GIS支持下豫東地區(qū)土壤野外采樣布點(diǎn)方法探索中采用規(guī)則網(wǎng)格法，同時(shí)采用表層采樣與剖面采樣相結(jié)合、網(wǎng)格點(diǎn)與類型控制點(diǎn)相結(jié)合的方法，最后得出采樣密度為每4 km2一個(gè)點(diǎn)，通過公式n=Z2p(1-p)/d2驗(yàn)證采樣點(diǎn)間距的合理性，結(jié)果顯示采樣點(diǎn)間距可行，均勻性比較好。同時(shí)通過與通許縣土壤圖做疊置分析表明，采樣點(diǎn)基本涵蓋了所有的土種類型，代表性好。李玲等[5]以鄭州市土壤肥力監(jiān)測(cè)為例，通過統(tǒng)計(jì)不同區(qū)域和不同土壤類型的采樣點(diǎn)數(shù)量，分析了采樣點(diǎn)的均勻性和代表性，其采樣方法也是以規(guī)則網(wǎng)格法為基礎(chǔ)，在布點(diǎn)中參考土壤圖和土地利用規(guī)劃圖及第二次土壤普查時(shí)的主剖面及農(nóng)化樣點(diǎn)位置(盡可能的落在同一田塊中)，同時(shí)兼顧土壤類型面積大小上的差異，結(jié)果表明，采樣點(diǎn)的均勻性和代表性較好。

表1是不同采樣布點(diǎn)方式下得出的相應(yīng)參數(shù)，結(jié)合附圖1—2可以看出，實(shí)際采樣點(diǎn)數(shù)在每個(gè)網(wǎng)格單元中多為1個(gè)。實(shí)際采樣點(diǎn)數(shù)量的標(biāo)準(zhǔn)偏差、變異系數(shù)分別為0.69和0.95，規(guī)則網(wǎng)格法得到的采樣點(diǎn)數(shù)量的標(biāo)準(zhǔn)偏差、變異系數(shù)分別為0，集中布點(diǎn)方式的標(biāo)準(zhǔn)偏差、變異系數(shù)分別為8.54和11.75。從3種方法布設(shè)的采樣點(diǎn)數(shù)量分布的變異系數(shù)來看，基于遙感影像的大均勻小隨機(jī)采樣方法獲得的采樣點(diǎn)分布比較均勻，十分接近于規(guī)則網(wǎng)格法。

第二種方法采用平均最鄰近距離工具分析測(cè)量樣點(diǎn)均勻性。該方法采用最近鄰點(diǎn)指數(shù)即測(cè)量距離與預(yù)期距離的比例(預(yù)期距離是由同樣數(shù)目的同類要素在同樣區(qū)域假定為隨機(jī)分布得來)分析采樣點(diǎn)的均勻性，最近鄰點(diǎn)指數(shù)見公式(1)：

(1)

表1　不同采樣布點(diǎn)方式下的統(tǒng)計(jì)信息

本文研究區(qū)域的總點(diǎn)數(shù)為101，面積為1.39 km2，根據(jù)最近鄰點(diǎn)指數(shù)分析結(jié)果顯示最近鄰距離值為106 m，帶入公式(1)，得到最近點(diǎn)指數(shù)值為1.9，說明研究區(qū)采樣點(diǎn)為均勻分布。

綜上分析可以看出，利用變異系數(shù)和最近鄰近指數(shù)法得到的結(jié)果是一致的，基于遙感影像的大均勻小隨機(jī)方法得到的采樣點(diǎn)分布均勻。

3.2　采樣點(diǎn)代表性分析

由于土地利用方式、坡度、坡向、高程等對(duì)土壤各種性狀的空間分布有重要的影響，因此本文通過坡度、坡向、土地利用方式對(duì)以遙感影像為基礎(chǔ)的大均勻小隨機(jī)采樣方法獲得的采樣點(diǎn)進(jìn)行了研究。

3.2.1采樣點(diǎn)在不同地形條件的分布情況

(1) 采樣點(diǎn)在不同坡度范圍的分布情況。根據(jù)測(cè)繪得到1∶2 000的地形圖生成DEM數(shù)字高程模型，在此基礎(chǔ)上按照3°，8°，15°，25°，35°與45°的臨界坡度分級(jí)系統(tǒng)重分類后得到了研究區(qū)坡度圖，統(tǒng)計(jì)不同坡度范圍的面積；將坡度圖與規(guī)則網(wǎng)格布點(diǎn)、實(shí)際采樣布點(diǎn)進(jìn)行疊加分析(附圖1—2)，提取采樣點(diǎn)的坡度信息，得出采樣點(diǎn)在各坡度范圍分布的統(tǒng)計(jì)特征(表2)，如果采樣點(diǎn)的個(gè)數(shù)比例與相應(yīng)的坡度范圍面積比例越接近，則說明采樣點(diǎn)代表性越好。從表2可以看出，研究區(qū)按坡度從低到高的順序，不同范圍面積比例依次為10%，22%，17%，14%，12%，11%和14%，坡度主要集中在3°～15°之間，無坡地最少。

規(guī)則網(wǎng)格法采樣點(diǎn)在不同坡度范圍內(nèi)分布的比例依次為12%，26%，18%，11%，12%，9%和12%，由此可以看出，規(guī)則網(wǎng)格布點(diǎn)法與研究區(qū)不同坡度范圍的面積比例基本相符，代表性很好，但是由于規(guī)則中心網(wǎng)格布點(diǎn)落在45°～90°的點(diǎn)有16個(gè)，如果某些地方坡度過陡，不便于采樣，為了兼顧實(shí)際采樣方案的可操作性，對(duì)規(guī)則網(wǎng)格布點(diǎn)進(jìn)行了適當(dāng)?shù)恼{(diào)整后得出大均勻小隨機(jī)采樣方案，通過分析得出采樣點(diǎn)在不同坡度范圍內(nèi)分布的比例依次為21%，43%，21%，6%，4%，2%以及3%，與研究區(qū)不同坡度范圍的面積比例相比，采樣點(diǎn)在<15°內(nèi)分布的比例大于研究區(qū)面積比例，在>15°內(nèi)的分布情況則相反，所以對(duì)坡度而言雖然此次采樣方案的可操作性好，但是其代表性較差。

表2　基于DEM導(dǎo)出的坡度信息統(tǒng)計(jì)量

(2) 采樣點(diǎn)在不同坡向范圍的分布情況。根據(jù)測(cè)繪得到1∶2 000的地形圖生成DEM數(shù)字高程模型，在此基礎(chǔ)上，按坡向分類標(biāo)準(zhǔn)得到了研究區(qū)坡向分布圖，將坡向圖與規(guī)則網(wǎng)格采樣點(diǎn)分布圖、實(shí)際采樣點(diǎn)分布圖進(jìn)行疊加(附圖3—4)，得出研究區(qū)域各類坡向的統(tǒng)計(jì)特征值(表3)。如果采樣點(diǎn)的個(gè)數(shù)比例與相應(yīng)的坡向范圍面積比例越接近，則說明采樣點(diǎn)代表性越好。從表3可以看出，研究區(qū)坡向等級(jí)1—10所占面積比例依次為4%，8%，16%，18%，13%，8%，6%，9%，11%和7%，大均勻小隨機(jī)采樣點(diǎn)在不同坡向范圍內(nèi)的比例依次為8%，4%，27%，21%，5%，4%，1%，11%，12%和7%，規(guī)則網(wǎng)格采樣點(diǎn)在不同坡向范圍內(nèi)的比例依次為2%，6%，17%，17%，17%，6%，1%，16%，12%和6%，可以看出，規(guī)則中心網(wǎng)格布點(diǎn)與該坡向范圍面積所占比例基本相似，但結(jié)合坡度圖(附圖1—2)可以看出，坡向在112°～157°之間，規(guī)則網(wǎng)格法采樣點(diǎn)的坡度都在45°以上，在實(shí)際采樣中難以實(shí)現(xiàn)。因此，基于所占比例，綜合實(shí)際情況和數(shù)據(jù)的有效性分析，說明大均勻小隨機(jī)采樣點(diǎn)代表性比較好。

3.2.2采樣點(diǎn)在不同土地利用方式的分布情況將土地利用圖與規(guī)則網(wǎng)格采樣點(diǎn)圖、實(shí)際采樣點(diǎn)圖進(jìn)行疊加，得到不同土地利用類型分布圖(附圖5—6)，面積比例以及采樣點(diǎn)在不同土地利用類型中的分布比例，分析結(jié)果見表4。如果采樣點(diǎn)的個(gè)數(shù)比例與相應(yīng)的土地利用類型面積比例越接近，則說明采樣點(diǎn)代表性越好。從表4可以看出，研究區(qū)域的土地利用類型主要有耕地、林地、草地、果園以及居民宅基地，所占面積比例依次為51%，36%，9%，1%和3%，研究區(qū)域的土地利用類型以耕地和林地(主要是灌木林地)為主。此外，從表2可以看出，規(guī)則網(wǎng)格采樣點(diǎn)分布在耕地、林地、草地、果園和居民宅基地的比例依次為53%，33%，11%，1%和2%，大均勻小隨機(jī)采樣點(diǎn)個(gè)數(shù)落在耕地、林地、草地的比例分布依次為77%，15%和8%，與土地利用類型面積比例相比，規(guī)則網(wǎng)格的個(gè)數(shù)比例更接近土地利用類型面積比例，但是，ArcGIS自動(dòng)生成的規(guī)則網(wǎng)格采樣點(diǎn)不一定都是符合要求的有效點(diǎn)。例如有3個(gè)點(diǎn)落到建筑區(qū)域，由于城鎮(zhèn)、住宅、道路、溝渠、墳?zāi)垢浇忍幦藶楦蓴_大，失去土壤的代表性，同時(shí)由于黃土高原地形破碎，還有部分采樣點(diǎn)落在溝谷陡坡上，在實(shí)際采樣過程中難于操作，有一定的局限性。而與研究區(qū)中不同土地利用類型面積比例相比，本文采用的大均勻小隨機(jī)采樣法得到的采樣點(diǎn)在草地中的分布比例更為接近，耕地中的較多而林地較少。通過對(duì)土地利用圖與坡度圖做疊加分析，得到基于坡度分級(jí)的土地利用圖(附圖7)，從表5和附圖7可以看出，林地主要分布在大于25°的陡坡上，比例為67%，其中在坡度為45°~90°之間分布最多，比例達(dá)到27%。耕地主要分布在0°～15°之間，比例為76%。以上結(jié)果看出林地在大于45°的坡上分布較多，采樣比較困難，因此，導(dǎo)致實(shí)際采樣點(diǎn)在林地中分布的比例較小。

表3　基于DEM導(dǎo)出的坡向信息統(tǒng)計(jì)量

表4　各類土地利用類型的信息統(tǒng)計(jì)量

表5　基于坡度分級(jí)的土地利用類型的面積及所占面積比例

4結(jié) 論

(1) 本研究區(qū)域采樣點(diǎn)數(shù)量分布的均勻性分析表明，以遙感影像為基礎(chǔ)，采用大均勻小隨機(jī)的布點(diǎn)方式，采樣點(diǎn)分布的變異系數(shù)接近于規(guī)則網(wǎng)格布點(diǎn)，表明此次布點(diǎn)比較均勻，最近鄰點(diǎn)指數(shù)也說明此次大均勻小隨機(jī)布點(diǎn)比較均勻。

(2) 從大均勻小隨機(jī)采樣點(diǎn)的代表性分析可知，整體來看，通過對(duì)規(guī)則網(wǎng)格中心法獲得的采樣點(diǎn)在不同地形和土地利用類型中分布的比例更接近相應(yīng)的坡度范圍、坡向范圍、土地利用的面積比例的原因分析出大均勻小隨機(jī)布點(diǎn)兼顧實(shí)際可操作性，并且大均勻小隨機(jī)采樣主要分布在0°～25°緩坡和緩陡坡的耕地和林地上，坡向范圍主要集中在22°～67°和67°～112°，所占比例基本與其面積比例一致，說明此次大均勻小隨機(jī)布點(diǎn)代表性比較好。

綜上所述，結(jié)果表明利用隨機(jī)布點(diǎn)與規(guī)則網(wǎng)格布點(diǎn)相結(jié)合的大均勻小隨機(jī)的方法進(jìn)行土壤采樣布點(diǎn)優(yōu)化在黃土高原丘陵溝壑區(qū)域具有較好的適用性，如果在此基礎(chǔ)上，考慮到土地利用、坡度、坡向的分布情況，可以更精確得進(jìn)行土壤采樣布點(diǎn)。同時(shí)還表明小流域或區(qū)域中土壤采樣均勻性與代表性需要運(yùn)用ArcGIS的空間分析方法，把數(shù)字地形分析的結(jié)果以及土地利用現(xiàn)狀圖運(yùn)用到大均勻小隨機(jī)采樣布點(diǎn)上，使土壤剖面研究從采樣點(diǎn)布設(shè)到分析的每一個(gè)環(huán)節(jié)都變得更為便捷，并且在樣點(diǎn)布置的合理性、取樣的科學(xué)性和數(shù)據(jù)分析結(jié)果的準(zhǔn)確可靠性等方面都得到了顯著的提升。因此，利用大均勻小隨機(jī)的方法，結(jié)合土地利用、坡度、坡向的特點(diǎn)進(jìn)行采樣布點(diǎn)，可以滿足均勻性和代表性的要求，又能減少樣品采集過程中的誤差，避免造成資金和人力上的浪費(fèi)。

[參考文獻(xiàn)]

[1]王栓全,岳宏昌,王偉.黃土丘陵溝壑區(qū)不同土地類型的土壤水分特性[J].干旱地區(qū)農(nóng)業(yè)研究,2009,27(6):93-96.

[2]鄺高明,朱清科,劉中奇,等.黃土丘陵溝壑區(qū)微地形對(duì)土壤水分及生物量的影響[J].水土保持研究,2012,19(3):74-77.

[3]陳天恩,董靜,陳立平,等.縣域農(nóng)田土壤采樣布局多目標(biāo)優(yōu)化分析模型[J].農(nóng)業(yè)工程學(xué)報(bào),2012,28(23):67-73.

[4]蘇曉燕,趙永存,楊浩,等.不同釆樣點(diǎn)數(shù)量下土壤有機(jī)質(zhì)含量空間預(yù)測(cè)方法對(duì)比[J].地學(xué)前緣,2012,18(6):34-40.

[5]李玲,陳偉強(qiáng),江輝,等.3S在土壤布點(diǎn)與采樣中的應(yīng)用[J].中國農(nóng)學(xué)通報(bào),2007,23(6):388-391.

[6]李翔,潘瑜春,趙春江,等.利用不同方法估測(cè)土壤有機(jī)質(zhì)及其對(duì)采樣數(shù)的敏感性分析[J].地理科學(xué),2007,27(5):689-694.

[7]李梅,張學(xué)雷,武繼承.GIS支持下豫東地區(qū)土壤野外采樣布點(diǎn)方法探索[J].土壤,2011,43(3):459-465.

[8]Peigné J, Vian J F, Cannavacciuolo M, et al. Soil sampling based on field spatial variability of soil microbial indicators[J]. European Journal of Soil Biology, 2009(45):488-495.

[9]Englund E, Weber D, Leviant N. The effects of sampling design parameters on block selection[J].Mathematical Geology, 1992,24(3):329-343.

[10]蘇曉燕.采樣設(shè)計(jì)對(duì)土壤有機(jī)質(zhì)含量空間預(yù)測(cè)精度的影響研究[D].南京：南京師范大學(xué),2012.

[11]申祥民,雷曉云,陳大春,等.不同布點(diǎn)方式的膜下滴灌棉田土壤水分的空間變異研究[J].新疆農(nóng)業(yè)大學(xué)學(xué)報(bào),2010,33(4):363-368.

[12]姚麗賢,周修沖,蔡永發(fā),等.不同采樣密度下土壤特性的空間變異特征及其推估精度研究[J].土壤,2004,36(5):538-542.

[13]Visschers R, Finke P A, De Gruijter J J. A soil sampling program for the Netherlands[J]. Geoderma,2007,139(1):60-72.

[14]史志華,朱華德,陳佳,等.小流域土壤水分空間異質(zhì)性及其與環(huán)境因子的關(guān)系[J].應(yīng)用生態(tài)學(xué)報(bào),2012,23(4):889-895.

[15]白天路,楊勤科,申佳.黃土高原丘陵溝壑小流域土壤水分垂直分布變異特征及影響因子[J].生態(tài)雜志,2009,28(12):2508-2514.

[16]閆健康.GIS支持的平原地區(qū)鄉(xiāng)(鎮(zhèn))域村莊布點(diǎn)規(guī)劃研究[D].鄭州：河南大學(xué),2007.

[17]劉瑞民,楊志峰,沈珍瑤,等.基于DEM的長(zhǎng)江上游土地利用分析[J].地理科學(xué)進(jìn)展,2006,25(1):102-108.

[18]湯國安,宋佳.基于DEM坡度圖制圖中坡度分級(jí)方法的比較研究[J].水土保持學(xué)報(bào),2006,20(2):157-160.

[19]趙牡丹,湯國安.黃土丘陵溝壑區(qū)不同坡度分級(jí)系統(tǒng)及地面坡譜對(duì)比[J].水土保持通報(bào),2002,22(4):33-36.

A Study on Image-based Soil Moisture Sampling Method in Hilly and Gully Area of Loess Plateau

WANG Xiaojun, WU Jiangtao, WANG Bing, WEN Fenxiang

(TheCollegeofEnvironmentalScienceandResources,ShanxiUniversity,Taiyuan,Shanxi030006,China)

Abstract：[Objective] To test the rationality of soil sampling design, we took Sanyanjing basin of Shouyang county in Shanxi Province as a case study, to analyze the uniformity and representativeness of soil water sampling. [Methods] Using statistics analysis software and ArcGIS, we took grid sampling points as the initial layout. And then combined with remote sensing images, the automatic generated regular grid sampling points were finely adjusted giving the consideration of land use types and operation, namely, using uniform sampling in large scale and random sampling in small scale method for soil sampling point layout.[Results] Both of the nearest neighbor index and coefficient of variation showed a uniform distribution of sampling points. Sampling points obtained by uniform sampling in large scale and random sampling in small scale method were mainly concentrated in the slope of 0°~15°, which is not representative according to the distribution of slope. While the selected sampling points were representative according to the distribution of aspect. The points obtained by uniform sampling in large scale and random sampling in small scale method were more representative in grassland and farmland, but fewer representatives in forest land, due to the woodland is mainly distributed on steep slope and difficult to sample.[Conclusion] The uniformity and representativeness of sampling points obtained by large uniform small sample random sampling method are quite good. This sampling method can be well applied in hilly and gully area of the Loess Plateau to obtain uniform and representative sampling points.

Keywords:GIS; uniform sampling in large scale and random sampling in small scale; grid sampling; uniformity; representativeness

文獻(xiàn)標(biāo)識(shí)碼：B

文章編號(hào)：1000-288X(2015)01-0181-06

中圖分類號(hào)：S159.3

通信作者：王兵(1979—),女(漢族),河北省邢臺(tái)市人,博士,講師,主要從事農(nóng)田生態(tài)系統(tǒng)水氮利用方面的研究。E-mail:wang_bing1996@163.com。

收稿日期：2014-02-27修回日期：2014-03-26

資助項(xiàng)目：國家自然科學(xué)青年基金項(xiàng)目“黃土高原荒溪近自然治理的土壤—植被水文生態(tài)作用研究”(41201277)； 山西省科技重大專項(xiàng)“晉北沙化土地防治的關(guān)鍵技術(shù)研究與試驗(yàn)示范”(20121101011)； 山西省自然基金項(xiàng)目“黃土殘?jiān)珳羡謪^(qū)小流域土壤水分空間異質(zhì)性研究”(2014011034-2)

第一作者：王曉軍(1968—),男(漢族),山西省運(yùn)城市人,博士,副教授,主要從事自然資源評(píng)價(jià)與管理研究。E-mail:xjwang@sxu.edu.cn。