侯建花, 周國(guó)模, 王國(guó)英, 莫路鋒

1 浙江農(nóng)林大學(xué)環(huán)境與資源學(xué)院， 浙江省森林生態(tài)系統(tǒng)碳循環(huán)與固碳減排重點(diǎn)實(shí)驗(yàn)室, 臨安 311300 2 浙江省景寧縣林業(yè)局, 景寧 323500 3 浙江農(nóng)林大學(xué)信息工程學(xué)院, 臨安 311300 4 西安交通大學(xué)電子與信息工程學(xué)院, 西安 710049

面向區(qū)域土壤碳通量估算的多向插值空間采樣策略

侯建花1,2, 周國(guó)模1,*, 王國(guó)英3,4, 莫路鋒1,3

1 浙江農(nóng)林大學(xué)環(huán)境與資源學(xué)院， 浙江省森林生態(tài)系統(tǒng)碳循環(huán)與固碳減排重點(diǎn)實(shí)驗(yàn)室, 臨安 311300 2 浙江省景寧縣林業(yè)局, 景寧 323500 3 浙江農(nóng)林大學(xué)信息工程學(xué)院, 臨安 311300 4 西安交通大學(xué)電子與信息工程學(xué)院, 西安 710049

區(qū)域土壤碳通量的準(zhǔn)確測(cè)量對(duì)陸地生態(tài)系統(tǒng)碳循環(huán)過(guò)程分析具有十分重要的作用。由于土壤碳通量空間異質(zhì)性強(qiáng)，采用隨機(jī)抽樣的方法無(wú)法對(duì)區(qū)域土壤碳通量進(jìn)行準(zhǔn)確估算，而大范圍的多點(diǎn)采樣則需要大量的人力及設(shè)備成本。基于一種自制的儀器，提出了一種遞增式采樣的多向插值采樣策略(MDI Multiple Directional Interpolation)：在設(shè)定初始采樣點(diǎn)的基礎(chǔ)上，通過(guò)對(duì)已有采樣點(diǎn)的測(cè)量，通過(guò)徑向插值的方法計(jì)算采樣點(diǎn)連線交點(diǎn)，將不同徑向計(jì)算值差異最大的點(diǎn)作為新增采樣點(diǎn)，以此逐步增加。通過(guò)對(duì)20 幅的50×50 網(wǎng)格區(qū)域仿真，結(jié)果表明(1)MDI布局策略能夠針對(duì)土壤碳通量的變化情況而反饋采樣點(diǎn)的疏密。(2)誤差分析得出采樣點(diǎn)數(shù)量(n=10)較少, MDI布局策略對(duì)碳通量的估算誤差比隨機(jī)布局策略低, 比平均布局策略稍高；隨采樣點(diǎn)增多，3種布局策略誤差均降低；采樣點(diǎn)數(shù)量n=40, MDI布局策略對(duì)碳通量的估算誤差(0.028)比平均布局策略的誤差(0.32)降低了12.5%，比隨機(jī)布局策略的誤差(0.04)降低了30.0%。MDI布局策略根據(jù)土壤碳通量的變化梯度合理分配采樣點(diǎn)，降低區(qū)域土壤碳通量監(jiān)測(cè)誤差。

多向插值; MDI采樣策略; 土壤碳通量估算

土壤碳通量是陸地生態(tài)系統(tǒng)碳循環(huán)中的重要組成部分[1]，它所釋放的CO2量約占大氣CO2總量的10%—20%，與植被冠層光合作用固碳量相當(dāng)[2-3]。準(zhǔn)確的測(cè)量區(qū)域土壤碳通量對(duì)研究陸地生態(tài)系統(tǒng)碳循環(huán)過(guò)程具有十分重要的作用[4]。現(xiàn)階段遙感和模型相結(jié)合是研究區(qū)域碳收支的重要手段，但結(jié)果不確定性較大，需要采用通量實(shí)測(cè)、模型和遙感數(shù)據(jù)融合的方法來(lái)提高模型結(jié)果的可靠性[5-8]。直接測(cè)定土壤碳通量能夠保證區(qū)域土壤碳通量的精確性。對(duì)土壤碳通量的主要測(cè)量方法為開(kāi)路式的通量測(cè)定，對(duì)于單點(diǎn)的測(cè)量目前認(rèn)為最準(zhǔn)確的是利用LI-8100儀器進(jìn)行測(cè)量。但是由于LI-8100等土壤碳通量測(cè)量?jī)x器昂貴，在長(zhǎng)期持續(xù)的測(cè)量過(guò)程需要耗費(fèi)大量人力，同時(shí)監(jiān)測(cè)點(diǎn)的時(shí)間同步性也很難做到一致，因此通過(guò)簡(jiǎn)單增加監(jiān)測(cè)點(diǎn)的方法進(jìn)行區(qū)域土壤碳通量的監(jiān)測(cè)顯得十分困難。無(wú)線傳感網(wǎng)具有成本低、覆蓋范圍廣、多點(diǎn)同步、持續(xù)監(jiān)測(cè)的優(yōu)點(diǎn)，可應(yīng)用于大范圍的區(qū)域監(jiān)測(cè)，作為一種新興手段已應(yīng)用在土壤碳通量的直接測(cè)定中[2]，其相對(duì)誤差已小于5%。

土壤屬性空間分布隨地形、氣候和植被等因素的變化而不同[9]，是一個(gè)復(fù)雜、變化且不均一的連續(xù)體[10]，具有明顯的空間異質(zhì)性[11]。少量監(jiān)測(cè)點(diǎn)的均值不能代表整個(gè)區(qū)域的土壤碳通量，對(duì)區(qū)域土壤通量的估算有很大的誤差。但是監(jiān)測(cè)點(diǎn)數(shù)量過(guò)多產(chǎn)生能源、財(cái)力消耗和數(shù)據(jù)冗余等問(wèn)題。因此如何通過(guò)盡量少的采樣點(diǎn)對(duì)土壤碳通量進(jìn)行測(cè)量以得到較準(zhǔn)確的區(qū)域土壤碳通量值，急需選擇一種合理的采樣策略滿足傳感器布置的要求。目前國(guó)內(nèi)外常用的采樣方法主要有主觀判斷采樣、簡(jiǎn)單隨機(jī)采樣、規(guī)則網(wǎng)格采樣、分區(qū)采樣及混合采樣等[12]。測(cè)量區(qū)域碳通量普遍運(yùn)用規(guī)則網(wǎng)格或隨機(jī)的采樣方式布設(shè)測(cè)量點(diǎn)，然后通過(guò)簡(jiǎn)單的統(tǒng)計(jì)方法計(jì)算區(qū)域碳通量。劉源月等人[1，13-15]在研究森林生態(tài)系統(tǒng)土壤碳通量過(guò)程中，使用簡(jiǎn)單隨機(jī)采樣方法將測(cè)量數(shù)據(jù)的均值作為區(qū)域土壤碳通量。這種方法不能準(zhǔn)確、全面地掌握土壤碳通量的空間分布信息[12]。

本文基于一種自制的儀器Lr100GE-6400，通過(guò)區(qū)域土壤碳通量空間分布仿真數(shù)據(jù)，提出了一種選擇性遞增的采樣方法多向插值布局策略(MDI)，主要研究以下兩個(gè)方面：(1)MDI采樣策略的原理和采樣過(guò)程；(2)MDI采樣策略相對(duì)隨機(jī)采樣、均勻采樣的優(yōu)勢(shì)。通過(guò)盡量少的采樣點(diǎn)對(duì)土壤碳通量進(jìn)行測(cè)量以得到較準(zhǔn)確的區(qū)域土壤碳通量值。

1 基于多向插值的空間采樣策略

在MDI布局策略中，根據(jù)已測(cè)量采樣點(diǎn)的位置及測(cè)量值形成交點(diǎn)，計(jì)算監(jiān)測(cè)區(qū)域中各交點(diǎn)在多個(gè)方向上的插值預(yù)測(cè)結(jié)果，選擇離散度高的交點(diǎn)進(jìn)行逐步增設(shè)采樣點(diǎn)。其中方向指監(jiān)測(cè)點(diǎn)的連線。而交點(diǎn)的離散度是各個(gè)方向上插值的方差。它代表了該點(diǎn)在各個(gè)方向插值的差異，離散度越大差異越大。因此根據(jù)離散度來(lái)確定采樣點(diǎn)的方式在一定程度上考慮了土壤碳通量空間異質(zhì)性。

基于多向插值離散度的空間采樣布局策略MDI步驟為：首先在研究區(qū)域選擇若干初始采樣點(diǎn)；然后依據(jù)空間插值離散度逐步增加采樣點(diǎn)。具體描述如下：

(2) 產(chǎn)生交點(diǎn) 監(jiān)測(cè)區(qū)域內(nèi)已布設(shè)的監(jiān)測(cè)點(diǎn)相互連接形成多條連線，連線相交形成交點(diǎn)。根據(jù)王豐川等[17]研究得到土壤碳通量的最近采樣距離為20 m，在20 m內(nèi)土壤碳通量的變化梯度很小，所以交點(diǎn)距離小于閾值ε(ε=20 m)時(shí)，近似認(rèn)為它們交于同一個(gè)公共交點(diǎn)。如圖2所示，假設(shè)已在5個(gè)位置部署了采樣點(diǎn)，已部署采樣點(diǎn)在監(jiān)測(cè)區(qū)域內(nèi)形成8個(gè)交點(diǎn)。

圖1 初始點(diǎn)的生成示意圖Fig.1 Generation of initial points

圖2 采樣點(diǎn)的生成Fig.2 Generation of sampling points

(3) 計(jì)算多向插值離散度 根據(jù)已測(cè)量的各采樣點(diǎn)的土壤碳通量數(shù)據(jù)，對(duì)各個(gè)方向的交點(diǎn)位置進(jìn)行插值估算，計(jì)算各交點(diǎn)的多向插值離散度。

高東等[16]對(duì)區(qū)域土壤碳通量的研究表明，土壤碳通量在空間上是連續(xù)變化的，因此在各個(gè)方向上可以通過(guò)插值方法估算未知點(diǎn)的碳通量。如圖3單向插值所示， A1、B1為區(qū)域中的兩個(gè)采樣點(diǎn)，已知A1和B1點(diǎn)的土壤碳通量分別為fluxΑ1和fluxB1，可以通過(guò)A1和B1點(diǎn)運(yùn)用插值方法來(lái)估算得到C點(diǎn)的土壤碳通量值fluxc，如式(1)所示，其中X是對(duì)應(yīng)點(diǎn)的相對(duì)橫坐標(biāo)。

(1)

采樣點(diǎn)在監(jiān)測(cè)區(qū)域中相互交叉形成多個(gè)方向，未測(cè)量交點(diǎn)的土壤碳通量根據(jù)多個(gè)方向的插值可以得到多個(gè)估算結(jié)果。如圖3多向插值所示，A1、B1、A2、B2、…、An、Bn為已進(jìn)行測(cè)量的位置，這些位置分別構(gòu)成L1、L2、…、Ln等多個(gè)方向，C點(diǎn)為多個(gè)方向的交點(diǎn)。則未測(cè)量交點(diǎn)C在L1、L2、…、Ln各方向上分別得到不同預(yù)測(cè)值fluxL1、fluxL2、…、fluxLn。計(jì)算出交點(diǎn)上預(yù)測(cè)值的離散度，選擇離散度高的點(diǎn)進(jìn)行采樣，通過(guò)該點(diǎn)的實(shí)際數(shù)據(jù)修正通量值提升各方向估值精度，使整個(gè)區(qū)域土壤碳通量的估算更準(zhǔn)確。多向插值離散度δ計(jì)算方法為(式2):

(2)

(4) 選擇新增采樣點(diǎn) 監(jiān)測(cè)區(qū)域內(nèi)已布設(shè)的監(jiān)測(cè)點(diǎn)形成多個(gè)方向，各個(gè)方向之間相交形成交點(diǎn)，從所有交點(diǎn)中選擇多向插值離散度最高的點(diǎn)作為新增采樣點(diǎn)。

重復(fù)步驟(2)、(3)和(4)，直到各交點(diǎn)的多向插值離散度均低于設(shè)定閾值為止。

圖3 多向插值離散度Fig.3 Dispersion of multiple directional interpolations

2 模擬實(shí)驗(yàn)與分析

通過(guò)模擬區(qū)域土壤碳通量數(shù)據(jù)，利用MDI空間采樣布局策略進(jìn)行20 次仿真實(shí)驗(yàn)，并對(duì)實(shí)驗(yàn)結(jié)果進(jìn)行分析。

2.1 實(shí)驗(yàn)數(shù)據(jù)

根據(jù)張義輝等統(tǒng)計(jì)出的區(qū)域土壤碳通量空間布局結(jié)果[18]，設(shè)定模擬區(qū)域土壤碳通量數(shù)據(jù)的特征(表1)。模擬實(shí)驗(yàn)區(qū)域?yàn)? km×1 km，分割成50 m×50 m的方格矩陣，在每個(gè)方格內(nèi)賦予土壤碳通量值，矩陣行列下標(biāo)(x,y)表示對(duì)應(yīng)模擬區(qū)域的位置，1≤x,y≤50。

表1 實(shí)驗(yàn)數(shù)據(jù)參數(shù)

圖4 土壤碳通量模擬圖Fig.4 Simulation map of soil carbon flux

為了保證實(shí)驗(yàn)?zāi)M數(shù)據(jù)的代表性，實(shí)驗(yàn)按照表1的規(guī)范模擬產(chǎn)生20幅不同的區(qū)域土壤碳通量分布數(shù)據(jù)，作為實(shí)驗(yàn)數(shù)據(jù)。圖4為其中一幅符合空間分布特征(平均值為3.56，標(biāo)準(zhǔn)差為2.54，變程為210， 塊金值為0.81，變異系數(shù)為45)的1 km×1 km區(qū)域土壤碳通量空間分布模擬圖。在圓形范圍1附近，土壤的碳通量值變化梯度相對(duì)較大；在圓形范圍2附近，土壤的碳通量值變化梯度相對(duì)較小。

2.2 實(shí)驗(yàn)方案

在仿真試驗(yàn)中，采樣點(diǎn)位置用模擬數(shù)據(jù)的坐標(biāo)表示，采樣點(diǎn)的測(cè)量值用該位置的土壤碳通量仿真值表示。使用傳統(tǒng)和地統(tǒng)計(jì)學(xué)的方法對(duì)土壤碳通量的空間分布進(jìn)行比較，運(yùn)用地統(tǒng)計(jì)學(xué)Kriging插值方法能更準(zhǔn)確地展現(xiàn)土壤碳通量的空間分布[12]。通過(guò)將采樣點(diǎn)的Kriging空間插值結(jié)果與整個(gè)區(qū)域的仿真值進(jìn)行對(duì)比，評(píng)價(jià)部署策略的性能。實(shí)驗(yàn)對(duì)傳統(tǒng)的主觀判斷采樣方法、分區(qū)和混合采樣和MDI布局策略3采樣方法進(jìn)行了比較。

基于20幅區(qū)域土壤碳通量分布的模擬數(shù)據(jù)，分別運(yùn)用MDI布局、隨機(jī)布局和均勻布局策略，進(jìn)行布設(shè)測(cè)量點(diǎn)。采樣點(diǎn)數(shù)量分別為區(qū)域土壤碳通量總數(shù)的0.4%、0.8%、1.2%、1.6%的，即采樣點(diǎn)數(shù)量n=10, 20, 30, 40(圖5)。在采樣點(diǎn)數(shù)量相同的條件下，比較MDI布局策略與隨機(jī)布局策略、均勻布局策略對(duì)區(qū)域土壤碳通量的估算效果。

圖5 MDI (Multiple Directional Interpolation)測(cè)量點(diǎn)布局Fig.5 MDI distribution of sampling points

2.3 實(shí)驗(yàn)結(jié)果

通過(guò)對(duì)20幅土壤碳通量空間模擬數(shù)據(jù)的布點(diǎn)實(shí)驗(yàn)結(jié)果進(jìn)行統(tǒng)計(jì)，得出MDI布局策略與隨機(jī)布局策略、均勻布局策略的區(qū)域碳通量估算誤差(圖6)以及MDI布局策略與隨機(jī)布局策略、均勻布局策略的估算誤差的標(biāo)準(zhǔn)差(圖7)。通過(guò)分析得出：

(1) 3種布局策略的誤差和誤差標(biāo)準(zhǔn)差與測(cè)量點(diǎn)的數(shù)量的關(guān)系整體趨勢(shì)是相同的：測(cè)量點(diǎn)數(shù)量越多，誤差和誤差標(biāo)準(zhǔn)差減小，趨勢(shì)逐漸變緩。

(2) 在相同的測(cè)量點(diǎn)數(shù)量n(n=10, 20, 30, 40)下，隨機(jī)布局策略的誤差和誤差標(biāo)準(zhǔn)差均比MDI布局策略的大。采樣點(diǎn)數(shù)量(n=40)時(shí)，MDI布局策略對(duì)碳通量的估算誤差(0.04)比隨機(jī)布局策略的估算誤差(0.028)低30.0%。MDI布局策略和隨機(jī)布局策略的誤差差異極顯著(P<0.01)(圖6)，MDI布局策略明顯優(yōu)于隨機(jī)布局策略。

(3) 在測(cè)量點(diǎn)數(shù)量n=10時(shí)，均勻布局策略的誤差和誤差標(biāo)準(zhǔn)差比MDI布局策略的低；隨著測(cè)量點(diǎn)數(shù)量n(n= 20, 30, 40)增大，均勻布局策略的誤差和誤差標(biāo)準(zhǔn)差比MDI布局策略的大。采樣點(diǎn)數(shù)量(n=40)時(shí)，MDI布局策略對(duì)碳通量的估算誤差(0.028)比平均布局策略的誤差(0.32)低12.5%。MDI布局策略和均勻布局策略的誤差差異極顯著(P<0.01)(圖6)，增加測(cè)量點(diǎn)數(shù)量后，MDI布局策略明顯優(yōu)于均勻布局策略。

圖6 碳通量估算誤差Fig.6 Estimation error of carbon flux

圖7 碳通量估算誤差的標(biāo)準(zhǔn)差Fig.7 Standard deviation of carbon flux estimation error

2.4 實(shí)驗(yàn)分析

對(duì)MDI布局策略的估算效果優(yōu)于均勻和隨機(jī)布局的原因進(jìn)行分析，圖8是對(duì)應(yīng)模擬區(qū)域中土壤碳通量采樣結(jié)果。測(cè)量點(diǎn)n=40時(shí)，MDI布局策略的估算誤差與隨機(jī)布局策略、均勻布局策略(圖9)的誤差對(duì)比得出，對(duì)應(yīng)區(qū)域土壤碳通量分布(圖4)范圍1中，均勻布設(shè)策略和隨機(jī)布設(shè)策略誤差比MDI布設(shè)策略的估算誤差高；對(duì)于范圍2， MDI布局策略比其它兩種的估算誤差略高。均勻和隨機(jī)兩種布局策略在范圍1中的較大的估算誤差導(dǎo)致最終估算誤差及誤差標(biāo)準(zhǔn)差較大。

圖8 MDI與隨機(jī)方法、均勻方法的采樣點(diǎn)分布對(duì)比(n=40)Fig.8 Contrast between MDI and random & uniform sampling points distribution

圖9 MDI與隨機(jī)方法、均勻方法的誤差分布對(duì)比Fig.9 Contrast between MDI and random & uniform distribution of error

局部測(cè)量點(diǎn)密度(表2)是導(dǎo)致誤差分布不同的原因。在土壤碳通量變化較復(fù)雜的范圍1處，較多測(cè)量點(diǎn)能夠準(zhǔn)確地得到該區(qū)域土壤碳通量的變化趨勢(shì)，而測(cè)量點(diǎn)不足是導(dǎo)致該區(qū)域土壤碳通量的變化趨勢(shì)不精確的原因。

表2 MDI(Multiple Directional Interpolation)與隨機(jī)方法、均勻方法的采樣點(diǎn)數(shù)量分布對(duì)比

對(duì)應(yīng)圖4中范圍1的區(qū)域，將圖8采樣分布和圖9誤差分布進(jìn)行對(duì)應(yīng)分析，MDI中測(cè)量點(diǎn)比較密集，高達(dá)21個(gè)(表2)；但是隨機(jī)布局策略和均勻布局策略(圖8)相對(duì)MDI的測(cè)量點(diǎn)布局要稀疏許多，測(cè)量點(diǎn)數(shù)量分別為16個(gè)和10個(gè)，尤其是隨機(jī)布設(shè)的測(cè)量點(diǎn)隨意性較大，導(dǎo)致此區(qū)域的測(cè)量點(diǎn)不足。而對(duì)應(yīng)圖5中范圍2處，MDI的布點(diǎn)比較稀疏，僅有2個(gè)(表2)；但是隨機(jī)布局策略和均勻布局策略相對(duì)MDI的布設(shè)要密集，其測(cè)量點(diǎn)數(shù)量均為5個(gè)。可見(jiàn)，MDI布局測(cè)量能夠根據(jù)土壤碳通量的變化梯度合理分配測(cè)量點(diǎn)。

3 結(jié)論

本文針對(duì)區(qū)域土壤碳通量估算的采樣點(diǎn)布局問(wèn)題提出了MDI空間采樣布局策略。使用MDI布局策略進(jìn)行模擬實(shí)驗(yàn)，驗(yàn)證了MDI布局比均勻布局和隨機(jī)布局能夠得到更高的估算精度。對(duì)實(shí)驗(yàn)結(jié)果的誤差進(jìn)行比較，采樣點(diǎn)數(shù)量(n=10)較少, MDI布局策略對(duì)碳通量的估算誤差比隨機(jī)布局策略低, 比平均布局策略稍高；隨測(cè)量點(diǎn)增多，3種布局策略誤差均降低；采樣點(diǎn)數(shù)量(n=40)增多時(shí), MDI布局策略對(duì)碳通量的估算誤差(0.028)比平均布局策略的誤差(0.32)低12.5%，比隨機(jī)布局策略的誤差(0.04)低30.0%。綜合比較3種布局方法， MDI布局策略從整體上對(duì)測(cè)量點(diǎn)數(shù)量進(jìn)行了合理的分配，使得整個(gè)區(qū)域內(nèi)誤差低。MDI部署策略在土壤碳通量變化梯度大的區(qū)域設(shè)置較多采樣點(diǎn)，能夠更好的擬合出此區(qū)域土壤碳通量的變化，很大程度上降低了誤差；在土壤碳通量變化梯度小的區(qū)域，MDI布局策略設(shè)置較少的采樣點(diǎn)，但對(duì)土壤碳通量的估算準(zhǔn)確度影響不大。通過(guò)實(shí)驗(yàn)得到MDI的布局策略使每個(gè)采樣點(diǎn)的作用最大化。

在真實(shí)條件下最佳(多向插值離散度最大)位置并不適合部署土壤測(cè)量設(shè)備(例如水池、石頭等非土壤表面)。此時(shí)，可以考慮采用兩種措施：依據(jù)土壤碳通量的局部相關(guān)性，選擇理論最佳部署位置周邊作為新增采樣點(diǎn)；或者在多向插值離散度相差不大的情況下，選用離散度相對(duì)較大的其它位置作為新增采樣點(diǎn)。但無(wú)論哪種方法都降低了MDI布局策略的性能，并且在理論最佳位置附近仍然不適合部署節(jié)點(diǎn)，最佳位置的多向插值離散度與其余交點(diǎn)的多向插值離散度接近時(shí)，不容易確定新增采樣點(diǎn)位置。下一步工作要進(jìn)行包含大面積非土壤表面的區(qū)域土壤碳通量空間采樣策略研究。

[1] 劉源月, 江洪, 邱忠平, 原煥英, 李雅紅. 亞熱帶典型森林生態(tài)系統(tǒng)土壤呼吸. 西南交通大學(xué)學(xué)報(bào), 2009, 44(4): 590-594.

[2] 徐林, 莫路鋒, 宣子蔚, 徐小軍, 王圣輝, 周國(guó)模. 基于WSN的土壤碳通量測(cè)量系統(tǒng). 林業(yè)科學(xué), 2013, 49(2): 122-126.

[3] Raich J W, Schlesinger W H. The global carbon dioxide flux in soil respiration and its relationship to vegetation and climate. Tellus, 1992, 44(2): 81-99.

[4] Fang C, Moncrieff J B, Gholz H L, Clark K L. Soil CO2efflux and its spatial variation in a florida slash pine plantation. Plant and Soil, 1998, 205(2): 135-146.

[5] 米娜, 于貴瑞, 王盤興, 溫學(xué)發(fā), 孫曉敏, 張雷明, 宋霞, 王樹(shù)森. 基于EALCO模型對(duì)中亞熱帶人工針葉林CO2通量季節(jié)變異的模擬. 植物生態(tài)學(xué)報(bào), 2007, 31(6): 1119-1131.

[6] Zhou L, Wang S Q, Ju W M, Xiong Z, Georg K, Chen J M, Shi H. Assessment of carbon dynamics of forest ecosystems in the Poyang lake basin responding to afforestation and future climate change. Journal of Resources and Ecology, 2013, 4(1): 11-19.

[7] Tate K R, Scott N A, Parshotam A, Brown L, Wilde R H, Giltrap D J, Trustrum N A, Gomez B, Ross D J. A multi-scale analysis of a terrestrial carbon budget: is New Zealand a source or sink of carbon. Agriculture, Ecosystems and Environment, 2000, 82(1/3): 229-246.

[8] Potter C S, Klooster S A, de Carvalho C R, Genovese V B, Torregrosa A, Dungan J, Bobo M, Coughlan J. Modeling seasonal and interannual variability in ecosystem carbon cycling for the Brazilian Amazon region. Journal of Geophysical Research, 2001, 106(D10): 10423-10446.

[9] 魏書(shū)精, 羅碧珍, 孫龍, 魏書(shū)威, 劉芳芳, 胡海清. 森林生態(tài)系統(tǒng)土壤呼吸時(shí)空異質(zhì)性及影響因子研究進(jìn)展. 生態(tài)環(huán)境學(xué)報(bào), 2013, 22(4): 689-704.

[10] Huggett R J. Soil chronosequences, soil development, and soil evolution: a critical review. Catena, 1998, 32(3/4): 155-172.

[11] 周濤, 史培軍, 惠大豐, 駱亦其. 中國(guó)土壤呼吸溫度敏感性空間格局的反演. 中國(guó)科學(xué)C輯: 生命科學(xué), 2009, 39(3): 315-322.

[12] 陳天恩, 陳立平, 王彥集, 郜允兵, 任仲山. 基于地統(tǒng)計(jì)的土壤養(yǎng)分采樣布局優(yōu)化. 農(nóng)業(yè)工程學(xué)報(bào), 2009, 25(增2): 49-55.

[13] 王超, 黃蓉, 楊智杰, 劉強(qiáng), 陳光水, 萬(wàn)曉華. 萬(wàn)木林保護(hù)區(qū)柑橘和錐栗園土壤呼吸的比較. 應(yīng)用生態(tài)學(xué)報(bào), 2012, 23(6): 1469-1475.

[14] 常宗強(qiáng), 馮起, 司建華, 李建林, 蘇永紅. 祁連山不同植被類型土壤碳貯量和碳通量. 生態(tài)學(xué)雜志, 2008, 27(5): 681-688.

[15] 周文嘉, 石兆勇, 王娓. 中國(guó)東部亞熱帶森林土壤呼吸的時(shí)空格局. 植物生態(tài)學(xué)報(bào), 2011, 35(7): 731-740.

[16] 高東, 魯紹偉, 繞良懿, 李佳, 張滕, 靳阿亮. 淮北平原四種土地利用類型非生長(zhǎng)季土壤呼吸速率. 農(nóng)業(yè)工程學(xué)報(bào), 2011, 27(4): 94-99.

[17] 王豐川, 劉加珍, 陳永金. 黃河三角洲濕地土壤呼吸及其環(huán)境因子分析. 人民黃河, 2013, 35(1): 81-84.

[18] 張義輝, 李洪建, 榮燕美, 嚴(yán)俊霞, 李君劍. 太原盆地土壤呼吸的空間異質(zhì)性. 生態(tài)學(xué)報(bào), 2010, 30(23): 6606-6612.

Spatial sampling layout strategy for regional soil carbon flux estimation based on multiple directional interpolation

HOU Jianhua1,2, ZHOU Guomo1,*, WANG Guoying3,4, MO Lufeng1,3

1ZhejiangProvincialKeyLaboratoryofCarbonCyclinginForestEcosystemsandCarbonSequestration,SchoolofEnvironmentalandResourceSciences,ZhejiangAgricultureandForestryUniversity,Lin′an311300,China2JingningForestryBureau,Jingning323500,China3SchoolofInformationEngineering,ZhejiangAgricultureandForestryUniversity,Lin′an311300,China4SchoolofElectronicandInformationEngineering,Xi′anJiaoTongUniversity,Xi′an710049,China

Soil carbon flux is a vital part of the terrestrial ecosystem carbon cycle. Carbon dioxide released from soil is 10% to 20% of the total amount of that in atmosphere, equal to carbon assimilation by photosynthesis in plant canopies. Thus accurate measurement of regional soil carbon flux is significant for research on carbon cycling in terrestrial ecosystems. Remote data combined with models to analyze regional carbon balance is important but unreliable. The data from direct measurement of soil carbon flux turns to be more accurate. During the measuring process, random sampling results in inaccurate estimates of the regional soil carbon flux because of strong spatial heterogeneity. Another method, multi-point sampling in a wide area, requires much more human labor and equipment costs. However, a wireless sensor net (WSN) is a new method that has been applied to the measurement of soil carbon flux. Its advantages are low expense, wide coverage, multi-point synchronization, continuous monitoring, and be applicability to a wide range of areas. By using a self-regulating instrument based on WSN, Lr100GE-6400, to measure regional soil carbon flux, this paper presents a new layout strategy, Multiple Directional Interpolation (MDI). The strategy is as follows: (a) Set up the original well-distributed sampling sites and measure the carbon flux of these sites. (b) Calculate the intersections of different connection lines between sampling points by the radial point interpolation method and regard the points with maximum radial differences as new sampling points. Consequently, points are gradually increased based on this method to achieve accurate estimates of regional soil carbon flux. MDI layout strategy considers the spatial heterogeneity of soil carbon flux. It determines the sampling points based on dispersion, so that MDI can provide more accurate and comprehensive spatial information. Based on simulations in 20 fields of the 50×50 mesh region, layout is set up by MDI, random layout strategy, and uniform layout strategy, respectively. Keeping the same sample points, the three strategies were compared, and results were as follows: (1) MDI layout strategy can reflect the density of sampling points in terms of the variation gradient of soil carbon flux, decreasing use of sampling points and increasing the accuracy of measurement. (2) Error analysis of the experiments verified that the MDI layout strategy is more accurate than the average and random layout strategies. If there are fewer sampling points (n=10), errors in the MDI layout strategy are less than the random and more than the average strategies; if sampling points increase, the three types of layout strategies tend to be more accurate. If there are 40 sampling points, errors in MDI (0.028) are 12.5% less than that of the average (0.32),and 30.0% less than that of the random strategy. In conclusion, MDI layout strategy generally allocates the sampling points more reasonably based on the variation gradient of soil carbon flux. More sampling points are selected in the region with a large variation gradient in soil carbon flux, which matches the variation of soil carbon flux accurately, while fewer are located in the region with a small variation gradient in soil carbon flux but these do not affect the accuracy of soil carbon flux measurement. This method maximizes the effect of every sampling point and minimizes errors.

multiple-directional Interpolation; MDI layout strategy; soil carbon flux estimation

國(guó)家自然科學(xué)基金重大項(xiàng)目(61190114/F0102); 國(guó)家自然科學(xué)基金(61303236); 國(guó)家林業(yè)局948項(xiàng)目(2013-4-71);浙江省科技計(jì)劃項(xiàng)目重大科技專項(xiàng)(2012C13011-1); 景寧畬族自治縣科技計(jì)劃項(xiàng)目(2014A05-5)

2014-01-29;

日期:2014-11-19

10.5846/stxb201401290209

*通訊作者Corresponding author.E-mail: zhougm@zafu.edu.cn

侯建花, 周國(guó)模, 王國(guó)英, 莫路鋒.面向區(qū)域土壤碳通量估算的多向插值空間采樣策略.生態(tài)學(xué)報(bào),2015,35(18):6070-6077.

Hou J H, Zhou G M, Wang G Y, Mo L F.Spatial sampling layout strategy for regional soil carbon flux estimation based on multiple directional interpolation.Acta Ecologica Sinica,2015,35(18):6070-6077.