Stand dynamics modelling approaches for multicohort management of eastern Canadian boreal forests

Abstract

Treeline dynamics in Siberia under changing climates as inferred from an individual-based model for Larix, supplementary material

Siberian boreal forests are expected to expand northwards in the course of global warming. However, processes of the treeline ecotone transition, as well astiming and related climate feedbacks are still not understood. Here, we present 'Larix Vegetation Simulator' LAVESI, an individual-based spatially-explicit model that can simulate Larix gmelinii (RUPR.) RUPR. stand dynamics in an attempt to improve our understanding about past and future treeline movements under changing climates. The relevant processes (growth, seed production and dispersal, establishment and mortality) are incorporated and adjusted to observation data mainly gained from the literature. Results of a local sensitivity analysis support the robustness of the model's parameterization by giving relatively small sensitivity values. We tested the model by simulating tree stands under modern climate across the whole Taymyr Peninsula, north-central Siberia (c. 64-80° N; 92-119° E). We find tree densities similar to observed forests in the northern to mid-treeline areas, but densities are overestimated in the southern parts of the simulated region. Finally, from a temperature-forcing experiment, we detect that the responses of tree stands lag the hypothetical warming by several decades, until the end of 21st century. With our simulation experiments we demonstrate that the newly-developed model captures the dynamics of the Siberian latitudinal treeline.

Critical evaluation of diameter increment modelling in the Great Lakes region

Simulations of forest stand dynamics in a modelling framework including Forest Vegetation Simulator (FVS) are diameter driven, thus the diameter or basal area increment model needs a special attention. This dissertation critically evaluates diameter or basal area increment models and modelling approaches in the context of the Great Lakes region of the United States and Canada. A set of related studies are presented that critically evaluate the sub-model for change in individual tree basal diameter used in the Forest Vegetation Simulator (FVS), a dominant forestry model in the Great Lakes region. Various historical implementations of the STEMS (Stand and Tree Evaluation and Modeling System) family of diameter increment models, including the current public release of the Lake States variant of FVS (LS-FVS), were tested for the 30 most common tree species using data from the Michigan Forest Inventory and Analysis (FIA) program. The results showed that current public release of the LS-FVS diameter increment model over-predicts 10-year diameter increment by 17% on average. Also the study affirms that a simple adjustment factor as a function of a single predictor, dbh (diameter at breast height) used in the past versions, provides an inadequate correction of model prediction bias. In order to re-engineer the basal diameter increment model, the historical, conceptual and philosophical differences among the individual tree increment model families and their modelling approaches were analyzed and discussed. Two underlying conceptual approaches toward diameter or basal area increment modelling have been often used: the potential-modifier (POTMOD) and composite (COMP) approaches, which are exemplified by the STEMS/TWIGS and Prognosis models, respectively. It is argued that both approaches essentially use a similar base function and neither is conceptually different from a biological perspective, even though they look different in their model forms. No matter what modelling approach is used, the base function is the foundation of an increment model. Two base functions – gamma and Box-Lucas – were identified as candidate base functions for forestry applications. The results of a comparative analysis of empirical fits showed that quality of fit is essentially similar, and both are sufficiently detailed and flexible for forestry applications. The choice of either base function in order to model diameter or basal area increment is dependent upon personal preference; however, the gamma base function may be preferred over the Box-Lucas, as it fits the periodic increment data in both a linear and nonlinear composite model form. Finally, the utility of site index as a predictor variable has been criticized, as it has been widely used in models for complex, mixed species forest stands though not well suited for this purpose. An alternative to site index in an increment model was explored, using site index and a combination of climate variables and Forest Ecosystem Classification (FEC) ecosites and data from the Province of Ontario, Canada. The results showed that a combination of climate and FEC ecosites variables can replace site index in the diameter increment model.

Hurricane disturbance and recovery of energy balance, CO2 fluxes and canopy structure in a mangrove forest of the Florida Everglades

Eddy covariance (EC) estimates of carbon dioxide (CO2) fluxes and energy balance are examined to investigate the functional responses of a mature mangrove forest to a disturbance generated by Hurricane Wilma on October 24, 2005 in the Florida Everglades. At the EC site, high winds from the hurricane caused nearly 100% defoliation in the upper canopy and widespread tree mortality. Soil temperatures down to -50 cm increased, and air temperature lapse rates within the forest canopy switched from statically stable to statically unstable conditions following the disturbance. Unstable conditions allowed more efficient transport of water vapor and CO2 from the surface up to the upper canopy layer. Significant increases in latent heat fluxes (LE) and nighttime net ecosystem exchange (NEE) were also observed and sensible heat fluxes (H) as a proportion of net radiation decreased significantly in response to the disturbance. Many of these impacts persisted through much of the study period through 2009. However, local albedo and MODIS (Moderate Resolution Imaging Spectro-radiometer) data (the Enhanced Vegetation Index) indicated a substantial proportion of active leaf area recovered before the EC measurements began 1 year after the storm. Observed changes in the vertical distribution and the degree of clumping in newly emerged leaves may have affected the energy balance. Direct comparisons of daytime NEE values from before the storm and after our measurements resumed did not show substantial or consistent differences that could be attributed to the disturbance. Regression analyses on seasonal time scales were required to differentiate the storm's impact on monthly average daytime NEE from the changes caused by interannual variability in other environmental drivers. The effects of the storm were apparent on annual time scales, and CO2 uptake remained approximately 250 g C m-2 yr-1 lower in 2009 compared to the average annual values measured in 2004-2005. Dry season CO2 uptake was relatively more affected by the disturbance than wet season values. Complex leaf regeneration dynamics on damaged trees during ecosystem recovery are hypothesized to lead to the variable dry versus wet season impacts on daytime NEE. In contrast, nighttime CO2 release (i.e., nighttime respiration) was consistently and significantly greater, possibly as a result of the enhanced decomposition of litter and coarse woody debris generated by the storm, and this effect was most apparent in the wet seasons compared to the dry seasons. The largest pre- and post-storm differences in NEE coincided roughly with the delayed peak in cumulative mortality of stems in 2007-2008. Across the hurricane-impacted region, cumulative tree mortality rates were also closely correlated with declines in peat surface elevation. Mangrove forest-atmosphere interactions are interpreted with respect to the damage and recovery of stand dynamics and soil accretion processes following the hurricane.

Chilling damage in a changing climate in coastal landscapes of the subtropical zone: a case study from south Florida

Freeze events significantly influence landscape structure and community composition along subtropical coastlines. This is particularly true in south Florida, where such disturbances have historically contributed to patch diversity within the mangrove forest, and have played a part in limiting its inland transgression. With projected increases in mean global temperatures, such instances are likely to become much less frequent in the region, contributing to a reduction in heterogeneity within the mangrove forest itself. To understand the process more clearly, we explored the dynamics of a Dwarf mangrove forest following two chilling events that produced freeze-like symptoms, i.e., leaf browning, desiccation, and mortality, and interpreted the resulting changes within the context of current winter temperatures and projected future scenarios. Structural effects from a 1996 chilling event were dramatic, with mortality and tissue damage concentrated among individuals comprising the Dwarf forest's low canopy. This disturbance promoted understory plant development and provided an opportunity for Laguncularia racemosa to share dominance with Rhizophora mangle. Mortality due to the less severe 2001 event was greatest in the understory, probably because recovery of the protective canopy following the earlier freeze was still incomplete. Stand dynamics were static over the same period in nearby unimpacted sites. The probability of reaching temperatures as low as those recorded at a nearby meteorological station (≤3 °C) under several warming scenarios was simulated by applying 1° incremental temperature increases to a model developed from a 42-year temperature record. According to the model, the frequency of similar chilling events decreased from once every 1.9 years at present to once every 3.4 and 32.5 years with 1 and 4 °C warming, respectively. The large decrease in the frequency of these events would eliminate an important mechanism that maintains Dwarf forest structure, and promotes compositional diversity.

Improved models for estimating temporal changes in carbon sequestration in above-ground biomass of mixed-species environmental plantings

Plantings of mixed native species (termed 'environmental plantings') are increasingly being established for carbon sequestration whilst providing additional environmental benefits such as biodiversity and water quality. In Australia, they are currently one of the most common forms of reforestation. Investment in establishing and maintaining such plantings relies on having a cost-effective modelling approach to providing unbiased estimates of biomass production and carbon sequestration rates. In Australia, the Full Carbon Accounting Model (FullCAM) is used for both national greenhouse gas accounting and project-scale sequestration activities. Prior to undertaking the work presented here, the FullCAM tree growth curve was not calibrated specifically for environmental plantings and generally under-estimated their biomass. Here we collected and analysed above-ground biomass data from 605 mixed-species environmental plantings, and tested the effects of several planting characteristics on growth rates. Plantings were then categorised based on significant differences in growth rates. Growth of plantings differed between temperate and tropical regions. Tropical plantings were relatively uniform in terms of planting methods and their growth was largely related to stand age, consistent with the un-calibrated growth curve. However, in temperate regions where plantings were more variable, key factors influencing growth were planting width, stand density and species-mix (proportion of individuals that were trees). These categories provided the basis for FullCAM calibration. Although the overall model efficiency was only 39-46%, there was nonetheless no significant bias when the model was applied to the various planting categories. Thus, modelled estimates of biomass accumulation will be reliable on average, but estimates at any particular location will be uncertain, with either under- or over-prediction possible. When compared with the un-calibrated yield curves, predictions using the new calibrations show that early growth is likely to be more rapid and total above-ground biomass may be higher for many plantings at maturity. This study has considerably improved understanding of the patterns of growth in different types of environmental plantings, and in modelling biomass accumulation in young (<25. years old) plantings. However, significant challenges remain to understand longer-term stand dynamics, particularly with temporal changes in stand density and species composition. © 2014.

Free atmospheric CO2 enrichment increased above ground biomass but did not affect symbiotic N2-fixation and soil carbon dynamics in a mixed deciduous stand in Wales

Through increases in net primary production (NPP), elevated CO2 is hypothesizes to increase the amount of plant litter entering the soil. The fate of this extra carbon on the forest floor or in mineral soil is currently not clear. Moreover, increased rates of NPP can be maintained only if forests can escape nitrogen limitation. In a Free atmospheric CO2 Enrichment (FACE) experiment near Bangor, Wales, 4 ambient CO2 and 4 FACE plots were planted with patches of Betula pendula, Alnus glutinosa and Fagus sylvatica on a former arable field. Four years after establishment, only a shallow L forest floor litter layer had formed due to intensive bioturbation. Total soil C and N contents increased irrespective of treatment and species as a result of afforestation. We could not detect an additional C sink in the soil, nor were soil C stabilization processes affected by FACE. We observed a decrease of leaf N content in Betula and Alnus under FACE, while the soil C/N ratio decreased regardless of CO2 treatment. The ratio of N taken up from the soil and by N2-fixation in Alnus was not affected by FACE. We infer that increased nitrogen use efficiency is the mechanism by which increased NPP is sustained under elevated CO2 at this site.

Cation dynamics in soils with different salinity levels growing irrigated rice

Salinity levels in soils of the Outer Coastal Plain of Rio Grande do Sul, Brazil, can be high, due to excess of Na in the irrigation water, evapotranspiration and soil development from marine sediments. The cultivation of irrigated rice could be an alternative, since ion uptake as well as leaching by the establishment of a water layer could mitigate the effects of soil salinity. This study aimed to evaluate the dynamics of basic cations in the solution of Albaqualf soils with different salinity levels growing irrigated rice. The plow layer contained exchangeable Na percentages (ESP) of 5.6, 9.0, 21.2 and 32.7 %. The plant stand, dry matter, Na, K and Ca + Mg uptake at full flowering and grain yield were evaluated. The levels of Na, K, Ca + Mg and electrical conductivity (EC) in the soil solution were also measured weekly during the rice cycle at four soil depths, in the water layer and irrigation water. The Na, K and Ca + Mg uptake by rice at full flowering was used to estimate ion depletion from the layer under root influence. Soil salinity induced a reduction in the rice stand, especially in the soil with ESP of 32.7 %, resulting in lower cation uptake and very low yield at that site. As observed in the water layer and irrigation water, the Na, K, Ca + Mg and EC levels in the soil solution decreased with time at depths of 5, 10 and 20 cm, regardless of the original soil salinity, showing that cation dynamics in the plow layer was determined by leaching and root uptake, rather than by the effect of evapoconcentration of basic cations in the soil surface layer.

Dynamics and vertical distribution of a Hrabeiella periglandulata (Annelida) population in South Moravia, Czech Republic

The objective of this work was to evaluate the density, dynamics and vertical distribution of a Hrabeiella periglandulata population in a forest soil at Brno, Czech Republic. From December 2003 to November 2004, two plots covered by mixed stands and two covered by coniferous stands were sampled monthly. Six soil cores per plot were taken down to 15 cm and subdivided into layers, which were subjected to wet funnel extraction. Missing in one of the coniferous stands H. periglandulata was abundant in the mixed stand with the highest soil pH. In this stand, monthly sampling continued until November 2005, with three additional samplings up to January 2007. Mean annual density was 2,672±1,534 individuals m-2. Population dynamics differed from those reported from Germany. Highest densities were reached in early summer, lowest between August and December. Due to aggregated horizontal distribution, differences between monthly values were often nonsignificant. No significant correlation with climatic data was found. Nevertheless, the observed dynamics corresponded to the climatic conditions, showing particularly the negative effect of drought. The population was evenly distributed in the sampled soil profile, only avoiding the organic layer. Except for a locality in Poland, this is the easternmost record of the species.

Modelling of the natural chlorine cycling in a coniferous stand: implications for chlorine-36 behaviour in a contaminated forest environment

Considered as one of the most available radionuclide in soileplant system, 36Cl is of potential concern for long-term management of radioactive wastes, due to its high mobility and its long half-life. To evaluate the risk of dispersion and accumulation of 36Cl in the biosphere as a consequence of a potential contamination, there is a need for an appropriate understanding of the chlorine cycling dynamics in the ecosystems. To date, a small number of studies have investigated the chlorine transfer in the ecosystem including the transformation of chloride to organic chlorine but, to our knowledge, none have modelled this cycle. In this study, a model involving inorganic as well as organic pools in soils has been developed and parameterised to describe the biogeochemical fate of chlorine in a pine forest. The model has been evaluated for stable chlorine by performing a range of sensitivity analyses and by comparing the simulated to the observed values. Finally a range of contamination scenarios, which differ in terms of external supply, exposure time and source, has been simulated to estimate the possible accumulation of 36Cl within the different compartments of the coniferous stand. The sensitivity study supports the relevancy of the model and its compartments, and has highlighted the chlorine transfers affecting the most the residence time of chlorine in the stand. Compared to observations, the model simulates realistic values for the chlorine content within the different forest compartments. For both atmospheric and underground contamination scenarios most of the chlorine can be found in its organic form in the soil. However, in case of an underground source, about two times less chlorine accumulates in the system and proportionally more chlorine leaves the system through drainage than through volatilisation.

Dynamics of entangled polymers

In this book chapter, we focus on the basic questions on dynamics of entangled polymers that faced the earliest researchers, and see how their answers stand up in light of what has been learned in the last 40 years. In the process, we hope to provide the novice with a basic understanding of the physical concepts without burdening the reader with the mathematics, or with detailed data that can easily be found in the books referred to above. At the same time, we hope to update the knowledgeable non-specialist with the latest developments in the area. For the expert, we hope to provoke new ways of thinking about an old, but incompletely solved, problem area.

Dynamics of soil exploration by fine roots down to a depth of 10 m throughout the entire rotation in Eucalyptus grandis plantations

Although highly weathered soils cover considerable areas in tropical regions, little is known about exploration by roots in deep soil layers. Intensively managed Eucalyptus plantations are simple forest ecosystems that can provide an insight into the belowground growth strategy of fast-growing tropical trees. Fast exploration of deep soil layers by eucalypt fine roots may contribute to achieving a gross primary production that is among the highest in the world for forests. Soil exploration by fine roots down to a depth of 10 m was studied throughout the complete cycle in Eucalyptus grandis plantations managed in short rotation. Intersects of fine roots, less than 1 mm in diameter, and medium-sized roots, 1-3 mm in diameter, were counted on trench walls in a chronosequence of 1-, 2-, 3.5-, and 6-year-old plantations on a sandy soil, as well as in an adjacent 6-year-old stand growing in a clayey soil. Two soil profiles were studied down to a depth of 10 m in each stand (down to 6 m at ages 1 and 2 years) and 4 soil profiles down to 1.5-3.0 m deep. The root intersects were counted on 224 m(2) of trench walls in 15 pits. Monitoring the soil water content showed that, after clear cutting, almost all the available water stored down to a depth of 7 m was taken up by tree roots within 1.1 year of planting. The soil space was explored intensively by fine roots down to a depth of 3 m from 1 year after planting, with an increase in anisotropy in the upper layers throughout the rotation. About 60% of fine root intersects were found at a depth of more than 1 m, irrespective of stand age. The root distribution was isotropic in deep soil layers and kriged maps showed fine root clumping. A considerable volume of soil was explored by fine roots in eucalypt plantations on deep tropical soils, which might prevent water and nutrient losses by deep drainage after canopy closure and contribute to maximizing resource uses.

A new probabilistic canopy dynamics model (SLCD) that is suitable for evergreen and deciduous forest ecosystems

There are strong uncertainties regarding LAI dynamics in forest ecosystems in response to climate change. While empirical growth & yield models (G&YMs) provide good estimations of tree growth at the stand level on a yearly to decennial scale, process-based models (PBMs) use LAI dynamics as a key variable for enabling the accurate prediction of tree growth over short time scales. Bridging the gap between PBMs and G&YMs could improve the prediction of forest growth and, therefore, carbon, water and nutrient fluxes by combining modeling approaches at the stand level.Our study aimed to estimate monthly changes of leaf area in response to climate variations from sparse measurements of foliage area and biomass. A leaf population probabilistic model (SLCD) was designed to simulate foliage renewal. The leaf population was distributed in monthly cohorts, and the total population size was limited depending on forest age and productivity. Foliage dynamics were driven by a foliation function and the probabilities ruling leaf aging or fall. Their formulation depends on the forest environment.The model was applied to three tree species growing under contrasting climates and soil types. In tropical Brazilian evergreen broadleaf eucalypt plantations, the phenology was described using 8 parameters. A multi-objective evolutionary algorithm method (MOEA) was used to fit the model parameters on litterfall and LAI data over an entire stand rotation. Field measurements from a second eucalypt stand were used to validate the model. Seasonal LAI changes were accurately rendered for both sites (R-2 = 0.898 adjustment, R-2 = 0.698 validation). Litterfall production was correctly simulated (R-2 = 0.562, R-2 = 0.4018 validation) and may be improved by using additional validation data in future work. In two French temperate deciduous forests (beech and oak), we adapted phenological sub-modules of the CASTANEA model to simulate canopy dynamics, and SLCD was validated using LAI measurements. The phenological patterns were simulated with good accuracy in the two cases studied. However, IA/max was not accurately simulated in the beech forest, and further improvement is required.Our probabilistic approach is expected to contribute to improving predictions of LAI dynamics. The model formalism is general and suitable to broadleaf forests for a large range of ecological conditions. (C) 2014 Elsevier B.V. All rights reserved.

From inverted pendulum to N-link chain: inertial sensors-based assessment of movement kinematics and dynamics for functional evaluation and rehabilitation

Despite several clinical tests that have been developed to qualitatively describe complex motor tasks by functional testing, these methods often depend on clinicians' interpretation, experience and training, which make the assessment results inconsistent, without the precision required to objectively assess the effect of the rehabilitative intervention. A more detailed characterization is required to fully capture the various aspects of motor control and performance during complex movements of lower and upper limbs. The need for cost-effective and clinically applicable instrumented tests would enable quantitative assessment of performance on a subject-specific basis, overcoming the limitations due to the lack of objectiveness related to individual judgment, and possibly disclosing subtle alterations that are not clearly visible to the observer. Postural motion measurements at additional locations, such as lower and upper limbs and trunk, may be necessary in order to obtain information about the inter-segmental coordination during different functional tests involved in clinical practice. With these considerations in mind, this Thesis aims: i) to suggest a novel quantitative assessment tool for the kinematics and dynamics evaluation of a multi-link kinematic chain during several functional motor tasks (i.e. squat, sit-to-stand, postural sway), using one single-axis accelerometer per segment, ii) to present a novel quantitative technique for the upper limb joint kinematics estimation, considering a 3-link kinematic chain during the Fugl-Meyer Motor Assessment and using one inertial measurement unit per segment. The suggested methods could have several positive feedbacks from clinical practice. The use of objective biomechanical measurements, provided by inertial sensor-based technique, may help clinicians to: i) objectively track changes in motor ability, ii) provide timely feedback about the effectiveness of administered rehabilitation interventions, iii) enable intervention strategies to be modified or changed if found to be ineffective, and iv) speed up the experimental sessions when several subjects are asked to perform different functional tests.