Corrections

CELL BIOLOGY

CELL BIOLOGY

Correction for “Cell fate regulation by gelsolin in human gynecologic cancers,” by Mohammad R. Abedini, Pei-Wen Wang, Yu-Fang Huang, Mingju Cao, Cheng-Yang Chou, Dar-Bin Shieh, and Benjamin K. Tsang, which appeared in issue 40, October 7, 2014, of Proc Natl Acad Sci USA (111:14442–14447; first published September 22, 2014; 10.1073/pnas.1401166111). The authors note that the affiliation for Yu-Fang Huang and Cheng-Yang Chou should instead appear as Department of Obstetrics and Gynecology, National Cheng Kung University Hospital, College of Medicine, National Cheng Kung University, Tainan 704, Taiwan. The corrected author and affiliation lines appear below. The online version has been corrected.

Correction for “Evidence for multiple roles for grainyheadlike 2 in the establishment and maintenance of human mucociliary airway epithelium,” by Xia Gao, Christopher M. Vockley, Florencia Pauli, Kimberly M. Newberry, Yan Xue, Scott H. Randell, Timothy E. Reddy, and Brigid L. M. Hogan, which appeared in issue 23, June 4, 2013, of Proc Natl Acad Sci USA (110:9356– 9361; first published May 20, 2013; 10.1073/pnas.1307589110). The authors note that the title appeared incorrectly. The title should instead appear as “Evidence for multiple roles for grainyhead-like 2 in the establishment and maintenance of human mucociliary airway epithelium.” The online version has been corrected.

Mohammad R. Abedinia,b,c, Pei-Wen Wangd, Yu-Fang Huange, Mingju Caoa,b, Cheng-Yang Choue, Dar-Bin Shiehd,f, and Benjamin K. Tsanga,b

www.pnas.org/cgi/doi/10.1073/pnas.1421315111

a

Departments of Obstetrics and Gynaecology and Cellular and Molecular Medicine, Interdisciplinary School of Health Sciences, University of Ottawa, Ottawa, ON, Canada K1H 8L6; bChronic Disease Program, Ottawa Hospital Research Institute, Ottawa, ON, Canada K1H 8L6; cCellular and Molecular Medicine Research Center, Department of Pharmacology, Birjand University of Medical Sciences, Birjand 97178, Iran; dInstitute of Basic Medical Science, Institute of Oral Medicine and Department of Stomatology, National Cheng Kung University Hospital, College of Medicine, National Cheng Kung University, Tainan 704, Taiwan; eDepartment of Obstetrics and Gynecology, National Cheng Kung University Hospital, College of Medicine, National Cheng Kung University, Tainan 704, Taiwan; and fAdvanced Optoelectronic Technology Center and Center for Micro/Nano Science and Technology, National Cheng Kung University, Tainan 704, Taiwan www.pnas.org/cgi/doi/10.1073/pnas.1421407111

ECOLOGY

Correction for “Biogeographic variation in evergreen conifer needle longevity and impacts on boreal forest carbon cycle projections,” by Peter B. Reich, Roy L. Rich, Ying-Ping Wang, and Jacek Oleksyn, which appeared in issue 38, September 23, 2014, of Proc Natl Acad Sci USA (111:13703–13708; first published September 15, 2014; 10.1073/pnas.1216054110). The authors note that on page 13705, left column, second full paragraph, lines 11 and 12, “foliage and roots represent a greater and smaller fraction, respectively, of total biomass in forests in increasingly cold environments (26)” should instead appear as “foliage and roots represent a smaller and greater fraction, respectively, of total biomass in forests in increasingly cold environments (26).” www.pnas.org/cgi/doi/10.1073/pnas.1421130111

17684 | PNAS | December 9, 2014 | vol. 111 | no. 49

www.pnas.org

SPECIAL FEATURE

Biogeographic variation in evergreen conifer needle longevity and impacts on boreal forest carbon cycle projections Peter B. Reicha,b,1, Roy L. Richa, Xingjie Luc, Ying-Ping Wangd, and Jacek Oleksyna,e a Department of Forest Resources, University of Minnesota, St. Paul, MN 55108; bHawkesbury Institute for the Environment, University of Western Sydney, Penrith, NSW 2751, Australia; cCollege of Global Change and Earth System Science, Beijing Normal University, Beijing 100875, China; dCommonwealth Scientific and Industrial Research Organisation Marine and Atmospheric Research, Aspendale, VIC 3195, Australia; and eInstitute of Dendrology, Polish Academy of Sciences, PL-62-035, Kornik, Poland

Leaf life span is an important plant trait associated with interspecific variation in leaf, organismal, and ecosystem processes. We hypothesized that intraspecific variation in gymnosperm needle traits with latitude reflects both selection and acclimation for traits adaptive to the associated temperature and moisture gradient. This hypothesis was supported, because across 127 sites along a 2,160-km gradient in North America individuals of Picea glauca, Picea mariana, Pinus banksiana, and Abies balsamea had longer needle life span and lower tissue nitrogen concentration with decreasing mean annual temperature. Similar patterns were noted for Pinus sylvestris across a north–south gradient in Europe. These differences highlight needle longevity as an adaptive feature important to ecological success of boreal conifers across broad climatic ranges. Additionally, differences in leaf life span directly affect annual foliage turnover rate, which along with needle physiology partially regulates carbon cycling through effects on gross primary production and net canopy carbon export. However, most, if not all, global land surface models parameterize needle longevity of boreal evergreen forests as if it were a constant. We incorporated temperature-dependent needle longevity and %nitrogen, and biomass allocation, into a land surface model, Community Atmosphere Biosphere Land Exchange, to assess their impacts on carbon cycling processes. Incorporating realistic parameterization of these variables improved predictions of canopy leaf area index and gross primary production compared with observations from flux sites. Finally, increasingly low foliage turnover and biomass fraction toward the cold far north indicate that a surprisingly small fraction of new biomass is allocated to foliage under such conditions. biogeography

such conditions (14, 15). However, although evidence that ecotypic variation occurs is abundant (7, 8, 10), comprehensive characterization of ecologically important intraspecific variation across biogeographic gradients is lacking (16). Global land surface models that simulate carbon cycling of the world’s terrestrial biomes face distinct challenges for both species-poor biomes, such as the boreal forest, and species-rich biomes, such as tropical rain forests. For the latter, the challenge is to sufficiently well parameterize leaf and canopy properties given enormous species diversity within and among sites (e.g., >500 species for 125 sites along a 2,000-km gradient. For each species, individuals in colder, more northern environments had longer needle life span, highlighting its importance to evergreen ecological success. Incorporating biogeography of needle longevity into a global model improved predictions of forest productivity and carbon cycling and identified specific problems for models that ignore such variability.

T

he boreal forest is an enormous terrestrial biome, characterized by long, cold winters, low rates of productivity and nutrient cycling (1–4), and low richness of the dominant tree species. All of these tree species have large geographic ranges (5), and ecological theory suggests that populations should vary both genotypically and phenotypically in response to such variation. In fact, it has been long recognized that terrestrial plants in general, and boreal conifers specifically, exhibit ecotypic population differentiation that reflects influence of site origin on a variety of chemical, morphological, physiological, and allometric traits (6–12). It is also well known that biogeographic variation in the environment influences the phenotypic traits plastically realized by species across their ranges. Both phenotypic and genotypic processes therefore contribute to observed variation in traits across environmental gradients. Hence, across the boreal biome, conditions that lead to slow growth, such as short growing season, cool summer temperatures, and low soil nutrient supply (1, 2, 13), should be reflected by shifts toward more conservative traits (such as long leaf life span and low nutrient concentrations) considered to be advantageous under www.pnas.org/cgi/doi/10.1073/pnas.1216054110

Author contributions: P.B.R. and R.L.R. designed research; R.L.R. performed research; P.B.R., X.L., Y.-P.W., and J.O. analyzed data; and P.B.R. wrote the paper. The authors declare no conflict of interest. This article is a PNAS Direct Submission. C.V. is a guest editor invited by the Editorial Board. 1

To whom correspondence should be addressed. Email: [email protected].

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10. 1073/pnas.1216054110/-/DCSupplemental.

PNAS | September 23, 2014 | vol. 111 | no. 38 | 13703–13708

ECOLOGY

Edited by Cyrille Violle, Centre National de la Recherche Scientifique, Montpellier, France, and accepted by the Editorial Board October 22, 2013 (received for review January 30, 2013)

NL relation to geographic climate variation? (ii) Are proportional shifts in NL for a given climate gradient similar among species? (iii) Do other leaf traits often associated interspecifically with leaf life span (e.g., refs. 15 and 20) covary with NL intraspecifically in a similar fashion? (iv) What are the consequences of geographic patterns of needle traits to ecosystem and regional carbon cycling? To improve our understanding of these issues, we measured needle traits of the four dominant evergreen conifers from Minnesota to northern Canada, compiled similar data for the dominant evergreen species in Eurasia (Pinus sylvestris), characterized the relations of needle traits with climate variables for all five species, and incorporated such findings into a land surface model [Community Atmosphere Biosphere Land Exchange (CABLE)] that simulated biome-wide carbon and nitrogen cycling. We sampled foliage from naturally grown trees (range of height 2.5–5 m) at 127 sites ranging from Minnesota to Ontario, Manitoba, Saskatchewan, Alberta, and the Northwest Territories, in central Canada. We sampled sunlit upper canopy branches of trees ranging from 2.5 to 5 m in height, to standardize tree size and age across sites, as well as the light environment of the sampled branches. The latter was of paramount concern, because needle longevity in the crown of a conifer increases systematically with shading (by as much as 50–60% for spruce and fir). The four species were present (and sampled) at 83 (Pinus banksiana), 50 (Picea mariana), 45 (Picea glauca), and 21 (Abies balsamea) of the sites, respectively. We also use a literature compilation for Pinus sylvestris across northern Fenno-Scandinavia to extend our analysis beyond North America. Because several climatic factors covary with latitude and are hypothesized to be the mechanism behind the latitudinal pattern, we focus our analyses and interpretation on climate, not latitude. Additionally, it is difficult to ascertain which of the related climate metrics (including annual and seasonal measures of temperature and precipitation) are responsible for the observed latitudinal patterns (Methods and Supporting Information). Because mean annual temperature (MAT) was identified in the model selection process as the climate metric that explained the most variance in NL, we use it in the results presented herein. However, we discuss the ways in which multiple aspects of climate, and associated soil resource availability, can combine to influence these results, below and in Supporting Information.

The biogeographic range in NL is particularly pronounced in the Picea species—ranging from ∼6 y at their southern margins (∼46°N) to 15 y at the most northern locations sampled (∼62– 63°N), which are not far from their northern limit in that part of central Canada. Abies has almost as strong a shift in NL per unit shift in MAT as the Picea species but has a narrower range, and thus NL varies from 5 y at its southern limit (∼46°N) to 9 y at its northern (∼55°N) limit. P. sylvestris shows a similar proportional shift (more than doubling NL from 3.5 to 7.0 y from 58 to 69°N) as the North America Picea species. In contrast, P. banksiana, which has as large a geographic range as the Picea species, has a much weaker dependence of NL on MAT—it varies from 4.0 y in Minnesota to 5.5 in northern Canada. It is not clear why P. banksiana has such a noisy and modest relation of NL to MAT, compared with the other NA species or its European congener. It is possible that further to the south in its range it is relegated by competition to particularly poor sites that support slow growth and therefore longer NL than would otherwise occur, but data to test this hypothesis are not available. Regardless, there is a strong signal that NL decreases with MAT for these five species that are among the most dominant species in much of North American and Eurasian boreal forests. Needle %N increases significantly with MAT in P. banksiana, P. glauca, and P. sylvestris, with nonsignificant relations in the other two species. For data pooled across all five species, there are significant (P < 0.0001) relationships between %N and both NL (Fig. S2) and MAT. This indicates that boreal evergreen forests dominated by individuals with long NL (owing to prevalence of long NL species such as Picea, or to adaptive response to NL to cold temperatures, or both) also have low needle %N. In addition to illuminating intraspecific ecological strategies, the biogeographic patterns documented here have important implications for ecosystem-scale function. We used the CABLE land surface model (21, 22) to explore the impact of biogeographic variation in needle traits on evergreen boreal forest C cycling (Methods and Table 1). This exploration is important because most, perhaps all, widely recognized global land surface models consider all boreal forest as a single plant functional type (PFT) with a single set of trait values, and land surface models in general are inadequate with respect to biomass allocation (23, 24). To our knowledge, leading models such as Biome-BioGeochemical Cycles (Biome-BCG), Community Land Model version 4 (CLM-4), Lund–Potsdam–Jena Dynamic Global Model (LPJ) (and its offspring), and others do not include a temperate dependence of needle turnover rate (or biomass allocation) within the boreal evergreen PFT. Because CABLE does not recognize species

Results For all five species, NL increases with increasing latitude (Fig. S1) and decreases with MAT (Fig. 1). The relationship is strongest for both Picea species and P. sylvestris and weakest in P. banksiana.

10 black spruce 5

2.5

Abies balsamea -4

-2

Pinus banksiana

Picea glauca 0

2

4

Needle life span (years)

Needle life span (years)

15

6 10

5

2.5

Picea mariana -4

-2

Pinus sylvestris 0

2

4

Mean annual temperature (oC)

6

-4 -2 0 2 4 6 Mean annual temperature (oC)

Fig. 1. Mean NL (years) for four North American and one European evergreen conifer species in relation to MAT along a large latitudinal gradient: A. balsamea (logNL = 0.860–0.041*MAT; P = 0.0083, R2 = 0.34, n = 19), P. glauca (logNL = 0.911–0.044*MAT; P < 0.0001, R2 = 0.65, n = 46), P. banksiana (logNL = 0.650−0.014*MAT; P = 0.0052, R2 = 0.10, n = 80), P. mariana (logNL = 0.932–0.056*MAT; P < 0.0001, R2 = 0.66, n = 51), and P. sylvestris (logNL = 0.792–0.040*MAT; P < 0.0001, R2 = 0.65, n = 79).

13704 | www.pnas.org/cgi/doi/10.1073/pnas.1216054110

Reich et al.

Allocation

Needle life span

%N

Default cALLcNLc%N cALLfNLc%N fALLcNLc%N fALLfNLc%N cALLfNLf%N fALLfNLf%N

Constant Constant Constant ALL(MAT) ALL(MAT) Constant ALL(MAT)

Constant(3.76) Constant(6.02) NL(MAT) Constant(6.02) NL(MAT) NL(MAT) NL(MAT)

Constant(1.23) Constant(1.23) Constant(1.23) Constant(1.23) Constant(1.23) %N(MAT) %N(MAT)

The default version of CABLE represents the prior version (22, 33). Allocation (ALL) was either the default constant or a continuous function of mean annual temperature (MAT) (25). Needle life span (NL) was either constant at the prior default value (3.76), constant at the mean value (6.02) for the observed data, or modeled as a continuous function of MAT. Needle %N was either constant at the prior default value (1.23) or a continuous function of NL, as adjusted using a canopy depth extinction coefficient.

differences, we used a generalized algorithm to modify NL with MAT in a fashion that averages responses across the species available (Methods). Structure, parameterizations, and functional routines of CABLE are roughly similar to many other global land surface models, such as those mentioned above. Of relevance to this study, its parameterizations of PFTs and their key traits are similar to many other global land surface models. Thus, the examination herein of how changes in parameters and functions alter model performance should be generally reflective and not exclusive to CABLE. Temperature sensitivity of biomass allocation is likely an important but poorly understood regulator of C cycling in boreal forests (25), and models at present do not predict allocation well, even ignoring biogeographic variation (23, 24). Standing biomass distribution to leaves, stems, and roots is a joint function of turnover rate and initial allocation of new biomass. Given this, and that our goal includes assessing ecosystem impacts of changes in foliage turnover rate, we jointly examine impacts of temperature-sensitive biomass distribution (26) on CABLE predictions. Our modification of the biomass distribution in CABLE is based on analyses that show that foliage and roots represent a greater and smaller fraction, respectively, of total biomass in forests in increasingly cold environments (26), which is consistent with independent data showing proportionally greater carbon allocation below ground for Picea in colder, drier habitats (25). Incorporating temperature-dependent biomass partitioning yields the result that allocation of new biomass (NPP) to needles must be markedly low in the cold far north (Fig. 2), given that NL increases in colder, higher latitudes (Fig. 1), but the fraction of standing biomass located in needles is low at low MAT (26) (Fig. S3). The CABLE model as previously published (21) (called default hereafter) assumed NL of 3.76 y, root life span of 18 y (this includes and averages across fine roots, coarse roots, and woody roots), and wood life span of 70 y, needle %N of 1.23%, and a constant ratio of allocation of new biomass to needles, stems, and roots for evergreen needle leaf forest (21). To assess impacts of incorporating realistic biogeographic variation of biomass allocation (ALL), NL, and needle %N into the model, we compare these default results to permutations (Table 1) that incorporate these variables either as constants (cALLcNLc%N) or as functions of MAT (fALLfNLf%N), in various combinations. The constant values for ALL and %N are from the default CABLE, the constant value of NL is either the default value or the average across all sites for our field data, and the functions that vary with MAT are from our observations (Figs. 1 and 2 and Figs. S2 and S3), although for %N this is indirect, because we predict %N from temperature-dependent NL, to retain their covariance. We map these predictions spatially (Figs. 3 and 4) and also Reich et al.

1 Mg/ha

0.8

50

0.6

100

0.4

200

0.2

300

0 -10

0 MAT (0C)

10

Fig. 2. Representation of the fraction of new total biomass allocated to foliage in relation to MAT for evergreen conifer forests, for four aboveground biomass classes (50–300 Mg/ha). These curves result from the MATdependency of needle life span (Fig. 1) and fractional standing biomass distribution to foliage (Fig. S1). The vertical dashed lines represent the range of data for which they were developed. The allocation fraction is mathematically derived as a function of life spans of leaf, wood, and root; biomass fraction of leaf, wood, and root; and their relations with MAT and total biomass. The idea of calculating allocation fraction in this fashion is from a semianalytic solution to accelerate biogeochemical model spin-up (39).

PNAS | September 23, 2014 | vol. 111 | no. 38 | 13705

SPECIAL FEATURE ECOLOGY

Version

compare predictions for four north–south climate zones to values of canopy leaf area index (LAI) and gross primary production (GPP) from flux sites published by Luyssaert et al. (3) (Fig. 5). In the figures, moving from top to bottom (Figs. 3 and 4) or left to right (Fig. 5) generally moves from the default model with constant input parameter values, to new versions with one or two modified parameters, to the model where values of allocation, NL, and needle %N are all biogeographically varied. The default version of CABLE successfully predicts the north– south direction of LAI and GPP (Figs. 3–5) and slightly overpredicts GPP while underpredicting LAI, except in the coldest zone (Fig. 5). If a more realistic mean NL (6.02 y) as estimated from our data is used as a constant NL value everywhere (cALLcNLc%N), the model still gets the direction of LAI and GPP correct, but overpredicts GPP and LAI. This directly results from model logic—with a constant allocation routine and a constant NL, greater NL results in greater LAI. Incorporating our general temperature-dependent relationship for NL in the model instead of a constant value (cALLfNLc%N) also overpredicts GPP and LAI (Fig. 5) but additionally also fails to capture the directional signatures for GPP and LAI—with GPP stable and LAI decreasing with MAT. This is because with a constant allocation routine a realistically longer NL in the north unrealistically leads to higher LAI (and consequently higher GPP than would otherwise occur), despite a shorter growing season and colder temperatures. Incorporating MAT-dependent biomass allocation (and with NL constant, 6.02 y, fALLcNLc%N) also overpredicts GPP and LAI, but gets their directional variation correct (they are all higher toward the warmer south). The permutation (fALLfNLc %N) that includes temperature-dependent NL as well as temperature-dependent biomass allocation makes predictions that are similar. It is likely that the four new CABLE versions presented above all overpredict GPP in part because the constant % N (1.23%) used in the model is higher than the empirically observed average (∼0.9%), and LAI and %N jointly influence productivity in the CABLE model, just as observed empirically for both maximum instantaneous GPP and annual NPP (4). The default CABLE predictions may approximate GPP (Fig. 5) accurately because unrealistically high %N offsets unrealistically low LAI.

Allocaon fracon to foliage

Table 1. Versions used to assess CABLE model performance for boreal evergreen conifer forests

conservative foliage traits at the highest latitudes (25, 27). Thus, patterns of NL and %N are consistent with the notion that under cold, nutrient-poor conditions and a short growing season, which limit realized carbon gain, a successful ecological strategy involves trading off potential for high carbon gain by having leaves with low maximum photosynthetic capacity but associated low respiratory costs and having canopies that require only 6–12% replacement each year (14). Although we present results solely in relation to MAT, it is unclear whether it is growing-season temperatures (or precipitation), average annual temperatures (or precipitation), or relative moisture and nutrient supply that are mechanistically responsible for the low growth and associated long needle life span in high-latitude, cold locations (Supporting Information). Similarly, it is likely that low temperatures and dry conditions both contribute to high carbon allocation below ground in colder forests (25, 26). For the study sites on either continent, MAT was correlated strongly with mean annual precipitation (MAP), as well as with seasonal climate metrics and with the ratio between MAP and potential evapotranspiration. MAT explained as much or more of the variance in needle traits than any other climate metric and was selected statistically in model selection, hence its

Fig. 3. Patterns of simulated GPP across the boreal forest zone in vegetation under seven model versions, representing different model permutations with panels corresponding to modifications described in Table 1.

The model permutation (cALLfNLf%N) in which both temperature-dependent NL and NL-dependent %N are incorporated in the model (but with the default allocation scheme) does improve predictions on average compared with the model with temperature-dependent NL but with the default constant %N, but the directionality of LAI and GPP are not captured. The model (fALLfNLf%N) that incorporates both of the temperature-dependent functions (for NL and biomass partitioning) and the NL-dependent %N seems to do well in matching not only the direction of LAI and GPP but also the values observed empirically, although least well in the coldest zone (Fig. 5). Discussion Given the strong and consistent inverse relationship between leaf life span, nitrogen, and leaf carbon dioxide exchange rates documented interspecifically, including for conifers (14, 15, 20), it is highly likely that evergreen conifers in cold, northern zones of the boreal forest have low photosynthesis and respiration rates compared with either conspecifics or unrelated taxa in warmer, southern habitats. A combination of direct and indirect effects of low temperatures, including low nutrient availability and a short growing season, likely contribute to selection for increasingly 13706 | www.pnas.org/cgi/doi/10.1073/pnas.1216054110

Fig. 4. Patterns of simulated LAI across the boreal forest zone in vegetation under seven model versions, representing different model permutations with panels corresponding to modifications described in Table 1.

Reich et al.

Fig. 5. Comparison of GPP (A) and LAI (B) predictions from the CABLE land surface to observed FLUXNET data from Luyssaert et al. (3). The comparisons are shown for the predictions under seven versions, representing different model permutations (Table 1).

use herein. Nonetheless, it is likely that the documented patterns of needle trait variation with MAT are a response to seasonal as well as annual aspects of the thermal and moisture regimes along the latitudinal gradients. Boreal forests play key roles in global carbon cycling, so developing accurate ecosystem understanding is critical for improved carbon cycle modeling. The region we studied here represents 29% of global forests and was found to have the strongest positive carbon feedback to future warming using earth system models (28). Empirical observations have found that needle life span and needle %N (this paper) and the fraction of total stand biomass distributed to foliage (26) all vary with MAT within this biome, but such patterns have not been accounted for by any of the earth system models that were used for estimating the sensitivity of terrestrial carbon balance to warming. Using CABLE, we showed that implementing one or two of these three model changes (Table 1) as shown in Figs. 3–5 will either minimally alter or worsen the model predictions (compared with FLUXNET observations) compared with simulations without any of those changes as shown in the default output. Only with all three changes implemented in the model (fALLfNLf%N, Table 1) does the performance of CABLE become quite similar to the default CABLE for GPP and improved for LAI (Fig. 5). Therefore, implementing ecological realism into CABLE, and perhaps other global land surface models, may not lead to better agreements with observations or independent estimates of GPP. However, given that the default CABLE model uses unrealistically low and unrealistically geographically stable values for Reich et al.

Methods Needle traits for P. banksiana, P. mariana, P. glauca, and A. balsamea were sampled from individual 2.5- to 5-m-tall trees from 15- to 30-y-old forests, after full leaf expansion and before senescence in the summers of 2005 and 2006. Sampling of canopy leaf traits was made at 127 sites encompassing an area from 43.2° to 67.0°N, and 84.4° to 117.7°W. Dominant upland tree species (including five deciduous angiosperms in addition to the four conifers) were sampled in zones of ∼0.5 ° latitude (∼111 km) delineated before travel. The first suitable site observed with at least three species was sampled and supplemental sites were added as feasible until as many of the available target species as possible were sampled from each zone. All target species present were sampled at each site. At least two individuals were sampled from each site. Sampling occurred within forest stands and >25 m from disturbed areas such as roads. Only healthy/vigorous trees were sampled. Lateral branches with southern exposure from the lower portion of the top 25% of the crown were cut using telescoping clippers. From those branches, prior year (1-y-old) needles were sampled and needle life span was evaluated. Branches with male cones, cone scars (in P. bansksiana), disease, or insect damage were avoided. These species produce needles once per year in spring. Thus, needle life span was calculated by counting all living annual cohorts (from terminal scar to terminal scar) up to and including the youngest cohort with >50% but 50% of the original live foliage was estimated to the nearest 10th of a cohort (see refs. 18 and 33). Because cohorts younger than this one tended to retain nearly all of their foliage, and cohorts older than this tended to retain few, this technique provides a reasonable estimate of both mean cohort and mean needle longevity. Needles were scanned for leaf area and killed and dried in a 65 °C drying oven within 12 h of their collection. Leaves or needles from individual trees were ground into a composite sample and analyzed for analyzed for %N at the University of California, Davis Stable Isotope Facility (Department of Plant Sciences) using an Elementar Vario EL Cube or Micro Cube elemental analyzer (Elementar Analysensysteme GmbH) interfaced to a PDZ Europa 20-20 isotope ratio mass spectrometer (Sercon Ltd.). CABLE is the Australian community land surface model (22, 34) and has been implemented into the Australian community earth system model (ACCESS) (34). CABLE simulates the exchange of momentum, heat, water, and CO2 exchange between land surface and lower atmosphere at hourly time steps, plant growth, litter fall, and biogeochemical cycling of carbon, nitrogen, and phosphorus at daily time steps. As in most global land surface models, 13 plant functional types are used to approximate the variation of vegetated surface on Earth and each plant functional type is defined by the unique set of parameters values that are estimated from calibration or published literature (see ref. 21). CABLE performs well compared with other global land surface models at site level (see ref. 35) and can also reproduce major global biogeochemical pools and fluxes quite well (see ref. 21). The model has recently been used to study the effects of correlations among key leaf traits on the simulated gross primary production (36). For this study, our goal was to explore how different assumptions about needle life span, needle N concentration, and the fraction of NPP allocation to foliage affect the simulated GPP, NPP, and canopy LAI of evergreen needle leaf forests in North America, Europe, and Asia. For some of the CABLE scenarios we generated temperature-dependent NL (logNL = 0.8246–

0.0320*MAT°C) based on regression of NL in relation to MAT that included all available data for all five species (n = 845). In CABLE we also generated foliage N (mg/g) as a function of NL (logN = 0.1356–0.2252*logNL) based on regressing N vs. NL, including all available data for all five species (n = 391). To account for the decline in needle N with canopy depth within the canopy, as observed for evergreen conifer forests (e.g., refs. 37 and 38), we assumed that needle N decreases exponentially with the cumulative canopy LAI (see ref. 39). This adjustment has a similar impact on whole-canopy N as adjusting N by weighting the canopy by proportion in different needle age classes and adjusting for differences in %N with needle age. We drove CABLE by reusing the three-hourly 1° by 1° meteorological forcings for 1990 from the Global Soil Water Project II (40, 41) year after year until the simulated annual GPP and NPP reached their steady-state values. All of the simulated results as reported here are steady-state values. Additionally, to analyze the sensitivities of simulated annual GPP, NPP, and canopy LAI to a 3 °C warming, we recalculated the parameter values (fraction of NPP allocated to needle, needle life span, and %N) that vary with MAT using MAT+3 if applicable for every land point within the domain of this study for each of the seven scenarios, and ran CABLE by adding 3 °C to the surface air temperature across every time step. CABLE was run by reusing the GSWPII meteorological forcing for 1990 until steady-state values of GPP, NPP, and canopy LAI were obtained. The sensitivities of annual GPP, NPP, and canopy LAI to 3 °C warming were then calculated using the results from the simulations here and previously without 3 °C warming.

1. McGuire AD, et al. (2002) Environmental variation, vegetation distribution, carbon dynamics, and water/energy exchange in high latitudes. J Veg Sci 13:301–314. 2. Van Cleve K, Yarie J, Erickson R, Dyrness CT (1993) Nitrogen mineralization and nitrification in successional ecosystems on the Tanana River floodplain, interior Alaska. Can J Res 23:970–978. 3. Luyssaert S, et al. (2007) Photosynthesis drives anomalies in net carbon-exchange of pine forests at different latitudes. Glob Change Biol 13(10):2110–2127. 4. Reich PB (2012) Key canopy traits drive forest productivity. Proc Biol Sci 279(1736): 2128–2134. 5. Burns RM, Honkala BH (1990) Silvics of North America: 1. Conifers; 2. Hardwoods. Agriculture Handbook 654 (US Dept Agric, Forest Service, Washington, DC), Vol 2. 6. Turesson G (1922) The genotypical response of the plant species to the habitat. Hereditas 3:211–350. 7. Clausen J, Keck DD, Hiesey WM (1948) Experimental studies on the nature of species. III. Environmental responses of climatic races of Achillea. Publication 581 (Carnegie Inst Washington, Washington, DC). 8. Langlet O (1959) A cline or not a cline – a question of Scots pine. Silvae Genet 8:13–22. 9. Wright JW, Bull WI (1963) Geographic variation in Scotch pine. Silvae Genet 12:1–25. 10. Oleksyn J, Tjoelker MG, Reich PB (1992) Growth and biomass partitioning of populations of European Pinus sylvestris L. under simulated 50° and 60° N daylengths: evidence for photoperiodic ecotypes. New Phytol 120:561–574. 11. Zhang JW, Marshall JD (1995) Variation in carbon isotope discrimination and photosynthetic gas exchange among populations of Pseudotsuga menziesii and Pinus ponderosa in different environments. Funct Ecol 9:402–412. 12. Jonas CS, Geber MA (1999) Variation among populations of Clarkia unguiculata (Onagraceae) along altitudinal and latitudinal gradients. Am J Bot 86(3):333–343. 13. Van Cleve K, Oliver L, Schlentner R, Viereck LA, Dyrness CT (1983) Productivity and nutrient cycling in taiga forest ecosystems. Can J Res 13:747–766. 14. Reich PB, Walters MB, Ellsworth DS (1992) Leaf lifespan in relation to leaf, plant and stand characteristics among diverse ecosystems. Ecol Monogr 62:365–392. 15. Wright IJ, et al. (2004) The worldwide leaf economics spectrum. Nature 428(6985):821–827. 16. Albert CH, et al. (2010) Intraspecific functional variability: Extent, structure and sources of variation. J Ecol 98:604–613. 17. Pravdin LF (1969) Scots Pine Variation, Interspecific Taxonomy and Selection (Izd. Nauka, Moscow, 1964) trans. Israel Program for Scientific Translations, Jerusalem. 18. Reich PB, Oleksyn J, Modrzynski J, Tjoelker MG (1996) Evidence that longer needle retention of spruce and pine populations at high elevations and high latitudes is largely a phenotypic response. Tree Physiol 16(7):643–647. 19. Schoettle AW (1990) The interaction between leaf longevity and shoot growth and foliar biomass per shoot in Pinus contorta at two elevations. Tree Physiol 7(1–4):209–214. 20. Reich PB, Walters MB, Ellsworth DS (1997) From tropics to tundra: Global convergence in plant functioning. Proc Natl Acad Sci USA 94(25):13730–13734. 21. Wang YP, Law RM, Pak B (2010) A global model of carbon, nitrogen and phosphorus cycles for the terrestrial biosphere. Biogeosciences 7:2261–2282.

22. Wang YP, et al. (2011) Diagnosing errors in a land surface model (CABLE) in the time and frequency domains. J Geophys Res 116:G01034. 23. Wolf A, Field CB, Berry JA (2011) Allometric growth and allocation in forests: A perspective from FLUXNET. Ecol Appl 21(5):1546–1556. 24. Wolf A, et al. (2011) Forest biomass allometry in global land surface models. Global Biogeochem Cycles 25:GB3015. 25. Vogel JG, et al. (2008) Carbon allocation in boreal black spruce forests across regions varying in soil temperature and precipitation. Glob Change Biol 14:1503–1516. 26. Reich PB, et al. (2014) Temperature drives global patterns in forest biomass distribution in leaves, stems, and roots. Proc Natl Acad Sci USA 111:13721–13726. 27. Kikuzawa K, Onoda Y, Wright IJ, Reich PB (2013) Mechanisms underlying global temperature-related patterns in leaf longevity. Glob Ecol Biogeogr 22:982–993. 28. Yoshikawa C, Kawamiya M, Kato T, Yamanaka Y, Matsuno T (2008) Geographical distribution of the feedback between future climate change and the carbon cycle. J Geophys Res 113:G03002. 29. Friedlingstein P, et al. (2006) Climate-carbon cycle feedback analysis, results from the C4MIP model intercomparison. J Clim 19:3337–3353. 30. Nemani RR, et al. (2003) Climate-driven increases in global terrestrial net primary production from 1982 to 1999. Science 300(5625):1560–1563. 31. Kimball JS, et al. (2007) Recent climate-driven increases in vegetation productivity for the western arctic: Evidence of an acceleration of the northern terrestrial carbon cycle. Earth Interact 11:1–30. 32. Reich PB, Oleksyn J (2008) Climate warming will reduce growth and survival of Scots pine except in the far north. Ecol Lett 11(6):588–597. 33. Perez Hardguindeguy N, et al. (2013) New handbook for standardised measurement of plant functional traits worldwide. Aust J Bot 61(3):167–234. 34. Kowalczyk EA, et al. (2013) The land surface model component of ACCESS: description and impact on the simulated surface climatology. Aust Meteorol Oceanogr J 63:65–82. 35. Abramowitz G, Pitman AJ, Gupta H, Kowalczyk EA, Wang YP (2007) Systematic bias in land surface models. J Hydrometeorol 8:989–1001. 36. Wang YP, et al. (2012) Correlations among leaf traits provide a significant constraint on the estimate of global gross primary production. Geophys Res Lett 39:L19405. 37. Dang Q-L, et al. (1997) Profiles of PAR, nitrogen and photosynthetic capacity in the boreal forest: Implications for scaling from leaf to canopy. J Geophys Res BOREAS Special Issue 102(D24):28845–28860. 38. Nippert JB, Marshall JD (2003) Sources of variation in ecophysiological parameters in Douglas-fir and grand fir canopies. Tree Physiol 23(9):591–601. 39. Lloyd J, et al. (2010) Optimisation of photosynthetic carbon gain and within-canopy gradients of associated foliar traits for Amazon forest trees. Biogeosciences 7:1833–1859. 40. Dirmeyer PA, et al. (2006) GSWP-2- multimodel analysis and implications for our perception of the land surface. Bull Am Meteorol Soc 87:1381–1397. 41. Xia JY, Luo YQ, Wang YP, Weng ES, Hararuk O (2012) A semi-analytic solution to accelerate spin-up of a coupled carbon and nitrogen land model to steady state. Geosci Model Dev 5:1259–1271.

13708 | www.pnas.org/cgi/doi/10.1073/pnas.1216054110

ACKNOWLEDGMENTS. P.B.R. thanks the Wilderness Research Foundation; the Institute on the Environment (University of Minnesota), and the U.S. National Science Foundation (DEB-1234162); Y.-P.W. thanks the Australian Department of Climate Change and Energy Efficiency for financial support. X.L. thanks the China Scholarship Council for financial support for his study in Australia. J.O.’s work was supported by National Science Center (Poland) Grant 2011/02/A/NZ9/00108 and Ministry of Science and Higher Education (Poland) Grant N N304 375738.

Reich et al.