Stand structural dynamics on pristine and managed boreal peatlands

The objectives of this study were to investigate the stand structure and succession dynamics in Scots pine (Pinus sylvestris L.) stands on pristine peatlands and in Scots pine and Norway spruce (Picea abies (L.) Karst.) dominated stands on drained peatlands. Furthermore, my focus was on characterising how the inherent and environmental factors and the intermediate thinnings modify the stand structure and succession. For pristine peatlands, the study was based on inventorial stand data, while for drained peatlands, longitudinal data from repeatedly measured stands were utilised. The studied sites covered the most common peatland site types in Finland. They were classified into two categories according to the ecohydrological properties related to microsite variation and nutrient levels within sites. Tree DBH and age distributions in relation to climate and site type were used to study the stand dynamics on pristine sites. On drained sites, the Weibull function was used to parameterise the DBH distributions and mixed linear models were constructed to characterise the impacts of different ecological factors on stand dynamics. On pristine peatlands, both climate and the ecohydrology of the site proved to be crucial factors determining the stand structure and its dynamics. Irrespective of the vegetation succession, enhanced site productivity and increased stand stocking they significantly affected the stand dynamics also on drained sites. On the most stocked sites on pristine peatlands the inter-tree competition seemed to also be a significant factor modifying stand dynamics. Tree age and size diversity increased with stand age, but levelled out in the long term. After drainage, the stand structural unevenness increased due to the regeneration and/or ingrowth of the trees. This increase was more pronounced on sparsely forested composite sites than on more fully stocked genuine forested sites in Scots pine stands, which further undergo the formation of birch and spruce undergrowth beneath the overstory as succession proceeds. At 20-30 years after drainage the structural heterogeneity started to decrease, indicating increased inter-tree competition, which increased the mortality of suppressed trees within stand. Peatland stands are more dynamic than anticipated and are generally not characterized by a balanced, self-perpetuating structure. On pristine sites, various successional pathways are possible, whereas on drained sites the succession has more uniform trend. Typically, stand succession proceeds without any distinct developmental stages on pristine peatlands, whereas on drained peatlands, at least three distinct stages could be identified. Thinnings had only little impact on the stand succession. The new information on stand dynamics may be utilised, e.g. in forest management planning to facilitate the allocation of the growth resources to the desired crop component by appropriate silvicultural treatments, as well as assist in assessing the effects of the climate change on the forested boreal peatlands.

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.

Kohti hyvää suometsien hoitoa : harvennusten ja kunnostusojitusten vaikutus ojitusaluemänniköiden puuntuotokseen ja metsänkasvatuksen taloustulokseen

This thesis examines the impacts of silvicultural activities on productivity and financial returns of Scots pine (Pinus sylvestris L.) stands on drained peatlands in Finland. The effects of ditch network maintenance operations (DNM) and thinnings, with different timings and intensities, were studied. Based on stand development simulations, the best regimes for different types of stands according to site type, climatic area, and stand silvicultural status were defined from the viewpoint of both wood production and financial profitability. Certain aspects affecting the management outcomes, such as the timing of the first thinning, were examined using data from thinning experiments. Long-term predictions of the impacts of different management regimes were carried out by simulating the development of well-representative model-stands which were composed from appropriate inventory data sets. The MOTTI stand simulator used to perform the simulations enables the predictions by utilizing specific models for drained peatland stands. In addition to natural stand dynamics, these models describe the effects of silvicultural treatments on the development of a given stand. The mean annual increment of merchantable wood (MAImerch) was used as the measure of wood productivity, and the financial feasibility of the regimes was compared using net present value (NPV) analysis. Silvicultural treatments, when applied to appropriately match stand condition, increased both the productivity and financial returns of stand management. Applying DNM resulted in a small increase in MAImerch. When thinning was introduced along with DNM, their combined effect on wood productivity was considerable. According to current operational practices, DNM is generally combined with thinning. In some cases, e.g., in sites of low productivity, the need for DNM may become apparent prior to the thinning stage. As for profitability, thinnings proved to be crucial. The regimes with heavy and late thinnings were generally more profitable than those with normal thinnings. Further, early thinning (relative to stand volume) lacked appeal when seeking a financially profitable removal from the first thinning. In young stands with an initially poor silvicultural condition, however, applying even a low-yielding first thinning considerably increased the NPV when compared to a regime with no thinning at all. Generally, the regimes resulting in the best profitability included heavier thinnings and fewer DNM and thinning treatments than did the regimes resulting in the best yield results. This study demonstrates considerable potential for profitable wood production-oriented management in pine stands on drained peatlands despite their challenging circumstances and long rotations. The results can be used for defining new and more site-specific silvicultural guidelines for various types of drained, pine-dominated peatland stands within the entire range of boreal conditions.

Landscape-scale soil erosion modelling and ecological restoration for a mountainous watershed in Sichuan, China

Sichuanissa Tiibetin ylängön metsäkato on pysähtynyt mutta eroosio-ongelmat jatkuvat Viikin tropiikki-instituutin tutkija Ping ZHOU kartoitti trooppisen metsänhoidon alaan kuuluvassa väitöskirjatyössään maaperän eroosioalttiutta ja sen riippuvuutta metsäkasvillisuudesta Jangtsen tärkeää sivuhaaraa Min-jokea ympäröivällä n. 7400 neliökilometrin suuruisella valuma-alueella Sichuanin Aba-piirikunnassa. Aineistonaan hän käytti muun muassa satelliittikartoitustietoja ja mittaustuloksia yli 600 maastokoealalta. Tutkimuksen nimi suomeksi on "Maaperän eroosion mallinnus ja vuoristoisen valuma-alueen ekologinen ennallistaminen Sichuanissa Kiinassa". Aikaisempien tutkimusten perusteella oli tiedossa että metsien häviäminen tällä alueella pysähtyi jo 1980-luvun alussa. Sen jälkeen on metsien pinta-ala hitaasti kasvanut etupäässä sen vuoksi, että teollinen puunhakkuu luonnonmetsissä kiellettiin kokonaan v. 1998 ja 25 astetta jyrkemmillä rinteillä myös maatalouden harjoittaminen on saatu lopetetuksi viljelijöille tarjottujen taloudellisten houkuttimien avulla. Täten myös pelto- ja laidunmaata on voitu ennallistaa metsäksi. Ping Zhou pystyi jakamaan 5700 metrin korkeuteen saakka kohoavan vuoristoalueen eroosioalttiudeltaan erilaisiin vyöhykkeisiin rinteen kaltevuuden, sademäärän, kasvipeitteen ja maalajin perusteella. Noin 15 prosentilla tutkitun valuma-alueen pinta-alasta, lähinnä Min-joen pääuomaa ympäröivillä jyrkillä rinteillä, eroosioriski oli suuri tai erittäin suuri. Eri tyyppisellä kasvillisuudella oli hyvin erilainen vaikutus eroosioalttiuteen, ja myös alueen sijainti vuoriston eri korkeuksilla vaikutti eroosioon. Säästyneet lähes luonnontilaiset havumetsät, joita on etupäässä vuoriston ylimmissä osissa 2600-4000 metrin korkeudella, edistävät tehokkaasti metsän luontaista uudistumista ja levittäytymistä vaurioituneille alueille. Säilyneiden metsien puulajikoostumus antoi tutkimuksessa mahdollisuuden ennustaa metsien tulevaa kehitystä koko tutkitulla valuma-alueella sen eri korkeusvyöhykkeissä ja eri maaperätyypeillä. Ennallistamisen kannalta ongelmallisimpia olivat alueet joilta metsäpeite oli lähinnä puiden teollisen hakkuun vuoksi kokonaan hävinnyt ja joilla maaperä yleisesti oli eroosion pahoin kuluttama. Näillä alueilla ei ole tehty juuri mitään uudistamis- tai ennallistamistoimenpiteitä. Niillä metsien ennallistaminen vaatii myös puiden tai pensaiden istuttamista. Tähän sopivia ovat erityisesti ilmakehän typpeä sitovat lajit, joista alueella kasvaa luontaisena mm. sama tyrnilaji joka esiintyy myös Suomessa. Työssä tutkittiin yli kahdeksankymmenen paikallisen luontaisen puulajin (joista peräti noin kolmannes on havupuulajeja) ekologisia ominaisuuksia ja soveltuvuutta metsien ennallistamiseen. Avainasemassa työn onnistumisen kannalta ovat nyt paikalliset asukkaat, joiden maankäytön muutokset ovat jo selvästi edistänet luonnonmetsän ennalleen palautumista. Suomen Akatemia rahoitti vuosina 2004-2006 VITRI:n tutkimushanketta, josta Ping Zhou'n väitöskirjatyö muodosti keskeisen osan. Kenttätyö Sichuanissa avasi mahdollisuuden hedelmälliseen monitieteiseen yhteistyöhön ja tutkijavaihtoon Kiinan tiedeakatemian alaisen Chengdun biologiainstituutin (CIB) kanssa; tämä tieteellinen kanssakäyminen jatkuu edelleen.

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.

Disturbance in boreal spruce forest - immediate dynamics from stand to understorey level

The immediate effects of two human-related vegetation disturbances, (1) green tree retention (GTR) patch felling and scarification by harrowing and (2) experimental understorey vegetation layer removal, were examined in boreal forest stands in Finland. Effects of GTR patch felling and scarification on tree uprootings, on coarse woody debris (CWD) and on epixylic plant community were followed in upland and in paludified forest types. Uprootings increased considerably during 2-3 years after the fellings and were more frequent (47%) in the paludified than in the upland forest (13%). Scarification reduced 68% of the CWD in the felling area. Cover and especially species richness of epixylics declined in the both areas during 1-2 years after the felling. The increasing size of GTR patch correlated positively with the species richness. Regeneration of understorey vegetation community and Vaccinium myrtillus and Vaccinium vitis-idaea after different removals of vegetation layers in an old-growth forest took four years. The regeneration occurred mainly by vegetative means and it was faster in the terms of species richness than in the cover. In the most severe treatment, recovery occurred merely by sexual reproduction. V. myrtillus recovered mainly by producing new shoots. V. vitis-idaea recovered faster than V. myrtillus, mainly by increasing length growth. For ecological reasons, use of larger GTR patches on paludified biotope would be recommendable. In felling areas, scarification by harrowing could be replaced with some other spot-wise method. After moderate intensity level disturbance, recovery occurs rapidly by vegetative regrowth of the dominating species. High level of intensity may prevent the recovery of vegetation community for years, while enabling also the genetic regeneration of the initial species. Local anthropogenic-related disturbances are currently increasing and they can interact during temporally short times, which should be taken in to account in the future forest management plans.

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.

Tree mortality and deadwood dynamics in late-successional boreal forests

Here I aimed at quantifying the main components of deadwood dynamics, i.e. tree mortality, deadwood pools, and their decomposition, in late-successional boreal forests. I focused on standing dead trees in three stand types dominated by Picea mariana and Abies balsamea in eastern Canada, and on standing and down dead trees in Picea abies-dominated stands in three areas in Northern Europe. Dead and living trees were measured on five sample plots of 1.6-ha size in each study area and stand type. Stem disks from dead trees were sampled to determine wood density and year of death, using dendrochronological methods. The results were applied to reconstruct past tree mortality and to model deadwood decay class dynamics. Site productivity, stand developmental stage, and the occurrence of episodic tree mortality influenced deadwood volume and quality. In all study areas tree mortality was continuous, leading to continuity in deadwood decay stage distribution. Episodic tree mortality due to either autogenic or allogenic causes influenced deadwood volume and quality in all but one study area. However, regardless of productivity and disturbance history deadwood was abundant, accounting for 20 53% of total wood volume in European study areas, and 15 27% of total standing volume in eastern Canada. Deadwood was a persistent structural component, since its expected residence time in early- and midstages of decay was 18 yr even in the area with the most rapid decomposition. The results indicated that in the absence of episodic tree mortality, stands may eventually develop to a steady state, in which deadwood volume fluctuates around an equilibrium state. However, in many forests deadwood is naturally variable, due to recurrent moderate-severity disturbances. This variability, the continuous tree mortality, and variation in rates of wood decomposition determine the dynamics and availability of deadwood as a habitat and carbon storage medium in boreal coniferous forest ecosystems.

Hybrid Simulation of Queuing-System Dynamics

This paper presents a simple hybrid computer technique to study the transient behaviour of queueing systems. This method is superior to stand-alone analog or digital solution because the hardware requirement is excessive for analog technique whereas computation time is appreciable in the latter case. By using a hybrid computer one can share the analog hardware thus requiring fewer integrators. The digital processor can store the values, play them back at required time instants and change the coefficients of differential equations. By speeding up the integration on the analog computer it is feasible to solve a large number of these equations very fast. Hybrid simulation is even superior to the analytic technique because in the latter case it is difficult to solve time-varying differential equations.

Litter dynamics and phenology of Melaleuca quinquenervia in south Florida

We monitored litterfall biomass at six different sites of melaleuca (Melaleuca quinquenervia (Cav.) S.T. Blake) forested wetlands in South Florida from July 1997 to June 1999. Annual litterfall of melaleuca varied between sites from 6.5 to 9.9 t dry wt ha(-1) yr(1) over the two-year period. Litterfall was significantly higher (p < 0.0001) in scasonally flooded habitats (9.3 t ha(-1) yr(1)) than in non-flooded (7.5 t ha(-1) yr(1)) and permanently flooded habitats (8.0 t ha(-1) yr(1)). Leaf fall was the major component forming 70% of the total litter, woody material 16%, and reproductive material 11%. Phenology of flowering and leaf flush was investigated by examination of the timing and duration of the fall of different plant parts in the litter traps, coupled with monthly field observations during the two-year study. In both years, flowering began in October and November, with peak flowers production around December, and was essentially completed by February and March. New shoot growth began in mid winter after peak flowering, and extended into the spring. Very little new growth was observed in melaleuca forests during the summer months, from May to August, in South Florida. In contrast, the fall of leaves and small wood was recorded in every month of the year, but generally increased during the dry season with higher levels observed from February to April. Also, no seasonality was recorded in the fall of seed capsules, which apparently resulted from the continual self-thinning of small branches and twigs inside the forest stand. In planning management for perennial weeds, it is important to determine the period during its annual growth cycle when the plant is most susceptible to control measures. These phenological data suggest that the appropriate time for melaleuca control in South Florida might be during late winter and early spring, when the plant is most active.

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.