Disjunct eddy covariance measurements of volatile organic compound fluxes using proton transfer reaction mass spectrometry

Volatile organic compounds (VOCs) are emitted into the atmosphere from natural and anthropogenic sources, vegetation being the dominant source on a global scale. Some of these reactive compounds are deemed major contributors or inhibitors to aerosol particle formation and growth, thus making VOC measurements essential for current climate change research. This thesis discusses ecosystem scale VOC fluxes measured above a boreal Scots pine dominated forest in southern Finland. The flux measurements were performed using the micrometeorological disjunct eddy covariance (DEC) method combined with proton transfer reaction mass spectrometry (PTR-MS), which is an online technique for measuring VOC concentrations. The measurement, calibration, and calculation procedures developed in this work proved to be well suited to long-term VOC concentration and flux measurements with PTR-MS. A new averaging approach based on running averaged covariance functions improved the determination of the lag time between wind and concentration measurements, which is a common challenge in DEC when measuring fluxes near the detection limit. The ecosystem scale emissions of methanol, acetaldehyde, and acetone were substantial. These three oxygenated VOCs made up about half of the total emissions, with the rest comprised of monoterpenes. Contrary to the traditional assumption that monoterpene emissions from Scots pine originate mainly as evaporation from specialized storage pools, the DEC measurements indicated a significant contribution from de novo biosynthesis to the ecosystem scale monoterpene emissions. This thesis offers practical guidelines for long-term DEC measurements with PTR-MS. In particular, the new averaging approach to the lag time determination seems useful in the automation of DEC flux calculations. Seasonal variation in the monoterpene biosynthesis and the detailed structure of a revised hybrid algorithm, describing both de novo and pool emissions, should be determined in further studies to improve biological realism in the modelling of monoterpene emissions from Scots pine forests. The increasing number of DEC measurements of oxygenated VOCs will probably enable better estimates of the role of these compounds in plant physiology and tropospheric chemistry. Keywords: disjunct eddy covariance, lag time determination, long-term flux measurements, proton transfer reaction mass spectrometry, Scots pine forests, volatile organic compounds

The effects of abiotic and biotic factors on somatic embryogenesis and seedlings of Pinus sylvestris (L.)

Somatic embryogenesis (SE) is an asexual form of plant propagation that occurs in nature and mimics many of the events of sexual reproduction. Pinus sylvestris (L.) is an important source of timber in Northern Eurasia but it is recalcitrant to somatic embryogenesis. Several factors important for the success of the P. sylvestris embryogenic cultures have not been thoroughly investigated. In this study, we examined the effects of parental genotypes on the SE in P. sylvestris, the involvement of the gaseous plant growth regulator, ethylene in SE, and also biotic effects on somatic embryos as well as on seedlings. We tested parental effects on immature embryo initiation for different media, storage periods, and on the maturation process. Maternal effects were found to be crucial for SE in the absence of paternal effects. No maternal-paternal interaction was observed at any stage of somatic embryo production. Additionally the role of ethylene at different developmental stages of SE was investigated. Two ACC synthase genes, PsACS1 and PsACS2, were isolated and characterized. PsACS1 was expressed during the proliferation stage in all tested genotypes, whereas PsACS2 was only expressed in somatic embryos of each genotype. Ethylene production in embryos at stage 3 was significantly higher than the other stages. In a parallel study, the response of somatic embryos to fungal elicitors was investigated. Three fungi, a mutualistic ectomycorrhizal (ECM) fungus (Suillus bovinus), a weak Scots pine pathogen (Heterobasidion parviporum) and a strong pathogen (H. annosum) were used. The gene expression patterns for embryos exposed to the H. parviporum elicitor were found to be similar to that documented for S. bovinus among the tested genes. By contrast somatic embryos exposed to the H. annosum elicitor had a different pattern of regulation which was marked by a delayed response, and in some cases death of the embryos. Furthermore, interaction without direct contact between P. sylvestris seedlings and microbes (mutualistic and pathogenic fungus, cyanobacterium) were investigated. Several novel genes expressed in seedlings treated with ECM fungus were isolated which suggested that physical contact is not necessary for elicitation of host responses. The results suggest that somatic embryos and seedlings of P. sylvestris are genetically well equipped to respond to fungal elicitor/exudates and could serve as a suitable model for reproducible molecular studies in conifer tree patho- and symbiotic systems.

Environmental and climatic dependences of stable isotope ratios in tree rings on different temporal scales

This work examines stable isotope ratios of carbon, oxygen and hydrogen in annual growth rings of trees. Isotopic composition in wood cellulose is used as a tool to study past climate. The method benefits from the accurate and precise dating provided by dendrochronology. In this study the origin, nature and the strength of climatic correlations are studied on different temporal scales and at different sites in Finland. The origin of carbon isotopic signal is in photosynthetic fractionation. The basic physical and chemical fractionations involved are reasonably well understood. This was confirmed by measuring instantaneous photosynthetic discrimination on Scots pine (Pinus sylvestris L.). The internal conductance of CO2 was recognized to have a significant impact on the observed fractionation, and further investigations are suggested to quantify its role in controlling the isotopic signal of photosynthates. Isotopic composition of the produced biomass can potentially be affected by variety of external factors that induce physiological changes in trees. Response of carbon isotopic signal in tree ring cellulose to changes in resource availability was assessed in a manipulation experiment. It showed that the signal was relatively stable despite of changes in water and nitrogen availability to the tree. Palaeoclimatic reconstructions are typically based on functions describing empirical relationship between isotopic and climatic parameters. These empirical relationships may change depending on the site conditions, species and timeframe studied. Annual variation in Scots pine tree ring carbon and oxygen isotopic composition was studied in northern and in central eastern Finland and annual variation in tree ring latewood carbon, oxygen and hydrogen isotopic ratio in Oak (Quercus robur L.) was studied in southern Finland. In all of the studied sites at least one of the studied isotope ratios was shown to record climate strongly enough to be used in climatic reconstructions. Using the observed relationships, four-century-long climate reconstructions from living Scots pine were created for northern and central eastern Finland. Also temporal stability of the relationships between three proxy indicators, tree ring growth and carbon and oxygen isotopic composition was studied during the four-hundred-year period. Isotope ratios measured from tree rings in Finland were shown to be sensitive indicators of climate. Increasing understanding of environmental controls and physiological mechanisms affecting tree ring isotopic composition will make possible more accurate interpretation of isotope data. This study also demonstrated that by measuring multiple isotopes and physical proxies from the same tree rings, additional information on tree physiology can be obtained. Thus isotopic ratios measured from tree ring cellulose provide means to improve the reliability of climate reconstructions.

Turvemaan männyn (Pinus sylvestris L.) ja kuusen (Picea abies [L.] Karst.) kanto- ja juuripuun ominaisuudet ja asetoniliukoisten uuteaineiden määrä

Tutkimus on osa Metsäklusteri Oy:n Future Biorefinery –tutkimusohjelmaa, jossa kartoitetaan mahdollisuuksia hyödyntää metsäteollisuuden raaka-aineita aiempaa tarkemmin ja uusissa tuotteissa. Tutkimuksen tavoitteena on selvittää männyn (Pinus sylvestris L.) ja kuusen (Picea abies [L.] Karst.) juurten ja kantopuun rakenne ja ominaisuudet. Tutkimuksessa selvitetään löytyykö männyn ja kuusen juurista reaktiopuuta ja määritetään asetoniliukoisten uuteaineiden osuus kantoja juuripuussa. Tutkimusaineistona oli viisi eri-ikäistä mäntyä ja viisi eri-ikäistä kuusta. Juuri- ja kantoaineisto kerättiin Metsäntutkimuslaitoksen toimesta Parkanon seudulta (62.017°N, 23.017°E) hakkuun jälkeen. Maanalaisista juurista otettiin näytteet kolmelta eri etäisyydeltä juurenniskaan nähden. Kummankaan lajin juurista ei löytynyt varsinaista reaktiopuuta, mutta joissakin näytteissä havaittiin lievää reaktiopuuta. Lievää reaktiopuuta löytyi enemmän männyn kuin kuusen juurista, eikä sitä löytynyt lainkaan kaikkein ohuimmista, noin 2 cm paksuisista juurenosista. Männyn kannoissa uuteaineprosentti oli korkeampi kuin kuusen. Männyn juurissa uuteaineprosentti kasvoi edettäessä kohti juuren kärkeä. Kuusella uuteaineprosentti laski aluksi, mutta lähellä juuren kärkeä taas kasvoi. Kuoren uuteainepitoisuus oli molemmilla puulajeilla korkeampi kuin puuaineen. Tutkimusaineisto oli suppea, eikä tutkimuksessa pyritty tilastolliseen yleistettävyyteen. Laajemmasta aineistosta tehdylle tutkimukselle on tarvetta, sillä turvekankailta saatavan puun tarjonta on Suomessa kasvussa, mutta juurten uuteaine- ja reaktiopuututkimuksia on tehty vain kivennäismailta kerätyistä aineistoista.

Climatic adaptation of Norway spruce (Picea abies (L.) Karsten) in Finland based on male flowering phenology.

Anthesis was studied at the canopy level in 10 Norway spruce stands from 9 localities in Finland from 1963 to 1974. Distributions of pollen catches were compared to the normal Gaussian distribution. The basis for the timing studies was the 50 per cent point of the anthesis-fitted normal distribution. Development up to this point was given in calendar days, in degree days (>5 °C) and in period units. The count of each parameter began on March 19 (included). Male flowering in Norway spruce stands was found to have more annual variation in quantity than in Scots pine stands studied earlier. Anthesis in spruce in northern Finland occurred at a later date than in the south. The heat sums needed for anthesis varied latitudinally less in spruce than in pine. The variation of pollen catches in spruce increased towards north-west as in the case of Scots pine. In the unprocessed data, calendar days were found to be the most accurate forecast of anthesis in Norway spruce both for a single year and for the majority of cases of stand averages over several years. Locally, the period unit could be a more accurate parameter for the stand average. However, on a calendar day basis, when annual deviations between expected and measured heat sums were converted to days, period units were narrowly superior to days. The geographical correlations respect to timing of flowering, calculated against distances measured along simulated post-glacial migration routes, were stronger than purely latitudinal correlations. Effects of the reinvasion of Norway spruce into Finland are thus still visible in spruce populations just as they were in Scots pine populations. The proportion of the average annual heat sum needed for spruce anthesis grew rapidly north of a latitude of ca. 63° and the heat sum needed for anthesis decreased only slighty towards the timberline. In light of flowering phenology, it seems probable that the northwesterly third of Finnish Norway spruce populations are incompletely adapted to the prevailing cold climate. A moderate warming of the climate would therefore be beneficial for Norway spruce. This accords roughly with the adaptive situation in Scots pine.

Assessment of insect occurrence in boreal forests based on satellite imagery and field measurements.

The presence/absence data of twenty-seven forest insect taxa (e.g. Retinia resinella, Formica spp., Pissodes spp., several scolytids) and recorded environmental variation were used to investigate the applicability of modelling insect occurrence based on satellite imagery. The sampling was based on 1800 sample plots (25 m by 25 m) placed along the sides of 30 equilateral triangles (side 1 km) in a fragmented forest area (approximately 100 km2) in Evo, S Finland. The triangles were overlaid on land use maps interpreted from satellite images (Landsat TM 30 m multispectral scanner imagery 1991) and digitized geological maps. Insect occurrence was explained using either environmental variables measured in the field or those interpreted from the land use and geological maps. The fit of logistic regression models varied between species, possibly because some species may be associated with the characteristics of single trees while other species with stand characteristics. The occurrence of certain insect species at least, especially those associated with Scots pine, could be relatively accurately assessed indirectly on the basis of satellite imagery and geological maps. Models based on both remotely sensed and geological data better predicted the distribution of forest insects except in the case of Xylechinus pilosus, Dryocoetes sp. and Trypodendron lineatum, where the differences were relatively small in favour of the models based on field measurements. The number of species was related to habitat compartment size and distance from the habitat edge calculated from the land use maps, but logistic regressions suggested that other environmental variables in general masked the effect of these variables in species occurrence at the present scale.

Understorey vegetation as an indicator of forest site potential in Southern Finland.

The relationship between site characteristics and understorey vegetation composition was analysed with quantitative methods, especially from the viewpoint of site quality estimation. Theoretical models were applied to an empirical data set collected from the upland forests of southern Finland comprising 104 sites dominated by Scots pine (Pinus sylvestris L.), and 165 sites dominated by Norway spruce (Picea abies (L.) Karsten). Site index H100 was used as an independent measure of site quality. A new model for the estimation of site quality at sites with a known understorey vegetation composition was introduced. It is based on the application of Bayes' theorem to the density function of site quality within the study area combined with the species-specific presence-absence response curves. The resulting posterior probability density function may be used for calculating an estimate for the site variable. Using this method, a jackknife estimate of site index H100 was calculated separately for pine- and spruce-dominated sites. The results indicated that the cross-validation root mean squared error (RMSEcv) of the estimates improved from 2.98 m down to 2.34 m relative to the "null" model (standard deviation of the sample distribution) in pine-dominated forests. In spruce-dominated forests RMSEcv decreased from 3.94 m down to 3.16 m. In order to assess these results, four other estimation methods based on understorey vegetation composition were applied to the same data set. The results showed that none of the methods was clearly superior to the others. In pine-dominated forests, RMSEcv varied between 2.34 and 2.47 m, and the corresponding range for spruce-dominated forests was from 3.13 to 3.57 m.

Heating values of mature trees.

The effective heating values of the above and below ground biomass components of mature Scots pine (Pinus sylvestris), Norway spruce (Picea abies), downy birch (Betula pubescens), silver birch (Betula pendula), grey alder (Alnus incana), black alder (Alnus glutinosa) and trembling aspen (Populus tremula) were studied. Each sample tree was divided into wood, bark and foliage components. Bomb calorimetry was used to determine the calorimetric heating values. The species is a significant factor in the heating value of individual tree components. The heating value of the wood proper is highest in conifers. Broad-leaved species have a higher heating value of bark than conifers. The species factor diminishes when the weighted heating value of crown, whole stems or stump-root-system are considered. The crown material has a higher heating value per unit weight in comparison with fuelwood from small-sized stems or wholetrees. The additional advantages of coniferous crown material are that it is a non-industrial biomass resource and is readily available. The variability of both the chemical composition and the heating value is small in any given tree component of any species. However, lignin, carbohydrate and extractive content were found to vary from one part of the tree to another and to correlate with the heating value.

Pre-harvest measurement of pine stands for sawing production planning.

To enhance the utilization of the wood, the sawmills are forced to place more emphasis on planning to master the whole production chain from the forest to the end product. One significant obstacle to integrating the forest-sawmill-market production chain is the lack of appropriate information about forest stands. Since the wood procurement point of view in forest planning systems has been almost totally disregarded there has been a great need to develop an easy and efficient pre-harvest measurement method, allowing separate measurement of stands prior to harvesting. The main purpose of this study was to develop a measurement method for pine stands which forest managers could use in describing the properties of the standing trees for sawing production planning. Study materials were collected from ten Scots pine stands (Pinus sylvestris) located in North Häme and South Pohjanmaa, in southern Finland. The data comprise test sawing data on 314 pine stems, dbh and height measures of all trees and measures of the quality parameters of pine sawlog stems in all ten study stands as well as the locations of all trees in six stands. The study was divided into four sub-studies which deal with pine quality prediction, construction of diameter and dead branch height distributions, sampling designs and applying height and crown height models. The final proposal for the pre-harvest measurement method is a synthesis of the individual sub-studies. Quality analysis resulted in choosing dbh, distance from stump height to the first dead branch (dead branch height), crown height and tree height as the most appropriate quality characteristics of Scots pine. Dbh and dead branch height are measured from each pine sample tree while height and crown height are derived from dbh measures by aid of mixed height and crown height models. Pine and spruce diameter distribution as well as dead branch height distribution are most effectively predicted by the kernel function. Roughly 25 sample trees seems to be appropriate in pure pine stands. In mixed stands the number of sample trees needs to be increased in proportion to the intensity of pines in order to attain the same level of accuracy.

Studies of fungi associated with Tomicus piniperda L. in Finland

The Scots pine bark beetle, Tomicus piniperda is a secondary colonizer of pine and other conifers. It is a native species from Europe and Asia that was recently introduced in North America. Although it is necessary to understand this insect's interactions with other organisms, few studies have focussed on its fungal associates. This study focused on the effect of latitude in the occurrence of fungi associated with T. piniperda. T. piniperda were collected from Pinus sylvestris in Northern (Rovaniemi) and Southern (Hyytiala) Finland. Both endo- and epi- mycota were isolated. The fungi were identified using a combination of morphological features and molecular data. The results revealed a great diversity of fungi species associated with T. piniperda, with a total of 3073 isolates representing 23 species. The most frequently isolated fungi in the bark beetles from Northern Finland were Beauvaria bassiana, Kuraishia sp. and Penicillium sp. whereas P. brevicompactum and Mortierella sp. were mostly observed in the South. Ophiostoma canum and O. minus were also observed. The number of isolates per insect in the north was 2.83 epi- and 2.38 for endo-mycota fungus. In the south, the number of isolates per insect was 4.1 for epi- and 3.5 for endo-mycota. Statistical analysis indicated that there was significant differences in fungal populations associated with the beetles in Southern and Northern Finland. There was however no significant difference between the epi- and endo-mycota fungal populations. The highest richness and diversity of the fungal species was observed in the South. However, the overall fungal diversity index analysis revealed that the mycobiota was undersampled.

Tree mortality, fire scar formation and regeneration 8 years after low-intensity fire in managed Pinus sylvestris stands

Fire is an important driver of the boreal forest ecosystem, and a useful tool for the restoration of degraded forests. However, we lack knowledge on the ecological processes initiated by prescribed fires, and whether they bring about the desired restoration effects. The purpose of this study was to investigate the impacts of low-intensity experimental prescribed fires on four ecological processes in young commercial Scots pine (Pinus sylvestris) stands eight years after the burning. The processes of interest were tree mortality, dead wood creation, regeneration and fire scar formation. These were inventoried in twelve study plots, which were 30 m x 30 m in size. The plots belonged to two different stand age classes: 30-35 years or 45 years old at the time of burning. The study was partly a follow-up of study plots researched by Sidoroff et al. (2007) one year after burning in 2003. Tree mortality increased from 183 stems ha-1 in 2003 to 259 stems ha-1 in 2010, corresponding to 15 % and 21 % of stem number respectively. Most mortality was experienced in the stands of the younger age class, in smaller diameter classes and among species other than Scots pine. By 2010, the average mortality of Scots pine per plot was 18%, but varied greatly ranging from 0% to 63% of stem number. Delayed mortality, i.e. mortality that occurred between 2 and 8 years after fire, seemed to become more important with increasing diameter. The input of dead wood also varied greatly between plots, from none to 72 m3 ha-1, averaging at 12 m3 ha-1. The amount of fire scarred trees per plot ranged from none to 20 %. Four out of twelve plots (43 %) did not have any fire scars. Scars were on average small: 95% of scars were less than 4 cm in width, and 75% less than 40 cm in length. Owing to the light nature of the fire, the remaining overstorey and thick organic layer, regeneration was poor overall. The abundance of pine and other seedlings indicated a viable seed source existed, but the seedlings failed to establish under dense canopy. The number of saplings ranged from 0 to 12 333 stems ha-1. The results of this study indicate that a low intensity fire does not necessarily initiate the ecological processes of tree mortality, dead wood creation and regeneration in the desired scale. Fire scars, which form the basis of fire dating in fire history studies, did not form in all cases.

Does the age of fine root carbon indicate the age of fine roots in boreal forests?

To test the reliability of the radiocarbon method for determining root age, we analyzed fine roots (originating from the years 1985 to 1993) from ingrowth cores with known maximum root age (1 to 6 years old). For this purpose, three Scots pine (Pinus sylvestris L.) stands were selected from boreal forests in Finland. We analyzed root 14C age by the radiocarbon method and compared it with the above-mentioned known maximum fine root age. In general, ages determined by the two methods (root 14C age and ingrowth core root maximum age) were in agreement with each other for roots of small diameter (<0.5mm). By contrast, in most of the samples of fine roots of larger diameter (1.5-2mm), the 14C age of root samples of 1987-89 exceeded the ingrowth core root maximum age by 1-10 years. This shows that these roots had received a large amount of older stored carbon from unknown sources in addition to atmospheric CO2 directly from photosynthesis. We conclude that the 14C signature of fine roots, especially those of larger diameter, may not always be indicative of root age, and that further studies are needed concerning the extent of possible root uptake of older carbon and its residence time in roots. Keywords: fine root age, Pinus sylvestris, radiocarbon, root carbon, ingrowth cores, tree ring

Colonization and development of Oribatida mite communities (Acari: Oribatida) in the forest reclaimed limestone mine dumps

In order to study the colonization and development of moss mites (Oribatida) communities in a Scots pine forest of a reclaimed limestone mine dump in Northern Poland, 3 plots from the dump were chosen. The selected plots differed in age, 5 years old, 35 and 50 years old. From a total of 30 samples 499 mites (Acari) were extracted in Tullgren funnel from which 262 were Oribatida. Abundance (N) was analyzed in all mites and after determining the species of both, juvenile and adult stages of oribatids, the following indices were analyzed: Abundance (N), Dominance (D), Species diversity (S), Species richness (s) and Shannon’s diversity index (H). Regarding to the results obtained; oribatid mites were dominant with the highest abundance in all assemblages (Plot 1: 139 Oribatida /299 Acari. Plot 2: 40/55 and Plot 3: 83/145). Tectocepheus velatus showed a very high dominance (45,99%) in plot 1; the highest value for Shannon’s diversity index belonged to plot 3. On the other hand, juvenile’s percentage was significantly higher than adult’s percentage, especially at plot 2 (95,02%). These results made us to conclude that the high abundance of oribatids in the youngest forest is due to T. velatus’s high abundance and that plot 3 is the best habitat for mites. Finally, the high occurrence of juvenile stages requires keeping on studying the area.

Rainfall infiltration and soil hydrological characteristics below ancient forest, planted forest, and grassland in a temperate northern climate

How rainfall infiltration rate and soil hydrological characteristics develop over time under forests of different ages in temperate regions is poorly understood. In this study, infiltration rate and soil hydrological characteristics were investigated under forests of different ages and under grassland. Soil hydraulic characteristics were measured at different scales under a 250 year old grazed grassland (GL), a six (6 yr) and 48 (48 yr) year old Scots pine (Pinus sylvestris) plantation, remnant 300 year old individual Scots pines (OT) and a 4000 year old Caledonian Forest (AF). In-situ field saturated hydraulic conductivity (Kfs) was measured and visible root:soil area was estimated from soil pits. Macroporosity, pore structure, and macropore connectivity were estimated from X-ray tomography of soil cores, and from water-release characteristics. At all scales the median values for Kfs, root fraction, macro-porosity and connectivity values tended to AF > OT > 48 yr > GL > 6 yr, indicating that infiltration rates and water storage increased with forest age. The remnant Caledonian Forest had a huge range of Kfs (12 to > 4922 mm h-1), with maximum Kfs values 7 to 15 times larger than 48-year-old Scots pine plantation, suggesting that undisturbed old forests, with high rainfall and minimal evapotranspiration in winter, may act as important areas for water storage and sinks for storm rainfall to infiltrate and transport to deeper soil layers via preferential flow. The importance of the development of soil hydrological characteristics under different aged forests is discussed.

Inter-annual variability of precipitation constrains the production response of boreal Pinus sylvestris to nitrogen fertilization

© 2015 Published by Elsevier B.V.Tree growth resources and the efficiency of resource-use for biomass production determine the productivity of forest ecosystems. In nutrient-limited forests, nitrogen (N)-fertilization increases foliage [N], which may increase photosynthetic rates, leaf area index (L), and thus light interception (IC). The product of such changes is a higher gross primary production and higher net primary production (NPP). However, fertilization may also alter carbohydrate partitioning from below- to aboveground, increasing aboveground NPP (ANPP). We analyzed effects of long-term N-fertilization on NPP, and that of long-term carbon storing organs (NPPS) in a Pinus sylvestris forest on sandy soil, a wide-ranging forest type in the boreal region. We based our analyses on a combination of destructive harvesting, consecutive mensuration, and optical measurements of canopy openness. After eight-year fertilization with a total of 70gNm^-2, ANPP was 27±7% higher in the fertilized (F) relative to the reference (R) stand, but although L increased relative to its pre-fertilization values, IC was not greater than in R. On the seventh year after the treatment initiation, the increase of ANPP was matched by the decrease of belowground NPP (78 vs. 92gCm^-2yr^-1; ~17% of NPP) and, given the similarity of IC, suggests that the main effect of N-fertilization was changed carbon partitioning rather than increased canopy photosynthesis. Annual NPPS increased linearly with growing season temperature (T) in both treatments, with an upward shift of 70.2gCm^-2yr^-1 by fertilization, which also caused greater amount of unexplained variation (r²=0.53 in R, 0.21 in F). Residuals of the NPPS-T relationship of F were related to growing season precipitation (P, r²=0.48), indicating that T constrains productivity at this site regardless of fertility, while P is important in determining productivity where N-limitation is alleviated. We estimated that, in a growing season average T (11.5±1.0°C; 33-year-mean), NPPS response to N-fertilization will be nullified with P 31mm less than the mean (325±85mm), and would double with P 109mm greater than the mean. These results suggest that inter-annual variation in climate, particularly in P, may help explaining the reported large variability in growth responses to fertilization of pine stands on sandy soils. Furthermore, forest management of long-rotation systems, such as those of boreal and northern temperate forests, must consider the efficiency of fertilization in terms of wood production in the context of changes in climate predicted for the region.