Biome-BGC模型模拟阔叶红松林碳水通量的参数敏感性检验和不确定性分析

生态过程模型的发展为研究者在长时间序列和区域尺度的研究提供了便利, 但模型模拟的准确性受到模型自身结构、模型参数估计合理性的影响。敏感性分析能够定量或定性筛选出对模型模拟结果影响较大的敏感参数, 是模型参数校准过程中的重要工具, 也是建模和应用的先决条件。该文以阔叶红松林为研究对象, 采用全局敏感性分析方法——傅里叶幅度灵敏度检验扩展法(EFAST)对Biome-BGC模型的生理生态参数进行了敏感性分析, 分别分析了红松(Pinus koraiensis)和阔叶树的净初级生产力(NPP)、蒸散(ET)对参数变化的敏感性。结果表明: (1)模拟红松NPP的不确定性高于阔叶树, 但二者的模拟ET的不确定性均较小。阔叶树的NPP和ET对生理生态参数的敏感性总体上都小于红松。(2)无论是红松、阔叶或其他植被类型, 模拟NPP均表现出对叶片碳氮比、细根碳氮比、比叶面积(SLA)和冠层截留系数的敏感性, 这4个参数的高敏感性主要是由模型自身结构所决定的, 与植被类型和研究地区的关系较小。对模拟ET而言, 细根与叶片碳分配比、新茎与新叶碳分配比和SLA均是影响红松和阔叶树ET的敏感参数, 但红松ET主要受参数与参数间的二阶或多阶交互作用的间接影响, 而阔叶树ET则主要是受到敏感参数直接效应的影响。(3)除了上述影响红松和阔叶树碳水通量的共性参数外, 诸如核酮糖-1,5-二磷酸羧化酶中叶氮含量、叶片与细根周转率、所有叶面积与投影叶面积之比等也是对模拟结果有影响的重要参数, 但是其敏感程度随物种不同和研究区不同而不同, 所以这类参数可以根据具体情况进行参数本地化, 对于其他不敏感参数则可以采用模型缺省值。     

Remote sensing estimation of gross primary productivity and its response to climate change in the upstream of Heihe River Basin

AimsQuantifying the gross primary productivity (GPP) of vegetation is of primary interest in studies of global carbon cycle. This study aims to optimize the MODIS GPP model for specific environments of a fragile waterhead ecosystem, by performing simulations of long-term (from 2001 to 2012) GPP with optimized MOD_17 model, and to analyze the response of GPP to the local climatic variations.Methods The original MODIS GPP products that underestimate GPP were validated against two years (2010-2011) of eddy covariance (EC) data at two sites (i.e. an alpine pasture site and a forest site, respectively) in the upstream of Heihe River Basin. Three comparative experiments were then conducted to analyze the effects of input parameters derived from three sources (i.e. meteorological, biome-specific, and fraction of absorbed photosynthetically active radiation (fPAR) parameters) on the model behavior. After refining the model-driven parameters, long-term GPPs of the study area were estimated using the optimized MOD_17 model, and the Least Absolute Deviation method was applied to analyze the partial correlations between interannual GPPs and climatic variables (temperature, precipitation and vapor pressure deficit (VPD)). Important findings The uncertainties in the original MODIS GPP products are attributable to biome-specific parameters, input data (e.g. meteorological and radiometry data) and vegetation maps. At the pasture site, the light use efficiency had the strongest impact on the GPP simulations. The refined fPAR calculated from the leaf area index (LAI) products of Global Land Surface Satellite (GLASS) greatly improved the GPP estimates, especially at the forest site. The GPPs from the optimized MOD_17 model well matched the EC data (R² = 0.90, root mean squared error (RMSE) = 1.114 g C·m^-2·d^-1 at the alpine pasture site; R² = 0.91, RMSE = 0.649 g C·m^-2·d^-1 at the forest site). The time series of GPPs displayed an up trend at an average rate of 9.58 g C·m^-2·a^-1 from 2001 to 2012. Examination of the partial correlations between interannual GPPs and climatic variables showed that the annual mean temperature and VPD generally had significant positive impacts on GPP, and the annual precipitation had a negative impact on GPP.     

Application of a rice simulation model in high temperature sensitivity study

以调试校正较高精度的ORYZA2000模型参数及高温敏感性模拟验证为目的, 为模型适应性和本地化提供依据, 利用江苏省9个试验点5个水稻(Oryza sativa)品种的田间观测数据及当地逐日气象数据, 采用ORYZA2000最新版本(V2.13)水稻生长模型, 首先挑选出5个试验点3个品种的观测数据进行模型参数适应性调试校正, 确定了水稻发育生长阶段的各项参数, 然后用该参数对独立样本的4个试验点2个水稻品种地上部分各器官生物量、叶面积指数动态变化过程及最终产量进行了动态模拟。通过t检验和质量评价指标对模拟结果进行了显著性检验。利用通过检验的模型及其参数在假设环境温度不同时间段的持续升高条件下, 开展了高温对水稻生物量及产量影响的模拟研究, 模拟结果的影响幅度与实际高温处理结果的影响幅度进行了比较。结果表明: 1)经过调试校正获得较高精度的水稻发育阶段各参数, 较准确地模拟了水稻生物量和叶面积指数的动态累积过程, 模拟值与观测值基本一致, 说明校正后参数的合理性和有效性; 2)调整参数后高温敏感性模拟结果表明, 孕穗期到开花期温度连续3天、5天、7天升高到35 ℃时, 总生物量、穗生物量和总产量与对照(CK)相比分别下降了12%-25%; 不同时间段连续升高到38 ℃时下降18%-31%; 不同时间段升高到41 ℃时, 各生物量与对照相比分别下降了20%-38%。模型模拟值与控制试验室的观测数据的下降幅度基本一致, 表明经过参数校正的ORYZA2000可以应用于水稻对气温升高响应的预测。     

Branching and metabolic exponents in seven woody plants

West、Brown和Enquist提出的植物分形网络模型(简称WBE模型)认为: 植物的分支指数(1/a, 1/b)决定植物的代谢指数, 当分支指数1/a、1/b分别为理论值2.0、3.0时, 代谢速率与个体大小的3/4次幂成正比, 但是恒定的3/4代谢指数并不能全面地反映植物的代谢情况。基于分支指数的协同变化, Price、Enquist和Savage对WBE模型进行扩展, 提出植物分支参数协同变化模型(简称PES模型)。该文借助于PES模型分析了7种木本植物的分支指数和代谢指数。结果表明: 物种间叶面积与叶生物量呈异速生长关系, 基于叶面积得到的分支指数1/a和代谢指数θ在物种间无显著差异, 基于叶生物量得到的分支指数1/a、1/b和代谢指数θ在物种间均存在显著差异, 但基于叶面积和叶生物量分别拟合出的整体分支指数1/a、1/b和代谢指数θ与理论值均无显著差异, 且用叶面积作为代谢速率的替代指标比用叶生物量分析得出的代谢指数与理论值更接近。今后研究应当关注植物叶面积与叶生物量的异速生长关系对植物代谢速率及相关功能特性的影响。     

Forest net primary productivity dynamics and driving forces in Jiangxi Province,China

Aims Subtropical forest ecosystem has great carbon sequestration capacity. Net primary productivity (NPP) plays a critical role in forest carbon cycle and is affected by a number of factors, including climate change, atmospheric composition, forest disturbance intensity and frequency, and forest age, etc. However, the contribution of these factors to the temporal-spatial dynamics of NPP is still not clear. Quantifying the main driving forces on the temporal-spatial dynamics of NPP for subtropical forest ecosystems is a critical foundation for understanding their carbon cycle.
Methods We utilized multi-sources dataset, including observed meteorological data, inversed annual maximum leaf area index (LAI), referenced NPP (simulated by Boreal Ecosystem Productivity Simulator (BEPS) model), forest age and forest types, land cover, digital elevation model (DEM), soil texture, CO₂ concentration and nitrogen deposition. We used the InTEC (integrated terrestrial ecosystem carbon-budget) model to simulate the NPP dynamics for forest ecosystems in Jiangxi Province during the period of 1901-2010. The effects of climate change, forest age, CO₂ concentration and nitrogen (N) deposition on forest NPP from 1970 to 2010 were discussed through designed scenarios.
Important findings (1) Validations by flux measurements and forest inventory data indicated that the InTEC model was able to capture the interannual and spatial variations of forest NPP. (2) The average forest NPP was 47.7 Tg C·a^-1 (± 4.2 Tg C·a^-1) during 1901-2010. The NPP in the 1970s, 1980s, 1990s and 2000s was 50.7, 48.8, 45.4, and 55.2 Tg C·a^-1, respectively. As forest regrows, NPP significantly increased for forests in Jiangxi Province in the 2000s, and exceed that in the 1970s for more than 60% of the forest area. (3) During 1970-2010, under the scenarios of disturbance and non-disturbance, the forest NPP were underestimated by 7.3 Tg C·a^-1 (14.5%) and overestimated by 3.6 Tg C·a^-1 (7.1%) compared to the scenarios of all disturbance and non-disturbance factors, respectively. Compared to the average NPP during 1970-2010, climate change decreased NPP by -2.0 Tg C·a^-1 (-4.7%), N deposition increased NPP by 4.5 Tg C·a^-1 (10.4%), CO₂ concentration change, and the integrated fertilization of CO₂ and N deposition increased NPP by 4.4 Tg C·a^-1 (10.3%) and 9.4 Tg C·a^-1 (21.8%), respectively.     

Spatial variation and controlling factors of temperature sensitivity of soil respiration in forest ecosystems across China

土壤呼吸的温度敏感性(Q₁₀)是陆地碳循环与气候系统间相互作用的关键参数。尽管已有大量关于不同类型森林Q₁₀季节和年际变化规律的研究, 但是对Q₁₀在区域尺度的空间变异特征及其影响因素仍认识不足, 已有结果缺乏一致结论。该研究通过整合已发表论文, 构建了中国森林生态系统年尺度Q₁₀数据集, 共包含399条记录、5种森林类型(落叶阔叶林(DBF)、落叶针叶林(DNF)、常绿阔叶林(EBF)、常绿针叶林(ENF)、混交林(MF))。分析了不同森林类型Q₁₀的空间变异特征及其与地理、气候和土壤因素的关系。结果显示, 1) Q₁₀介于1.09到6.24之间, 平均值(±标准误差)为2.37 (± 0.04), 且在不同森林类型之间无显著差异; 2)当考虑所有森林类型时, Q₁₀随纬度、海拔、土壤有机碳含量(SOC)和土壤全氮含量(TN)的增加而增大, 随经度、年平均气温(MAT)、平均年降水量(MAP)的增加而减小。气候(MAT、MAP)和土壤(SOC、TN)因素间存在相互作用, 共同解释了33%的Q₁₀空间变异, 其中MAT和SOC是Q₁₀空间变异的主要驱动因素; 3)不同类型森林Q₁₀对气候和土壤因素的响应存在差异。在DNF中Q₁₀随MAP的增加而减小, 而其他类型森林中Q₁₀与MAP无显著相关性; 在EBF、DBF、ENF中Q₁₀随TN的增加而增大, 但Q₁₀对TN的敏感性在EBF中最高, 在ENF中最低。这些结果表明, 尽管Q₁₀有一定的集中分布趋势, 但仍有较大范围的空间变异, 在进行碳收支估算时应注意尺度问题。Q₁₀的主要驱动因素和Q₁₀对环境因素的响应随森林类型而变化, 在气候变化情景下, 不同森林类型间Q₁₀可能发生分异。因此, 未来的碳循环-气候模型还应考虑不同类型森林碳循环关键参数对气候变化的响应差异。     

Application and comparison of remote sensing GPP models with multi-site data in China

Aims Estimation of gross primary productivity (GPP) of vegetation at the global and regional scales is important for understanding the carbon cycle of terrestrial ecosystems. Due to the heterogeneous nature of land surface, measurements at the site level cannot be directly up-scaled to the regional scale. Remote sensing has been widely used as a tool for up-saling GPP by integrating the land surface observations with spatial vegetation patterns. Although there have been many models based on light use efficiency and remote sensing data for simulating terrestrial ecosystem GPP, those models depend much on meteorological data; use of different sources of meteorological datasets often results in divergent outputs, leading to uncertainties in the simulation results. In this study, we examines the feasibility of using two GPP models driven by remote sensing data for estimating regional GPP across different vegetation types. Methods Two GPP models were tested in this study, including the Temperature and Greenness Model (TG) and the Vegetation Index Model (VI), based on remote sensing data and flux data from the China flux network (ChinaFLUX) for different vegatation types for the period 2003-2005. The study sites consist of eight ecological stations located in Xilingol (grassland), Changbaishan (mixed broadleaf-conifer forest), Haibei (shrubland), Yucheng (cropland), Damxung (alpine meadow), Qianyanzhou (evergreen needle-leaved forest), Dinghushan (evergreen broad-leaved forest), and Xishuangbanna (evergreen broad-leaved forest), respectively. Important findings All the remote sensing parameters employed by the TG and VI models had good relationships with the observed GPP, with the values of coefficient of determination, R², exceeding 0.67 for majority of the study sites. However, the root mean square errors (RMSEs) varied greatly among the study sites: the RMSE of TG ranged from 0.29 to 6.40 g·m^-2·d^-1, and that of VI ranged from 0.31 to 7.09 g·m^-2·d^-1, respectively. The photosynthetic conversion coefficients m and a can be up-scaled to a regional scale based on their relationships with the annual average nighttime land surface temperature (LST), with 79% variations in m and 58% of variations in a being explainable in the up-scaling. The correlations between the simulated outputs of both TG and VI and the measured values were mostly high, with the values of correlation coefficient, r, ranging from 0.06 in the TG model and 0.13 in the VI model at the Xishuangbanna site, to 0.94 in the TG model and 0.89 in the VI model at the Haibei site. In general, the TG model performed better than the VI model, especially at sites with high elevation and that are mainly limited by temperature. Both models had potential to be applied at a regional scale in China.     

Effects of nitrogen addition on soil respiration in shrublands in Mt. Dongling,Beijing, China

Aims Soil respiration from terrestrial ecosystems is an important component of terrestrial carbon budgets. Compared to forests, natural or semi-natural shrublands are mostly distributed in nutrient-poor sites, and usually considered to be relatively vulnerable to environmental changes. Increased nitrogen (N) input to ecosystems may remarkably influence soil respiration in shrublands. So far the effects of N deposition on shrubland soil respiration are poorly understood. The aim of this study is to investigate the soil respiration of Vitex negundo var. heterophylla and Spiraea salicifolia shrublands and their response to N deposition. Methods We carried out a N enrichment experiment in V. negundo var. heterophylla and S. salicifolia shrublands in Mt. Dongling, Beijing, with four N addition levels (N₀, control, 0; N₁, low N, 20 kg N·hm^-2·a^-1; N₂, medium N, 50 kg N·hm^-2·a^-1 and N₃, high N, 100 kg N·hm^-2·a^-1). Respiration was measured from 2012-2013 within all treatments.Important findings Under natural conditions, annual total and heterotrophic respiration were 5.91 and 4.23, 5.76 and 3.53 t C·hm^-2·a^-1 for the V. negundo var. heterophylla and S. salicifolia shrublands, respectively and both were not affected by short-term N addition. In both shrubland types, soil respiration rate exhibited significant exponential relationships with soil temperature. Temperature sensitivity (Q₁₀) of total soil respiration in V. negundo var. heterophylla and S. salicifolia shrublands ranged from 1.44 to 1.58 and 1.43 to 1.98, and Q₁₀ of heterotrophic soil respiration ranged from 1.38 to 2.11 and 1.49 to 1.88, respectively. Short-term N addition decreased only autotrophic respiration rate during the growing season, but had no significant effects on total and heterotrophic soil respiration in V. negundo var. heterophylla shrubland. In contrast, N addition enhanced the heterotrophic soil respiration rate and did not influence autotrophic and total soil respiration in S. salicifolia shrubland.     

Relationship between photochemical reflectance index with multi-angle hyper-spectrum and light use efficiency in urban green-land ecosystems

Aims Light-use efficiency (LUE) is one of critical parameters in the terrestrial ecosystem production studies. Accurate determination of LUE is very important for LUE models to simulate gross primary productivity (GPP) at regional and global scales. We used eddy covariance technique measurement and tower-based, multi-angular spectro-radiometer observations in autumn 2012 to explore the relationship between bidirectional reflectance distribution function (BRDF) corrected photochemical reflectance index (PRI) and LUE in different phenology and environment conditions in urban green-land ecosystems. Methods Using the eddy covariance technique, we estimated the temporal changes in GPP during the autumn 2012 over Beijing Olympic Forest Park. LUE was calculated as the ratio of GPP to the difference between incoming photosynthetically active radiation (PAR) and PAR reflected from the canopy. Daily PRI values were averaged from the BRDF using semi-empirical kernel driven models. The absolute greenness index (2G_RB) was made by webcam at a constant view zenith and view azimuth angle at solar noon. The logistic function was used to fit the time series of the greenness index. The onset of phonological stages was defined as the point when the curvature reached its maximum value. Important findings Webcamera-observed greenness index (2G_RB) showed a decreasing trend. There was a highly significant relationship between 2G_RB and air temperature (R² = 0.60, p < 0.001). This demonstrates that air temperature is the main driving factor to determine the phenology. PRI estimated from multi-angle hyper-spectrum can estimate LUE in urban green-land ecosystems in vigorous photosynthetic period. The correlation was the strongest (R² = 0.70, p < 0.001) in the peak photosynthetic period. PRI relates better to LUE under high temperature (>15 °C) with high vapour pressure deficit (VPD) (>700 Pa) and high PAR (>300 μmol·m^-2·s^-1). The LUE was up-scaled to landscape/regional scales based on these relationships and phenology. It can also be used for the estimation of GPP of urban green-land with high accuracy.     

Photosynthetic characteristics and photosynthesis-light response curve models of summer maize

Aim The photosynthesis-light response curve is the most commonly used method to explore the relationship between photosynthetically active radiation and the net photosynthetic rate, because it is more effective to reflect the plant photosynthetic characteristics. And it is very meaningful for researchers to choose a suitable summer corn (Zea mays) photosynthetic model and optimal the models of photosynthesis-light response curve by analyzing the differences between simulation and observation in each growth period of some plant. So the object of this paper was to propose some useful suggestions for the choice of summer corn photosynthetic modes and the optimization of the photosynthetic light response curves model in further.Methods In this paper, three typical photosynthetic models were used to fitting the photosynthetic light response curve for upper leaf, leaf at ear of grain and lower leaf of summer corn during bell and milk period. And then the fitting degree of each model was compared to the measured data. Photosynthetic active radiation was divided into three parts, and the fitting residual errors of these three models were analyzed individually.Important findings The photosynthetic characteristic parameters such as maximum net photosynthetic rate (P_nmax), saturated light intensity (I_sat) and dark respiration rate (R_d) decreased constantly with a top-down leaf position and the parameters at milk stage were generally lower than the bell stage. Each growth period and leaf position could fit the curve, but some deviation exists for the P_nmax and I_sat in the rectangular hyperbolic model and the non-rectangular hyperbolic model. The results of curve fitting residual showed that the simulation values from Ye Zi-Piao model were closest to the actual values, and especially for the high photosynthetically active radiation section.     

Shrub encroachment effect on the evapotranspiration and its component—A numerical simulation study of a shrub encroachment grassland in Nei Mongol,China

Aims Shrub encroachment is a common global change phenomenon occurring in arid and semi-arid regions. Due to the difficulty of partitioning evapotranspiration into shrub plants, grass plants and soil in the field, there are few studies focusing on shrub encroachment effect on the evapotranspiration and its component in China. This study aims to illustrate shrub encroachment effect on evapotranspiration by the numerical modeling method. Methods A two-source model was applied and calibrated with the measured evapotranspiration (ET) by the Bowen ratio system to simulate evapotranspiration and its component in a shrub encroachment grassland in Nei Mongol, China. Based on the calibrated model and previous shrub encroachment investigation, we set three scenarios of shrub encroachment characterized by relative shrub coverage of 5%, 15% and 30%, respectively, and quantified their effects caused by shrub encroachment through localized and calibrated two-source model.Important findings The two-source model can well reconstruct the evapotranspiration characteristics of a shrub encroachment grassland. Sensitivity analysis of the model shows that errors for the input variables and parameters have small influence on the result of partitioning evapotranspiration. The result shows that shrub encroachment has relatively small influence on the total amount of ET, but it has clear influence on the proportion of the components of evapotranspiration (E/ET). With shrub coverage increasing from 5% to 15% and then 30%, the evapotranspiration decreased from 182.97 to 180.38 and 176.72 W·m^-2, decreasing amplitude values of 0.34% and 0.44%, respectively. On average, E/ET rises from 52.9% to 53.9% and 55.5%, increasing amplitude values to 2.04% and 3.25%. Data analysis indicates that shrub encroachment results in smaller soil moisture changes, but clear changes of ecosystem structure (decreasing ecosystem leaf area index while increasing vegetation height) which lead to the decrease of transpiration fraction through decreasing canopy conductance. The research highlights that, with the shrub encroachment, more water will be consumed as soil evaporation which is often regarded as invalid part of evapotranspiration and thus resulting in the decrease of water use efficiency.     

Dynamics and regulations of ecosystem light use efficiency in a broad-leaved Korean pine mixed forest,Changbai Mountain

Aims Ecosystem light use efficiency (LUE) reflects the ability of CO₂ uptake and light utilization via photosynthesis, which is a key parameter in ecosystem models to evaluate ecosystem productivity. The objectives of this study were to: (1) compare the differences of LUE derived from different methods; (2) elucidate the seasonal dynamics of LUE and its regulatory factors; and (3) evaluate the maximum LUE (LUE_max) and its variability based on eddy-covariance (EC) flux.Methods Using the flux data from an EC tower during 2003-2005 at a broad-leaved Korean pine (Pinus koraiensis) mixed forest, Changbai Mountain, two types of LUE indicators were generated from: 1) the apparent quantum yield (ε) estimated with rectangular hyperbolic curve, and 2) the ecological light use efficiency (LUE_eco) calculated as the ratio between gross ecosystem productivity (GEP) and photosynthetically-active radiation (Q).Important findings The seasonal variation of ε and LUE_eco appeared a unimodal pattern within a year, with the variations significantly dominated by soil surface temperature and Normalized Difference Vegetation Index (NDVI). A positive correlation between GEP and LUE was found for both ε and LUE_eco, with the effect of Q on LUE relatively weak. The increase in diffusion radiation appeared favorable for enhanced LUE. Generally, there was a significant positive relationship between ε and LUE_eco, while ε was higher than LUE_eco, especially during the mid-season. The annual maximum value of ε and LUE_eco was (0.087 ± 0.003) and (0.040 ± 0.002) μmol CO₂·μmol photon^-1 over the three years, respectively. The interannual variability of LUE_max for ε and LUE_eco was 4.17% and 4.25%, respectively, with a maximum difference of >8%, likely resulted from considerable uncertainty in model simulations. Our results indicated that the inversion and optimization of maximum LUE should be taken seriously in the application of LUE models.     

光合电子流对光响应的机理模型及其应用

光合电子流对光响应的机理可以揭示植物光合电子流与光强、植物捕光色素分子物理特性之间的关系。该文讨论了光合电子流对光响应的机理模型的特性以及捕光色素分子的物理性质, 并利用此模型拟合了山莴苣(Lagedium sibiricum)、一年蓬(Erigeron annuus)和紫菀(Aster tataricus)的光合电子流对光响应的曲线。由此模型不仅可以得到植物的最大光合电子流、饱和光强、初始斜率等参数, 还可以获得捕光色素分子有效光能吸收截面和处于最低激发态的捕光色素分子数对光的响应关系。结果表明: 随光强的增加, 山莴苣的捕光色素分子的有效光能吸收截面下降最快, 紫苑的下降速度最慢; 山莴苣处于最低激发态的捕光色素分子数增长速度最快, 紫苑的增长速度最小。捕光色素分子的有效光能吸收截面随光强增加而下降、处于最低激发态的捕光色素分子数随光强增加而增加的特性将减少其光能的吸收和激子的传递, 因而有利于减少强光对植物产生的光伤害。     

菌根共生网络嵌套性判定的零模型选择

零模型是判定网络嵌套性的重要依据, 菌根共生关系网络经常出现高度非对称性, 该文通过探究矩阵非对称变化对基于不同零模型构建方法的网络嵌套性的影响, 试图为非对称网络零模型的选择提供依据。结果表明: 不同零模型保守性不同, 增加限定条件减少零模型构建过程中的自由空间, 高度限定条件易导致第II类错误。高度非对称网络会增加基于完全随机(r00)零模型的矩阵温度(NT)偏离、降低配对重叠度(NODF)偏离, 标准化指数z-score值显示网络非对称增加后有助于NT和NODF显著性判定。行或列限定对非对称网络嵌套性判定的影响存在差异, 列限定(c0)的网络嵌套性判定对网络非对称性变化的响应规律与r00零模型的响应趋势基本一致, 具有更低的嵌套性偏离和标准差值。行限定(r0, 包括行列限定(backtrack))零模型NT值和NT偏移随矩阵非对称性的变化保持稳定, 较之c0零模型在高度非对称网络中呈现更低的NODF偏离值。选用完全随机和限定零模型相结合的方法, 有助于更加准确判断非对称网络是否具有嵌套结构。高度非对称网络嵌套性判定中对行属性特征比较敏感, 不同非对称性网络间嵌套性水平相比较时选用r0零模型要优于r00和c0零模型。     

Soil C mineralization and temperature sensitivity in alpine grasslands of the Qinghai-Xizang Plateau

青藏高原具有独特的海拔、气候和生态系统类型, 弄清其土壤有机质分解及其温度敏感性对于揭示青藏高原土壤碳储量变化及其碳汇功能具有重要意义。该文利用青藏高原西北部草地的11个封育-自由放牧成对草地, 通过测定不同温度(5、10、15、20和25 ℃)培养下的土壤碳矿化速率, 探讨了土地利用方式对该地区土壤碳矿化及其温度敏感性的影响。实验结果表明: 温度对青藏高原高寒草地的土壤碳矿化具有显著影响, 温度越高土壤碳矿化量越大。从东至西, 土壤碳矿化量逐渐降低。草地土壤碳矿化量与土壤有机碳和土壤全氮含量显著正相关; 即土壤有机碳和土壤全氮含量越高, 土壤碳矿化量就越高。土地利用方式对土壤碳矿化的温度敏感性(Q₁₀)无显著影响, Q₁₀值变化范围为1.4-2.4; 其中, 放牧草地Q₁₀的平均值为1.83, 封育草地Q₁₀的平均值为1.86。此外, Q₁₀与土壤有机碳和土壤全氮含量无显著的相关关系, 也无明显的空间格局。放牧和封育对青藏高原高寒草地土壤碳矿化的温度敏感性无显著影响, 为深入分析青藏高原土壤碳汇功能及其对未来气温升高的响应提供了重要的理论依据。     

Modeling study on photosynthetic-light response curves of a C4 plant,maize

Aims A light response curve can reflect a plant’s ability to utilize light, which is also a key tool in determining the relationship between photosynthetic capacity and environmental factors; however the model accuracies concerning the light response curve remain elusive. The objectives of this study were to compare and assess the model accuracies related to a light response curve and the effects of drought. Methods A field rain shelter was used to control the soil water conditions. To obtain photosynthesis parameters from the light response curve and the drought effects, the relevant models (including the rectangular model, non-rectangular hyperbolic model, modified rectangular hyperbolic model, exponential model, quadratic function model, and a newly modified model) were applied to fit the light response curves. The validity of each model was tested by analyzing the differences between the fitted values obtained by the models and the measured values. Important findings The newly modified model has been proved to performing relatively better in accurately describing the light response curve patterns, and credibly obtaining the crucial photosynthetic parameters such as the maximum net photosynthetic rate, light saturation point, light compensation point, and dark respiration rate, especially under high radiation conditions.     

中国典型陆地生态系统波文比特征及影响因素

波文比(β)是陆面过程中的重要参数, 影响着地表和大气间的能量交换, 明确β的空间变异规律和影响因素有助于对地表能量平衡和气候间反馈关系的预测。该研究收集了在中国不同生态系统类型开展的用涡度相关法(EC)测量地表能量平衡的公开发表文献, 构建了β和气象环境因子数据库, 分析了β在生态系统之间的差异、空间变异特征及影响因素。主要结果: (1)所有生态系统β平均值为0.95 ± 0.64, 变异系数67%, 偏度1.58, 峰度3.07, 整体服从对数正态分布, β平均值最高为灌木生态系统(1.26), 最低为湿地生态系统(0.49)。(2) β在生态系统类型间差异显著: 森林和湿地生态系统β无显著差异, 灌木生态系统β >草地生态系统 β >森林和湿地生态系统 β, 农田生态系统β介于草地生态系统与森林和湿地生态系统之间。(3) β随着纬度的增加而增加, 不随经度和海拔变化。纬度每增加1°,β增加0.038。(4) β随着年降水量(MAP)、年平均气温(MAT)、净辐射(R_n)、当年降水量(PPT)、当年平均气温(T_a)和叶面积指数(LAI)的增加而降低。(5)不同生态系统中β对生物和非生物因素的响应存在显著差异: 草地、森林和灌木生态系统的β对生物和非生物因素变化较为敏感, 而农田和湿地生态系统的β与所有生物和非生物因素均无显著相关关系。(6) MAP和R_n是β变化的直接影响因素, LAI通过影响R_n间接影响β。结果表明了植被类型与气候因素之间具有交互作用, 能量分配最主要的影响因素是降水, 叶面积对能量分配的调节作用并不显著。     

Exploring the influence of soil types underneath the canopy in winter wheat leaf area index remote estimating

Aims Remote sensing is an effective and nondestructive way to retrieve leaf area index (LAI) from plot, regional and global range. Soil background is one of the confounding factors limiting remotely estimating LAI. And soil type contains a large proportion of soil background information, which can influence the optical properties of vegetation canopy and soil. However, our knowledge on the effects stemmed from soil types underneath the canopy on LAI remote estimating have been in shortage. Thus, this study aims to explore the influences of soil types underneath the canopy on winter wheat LAI remote estimating. Methods We analyzed the sensitivity variation of eight spectral indices, named normalized difference vegetation index (NDVI), modified soil-adjusted vegetation index (MSAVI), modified chlorophyll absorption ratio index 2 (MCARI₂), red edge inflection point (REIP), red edge amplitude (Dr), red edge area (SDr), red edge symmetry (RES), normalized difference spectral index (NDSI), to LAI in different soil types, and then we identified some spectral intervals or parameters that were insensitive to soil type variations underneath the canopy. We also compared the accuracy of two commonly used regression models, partial least squares regression (PLSR) and random forest regression (RFR), in estimating LAI for different soil types. We also explored the problems arising from applying the regression model developed in single soil type area to complex soil types area in retrieving LAI. Important findings This paper demonstrates the effects of soil types underneath the canopy on LAI retrieving. 1) The sensitivity of spectral indices to LAI is significantly different due to the soil type variation, but REIP has the least effects from soil type variation among the eight spectral indices. Meanwhile, the band selection algorithm of lambda-by-lambda not only chooses the most sensitive spectral interval for LAI, but also provides a feasible way to construct the spectral index that exhibits strong resistances to the effects of soil types underneath the canopy. 2) The accuracy of LAI estimation by regression models differs under soil type considered or not. So we suggest that in small scale researches, especially in a field scale, the ability of regression models in explaining variables is the priority consideration, while the PLSR is superior to RFR in this respect. Under the premise of unknown priori knowledge of land surfaces, the RFR is more suitable for retrieving LAI than PLSR, but land surface priori knowledge is still necessary. These findings provide the theoretical basis and methods for developing remotely sensing estimating LAI models adapted to various land surfaces. Further analysis is needed in applying the findings in more crop types, cultivars and growth stages.     

Leaf decomposition and nutrient release of dominant species in the forest and lake in the Jiuzhaigou National Nature Reserve,China

AimsLitter decomposition is an important ecological process in nutrient cycling and productivity of ecosystems. Our objective is to quantify the differences of litter decomposition and nutrient release (N and P) under the forest and in an alpine lake among the dominant tree species in the Jiuzhaigou National Nature Reserve.
Methods Fresh leaf litters of Abies ernestii, Pinus tabulaeformis, Betula albo-sinensis, and Salix cupularis were collected and placed in bags under the forest and in an alpine lake for a year.
Important findings The mass remaining ratio (MR) of the leaf litters was well predicted with Olson’s decay model (r > 0.93, p < 0.01). The time for 99% decomposition was the shortest for S. cupularis (6.80 a), followed by B. albo-sinensis (10.34 a), A. ernestii (18.88 a), and P. tabulaeformis (27.21 a). These values were 1.48-, 1.55-, 1.80-, and 1.65-folds of the corresponding values in the lake, respectively. Both MR and nitrogen remaining ratio (NR) had significantly negative correlations with the leaf initial N concentration, but significantly positive correlations with the initial C:N. The nutrient release was significantly different among the four species and between the two sites (i.e., forest and alpine lake). The N release of S. cupularis was consistent between forest and the lake (i.e. directly released in the beginning of decomposition), while other species had an obvious N enrichment process before it released. The release of P among was similar among the four species and between the two sites, with a release—enrichment—release pattern. Overall, the leaf litter decomposition appeared as an intricate process that was affected by the litter chemistry and and the environment. The fast litter decomposition in the lake may have a profound influence on the water quanlity in the Jiuzhaigou National Nature Reserve.     

Simulation on the light-response curves of electron transport rate of Quercus variabilis and Robinia pseudoacacia leaves in the Xiaolangdi area,China

为了比较直角双曲线模型、非直角双曲线模型与叶子飘模型的优缺点, 研究阴生叶和阳生叶电子传递速率的差异, 探讨环境/生物因素对电子传递速率等参数的影响, 该文采用LI-6400XT荧光测定系统对黄河小浪底栓皮栎(Quercus variabilis)、刺槐(Robinia pseudoacacia)叶片电子传递速率-光响应(J-I)曲线进行了测定, 利用直角双曲线模型、非直角双曲线模型和叶子飘模型对J-I曲线进行了拟合。结果表明, 3种模型对叶片J-I曲线拟合的决定系数(R ²)在0.96以上, 叶子飘模型的R ²最高(> 0.99)。直角双曲线模型和非直角双曲线模型无法模拟植物叶片光系统II动力学下调现象, 且不能得出饱和光强(I_sat); 直角双曲线模型对最大电子传递速率(J_max)的模拟明显大于实测值; 叶子飘模型能很好地模拟光系统II动力学下调现象, 得出的J_max和I_sat均最接近实测值。对阴生叶和阳生叶J-I曲线研究发现, 栓皮栎、刺槐阴生叶的J_max分别低于阳生叶25.0%和18.0%, 阳生叶的I_sat分别高于阴生叶26.0%和10.1%。栓皮栎和刺槐J_max与气温显著正相关; 刺槐I_sat与气温、土壤水分含量和净光合速率具有显著的正相关关系; 栓皮栎和刺槐J-I曲线初始斜率α均与净光合速率呈显著负相关关系。