Incorporating shape constraints in generalized additive modelling of the height-diameter relationship for Norway spruce

Background： Measurements of tree heights and diameters are essential in forest assessment and modelling. Tree heights are used for estimating timber volume, site index and other important variables related to forest growth and yield, succession and carbon budget models. However, the diameter at breast height （dbh） can be more accurately obtained and at lower cost, than total tree height. Hence, generalized height-diameter （h-d） models that predict tree height from dbh, age and other covariates are needed. For a more flexible but biologically plausible estimation of covariate effects we use shape constrained generalized additive models as an extension of existing h-d model approaches. We use causal site parameters such as index of aridity to enhance the generality and causality of the models and to enable predictions under projected changeable climatic conditions. Methods： We develop unconstrained generalized additive models （GAM） and shape constrained generalized additive models （SCAM） for investigating the possible effects of tree-specific parameters such as tree age, relative diameter at breast height, and site-specific parameters such as index of aridity and sum of daily mean temperature during vegetation period, on the h-d relationship of forests in Lower Saxony, Germany. Results： Some of the derived effects, e.g. effects of age, index of aridity and sum of daily mean temperature have significantly non-linear pattern. The need for using SCAM results from the fact that some of the model effects show partially implausible patterns especially at the boundaries of data ranges. The derived model predicts monotonically increasing levels of tree height with increasing age and temperature sum and decreasing aridity and social rank of a tree within a stand, The definition of constraints leads only to marginal or minor decline in the model statistics like AIC An observed structured spatial trend in tree height is modelled via 2-dimensional surface fitting. Conclusions： We demonstrate that the SCAM approach allows optimal regression modelling flexibility similar to the standard GAM but with the additional possibility of defining specific constraints for the model effects. The longitudinal character of the model allows for tree height imputation for the current status of forests but also for future tree height prediction.

Longitudinal height-diameter curves for Norway spruce, Scots pine and silver birch in Norway based on shape constraint additive regression models

Background： Generalized height-diameter curves based on a re-parameterized version of the Korf function for Norway spruce （Piceo abies （L.） Karst.）, Scots pine （Pinus sylvestris L.） and silver birch （Betula pendula Roth） in Norwa are presented. The Norwegian National Forest Inventory （NFI） is used as data base for estimating the model parameters. The derived models are developed to enable spatially explicit and site sensitive tree height imputatio in forest inventories as well as future tree height predictions in growth and yield scenario simulations. Methods： Generalized additive mixed models （gamm） are employed to detect and quantify potentially non-linear effects of predictor variables. In doing so the quadratic mean diameter serves as longitudinal covariate since stand ag as measured in the NFI, shows only a weak correlation with a stands developmental status in Norwegian forests. Additionally the models can be locally calibrated by predicting random effects if measured height-diameter pairs are available. Based on the model selection of non-constraint models, shape constraint additive models （scare） were fit tc incorporate expert knowledge and intrinsic relationships by enforcing certain effect patterns like monotonicity. Results： Model comparisons demonstrate that the shape constraints lead to only marginal differences in statistical characteristics but ensure reasonable model predictions. Under constant constraints the developed models predict increasing tree heights with decreasing altitude, increasing soil depth and increasing competition pressure of a tree. / two-dimensional spatially structured effect of UTM-coordinates accounts for the potential effects of large scale spatial correlated covariates, which were not at our disposal. The main result of modelling the spatially structured effect is lower tree height prediction for coastal sites and with increasing latitude. The quadratic mean diameter affects both the level and the slope of the height-diameter curve and both effects are positive. Conclusions： In this investigation it is assumed that model effects in additive modelling of height-diameter curves which are unfeasible and too wiggly from an expert point of view are a result of quantitatively or qualitatively limited data bases. However, this problem can be regarded not to be specific to our investigation but more general since growth and yield data that are balanced over the whole data range with respect to all combinations of predictor variables are exceptional cases. Hence, scare may provide methodological improvements in several applications by combining the flexibility of additive models with expert knowledge.

Constrained Control of Interpolating Surfaces by Parameters

In order to meet the needs of practical design, an interpolation technique is employed to constrain the shape of surfaces. The method of preserving positivity on the interpolation surface and constraint on interpolating data is also developed. The advantage of this new method is that it can be used to constrain the shape of an interpolating surface only by selecting suitable parameters, and numerical examples are presented to show the performance of the method.

First application of large reactivity measurement through rod drop based on three-dimensional space–time dynamics

Abstract Reactivity measurement is an essential part of a zero-power physics test,which is critical to reactor design and development.The rod drop experimental technique is used to measure the control rod worth in a zero-power physics test.The conventional rod drop experimental technique is limited by the spatial effect and the difference between the calculated static reactivity and measured dynamic reactivity;thus,the method must be improved.In this study,a modified rod drop experimental technique that constrains the detector neutron flux shape function based on three-dimensional space–time dynamics to reduce the reactivity perturbation and a new method for calculating the detector neutron flux shape function are proposed.Correction factors were determined using Monte Carlo N-particle transport code and transient analysis code for a pressurized water reactor at the Ulsan National Institute of Science and Technology and Xi’an Jiaotong University,and a large reactivity of over 2000 pcm was measured using the modified technique.This research evaluated the modified technique accuracy,studied the influence of the correction factors on the modification,and investigated the effect of constraining the shape function on the reactivity perturbation reduction caused by the difference between the calculated neutron flux and true value,using the new method to calculate the shape function of the detector neutron flux and avoiding the neutron detector response function(weighting factor)calculation.