Accuracy, realism and general applicability of European forest models.

Forest models are instrumental for understanding and projecting the impact of climate change on forests. A considerable number of forest models have been developed in the last decades. However, few systematic and comprehensive model comparisons have been performed in Europe that combine an evaluation of modelled carbon and water fluxes and forest structure. We evaluate 13 widely used, state-of-the-art, stand-scale forest models against field measurements of forest structure and eddy-covariance data of carbon and water fluxes over multiple decades across an environmental gradient at nine typical European forest stands. We test the models' performance in three dimensions: accuracy of local predictions (agreement of modelled and observed annual data), realism of environmental responses (agreement of modelled and observed responses of daily gross primary productivity to temperature, radiation and vapour pressure deficit) and general applicability (proportion of European tree species covered). We find that multiple models are available that excel according to our three dimensions of model performance. For the accuracy of local predictions, variables related to forest structure have lower random and systematic errors than annual carbon and water flux variables. Moreover, the multi-model ensemble mean provided overall more realistic daily productivity responses to environmental drivers across all sites than any single individual model. The general applicability of the models is high, as almost all models are currently able to cover Europe's common tree species. We show that forest models complement each other in their response to environmental drivers and that there are several cases in which individual models outperform the model ensemble. Our framework provides a first step to capturing essential differences between forest models that go beyond the most commonly used accuracy of predictions. Overall, this study provides a point of reference for future model work aimed at predicting climate impacts and supporting climate mitigation and adaptation measures in forests.

Linking canopy-scale mesophyll conductance and phloem sugar δ13 C using empirical and modelling approaches.

Interpreting phloem carbohydrate or xylem tissue carbon isotopic composition as measures of water-use efficiency or past tree productivity requires in-depth knowledge of the factors altering the isotopic composition within the pathway from ambient air to phloem contents and tree ring. One of least understood of these factors is mesophyll conductance (gm ). We formulated a dynamic model describing the leaf photosynthetic pathway including seven alternative gm descriptions and a simple transport of sugars from foliage down the trunk. We parameterised the model for a boreal Scots pine stand and compared simulated gm responses with weather variations. We further compared the simulated δ13 C of new photosynthates among the different gm descriptions and against measured phloem sugar δ13 C. Simulated gm estimates of the seven descriptions varied according to weather conditions, resulting in varying estimates of phloem δ13 C during cold/moist and warm/dry periods. The model succeeded in predicting a drought response and a postdrought release in phloem sugar δ13 C indicating suitability of the model for inverse prediction of leaf processes from phloem isotopic composition. We suggest short-interval phloem sampling during and after extreme weather conditions to distinguish between mesophyll conductance drivers for future model development.

A new method to estimate branch biomass from terrestrial laser scanning data by bridging tree structure models.

BACKGROUND AND AIMS: Branch biomass and other attributes are important for estimating the carbon budget of forest stands and characterizing crown structure. As destructive measuring is time-consuming and labour-intensive, terrestrial laser scanning (TLS) as a solution has been used to estimate branch biomass quickly and non-destructively. However, branch information extraction from TLS data alone is challenging due to occlusion and other defects, especially for estimating individual branch attributes in coniferous trees. METHODS: This study presents a method, entitled TSMtls, to estimate individual branch biomass non-destructively and accurately by combining tree structure models and TLS data. The TSMtls method constructs the stem-taper curve from TLS data, then uses tree structure models to determine the number, basal area and biomass of individual branches at whorl level. We estimated the tree structural model parameters from 122 destructively measured Scots pine (Pinus sylvestris) trees and tested the method on six Scots pine trees that were first TLS-scanned and later destructively measured. Additionally, we estimated the branch biomass using other TLS-based approaches for comparison. KEY RESULTS: Tree-level branch biomass estimates derived from TSMtls showed the best agreement with the destructive measurements [coefficient of variation of root mean square error (CV-RMSE) = 9.66 % and concordance correlation coefficient (CCC) = 0.99], outperforming the other TLS-based approaches (CV-RMSE 12.97-57.45 % and CCC 0.43-0.98 ). Whorl-level individual branch attributes estimates produced from TSMtls showed more accurate results than those produced from TLS data directly. CONCLUSIONS: The results showed that the TSMtls method proposed in this study holds promise for extension to more species and larger areas.

Estimating canopy gross primary production by combining phloem stable isotopes with canopy and mesophyll conductances.

Gross primary production (GPP) is a key component of the forest carbon cycle. However, our knowledge of GPP at the stand scale remains uncertain, because estimates derived from eddy covariance (EC) rely on semi-empirical modelling and the assumptions of the EC technique are sometimes not fully met. We propose using the sap flux/isotope method as an alternative way to estimate canopy GPP, termed GPPiso/SF , at the stand scale and at daily resolution. It is based on canopy conductance inferred from sap flux and intrinsic water-use efficiency estimated from the stable carbon isotope composition of phloem contents. The GPPiso/SF estimate was further corrected for seasonal variations in photosynthetic capacity and mesophyll conductance. We compared our estimate of GPPiso/SF to the GPP derived from PRELES, a model parameterized with EC data. The comparisons were performed in a highly instrumented, boreal Scots pine forest in northern Sweden, including a nitrogen fertilized and a reference plot. The resulting annual and daily GPPiso/SF estimates agreed well with PRELES, in the fertilized plot and the reference plot. We discuss the GPPiso/SF method as an alternative which can be widely applied without terrain restrictions, where the assumptions of EC are not met.

Extending the range of applicability of the semi-empirical ecosystem flux model PRELES for varying forest types and climate.

Applications of ecosystem flux models on large geographical scales are often limited by model complexity and data availability. Here we calibrated and evaluated a semi-empirical ecosystem flux model, PREdict Light-use efficiency, Evapotranspiration and Soil water (PRELES), for various forest types and climate conditions, based on eddy covariance data from 55 sites. A Bayesian approach was adopted for model calibration and uncertainty quantification. We applied the site-specific calibrations and multisite calibrations to nine plant functional types (PFTs) to obtain the site-specific and PFT-specific parameter vectors for PRELES. A systematically designed cross-validation was implemented to evaluate calibration strategies and the risks in extrapolation. The combination of plant physiological traits and climate patterns generated significant variation in vegetation responses and model parameters across but not within PFTs, implying that applying the model without PFT-specific parameters is risky. But within PFT, the multisite calibrations performed as accurately as the site-specific calibrations in predicting gross primary production (GPP) and evapotranspiration (ET). Moreover, the variations among sites within one PFT could be effectively simulated by simply adjusting the parameter of potential light-use efficiency (LUE), implying significant convergence of simulated vegetation processes within PFT. The hierarchical modelling of PRELES provides a compromise between satellite-driven LUE and physiologically oriented approaches for extrapolating the geographical variation of ecosystem productivity. Although measurement errors of eddy covariance and remotely sensed data propagated a substantial proportion of uncertainty or potential biases, the results illustrated that PRELES could reliably capture daily variations of GPP and ET for contrasting forest types on large geographical scales if PFT-specific parameterizations were applied.

Plant respiration: Controlled by photosynthesis or biomass?

Two simplifying hypotheses have been proposed for whole-plant respiration. One links respiration to photosynthesis; the other to biomass. Using a first-principles carbon balance model with a prescribed live woody biomass turnover, applied at a forest research site where multidecadal measurements are available for comparison, we show that if turnover is fast the accumulation of respiring biomass is low and respiration depends primarily on photosynthesis; while if turnover is slow the accumulation of respiring biomass is high and respiration depends primarily on biomass. But the first scenario is inconsistent with evidence for substantial carry-over of fixed carbon between years, while the second implies far too great an increase in respiration during stand development-leading to depleted carbohydrate reserves and an unrealistically high mortality risk. These two mutually incompatible hypotheses are thus both incorrect. Respiration is not linearly related either to photosynthesis or to biomass, but it is more strongly controlled by recent photosynthates (and reserve availability) than by total biomass.

Early snowmelt significantly enhances boreal springtime carbon uptake.

We determine the annual timing of spring recovery from space-borne microwave radiometer observations across northern hemisphere boreal evergreen forests for 1979-2014. We find a trend of advanced spring recovery of carbon uptake for this period, with a total average shift of 8.1 d (2.3 d/decade). We use this trend to estimate the corresponding changes in gross primary production (GPP) by applying in situ carbon flux observations. Micrometeorological CO2 measurements at four sites in northern Europe and North America indicate that such an advance in spring recovery would have increased the January-June GPP sum by 29 g⋅C⋅m-2 [8.4 g⋅C⋅m-2 (3.7%)/decade]. We find this sensitivity of the measured springtime GPP to the spring recovery to be in accordance with the corresponding sensitivity derived from simulations with a land ecosystem model coupled to a global circulation model. The model-predicted increase in springtime cumulative GPP was 0.035 Pg/decade [15.5 g⋅C⋅m-2 (6.8%)/decade] for Eurasian forests and 0.017 Pg/decade for forests in North America [9.8 g⋅C⋅m-2 (4.4%)/decade]. This change in the springtime sum of GPP related to the timing of spring snowmelt is quantified here for boreal evergreen forests.

New insights into the covariation of stomatal, mesophyll and hydraulic conductances from optimization models incorporating nonstomatal limitations to photosynthesis.

Optimization models of stomatal conductance (gs ) attempt to explain observed stomatal behaviour in terms of cost--benefit tradeoffs. While the benefit of stomatal opening through increased CO2 uptake is clear, currently the nature of the associated cost(s) remains unclear. We explored the hypothesis that gs maximizes leaf photosynthesis, where the cost of stomatal opening arises from nonstomatal reductions in photosynthesis induced by leaf water stress. We analytically solved two cases, CAP and MES, in which reduced leaf water potential leads to reductions in carboxylation capacity (CAP) and mesophyll conductance (gm ) (MES). Both CAP and MES predict the same one-parameter relationship between the intercellular : atmospheric CO2 concentration ratio (ci /ca ) and vapour pressure deficit (VPD, D), viz. ci /ca  ≈ ξ/(ξ + âˆšD), as that obtained from previous optimization models, with the novel feature that the parameter ξ is determined unambiguously as a function of a small number of photosynthetic and hydraulic variables. These include soil-to-leaf hydraulic conductance, implying a stomatal closure response to drought. MES also predicts that gs /gm is closely related to ci /ca and is similarly conservative. These results are consistent with observations, give rise to new testable predictions, and offer new insights into the covariation of stomatal, mesophyll and hydraulic conductances.

CASSIA--a dynamic model for predicting intra-annual sink demand and interannual growth variation in Scots pine.

The control of tree growth vs environment by carbon sources or sinks remains unresolved although it is widely studied. This study investigates growth of tree components and carbon sink-source dynamics at different temporal scales. We constructed a dynamic growth model 'carbon allocation sink source interaction' (CASSIA) that calculates tree-level carbon balance from photosynthesis, respiration, phenology and temperature-driven potential structural growth of tree organs and dynamics of stored nonstructural carbon (NSC) and their modifying influence on growth. With the model, we tested hypotheses that sink demand explains the intra-annual growth dynamics of the meristems, and that the source supply is further needed to explain year-to-year growth variation. The predicted intra-annual dimensional growth of shoots and needles and the number of cells in xylogenesis phases corresponded with measurements, whereas NSC hardly limited the growth, supporting the first hypothesis. Delayed GPP influence on potential growth was necessary for simulating the yearly growth variation, indicating also at least an indirect source limitation. CASSIA combines seasonal growth and carbon balance dynamics with long-term source dynamics affecting growth and thus provides a first step to understanding the complex processes regulating intra- and interannual growth and sink-source dynamics.

Assimilate transport in phloem sets conditions for leaf gas exchange.

Carbon uptake and transpiration in plant leaves occurs through stomata that open and close. Stomatal action is usually considered a response to environmental driving factors. Here we show that leaf gas exchange is more strongly related to whole tree level transport of assimilates than previously thought, and that transport of assimilates is a restriction of stomatal opening comparable with hydraulic limitation. Assimilate transport in the phloem requires that osmotic pressure at phloem loading sites in leaves exceeds the drop in hydrostatic pressure that is due to transpiration. Assimilate transport thus competes with transpiration for water. Excess sugar loading, however, may block the assimilate transport because of viscosity build-up in phloem sap. Therefore, for given conditions, there is a stomatal opening that maximizes phloem transport if we assume that sugar loading is proportional to photosynthetic rate. Here we show that such opening produces the observed behaviour of leaf gas exchange. Our approach connects stomatal regulation directly with sink activity, plant structure and soil water availability as they all influence assimilate transport. It produces similar behaviour as the optimal stomatal control approach, but does not require determination of marginal cost of water parameter.

Duration of shoot elongation in Scots pine varies within the crown and between years.

BACKGROUND AND AIMS: Shoot elongation in boreal and temperate trees typically follows a sigmoid pattern where the onset and cessation of growth are related to accumulated effective temperature (thermal time). Previous studies on leader shoots suggest that while the maximum daily growth rate depends on the availability of resources to the shoot, the duration of the growth period may be an adaptation to long-term temperature conditions. However, other results indicate that the growth period may be longer in faster growing lateral shoots with higher availability of resources. This study investigates the interactions between the rate of elongation and the duration of the growth period in units of thermal time in lateral shoots of Scots pine (Pinus sylvestris). METHODS: Length development of 202 lateral shoots were measured approximately three times per week during seven growing seasons in 2-5 trees per year in a mature stand and in three trees during one growing season in a sapling stand. A dynamic shoot growth model was adapted for the analysis to determine (1) the maximum growth rate and (2) the thermal time reached at growth completion. The relationship between those two parameters and its variation between trees and years was analysed using linear mixed models. KEY RESULTS: The shoots with higher maximum growth rate within a crown continued to grow for a longer period in any one year. Higher July-August temperature of the previous summer implied a higher requirement of thermal time for growth completion. CONCLUSIONS: The results provide evidence that the requirement of thermal time for completion of lateral shoot extension in Scots pine may interact with resource availability to the shoot both from year to year and among shoots in a crown each year. If growing season temperatures rise in the future, this will affect not only the rate of shoot growth but its duration also.

Quantification of forest carbon flux and stock uncertainties under climate change and their use in regionally explicit decision making: Case study in Finland.

Uncertainties are essential, yet often neglected, information for evaluating the reliability in forest carbon balance projections used in national and regional policy planning. We analysed uncertainties in the forest net biome exchange (NBE) and carbon stocks under multiple management and climate scenarios with a process-based ecosystem model. Sampled forest initial state values, model parameters, harvest levels and global climate models (GCMs) served as inputs in Monte Carlo simulations, which covered forests of the 18 regions of mainland Finland over the period 2015-2050. Under individual scenarios, the results revealed time- and region-dependent variability in the magnitude of uncertainty and mean values of the NBE projections. The main sources of uncertainty varied with time, by region and by the amount of harvested wood. Combinations of uncertainties in the representative concentration pathways scenarios, GCMs, forest initial values and model parameters were the main sources of uncertainty at the beginning, while the harvest scenarios dominated by the end of the simulation period, combined with GCMs and climate scenarios especially in the north. Our regionally explicit uncertainty analysis was found a useful approach to reveal the variability in the regional potentials to reach a policy related, future target level of NBE, which is important information when planning realistic and regionally fair national policy actions.

Effect of forest management choices on carbon sequestration and biodiversity at national scale.

Forest management methods and harvest intensities influence wood production, carbon sequestration and biodiversity. We devised different management scenarios by means of stakeholder analysis and incorporated them in the forest growth simulator PREBAS. To analyse impacts of harvest intensity, we used constraints on total harvest: business as usual, low harvest, intensive harvest and no harvest. We carried out simulations on a wall-to-wall grid in Finland until 2050. Our objectives were to (1) test how the management scenarios differed in their projections, (2) analyse the potential wood production, carbon sequestration and biodiversity under the different harvest levels, and (3) compare different options of allocating the scenarios and protected areas. Harvest level was key to carbon stocks and fluxes regardless of management actions and moderate changes in proportion of strictly protected forest. In contrast, biodiversity was more dependent on other management variables than harvesting levels, and relatively independent of carbon stocks and fluxes.

Role of data uncertainty when identifying important areas for biodiversity and carbon in boreal forests.

Forest conservation plays a central role in meeting national and international biodiversity and climate targets. Biodiversity and carbon values within forests are often estimated with models, introducing uncertainty to decision making on which forest stands to protect. Here, we explore how uncertainties in forest variable estimates affect modelled biodiversity and carbon patterns, and how this in turn introduces variability in the selection of new protected areas. We find that both biodiversity and carbon patterns were sensitive to alterations in forest attributes. Uncertainty in features that were rare and/or had dissimilar distributions with other features introduced most variation to conservation plans. The most critical data uncertainty also depended on what fraction of the landscape was being protected. Forests of highest conservation value were more robust to data uncertainties than forests of lesser conservation value. Identifying critical sources of model uncertainty helps to effectively reduce errors in conservation decisions.

Role of land cover in Finland's greenhouse gas emissions.

We present regionally aggregated emissions of greenhouse gases (GHG) from five land cover categories in Finland: artificial surfaces, arable land, forest, waterbodies, and wetlands. Carbon (C) sequestration to managed forests and unmanaged wetlands was also assessed. Models FRES and ALas were applied for emissions (CO2, CH4, N2O) from artificial surfaces and agriculture, and PREBAS for forest growth and C balance. Empirical emission coefficients were used to estimate emissions from drained forested peatland (CH4, N2O), cropland (CO2), waterbodies (CH4, CO2), peat production sites and undrained mires (CH4, CO2, N2O). We calculated gross emissions of 147.2 ± 6.8 TgCO2eq yr-1 for 18 administrative units covering mainland Finland, using data representative of the period 2017-2025. Emissions from energy production, industrial processes, road traffic and other sources in artificial surfaces amounted to 45.7 ± 2.0 TgCO2eq yr-1. The loss of C in forest harvesting was the largest emission source in the LULUCF sector, in total 59.8 ± 3.3 TgCO2eq yr-1. Emissions from domestic livestock production, field cultivation and organic soils added up to 12.2 ± 3.5 TgCO2eq yr-1 from arable land. Rivers and lakes (13.4 ± 1.9 TgCO2eq yr-1) as well as undrained mires and peat production sites (14.7 ± 1.8 TgCO2eq yr-1) increased the total GHG fluxes. The C sequestration from the atmosphere was 93.2 ± 13.7 TgCO2eq yr-1. with the main sink in forest on mineral soil (79.9 ± 12.2 TgCO2eq yr-1). All sinks compensated 63% of total emissions and thus the net emissions were 53.9 ± 15.3 TgCO2eq yr-1, or a net GHG flux per capita of 9.8 MgCO2eq yr-1.

Modelling the regional potential for reaching carbon neutrality in Finland: Sustainable forestry, energy use and biodiversity protection.

The EU aims at reaching carbon neutrality by 2050 and Finland by 2035. We integrated results of three spatially distributed model systems (FRES, PREBAS, Zonation) to evaluate the potential to reach this goal at both national and regional scale in Finland, by simultaneously considering protection targets of the EU biodiversity (BD) strategy. Modelling of both anthropogenic emissions and forestry measures were carried out, and forested areas important for BD protection were identified based on spatial prioritization. We used scenarios until 2050 based on mitigation measures of the national climate and energy strategy, forestry policies and predicted climate change, and evaluated how implementation of these scenarios would affect greenhouse gas fluxes, carbon storages, and the possibility to reach the carbon neutrality target. Potential new forested areas for BD protection according to the EU 10% protection target provided a significant carbon storage (426-452 TgC) and sequestration potential (- 12 to - 17.5 TgCO2eq a-1) by 2050, indicating complementarity of emission mitigation and conservation measures. The results of the study can be utilized for integrating climate and BD policies, accounting of ecosystem services for climate regulation, and delimitation of areas for conservation.

Modeling forest stand dynamics from optimal balances of carbon and nitrogen.

We formulate a dynamic evolutionary optimization problem to predict the optimal pattern by which carbon (C) and nitrogen (N) are co-allocated to fine-root, leaf, and wood production, with the objective of maximizing height growth rate, year by year, in an even-aged stand. Height growth is maximized with respect to two adaptive traits, leaf N concentration and the ratio of fine-root mass to sapwood cross-sectional area. Constraints on the optimization include pipe-model structure, the C cost of N acquisition, and agreement between the C and N balances. The latter is determined by two models of height growth rate, one derived from the C balance and the other from the N balance; agreement is defined by identical growth rates. Predicted time-courses of maximized height growth rate accord with general observations. Across an N gradient, higher N availability leads to greater N utilization and net primary productivity, larger trees, and greater stocks of leaf and live wood biomass, with declining gains as a result of saturation effects at high N availability. Fine-root biomass is greatest at intermediate N availability. Predicted leaf and fine-root stocks agree with data from coniferous stands across Finland. Optimal C-allocation patterns agree with published observations and model analyses.

A model bridging waterlogging, stomatal behavior and water use in trees in drained peatland.

Waterlogging causes hypoxic or anoxic conditions in soils, which lead to decreases in root and stomatal hydraulic conductance. Although these effects have been observed in a variety of plant species, they have not been quantified continuously over a range of water table depths (WTD) or soil water contents (SWC). To provide a quantitative theoretical framework for tackling this issue, we hypothesized similar mathematical descriptions of waterlogging and drought effects on whole-tree hydraulics and constructed a hierarchical model by connecting optimal stomata and soil-to-leaf hydraulic conductance models. In the model, the soil-to-root conductance is non-monotonic with WTD to reflect both the limitations by water under low SWC and by hypoxic effects associated with inhibited oxygen diffusion under high SWC. The model was parameterized using priors from literature and data collected over four growing seasons from Scots pine (Pinus sylvestris L.) trees grown in a drained peatland in Finland. Two reference models (RMs) were compared with the new model, RM1 with no belowground hydraulics and RM2 with no waterlogging effects. The new model was more accurate than the RMs in predicting transpiration rate (fitted slope of measured against modeled transpiration rate = 0.991 vs 0.979 (RM1) and 0.984 (RM2), R2 = 0.801 vs 0.665 (RM1) and 0.776 (RM2)). Particularly, RM2's overestimation of transpiration rate under shallow water table conditions (fitted slope = 0.908, R2 = 0.697) was considerably reduced by the new model (fitted slope = 0.956, R2 = 0.711). The limits and potential improvements of the model are discussed.

Sources and sinks of greenhouse gases in the landscape: Approach for spatially explicit estimates.

Climate change mitigation is a global response that requires actions at the local level. Quantifying local sources and sinks of greenhouse gases (GHG) facilitate evaluating mitigation options. We present an approach to collate spatially explicit estimated fluxes of GHGs (carbon dioxide, methane and nitrous oxide) for main land use sectors in the landscape, to aggregate, and to calculate the net emissions of an entire region. Our procedure was developed and tested in a large river basin in Finland, providing information from intensively studied eLTER research sites. To evaluate the full GHG balance, fluxes from natural ecosystems (lakes, rivers, and undrained mires) were included together with fluxes from anthropogenic activities, agriculture and forestry. We quantified the fluxes based on calculations with an anthropogenic emissions model (FRES) and a forest growth and carbon balance model (PREBAS), as well as on emission coefficients from the literature regarding emissions from lakes, rivers, undrained mires, peat extraction sites and cropland. Spatial data sources included CORINE land use data, soil map, lake and river shorelines, national forest inventory data, and statistical data on anthropogenic activities. Emission uncertainties were evaluated with Monte Carlo simulations. Artificial surfaces were the most emission intensive land-cover class. Lakes and rivers were about as emission intensive as arable land. Forests were the dominant land cover in the region (66%), and the C sink of the forests decreased the total emissions of the region by 72%. The region's net emissions amounted to 4.37 ± 1.43 Tg CO2-eq yr-1, corresponding to a net emission intensity 0.16 Gg CO2-eq km-2 yr-1, and estimated per capita net emissions of 5.6 Mg CO2-eq yr-1. Our landscape approach opens opportunities to examine the sensitivities of important GHG fluxes to changes in land use and climate, management actions, and mitigation of anthropogenic emissions.

Modeling acclimation of photosynthesis to temperature in evergreen conifer forests.

• In this study, we used a canopy photosynthesis model which describes changes in photosynthetic capacity with slow temperature-dependent acclimations. • A flux-partitioning algorithm was applied to fit the photosynthesis model to net ecosystem exchange data for 12 evergreen coniferous forests from northern temperate and boreal regions. • The model accounted for much of the variation in photosynthetic production, with modeling efficiencies (mean > 67%) similar to those of more complex models. The parameter describing the rate of acclimation was larger at the northern sites, leading to a slower acclimation of photosynthesis to temperature. The response of the rates of photosynthesis to air temperature in spring was delayed up to several days at the coldest sites. Overall photosynthesis acclimation processes were slower at colder, northern locations than at warmer, more southern, and more maritime sites. • Consequently, slow changes in photosynthetic capacity were essential to explaining variations of photosynthesis for colder boreal forests (i.e. where acclimation of photosynthesis to temperature was slower), whereas the importance of these processes was minor in warmer conifer evergreen forests.