The structure of Norway spruce (Picea abies [L.] Karst.) stems in relation to wood properties of sawn timber

An important challenge in forest industry is to get the appropriate raw material out from the forests to the wood processing industry. Growth and stem reconstruction simulators are therefore increasingly integrated in industrial conversion simulators, for linking the properties of wooden products to the three-dimensional structure of stems and their growing conditions. Static simulators predict the wood properties from stem dimensions at the end of a growth simulation period, whereas in dynamic approaches, the structural components, e.g. branches, are incremented along with the growth processes. The dynamic approach can be applied to stem reconstruction by predicting the three-dimensional stem structure from external tree variables (i.e. age, height) as a result of growth to the current state. In this study, a dynamic growth simulator, PipeQual, and a stem reconstruction simulator, RetroSTEM, are adapted to Norway spruce (Picea abies [L.] Karst.) to predict the three-dimensional structure of stems (tapers, branchiness, wood basic density) over time such that both simulators can be integrated in a sawing simulator. The parameterisation of the PipeQual and RetroSTEM simulators for Norway spruce relied on the theoretically based description of tree structure developing in the growth process and following certain conservative structural regularities while allowing for plasticity in the crown development. The crown expressed both regularity and plasticity in its development, as the vertical foliage density peaked regularly at about 5 m from the stem apex, varying below that with tree age and dominance position (Study I). Conservative stem structure was characterized in terms of (1) the pipe ratios between foliage mass and branch and stem cross-sectional areas at crown base, (2) the allometric relationship between foliage mass and crown length, (3) mean branch length relative to crown length and (4) form coefficients in branches and stem (Study II). The pipe ratio between branch and stem cross-sectional area at crown base, and mean branch length relative to the crown length may differ in trees before and after canopy closure, but the variation should be further analysed in stands of different ages and densities with varying site fertilities and climates. The predictions of the PipeQual and RetroSTEM simulators were evaluated by comparing the simulated values to measured ones (Study III, IV). Both simulators predicted stem taper and branch diameter at the individual tree level with a small bias. RetroSTEM predictions of wood density were accurate. For focusing on even more accurate predictions of stem diameters and branchiness along the stem, both simulators should be further improved by revising the following aspects in the simulators: the relationship between foliage and stem sapwood area in the upper stem, the error source in branch sizes, the crown base development and the height growth models in RetroSTEM. In Study V, the RetroSTEM simulator was integrated in the InnoSIM sawing simulator, and according to the pilot simulations, this turned out to be an efficient tool for readily producing stand scale information about stem sizes and structure when approximating the available assortments of wood products.

Improvement of mango tree architecture

This paper discusses the role of mango canopy architecture in mango productivity and orchard management and considers potential increases on production of high density orchards through improved canopy architecture. Lower tree height, reduced vigour and smaller more open canopies are recognised as important aspects of high density orchards. However, vigour management, light relations, flowering and crop load of high density orchards needs to be better understood if we are to developed sustainable highly productive canopy training and pruning systems that are easy to maintain at high planting densities.

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.

Productive coexistence and gain in agroforestry systems.

Productive coexistence and coexistence gain of populations were studied using nine years' data from field experiments of Taxodium ascendens-intercrop systems in Lixiahe, Jiangsu Province, China. A theoretical framework for productive coexistence in agroforestry was developed. Interaction patterns between trees and intercrops were presented within this framework. A model framework was developed to describe the coexistence gain and interaction of populations in T. ascendens-intercrop systems. Facilitation and resource sharing were identified as main contribution to the advantage of species combination in agroforestry. The model of population interaction developed in the present study was accepted for describing the interaction of populations in T. ascendens-intercrop systems, because it explained a high proportion of the variance of experimental data and fitted well the observations in most intercropping types. The model developed in the present study provides flexibility for describing different patterns of intra- and inter-specific interactions. Model coefficients were applied to the determination of the ecological compatibility of species. Managed T. ascendens-intercrop systems were advantageous as compared to a monoculture of trees or arable crops. In T. ascendens stands up to the age of three, arable crops contributed about 50-80 % of the total biomass yield of agroforestry. The diameter or height growth of T. ascendens was not significantly influenced by intercrops, indicating that intercropping under trees produced extra yields but did not depress the tree growth. When the trees were young (during the first three years), T. ascendens did not depress the crop yields, and a land equivalent ratio greater than unity was obtained together with a high yield of both components. The diameter and height of the trees were similar in four spacing configurations with an equal number of trees per hectare up to the age of eight, but wider between-rows open range were beneficial for the intercrops. The relationship between open-ranges and species coexistence was also analysed and the distribution of soil nutrients studied.

川西亚高山针叶林主要针叶树种年轮生态学研究

米亚罗地区是四川西部较为典型的亚高山针叶林区域之一。为建立该地区主要针叶树种岷江冷杉、云杉、紫果云杉和红杉的年轮宽度年表资料，了解不同海拔高度岷江冷杉原始林和不同恢复过程的人工针叶林及次生混交林树木径向生长规律，结合样地调查，用生长锥钻取了树木芯样做年轮生态学分析。芯样经过标准化程序固定和打磨抛光后，用WinDENDRO图像分析系统测量年轮宽度序列，用COFECHA程序交叉定年和控制测量数据质量，用ARSTAN程序建立了4个主要针叶树种的地区年表和不同海拔高度岷江冷杉林及人工针叶林和次生混交林针叶树的样地年表。 4个主要针叶树种年轮宽度年表的平均敏感度低于0.2，而其晚材宽度年表都具相对较高的平均敏感度。早材宽度与年轮总宽度标准化年表间的相关系数均在0.9以上；晚材宽度与年轮总宽度标准化年表间的相关系数则种间差异较大，红杉的最高，岷江冷杉的最低。岷江冷杉晚材宽度与年轮总宽度的相关性从1970年以后明显下降，而其他种的相关系数则随时间变化较小。树种之间标准化年表显著正相关，而云杉与紫果云杉和红杉与岷江冷杉之间相关系数明显较高。年表序列的第1主分量表达了4个树种树木共同径向生长变化格局；第2至第4主分量分别表达了云杉属和冷杉属、常绿针叶树种和落叶针叶树种以及云杉和紫果云杉树木径向生长变化差异。 不同海拔高度的8个岷江冷杉样地年轮宽度年表序列敏感度大体上随海拔高度升高而降低。各样地早材宽度与年轮总宽度年表之间的相关系数均在0.9以上，且随海拔高度变化不大；晚材宽度与年轮总宽度之间的相关系数随海拔高度的变化较大，并有随海拔升高而降低的趋势。样地年表序列之间相关系数差异很大，高海拔样地年表间多为显著正相关；低海拔样地年表间的相关系数变化不一；高海拔和低海拔样地年表之间相关性较差，且多不显著。样地年表的第1主分量能解释年表序列总方差的37.5%，反映了不同海拔高度岷江冷杉林木共同的径向生长变化格局；第2和第3主分量分别解释总方差的24.5%和18.2%，表现出明显的高海拔和低海拔样地树木间不同的径向生长变化，除一些样地例外，它们一般与低海拔样地年表有正相关，与高海拔样地年表有负相关。在那些另外的样地，海拔以外的其他因素可能也影响了树木径向生长变化。不同海拔高度样地林木的生长抑制和生长释放频率在不同时期表现出较大的差异，表明了不同的干扰历史和林木补充时间。 人工针叶林和次生混交林各样地林木早材宽度与其年轮总宽度年表之间相关系数均高达0.9以上；晚材宽度与年轮总宽度年表之间也都显著正相关，但人工针叶林样地的明显较高。样地年表序列之间的相关关系表现为，林分起源和经营管理相似的样地年表之间的相关系数明显较高，如人工针叶林与人工针叶林尽管树种不同，但样地年表之间显著正相关，而与次生混交林样地年表间关系不显著；反之亦然。综合比较各项生长参数及不同时期的树木径向生长速率，人工针叶林树木的胸径增长至少在40年以内是优于次生混交林的同种（或不同种）针叶树的。不同样地林木生长释放和生长抑制及人工针叶林树木胸高断面积增长分析表明，除严重的人为干扰外，林分郁闭后林木密度过大是造成高频率生长抑制的主要原因，在林分发育的适当时期通过抚育间伐等措施调控林分密度，是保证林木胸高断面积在一定时期内保持较高的连年增长的关键。日本落叶松作为引进的树种，在海拔3100 m左右种植表现良好，近30年来各项生长指标均高于林龄相近的云杉人工林，因此，适当用其作为川西亚高山针叶林采伐迹地快速恢复是合理的。 Miyaluo area is one of the typical regions covered by subalpine coniferous forests in western Sichuan province of southwestern China. To develop the regional tree-ring width chronology series for the dominant conifers such as Abies faxoniana, Picea asperata, P. purpurea and Larix potaninii, and to understand the radial growth patterns of conifers in Abies faxoniana natural forest stands at different altitudes, and in coniferous plantations and natural regenerated mixed stands in their different restoring processes as well, increment cores were sampled in the field together with conventional plots investigations for dendroecological analyses. After the increment cores being prepared according to standard procedures, the ring widths (total-ring and intra-ring widths) were measured with a WinDENDRO image-analysis system, and the measured tree-ring sequences were crossdated and quality-controlled with the software COFECHA. Using the software ARSTAN, we developed tree-ring width based chronology series of the four dominant conifers, eight site-specific Abies faxoniana chronologies, and seven site-specific chronologies of conifers in coniferous plantations and natural regenerated mixed stands. Mean sensitivities for total ring width chronologies of the four sampled dominant conifers were all below 0.2, while those for the latewood width chronologies of the same species were relatively much higher. Correlation coefficients between standard earlywood and total ring width chronologies of the four conifers were all above 0.9, but those between standard latewood and total ring width chronologies exhibited differences among species, with the coefficient of Larix potaninii the highest and that of Abies faxoniana the lowest. Correlation coefficients between latewood and total ring width of A. faxoniana obviously decreased from 1920-1970 for successive 50-year segments with 10-years lag analyses, though the same for the other three species changed unnoticeably with time. Tree-ring standard chronologies among species showed significant positive correlations, with the correlation coefficients between chronologies of Picea asperata and P. purpurea, and of Larix potaninii and Abies faxoniana relatively much higher. The first principal component of tree-ring chronologies represented the common radial growth patterns of the four conifers in Miyaluo area. The second, third and fourth PCs expressed the differences in radial growth responses for the genus Picea and Abies, for the evergreen and deciduous confers, and for the two species of the genus Picea, respectively. In general, mean sensitivities of the eight Abies faxoniana site-specific tree-ring width chronologies decreased with increasing altitude. The correlation coefficients between earlywood and total ring width chronologies for all sites reached 0.9, which did not change much with altitude; but those between latewood and total ring width chronologies diversified, with a decreasing tendency from lower altitudinal sites to higher altitudinal sites. Correlation coefficients among site chronologies varied considerably, with significant positive correlations among higher site chronologies, mixed results among lower site chronologies, and poor and insignificant correlations between chronologies of higher site and lower site. The first PC, which represents 37.5% of the total variance, reflected a common radial growth response at sites of different altitudes, and it showed a tendency of explaining more variance with increasing altitude. The second and the third PCs contributed to 24.5% and 18.2% of the total variance, respectively, exhibiting distinctive differences in radial growth responses between low- and high-altitudinal sites. With some exceptions, the radial growth represented by the second and third PCs had a positive correlation with that at the low-altitudinal sites and a negative correlation with that at the high-altitudinal sites. For those exceptional sites, factors other than altitude might also play a role in affecting tree-ring growth variations. Trees in stands of different altitudes showed great differences in frequencies of growth suppressions and releases through times, suggesting different disturbance histories and periods when trees recruiting to the canopy. Correlation coefficients between earlywood and total ring width chronologies for all sites of coniferous plantations and natural regenerated mixed stands were also above 0.9; and the same between latewood width and total ring width chronologies all positively correlated, too, with the coefficients of the coniferous plantations obviously much higher. Correlations among site chronologies showed that the coefficients among sites with similar stand origin and management regimes were much higher than those among sites with different stand origin and management regimes. For example, significant positive correlations were found for chronologies among different coniferous plantations, irrespective of species differences; while insignificant correlations between chronologies of the same conifer from a coniferous plantation and a natural regenerated mixed stand, and vise versa. Integrative comparisons of different tree growth parameters and radial growth rates at different stages indicated that the diameter at breast height (DBH) increments for trees in coniferous plantations were faster than those for trees of the same (or different) species in the natural regenerated mixed stands, at least within their first 40 years of stand development. Analyses of growth releases and suppressions, and basal area increments of trees in different stands demonstrated that over-dense individuals after canopy closure was the main factor resulting in high frequencies of radial growth suppressions, with some exceptions of severe man-made disturbances. Therefore, to ensure a continuous basal area current annual increment in certain periods, tree density controlling through thinning in due time during the stand development process are necessary. It should be mentioned that, as an introduced conifer to Miyaluo area, Larix kaempferi grew quite well at altitude of ca. 3100 m after planting in 1970s. In their near 30 years of stand development, Larix kaempferi trees exhibited faster growth in various parameters than Picea asperata trees of the similar stand age did. Thus we think it reasonable to use Larix kaempferi as a fast restoring species at appropriate sites of cutting blanks of subalpine coniferous forests in western Sichuan.

Effets de facteurs pédoclimatiques sur la composition et l’anatomie du bois de cinq cultivars de saule destinés aux biocarburants

En 2011, cinq (5) cultivars de saules ont été sélectionnés pour leur rendement en biomasse. Ils ont été plantés sur quatre sites de la province du Québec et ont été maintenus selon le protocole de la culture intensive sur courtes rotations (CICR) afin de déterminer leur potentiel pour la bioénergie dans des environnements contrastés. La composition et l’anatomie du bois de ces cultivars ont été caractérisées et comparés en fonction des conditions environnementales caractéristiques de chaque site. La hauteur et le diamètre à la base des plantes diffèrent selon les sites. Ainsi, les cultivars répondent de façon spécifique aux conditions pédoclimatiques dans lesquelles ils sont cultivés. L’effet de l’environnement n’a pas été mis en évidence sur la teneur en lignine des cultivars. Cependant, un effet génotypique a pu être constaté soulignant l’importance de la sélectivité des cultivars. La densité du bois a étonnamment conservé la même hiérarchie génotypique entre les sites. À l’opposé, l’anatomie du bois présente des différences notamment au niveau des caractéristiques des fibres et des vaisseaux. Une forte teneur en polyphénols ainsi que des fibres moins larges et des vaisseaux plus nombreux ont été observés sur le site dont le bois est le plus dense supposant l’effet probable d’un stress abiotique. De plus, deux fois plus de fibres gélatineuses, fibres riches en cellulose, ont été identifiées sur ce site montrant un intérêt pour la production de bioéthanol.

Replacing Norway spruce with European beech: a comparison of biomass and Net Primary Production patterns in young stands

An alteration of species composition in temperate forests – both managed and natural - is one of the expected effects of environmental change. Present forest tree species ranges will be altered by changing environmental conditions. By a combination of continuous and destructive sampling, we compared biomass stocks and annual NPP in naturally regenerated stands of Norway spruce and European beech. We purposely selected a site where future environmental conditions are predicted to favour beech over presently dominant spruce. We found no difference in overall productivity, but biomass allocation differed significantly between the two species. Beech allocated more assimilates to stem and roots than spruce. There was no significant difference between the species in NPP of the fast turnover biomass pool comprising foliage and fine roots. Maximum height growth occurred about a month earlier than in spruce, potentially changing the timing of carbon (C) flow into the soil pools. We show that the replacement of spruce by beech will result in changes in forest biomass allocation and in alterations of belowground C cycle. Such changes will affect forest ecosystem function by modifying the magnitude and timing of certain C fluxes, but also by potentially changing the species composition of forest biota dependent on them.

THE ENDOPARASITE PILOSTYLES ULEI (APODANTHACEAE - CUCURBITALES) INFLUENCES WOOD STRUCTURE IN THREE HOST SPECIES OF MIMOSA

Pilostyles species (Apodanthaceae) are endoparasites in stems of the plant family Fabaceae. The body comprises masses of parenchyma in the host bark and cortex, with sinkers, comprising groups of twisted tracheal elements surrounded by parenchyma that enter the secondary xylem of the host plant. Here we report for the first time the effects of Pilostyles parasitism on host secondary xylem. We obtained healthy and parasitized stems from Mimosa foliolosa, M. maguirei and M. setosa and compared vessel element length, fiber length, vessel diameter and vessel frequency, measured through digital imaging. Also, tree height and girth were compared between healthy and parasitized M. setosa. When parasitized, plant size, vessel diameter, vessel element length and fiber length are all less than in healthy plants. Also, vessel frequency is greater and vessels are narrower in parasitized stems. These responses to parasitism are similar to those observed in stressed plants. Thus, hosts respond to the parasite by changing its wood micromorphology in favour of increased hydraulic safety.

Evaluation of Trichospilus diatraeae (Hymenoptera: Eulophidae) as parasitoid of the eucalyptus defoliator Euselasia eucerus (Lepidoptera: Riodinidae).

Lepidopteran pests have occurred in eucalyptus plantations in Brazil since 1948, reaching high population levels, reducing tree growth and causing considerable losses in wood production. The control of defoliating caterpillars in eucalyptus forests is complex, mainly due to the large extent of forest plantations and tree height, necessitating the aerial application of chemical or biological insecticides. Due to this complexity, alternative control methods have been proposed, for instance biological control through the use of parasitoids. Trichospilus diatraeae Cherian & Margabandhu, 1942 is a gregarious pupal parasitoid that preferentially attacks species of the order Lepidoptera. This is the first report of T. diatraeae parasitizing pupae of the eucalyptus defoliator Euselasia eucerus Hewitson in Brazil. This parasitoid offers new perspectives for the biological control programmes of this species in eucalyptus plantations in Brazil.

Growth and carbon isotopes of Mediterranean trees reveal contrasting responses to increased carbon dioxide and drought

Forest dynamics will depend upon the physiological performance of individual tree species under more stressful conditions caused by climate change. In order to compare the idiosyncratic responses of Mediterranean tree species (Quercus faginea, Pinus nigra, Juniperus thurifera) coexisting in forests of central Spain, we evaluated the temporal changes in secondary growth (basal area increment; BAI) and intrinsic water-use efficiency (iWUE) during the last four decades, determined how coexisting species are responding to increases in atmospheric CO2 concentrations (Ca) and drought stress, and assessed the relationship among iWUE and growth during climatically contrasting years. All species increased their iWUE (ca. +15 to +21 %) between the 1970s and the 2000s. This increase was positively related to Ca for J. thurifera and to higher Ca and drought for Q. faginea and P. nigra. During climatically favourable years the study species either increased or maintained their growth at rising iWUE, suggesting a higher CO2 uptake. However, during unfavourable climatic years Q. faginea and especially P. nigra showed sharp declines in growth at enhanced iWUE, likely caused by a reduced stomatal conductance to save water under stressful dry conditions. In contrast, J. thurifera showed enhanced growth also during unfavourable years at increased iWUE, denoting a beneficial effect of Ca even under climatically harsh conditions. Our results reveal significant inter-specific differences in growth driven by alternative physiological responses to increasing drought stress. Thus, forest composition in the Mediterranean region might be altered due to contrasting capacities of coexisting tree species to withstand increasingly stressful conditions. © 2013 Springer-Verlag Berlin Heidelberg.

Modelagem do crescimento e produção em povoamentos desbastados de Eucalyptus grandis

The aim of the study was to develop a system of growth and yield models for thinned stands of Eucalyptus spp.; and to assess the behavior of the growth in scenarios with 10% decrease or increase in rainfall. The probability distribution functions Weibull 2 and 3 parameters and Johnson SB for different methods were fitted. Correlation between the fitted parameters with age was evaluated. Dominant height growth behavior was evaluated to check if thinned stand changes its growth when compared to a non-thinned stands. The stand variables dominant height and basal area were projected and simultaneously predicted and projected, respectively. Individual tree equations were fitted, which were fitted as functions of stand level variables in order to decrease the error propagation. R software was used to fit all the proposed models and consequently all the fitted models were evaluated by their parameters significance (F-test) and graphs of predicted values in relation to the observed values around the 1:1 line. Thus, the prognosis system was made by two ways, first one using the full data set, and for the second one the dataset was restricted at age 7.5. Increase and decrease in 20% of rainfall were assessed by updating the site index function. Method of moments was the most precise to describe the diameter distribution for every age in eucalyptus stands for Johnson SB and Weibull 2 parameters pdfs. When observed for each pdf the correlation for their fitted parameters with age, we noticed that shape parameters for a thinned stand were no longer correlated with age, differently of non-thinned stands. Thus, thinning effect was accounted in the basal area prediction and projection modeling. This result emphasized the necessity of applying the Parameter Recovery method in order to assess differences and capture the right pattern for thinned and non-thinned stands in the future. Dominant height was not influenced by thinning intensity. Therefore the fitted Chapman-Richards model did not account for a stand being thinned or not. All the fitted equations behaved with good precision, no matter using full or precocious dataset. The prognosis system using full and/or precocious date set was evaluated for when using Parameter Recovery method for Sb and Weibull pdfs, and by then, graphical analysis and precision statistics showed appropriated results. Finally, the increase or decrease in rainfall regime were observed for eucalyptus stand yields and we may notice how important is to observe this effect, since the growth pattern is strictly affected by water.

Measured and modeled interactive effects of potassium deficiency and water deficit on gross primary productivity and light-use efficiency in Eucalyptus grandis plantations

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Allometric equations for estimating tree biomass in restored mixed-species Atlantic Forest stands

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Wood and leaf anatomy of Copaifera langsdorffii dwarf trees

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Response of nutrient cycles of an old-growth montane forest in Ecuador to experimental low-level nutrient amendments

Atmospheric nitrogen (N) and phosphorus (P) depositions are expected to increase in the tropicsrnas a consequence of increasing human activities in the next decades. Furthermore, a possiblernshortened El Niño Southern Oscillation cycle might come along with more frequent calcium (Ca)rndepositions on the eastern slope of the Ecuadorian Andes originating from Saharan dust. It isrncrucial to understand the response of the old-growth montane forest in Ecuador to increasedrnnutrient deposition to predict the further development of this megadiverse ecosystem.rnI studied experimental additions of N, P, N+P and Ca to the forest and an untreatedrncontrol, all in a fourfold replicated randomized block design. These experiments were conductedrnin the framework of a collaborative research effort, the NUtrient Manipulation EXperimentrn(NUMEX). I collected litter leachate, mineral soil solution (0.15 and 0.30 m depths), throughfallrnand fine litterfall samples and determined N, P and Ca concentrations and fluxes. This approachrnalso allowed me to assess whether N, P and/or Ca are limiting nutrients for forest growth.rnFurthermore, I evaluated the response of fine root biomass, leaf area index, leaf area and specificrnleaf area, tree diameter growth and basal area increment contributed from a cooperating group inrnthe Ca applied and control treatments.rnDuring the observation period of 16 months after the first fertilizer application, less thanrn10, 1 and 5% of the applied N, P and Ca, respectively, leached below the organic layer whichrncontained almost all roots but no significant leaching losses occurred to the deeper mineral soil.rnDeposited N, P and Ca from the atmosphere in dry and wet form were, on balance, retained in therncanopy in the control treatment. Retention of N, P and Ca in the canopy in their respectiverntreatments was reduced resulting in higher concentrations and fluxes of N, P and Ca inrnthroughfall and litterfall. Up to 2.5% of the applied N and 2% of the applied P and Ca werernrecycled to the soil with throughfall. Fluxes of N, P and Ca in throughfall+litterfall were higher inrnthe fertilized treatments than in the control; up to 20, 5 and 25% of the applied N, P and Ca,rnrespectively, were recycled to the soil with throughfall+litterfall.rnIn the Ca-applied plots, fine root biomass decreased significantly. Also the leaf area of thernfour most common tree species tended to decrease and the specific leaf area increasedrnsignificantly in Graffenrieda emarginata Triana, the most common tree species in the study area.rnThese changes are known plant responses to reduced nutrient stress. Reduced aluminium (Al)rntoxicity as an explanation of the Ca effect was unlikely, because of almost complete organocomplexationrnof Al and molar Ca:Al concentration ratios in solution above the toxicity threshold.rnThe results suggest that N, P and Ca co-limit the forest ecosystem functioning in thernnorthern Andean montane forests in line with recent assumptions in which different ecosystemrncompartments and even different phenological stages may show different nutrient limitationsrn(Kaspari et al. 2008). I conclude that (1) the expected elevated N and P deposition will bernretained in the ecosystem, at least in the short term and hence, quality of river water will not bernendangered and (2) increased Ca input will reduce nutrient stress of the forest.