Re-assessment of the climatic controls on the carbon and water fluxes of a boreal aspen forest over 1996-2016: Changing sensitivity to long-term climatic conditions

Recent evidence suggests that the relationships between climate and boreal tree growth are generally non-stationary; however, it remains uncertain whether the relationships between climate and carbon (C) fluxes of boreal forests are stationary or have changed over recent decades. In this study, we used continuous eddy-covariance and microclimate data over 21 years (1996-2016) from a 100-year-old trembling aspen stand in central Saskatchewan, Canada to assess the relationships between climate and ecosystem C and water fluxes. Over the study period, the most striking climatic event was a severe, 3-year drought (2001-2003). Gross ecosystem production (GEP) showed larger interannual variability than ecosystem respiration (R_e ) over 1996-2016, but R_e was the dominant component contributing to the interannual variation in net ecosystem production (NEP) during post-drought years. The interannual variations in evapotranspiration (ET) and C fluxes were primarily driven by temperature and secondarily by water availability. Two-factor linear models combining precipitation and temperature performed well in explaining the interannual variation in C and water fluxes (R²  > .5). The temperature sensitivities of all three C fluxes (NEP, GEP and R_e ) declined over the study period (p < .05), and, as a result, the phenological controls on annual NEP weakened. The decreasing temperature sensitivity of the C fluxes may reflect changes in forest structure, related to the over-maturity of the aspen stand at 100 years of age, and exacerbated by high tree mortality following the severe 2001-2003 drought. These results may provide an early warning signal of driver shift or even an abrupt status shift of aspen forest dynamics. They may also imply a universal weakening in the relationship between temperature and GEP as forests become over-mature, associated with the structural and compositional changes that accompany forest ageing.

Spatial and temporal variations in global soil respiration and their relationships with climate and land cover

Soil respiration (R _s) represents the largest flux of CO₂ from terrestrial ecosystems to the atmosphere, but its spatial and temporal changes as well as the driving forces are not well understood. We derived a product of annual global R _s from 2000 to 2014 at 1 km by 1 km spatial resolution using remote sensing data and biome-specific statistical models. Different from the existing view that climate change dominated changes in R _s, we showed that land-cover change played a more important role in regulating R _s changes in temperate and boreal regions during 2000-2014. Significant changes in R _s occurred more frequently in areas with significant changes in short vegetation cover (i.e., all vegetation shorter than 5 m in height) than in areas with significant climate change. These results contribute to our understanding of global R _s patterns and highlight the importance of land-cover change in driving global and regional R _s changes.

Divergent long-term trends and interannual variation in ecosystem resource use efficiencies of a southern boreal old black spruce forest 1999-2017

Long-term trends in ecosystem resource use efficiencies (RUEs) and their controlling factors are key pieces of information for understanding how an ecosystem responds to climate change. We used continuous eddy covariance and microclimate data over the period 1999-2017 from a 120-year-old black spruce stand in central Saskatchewan, Canada, to assess interannual variability, long-term trends, and key controlling factors of gross ecosystem production (GEP) and the RUEs of carbon (CUE = net primary production [NPP]/GEP), light (LUE = GEP/absorbed photosynthetic radiation [APAR]), and water (WUE = GEP/evapotranspiration [E]). At this site, annual GEP has shown an increasing trend over the 19 years (p < 0.01), which may be attributed to rising atmospheric CO₂ concentration. Interannual variability in GEP, aside from its increasing trend, was most strongly related to spring temperatures. Associated with the significant increase in annual GEP were relatively small changes in NPP, APAR, and E, so that annual CUE showed a decreasing trend and annual LUE and WUE showed increasing trends over the 19 years. The long-term trends in the RUEs were related to the increasing CO₂ concentration. Further analysis of detrended RUEs showed that their interannual variation was impacted most strongly by air temperature. Two-factor linear models combining CO₂ concentration and air temperature performed well (R² ~0.60) in simulating annual RUEs. LUE and WUE were positively correlated both annually and seasonally, while LUE and CUE were mostly negatively correlated. Our results showed divergent long-term trends among CUE, LUE, and WUE and highlighted the need to account for the combined effects of climatic controls and the 'CO₂ fertilization effect' on long-term variations in RUEs. Since most RUE-based models rely primarily on one resource limitation, the observed patterns of relative change among the three RUEs may have important implications for RUE-based modeling of C fluxes.

Application of biochar and nitrogen influences fluxes of CO2, CH4 and N2O in a forest soil

Nitrogen (N) fertilization of forests for increasing carbon sequestration and wood volume is expected to influence soil greenhouse gas (GHG) emissions, especially to increase N₂O emissions. As biochar application is known to affect soil GHG emissions, we investigated the effect of biochar application, with and without N fertilization, to a forest soil on GHG emissions in a controlled laboratory study. We found that biochar application at high (10%) application rates increased CO₂ and N₂O emissions when applied without urea-N fertilizer. At both low (1%) and high biochar (10%) application rates CH₄ consumption was reduced when applied without urea-N fertilizer. Biochar application with urea-N fertilization did not increase CO₂ emissions compared to biochar amended soil without fertilizer. In terms of CO₂-eq, the net change in GHG emissions was mainly controlled by CO₂ emissions, regardless of treatment, with CH₄ and N₂O together accounting for less than 1.5% of the total emissions.

Nitrogen enrichment potential of biochar in relation to pyrolysis temperature and feedstock quality

Nitrogen (N) enrichment of biochar from both inorganic and organic waste N sources has the potential to add economic and environmental value through its use as a slow release N fertilizer. We investigated the sorption of N by, and its release from, biochar made at pyrolysis temperatures of 400, 500 and 600 °C from three feedstocks: poultry litter (PL with a carbon (C) to N ratio (C:N) of 14), softwood chips of spruce-pine-fir (SPF with a C:N of 470), and a 50:50 mixture of PL and SPF (PL/SPF). The prepared biochars were enriched with ammonium nitrate (AN) and urea ammonium nitrate (UAN). PL biochars had the lowest C content (50-56% C), but the highest pH (9.3-9.9), electrical conductivity (EC, 780-960 dS m(-1)), cation exchange capacity (CEC, 40-46 cmol kg(-1)), and N content (3.3-4.5%). While N content and hydrogen (H) to C atomic ratio (H:C) decreased with increasing pyrolysis temperature irrespective of the feedstock used, both pH and EC slightly increased with pyrolysis temperature for all feedstocks. The PL and SPF biochars showed similar H:C and also similar N sorption and N release at all pyrolysis temperatures. These biochars sorbed up to 5% N by mass, irrespective of the source of N. However, PL/SPF biochar performed poorly in sorbing N from either AN or UAN. Biochar H:C was found to be unrelated to N sorption rates, suggesting that physical adsorption on active surfaces was the main mechanism of N sorption in these biochars. There were minor differences between N sorbed from NO3-N and NH4-N among different biochars. Very small amounts of sorbed N (0.2-0.4 mg N g(-1) biochar) was released when extracted with 1 M KCl solution, indicating that the retained N was strongly held in complex bonds, more so for NH4-N because the release of NO3-N was 3-4 times greater than that of NH4-N. NH4-N sorption far exceeded the effective CEC of the biochars, thereby suggesting that most of the sorption may be due to physical entrapment of NH4(+) in biochar pores. The results of this study suggest that biochar can be used to remove excess N from poultry and dairy manure and be a good mitigation option for reducing N leaching and gaseous losses.

Looking deeper into the soil: biophysical controls and seasonal lags of soil CO2 production and efflux

We seek to understand how biophysical factors such as soil temperature (Ts), soil moisture (theta), and gross primary production (GPP) influence CO2 fluxes across terrestrial ecosystems. Recent advancements in automated measurements and remote-sensing approaches have provided time series in which lags and relationships among variables can be explored. The purpose of this study is to present new applications of continuous measurements of soil CO2 efflux (F0) and soil CO2 concentrations measurements. Here we explore how variation in Ts, theta, and GPP (derived from NASA's moderate-resolution imaging spectroradiometer [MODIS]) influence F0 and soil CO2 production (Ps). We focused on seasonal variation and used continuous measurements at a daily timescale across four vegetation types at 13 study sites to quantify: (1) differences in seasonal lags between soil CO2 fluxes and Ts, theta, and GPP and (2) interactions and relationships between CO2 fluxes with Ts, theta, and GPP. Mean annual Ts did not explain annual F0 and Ps among vegetation types, but GPP explained 73% and 30% of the variation, respectively. We found evidence that lags between soil CO2 fluxes and Ts or GPP provide insights into the role of plant phenology and information relevant about possible timing of controls of autotrophic and heterotrophic processes. The influences of biophysical factors that regulate daily F0 and Ps are different among vegetation types, but GPP is a dominant variable for explaining soil CO2 fluxes. The emergence of long-term automated soil CO2 flux measurement networks provides a unique opportunity for extended investigations into F0 and Ps processes in the near future.

Incorporating diffuse photosynthetically active radiation in a single-leaf model of canopy photosynthesis for a 56-year-old Douglas-fir forest

A simple top-down model of canopy photosynthesis (P) was developed and tested in this study. The model (referred to as the Q(e)-MM model) is P = alphaQ (e) P (max)/(alphaQ ( e ) + P (max)), alpha and P (max) are quantum-use efficiency and potential P, respectively. Q (e) is given by Q (d) (0) + kQ (b) (0), where Q (d) (0) and Q (b) (0) are the diffuse and direct photosynthetically active radiation (PAR) incident on the canopy, respectively. Q (e) can be considered to be the effective incident PAR contributing to P and k is a measure of the contribution of Q (b) (0) to Q (e). When k = 1, the Q(e)-MM model becomes the regular Michaelis-Menten type model of P (referred to as the MM model). A major objective of this study was to determine how well the Q(e)-MM model could estimate P of a 56-year-old coastal Douglas-fir stand. To this end, we parameterized the Q(e)-MM model using five and half years of eddy-covariance measurements of CO(2) flux above the Douglas-fir stand. The Q(e)-MM model, with the incorporation of a function of air temperature, accounted for 74% of the variance in over 34,000 half-hourly P measurements. P estimated using the Q(e)-MM model had no systematic errors with respect to Q (d) (0). Although the Q(e)-MM model has only one more parameter than the MM model, it accounted for 30% more variance in P than the latter when total incident PAR exceeded 900 micromol m(-2) s(-1). On average, k was found to be 0.22. We show that this small value of k reflects the significant effect of the scattering of the solar beam and the fraction of light-limited sunlit leaves. We also show that the success of the Q(e)-MM model was due to the fact that a large fraction of the sunlit leaves were light-limited as a result of their orientation to the solar beam.

Effects of mutual shading of tree crowns on prediction of photosynthetic light-use efficiency in a coastal Douglas-fir forest

Gross primary production (GPP) is often expressed as the product of absorbed photosynthetically active radiation and the efficiency (epsilon) with which a plant community uses absorbed radiation in biomass production. Light-use efficiency is affected by environmental stresses, and varies diurnally and seasonally. Uncertainty about epsilon can be a serious limitation when modeling GPP. An important determinant of epsilon is the amount and type of solar radiation incident on a canopy, because an abundance of light can trigger a photo-protective reaction, diminishing GPP. The radiation regime in a forest canopy is determined by the predominant sky conditions and by mutual shading of tree crowns. Shading effects, producing shifts in the amount of incident direct and diffuse solar radiation, have been largely ignored, however, because they depend on forest structure and are difficult to measure. We describe a new approach for estimating changes in mutual canopy shading throughout the day and year based on a canopy structure model derived from light detection and ranging (LiDAR). Proportions of canopy shading were then combined with eddy covariance data to assess the explanatory power for variance in epsilon by regression tree analysis over half-hourly, daily and weekly time scales. The approach explained between 75 and 97% of variance in epsilon, representing an increase of between 5 and 16% compared with models driven solely by meteorological variables.

Biophysical controls on rhizospheric and heterotrophic components of soil respiration in a boreal black spruce stand

We conducted a root-exclusion experiment in a 125-year-old boreal black spruce (Picea mariana (Mill.) BSP) stand in 2004 to quantify the physical and biological controls on temporal dynamics of the rhizospheric (R(r)) and heterotrophic (R(h)) components of soil respiration (R(s)). Annual R(r), R(h) and estimated moss respiration were 285, 269 and 57 g C m(-2) year(-1), respectively, which accounted for 47, 44 and 9% of R(s) (611 g C m(-2) year(-1)), respectively. A gradual transition from R(h)-dominated (winter, spring and fall) to R(r)-dominated (summer) respiration was observed during the year. Soil thawing in spring and the subsequent increase in soil water content (theta) induced a small and sustained increase in R(h) but had no effect on R(r). During the remainder of the growing season, no effect of theta was observed on either component of R(s). Both components increased exponentially with soil temperature (T(s)) during the growing season, but R(r) showed greater temperature sensitivity than R(h) (Q(10) of 4.0 and 3.0, respectively). Temperature-normalized variations in R(r) were highly correlated with eddy covariance estimates of gross ecosystem photosynthesis, and the correlation was greatest when R(r) was lagged by 24 days. Within diurnal cycles, variations in T(s) were highly coupled to variations in R(h) but were significantly decoupled from R(r). The patterns observed at both time scales strongly suggest that the flow of photosynthates to the rhizosphere is a key driver of belowground respiration processes but that photosynthate supply may control these processes in several ways.

Measurement of CO2 exchange between Boreal forest and the atmosphere

The Boreal forest is the world's second largest forested biome occupying the circumpolar region between 50 degrees N and 70 degrees N. This heterogeneous biome stores about 25% of all terrestrial carbon. We have reviewed EC measurements of CO2 exchange between the atmosphere and Boreal forests, and assessed progress in understanding the controlling processes. We have assessed net ecosystem productivity, the net balance between net primary productivity and heterotrophic respiration, measured using the EC method, for 38 Boreal forest sites. Gross ecosystem productivity has been estimated by adding day-time EC-measured CO2 fluxes to respiration estimated from night-time relationships between respiration and temperature. Maximum midday values of gross ecosystem productivity vary from 33 pmol m(-2) s(-1) for aspen to 6 micromol m(-2) s(-1) for larch stands. Long-term EC flux measurements, ongoing at nine Boreal sites, have shown the strong impact of spring weather and growing season water balance on annual net ecosystem productivity. Estimation of net biome production, incorporating the effects of disturbance resulting from forest fires and logging, has progressed significantly in recent years. After disturbance, summer measurements in Boreal chronosequences suggest that it takes about 10 years before growing season carbon uptake offsets the decomposition emissions. Small-scale exchange rate measurements using chambers and manipulative experiments such as stem girdling and soil heating help to understand the processes and mechanisms playing major roles in the carbon balance of terrestrial ecosystems. Aircraft EC flux measurements, convective boundary layer carbon budgets, and (13)C/12C changes in the atmosphere play an important role in validating estimates of regional carbon exchange based on scaled up EC measurements. Atmospheric inverse models are an important approach to studying regional and global carbon balance but need further improvement to yield reliable quantitative results.