Simulation of tree-ring widths with a model for primary production, carbon allocation, and growth

We present a simple, generic model of annual tree growth, called "T". This model accepts input from a first-principles light-use efficiency model (the "P" model). The P model provides values for gross primary production (GPP) per unit of absorbed photosynthetically active radiation (PAR). Absorbed PAR is estimated from the current leaf area. GPP is allocated to foliage, transport tissue, and fine-root production and respiration in such a way as to satisfy well-understood dimensional and functional relationships. Our approach thereby integrates two modelling approaches separately developed in the global carbon-cycle and forest-science literature. The T model can represent both ontogenetic effects (the impact of ageing) and the effects of environmental variations and trends (climate and CO2) on growth. Driven by local climate records, the model was applied to simulate ring widths during the period 1958–2006 for multiple trees of Pinus koraiensis from the Changbai Mountains in northeastern China. Each tree was initialised at its actual diameter at the time when local climate records started. The model produces realistic simulations of the interannual variability in ring width for different age cohorts (young, mature, and old). Both the simulations and observations show a significant positive response of tree-ring width to growing-season total photosynthetically active radiation (PAR0) and the ratio of actual to potential evapotranspiration (α), and a significant negative response to mean annual temperature (MAT). The slopes of the simulated and observed relationships with PAR0 and α are similar; the negative response to MAT is underestimated by the model. Comparison of simulations with fixed and changing atmospheric CO2 concentration shows that CO2 fertilisation over the past 50 years is too small to be distinguished in the ring-width data, given ontogenetic trends and interannual variability in climate.

Nonlinear identification using orthogonal forward regression with nested optimal regularization

An efficient data based-modeling algorithm for nonlinear system identification is introduced for radial basis function (RBF) neural networks with the aim of maximizing generalization capability based on the concept of leave-one-out (LOO) cross validation. Each of the RBF kernels has its own kernel width parameter and the basic idea is to optimize the multiple pairs of regularization parameters and kernel widths, each of which is associated with a kernel, one at a time within the orthogonal forward regression (OFR) procedure. Thus, each OFR step consists of one model term selection based on the LOO mean square error (LOOMSE), followed by the optimization of the associated kernel width and regularization parameter, also based on the LOOMSE. Since like our previous state-of-the-art local regularization assisted orthogonal least squares (LROLS) algorithm, the same LOOMSE is adopted for model selection, our proposed new OFR algorithm is also capable of producing a very sparse RBF model with excellent generalization performance. Unlike our previous LROLS algorithm which requires an additional iterative loop to optimize the regularization parameters as well as an additional procedure to optimize the kernel width, the proposed new OFR algorithm optimizes both the kernel widths and regularization parameters within the single OFR procedure, and consequently the required computational complexity is dramatically reduced. Nonlinear system identification examples are included to demonstrate the effectiveness of this new approach in comparison to the well-known approaches of support vector machine and least absolute shrinkage and selection operator as well as the LROLS algorithm.

Estimation of Gaussian process regression model using probability distance measures

A new class of parameter estimation algorithms is introduced for Gaussian process regression (GPR) models. It is shown that the integration of the GPR model with probability distance measures of (i) the integrated square error and (ii) Kullback–Leibler (K–L) divergence are analytically tractable. An efficient coordinate descent algorithm is proposed to iteratively estimate the kernel width using golden section search which includes a fast gradient descent algorithm as an inner loop to estimate the noise variance. Numerical examples are included to demonstrate the effectiveness of the new identification approaches.

Tensor regression based on linked multiway parameter analysis

Classical regression methods take vectors as covariates and estimate the corresponding vectors of regression parameters. When addressing regression problems on covariates of more complex form such as multi-dimensional arrays (i.e. tensors), traditional computational models can be severely compromised by ultrahigh dimensionality as well as complex structure. By exploiting the special structure of tensor covariates, the tensor regression model provides a promising solution to reduce the model’s dimensionality to a manageable level, thus leading to efficient estimation. Most of the existing tensor-based methods independently estimate each individual regression problem based on tensor decomposition which allows the simultaneous projections of an input tensor to more than one direction along each mode. As a matter of fact, multi-dimensional data are collected under the same or very similar conditions, so that data share some common latent components but can also have their own independent parameters for each regression task. Therefore, it is beneficial to analyse regression parameters among all the regressions in a linked way. In this paper, we propose a tensor regression model based on Tucker Decomposition, which identifies not only the common components of parameters across all the regression tasks, but also independent factors contributing to each particular regression task simultaneously. Under this paradigm, the number of independent parameters along each mode is constrained by a sparsity-preserving regulariser. Linked multiway parameter analysis and sparsity modeling further reduce the total number of parameters, with lower memory cost than their tensor-based counterparts. The effectiveness of the new method is demonstrated on real data sets.

A model analysis of climate and CO2 controls on tree growth in a semi-arid woodland

We used a light-use efficiency model of photosynthesis coupled with a dynamic carbon allocation and tree-growth model to simulate annual growth of the gymnosperm Callitris columellaris in the semi-arid Great Western Woodlands, Western Australia, over the past 100 years. Parameter values were derived from independent observations except for sapwood specific respiration rate, fine-root turnover time, fine-root specific respiration rate and the ratio of fine-root mass to foliage area, which were estimated by Bayesian optimization. The model reproduced the general pattern of interannual variability in radial growth (tree-ring width), including the response to the shift in precipitation regimes that occurred in the 1960s. Simulated and observed responses to climate were consistent. Both showed a significant positive response of tree-ring width to total photosynthetically active radiation received and to the ratio of modeled actual to equilibrium evapotranspiration, and a significant negative response to vapour pressure deficit. However, the simulations showed an enhancement of radial growth in response to increasing atmospheric CO2 concentration (ppm) ([CO2]) during recent decades that is not present in the observations. The discrepancy disappeared when the model was recalibrated on successive 30-year windows. Then the ratio of fine-root mass to foliage area increases by 14% (from 0.127 to 0.144 kg C m-2) as [CO2] increased while the other three estimated parameters remained constant. The absence of a signal of increasing [CO2] has been noted in many tree-ring records, despite the enhancement of photosynthetic rates and water-use efficiency resulting from increasing [CO2]. Our simulations suggest that this behaviour could be explained as a consequence of a shift towards below-ground carbon allocation.

Long-term conditioning of soil by plantation eucalypts and pines does not affect growth of the native jarrah tree

Plant species can condition the physico-chemical and biological properties of soil in ways that modify plant growth via plant–soil feedback (PSF). Plant growth can be positively affected, negatively affected or neutrally affected by soil conditioning by the same or other plant species. Soil conditioning by other plant species has particular relevance to ecological restoration of historic ecosystems because sites set aside for restoration are often conditioned by other, potentially non-native, plant species. We investigated changes in properties of jarrah forest soils after long-term (35 years) conditioning by pines (Pinus radiata), Sydney blue gums (Eucalyptus saligna), both non-native, plantation trees, and jarrah (Eucalyptus marginata; dominant native tree). Then, we tested the influence of the conditioned soils on the growth of jarrah seedlings. Blue gums and pines similarly conditioned the physico-chemical properties of soils, which differed from soil conditioning caused by jarrah. Especially important were the differences in conditioning of the properties C:N ratio, pH, and available K. The two eucalypt species similarly conditioned the biological properties of soil (i.e. community level physiological profile, numbers of fungal-feeding nematodes, omnivorous nematodes, and nematode channel ratio), and these differed from conditioning caused by pines. Species-specific conditioning of soil did not translate into differences in the amounts of biomass produced by jarrah seedlings and a neutral PSF was observed. In summary, we found that decades of soil conditioning by non-native plantation trees did not influence the growth of jarrah seedlings and will therefore not limit restoration of jarrah following the removal of the plantation trees.

Struggles over family land? Tree crops, land and labour in Ghana's Brong-Ahafo region

Agricultural land use in much of Brong-Ahafo region, Ghana has been shifting from the production of food crops towards increased cashew nut cultivation in recent years. This article explores everyday, less visible, gendered and generational struggles over family farms in West Africa, based on qualitative, participatory research in a rural community that is becoming increasingly integrated into the global capitalist system. As a tree crop, cashew was regarded as an individual man's property to be passed on to his wife and children rather than to extended family members, which differed from the communal land tenure arrangements governing food crop cultivation. The tendency for land, cash crops and income to be controlled by men, despite women's and young people's significant labour contributions to family farms, and for women to rely on food crop production for their main source of income and for household food security, means that women and girls are more likely to lose out when cashew plantations are expanded to the detriment of land for food crops. Intergenerational tensions emerged when young people felt that their parents and elders were neglecting their views and concerns. The research provides important insights into gendered and generational power relations regarding land access, property rights and intra-household decision-making processes. Greater dialogue between genders and generations may help to tackle unequal power relations and lead to shared decision-making processes that build the resilience of rural communities.

Leaf trapping and retention of particles by holm oak and other common tree species in Mediterranean urban environments

Holm oak (Quercus ilex), a widespread urban street tree in the Mediterranean region, is widely used as biomonitor of persistent atmospheric pollutants, especially particulate-bound metals. By using lab- and field-based experimental approaches, we compared the leaf-level capacity for particles’ capture and retention between Q. ilex and other common Mediterranean urban trees: Quercus cerris, Platanus × hispanica, Tilia cordata and Olea europaea. All applied methods were effective in quantifying particulate capture and retention, although not univocal in ranking species performances. Distinctive morphological features of leaves led to differences in species’ ability to trap and retain particles of different size classes and to accumulate metals after exposure to traffic in an urban street. Overall, P. × hispanica and T. cordata showed the largest capture potential per unit leaf area for most model particles (Na+ and powder particles), and street-level Cu and Pb, while Q. ilex acted intermediately. After wash-off experiments, P. × hispanica leaves had the greatest retention capacity among the tested species and O. europaea the lowest. We concluded that the Platanus planting could be considered in Mediterranean urban environments due to its efficiency in accumulating and retaining airborne particulates; however, with atmospheric pollution being typically higher in winter, the evergreen Q. ilex represents a better year-round choice to mitigate the impact of airborne particulate pollutants.

Probabilistic wind power forecasts with an inverse power curve transformation and censored regression

Forecasting wind power is an important part of a successful integration of wind power into the power grid. Forecasts with lead times longer than 6 h are generally made by using statistical methods to post-process forecasts from numerical weather prediction systems. Two major problems that complicate this approach are the non-linear relationship between wind speed and power production and the limited range of power production between zero and nominal power of the turbine. In practice, these problems are often tackled by using non-linear non-parametric regression models. However, such an approach ignores valuable and readily available information: the power curve of the turbine's manufacturer. Much of the non-linearity can be directly accounted for by transforming the observed power production into wind speed via the inverse power curve so that simpler linear regression models can be used. Furthermore, the fact that the transformed power production has a limited range can be taken care of by employing censored regression models. In this study, we evaluate quantile forecasts from a range of methods: (i) using parametric and non-parametric models, (ii) with and without the proposed inverse power curve transformation and (iii) with and without censoring. The results show that with our inverse (power-to-wind) transformation, simpler linear regression models with censoring perform equally or better than non-linear models with or without the frequently used wind-to-power transformation.

Tests of sunspot number sequences: 3. Effects of regression procedures on the calibration of historic sunspot data

We use sunspot group observations from the Royal Greenwich Observatory (RGO) to investigate the effects of intercalibrating data from observers with different visual acuities. The tests are made by counting the number of groups RB above a variable cut-off threshold of observed total whole-spot area (uncorrected for foreshortening) to simulate what a lower acuity observer would have seen. The synthesised annual means of RB are then re-scaled to the full observed RGO group number RA using a variety of regression techniques. It is found that a very high correlation between RA and RB (rAB > 0.98) does not prevent large errors in the intercalibration (for example sunspot maximum values can be over 30 % too large even for such levels of rAB). In generating the backbone sunspot number (RBB), Svalgaard and Schatten (2015, this issue) force regression fits to pass through the scatter plot origin which generates unreliable fits (the residuals do not form a normal distribution) and causes sunspot cycle amplitudes to be exaggerated in the intercalibrated data. It is demonstrated that the use of Quantile-Quantile (“Q  Q”) plots to test for a normal distribution is a useful indicator of erroneous and misleading regression fits. Ordinary least squares linear fits, not forced to pass through the origin, are sometimes reliable (although the optimum method used is shown to be different when matching peak and average sunspot group numbers). However, other fits are only reliable if non-linear regression is used. From these results it is entirely possible that the inflation of solar cycle amplitudes in the backbone group sunspot number as one goes back in time, relative to related solar-terrestrial parameters, is entirely caused by the use of inappropriate and non-robust regression techniques to calibrate the sunspot data.

Effects of selective logging on tropical forest tree growth

We combined measurements of tree growth and carbon dioxide exchange to investigate the effects of selective logging on the Aboveground Live Biomass (AGLB) of a tropical rain forest in the Amazon. Most of the measurements began at least 10 months before logging and continued at least 36 months after logging. The logging removed similar to 15% of the trees with Diameter at Breast Height (DBH) greater than 35 cm, which resulted in an instantaneous 10% reduction in AGLB. Both wood production and mortality increased following logging, while Gross Primary Production (GPP) was unchanged. The ratio of wood production to GPP (the wood Carbon Use Efficiency or wood CUE) more than doubled following logging. Small trees (10 cm < DBH < 35 cm) accounted for most of the enhanced wood production. Medium trees (35 cm < DBH < 55 cm) that were within 30 m of canopy gaps created by the logging also showed increased growth. The patterns of enhanced growth are most consistent with logging-induced increases in light availability. The AGLB continued to decline over the study, as mortality outpaced wood production. Wood CUE and mortality remained elevated throughout the 3 years of postlogging measurements. The future trajectory of AGLB and the forest`s carbon balance are uncertain, and will depend on how long it takes for heterotrophic respiration, mortality, and CUE to return to prelogging levels.

Effects of abscisic acid, ethylene and sugars on the mobilization of storage proteins and carbohydrates in seeds of the tropical tree Sesbania virgata (Leguminosae)

Endospermic legumes are abundant in tropical forests and their establishment is closely related to the mobilization of cell-wall storage polysaccharides. Endosperm cells also store large numbers of protein bodies that play an important role as a nitrogen reserve in this seed. In this work, a systems approach was adopted to evaluate some of the changes in carbohydrates and hormones during the development of seedlings of the rain forest tree Sesbania virgata during the period of establishment. Seeds imbibed abscisic acid (ABA), glucose and sucrose in an atmosphere of ethylene, and the effects of these compounds on the protein contents, alpha-galactosidase activity and endogenous production of ABA and ethylene by the seeds were observed. The presence of exogenous ABA retarded the degradation of storage protein in the endosperm and decreased alpha-galactosidase activity in the same tissue during galactomannan degradation, suggesting that ABA represses enzyme action. On the other hand, exogenous ethylene increased alpha-galactosidase activity in both the endosperm and testa during galactomannan degradation, suggesting an inducing effect of this hormone on the hydrolytic enzymes. Furthermore, the detection of endogenous ABA and ethylene production during the period of storage mobilization and the changes observed in the production of these endogenous hormones in the presence of glucose and sucrose, suggested a correlation between the signalling pathway of these hormones and the sugars. These findings suggest that ABA, ethylene and sugars play a role in the control of the hydrolytic enzyme activities in seeds of S. virgata, controlling the process of storage degradation. This is thought to ensure a balanced flow of the carbon and nitrogen for seedling development.

Cell wall polysaccharides from cell suspension cultures of the Atlantic Forest tree Rudgea jasminoides (Rubiaceae)

Rudgea jasminoides (Rubiaceae) is a tropical tree species native of the Atlantic Forest in the south of Brazil. Previous studies with leaf cell walls of R. jasminoides showed a different proportion of cross-linked glycans compared to what is usually reported for eudicots. However, due to the difficulties of working with whole plant organs, cell suspensions of R. jasminoides, consisting of predominantly undifferentiated cells with mainly primary cell walls, were used to examine cell walls and extracellular soluble polysaccharides (EP) released into the culture medium. Sugar composition and linkage analysis showed homogalacturonans, xylogalacturonans and arabinogalactans to be the predominant EP. In the cell wall, homogalacturonans and arabinogalactans are the major pectins, and xyloglucans and xylans are the major cross-linking glycans. The presence of xylogalacturonans in the R. jasminoides cell cultures seems to be related to the occurrence of a homogeneous cell suspension with loosely attached cells. Although all alkali extractions from the cell walls yielded amounts of xyloglucan that exceed those of the xylans, the latter was found in a proportion that is higher than what has been usually reported for primary cell walls of most eudicots. The xyloglucan from cell walls of cell suspension cultures of R. jasminoides has low fucosylation levels and high proportion of galactosyl residues, a branching pattern commonly found in storage cell-wall xyloglucans.