吉林省长白山区胡桃楸天然次生混交林立地指数模型

依据吉林省长白山胡桃楸天然次生混交林12个典型区域的坡向、坡位分布,结合胡桃楸龄级,布设样地109块,以年龄-树高数据为材料,编制了吉林省东部山区胡桃楸天然次生混交林立地指数表.结果表明: 应用3倍树高标准差法剔除异常数据,最终得到80块379对胡桃楸不同龄级的年龄与树高数据,运用7个数学模型进行比较,确定了单形立地指数导向曲线方程式: H=28.0181(1-e^-0.0341A)^1.0503,决定系数R²为0.85,残差平方和为0.89,基准年龄为50年,指数级距为2 m,得到了16～28 m的7条立地指数曲线.根据不同立地类型,建立了多形立地指数模型,幂函数拟合效果较好,即H=B₀t^B1,其中B₀=10^-9e^0.9213SI,B₁=-0.2139SI+5.6183,最终得到了多形立地指数模型方程式H=10^-9e^0.9213SIt^{-0.2139SI+5.6183}.指数级差法和落点检验法对2种立地指数模型的精度检验表明,多形立地指数的精确度较高,能够准确地反映胡桃楸不同立地类型间树高生长的差异,可用于吉林省长白山区胡桃楸天然次生混交林立地质量评价.     

极端干旱区多枝柽柳叶片气孔导度的环境响应模拟

气孔通过调节植物体水分散失和CO₂吸收在植物适应环境变化和环境胁迫中发挥重要作用。该文在对极端干旱区多枝柽柳(Tamarix ramosissma)叶片气体交换参数观测的基础上, 引入诊断函数f(H)对BWB模型和BBL模型提出的气孔导度(g_s)模型中的空气湿度(h_s或D_s)进行了评价, 并将评价结果引入叶子飘和于强推导出的g_s机理模型。结果表明: (1) BWB和BBL模型对h_s (或D_s)的模拟效果存在很大差异: BWB模型拟合效果较好(R² = 0.5354), BBL模型的结果显著但效果较差(R² = 0.1103)。试验结果显示: 随h_s(或D_s)的增大, g_s呈先增大后减小的趋势, 可用Gauss模型进行拟合, R²分别为0.593和0.258, 说明g_s与h_s的关系要比D_s更密切; (2)叶子飘和于强给出的简化模型(Simple模型)和该文给出的指数模型(Gauss-h模型)均具有较好的模拟效果(R²分别为0.8707和0.8286), η值分别为0.1245和0.0171, 其值均介于0-1之间; (3)模型验证中Gauss-h模型较Simple模型明显低估了观测值, 当观测条件无限趋近于Simple模型的假设时, Simple模型的拟合效果可得到显著提高(R² = 0.9606)。     

不同地形条件下植被盖度信息提取技术研究

为系统地研究特定区域的植被盖度信息提取技术, 在不同的地形条件下, 比较了目前流行的多种高光谱遥感植被盖度提取方法。结果表明: 最优高光谱归一化植被指数(NDVI₁)的建模和验证精度均高于其他两种归一化植被指数(NDVI), 直接采用NDVI建立的回归模型对研究区植被盖度的估测能力低于像元二分模型; 阴坡的最佳模型为基于一阶微分的偏最小二乘回归模型(PLSR模型), 其建模决定系数(R²)为0.810, 均方根误差(RMSE)为6.29, 验证R²为0.773, RMSE为8.85; 阳坡的最佳模型为基于二阶微分的PLSR模型, 其建模R²为0.823, RMSE为6.04, 验证R²为0.801, RMSE为7.35; 平原的最佳模型为全受限的线性光谱混合分解模型(FCLS), 其验证R²为0.852, RMSE为5.86。     

黑龙江省红松人工林枝条分布数量模拟

基于黑龙江省佳木斯市孟家岗林场的12块样地65株人工红松解析木的955个枝解析数据,以Poisson回归模型和负二项回归模型作为备选模型,构建了人工红松二级枝条数量分布模型,并采用AIC、Pseudo-R²、均方根误差(RMSE)和Vuong检验对模型的拟合优度进行比较.结果表明: 每轮一级枝条分布数量集中在3~5个,均值为4个,一级枝条分布数量与人工红松自身的枝条属性相关.一级标准枝上二级枝条分布的离散程度较大,利用全部子回归技术构建二级枝条分布数量模型,最终选择以负二项回归模型为基础的E(Y)=exp(β₀+β₁lnRDINC+β₂RDINC²+β₃HT/DBH+β₄CL+β₅DBH)作为二级枝条分布数量最优预测模型(β为参数;RDINC为相对着枝深度;HT为树高;DBH为胸径;CL为冠长).最优模型的Pseudo-R²为0.79,平均偏差接近于0,平均绝对偏差<7.对于所建立的模型,lnRDINC、CL和DBH的参数为正值,RDINC²和HT/DBH的为负值,随着RDINC增大,在树冠内二级枝条分布数量存在最大值.总的来说,所建立的人工红松二级枝条分布数量模型的预测精度为96.4%,可以很好地预估该研究区域人工红松二级枝条分布数量,为以后枝条的光合作用和生物量的研究提供了理论基础.     

鮸形态性状与体质量的关系

为了研究鮸的形态性状与体质量之间的关系,测量了524尾12月龄左右鮸的体质量(Y, g)和8个形态性状(cm),包括全长(X₁)、体长(X₂)、头长(X₃)、躯干长(X₄)、尾长(X₅)、尾柄长(X₆)、尾柄高(X₇)和体高(X₈).采用相关分析、回归分析、通径分析以及模型拟合等方法对其进行分析.结果表明: 鮸所有性状两两之间存在极显著正相关关系,相关系数在0.834～0.941.经共线性诊断排除了体长这一具有严重共线性的性状后,采用逐步引入-剔除法进行多元回归并进行显著性检验,成功筛选出全长(X₁)、躯干长(X₄)、尾柄高(X₇)和体高(X₈)4个统计检验极显著的形态性状,并建立多元回归方程:Y ＝-444.188 +7.943X₁ +12.861X₄ +38.254X₇ + 42.722X₈.通径分析结果显示,4个形态性状中,体高对体质量的直接作用最大(通径系数=0.351),全长次之(通径系数=0.335).对全长和体高与体质量之间的函数关系进行曲线模型拟合得出,在6个拟合模型中,幂函数模型拟合效果最好,方程为:y=0.013 x^2.891(全长); y=2.028 x^2.751(体高).经适合性检验得知,全长比体高通过幂函数对体质量的预测效果更好.研究结果可为鮸的形态性状在选择育种中的有效利用提供重要参考.     

温带-亚热带栓皮栎种子形态的变异及其与环境因子的关系

在区域尺度上(25.14°-40.25° N, 99.87°-122.07° E), 采集20个栓皮栎(Quercus variabilis)种群的种子样品, 测定种宽、种长, 并计算宽长比形态指标, 探讨了区域尺度上种子形态变异特点及其与环境因子的关系。结果表明, 栓皮栎种子的宽度和长度变化幅度分别为1.21-2.18 cm和1.20-2.96 cm; 宽长比的变化幅度为0.57-1.10。栓皮栎的种宽与种长呈显著的正相关关系。单因素方差分析表明, 种宽、种长和宽长比在种群间差异显著(p < 0.001)。种长与等效纬度呈负线性(R² = 0.18; p = 0.05), 与经度呈凸型的变异关系(R² = 0.43; p = 0.009)。种宽和种长与最热月平均气温呈极显著正相关(R² = 0.35; p = 0.006; R² = 0.30; p = 0.012), 而与最湿季降水量呈显著负相关(R² = 0.28, p = 0.019; R² = 0.24, p = 0.017)。种子宽长比没有明显的变化趋势, 大致趋于恒定(0.88 ± 0.08)。     

开放式空气CO2浓度增高试验中的NO和NO2地气交换观测研究

介绍了农田FACE(free-air CO₂ enrichment)试验中的NO和NO₂地气交换观测方法,即静态暗箱采样—NO和NO₂化学发光分析法,并对观测结果进行了分析讨论.此观测方法简单、易于操作,并可获得可靠的NO和NO₂净交换通量观测结果.在稻麦轮作农田的旱地阶段,无论FACE还是对照处理,NO主要表现为地面净排放,NO₂主要表现为地面净吸收.逐日的NO净排放不依赖于土壤温度,但却与土壤含水量呈线性负相关(R²=0.82,P<0.001).NO₂净吸收具有明显的季节变化特征,逐日的净吸收通量随土壤温度和土壤含水量的变化可分别用抛物线方程拟合(温度:R²=0.74,P<0.001;含水量:R²=0.69,P<0.001).大气CO₂浓度升高200±40μmol·mol^-1使NO净排放减弱19%(t检验P=0.096),NO₂净吸收减弱10%(t检验P=0.26),这主要是植物生长受到促进的缘故.     

内蒙古植被降水利用效率的时空格局及其驱动因素

植被降水利用效率(precipitation-use efficiency, PUE)是评价干旱、半干旱地区植被生产力对降水量时空动态响应特征的重要指标。该研究利用光能利用率CASA (Carnegie-Ames-Stanford Approach)模型估算了2001-2010年内蒙古地区植被净初级生产力(net primary productivity, NPP), 结合降水量的空间插值数据, 分析了近10年内蒙古地区植被PUE的空间分布、主要植被类型的PUE,及其时空格局的驱动因素。结果表明: 2001-2010年内蒙古地区所有植被的平均PUE为0.94 g C·m^-2·mm^-1, 且在105-120° E地带性规律明显,PUE上升速率为每10° 0.55 g C·m^-2·mm^-1。各植被类型间PUE差别较大, 其中灌丛PUE最高, 荒漠PUE最低。在不同的降水量区域, 植被PUE的空间分布与气候因子的关系有较大差别, 0-75 mm降水量区间内, PUE随降水量、气温的升高显著下降(R² = 0.226, p < 0.05); 175-300 mm降水量区间内, 植被 PUE的空间变化与降水量和气温呈极显著相关关系(R² = 0.878, p < 0.001), 且随降水量的增加显著上升( R² = 0.94, p < 0.001), 变化速率约为每100 mm降水0.57 g C·m ^-2·mm^-1; 在降水量大于475 mm的区域, 植被PUE的空间分布与降水量、气温的相关性显著(R² = 0.19, p < 0.05), 且随着气温的上升、降水量的下降而增加, 其中气温的贡献是降水量的8.61倍。在不同的降水量区域, 植被 PUE的年际波动与气候因子的关系也有较大差别, 对于年降水量0-220 mm的地区, PUE的年际波动与降水量呈正相关性、与气温呈负相关性; 在年降水量为220-310 mm的地区, PUE的年际波动主要受降水量的控制, 受气温影响较小; 在年降水量>310 mm的地区,PUE的年际波动与降水量、气温均呈正相关关系, 但在降水量越高的地区, PUE的年际波动与降水量的相关性越弱, 与气温的相关性越强。植被覆盖度与PUE的空间分布极显著相关(R² = 0.73, p < 0.001), 且与 PUE的年际波动也存在线性相关关系(R² = 0.11, p < 0.001); 叶面积指数( LAI)与PUE的年际波动呈线性相关关系(R² = 0.42, p < 0.001), 而当 LAI < 3.15时, PUE的空间分布随LAI增加而呈线性增加。     

黄土高原北部典型灌丛枝条生物量估算模型

于2015年8月末在陕西神木县六道沟小流域采集200个柠条和210个沙柳枝条,测定枝条的基径(D)、长度(H)、含水量(W₀)、鲜质量(W_F)和干质量(W),选用指数函数和异速生长方程建立了4种由枝条形态指标估算枝条生物量的简易模型,并对模型的拟合效果进行验证. 结果表明: 对于柠条和沙柳灌丛,基于D、H二者组合变量(D²H)的异速生长方程是估算枝条生物量的最优模型,该模型经线性转化后可以消除生物量数据的异方差性,且拟合效果最优,决定系数(R²)最大,平均误差(ME)、平均绝对误差(MAE)、总相对误差(TRE)、平均系统误差(MSE)和平均绝对百分误差(MPSE)整体上最小,基本满足生态学研究的精度要求.     

刈割、围封、放牧三种利用方式下草原生态系统的多功能性与植物物种多样性之间的关系

随着全球变化对生物多样性的影响不断加剧, 生物多样性与生态系统功能之间相互关系(BEF)的研究显得极为重要。过去的20多年, BEF的研究大多集中在对物种多样性与单一或少数生态系统功能之间关系的探讨, 但生态系统最为重要的价值是同时维持多种服务和功能的能力, 基于此, 该文首次在国内引入近年来不断完善的生态系统多功能性(multifunctionality)的概念, 并对目前主流的评价方法进行了改进, 从而对内蒙古三种利用方式(刈割、围封、放牧)下的草地群落进行了多功能性评价, 并探讨了多功能性与物种多样性之间的关系。结果显示本研究改进的方法和目前主流方法评价得出的多功能性指数在样方和样地尺度上都有很高的相关性(R² = 0.6956, p < 0.0001; R² = 0.9231, p < 0.0001), 表明该文作者改进后的方法是可靠的。重度放牧的草地群落物种多样性水平最低, 绝大多数土壤功能指标较差, 表现出退化特征; 7年的围封和刈割群落均有较高的物种多样性水平和改善的土壤功能指标; 三者的多功能性指数为刈割(0.2178) >围封(0.0704) >放牧(-0.8031)。植被样方主要沿水肥梯度分布; 多样性指数中, 均匀度指数(Pielou index)和丰富度指数(Margelf index)对多功能性的影响作用最大, 均为样方尺度(R² = 0.1871, p < 0.0001; R² = 0.1601, p < 0.0001)小于样地尺度(R² = 0.5921, p = 0.0093; R² = 0.7499, p = 0.0007), 有尺度依赖性; 多功能性在样方和样地尺度上均与物种均匀度呈线性正相关关系, 而与物种丰富度呈单峰曲线关系。该文研究结果表明, 相对于重度放牧和围封, 刈割更有利于维持该地区生态系统的多功能性; 物种丰富度适中且物种分布均匀的生态系统可能有更好的多功能性。     

华北地区胡桃楸林分布规律及群落构建机制分析

探究植物分布规律和群落构建机制是揭示植物群落空间分布、群落物种多样性的形成、发展及其影响因素的重要途径。该文以华北地区胡桃楸(Juglans mandshurica)林为研究对象, 基于野外84个样方的调查数据, 通过径级分析、典范对应分析等方法研究了胡桃楸林的空间分布规律, 并通过亲缘关系指数计算、植物功能性状等方法研究了胡桃楸林的物种共存机制。结果表明, 华北地区胡桃楸的胸径相对较小(平均5.36 cm), 种群年龄较低; 海拔、坡度、坡位和人为干扰程度是影响华北地区胡桃楸分布的主要因子。研究区域内的胡桃楸大部分生长于山体中下部海拔较低的缓坡, 且人为干扰相对较少的区域, 各分布区域中胡桃楸林的分布规律各异。胡桃楸林物种的构建由生态位机制主导, 其中河北、陕西、天津的胡桃楸林物种共存过程主要受负密度制约的影响, 北京、山西的胡桃楸林的物种共存过程主要受环境选择驱动。     

青藏高原东北部祁连圆柏径向生长对不同类型干旱的响应

为研究祁连圆柏径向生长对不同时期(生长季前2—4月和生长季5—7月)气候因子的响应及面对不同类型(高温、缺水、高温+缺水)干旱事件的弹性(抵抗力和恢复力)变化,利用青藏高原东北部17个采样点的祁连圆柏树轮宽度资料,分析径向生长与不同时期气候因子的相关性,探究高低海拔祁连圆柏面对各类干旱事件的弹性差异。结果表明:祁连圆柏径向生长与干旱指数呈显著正相关,与生长季温度呈负相关(P<0.1)。祁连圆柏面对不同时期干旱事件的弹性存在显著差异,生长季前发生的干旱事件中,低海拔圆柏的抵抗力比高海拔增高2.3%,恢复力降低25.1%;生长季干旱事件中,低海拔圆柏的抵抗力比高海拔降低23.7%,恢复力增高107.1%。祁连圆柏面对缺水型干旱时恢复力更强,均值达到1.68,而面对高温型干旱时祁连圆柏的抵抗力更强,均值达到1.43。未来我国西部高山祁连圆柏,尤其是处于低海拔区的,其径向生长受到全球变暖造成的极端干旱事件的影响会更加显著。     

新疆天山峡谷不同坡向野核桃种群分布格局

在天山峡谷野核桃分布的4条沟谷的阳坡和阴坡设置8个4000 m²样地,应用相邻格子法对野核桃每木调查,在不同尺度(面积)上采用方差/均值比率法和偏离指数(C_x)、负二项参数(K)、丛生指数(I)、平均拥挤度(m^*)、聚块性指数(PAI)、聚集指数(C_a)等聚集强度指数对野核桃种群分布格局进行分析.结果表明: 随尺度增大,野核桃种群由集群分布向随机分布过渡;在25 m²尺度上,种群整体呈集群分布,但不同样地种群聚集强度存在较大差异,西沟阴坡种群聚集程度最大,南沟阴坡种群聚集程度最小;在50～400 m²尺度上,种群整体呈随机分布.在25 m²尺度上,幼树、小树和壮年树均呈集群分布,壮年树聚集强度最大,小树聚集强度最小;当取样尺度达到200～400 m²时,幼树、小树、壮年树和大树均呈随机分布.野核桃种群格局规模为25 m²,格局强度为2.49～9.38 m²;种群组分中小树格局强度最大,壮年树格局强度最小.取样尺度与种群斑块大小接近,采用适合的测定方法并应用多指数综合判定,一般可获得可靠的结论.     

亚热带地区常绿阔叶林SPAC系统水分的氢氧稳定同位素特征

森林的土壤-植物-大气连续体(SPAC)是陆地重要的水循环连续界面过程。本研究通过分析亚热带常绿阔叶林的降水、大气水汽、土壤水、叶片水的同位素组成,探讨森林SPAC系统水分的氢氧同位素组成特征以及植物蒸腾与叶片性状和环境因子的关系。结果表明: 研究区大气降水、土壤水、竹柏枝条水、竹柏叶片水和大气水汽的δD-δ¹⁸O线性回归方程分别为: δD_P=7.97δ¹⁸O_P+12.68(R²=0.97)、δD_S=4.29δ¹⁸O_S-18.62(R²=0.81)、δD_B=3.31δ¹⁸O_B-29.73(R²=0.49)、δD_L=1.49δ¹⁸O_L-10.09(R²=0.81)、δD_V=3.89δ¹⁸O_V-51.29(R²=0.46)。在降水→土壤水→植物水的界面水输送过程中,氢氧同位素逐渐富集,而从土壤蒸发和从植物蒸腾的水汽同位素贫化。在降水和蒸发作用的影响下,土壤水同位素随深度增加有贫化的趋势,而且整体上旱季土壤水同位素比雨季富集。观测期间,枝条水同位素比土壤水略微富集,说明水分在植物体内运输过程中存在受到蒸腾富集作用的可能性。旱季,乔木的枝条水同位素比灌木贫化,说明根系分布更深的乔木植物更倾向于利用深层土壤水。由于在叶片性状、蒸腾速率以及对环境因子的响应程度等方面存在差异,不同植物的叶片水同位素组成随叶龄增长的变化特征有所不同。雨季的环境条件更有利于叶片蒸腾,使雨季的叶片水同位素比旱季富集。叶片水同位素组成与植物叶片含水量呈正相关关系,与相对湿度呈负相关关系,综合反映了植物应对环境变化的水分调控功能。     

不同去趋势方法对基于Dendrometer数据的茎干水分动态分析的影响——以白扦为例

树木茎干半径变化记录仪(Dendrometer)监测的高精度数据不仅包括木质部的年内径向生长过程, 还包含由茎干水分的消耗和补充引起的可逆变化。然而, 不同的年内生长去趋势方法获得的茎干水分波动之间的差异性仍缺乏对比研究。基于芦芽山北坡针叶林下限白扦(Picea meyeri) 2015年生长季的茎干半径变化和环境因子的实时监测数据, 使用Gompertz生长模型(GPZ)、线性生长模型(LG)、零生长模型(ZG)、日值法(D)和茎干循环法(SC)模拟并去除茎干年内的生长趋势, 然后提取5种不同类型树木水分缺乏引起的茎干收缩(TWD_GPZ、TWD_LG、TWD_ZG、TWD_D和TWD_SC)以表征茎干水分亏缺, 并进一步对比分析了不同茎干水分亏缺序列对环境中水分状况的响应特征。研究发现: (1)不同去趋势方法计算的茎干水分亏缺的趋势和幅度有所差异, 可聚类为3组: TWD_LG和TWD_ZG、TWD_GPZ以及TWD_D和TWD_SC。同组或聚类距离接近的序列在生长季内每个月份都展现出显著的相关性。然而, TWD_LG、TWD_ZG和TWD_GPZ与TWD_D和TWD_SC在8月份相关性较弱。(2) TWD_D和TWD_SC与空气饱和水汽压差(VPD)的正相关关系比TWD_GPZ、TWD_LG和TWD_ZG更加稳定, 且具有更大的相关系数。5种茎干水分亏缺序列和土壤含水量(SWC)的关系在生长季内变化很大。(3)不同去趋势方法的茎干水分亏缺都随着水分胁迫程度(VPD/SWC)升高而显著增长。当胁迫程度较低时, TWD_SC对VPD/SWC的变化最为敏感(R² = 0.39, p < 0.001), 但是与TWD_ZG差别不大(R ² = 0.37, p < 0.001); 当胁迫程度较高时, TWD_ZG对VPD/SWC的敏感性最高(R ² = 0.59, p < 0.001)。综合对比来看, 零生长模型是比较适合研究区白扦生长季内茎干水分波动的去趋势方法, 其可为干旱胁迫条件下预测研究区树木的茎干水分动态及特征提供科学 依据。     

Algorithm for application of Fourier analysis for biorhythmic baselines of pharmacodynamic indirect response models

The change of an indirect pharmacological response R(t) can be described by a periodic time-dependent production rate k_in (t) and a first-order loss constant k_out. If k_in(t) follows some biological rhythm (e.g., circadian), then the response R(t) also displays a periodic behavior. A new approach for describing the input function in indirect response models with biorhythmic baselines of physiologic substances is introduced. The present approach uses the baseline (placebo) response R_b(t) to recover the equation for k_in(t). Fourier analysis provides an approximate equation for R_b(t) that consists of terms (usually two or three) of the Fourier series (harmonics) that contribute most to the overall sum. The model differential equation is solved backward for k_in(t), yielding the equation involving R_b(t). A computer program was developed to perform the square L²-norm approximation technique. Fourier analysis was also performed based on nonlinear regression. Cortisol suppression after inhalation of fluticasone propionate (FP) was modeled based on the inhibition of the secretion rate k_in(t) using ADAPT II. The pharmacodynamic parameters k_out and IC₅₀ were estimated from the model equation with k_in(t) derived by the new approach. The proposed method of describing the input function needs no assumption about the behavior of k_in(t), is as efficient as methods used previously, and is more flexible in describing the baseline data than the nonlinear regression method. (Chronobiology International, 17(1), 77-93, 2000)     

帽儿山阔叶混交林天然更新幼苗幼树地径-树高模型

基于黑龙江省东北林业大学帽儿山实验林场48块天然阔叶林幼苗幼树调查数据,在8个备选模型中为主要更新树种选择最佳地径(D₀)-树高(H)模型作为基本模型,通过再参数化在基础模型中引入林分因子,并建立样地水平混合效应模型,最后分别对基础模型和混合效应模型进行独立样本检验.结果表明:各树种幼苗幼树的地径-树高关系存在明显的正相关,幂函数或包含幂函数的模型能较好地拟合幼苗幼树地径和树高的关系;基础模型中引入林分因子[林分优势高(H_T)、林分平均胸径(D_g)、林分胸高断面积(BA)]能提高模型的拟合效果,各树种剩余均方根误差(RMSE)下降1.3%～7.4%(平均3.8%),但调整后的决定系数(R_a²)仅仅提高0.1%～1.1%(平均0.6%),赤池信息准则(AIC)下降3.2%～35.2%(平均下降11.4%).对春榆、椴树、水曲柳等10个树种建立混合效应模型,混合效应模型的R_a²比基础模型有所提高,增幅为0.5%～3.5%(平均增加2.2%);RMSE和AIC比基础模型的小,RMSE下降的幅度很大,为3.9%～20.3%,平均下降13.9%,AIC减少4.0%～44.4%(平均减少22.3%).模型检验结果显示,相对于基础模型,混合效应模型的平均绝对误差(MAE)减小0.0001～0.46 m,平均减小0.08 m;平均预测误差百分比(MPSE)降幅较大,为0.1%～6.2%,平均降幅2.0%.说明混合效应模型既能提高模型的拟合效果,又能提高模型的预测能力.本研究构建的阔叶混交林主要更新树种幼苗幼树地径-树高模型为天然阔叶林结构分析和林分生长预测提供了参考.     

Comparative molecular field analysis and QSAR on substrates binding to cytochrome p450 2D6

In this study, we utilized comparative molecular field analysis (CoMFA) to gain a better understanding of the steric and electrostatic features of the cytochrome P450 2D6 (CYP2D6) active site. The training set consists of 24 substrates with reported K_M values from liver microsomal CYP2D6 spanning an activity range of almost three log units. The low energy conformers were fit by root mean square (RMS) to minaprine at the site of metabolism and to the protonated nitrogen. In this manner, we constructed two CoMFA models, one model with a distance constraint and another without. The model with the distance parameter (non-crossvalidated R²=0.99) was approximately equal to the CoMFA without a distance parameter (non-cross-validated R²=0.98). Validation of our CoMFA was accomplished by predicting the K_M values of 15 diverse CYP2D6 substrates not in the original training set resulting in a predictive R²=0.62. Finally, we also pursued correlations of pK_a and log P with CYP2D6 substrate K_M in an effort to investigate other physicochemical properties.