王俊波 漆喜林 李戰(zhàn)剛 高海銀 王耀軍 郭 峰 唐翠平

(1.陜西省防護(hù)林建設(shè)工作站，陜西 西安 710082；2.定邊縣林業(yè)工作站，陜西 定邊 718600；3.吳起縣種苗站，陜西 吳起 717600；4.西南林業(yè)大學(xué)林學(xué)院，云南 昆明 650224)

黃土高原油松正常生長(zhǎng)的最低需水量估計(jì)

王俊波1漆喜林1李戰(zhàn)剛1高海銀2王耀軍3郭 峰4唐翠平4

(1.陜西省防護(hù)林建設(shè)工作站，陜西 西安 710082；2.定邊縣林業(yè)工作站，陜西 定邊 718600；3.吳起縣種苗站，陜西 吳起 717600；4.西南林業(yè)大學(xué)林學(xué)院，云南 昆明 650224)

以13年生油松人工林為研究對(duì)象，采用樹高生長(zhǎng)模擬、逐步回歸、生長(zhǎng)指數(shù)等方法，考慮降水和整地蓄水作用，估計(jì)黃土高原地區(qū)油松正常生長(zhǎng)的最低需水量，結(jié)果表明：油松正常生長(zhǎng)的最低需水量為540 mm，并據(jù)此計(jì)算出不同降雨量下的集水區(qū)大小。

油松；生長(zhǎng)指數(shù)；蓄水量；需水量；黃土高原

油松(PinustabulaeformisCarr.)是我國(guó)北方地區(qū)的主要造林樹種之一，在暖溫帶及中溫帶的濕潤(rùn)、半濕潤(rùn)、半干旱、干旱區(qū)均有人工栽培。在這一廣大的分布區(qū)內(nèi)，水熱條件差別造成了林地生產(chǎn)力的地域性分異[1-5]。其中值得注意的是在黃土高原地區(qū)，干旱少雨和土壤干化使其形成了大面積的低質(zhì)、低效林[6-7]。因此，油松生長(zhǎng)狀況與水分條件的關(guān)系研究引起了眾多學(xué)者的興趣。王梅等、焦醒等的研究表明，陰坡油松的生長(zhǎng)量顯著大于陽(yáng)坡，其原因在于陰坡的水分條件優(yōu)于陽(yáng)坡[8-9]。靳天恩等通過回歸分析，建立了油松樹高生長(zhǎng)量與降水量、蒸發(fā)量、平均溫度等因子的回歸關(guān)系模型[10]。丁曉剛等通過逐步回歸分析，篩選出影響油松樹高生長(zhǎng)的主導(dǎo)因子，即降水量[11]。但在自然狀態(tài)下，黃土高原地區(qū)油松正常生長(zhǎng)的最低需水量到底是多少，明確回答這一問題，對(duì)當(dāng)?shù)赜退稍炝衷O(shè)計(jì)及油松林地水分管理和密度調(diào)控具有一定的理論意義和參考價(jià)值。然而，根據(jù)目前的研究結(jié)果尚無(wú)法對(duì)這一問題做出明確回答。針對(duì)這一情況，本文以樹高生長(zhǎng)指數(shù)為因變量，以相應(yīng)年份的氣象資料為自變量，通過逐步篩選建立樹高生長(zhǎng)指數(shù)與降水量的回歸模型，根據(jù)回歸模型以及整地蓄水作用估計(jì)出油松正常生長(zhǎng)的最低需水量。

1 研究區(qū)概況

調(diào)查地點(diǎn)位于陜西省榆林市榆陽(yáng)區(qū)，地貌為黃土丘陵，以覆沙黃土為主；氣候?qū)贉貛О敫珊荡箨懶约撅L(fēng)氣候，年平均氣溫8.1 ℃，極端最低氣溫-32.7 ℃，極端最高氣溫38.6 ℃；年均降水量438 mm，降水主要集中在7—9月，約占全年降水量的70%；年均風(fēng)速5.3 m/s，最大風(fēng)速28.0 m/s，風(fēng)向以西北為主。樣地設(shè)在山坡中部，海拔 1 142 ～ 1 150 m。在所調(diào)查的3塊樣地中，造林技術(shù)、經(jīng)營(yíng)措施及立地條件一致。其中，造林時(shí)間為1998年10月，株行距4 m×4 m；魚鱗坑整地，規(guī)格為直徑80 cm、深度30 cm。調(diào)查于2011年7月進(jìn)行，林分年齡13 a。

2 研究方法

樣地選擇采用典型抽樣法，即堅(jiān)持代表性、典型性和可比性原則，在全面踏查的基礎(chǔ)上選擇具有代表性和典型性的地段設(shè)置樣地。樣帶面積為100 m×8 m，樣地調(diào)查采用每木檢尺法，即逐株測(cè)定胸徑、樹高、冠幅生長(zhǎng)量，每塊樣帶調(diào)查株數(shù)>50株。每塊樣地選擇3株平均標(biāo)準(zhǔn)木測(cè)定逐年樹高生長(zhǎng)量，標(biāo)準(zhǔn)木平均樹高為4.15～4.43 m、胸徑為9.15～9.46 cm、冠幅為2.83～3.41 m。每塊樣地按對(duì)角線采集3個(gè)土樣，方法為“環(huán)刀法”和“鋁盒法”，采樣深度30～50 cm，土壤含水率及容重的測(cè)定采用“烘箱法”。

采用Excel及SPSS 16.0進(jìn)行數(shù)據(jù)分析。其中，影響生長(zhǎng)的主導(dǎo)氣象因子篩選采用逐步回歸法；樹高生長(zhǎng)過程采用邏輯斯蒂方程與自然對(duì)數(shù)函數(shù)相結(jié)合的四參數(shù)理論生長(zhǎng)復(fù)合模型進(jìn)行擬合，表達(dá)式如下：

(1)

式中：lny為累積生長(zhǎng)量的自然對(duì)數(shù)；t為生長(zhǎng)時(shí)間(年齡)；b0為常數(shù)；b1、b2、b3均為待定參數(shù)，采用麥夸特算法估計(jì)[12]。

3 結(jié)果與分析

3.1 樹高生長(zhǎng)過程擬合

根據(jù)模型(1)對(duì)各樣地標(biāo)準(zhǔn)木樹高生長(zhǎng)過程進(jìn)行回歸分析，結(jié)果見表1。

表1 樹高生長(zhǎng)過程擬合及其顯著性檢驗(yàn)

由表1可知，樹高生長(zhǎng)量與年齡之間的回歸關(guān)系均達(dá)到極顯著水平(p<0.01)，決定系數(shù)(R2)為0.990 3～0.998 5，表明該模型可解釋樹高生長(zhǎng)過程變異的99.0%以上。因此，采用該模型可對(duì)不同年齡的理論生長(zhǎng)量進(jìn)行預(yù)測(cè)。

3.2 樹高生長(zhǎng)指數(shù)估計(jì)

生長(zhǎng)指數(shù)是不同年齡實(shí)際生長(zhǎng)量與同期理論生長(zhǎng)量之比，可以反映環(huán)境因素變動(dòng)對(duì)樹木生長(zhǎng)的影響[13]。根據(jù)表1中的模型預(yù)測(cè)不同年齡的理論生長(zhǎng)量，估計(jì)不同年份油松樹高生長(zhǎng)指數(shù)，結(jié)果見表2。

表2 樹高生長(zhǎng)指數(shù)及主要?dú)庀笠蜃咏y(tǒng)計(jì)

由表2可知，2002、2003、2007、2008年的樹高生長(zhǎng)指數(shù)大于1.0，實(shí)際生長(zhǎng)量大于理論生長(zhǎng)量；2004、2005、2006、2009、2010年的樹高生長(zhǎng)指數(shù)小于1.0，實(shí)際生長(zhǎng)量小于理論生長(zhǎng)量。

3.3 影響樹高生長(zhǎng)指數(shù)的主導(dǎo)因子篩選

由于調(diào)查對(duì)象位于同一地點(diǎn)，且造林技術(shù)、經(jīng)營(yíng)措施、立地條件相同，因此樹高生長(zhǎng)指數(shù)年際之間的變化歸因于氣象因子的波動(dòng)。為此，以不同年份的樹高生長(zhǎng)指數(shù)為因變量、以相應(yīng)年份的氣象因子為自變量(表2)，采用逐步回歸法篩選影響油松生長(zhǎng)的主導(dǎo)因子，結(jié)果如下：

Ih= 0.497 99+0.001 24x

(R=0.895 0，p<0.001)

(2)

式中：Ih為樹高生長(zhǎng)指數(shù)；x為年降水量。

由此表明：在影響油松生長(zhǎng)的氣象因子中，年降水量是主導(dǎo)因子或限制因子。而且，油松樹高生長(zhǎng)指數(shù)與年降水量之間存在極顯著的回歸關(guān)系，年降水量可解釋油松樹高生長(zhǎng)指數(shù)變異的80.10%(決定系數(shù)R2)。

3.4 油松正常生長(zhǎng)的最低需水量估計(jì)

根據(jù)回歸模型(2)估計(jì)，當(dāng)樹高生長(zhǎng)指數(shù)為1.0(實(shí)際生長(zhǎng)量與理論生長(zhǎng)量相等)時(shí)，油松生長(zhǎng)所需降水量為405.0 mm。與表2比較可以發(fā)現(xiàn)：降水量為405.0 mm時(shí)，樹高生長(zhǎng)指數(shù)并不能穩(wěn)定通過1.0。為此，以樹高生長(zhǎng)指數(shù)為1.05作為油松正常生長(zhǎng)的指標(biāo)，即實(shí)際生長(zhǎng)量超過理論生長(zhǎng)量的5%，則所需的降水量為445.0 mm。

另一方面，根據(jù)回歸方程僅估計(jì)了油松正常生長(zhǎng)所需的天然降水量，忽略了通過整地?cái)r蓄的地表徑流部分。為了比較準(zhǔn)確地估計(jì)油松正常生長(zhǎng)所需的最低需水量，必須考慮整地?cái)r蓄的地表徑流數(shù)量。為此，測(cè)定了集水區(qū)內(nèi)、外的土壤含水率。結(jié)果表明：集水區(qū)內(nèi)土壤含水率為13.23%，集水區(qū)外為11.10%，兩者存在極顯著差異(p=0.002)。在此情況下，根據(jù)土壤含水率、集水區(qū)規(guī)格、土壤容重即可估計(jì)出整地?cái)r蓄的地表徑流數(shù)量。其中，集水面半徑(r)為40 cm、深度(h)為30 cm，則土壤蓄水體積為0.075 4 m3；水分儲(chǔ)量=蓄水體積×土壤容重×土壤含水率。通過測(cè)定得知土壤容重為 1 500 kg/m3，則集水區(qū)內(nèi)土壤的水分儲(chǔ)量為14.96 kg、集水區(qū)外土壤的水分儲(chǔ)量為12.55 kg，即通過整地在集水區(qū)可多攔蓄地表徑流2.41 kg。由于集水區(qū)的面積為 0.251 2 m2，因此，通過整地可提高蓄水量9.59 kg/m2，相當(dāng)于9.59 mm降水。然而，本次含水率測(cè)定在2011年7月2日進(jìn)行，上述推算結(jié)果屬一次性的蓄水量狀況。根據(jù)多年的測(cè)定結(jié)果，榆陽(yáng)區(qū)6月份的平均降水量約為40.0 mm，而年均降水量400.0 mm左右[14-15]。據(jù)此估計(jì)，集水區(qū)全年攔蓄的地表徑流為95.90 mm。

由上述分析可知：油松正常生長(zhǎng)(樹高生長(zhǎng)指數(shù)大于1.05)的最低需水量=天然降水+整地蓄水=445.0+96.0，約540 mm。

4 結(jié)論與討論

在研究區(qū)，影響油松生長(zhǎng)的主導(dǎo)因子是水分，降水不足限制了油松的正常生長(zhǎng)。根據(jù)生長(zhǎng)指數(shù)與降水量的回歸模型及整地蓄水作用估計(jì)，13年生油松林正常生長(zhǎng)的最低需水量為540 mm。在本文統(tǒng)計(jì)的10年期間，降水較多的2002年、2003年、2007年、2008年，天然降水與整地蓄水之和介于532.2～622.6 mm，樹高生長(zhǎng)指數(shù)平均為1.096，油松生長(zhǎng)正常；降水中等的2004年、2009年，天然降水與整地蓄水之和分別約為516.0、516.8 mm，樹高生長(zhǎng)指數(shù)平均為0.912，油松生長(zhǎng)受到一定的抑制；降水較少的2005年、2006年、2010年，天然降水與整地蓄水之和介于344.7～459.9 mm，樹高生長(zhǎng)指數(shù)平均為0.799，油松生長(zhǎng)受到明顯抑制。

上述研究結(jié)果還說明，通過人工整地起到了良好的蓄水、保水作用。雖說該研究只是對(duì)特定地區(qū)及特定對(duì)象的最低需水量進(jìn)行了估計(jì)，但其結(jié)果為探討集水區(qū)大小提供了線索。根據(jù)公式“需水量-降水量=降水量×集水區(qū)面積”，即可求出集水區(qū)面積或半圓形(如魚鱗坑)集水區(qū)的半徑、矩形集水面的邊長(zhǎng)。在該研究中，油松正常生長(zhǎng)的最低需水量為540 mm、年均降水量為400 mm，則集水區(qū)面積不小于0.35 m2、半圓形集水區(qū)的半徑不小于0.50 m。同理，當(dāng)降水量分別為300、500 mm時(shí)，集水區(qū)面積不小于0.80、0.08 m2，半圓形集水區(qū)的半徑不小于0.75、0.25 m。當(dāng)然，油松正常生長(zhǎng)的最低需水量隨著當(dāng)?shù)厮譅顩r、林分年齡、林分密度等因素的改變而變化[16-18]，在實(shí)際應(yīng)用中必須考慮。根據(jù)長(zhǎng)期的實(shí)踐經(jīng)驗(yàn)，陜西黃土高原和毛烏素沙地油松適宜的造林密度介于625～ 1 111 株/hm2，并以稀植喬木、喬冠草結(jié)合為原則，從而形成與地帶性植被類型接近的人工群落。

[1] 馬欽彥.中國(guó)油松生物量的研究[J].北京林業(yè)大學(xué)學(xué)報(bào)，1989，11(4)：1-10.

[2] 肖瑜.陜西省不同氣候區(qū)域油松人工林生物量和生產(chǎn)力的比較研究[J].植物生態(tài)學(xué)與地植物學(xué)學(xué)報(bào)，1990，14(3)： 237-245.

[3] 吳剛，馮宗煒.中國(guó)油松林群落特征及生物量的研究[J].生態(tài)學(xué)報(bào)，1994，14(4)：415-422.

[4] 羅天祥，李文華，趙士洞.中國(guó)油松林生產(chǎn)力格局與模擬[J].應(yīng)用生態(tài)學(xué)報(bào)，1999，10(3)：257-261.

[5] 李軍，王學(xué)春，邵明安，等.黃土高原3個(gè)不同降水量地點(diǎn)油松林地水分生產(chǎn)力與土壤干燥化效應(yīng)模擬[J].林業(yè)科學(xué)，2010，46(11)：25-35.

[6] 吳欽孝，楊文治.黃土高原植被建設(shè)與持續(xù)發(fā)展[M].北京：科學(xué)出版社，1998.

[7] 程積民，萬(wàn)惠娥.中國(guó)黃土高原植被建設(shè)與水土保持[M].北京：中國(guó)林業(yè)出版社，2002.

[8] 王梅，張文輝.不同坡向人工油松林生長(zhǎng)狀況與林下物種多樣性分析[J].西北植物學(xué)報(bào)，2009，29(8)：1678-1683.

[9] 焦醒，劉廣全.陜西黃土高原油松生長(zhǎng)狀況及其影響因子分析[J].西北植物學(xué)報(bào)，2009，29(5)：1026-1032.

[10] 靳天恩，馬彥紅.影響人工油松林生長(zhǎng)的相關(guān)因子的研究[J].防護(hù)林科技，1999(3)：12-14.

[11] 丁曉綱，李吉躍，哈什格日樂.毛烏素沙地氣候因子對(duì)樟子松、油松生長(zhǎng)的影響[J].河北林果研究，2005，20(4)：309-313.

[12] 唐啟義，馮明光.實(shí)用統(tǒng)計(jì)分析及其DPS數(shù)據(jù)處理系統(tǒng)[M].北京：科學(xué)出版社，2002.

[13] 李根前，趙一慶，唐德瑞.毛烏素沙地中國(guó)沙棘生長(zhǎng)過程與水熱條件的關(guān)系[J].西北林學(xué)院學(xué)報(bào)，1999，14(1)：10-15.

[14] 榆林地區(qū)農(nóng)業(yè)計(jì)劃委員會(huì).陜西省榆林地區(qū)農(nóng)業(yè)區(qū)劃[R].榆林：榆林地區(qū)農(nóng)業(yè)計(jì)劃委員會(huì)，1987.

[15] 任德存.走進(jìn)榆林[M].西安：陜西旅游出版社，2005.

[16] 劉淑明，孫長(zhǎng)忠，孫丙寅.油松和側(cè)柏人工林適宜灌溉定額的研究[J].林業(yè)科學(xué)，2004，40(6)：85-90.

[17] 張曉明，于新曉，武思宏，等.黃土高原主要造林樹種需水定額計(jì)算與分析[J].水土保持研究，2008，15(1)：36-40.

[18] 張建軍，賀維，納磊.黃土區(qū)刺槐和油松水土保持林合理密度的研究[J].中國(guó)水土保持科學(xué)，2007，5(2)：55-59.

(責(zé)任編輯 趙粉俠)

Estimation of Minimum Water Requirement for Normal Growth ofPinustabulaeformisin Loess Plateau

WANG Jun-bo1, QI Xi-lin1, LI Zhan-gang1, GAO Hai-yin2, WANG Yao-jun3, GUO Feng4, TANG Cui-ping4

(1.Shelterbelt Construction Station of Shaanxi Province, Xi'an Shaanxi 710082, China；2. Forestry Station of Dingbian County, Dingbian Shaanxi 718600, China；3. Forest Seed & Sapling Station of Wuqi County, Wuqi Shaanxi 717600, China；4. College of Forestry, Southwest Forestry University, Kunming Yunnan 650224, China)

The 13 a oldPinustabulaeformisplantation was taken as study object to determine the minimum water requirement for its normal growth under the natural conditions in the Loess Plateau by means of height growth simulation, stepwise regression, height growth index, and by integrating the precipitation with water storage capacity of the soil after site preparation. The results showed that the minimum water requirement for normal growth ofP.tabulaeformisin the Loess Plateau was 540 mm. And the area of water collection under different precipitation scale was calculated on the basis of minimum water requirement.

Pinustabulaeformis; growth index; capacity of water storage; water requirement; Loess Plateau

2013-12-06

陜西荒漠化地區(qū)針葉樹造林關(guān)鍵技術(shù)研究集成與示范推廣項(xiàng)目資助。

10.3969/j.issn.2095-1914.2014.02.005

S715.5

A

2095-1914(2014)02-0025-03

第1作者：王俊波(1963—)，男，碩士，高級(jí)工程師。研究方向：防護(hù)林、防沙治沙等工程管理及科技推廣。Email：wanjj808@126.com。