Simulated projections of boreal forest peatland ecosystem productivity are sensitive to observed seasonality in leaf physiology

We quantified seasonal CO2 assimilation capacities for seven dominant vascular species in a wet boreal forest peatland then applied data to a land surface model parametrized to the site (ELM-SPRUCE) to test if seasonality in photosynthetic parameters results in differences in simulated plant responses to elevated CO2 and temperature. We collected seasonal leaf-level gas exchange, nutrient content and stand allometric data from the field-layer community (i.e., Maianthemum trifolium (L.) Sloboda), understory shrubs (Rhododendron groenlandicum (Oeder) Kron and Judd, Chamaedaphne calyculata (L.) Moench., Kalmia polifolia Wangenh. and Vaccinium angustifolium Alton.) and overstory trees (Picea mariana (Mill.) B.S.P. and Larix laricina (Du Roi) K. Koch). We found significant interspecific seasonal differences in specific leaf area, nitrogen content (by area; Na) and photosynthetic parameters (i.e., maximum rates of Rubisco carboxylation (Vcmax25°C), electron transport (Jmax25°C) and dark respiration (Rd25°C)), but minimal correlation between foliar Na and Vcmax25°C, Jmax25°C or Rd25°C, which illustrates that nitrogen alone is not a good correlate for physiological processes such as Rubisco activity that can change seasonally in this system. ELM-SPRUCE was sensitive to the introduction of observed interspecific seasonality in Vcmax25°C, Jmax25°C and Rd25°C, leading to simulated enhancement of net primary production (NPP) using seasonally dynamic parameters as compared with use of static parameters. This pattern was particularly pronounced under simulations with higher temperature and elevated CO2, suggesting a key hypothesis to address with future empirical or observational studies as climate changes. Inclusion of species-specific seasonal photosynthetic parameters should improve estimates of boreal ecosystem-level NPP, especially if impacts of seasonal physiological ontogeny can be separated from seasonal thermal acclimation.