联合GF-6和Sentinel-2红边波段的森林地上生物量反演

光谱反射率能反映地物差异,是森林地上生物量（Aboveground Biomass,AGB）遥感反演的理论基础。红边波段处于近红外与红光波段交界处快速变化的区域,能对植被冠层结构和叶绿素含量的微小变化做出快速反应,对植被生长状况较敏感。研究以GF-6和Sentinel-2多光谱影像作为数据源,结合野外调查AGB数据,构建落叶松和樟子松AGB线性和非线性估测模型,通过比较模型精度选择最优模型进行森林AGB反演和空间分布制图。结果表明：GF-6和Sentinel-2影像红边波段反射率与落叶松、樟子松AGB均呈显著相关（P<0.05）,红边波段对AGB估测较敏感。多变量估测模型整体估测效果优于单变量模型,所有模型中多元线性回归模型取得了最优的决定系数（落叶松R²=0.66,樟子松R²=0.65）和最低的均方根误差（落叶松RMSE=31.45 t/hm²,樟子松RMSE=54.77 t/hm²）。相比单个数据源,联合GF-6和Sentinel-2影像构建的多元线性回归模型估测效果得到了显著提升,模型RMSE对于落叶松和樟子松AGB估测分别最大降低了22.9%和11.2%。增加红边波段进行AGB估测能显著提高模型估测精度,三组数据源分别加入红边波段信息后进行建模,模型RMSE得到了显著降低。GF-6拥有800 km观测幅宽和高效的重访周期,可以快速地提供大尺度时间序列数据,在森林地上生物量反演和动态监测方面有着很大潜力。     

南方红壤典型水土流失区马尾松林地上林木碳储量的遥感监测——以长汀县河田镇为例

森林碳储量动态变化对揭示区域水土流失治理成效具有重要指示意义。以长汀县河田镇为例,2017年随机设置34个马尾松林样本作为建模集,分别与同期Landsat影像的原始波段、植被指数及主成分因子进行回归分析,构建马尾松（（Pinus massoniana））林地上林木碳储量的最佳反演模型,基于伪不变特征原理的线性归一化法实现该模型在2003、2010年影像上的适用性校正转换,实现研究区2003、2010、2017年马尾松林地上林木碳储量的反演及时空分异特征的研究。结果表明：研究区2017年马尾松林地上林木碳储量最佳遥感反演模型是以绿色植被指数（GNDVI）为自变量构建的指数模型：C₂₀₁₇=0.006e^{14.357GNDVI₂₀₁₇},该模型拟合的决定系数为0.57,平均相对精度为82.19%;2003年、2010年马尾松林地上林木碳储量遥感估测模型为：C₂₀₀₃=0.006e^{（16.4086GNDVI₂₀₀₃+1.1428）}、C₂₀₁₀=0.006e^{（15.1677GNDVI₂₀₁₀+1.5821）},两期校正模型的决定系数均在0.85以上;2003、2010及2017年碳储量分别为8.24 t/hm²、11.34t/hm²、16.14 t/hm²,整体呈上升趋势;地上林木碳储量随海拔、坡度的升高而增加,向阳坡地上林木碳储量高于背阴坡;碳储量增长率随海拔、坡度的升高而降低,背阴坡碳储量增长率高于向阳坡。     

普洱季风常绿阔叶林乔木物种多样性高分辨率遥感估测

季风常绿阔叶林是我国南亚热带典型的地带性植被,也是云南省普洱地区重要森林类型。季风常绿阔叶林乔木物种多样性遥感估测对研究区域尺度生物多样性格局及其规律具有重要作用。根据光谱异质性假说和环境异质性假说,首先使用1m空间分辨率的机载高光谱数据和激光雷达数据提取了光谱多样性特征和垂直结构特征。然后利用基于随机森林算法的递归特征消除方法选择对研究区森林乔木物种多样性指数具有较好解释能力的遥感特征,并对Shannon-Winner物种多样性指数进行建模、制图。研究结果表明：（1）基于机载LiDAR数据提取的垂直结构特征和机载高光谱数据提取的光谱多样性特征均对研究区森林乔木物种多样性具有较好的解释能力,随机森林模型估测结果分别为R²=0.48,RMSE=0.46和R²=0.5,RMSE=0.45;两种数据源融合可以进一步提高遥感数据的森林乔木物种多样性估测精度,随机森林估测模型R²和RMSE分别为0.69和0.37。（2）机载激光雷达数据对研究区针阔混交林乔木物种多样性的估测能力优于机载高光谱数据。（3）机器学习方法有助于从高维遥感数据特征中选择适合于森林乔木物种多样性建模的少量特征。该研究在云南普洱开展对季风常绿阔叶林的遥感估测研究,可为森林生物多样性调查提供补充手段,有助于森林生物多样性大尺度、长期动态监测。     

Estimates of subtropical forest biomass based on airborne LiDAR and Landsat 8 OLI data

 快速、定量、精确地估算区域森林生物量一直是森林生态功能评价以及碳储量研究的重要问题。该研究基于机载激光雷达(LiDAR)点云与Landsat 8 OLI多光谱数据, 借助江苏省常熟市虞山地区55块调查样地数据, 首先提取并分析了87个特征变量(53个OLI特征变量, 34个LiDAR特征变量)与森林地上、地下生物量的Pearson’s相关系数以进行变量优选, 然后利用多元逐步回归法建立森林生物量估算模型(OLI生物量估算模型和LiDAR生物量估算模型), 并与基于两种数据建立的综合生物量估算模型的结果进行比较, 讨论预测结果及其精确性。结果表明: 3种模型(OLI模型、LiDAR模型和综合模型)在所有样地无区分分析时, 地上和地下生物量的估算精度均达到0.4以上, 基于不同森林类型(针叶林、阔叶林、混交林)分析时地上和地下生物量的估算精度均有明显提高, 达到0.67及以上。利用分森林类型模型估算生物量, 综合生物量估算模型精度(地上生物量: R²为0.88; 地下生物量: R²为0.92)优于OLI生物量估算模型(地上生物量: R²为0.73; 地下生物量: R²为0.81)和LiDAR生物量估算模型(地上生物量: R²为0.86; 地下生物量: R²为0.83)。     

西藏自治区草地地上生物量遥感反演研究

草地地上生物量（Aboveground Biomass,AGB）是反映草地生态系统功能和质量的关键指标,大尺度地准确估算草地AGB对草地生态系统经营管理至关重要。研究以MODIS影像为数据源,提取反射率、植被指数和植被产品三种不同类型的特征变量,结合野外实测样地草地AGB数据,构建以多元线性逐步回归为代表的参数模型以及随机森林、支持向量机和kNN等非参数模型进行西藏自治区草地AGB估测及空间分布制图。结果表明：（1）多元线性逐步回归、随机森林、支持向量机和kNN模型在加入植被产品特征变量后,RMSE分别降低了15.8%、13.5%、4.1%和17.3%,表明植被产品作为建模变量用于草地AGB估测可有效提高模型精度;（2）三种组合变量构建的草地AGB估测模型中,反射率、植被指数、植被产品组合构建的模型效果最佳,其中kNN模型估测精度最高,R²达到0.60,RMSE和MAE分别为0.43、0.34 t/hm²;（3）草地AGB空间分布呈现出西北地区较低、中部地区较高且分布形态较破碎和东部地区较高的变化特征;（4）利用MODIS植被产品结合kNN模型的预测值与草地实测的AGB空间分布趋势基本一致。综上,MODIS植被产品结合kNN模型可作为大尺度区域草地AGB遥感估测的有效参考。     

“山水林田湖草生命共同体”要素关联性分析——以长白山地区为例

对“山水林田湖草生命共同体”各要素相互关系的深入理解是对其实行保护、系统治理的重要基础。本研究以长白山地区为研究对象,利用空间叠置分析法定量分析该地区“山水林田湖草生命共同体”各要素间的关联性,解构各要素的关系以及相互影响程度,并量化其依存关系。研究结果显示：研究区内各市人口数量、城镇面积与耕地面积呈显著正相关关系（人口数量：R²=0.789,P < 0.05; 城镇面积：R²=0.863,P < 0.05）;耕地的分布与地表水体的分布具有较高的空间相关性（r=0.812）,两者的数量关系也呈显著正相关（R²=0.96,P < 0.01）;区内河流分布、土壤类型与地形、地貌密切相关;区内各市水土流失面积与其森林覆盖率呈显著线性负相关关系（R²=0.824,P < 0.05）。由此可见,6个要素之间存在很明显的依存关系。因此,长白山地区“山水林田湖草生命共同体”的保护与修复应该重视各要素之间的联系,针对其系统性和整体性特征实施相应的保护和治理。     

滨海湿地植被地上生物量遥感反演

盐沼植被生物量是滨海湿地生态系统碳循环研究的重要参数,是湿地生态系统健康评价、资源可持续利用的关键指标,开展盐沼植被地上生物量监测方法研究具有重要意义。目前,遥感技术在湿地生物量监测领域已经得到广泛应用,但反演方法仍以统计模型为主,模型构建需要实测数据支撑,时空拓展性不强。选择江苏盐城丹顶鹤保护区为研究区,基于冠层辐射(PROSAIL)传输模型,通过局部和全局敏感性分析,对模型参数本地化,构建了互花米草地上生物量半经验反演模型,应用于Landsat 8 OLI遥感影像,获得了互花米草地上生物量的时空分布。研究结果表明,利用PROSAIL模型模拟互花米草冠层反射率,叶面积指数(LAI)、叶片干物质含量(Cm)、叶倾角分布参数(LIDF)、等效水厚度(Cw)、叶绿素含量(Cab)、叶片结构参数(N)为高敏感性参数,类胡萝卜素含量(Car)、土壤参数(P_soil)为低敏感性参数;利用不同时刻的遥感影像反演了地上生物量,遥感反演结果与实测数据对比,拟合度R²为0.83,均方根误差(RMSE)为0.43kg/m²,平均相对误差(MRE)为15.7%,精度较高,模型具有较好的时空普适性。研究发展了盐沼植被地上生物量遥感反演方法,解决了以往过于依赖现场实测数据构建反演模型的局限性,该方法可以为研究滨海湿地生态系统碳循环以及准确估算其碳汇潜力提供技术支持。     

基于多端元光谱混合分析方法的大兴安岭火后植被盖度恢复研究

火干扰是北方针叶林结构、功能及动态的主要调节因子之一。研究火后植被恢复对理解火干扰和生态系统的交互作用具有重要意义。火烧迹地通常由植被与基质混合组成,在中低分辨率（ > 10 m）遥感影像中表现为混合像元,因此研究亚像元尺度上植被的恢复是精确量化植被恢复的关键。本研究以2000年大兴安岭呼中自然保护区中8700 hm²火烧迹地为研究区,以两期（2014年6月1日和2010年6月22日）中分辨率Landsat ETM+影像（30 m）为基础数据,比较多端元光谱混合分析（Multiple Endmember Spectral Mixture Analysis,MESMA）和归一化植被指数（Normalized Difference Vegetation Index,NDVI）获得的植被盖度,以高分辨率（2 m）WorldView-2影像（2014年7月1日）为验证数据,对两种方法计算的植被盖度精度进行比较。结果表明,MESMA方法获得的植被盖度（R²=0.691）与传统的NDVI获得的植被盖度（R²=0.700）精度无统计差异,中烈度下获得的植被覆盖精度高于低、高火烧烈度。为验证同一端元能否运用到不同时相的Landsat影像中,本研究将从2014年影像中获取的最佳端元运用到2010年影像中获得植被盖度图,结果表明2014年与2010年得到的RMSE（均方根误差）均值分别为0.0015和0.0065,说明最佳端元可用于不同时相的影像分解。本研究表明MESMA方法可有效监测北方针叶林中火后植被盖度恢复,并可运用于时间序列遥感影像监测植被恢复动态。     

西藏色季拉山急尖长苞冷杉原始林粗木质残体空间格局分析

粗木质残体（Coarse woody debris,CWD）的空间格局反映了森林群落的死亡格局和干扰格局,在一定程度上体现了群落内林木的死亡过程。采用相邻网格法对色季拉山急尖长苞冷杉（Abies georgei var.smithii）原始林1 hm²固定样地内CWD进行调查,从CWD类型、腐烂等级、径级3个方面对CWD空间分布格局进行分析。结果表明：样地内CWD总密度为582株/hm²,倒木占55.33%,是CWD的主要输入形式。CWD密度在腐烂等级上的分布可用多项式拟合（R²=0.9973）,在径级上的分布可用指数衰减模型拟合（R²=0.9746）,且在不同类型、腐烂等级及径级上的分布差异较大。在50 m尺度内,CWD整体表现为小尺度的集群分布和中、大尺度的随机分布。在3种CWD分类中,仅有大枯枝、Ⅰ级腐烂、径级ⅠCWD在小尺度或中尺度表现为较强的集群分布,其余则均以随机分布为主,只是在个别尺度达到或接近集群分布。不同类型CWD间整体关联不显著,只有枯立木与大枯枝在0-21 m尺度内达到显著负关联。CWD空间分布格局是急尖长苞冷杉原始林的重要结构特征,在很大程度上决定着林下植物群落及林型自然更新格局。     

长白落叶松林龄序列上的生物量及碳储量分配规律

由于多年来的过量采伐和重采轻育,伊春东折棱河林场人工长白落叶松林分质量普遍下降,森林生态功能严重衰退。结合对该研究地同一立地类型的人工长白落叶松林(Larix ologensis)林木各组分生物量垂直分配规律的分析,研究了其生物量在年龄序列上的分布及分配规律,为提高其林分生物量及碳储量采取相应的抚育管理措施提供一定的理论基础。结果表明,处于中龄、近熟及成熟林中的林木树干、树皮及活枝生物量所占比例受年龄影响较小,而叶生物量随林龄增大呈现明显递减变化;不同年龄长白落叶松的垂直分布规律基本一致:其树皮及树干生物量随树高增大呈现递减规律,其活枝及叶生物量主要集中分布于树冠中部,而其死枝生物量未呈现明显分布规律;长白落叶松根系生物量随着林分年龄的增大,其粗根、中根及细根所占比例呈现递减规律,而其大根所占比例随年龄的增大基本呈增大趋势。通过统计分析得出,长白落叶松生物量与林分蓄积的最优模型为:W=0.4909M+9.6624(R²=0.8893),进而估算得出:研究区域幼龄长白落叶松林分生物量为1273.72 t/hm²,碳储量为656.98 t/hm²;中龄长白落叶松林分生物量为15480.13 t/hm²,碳储量为7984.65 t/hm²;近熟、成熟龄长白落叶松林分生物量为7684.41 t/hm²,碳储量为3963.62 t/hm²。随林分结构的改善以及中龄、近熟及成熟林分的不断增加,生物量及碳储量会相应增加。     

基于HJ1B和ALOS/PALSAR数据的森林地上生物量遥感估算

森林地上生物量的精确估算能够减小碳储量估算的不确定性。为了探寻一种有效地提高森林生物量估算精度的方法,探讨了基于遥感物理模型和经验统计模型估算山地森林地上生物量的方法。首先,基于Li-Strahler几何光学模型和多元前向模式(MFM)进行模型模拟,结合查找表算法(LUT)从多光谱图像HJ1B估算贺兰山研究区的森林地上生物量。其次,采用统计方法建立了2种回归模型:(1)多光谱图像HJ1B进行混合像元分解(SMA),并与雷达图像ALOS/PALSAR进行图像融合建立生物量回归模型;(2)雷达图像ALOS/PALSAR后向散射系数和实测生物量建立了生物量回归模型。用实测数据对3种算法估算结果进行精度验证。研究结果表明:采用几何光学模型和MFM算法估算的森林地上生物量精度最好(决定系数R2=0.61,均方根误差RMSE=8.33 t/hm2,P0.001),其估算地上生物量与实测值一致性较好,估算生物量精度略优于SMA估算的精度(R2=0.60,RMSE=9.417 t/hm2);ALOS/PALSAR多元回归估算的精度最差(R2=0.39,RMSE=14.89 t/hm2)。由此可见,采用几何光学模型和混合像元分解SMA适合估算森林地上生物量,利用这2种方法进行森林地上生物量遥感监测研究具有一定的应用潜力。     

星载大光斑LiDAR与HJ-1A高光谱数据联合估测区域森林地上生物量

以吉林省汪清林业局经营区为研究区,利用HJ-1A/HSI高光谱数据和ICESat-GLAS波形数据,估测区域森林地上生物量。从平滑后的GLAS波形数据中提取波形长度W和地形坡度参数TS,建立GLAS森林最大树高估测模型;从GLAS波形数据中提取能量参数I(植被回波能量Ev和回波总能量E之比),建立GLAS森林郁闭度估测模型;利用GLAS估测的森林最大树高和森林郁闭度联合建立森林地上生物量模型。由于GLAS呈离散条带状分布,无法实现区域估测,因此研究将GLAS波形数据与HJ-1A/HSI高光谱数据联合,基于支持向量回归机算法实现森林地上生物量区域估测,得到研究区森林地上生物量分布图。研究结果显示,基于W和TS建立的GLAS森林最大树高估测模型的adj.R~2=0.78,RMSE=2.51m,模型验证的adj.R~2=0.85,RMSE=1.67m。地形坡度参数TS能够有效的降低地形坡度的影响;当林下植被高度为2m时,得到的基于参数I建立的GLAS森林郁闭度估测模型效果最好,模型的adj.R~2=0.64,RMSE=0.13,模型验证的adj.R~2=0.65,RMSE=0.12。利用森林最大树高和森林郁闭度建立的森林地上生物量模型的adj.R~2=0.62,RMSE=10.88 t/hm~2,模型验证的adj.R~2=0.60,RMSE=11.52 t/hm~2。基于支持向量回归机算法,利用HJ-1A/HSI和GLAS数据建立的森林地上生物量SVR模型,生成了森林地上生物量分布图,利用野外数据对得到的分布图进行验证,验证结果显示森林地上生物量估测值与实测值存在很强的线性关系(adj.R~2=0.62,RMSE=11.11 t/hm~2),能够满足林业应用的需要。因此联合ICESat-GLAS波形数据与HJ-1A高光谱数据,能够提高区域森林地上生物量的估测精度。     

森林生物量遥感降尺度研究

森林生物量是评价全球碳氧平衡、气候变化的重要指标。目前已有基于星载激光雷达数据的全球森林生物量产品,但空间分辨率较低,不能很好地满足小区域森林调查和动态监测的需要。针对这一现状,以美国马里兰州两个森林分布状况不同的区域为研究区,基于CMS(Carbon Monitoring System)30 m分辨率和GEOCARBON 1 km分辨率森林地上生物量产品以及TM等数据源,通过升尺度模拟低分辨率生物量数据和直接使用低分辨率产品两种方式,分别尝试建立了多光谱地表参数和低分辨率森林地上生物量之间的统计关系,以此作为降尺度模型实现了森林地上生物量空间分辨率从1 km到30 m的转换,并对降尺度结果进行精度评价和误差分析。结果表明:模拟数据降尺度后的30 m分辨率森林地上生物量空间分布和CMS森林地上生物量分布状况大致相同,RMSE=59.2—65.5 Mg/hm~2,相关系数约为0.7;其降尺度结果优于GEOCARBON产品直接降尺度结果RMSE=75.3—79.9 Mg/hm~2;相较于线性模型,非线性模型能更好地呈现森林地上生物量和地表参数间的关系;总体上,降尺度生物量呈现高值区低估,低值区高估的现象。     

大兴安岭5种典型林型森林生物碳储量

森林生态系统是陆地生态系统的重要碳库,森林生态系统的生物碳储量作为森林生态系统碳库的重要组成部分,对全球碳循环与碳平衡产生重要作用。以大兴安岭5种典型林型为研究对象,结合森林资源清查资料,采用地理信息技术(GIS),将5种林型分龄组分别对乔木层、林下的灌木层、草本层和凋落物层各组分的单位面积生物量、含碳率和生物碳储量进行测定和计量估算,并从林分水平上,采用分龄组的方法,计量估算了生物碳储量。结果表明:大兴安岭5种典型林型不同龄组的生物碳储量分别为:兴安落叶松幼龄林、中龄林、近熟林和成熟林的生物碳储量分别为15.20、50.96、95.80t/hm2和109.33t/hm2;白桦幼龄林、中龄林、近熟林和成熟林的生物碳储量分别为15.36、30.67、41.62t/hm2和64.35t/hm2;樟子松幼龄林、中龄林、近熟林和成熟林的生物碳储量分别为29.89、59.92、90.01t/hm2和117.08t/hm2;蒙古栎幼龄林、中龄林、近熟林和成熟林的生物碳储量分别为11.17、11.90、34.94t/hm2和59.49t/hm2;山杨幼龄林、中龄林、近熟林和成熟林的生物碳储量分别为21.81、28.58、42.84t/hm2和64.39t/hm2。研究发现:5种典型林型不同龄组的森林生物碳储量均随着林龄(幼龄林、中龄林、近熟林和成熟林)的增长而增加,但不同林型的碳汇功能存在差异,同一种林型在不同林龄的生物碳储量增幅差异亦较大。尤其是大兴安岭目前林分质量比较差,幼龄林和中龄林所占的比重较大,若能对现有林分加以更好地抚育和管理,该区森林植被仍具有较大的碳汇潜力,碳汇功能将进一步增强,大兴安岭在国家的生态功能区建设中将发挥更重要的碳汇功能,对此提出了森林生态系统碳增汇管理策略与管理路径。研究结果为正确认识森林生物碳储量对区域碳平衡及生态环境的影响具有重要意义,以及在未来营林、造林活动中充分发挥人工林碳汇效应提供参考依据。     

广西马山岩溶次生林群落生物量和碳储量

岩溶植被在岩溶生态系统碳循环和全球碳平衡中具有重要的作用。通过对马山县岩溶次生林年龄序列(幼龄林、中龄林和老龄林)3个演替阶段9个样地(20 m×50 m)的系统取样调查,研究了停止人为干扰后岩溶次生林生物量和碳储量的变化。结果表明:沿幼林、中林和老林群落的顺向演替发展,群落生物量显著增加(P0.05),从幼林群落的48.17 t/hm2、到中林群落113.47 t/hm2,再到老林群落242.59 t/hm2。老林生态系统的碳储量较高,平均为236.69 t/hm2,中林和幼林较低且非常相近,分别为225.17 t/hm2和224.76 t/hm2,各次生林生态系统的碳储量差异不显著(P0.05)。土壤碳储量的大小顺序为幼林(198.44 t/hm2)中林(167.39 t/hm2)老林(113.43 t/hm2)。沿群落正向演替,各次生林生态系统中植物碳储量和土壤碳储量的比例发生明显的变化。幼林的土壤碳储量占生态系统碳储量的88.29%,植物碳储量只占11.71%;中林相应为74.34%和25.66%;而老林为47.92%和52.08%。可见,随着岩溶植被的正向演替,土壤碳转变为植物碳的趋势十分明显,这是岩溶森林不同于酸性土森林的一个显著特征。     

基于森林群落和光谱曲线特征分异的福建省森林生物量遥感反演

基于福建省Landsat8 OLI影像,利用混合像元分解模型筛选出"纯净"的植被像元,提取296个调查样地对应植被像元的红光和近红外波段的中心波长(分别CWR和CWNIR)及其对应的反射率(分别R和NIR),构建以(NIR-R)/(CWNIR-CWR)为特征指数的叶生物量回归模型。然后根据针叶林、阔叶林及针阔混交林叶生物量与干、枝、叶所组成的地上生物量的关系方程,结合福建省植被覆盖分类数据,估测了整个福建省针叶林、阔叶林、混交林的地上生物量,并绘制了福建省地上生物量分布图。结果表明:红光和近红外两个波段反射率和其中心波长所组成的斜率与叶生物量相关性显著,与针叶林、阔叶林、混交林叶生物量的精度分别达到70.55%、68.89%、51.75%,采用这种方法对福建省叶生物量和地上总生物量进行估算,并进行精度验证,其中,针叶林、阔叶林、混交林叶物量的模型误差(RMSE)分别达到29.2467 t/hm~2(R~2=66.64%)、14.0258 t/hm~2(R~2=61.13%)、10.1788 t/hm~2(R~2=55.43%),地上总生物量的模型精度分别达到49.8315 t/hm~2(R~2=54.65%)、45.1820 t/hm~2(R~2=49.01%)、41.5131 t/hm~2(R~2=38.79%),这说明,采用红光波段和近红外波段与其中心波长所组成的斜率估测森林叶生物量,进而估算其地上总生物量的方法是可行的。     

基于地统计学丰林自然保护区森林生物量估测及空间格局分析

摘要： 基于丰林保护区1997年样地调查数据,根据一元生物量估测模型,计算样地生物量,在此基础上,利用ArcGIS地统计插值方法得到整个研究区森林生物量分布,并从林分结构（林型、林龄组）和地形因子（海拔、坡度、坡向）两个方面对保护区森林生物量空间格局进行了分析。结果表明,利用地统计插值得到区域水平的森林生物量是可行的,保护区森林平均生物量水平为171.5t/hm2,总生物量为3.08Tg(1Tg=1012g);不同林分结构（林型、林龄组）有不同的生物量水平;地形因子与生物量有显著的相关性,并得到它们之间的回归方程,为保护区森林生态系统的可持续经营提供了科学依据。     

基于生物量回归方程估算黔中喀斯特常绿落叶阔叶混交林木本植物的根系生物量

常规根系生物量研究方法在我国西南喀斯特森林地区实施困难,根系挖掘法所得研究结果不确定性高,导致目前根系生物量数据匮乏。选择贵州中部喀斯特常绿落叶阔叶混交林为对象,建立常规的根系生物量回归方程,结合群落调查数据,以期研究该森林木本植物的根系生物量特征及其空间分布格局。利用106株乔木、34株灌木和34株藤本标准木根系数据,构建了5种优势乔木(安顺润楠Machilus cavaleriei、化香树Platycarya strobilacea、云贵鹅耳枥Carpinus pubescens、云南鼠刺Itea yunnanensis和窄叶石栎Lithocarpus confinis)、3种优势灌木(刺异叶花椒Zanthoxylum dimorphophyllum、倒卵叶旌节花Stachyurus obovatus和异叶鼠李Rhamnus heterophylla)和2种优势藤本(藤黄檀Dalbergia hancai Benth和小果蔷薇Rosa cymosa)以及乔木通用、灌木通用和藤本通用共13个根系生物量回归方程。利用这些方程计算得到该喀斯特森林木本植物总根系生物量为22.72Mg/hm2。乔木根系生物量(22.57 Mg/hm2)远高于灌木和藤本,占森林总根系生物量的99.30%。5个优势乔木树种的根系生物量(19.67 Mg/hm2)占森林总根系生物量的86.54%。物种根系发达程度是影响根系生物量空间分布格局的重要因素。研究可为喀斯特地区植被地下生物量与碳储量的全面估算提供一个新途径。