Predicting mature cork biomass with t years of growth from one measurement taken at any other age

Cork is a natural product that is extracted from the outer bark of the cork oak tree. According to Portuguese legislation, the interval between two consecutive cork extractions on the same tree must be equal or greater than 9 years. Although the majority of the cork oak stands are debarked at the end of this period, this rotation may not be the optimum in many cases. The existing models for cork weight prediction can only be used for trees debarked at a 9 years or, in one model, at a 10-years rotation period, since the data used for its development was limited to these growth periods. The development of a method that allows for the prediction of the mature cork biomass with t years of growth, based in one measurement taken at any other age, was the main objective of this work. The method is based on the knowledge that the density of the cork tissue is constant between the inner and outer cork rings, being significantly different from density of the cork back. It can be implemented using two different equations that were developed during this work: a model to estimate cork biomass with 9 years of age and a model to estimate the cork back weight proportion at 9 years of age. For the first model, four different alternative models were developed, considering different levels of information collected during forest inventory. The model to estimate the cork back weight proportion leads to the biomass of cork tissue. Cork biomass at t years is obtained by decreasing or increasing the biomass of cork tissue according to the difference in cork thickness between t and 9 years of growth. The proposed method was evaluated by comparing the observed and estimated values of cork biomass from an independent data set with corks with 9, 10 and 11 years of age. The results showed similar precision for corks with 9, 10 or 11 years of age. As expected the precision of the predictions increases when the model to estimate cork biomass with 9 years of age uses more information. The presented method should be an important tool for cork oak stand management, for the prediction of the evolution of carbon stocks in cork oak stands, and will allow analyzing the impact in cork biomass production of decreasing or extending the interval between two consecutive cork extractions.     

A spatial distribution model of cork oak (Quercus suber) in southwestern Spain: A suitable tool for reforestation

Cork oak (Quercus suber) is an evergreen tree characterized by a thick bark, which grows in Mediterranean schlerophyllous forests. It is most prevalent in the southwestern Iberian Peninsula. Despite the potential of the province of Huelva (southwestern Spain) to maintain mature forests of cork oak, the tree has been severely depleted and most forests have either disappeared or are seriously threatened. This paper presents a spatial distribution model of cork oak for the province of Huelva with a view to determining the optimal areas for reforestation. The model draws on all available digital cartographic information with respect to cork oak distribution: topographic data (altitude, slope and orientation) were obtained from a Digital Terrain Model (20 m scale); rainfall, temperature and PET models were based on data collected from a network of meteorological stations; litologic data derive from the litologic map of Huelva (1:100,000). The result of this work is a mesh of points at a resolution of 100 m, sufficient to meet the needs of any kind of reforestation or management programmes in the area studied. Each point of this mesh contains the corresponding values for bioclimatic, topographic and litologic variables in a georeferenced data matrix. The independent variables responsible for cork oak distribution (binary dependent variables) were then identified by means of binary logistic regression analysis. North-facing slopes, abundant annual rainfall and litology were the main explaining variables. The spatial distribution model was produced by applying the formula obtained to spatial analysis software. This model is proposed as a basis for future reforestation plans, especially in those areas most affected by forest fires.     

Effects of treeshelters and polyethylene mulch sheets on survival and growth of cork oak (Quercus suber L.) seedlings planted in northwestern Tunisia

Success of plantation establishment depends on the ability of seedlings to withstand animal browsing and competition from vegetation. This study in northwestern Tunisia examined the influence of browsing protection and vegetation control on seedling mortality, diameter, and height growth components of an evergreen sclerophyllous oak species, cork oak (Quercus suber L.) for the first 2 years after planting (2005 and 2006). Browsing protection was controlled by varying both treeshelter type (non-vented, vented, and control) and height (1.2-m and 1.8-m tall). Competing vegetation was controlled through increasing sizes of square-shaped black polyethylene mulch sheets: 0 m, 0.7 m, 1 m and 1.3 m large. Seedling mortality rate was relatively low during the first year (9%) and was not affected by any of the studied factors. During the following year, mortality was lower for seedlings in both types (non-vented and vented) of 1.2-m tall treeshelters than for unsheltered ones. Mortality was only slightly higher for 1.8-m tall shelters than for 1.2-m shelters. Compared to unsheltered seedlings, basal diameter was reduced during both seasons for seedlings inside non-vented treeshelters (1.2-m and 1.8-m tall), whereas for vented shelters it was reduced during the first year and then increased during the second year. Height growth was positively affected by tree shelters during both years. Seedlings inside non-vented 1.8-m tall shelters were the tallest having an average height of up to 140 cm at the end of the second year, while the unsheltered seedlings had only an average height of up to 42 cm. The greater height of the sheltered seedlings was attributed to the increase of: (1) the number of shoot flushes yearly established on the main stem (up to three), and (2) the length of all the growth units produced during each flush. Seedlings inside vented shelters had balanced growth with a low proportion of seedlings unable to support themselves without a stake (7%). However, non-vented shelters had unbalanced height and diameter growths resulting in both a higher height-to-diameter ratio and in a higher proportion of seedlings having stability problems (47%). Polyethylene mulch sheets did not affect early survival, seedling basal diameter, or number of shoot flushes established yearly, but did slightly improve seedling height at the end of the second season, mostly due to the significant improvement in shoot elongation of the first flush. Results suggest that both 1.8 m treeshelters and plastic mulches may enhance the growth of cork oak seedlings planted on harsh sites in northwestern Tunisia.     

Cork oak (<Emphasis Type="Italic">Quercus suber</Emphasis> L.) wood growth and vessel characteristics variations in relation to climate and cork harvesting

Variations in tree ring growth of Quercus suber L. were analysed using dendrochronological techniques on cork oak discs from trees harvested in the cork producing region of Alentejo, Portugal. A tree-ring chronology containing a strong common signal and covering the period from 1970 to 1995 was build for ca. 30-year-old cork oaks never submitted to cork harvesting using 14 trees that crossdated satisfactorily out of 30 sampled trees. The tree ring indices correlated positively with September temperature (r = 0.48, P < 0.05) and very strongly with precipitation totals from previous October until current February (r = 0.82, P < 0.001) showing that the water stored in the soil during the autumn and winter months prior to the growing season has a primordial effect on the growth of the given season. The effects of cork harvesting were analysed by comparing mean ring width, mean annual vessel area, vessel density (n°vessels/mm²), and vessel coverage (percentage of transverse surface occupied by vessels) between three mature cork oak trees and three young trees, for the period from 1987 to 1996, corresponding to the growth between two consecutive cork removals in the case of mature trees. In 1988, 1989 and 1996 (corresponding to the first and second years after cork removal, and 1996 to a year of cork removal), the ratios between ring widths of young versus mature trees was twice that for the rest of the period. However, an effect of cork removal indicated by eventual alterations in vessel size and distribution in the wood rings corresponding to the years 1988, 1989 and 1996 in the mature cork oaks was not observed.     

Site quality evaluation by classification tree: an application to cork quality in Sardinia

Cork harvesting and stopper production represent a major forest industry in Sardinia (Italy). The target of the present investigation was to evaluate the classification tree as a tool to discover possible relationships between microsite characteristics and cork quality. Seven main cork oak (Quercus suber) producing areas have been identified in Sardinia, for a total of more than 122,000 ha. Sixty-three sample trees, distributed among different geographical locations and microsite conditions, were selected. A soil profile near each sample tree was described, soil samples were collected and analysed. After debarking, cork quality of each sample tree was graded by an independent panel of experts. Microsites where trees had more than 50% of the extracted cork graded in the best quality class, according to the official quality standard in Italy, were labelled as prime microsites, the others as nonprime microsites. Relationships between a binary dummy variable (0 for nonprime microsites, 1 for prime microsites) and site factors were investigated using classification tree analysis to select the relevant variables and to define the classification scheme. Prime quality microsites for cork production proved to be characterised by elevation, soil phosphorus content and sandiness. Results have been compared with those of the more conventional parametric approach by logistic regression. The work demonstrates the advantages of the classification tree method. The model may be appropriate for classifications at landscape and stand mapping levels, where it is possible to sample a number of microsites and to evaluate distributional characteristics of model output, while its precision is only indicative when estimating the prime quality of single microsites.     

Population differences in juvenile survival under increasing drought are mediated by seed size in cork oak (Quercus suber L.)

Many studies have focused on the physiological and morphological responses of Quercus suber to high temperatures and summer drought. Nevertheless, our understanding of the potential of this species to cope with climate change is incomplete. An increase in severity, length and frequency of summer droughts is expected in the Iberian Peninsula over the 21st century. We investigated the potential of cork oak to adapt to climate change in a 4-year study comparing seedling survival from thirteen Spanish populations in a common garden. Acorn size was evaluated as a possible adaptive trait enhancing stress resistance during establishment. Populations originating from sites with the driest summers exhibited the highest survival rates under dry conditions. These populations were characterized by bigger acorns, suggesting selection for this trait across sites experiencing drought. Our findings reveal that northern populations are not well adapted to cope with increasing drought but suggest that they might cope well with the moderate drought increase expected for these zones. However, continental populations are intermediately adapted to dry conditions, while drastic and fast increases in summer droughts are expected to occur in these regions. This extensive and quick change will provide scant chances to adapt making the populations of cork oak in these regions particularly vulnerable to the future climatic conditions.     

Fire-driven oak regeneration in French Mediterranean ecosystems

Regeneration by seeds for cork oak (Quercus suber) and companion oaks (holm oak Quercus ilex and downy oak Quercus pubescens) is likely to be poor in the fire-prone Maures massif (southern France) but the causes are poorly known. Our objective was to assess the effective recruitment for these three oak species and their temporal pattern of recruitment, in order to determine the main limitation factors and the regeneration window of each species. We studied oak recruits (height <3 m) in naturally regenerated populations according to a gradient of fire recurrence and in five main vegetation types including shrublands and mixed mature woodlands. Fire recurrence was the main explanatory factor of oak recruitment, either directly or through vegetation type and microsite characteristics. The results indicate nil to low recruitment for holm oak and downy oak in shrublands, especially those recurrently burned and dominated by Cistus species. Cork oak recruited better than the other oaks in medium and high shrublands dominated by Erica arborea. In contrast, recruitment was high for holm and downy oak in mixed oak stands and mixed pine-oak stands that have not burned for decades. Microsite conditions such as coverage by litter and shrubs influenced oak recruitment, whereas landscape configuration and stand basal area had no influence. Our results suggest that strategic shrub-clearing, oak planting and protection of mixed oak woodlands as seed sources would help maintaining oak populations in the woodland–shrubland mosaic.     

Opportunity for hybridization between two oak species in mixed stands as monitored by the timing and intensity of pollen production

The opportunity of cross-pollination in mixed stands of two oak species (cork oak and holm oak) was studied by characterizing individual phenologies of flowering. In the spring of 1998 at one stand consisting of 64 marked trees, there was a period of 19 days when maximal pollen release in one species and stigma receptivity in the other occurred simultaneously, enabling interspecific gene flow in either direction. This happened in spite of an average time separation of 22 days between the two species, reflecting a considerable intraspecific variation in the timings of flowering. Flowering intensities (as estimated from male flower abundance) were high, but fruiting intensities were comparatively low. Shortly after pollination, considerable abortion of female flowers and early fruits was recorded. In 2000, the interspecific overlap of phenologies was drastically reduced due to a delay in cork oak flowering. On the other hand, the individual timings were repeatable for most trees, at least in holm oak. Two other mixed stands were subject of parallel studies, with similar results in all traits except for a less dramatic reduction in fruiting intensities. In spite of the high opportunity for cross-pollination in 1998, and given the lack of hybrids among the progenies from the subsequent fruiting season [Oliveira, P., Custódio, A.C., Branco, C., Reforço, I., Rodrigues, F., Varela, M.C., Meierrose, C., 2003. Hybrids between cork oak and holm oak: isoenzyme analysis. Forest Genet. 10, 283–298], it can be concluded that the prerequisite of cross-pollination is clearly insufficient for hybridization to succeed. Post-pollination processes must play an important role in the maintenance of reproductive isolation between the two species.     

Litter flammability in oak woodlands and shrublands of southeastern France

Characterizing the flammability of litter fuels is of major importance for assessing wildland fire ignition hazard. Here we compared the flammability of litter within a mosaic of Quercus suber (cork oak) woodlands and shrublands in a Mediterranean fire-prone area (Maures massif, southeastern France) to test whether the characteristics and the flammability of litter vary with the vegetation types. We tested experimentally the ignitability, the sustainability, the combustibility and the consumability of undisturbed (=non-reconstructed) litter samples with a point-source mode of ignition. Although the frequency of ignition was similar between all the vegetation types, we distinguished four groups having litter of specific composition and flammability: low and sparse shrublands dominated by Cistus species, medium shrublands with cork oak, high Erica shrublands with sparse cork oak woodlands, and mixed mature oak woodlands with Q. suber, Q. ilex and Q. pubescens. As these vegetation types corresponded to a specific range of past fire recurrence, we also tested the effect of the number of fires and the time since the last fire on litter flammability. Litters of plots recurrently burned had low ability to propagate flames and low flame sustainability. We discuss how the recent fire history can modify vegetation and litter flammability, and thus the fire ignition hazard.     

Fire-vegetation interplay in a mosaic structure of Quercus suber woodlands and Mediterranean maquis under recurrent fires

We examined the effects of fire recurrence on a mosaic structure of Quercus suber and Erica-Cistus shrubland communities of Southeastern France in order to improve the understanding of fire-vegetation interplay. Plant communities that were similar in 1959 (woodlands on shrublands called maquis on acidic soils) were compared along a gradient of fire recurrence, from 0 to 4 fires, with different time intervals between fires. The results showed that understory cover increased roughly with fire recurrence, whereas tree height, cover, density, stand basal area and litter depth decreased. Different fire recurrences along the past decades led to different vegetation types. High fire recurrence corresponded to maquis and sparse cork-oak woodlands while pure oak woodlands established in the absence of fire during the same period. In all, the diversity of tree diameter decreased with fire recurrence. High fire recurrence (3 or 4 fires in 50 years) led to a simple vertical structure of vegetation with a mono-layered shrub cover and few Quercus suber trees. In contrast, spatial connections between plants were maximal at longer time intervals, leading to a multi-layered vegetation. We finally discuss the potential implications of past fires on the behavior of future fires in the perspective of a sustainable management of these Mediterranean ecosystems.     

Composition and structure of natural mixed-oak stands in northern and central Portugal

Stand composition and structure of natural mixed-oak stands of common-oak (Quercus robur L.) and pyrenean-oak (Quercus pyrenaica Willd.) were studied. Diverse compositional and structural elements in early and late successional stand stages were analysed. The study was conducted in north and central Portugal where different natural mixed oak forests types are located. The following mixed-oak forest types involving common-oak and pyrenean-oak were studied: common-oak & other hardwoods; common-oak & cork-oak (Quercus suber L.); ash (Fraxinus angustifolia Vahl) & pyrenean-oak; and pyrenean-oak & madrone (Arbutus unedo L.). Measurements were made in early and late successional stand stages on the different mixed oak forest types. Different stand characteristics and indices were used to describe and compare stand structure and composition. The study showed changes in species diversity and stand structure. Most tree species in mature stands are present in early stages but with higher abundance. Shannon diversity index may change between 0.798 and 1.915. Significant differences on species diversity and abundance were found depending on the forest type and successional stage. Mature mixed-oak forests have high species diversity with an abundance of small to medium tree size species. Species distribution and diameter differentiation indices range from 0.30 to 0.70 and 0.52 to 0.82, respectively, revealing significant structural complexity. The average number of standing and downed dead trees was 265 and 83 trees ha⁻¹ for early and late stage, respectively, with 6.9 and 65.4 m³ ha⁻¹. Higher values of stand diversity index were 41 and 53 in more complex and developed forests. Later stand stages have complex structure, with a wider range of tree diameter distribution and higher degree of irregularity.     

Differences in ecophysiological responses to summer drought between seedlings of three deciduous oak species

Seasonal variations in water status and gas exchange were assessed in 2-year-old seedlings of Quercus pubescens Willd., Quercus ithaburensis Decaisne subsp. macrolepis (Kotsch) Hedge & Yalt, and Quercus frainetto Ten. during the 2004 growing season characterized by the typical summer drought conditions observed in Greece. The seedlings were grown under field conditions. After establishment half of the seedlings were well watered, while the rest received no irrigation other than natural rainfall. Measurements were conducted from April to October.     

Development of oak and beech seedlings planted under varying shelterwood densities and with different site preparation methods in southern Sweden

Effects of site preparation, shelterwood density and planting depth on the survival and growth of planted beech and oak seedlings were studied. Experiments were performed in one oak and one beech stand in southern Sweden. Two areas with different densities of shelterwood and one clearcut were established in each stand. Growth, damage and survival of the planted seedlings were observed for three years. Soil water potential was recorded weekly and radiation and soil temperatures were recorded continuously during the growing seasons.Neither site preparation methods nor planting depth affected oak seedling growth, when planting was carried out on fresh clearcuts or in shelterwoods, while growth of beech seedlings was positively affected by mounding. Growth of oak seedlings was inhibited by the shelterwood treatments. In beech seedlings, growth was lowest in dense shelterwood, while there was no difference in growth between seedlings on the clearcut area and in the shelterwood of low density. Differences in growth may be explained by differences in radiation and soil water potential.When planting was carried out on a one-year-old clearcut, site preparation improved the subsequent growth of oak and beech seedlings.     

Environmental and plant-based controls of water use in a Mediterranean oak stand

In this work, new information is reported on water relations of the Mediterranean oak species Quercus pyrenaica based on environmental and physiological measurements carried out during the growing seasons of 2006 and 2007. The interest in this species has increased due to its use in reforestation programs and its impacts on the water resources due to the extensive spontaneous afforestation after the abandonment of forest (firewood, charcoal, livestock, etc.) and agricultural activities, in general in some areas in the Mediterranean region and in particular, in the studied area (Sistema Central range, Spain). The objectives were to evaluate the long-term water use of this stand and its limit and to analyse the specific traits to cope with summer drought, especially the use of stem water storage and deep soil water reserves. Tree water stress associated with depletion of soil water reserves was not observed since the oak trees appeared to avoid a marked water stress using water reserves from deeper soil layers as summer drought progresses. The contribution of mean daily stem water storage to transpiration was low (4%), although it could be greater under dryer conditions. Only at the end of summer of 2006, the transpiration and canopy conductance were reduced due to soil drought. Despite the absence of marked water stress an upper limit was found in transpiration (slightly higher than 3 mm day⁻¹). The heavy use of soil water resources by this species (75% of available soil water in this study) should be considered when evaluating the impact of spontaneous afforestation and reforestation programs on water resources.     

Drought matters - Declining precipitation influences growth of Fagus sylvatica L. and Quercus robur L. in north-eastern Germany

For north-eastern Germany regional climate models project rising temperatures in combination with decreasing summer and increasing winter precipitation. The resulting overall drier conditions during the growing season will considerably impact forest growth there. We evaluate the consequences of increasing drought on the growth of the two locally most important broadleaf tree species common beech (Fagus sylvatica L.) and pedunculate oak (Quercus robur L.). Three mixed forests of beech and oak were sampled along a west-east gradient of declining precipitation. In total we used 257 ring-width samples from 133 trees to build six species and site specific chronologies. Additionally, we modelled the soil water budget for each site. We performed continuous and discontinuous (pointer year) analysis of climate-tree-growth relationships with particular emphasis on inter-annual-variations and their dependence upon climatic factors (temperature, precipitation, soil moisture) and on the stability of the obtained relationships. Results of climate-growth correlations together with pointer year analysis indicate a strong dependency of growth of both species from water availability, especially during early summer (June and July). General correlation pattern between growth and climate are similar for both species, but climate sensitivity of beech is generally higher. We identified drought as the main driver of negative growth depressions in both species. Increasing drought stress along the gradient is expressed in higher correlations to climatic variables, higher sensitivity (variance) of growth, and a higher number of negative pointer years for both species. For beech we also found a significant trend of decreasing average growth rates along the gradient. Growth superiority of beech compared to oak declines with decreasing precipitation. The relationships were generally stable throughout the 20th century. A rise of sensitivity together with a higher frequency of negative pointer years during the last decades suggests that increasing climatic variability together with rising temperatures might be influencing growth of Fagus at the more humid sites. If we substitute space by time it seems that already small changes in precipitation regime can have considerable impact, especially on the growth of beech. Other, more drought tolerant species like oak might gain competitive advantages under the projected climatic changes.     

Pisolithus tinctorius increases the size of nursery grown red oak seedlings

Inoculation of nursery beds with vegetative mycelia of P. tinctorius (Pt) significantly increased the number of salable 1–0 red oak seedlings produced in two of three years and the average stem height and diameter each year. The percentage of Pt mycorrhizal primary lateral roots on seedlings in the inoculated plots was significantly correlated with stem diameter and weight of lateral roots. The high rates of colonization of short roots by Pt in inoculated plots suggested the test fungus could successfully compete with indigenous fungi in forming mycorrhizae.     

Site quality of pedunculate oak (<Emphasis Type="Italic">Quercus robur</Emphasis> L.) stands in Galicia (northwest Spain)

Accurate estimates of forest productivity are required for sustainable forest management. Height of dominant and codominant trees at a reference age is often used as a measure of site productivity. Eight algebraic difference equations based on the models proposed by Sloboda, Bertalanffy-Richards, Korf, Hossfeld, and McDill-Amateis were tested on the transformed, longitudinal data structure that considers all possible growth intervals. Autocorrelation was modelled by expanding the error term as an autoregressive process according to this data structure. Generalized nonlinear least squares methods were used for model fitting. Several numerical and graphic analyses were used to compare the different candidate models. The relative error in dominant height prediction was used to select 60 years as the best reference age. An algebraic difference equation based on the model proposed by Korf provided the best compromise between biological and statistical aspects and produced the most adequate site index curves. This model is therefore recommended for height growth prediction and site classification of pedunculate oak stands in Galicia. The model is polymorphic and base-age invariant, having only one asymptote. Predictions of height for age intervals between t₁ and t₂ of more than 15 years should be considered with caution because of the associated critical error.     

Growth,nutrient, water relations,and gas exchange in a holm oak plantation in response to irrigation and fertilization

Eighty 6-years-old Quercus ilex L. subsp. ballota seedlings planted on a former agricultural land were subjected during two growing seasons to one of four treatments: fertilization and irrigation, irrigation, fertilization, and control. Seasonal and between-treatment variations on water relations, gas exchange parameters, growth and nutrient status were analyzed. Water potential () was related to climatic conditions. Thus, the frequent rain storms during the summer allow seedlings to maintain relatively high values, joined to moderate photosynthetic activity. Differences on , photosynthesis and stomatal conductance due to irrigation were shown at the onset and end of growth, related with lower water availability. Fertilization had a greater effect on growth than irrigation. Both fertilization and irrigation positively affected the relative increment on leaf nutrient concentration at the end of the second year. The results indicate that water availability was not a limited factor, thus irrigation is not justified; while fertilization improve growth.     

Effects of litter quality and parent material on organic matter characteristics and N-dynamics in Luxembourg beech and hornbeam forests

To test effects of litter quality and soil conditions on N-dynamics, we selected seven forests in Luxembourg dominated by beech (Fagus sylvatica, L.) and hornbeam (Carpinus betulus L.), and located on acid loam, decalcified marl or limestone, and measured organic matter characteristics, microbial C and N and net N-mineralization in a laboratory incubation experiment. Organic layer characteristics were significantly affected by species, with lower litter decay and higher accumulation under the less palatable beech, even on limestone. However, beech and hornbeam did not show any differences in N-cycling at all. Instead of species, N-cycling was affected by site conditions, albeit different than expected. Microbial N generally increased from acid loam to limestone, but acid loam showed higher net N-mineralization, especially in the organic layer. Also, acid loam showed high instead of low efficiency of N-mineralization per unit microbe, in both organic layer and mineral topsoil. In addition, acid loam showed net consumption of DOC instead of release in both soil layers, which suggests that not N, but C was a limiting factor to decomposition. In contrast, limestone showed low net N-mineralization in the organic layer, despite high mass and well-decomposed organic matter, and low efficiency of N-mineralization per unit microbe in both organic layer and mineral topsoil. DOC was net released instead of consumed, which supported that not C, but N was a limiting factor. The general lack of differences in net N-cycling between species, but relatively clear site effects, is discussed in relation to different microbial strategies. Acid soil may have high net N-release despite low biological activity, because N-requirements of fungi are also low, while in calcareous soil, high bacterial N-demand may counteract high gross N-release. Thus, species producing litter that decomposes rapidly may be planted to improve soil conditions and plant biodiversity, but litter quality effects on N-availability may be less important than soil conditions.     

Ungulate herbivory of regenerating conifers in relation to foliar nutrition and terpenoid production

Ungulate browsing greatly influences regeneration dynamics of some forest ecosystems, yet the relationship between browse susceptibility and foliar chemistry of forest tree seedlings is not well understood. We applied field fertilization (15N-9P-10K controlled-release fertilizer at 0, 20, 40, and 60 g per seedling) and investigated how subsequent changes in terpenoid production and foliar nutrition influence ungulate browse preference for Douglas-fir (Pseudotsugamenziesii (Mirb.) Franco), western hemlock (Tsuga heterophylla Raf. Sarg.), and western redcedar (Thujaplicata Donn ex D. Don) seedlings across four sites in northwestern Oregon, USA. Fertilization increased foliar N concentration of all three species, but above-ground growth of only Douglas-fir and western hemlock. Foliar monoterpene concentrations for western hemlock and western redcedar also increased at higher fertilization rates, while Douglas-fir monoterpene production was not affected by fertilization. Regardless of monoterpene levels, ungulate browse preference was greater for fertilized western hemlock seedlings. The opposite response, however, was observed for western redcedar at two of the four study sites where the likelihood of browse was greater for non-fertilized than fertilized seedlings. Differences in browse preference among species may depend on the type and amount of individual monoterpenes manufactured in response to fertilization. Western redcedar produce α- and β-thujone, oxygenated monoterpenes known to promote gastroenteritis and possibly inhibit microbial rumen activity, which were absent from Douglas-fir or western hemlock. Higher concentrations of α- and β-thujone associated with increasing fertilizer rate provide a plausible explanation as to why ungulates preferred non-fertilized western redcedar. Our results illustrate species-specific adaptation in browse avoidance and selective ungulate browsing behavior of individual trees as linked to foliar chemistry.