劉宏偉,高　菲,余鐘波,高風(fēng)華,5,向　龍

(1.南京水利科學(xué)研究院水文水資源研究所,江蘇 南京　210029;2.南京水利科學(xué)研究院水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室,江蘇 南京　210029;3.江蘇省水利廳,江蘇 南京　210029; 4.河海大學(xué)水文水資源學(xué)院,江蘇 南京　210098;5.南京工程學(xué)院環(huán)境工程學(xué)院,江蘇 南京　211167)

?

濕潤地區(qū)坡面土壤含水率時(shí)空變異性研究

劉宏偉1,2,高菲3,余鐘波4,高風(fēng)華4,5,向龍4

(1.南京水利科學(xué)研究院水文水資源研究所,江蘇 南京210029;2.南京水利科學(xué)研究院水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室,江蘇 南京210029;3.江蘇省水利廳,江蘇 南京210029; 4.河海大學(xué)水文水資源學(xué)院,江蘇 南京210098;5.南京工程學(xué)院環(huán)境工程學(xué)院,江蘇 南京211167)

通過野外觀測數(shù)據(jù)分析,研究濕潤地區(qū)坡面尺度土壤含水率的時(shí)空變異性。利用頻域反射(FDR)傳感器在1年內(nèi)以1次/h的觀測密度連續(xù)監(jiān)測近地表24個(gè)位置及2個(gè)剖面上12個(gè)位置的土壤含水率。利用變異系數(shù)和半方差圖方法,分析了研究區(qū)域內(nèi)不同空間位置上土壤含水率的時(shí)間變異性和不同季節(jié)土壤含水率的空間變異性。結(jié)果表明,在太湖流域地表附近土壤含水率呈中等或者中等偏弱變異性;土壤性質(zhì)、地形、植被、氣象和人類活動(dòng)等因素影響了土壤水分的時(shí)空變異性;不同要素在時(shí)間變異性和空間變異性上影響力不同,并且不同季節(jié)有不同的主導(dǎo)影響因素;不同季節(jié)土壤含水率半方差具有相似的模式,這與各因素作用的強(qiáng)弱有關(guān)。

土壤含水率;時(shí)空變異性;變異系數(shù);半方差圖;濕潤地區(qū);坡面

土壤水運(yùn)動(dòng)是水文循環(huán)中最重要的組成部分,土壤水分分布與運(yùn)動(dòng)規(guī)律的研究是水資源利用和環(huán)境整治的主要研究內(nèi)容之一。土壤水分運(yùn)動(dòng)與下滲過程、蒸散發(fā)過程以及地下水分運(yùn)動(dòng)過程聯(lián)系緊密[1-3],并對地表與地下徑流形成、溶質(zhì)遷移、土壤-大氣之間相互作用等各環(huán)節(jié)都有著顯著的影響[2-7]。土壤水分不論在空間上還是時(shí)間上,都有著很強(qiáng)的變異性[8-9]。土壤水分變異性的存在增加了水文過程及溶質(zhì)遷移過程研究的復(fù)雜程度。

土壤含水率在水文過程和溶質(zhì)遷移中具有重要作用,近年來關(guān)于土壤水分國內(nèi)外學(xué)者都進(jìn)行了大量研究。胡偉等[7]在黃土高原退耕坡地地區(qū)研究了土壤水分空間變異性及尺度影響,發(fā)現(xiàn)不同研究尺度上影響土壤含水率的主導(dǎo)因素和過程是不同的。Andrew等[8]研究了澳大利亞小流域中土壤含水率的空間相關(guān)性及其與流域水文過程的關(guān)系,認(rèn)為流域上土壤含水率在空間和時(shí)間上均有相關(guān)關(guān)系并受地形影響明顯。汪星等[10]在黃土高原半干旱區(qū)開展了山地密植棗林土壤水分特性研究,分析了該區(qū)域土壤水月變化規(guī)律及其與植被的關(guān)系。王云強(qiáng)等[11]研究了黃土高原區(qū)域尺度上土壤水分空間變異性規(guī)律,分析了變異系數(shù)和空間異質(zhì)比在0～500 cm垂向分布的規(guī)律性。一些學(xué)者[9, 12-14]比較了黃土區(qū)草地和農(nóng)地兩種土地利用類型條件下土壤含水率變異性,認(rèn)為土地利用方式對水平和垂向含水率分布有明顯影響。Lin等[15-16]在美國賓州的研究顯示,對土壤含水率空間變化和時(shí)間變化的研究有助于了解降雨發(fā)生時(shí)流域地表徑流的發(fā)生機(jī)制和地下優(yōu)勢流的產(chǎn)生現(xiàn)象,對準(zhǔn)確模擬和預(yù)測徑流過程十分重要。

從根本上影響土壤水變化的因素在時(shí)間和空間上是同時(shí)變化的,因此研究土壤水分變化特性必須同時(shí)從時(shí)空兩方面入手。目前國內(nèi)此類研究多見于對黃土高原半干旱區(qū)的研究,鮮有對濕潤地區(qū)土壤含水率方面探討。本文選取太湖流域西部宜興梅林小流域觀測數(shù)據(jù),比較不同季節(jié)土壤含水率的時(shí)間和空間變化,掌握太湖流域濕潤地區(qū)土壤水變化規(guī)律和主要影響因素,以幫助加深對濕潤地區(qū)產(chǎn)流、產(chǎn)污機(jī)制的理解。通過探求濕潤地區(qū)土壤含水率時(shí)空變異性規(guī)律和影響因素,為研究各因素如何定量影響水文循環(huán)過程打好基礎(chǔ),進(jìn)而為流域水文、溶質(zhì)模擬以及非點(diǎn)源污染遷移規(guī)律的研究和分布式水文、溶質(zhì)模擬的改進(jìn)提供科學(xué)依據(jù)。

1　研究區(qū)域和方法

1.1研究區(qū)域概況

研究區(qū)域梅林小流域地處宜興市東南(31°20′N, 119°51′E),面積約0.7 km2,位于太湖以西約9 km處,屬于封閉型小流域(圖1)。多年平均降雨量約1 150 mm,平均氣溫15.5℃。土地利用類型主要有菜地、旱地、板栗林、竹林、茶園、梨園和水稻田等。流域內(nèi)主要土壤類型為紅壤和黃棕壤,并存在少量白土。土壤水分長期觀測設(shè)備設(shè)置在流域內(nèi)一處山坡上,坡度約19°,高程變化在15.5～23.8 m之間(圖2)。觀測場處于流域內(nèi)靠近上游東向山坡的下部,坡面以下處有水塘,常年積水。坡面種植板栗樹,植被覆蓋率高,根系生長到地面下0.5 m以上。夏季地面有雜草生長。試驗(yàn)區(qū)坡長18 m,面積81 m2。土壤粒徑組成及其他特性見表1。在坡頂、中段及坡腳分樹下和樹間分別取3個(gè)樣品(表1)。

圖1　梅林流域位置

圖2　土壤水傳感器平面和剖面布置

試驗(yàn)區(qū)土壤取樣位置不同土壤質(zhì)地比例/%D<0.002mm0.002mm≤D≤0.05mmD>0.05mm干容重/(g·cm-3)飽和水力傳導(dǎo)度/(cm·d-1)坡頂10.967.022.11.3246.0坡中13.269.317.61.3731.9坡腳12.970.117.11.3437.5

注：D為土壤顆粒直徑。

1.2觀測布置和分析方法

1.2.1野外觀測及試驗(yàn)方法

在選定的試驗(yàn)區(qū)中,布置有34個(gè)頻域反射(FDR)土壤水分傳感器,其中，地表傳感器布置在沿山坡東線、中線和西線，分布在8個(gè)不同高度,同一高度3個(gè)傳感器間距1.5 m(如圖2中“FDR傳感器”);剖面?zhèn)鞲衅鞑贾迷贏、B兩個(gè)垂向剖面中(如圖2“測量剖面”),布置深度見圖2,具體為A剖面?zhèn)鞲衅魃疃葹榈乇怼⒌乇硪韵?.3 m、0.5 m、2.0 m和3.0 m 5個(gè)傳感器,B剖面?zhèn)鞲衅魃疃葹榈乇?、地表以?.3 m、0.5 m、1.0 m、1.5 m和2.0 m。試驗(yàn)場坡面地形中線位于坡脊,西、東兩線位于坡溝。坡脊上種有板栗樹,坡溝主要是樹間空地。土壤水觀測每1h測量1次電壓(單位mV)數(shù)據(jù),根據(jù)率定的“電壓-體積含水率曲線”換算體積土壤含水率。分析采用的數(shù)據(jù)從2006年11月16日到2007年11月7日,利用數(shù)據(jù)采集器逐小時(shí)記錄數(shù)據(jù)。土壤機(jī)械組成采用LS13320全自動(dòng)流光粒度分析儀測定,土壤容重采用環(huán)刀法測定,土粒密度采用比重瓶法測定[17]。雨量觀測使用翻斗式雨量計(jì)觀測。

1.2.2分析方法

利用不同季節(jié)土壤含水率變異系數(shù)和半方差圖來分析土壤含水率的時(shí)間變異性,以及不同季節(jié)的空間變異性,比較各種因素對土壤含水率變異性的影響[18-19]。在分析土壤含水率觀測值不同空間位置的時(shí)間變異性和不同時(shí)間的空間變異性時(shí),參考了變異系數(shù)Cv[20]。當(dāng)Cv≤0.1時(shí)為弱變異性,當(dāng)0.1

半方差圖法是用來比較土壤水分的空間連續(xù)特性和區(qū)域變異特性的。假定土壤屬性在空間上二階平穩(wěn),用半方差函數(shù)來表示其空間相關(guān)結(jié)構(gòu)。半方差常采用式(1)計(jì)算:

(1)

式中:h為位差;N(h)為位差為h時(shí)的樣本對的數(shù)量;θ(xi)為xi位置的土壤含水率;n為樣本數(shù)量。半方差圖中半方差值隨步長增加到一個(gè)相對穩(wěn)定的水平上時(shí),對應(yīng)的半方差值稱為基臺值,這是區(qū)域化變量總體特征的體現(xiàn)。變程是指變異函數(shù)在有限步長上達(dá)到基臺值時(shí)對應(yīng)的步長,也稱自相關(guān)距離,是空間自相關(guān)性的最大距離。本文中半方差的計(jì)算使用地學(xué)統(tǒng)計(jì)程序庫(geostatistical software library,GSLIB)[19],分別從總體變異性、土壤特性影響、地形影響、植被影響、氣象因素影響等方面進(jìn)行分析。

2　結(jié)　果

2.1地表附近土壤含水率變異性

2.1.1沿坡長方向不同位置土壤含水率時(shí)間變異性

對沿坡長8個(gè)不同高度的平均土壤含水率時(shí)間序列進(jìn)行分析(圖3,表2,其中1～8位置為圖2中所標(biāo)注的不同高度),其變異系數(shù)均在0.20～0.26之間,變異性屬中等偏弱。沿坡長方向,體積含水率年平均值總體趨勢為上部較小,下部較大,而實(shí)測數(shù)據(jù)的時(shí)間序列變異系數(shù)則變化不大。坡頂和坡長中段6號位置的變異系數(shù)略高于其他位置的測量值,其時(shí)間變異性略強(qiáng)于其他位置。

圖3　不同季節(jié)坡面含水率時(shí)間序列變異系數(shù)

沿坡位置x/%σXCv124.815.920.24229.606.650.22331.486.240.20435.587.150.20535.087.280.21632.518.330.26730.476.790.22831.466.360.20

2.1.2垂直坡長方向土壤含水率時(shí)間變異性

表3為3條沿線上平均土壤含水率時(shí)間序列變異系數(shù)等統(tǒng)計(jì)數(shù)據(jù)。中線的土壤含水率平均值高于西東兩線,而最小值低于左右兩邊。3條測線土壤含水率的最大值非常接近。時(shí)間序列的變異系數(shù)顯示,3條線土壤含水率平均值均為中等偏弱變異性,中線變異性相對東西線較強(qiáng)。中線的全年土壤含水率平均值和最小值較低,變異系數(shù)較大,與東西兩線相比,體現(xiàn)出明顯的坡脊與坡溝土壤含水率變化的差異。根據(jù)表4不同季節(jié)的統(tǒng)計(jì)結(jié)果來看,冬春夏3季的土壤含水率時(shí)間變異系數(shù)依次略有升高,差別不太大,而秋季則有明顯上升,土壤含水率平均值則有與之相應(yīng)的相反結(jié)果。

表3　垂直坡長土壤含水率時(shí)間序列統(tǒng)計(jì)結(jié)果

表4　垂直坡長土壤含水率季節(jié)序列統(tǒng)計(jì)結(jié)果

2.1.3不同季節(jié)地表土壤含水率空間變異性

從不同季節(jié)平均土壤含水率在空間上的變異性來看,多屬于中等或者中等偏弱的變異性。坡面上3條測線沿坡長方向的空間變異系數(shù)較小,其中西線與東線變異性小于中線(圖4(b))。不同高度上垂直坡長方向的土壤含水率的空間變異性高低相間,呈波浪型變化(圖4(a)),并且坡頂和坡腳兩處垂直坡向的變異系數(shù)高于中間各個(gè)高度。從不同季節(jié)的空間變異性來看,春、夏、秋3季的變化規(guī)律相似,并且與冬季空間變異系數(shù)相比變化略大。

圖4　不同季節(jié)土壤含水率變異系數(shù)

圖5　不同季節(jié)土壤含水率半方差擬合

根據(jù)沿坡長方向半方差計(jì)算結(jié)果,可以發(fā)現(xiàn)春、夏、秋3個(gè)季節(jié)的半方差圖模式相似,可以用球狀模型擬合,冬季的變程明顯大于其他3個(gè)季節(jié),用高斯模型可以很好擬合(圖5)。春、夏、秋3季半方差變程約為8 m,而冬季的半方差在16 m以內(nèi)都沒有達(dá)到基臺值。這說明在春、夏、秋季沿坡長方向土壤含水率的相關(guān)距離要小于冬季。

2.2垂向剖面土壤含水率變異性

2.2.1垂直剖面土壤含水率時(shí)間序列變異性

沿垂直方向的土壤含水率時(shí)間變異性從淺到深有明顯的從大到小的變化趨勢(圖6)。地表土壤含水率的時(shí)間變異性為中等,而較深的土壤變異性逐漸減弱,逐漸變成弱變異性。變異系數(shù)沿深度變化存在轉(zhuǎn)折點(diǎn),位于0.5 ～ 1.0 m深度。這種趨勢的分布與降水、蒸發(fā)、植被生長和根系深度、地下水運(yùn)動(dòng)有著密切的關(guān)系。

圖6　A、B剖面土壤含水率時(shí)間序列變異系數(shù)

2.2.2不同季節(jié)垂向剖面土壤含水率空間變異性

從不同季節(jié)土壤含水率垂向剖面分布情況來看,可以將剖面分成3個(gè)部分(圖7)。第一部分為近地表層,包括A剖面的地表到0.5 m深度和B剖面上從地表到1.5 m與2.0 m之間。這部分土壤含水率在不同季節(jié)變化明顯,受前面分析的各個(gè)因素影響明顯。第二部分為穩(wěn)定層,包括A剖面0.5～1.5 m之間和B剖面2.0～3.0 m之間。這部分土壤含水率時(shí)間變化非常小,對地表水和地下水變化響應(yīng)均不明顯。第三部分為近潛水層,包括A剖面2.0 m上下,B剖面3.0 m上下。這部分受變動(dòng)地下潛水位的影響,當(dāng)潛水位在這個(gè)層間隨時(shí)間變動(dòng)時(shí)會造成土壤含水率的變化。A剖面位于坡腳,而B剖面位于坡頂。與坡腳剖面相比,坡頂剖面含水率變異性明顯較強(qiáng)。

圖7　季節(jié)平均土壤含水率沿垂向分布

3　討　論

濕潤丘陵地區(qū)土壤含水率不論在時(shí)間還是空間上均有中等或者中等偏弱的變異性。不同空間位置上的時(shí)間變異系數(shù)和不同時(shí)間段的空間平均土壤含水率變異系數(shù),多數(shù)情況下在0.13～0.26之間,少數(shù)情況出現(xiàn)高于0.4或者低于0.1。

3.1地表土壤水時(shí)間和空間變異性的影響因素

從土壤特性影響方面來看,坡頂土壤含沙量高(表1),平均飽和水力傳導(dǎo)度較高,透水性好于其他位置,在土壤含水率時(shí)間變化過程中,坡頂吸收和釋放水分快,因此坡頂含水率的時(shí)間變異性與其他位置相比略大。在相同高度的坡面上,樹下和樹間空地的土壤性質(zhì)略有差別[21],但總體差別不大。實(shí)測數(shù)據(jù)反映出3條測線上的飽和土壤含水率非常相似(表3),土壤特性對垂直坡長方向上土壤含水率的變異性影響不大。

從地形影響方面來看，降雨發(fā)生以后,坡頂水分在重力作用下將會向坡腳方向匯集,使得坡面上部土壤水分時(shí)間序列變異系數(shù)略高于中下部。位于西線和東線沿坡長方向的空間變異系數(shù)明顯小于中線(圖4(b))。這樣的現(xiàn)象首先是由于實(shí)驗(yàn)區(qū)局部微地形的影響,西東兩線位于坡溝,中線位于坡脊。相比之下,坡脊排水容易,而坡溝則更有利保持水分,所以坡脊土壤水分變化相比之下更劇烈。

從季節(jié)和天氣影響方面來看,研究區(qū)域全年都有降水分布,在夏季最為集中。從不同季節(jié)的時(shí)間系列變異系數(shù)來看,冬季變異性最弱,因?yàn)榻涤旰驼舭l(fā)量都較小。8—10月的秋季氣溫高,蒸發(fā)與植物散發(fā)量大,土壤含水率均值顯然低于其他3季,3個(gè)月觀測值的時(shí)間變異性也明顯高于其他3季,這主要受氣象因素影響。夏季坡面不同高度測量值的時(shí)間變異系數(shù)差別較大,總體呈上高下低。造成這種現(xiàn)象,主要是因?yàn)橄募鹃L時(shí)間降雨,下部測點(diǎn)很長時(shí)間處于接近飽和或者田間持水量的狀態(tài),坡上處于這種狀態(tài)的時(shí)間則短得多,因此其時(shí)間變異性會小于坡面上部的位置。即降水發(fā)生多的時(shí)段,坡頂土壤含水率變異性最大;雨期結(jié)束后,坡頂變異性逐漸接近坡腳的狀態(tài)。Famiglietti等[22]在沿坡長方向7個(gè)月的采樣分析也得出類似的觀測結(jié)果。

植被的影響在坡面不同高度上垂直坡長方向的空間變異性也有明顯的體現(xiàn)。如圖4(a)所示,變異系數(shù)呈波浪形變化,有板栗樹生長的高度垂直坡長方向的空間變異性較強(qiáng),無植被生長的高度位置則較弱。中線測量傳感器沿坡面是按樹下—樹間—樹下間隔布置,變異系數(shù)波浪形變化規(guī)律與之相應(yīng)。圖4(b)顯示夏季中線變異性最強(qiáng),春秋次之,冬季最弱,這與植物生長的影響關(guān)系密切。植被生長會導(dǎo)致土壤持水能力和蒸散發(fā)能力發(fā)生變化,從而導(dǎo)致土壤含水率變異性較強(qiáng)。

此外,人類活動(dòng)也對土壤水的變化有著明顯的影響。坡面上6號位置測值的時(shí)間序列變異系數(shù)明顯高于周圍(表2),主要原因是測區(qū)的人類活動(dòng)造成擾動(dòng),導(dǎo)致變異性增強(qiáng)。

與冬季相比,春、夏、秋3季降雨多且頻繁,植被生長也較旺盛,因此植被、氣象因素的影響強(qiáng)于地形因素引起的土壤水分沿山坡再分布的作用,而冬季這種土壤水分再分布的運(yùn)動(dòng)影響明顯得多。土壤水分再分布作用聯(lián)系著坡面上很長距離上的土壤含水率。冬季其他因素影響減弱,使這種再分布作用成為主要影響因素之一,從而造成不同時(shí)間段沿坡面方向土壤含水率相關(guān)距離,冬季大于其他季節(jié)的現(xiàn)象。因此,從不同季節(jié)土壤含水率變化模式可以識別出植被、氣象、地形等不同影響因素的不同作用效果。若要進(jìn)一步弄清各影響因素的作用機(jī)理和定量規(guī)律,則有必要開展更加可控的針對性實(shí)驗(yàn)。

3.2垂直剖面土壤含水率變異性

首先,越上層的土壤接受降水補(bǔ)給相對越容易,土壤水分對降水的響應(yīng)越明顯,造成土壤水分時(shí)間序列變化劇烈。其次,在無雨期越接近地表的土層受太陽輻射、氣溫等氣象因素影響越大,地表的蒸發(fā)過程與深層土壤相比也更為活躍,這造成土壤含水率變異性增大。第三,植被根系的生長影響較大,A剖面位于坡腳的樹間區(qū),地表植被為較矮小的雜草和茶樹,其根系生長深度一般為0.3～0.5 m。在這個(gè)范圍之內(nèi),植物生長對土壤水分影響明顯,根系對水分的吸收造成了0～ 0.5 m一段內(nèi)土壤水分的較強(qiáng)變異性。B剖面布置的位置為板栗樹附近。板栗樹的根系比雜草和茶樹深得多,所以從圖6(b)可以發(fā)現(xiàn),土壤含水率變異性較高的區(qū)域比A剖面(圖6(a))要深一些。由于試驗(yàn)區(qū)位于濕潤地區(qū),地下潛水埋深很淺,所以觀測剖面均受到變動(dòng)地下水位的影響。從圖6可以發(fā)現(xiàn),A、B兩剖面最深處都有變異性反而增加的現(xiàn)象,A剖面為1.5 m,B剖面為3.0 m。這一現(xiàn)象主要是因?yàn)樵诓糠重S水時(shí)間里,潛水面達(dá)到測點(diǎn)以上,而缺水時(shí)間里,潛水面則退下去,這樣就造成了土壤含水率時(shí)間的變動(dòng)性有所增加。

4　結(jié)　論

a. 濕潤丘陵地區(qū)土壤含水率不論在時(shí)間還是空間上均有中等或者中等偏弱的變異性。時(shí)間變異性總體趨勢為冬季變異性最弱,春、夏季比冬季強(qiáng),秋季最強(qiáng)。空間變異性總體趨勢為坡頂變異性強(qiáng)于坡腳;坡脊變異性強(qiáng)于坡溝;垂直剖面上,地表層和受變動(dòng)地下水位影響大的深層,土壤含水率變異性強(qiáng)于中等深度土層。

b. 土壤含水率變化受到土壤特性、地形、植被、氣象和人類活動(dòng)等因素的影響明顯。各個(gè)影響因子在不同時(shí)間和空間條件有著不同的影響強(qiáng)度。在小尺度的山坡試驗(yàn)區(qū)上,不同空間位置不同的土壤特性主要影響各空間位置上土壤含水率在時(shí)間上的變異性。地形因素對沿坡長、垂直坡長和剖面3個(gè)維度均有影響。主要由于細(xì)節(jié)地形的差異影響匯流路徑和土壤水分再分布的方式影響土壤水分變異性。植被生長在不同季節(jié)有不同的影響,春、夏、秋強(qiáng),冬季弱。同時(shí),植被生長旺盛的位置土壤水分變異性越強(qiáng)。氣象因素包括降雨、蒸發(fā)等,是土壤含水率變化的最主要的驅(qū)動(dòng)因素。雨多、氣溫高、植被生長旺盛的時(shí)間段氣象因素影響最強(qiáng)烈。人類活動(dòng)對地表土壤含水率的影響十分明顯并且隨機(jī)性強(qiáng)。

c. 沿坡長方向地表土壤含水率的空間變異性按不同時(shí)間段有不同的變異模式。當(dāng)植被、氣候等外部因素影響較顯著時(shí),土壤含水率半方差具有相似的模式,如試驗(yàn)采樣期內(nèi)的春、夏、秋季半方差圖十分相似,都具有相似的基臺值和變程。當(dāng)外部因素影響不顯著,地形等因素引起的土壤水分再分布占主要影響地位時(shí),則會出現(xiàn)另一種模式,冬季的半方差圖變程超過16.0 m,變異性強(qiáng)于其他季節(jié)。

[ 1 ] IGNACIO R I,GREGOR K.V,RICCARDO R.On the spatial organization of soil moisture fields[J].Geophysical Research Letters,1995,22(20)：2757-2760.

[ 2 ] YU Z,CARLSON T N,BARRON E J,et al.On evaluating the spatial-temporal variation of soil moisture in the Susquehanna River Basin[J].Water Resources Research,2001,34:1313-1326.

[ 3 ] DONG W,YU Z,WEBER D.Simulation on soil water variation in arid regions[J].Journal of Hydrology,2003,275:162-181

[ 4 ] ANDREW W,WESTERN,ZHOU Senlin,et al.Spatial correlation of soil moisture in small catchments and its relationship to dominant spatial hydrological processes[J].Journal of Hydrology,2004,286:113-134.

[ 5 ] 劉宏偉,余鐘波,崔廣柏.濕潤地區(qū)土壤水分對降雨的響應(yīng)模式研究[J].水利學(xué)報(bào),2009,40(7):822-829.(LIU Hongwei,YU Zhongbo,CUI Guangbo.Pattern of soil moisture responding to precipitation in humid area[J].Journal of Hydraulic Engineering,2009,40(7):822-829.(in Chinese))

[ 6 ] 劉宏偉,余鐘波.用Hargreaves法與Penman-Monteith法計(jì)算ET0:以太湖流域的應(yīng)用為例[J].水資源保護(hù),2010,26(1):6-8.(LIU Hongwei,YU Zhongbo.Application of Hargreaves and Penman-Monteith Equation to estimatingET0: a case study in Taihu Basin[J].Water Resources Protection,2010,26(1):6-8.(in Chinese))

[ 7 ] 胡偉,邵明安,王全九.黃土高原退耕土壤水分空間變異性研究[J].水科學(xué)進(jìn)展,2006,17(1):74-81.(HU Wei,SHAO Ming’an,WANG Quanjiu.Study on spatial variability of soil moisture on the recultivated slope-land on the Loess Plateau[J].Advances in Water Science,2006,17(1):74-81.(in Chinese))

[ 8 ] ANDREW W,RODGER B G,GUNTER B,et al.Observed spatial organization of soil moisture and its relation to terrain indices[J].Water Resources Research,1999,35(3): 797-810.

[ 9 ] 王軍,傅伯杰,邱揚(yáng),等.用空間內(nèi)插法研究黃土丘陵小流域土壤水分時(shí)空分布特征[J].自然科學(xué)進(jìn)展,2002,12(4): 430-433.(WANG Jun,FU Bojie,QIU Yang,et al.Study on soil moisture of a small watershed on the Loess Plateau with spatial interpolative methods[J].Advances in Natural Science,2002,12(4):430-433.(in Chinese))

[10] 汪星,周玉紅,汪有科,等.黃土高原半干旱區(qū)山地密植棗林土壤水分特性研究[J].水利學(xué)報(bào),2015,46(3):263-270.(WANG Xing,ZHOU Yuhong,WANG Youke,et al.Soil water characteristic of a dense jujube plantation in the semi-arid hilly regions of the Loess Plateau in China[J].Journal of Hydraulic Engineering,2015,46(3):263-270.(in Chinese))

[11] 王云強(qiáng),邵明安,劉志鵬.黃土高原區(qū)域尺度土壤水分空間變異性[J].水科學(xué)進(jìn)展,2012,23(3):310-316.(WANG Yunqiang,SHAO Ming’an,LIU Zhipeng.Spatial variability of soil moisture at a regional scale in the Loess Plateau[J].Advances in Water Science,2012,23(3):310-316.(in Chinese))

[12] 王云強(qiáng),張興昌,從偉,等.黃土區(qū)不同土地利用方式坡面土壤含水率的空間變異性研究[J].農(nóng)業(yè)工程學(xué)報(bào),2006,22(12):65-71.(WANG Yunqiang,ZHANG Xingchang,CONG Wei,et al.Spatial variability of soil moisture on slope-land under different land uses on the Loess Plateau[J].Transactions of the CSAE,2006,22(12): 65-71.(in Chinese))

[13] 王軍,傅伯杰.黃土丘陵小流域土地利用結(jié)構(gòu)對土壤水分時(shí)空分布的影響[J].地理學(xué)報(bào),2000,55(1):84-91.(WANG Jun,FU Bojie.The impact of land use on spatial and temporal distribution of soil moisture on the Loess Plateau[J].Acta Geographica Sinica,2000,55(1):84-91.(in Chinese))

[14] QIU Yang,FU Bojie,WANG Jun,et al.Soil moisture variation in relation to topography and land use in a hillslope catchment of the Loess Plateau,China[J].Journal of Hydrology,2001,240(3/4): 243-246.

[15] LIN H S,KOGELMANN W,WALKER C,et al.Soil moisture patterns in a forested catchment:a hydropedological perspective[J].Geoderma,2006,131:345-368.

[16] LIN H,ZHOU X.Evidence of subsurface preferential flow using soil hydrologic monitoring in the Shale Hills catchment[J].European Journal of Soil Science,2008,59:34-49.

[17] 魯如坤.土壤農(nóng)業(yè)化學(xué)分析方法[M].北京：中國農(nóng)業(yè)出版社,2000:266-295.

[18] 馮益明,唐守正,李增元.空間統(tǒng)計(jì)分析在林業(yè)中的應(yīng)用[J].林業(yè)科學(xué)，2004,40(3):149-155.(FENG Yiming,TANG Shouzheng,LI Zengyuan.Application of spatial statistic analysis in forestry[J].Scientia Silvae Sinicae,2004,40(3):149-155.(in Chinese))

[19] DEUTSCH C V,JOURNAL A G.GSLIB geostatistical software library and user’s guide[M].2nd Edtion.Oxford:Oxford University Press,1998.

[20] 王孟本,李宏建.晉西北黃土區(qū)人工林土壤水分動(dòng)態(tài)的定量研究[J].生態(tài)學(xué)報(bào),1995,15(2):178-184.(WANG Mengben,LI Hongjian.Quantitative study on the soil water dynamics of various forest plantations in the Loess Plateau region in northwestern Shanxi[J].Acta Ecol Sin,1995,15(2):178-184.(in Chinese))

[21] YIN J,Young M,YU Z.Effects of paleoclimate and time-varying canopy structures on paleowater fluxes[J].Journal of Geophysical Research,2008,113(D6):304-312.doi:10.1029/2007JD009010.

[22] FAMIGLIETTI J S,RUDNICKI J W,RODELL M.Variability in surface moisture content along a hillslope transect: Rattlesnake Hill,Texas[J].Journal of Hydrology,1998,210:259-281.

Study on temporal-spatial variability of soil moisture content on hillslope in a humid area

LIU Hongwei1, 2, GAO Fei3, YU Zhongbo4, GAO Fenghua4, 5, XIANG Long4

(1.HydrologyandWaterResourcesDepartment,NanjingHydraulicResearchInstitute,Nanjing210029,China;2.StateKeyLaboratoryofHydrology-WaterResourcesandHydraulicEngineering,>NanjingHydraulicResearchInstitute,Nanjing210029,China;3.WaterResourcesDepartmentofJiangsuProvince,Nanjing210029,China;4.CollegeofHydrologyandWaterResources,HohaiUniversity,Nanjing210098,China;5.SchoolofEnvironmentalEngineering,NanjingInstituteofTechnology,Nanjing211167,China)

The temporal-spatial variability of soil moisture content on a hillslope in a humid area was analyzed based on field monitoring data. Hourly data over a year were collected with frequency domain reflectometry (FDR) sensors. Twenty-four sensors near the soil surface and 12 sensors in two profiles were set up to monitor the soil moisture content continuously. The coefficient of variation and semi-variogram were used to analyze the soil moisture content’s temporal variability in different locations and the spatial variability in different seasons in the study area. The results show that, in the Taihu Basin, the surface soil moisture content exhibited a medium or low level of variability. Soil characteristics, topography, vegetation, meteorological factors, and human activities influenced the temporal-spatial variability of soil moisture. These factors had different effects on the temporal-spatial variability, and the dominant factors were different in each of the four seasons. The semi-variograms for the soil moisture contents in different seasons had similar patterns due to different effects of the factors.

soil moisture content; temporal-spatial variability; coefficient of variation; semi-variogram; humid area; hillslope

10.3880/j.issn.1004-6933.2016.05.004

國家自然科學(xué)基金青年項(xiàng)目(51209139)

劉宏偉(1982—)，男，高級工程師，博士，主要從事水文實(shí)驗(yàn)、模擬以及防洪減災(zāi)等研究。E-mail: hwliu@nhri.cn

P33

A

1004-6933(2016)05-0017-07

2015-10-15編輯：彭桃英)