黄土高原年均降水量空间插值及其方法比较

采用不同空间插值方法(普通克里格法OK、反距离权重法IDW和径向基函数法RBF),对黄土高原地区85个气象站点1957—2013年逐日降水数据进行年均降水量空间插值分析,采用平均绝对误差(MAE)、均方根误差(RMSE)、Pearson相关系数(R)以及准确度(AC)对插值结果进行对比,研究不同方法对该地区年均降水量插值结果的影响.结果表明:研究期间,黄土高原地区年均降水量的Moran I值为0.67,各站点年均降水量存在较强的空间自相关性.63个预测站点和22个检验站点的验证结果显示,3种方法的实测值和预测值均显著相关,但反距离权重法和径向基函数法的平均绝对误差(MAEIDW=38.98,MAERBF=34.61)和均方根误差(RMSEIDW=51.49,RMSERBF=43.79)均高于普通克里格法.通过普通克里格法的不同半变异函数模型(环形模型、球体模型、指数模型、高斯模型)对比发现,环形模型的平均绝对误差最小(MAE=32.34),准确度最高(AC=0.976),指数模型平均绝对误差最大(MAE=33.24).综合对比不同插值方法实测值与估算值的验证结果表明,采用半变异函数为环形模型的普通克里格法是进行黄土高原年均降水量插值较好的方法.     

某石油泄漏场地土壤二甲苯空间分布的三维插值

为精确界定污染场地中污染物三维空间分布范围,本文以中国某石油储罐泄露污染场地为研究对象,采用最近邻点法、反距离加权法及克里金法3种三维空间插值模型,对场地中特征污染物二甲苯进行三维空间插值计算及污染分布表征,比较不同插值模型的计算精度,揭示场地中实际受污染土方量。结果表明,反距离加权法预测平均误差为4.94,均方根误差为14.11,其精度明显高于另外2种模型,界定的污染分布范围较符合实际污染情况,通过该方法进行插值后计算受污染土方量为666.2 m~3。研究结果对该类型污染场地污染物空间分布表征和修复治理决策制定提供科学依据。     

样本量对云南松幼苗生物量模型构建及预估精度的影响

为探究不同样本量对生物量模型构建及建模精度的影响，以实际调查的20个家系，共615株云南松幼苗为例，通过编写计算机程序建立进行简单随机抽样，构建不同样本量云南松幼苗各器官及单株生物量异速生长方程。利用决定系数（R²）、估计值标准误（SEE）、均方根误差（RMSE）、总相对误差（RS）及平均误差绝对值（MAB），对模型拟合优度与精度进行比较分析。结果表明：幂函数方程可较好地用于估测云南松幼苗生物量；随着样本量的增加模型精度评估指数MAB呈幂函数形式逐渐减小；当样本数量小于200时MAB较为敏感，模型精度较差，样本量大于200时，其精度随样本量逐渐增加，但变化幅度逐步减小并趋于稳定。因此，根据MAB的变化趋势，样本量达到200时可以构建精度较高且稳定模型。     

中国沿海鼠尾藻遗传多样性分析

鼠尾藻是我国沿海重要经济海藻,本文采用AFLP分子标记技术对我国辽宁、山东和浙江沿海的野生鼠尾藻共120个个体进行了遗传多样性分析.5对引物共得到340个位点,多态性位点比例分别为旅顺盐场群体82.31%、大连长海群体78.97%、山东海阳群体91.03%和浙江洞头群体94.87%,平均杂合度分别为0.2088、0.1862、0.2505和0.2736.Shannon多样性指教分别为0.3273、0.2974、0.3911和0.4236,群体间的遗传距离在0.0301～0.1437之间.分析结果表明,4个群体的遗传多样性水平均出现了一定程度的降低.     

基于地理加权回归克里格的日平均气温插值

气温是大量农业、水文、气候、生态模型的输入变量.在地形复杂的区域,考虑气温与环境变量的线性回归关系和残差的自相关性的方法(如回归克里格法,regression Kriging,RK)是目前气温插值的主要方法.但此类方法多使用基于普通最小二乘的全局回归技术,没有顾及回归关系的空间非平稳性.地理加权回归克里格(geographically weighted regression-Kriging,GWRK)是一种既能顾及回归关系的空间非平稳性、又能考虑残差的自相关性的一种插值方法.本文用RK和GWRK对海南岛2013年12月18日的日平均气温进行插值并进行比较研究.依相关性分析和逐步回归分析的结果,采用RK1(以海拔为辅助变量)、GWRK1(以海拔为辅助变量)、RK2(以纬度、海拔、海陆距离为辅助变量)和GWRK2(以海拔、海陆距离为辅助变量)4种模型进行研究,并用80个验证站评估4种模型的精度.结果表明:GWRK1模型的最大正误差、最大负误差、平均绝对误差、均方根误差均最接近于0.从最大正误差、平均绝对误差、均方根误差3个指标看,考虑更多辅助变量的RK2、GWRK2模型反而不及只考虑海拔的RK1、GWRK1模型,表明RK2、GWRK2模型中辅助变量之间的相关性对插值结果有较大影响.     

基于辅助环境变量的土壤有机碳空间插值——以黄土丘陵区小流域为例

利用普通克里格法(OK)、反距离加权法(IDW)、径向基函数法(RBF)、基于土地利用类型修正的普通克里格法(OK_LU)4种插值方法,对黄土丘陵羊圈沟小流域的土壤有机碳含量进行空间插值。预测结果的准确性通过Pearson相关系数(R),平均绝对误差(MAE),均方根误差(RMSE),准确度(AC)来评价。研究结果表明:(1)在前3种常规空间插值方法中,OK对刻画区域土壤有机碳的空间分布趋势效果最佳,其预测MAE值和RMSE值均为最小,Pearson相关系数(R)和准确度(AC)最大,说明其预测结果的准确性最好、预测的极端误差也最小;其次为RBF;IDW预测的效果最差。(2)OK_LU在空间特征表达方面能够更好地反映复杂地形区的局部变异,其插值结果的精度相比OK有一定程度的提高,其平均绝对误差(MAE)从0.900%降到了0.567%,均方根误差(RMSE)从1.101%降到了0.777%,Pearson相关系数(R)从0.4026提高到0.5589,准确度(AC)从0.9081提高到0.9505。综合比较,在黄土丘陵地区,OK_LU能使插值结果的精度有较大提高,是土壤有机碳空间制图的有效途径。     

鼠尾藻生长与生殖的权衡

2009年7月初至2009年9月中旬对烟台长岛县望福礁附近水域鼠尾藻(Sargassum thunbergii)进行了生态学调查。实验设立固定采样地点通过随机采样以及现场标记测量的方法观测了鼠尾藻的生长与生殖情况。结果表明:(1)鼠尾藻生殖始于7月初,8月初为高峰期,此间87%的藻体产生生殖托;(2)藻体长度对成熟比例和生殖力均有显著的影响,藻体成熟具有大小依赖性;(3)营养组织生物量大的鼠尾藻群体生殖组织生物量也相对较大,生殖期间营养组织生长和生殖组织生长呈异速关系;(4)成熟藻体在有性生殖高峰期过后,生长率下降,而始终未成熟的藻体则持续生长(生长率0)。以上结果表明,在鼠尾藻的生活史中存在着生长与生殖的权衡。     

穗状狐尾藻与不同生长期苦草种间竞争研究

水生植物之间的竞争作用在水体生态系统恢复过程中具有重要作用。本研究采用取代系列实验方法,对穗状狐尾藻(Myriophyllum spicatum L.)与不同生长期(幼苗期和成株期)苦草(Vallisneria natans(Lour.)Hara)的种间竞争进行了研究。结果显示,对于不同生长期的苦草,穗状狐尾藻均具有较强的竞争优势,成株期苦草的竞争能力强于其幼苗期;苦草成株组中的穗状狐尾藻平均株高比苦草幼苗组中的穗状狐尾藻平均株高增加了128.6%(P0.05),但苦草成株组中的穗状狐尾藻植株平均干重比苦草幼苗组中的穗状狐尾藻植株平均干重降低了62.8%(P0.05);苦草幼苗组中的穗状狐尾藻分枝较多、株丛大,对苦草幼苗形成了较大的遮光作用;苦草成株组中的穗状狐尾藻分枝少、植株较长,顶端漂浮于水面生长,表明受苦草的种间竞争压力增大,但仍小于穗状狐尾藻的种内竞争;另外,苦草幼苗组中的苦草根长比苦草成株组中的苦草根长平均长28.6%(P0.05),表明在与穗状狐尾藻竞争过程中,苦草幼苗更趋于竞争地下资源。本研究结果说明穗状狐尾藻竞争能力强于苦草,苦草幼苗受穗状狐尾藻的竞争影响较大。     

不同放牧模式下高原鼠兔洞口的空间分布格局变化

2019年4-11月,在四川省若尔盖县选取全年禁牧、冬季放牧、全年放牧3种放牧模式样地,对高原鼠兔（Ochotona curzoniae）洞口的数量、空间分布格局变化及影响因子进行了研究。计算平均最近邻指数和S²/m比值来分析样方（10 m×10 m）和样地水平下的空间分布格局,通过非参数检验评估空间分布格局在不同放牧模式间和月份间的差异显著性。结果显示：（1）高原鼠兔的洞口总数和有效洞口数在不同放牧模式间均存在显著差异,以冬季放牧模式的洞口数量最多;冬季放牧和全年禁牧模式的洞口数量在月份间有显著变化;（2）无论样方和样地水平,洞口空间分布格局在不同放牧模式间都存在显著差异;在同一放牧模式下,所有洞口和有效洞口的分布型存在差异;（3）植被高度和盖度与高原鼠兔的洞口数量和洞口分布格局相关;（4）洞口数量与平均最近邻指数有显著正相关关系。研究结果说明：不同放牧模式会影响高原鼠兔的空间利用模式及种群动态;高原鼠兔洞口数量及分布型受到植被高度和盖度的影响。     

中高纬度地区TRMM卫星降雨数据的精度评价

卫星降雨产品在降雨空间格局分析中扮演着重要角色,其中TRMM 3B42（3B42）是重要的产品之一,具有适宜的空间和时间尺度,已成为生态、气候、水文模型的重要驱动因子.当前对3B42 V7新版本数据的精度和误差认识仍然有限,尤其缺少中高纬度地区的数据精度评价研究,制约其在生态、气候和水文等领域的应用.本文基于内蒙古自治区1998—2012年3B42 V7产品和53个气象站点数据,评估其在中高纬度地区的精度,揭示了各精度评价指标的空间分布特征及其影响因子.总体精度评价结果表明： 3B42在年、月、日尺度上均高估降雨量,且误差随时间尺度增大而增加,日降雨量的平均误差和平均绝对误差分别为-0.06和0.88 mm;降雨事件侦测能力表现良好,公正预兆评分达到0.23,但对降雨事件总频率存在微弱高估.数据精度受海拔和多年平均降雨量影响,随海拔升高,降雨绝对误差减小,降雨事件探测能力减弱,多年平均降雨量对数据精度的影响则相反.     

Bioclimatic and vegetation mapping of a topographically complex oceanic island applying different interpolation techniques

Different spatial interpolation techniques have been applied to construct objective bioclimatic maps of La Palma, Canary Islands. Interpolation of climatic data on this topographically complex island with strong elevation and climatic gradients represents a challenge. Furthermore, meteorological stations are not evenly distributed over the island, with few stations at high elevations. We carried out spatial interpolations of the compensated thermicity index (Itc) and the annual ombrothermic Index (Io), in order to obtain appropriate bioclimatic maps by using automatic interpolation procedures, and to establish their relation to potential vegetation units for constructing a climatophilous potential natural vegetation map (CPNV). For this purpose, we used five interpolation techniques implemented in a GIS: inverse distance weighting (IDW), ordinary kriging (OK), ordinary cokriging (OCK), multiple linear regression (MLR) and MLR followed by ordinary kriging of the regression residuals. Two topographic variables (elevation and aspect), derived from a high-resolution digital elevation model (DEM), were included in OCK and MLR. The accuracy of the interpolation techniques was examined by the results of the error statistics of test data derived from comparison of the predicted and measured values. Best results for both bioclimatic indices were obtained with the MLR method with interpolation of the residuals showing the highest R ² of the regression between observed and predicted values and lowest values of root mean square errors. MLR with correction of interpolated residuals is an attractive interpolation method for bioclimatic mapping on this oceanic island since it permits one to fully account for easily available geographic information but also takes into account local variation of climatic data.     

基于地理加权回归克里格模型的帽儿山地区森林碳储量空间分布

森林碳储量对于全球气候变化具有重要影响,以往的模型估算未考虑到模型残差的空间相关性和碳储量数据的非平稳性,影响模型的预测精度.本研究基于东北林业大学帽儿山实验林场的ETM+遥感影像数据和193块固定样地,利用地理加权克里格回归(GWRK)建立森林碳储量与遥感和地形因子的回归模型,同时对比最小二乘模型(OLS)、地理加权回归模型(GWR)的预测精度.结果表明: 对于帽儿山地区的森林碳储量估算,GWRK的平均绝对误差(MAE)、均方根误差(RMSE)低于OLS模型和GWR模型,GWRK模型的平均误差(ME)低于GWR模型,与OLS模型相近.GWRK模型的预测精度为83.2%,较OLS模型(73.7%)和GWR模型(77.3%)分别提高6%和10%,拟合精度明显提高,说明GWRK模型是森林碳储量估算的有效方法.利用GWRK模型预测的研究区森林碳储量平均值为70.31 t·hm^-2,在海拔较高的地区,森林碳储量值相对较高,说明海拔对其有较大影响.     

基于带宽优选地理加权回归模型的深圳市植被碳储量反演

植被碳储量估测是自然资源监测的重要内容,遥感技术结合地面样地进行反演可以获得区域范围内植被碳储量的空间连续分布,弥补了传统人工抽样调查估测的不足。然而,现有的参数和非参数遥感估测模型大多忽略了样地数据的变异与空间自相关关系。研究以Landsat 8 OLI影像为数据源提取遥感变量,结合植被碳储量实测调查数据,利用最小信息准则(AICc)、最大空间自相关距离(MSAD)和交叉验证(CV)分别确定最优带宽,组合Gaussian、Bi-square和Exponential核函数构建地理加权回归(GWR)模型估算深圳市植被碳储量,并与多元线性回归(MLR)进行比较,选择最优模型绘制深圳市植被碳储量空间分布图。研究结果表明,GWR模型整体精度优于MLR模型,GWR模型的决定系数(R~2)均高于MLR模型,且均方根误差(RMSE)和平均绝对误差(MAE)显著降低。带宽和核函数的选择对GWR模型估测结果产生了显著影响。以CV确定带宽、Exponential为核函数组合构建的GWR模型效果最佳,其R~2为0.697,RMSE为10.437 Mg C/hm~2,相比其它模型精度上升了13.87%—32....     

Ensemble Learning for Spatial Interpolation of Soil Potassium Content Based on Environmental Information

One important method to obtain the continuous surfaces of soil properties from point samples is spatial interpolation. In this paper, we propose a method that combines ensemble learning with ancillary environmental information for improved interpolation of soil properties (hereafter, EL-SP). First, we calculated the trend value for soil potassium contents at the Qinghai Lake region in China based on measured values. Then, based on soil types, geology types, land use types, and slope data, the remaining residual was simulated with the ensemble learning model. Next, the EL-SP method was applied to interpolate soil potassium contents at the study site. To evaluate the utility of the EL-SP method, we compared its performance with other interpolation methods including universal kriging, inverse distance weighting, ordinary kriging, and ordinary kriging combined geographic information. Results show that EL-SP had a lower mean absolute error and root mean square error than the data produced by the other models tested in this paper. Notably, the EL-SP maps can describe more locally detailed information and more accurate spatial patterns for soil potassium content than the other methods because of the combined use of different types of environmental information; these maps are capable of showing abrupt boundary information for soil potassium content. Furthermore, the EL-SP method not only reduces prediction errors, but it also compliments other environmental information, which makes the spatial interpolation of soil potassium content more reasonable and useful.     

Topographic variations of soil properties and stand structure in aPinus thunbergii plantation

Soil properties and above-and belowground forest structure were studied across various topographies in a 20-year-oldPinus thunbergii Parl. plantation on Mt Tanakami, Japan. The soil properties and stand structure varied greatly with slope position fromridge top to valley floor. Soil thickness, fine soil content and soil moisture content were greater in lower slope positions. The amount of organic carbon in the forest floor was greater in upper slope positions. The organic carbon content in the mineral soil was slightly greater in lower slope positions. These changes in soil properties suggested an upslope decrease in decomposition rate and water and/or nutrient availability. The aboveground structure ofP. thunbergii was more developed at lower slope positions. The mean stem diameter, height and volume ofP. thunbergii increased downslope with decreasing tree density. However, fine root biomass increased greatly upslope. This inverse relationship between tree height and fine root biomass indicated morphological plasticity ofP. thunbergii in exploiting environmental heterogeneity. Variations in soil-plant interactions in the stand along various topographies caused spatial heterogeneity in the accumulation pattern of organic matter in plants and the soil.     

Peroxynitrite-scavenging constituents from the brown alga<Emphasis Type="Italic">Sargassum thunbergii</Emphasis>

Peroxynitrite formationin vivo is implicated in numerous human diseases and there is considerable interest in the use of antioxidants and natural products for their treatment. The three components (1–3) isolated fromSargassum thunbergii as well as the organic solventsoluble fractions and the aqueous layer ofS. thunbergii were evaluated for their potential to scavenge authentic ONOO and ONOO derived from 3-morpholinosydnonimine (SIN-1). The antioxidant activity of the individual fractions was in the order of 85% aquaous (aq.) MeOH>n-BuOH>n-hexane> H₂O. The three known compounds sargahydroquinoic acid (1), sargaquinoic acid (2) and sargachromenol (3) showed peroxynitrite-scavenging activities comparable to those of L-ascorbic acid and penicillamine. These results showed a possible antioxidant activity in major constituents ofS. thunbergii.     

Size-dependence of reproductive allocation of Sargassum thunbergii (Sargassaceae, Phaeophyta) in Bohai Bay, China

Reproductive allocation (the proportion of biomass allocated to reproductive tissue, RA) in Sargassum thunbergii was studied in Yantai, Shandong Peninsula, southern coast of Bohai Bay, China. Annual reproduction initiated in mid-June and peaked in mid-July (90 ± 8% fertile thalli and 75 ± 6% RA mean). Both RA and percentage of fertile thalli exhibited significant temporal variations during reproduction. Sterile thalli were only observed in small length hierarchies at peak reproduction and mean values of RA showed a significantly hierarchical variation, suggesting that the size of thalli played an important role in reproduction and RA were size-dependent. Numerous receptacles were produced along the lateral branches during the reproductive period. A distinct seasonal pattern was observed wherein the presence of lateral branches was followed by the onset of reproduction. RA was positively correlated with the number of lateral branches, as well as the total length of lateral branches. In addition, fertile thalli decayed quickly after peak reproduction. The probability of decay was evidently higher for fertile than for sterile thalli, because all surviving thalli were sterile and short during July–August. So, a trade-off between reproduction and survival may exist at individual levels in S. thunbergii.     

双色蜡蘑对黑松幼苗生长量及其根系形态的影响

为研究双色蜡蘑(Laccaria bicolor)对黑松(Pinus thunbergii)幼苗生长及其根系形态的影响,在营养杯培育条件下,用双色蜡蘑液体菌剂对黑松幼苗进行接种处理,接种第15、30、60、90、120天时取样,比较接种和未接种黑松幼苗的生物量、根系形态及根系分形维数的差异。结果表明:双色蜡蘑在黑松幼苗地上植株、地下根系的生长及其生物量的积累方面都表现出明显的促进作用。接种双色蜡蘑也显著改善了根系总长度、分支数、表面积、体积等参数和根系分形维数,并对地下根系生长的促进作用时间早于地上部分,且效应显著高于地上部分。接种双色蜡蘑第15~30天时对地上部分基本无影响,但对根系促进作用明显,而地上部分在第60天时开始表现出显著的生长效应。研究发现,双色蜡蘑能够成功定殖于黑松根部,促进黑松幼苗生长及其生物量的积累,同时显著促进根系总长度、分支数、表面积和体积增加,并使根系分形维数增大,表现出明显的促生作用,且对根系发育的显著促进作用早于地上部分。     

深圳市前海桂湾公园土壤蒸发及植被蒸腾模拟

蒸散发（ET）是连接水文循环和地表能量平衡的重要组成部分,准确的估算土壤蒸发（LE_s）和植被蒸腾（LE_v）对城市水资源分配和管理具有重要意义。针对城市林地蒸散发双源模型的缺失,本研究基于MOD16模型提出了应用于城市林地区域改进的MOD16双源模型。改进的MOD16双源模型更加准确的描述城市林地区域复杂下垫面的能量分配过程。结合高时空分辨率的Sentinel-2遥感卫星影像,对2020年深圳市前海桂湾公园20个无云日的LE_s和LE_v进行反演研究。用Shuttleworth-Wallace （S-W）模型和双作物系数法（FAO dual-K_c）验证模型性能,同时分析改进的MOD16模型对输入变量的敏感度,结果表明：改进的MOD16模型与S-W模型和FAO dual-K_c之间的均方根误差（RMSE）分别为21.39 W/m²和20.41 W/m²,平均绝对误差（MAE）分别为18.81 W/m²和19.05 W/m²,R²分别为0.801和0.634。改进的MOD16模型可以应用在城市林地LE_s和LE_v估算中。桂湾公园地区总蒸散月均值在85-165 W/m²之间。研究区域的LE_s和LE_v呈现出明显的季节变化,春夏季蒸散量升高,秋冬季蒸散量降低,LE_s的变化范围为0-50 W/m²,LE_v的变化范围为0-120 W/m²。LE_v高值在研究区呈现零散分布,多集中在前海桂湾沿岸,LE_s高值集中在研究区的西北、东北和东南部。敏感性分析结果表明,改进后的MOD16模型对植被覆盖、太阳辐射和湿度更敏感。其中,LE_v模拟结果受净辐射和植被覆盖影响最大,LE_s的模拟结果受湿度和植被覆盖的影响最大。因此,应用改进的MOD16模型时需要优先确保这些参数的准确输入。修正后的MOD16模型极大地提高了高分辨率、小尺度区域城市林地LE_s和LE_v模拟的准确性,为准确获取植被耗水量信息从而科学指导城市林地灌溉、解决城市水资源分配与管理提供有力工具。     

The effect of elevated CO2 on the production and respiration of a Sargassum thunbergii community: A mesocosm study

Approximately one‐third of anthropogenic carbon dioxide is absorbed into the ocean and causes it to become more acidic. The Intergovernmental Panel on Climate Change (IPCC) suggested that the surface ocean pH, by the year 2100, would drop by a further 0.3 and 0.4 pH units under RCP (Representative Concentration Pathway) 6.0 and 8.5 climate scenarios. The macroalgae communities that consisted of Sargassum thunbergii and naturally attached epibionts were exposed to fluctuations of ambient and manipulated pH (0.3–0.4 units below ambient pH). The production and respiration in S. thunbergii communities were calculated from dissolved oxygen time‐series recorded with optical dissolved oxygen sensors. The pH, irradiance, and dissolved oxygen occurred in parallel with diurnal (day/night) patterns. According to net mesocosm production – photosynthetically active radiation (PAR) model, the saturation and compensation PAR, the mean maximum gross mesocosm production (GMP), and daily mesocosm respiration were higher in the CO₂ enrichment, than in the ambient condition, while the mean of photosynthetic coefficient was similar. In conclusion, elevated CO₂ stimulated oxygen production and consumption of S. thunbergii communities in the mesocosm. Furthermore, the sensitivity of the GMP of the S. thunbergii community to irradiance was reduced, and achieved maximum production rate at higher PAR. These positive responses to CO₂ enrichment suggest that S. thunbergii communities may thrive in under high CO₂ conditions.