Forest biomass and volume estimation using airborne LiDAR in a cool-temperate forest of northern Hokkaido,Japan

Trees are recognized as a carbon reservoir, and precise and convenient methods for forest biomass estimation are required for adequate carbon management. Airborne light detection and ranging (LiDAR) is considered to be one of the solutions for large-scale forest biomass evaluation. To clarify the relationship between mean canopy height determined by airborne LiDAR and forest timber volume and biomass of cool-temperate forests in northern Hokkaido, Japan, we conducted LiDAR observations covering the total area of the Teshio Experimental Forest (225 km²) of Hokkaido University and compared the results with ground surveys and previous studies. Timber volume and aboveground tree carbon content of the studied forest stands ranged from 101.43 to 480.40 m³ ha^–1 and from 30.78 to 180.54 MgC ha^–1, respectively. The LiDAR mean canopy height explained the variation among stands well (volume: r² = 0.80, RMSE = 55.04 m³ ha^–1; aboveground tree carbon content: r² = 0.78, RMSE = 19.10 MgC ha^–1) when one simple linear regression equation was used for all types (hardwood, coniferous, and mixed) of forest stands. The determination of a regression equation for each forest type did not improve the prediction power for hardwood (volume: r² = 0.84, RMSE = 62.66 m³ ha^–1; aboveground tree carbon content: r² = 0.76, RMSE = 27.05 MgC ha^–1) or coniferous forests (volume: r² = 0.75, RMSE = 51.07 m³ ha^–1; aboveground tree carbon content: r² = 0.58, RMSE = 19.00 MgC ha^–1). Thus, the combined regression equation that includes three forest types appears to be adequate for practical application to large-scale forest biomass estimation.     

Biomass storage and stand structure in a conservation unit in the Atlantic Rainforest—The role of big trees

This study represents a small-scale approach to forest structure and biomass in the Atlantic Rainforest in Brazil and provides information on an ecosystem in which there still is a lack of data in this regard.The project was carried out in the National Park “Serra dos Orgãos” in the state of Rio de Janeiro, which is one of the largest remnants of continuous forest in this area. This forest is marked by a mosaic of forest types differing in tree composition and structure. Within this heterogeneous habitat the stand structure in three investigation plots was assessed to estimate the above-ground dry biomass (AGB) for all trees with a dbh ≥ 5 cm.This study indicates the structural diversity of the Atlantic Rainforest. Trees with a dbh > 30 cm were represented by 6% of all sampled individuals (18 out of 318 trees), but contributed 72% of total estimated AGB. The results suggest that big trees in the Atlantic Rainforest may contribute more into total AGB as reported for other tropical rainforests. Small-scale structural approaches like this study are able to form an initiating framework of more detailed results and help to improve estimates on biomass amounts and therefore on carbon storage capacity.     

Toward consistent measurements of carbon accumulation: A multi-site assessment of biomass and basal area increment across Europe

The use of tree-ring data in carbon cycle research has so far been limited because traditional study designs are not geared toward quantifying forest carbon accumulation. Existing studies that assessed biomass increment from tree rings were often confined to individual sites and used inconsistent sampling schemes. We applied a consistent biomass-oriented sampling design at five managed forest sites located in different climate zones to assess the annual carbon accumulation in above-ground woody tissues (i.e. stems and branches) and its climate response. Radial growth and biometric measurements were combined to reconstruct the annual biomass increment in individual trees and upscaled to the site level. In addition to this, we estimated that 32–60 trees are required at these five sites to robustly quantify carbon accumulation rates. Tree dimensions and growth rates varied considerably among sites as a function of differing stand density, climatic limitations, and management interventions. Accordingly, mean site-level carbon accumulation rates between 65 g C m⁻² y⁻¹ and 225 g C m⁻² y⁻¹ were reconstructed for the 1970–2009 period. A comparison of biomass increment with the widely used basal area increment (BAI) revealed very similar growth trends but emphasized the merits of biomass assessments due to species-specific BAI/biomass relationship. Our study illustrates the benefits and challenges of combining tree-ring data with biometric measurements and promotes the consistent application of a standardized sampling protocol across large spatial scales. It is thus viewed as a conceptual basis for future use of tree-ring data to approach research questions related to forest productivity and the terrestrial carbon balance.     

Estimating biomass and carbon mitigation of temperate coniferous forests using spectral modeling and field inventory data

Realizing the importance of forest carbon monitoring and reporting in climate change, the present study was conducted to derive spectrally modeled aboveground biomass and mitigation using Landsat data in combination with sampled field inventory data in the coniferous forests of Western Himalaya. After conducting preliminary survey in 2009, 90 quadrats (45 each for calibration and validation) of 0.1 ha were laid in six forest types for recording field inventory data viz. diameter at breast height, height, slope and aspect. Biomass carbon (Mg ha^− 1) was worked out for different forest types and crown density classes (open with 10–40% crown density and closed with > 40% crown density) using recommended volume equations, ratios and factors. Biomass carbon map (aboveground + belowground) was generated for the entire region using geospatial techniques. Normalized difference vegetation index (NDVI) was generated and spectral values were extracted to establish relation (R² = 0.72, p < 0.01) with the field inventory data. The model developed was validated (R² = 0.73, p < 0.01) with 45 sample observations not used earlier for predicting and generating biomass carbon map (2009) for the entire region. The data from field based inventory indicates highest total biomass carbon (171.40, σ ± 23.19) Mg ha^− 1 for Fir–Spruce (closed) which has relatively more mature girth classes and low tree density. This value was found to be significantly higher than other forest types. Lowest biomass carbon was observed for Blue Pine (open) (37.15, σ ± 11.82) Mg ha^− 1. The NDVI values for the entire region ranged from 0 to 0.62 and consequently the spectrally derived aboveground biomass carbon varied from 0 to 600 Mg ha^− 1. The study demonstrates the application of mapping, spectral responses and sampled field inventory for type wise assessment of carbon mitigation in temperate coniferous forests of Himalayas.     

Loss of aboveground forest biomass and landscape biomass variability in Missouri,US

Disturbance regimes and forests have changed over time in the eastern United States. We examined effects of historical disturbance (circa 1813 to 1850) compared to current disturbance (circa 2004 to 2008) on aboveground, live tree biomass (for trees with diameters ≥13 cm) and landscape variation of biomass in forests of the Ozarks and Plains landscapes in Missouri, USA. We simulated 10,000 one-hectare plots using random diameters generated from parameters of diameter distributions limited to diameters ≥13 cm and random densities generated from density estimates. Area-weighted mean biomass density (Mg/ha) for historical forests averaged 116 Mg/ha, ranging from 54 Mg/ha to 357 Mg/ha by small scale ecological subsections within Missouri landscapes. Area-weighted mean biomass density for current forests averaged 82 Mg/ha, ranging from 66 Mg/ha to 144 Mg/ha by ecological subsection for currently forested land. Biomass density of current forest was greater than historical biomass density for only 2 of 23 ecological subsections. Current carbon sequestration of 292 TgC on 7 million ha of forested land is less than half of the estimated historical total carbon sequestration of 693 TgC on 12 million ha. Cumulative tree cutting disturbances over time have produced forests that have less aboveground tree biomass and are uniform in biomass compared to estimates of historical biomass, which varied across Missouri landscapes. With continued relatively low rates of forest disturbance, current biomass per ha will likely increase to historical levels as the most competitive trees become larger in size and mean number of trees per ha decreases due to competition and self-thinning. Restoration of large diameter structure and forested extent of upland woodlands and floodplain forests could fulfill multiple conservation objectives, including carbon sequestration.     

Carbon profiles of typical forest types across Europe assessed with CO2FIX

This paper presents for 16 typical forest types across Europe a standard carbon sequestration profile. The study was carried out with the model CO2FIX which was parameterised with local yield table data and additional required parameters. CO2FIX quantifies the carbon of the forest ecosystem–soil–wood products chain at the stand level. To avoid misleading results annual net sequestration rates are not presented here, because these strongly fluctuate in time. Therefore, only its advancing mean is presented as a more reliable indicator. This avoids a great deal of uncertainty for policy makers. The variation between forest types is large, but mean sequestration rates mostly peak after some 38 years (with a net source lasting up to 15 years after afforestation) at an average value of 2.98 Mg C ha⁻¹ per year (ranging between forest types from 4.1 to 1.15). After 200 years, the net sequestration rate saturates to a value of 0.8 Mg C ha⁻¹ per year (ranging from 1.4 to 0.13). The long-term mean carbon stock in tree biomass and products amounts on average to 114 Mg C ha⁻¹ (ranging from 52 to 196).     

Extremely low fine root biomass in Larix sibirica forests at the southern drought limit of the boreal forest

Mongolia's Larix sibirica forests at the southern fringe of the Eurosiberian boreal forest belt are exposed not only to very low winter temperatures, but also to frequent summer droughts. It is not completely known how Siberian larch adapts to these stressors. We examined whether (i) these forests differ in their fine root bio- and necromass from more humid boreal forests further in the North and (ii) inter-annual fluctuations in fine root biomass are related to tree vitality. In two exceptionally dry summers, we found only 4–5 g DM m^?2 of fine root biomass (in 0–20 cm depth), which is far less than typical conifer fine root biomass figures from boreal forests (c. 200–400 g m^?2) and the lowest forest fine root biomass reported worldwide; in a moist summer, fine root biomass was 20 fold higher. In contrast to fine root biomass, both necromass and non-tree root mass were high in all three years. From the large increase of fine root biomass in the moist summer and the generally high root necromass, we conclude that drought-induced fine root dieback was the likely cause of the very small amount of live root mass in the dry summers. Larch fine roots seem to be more drought-sensitive than shoots, since marked needle loss did not occur under the extreme conditions.     

A power-driven increment borer for sampling high-density tropical wood

High-density hardwood trees with large diameters have been found to damage manually operated increment borers, thus limiting their use in the tropics. Therefore, we herein report a new, low-cost gasoline-powered sampling system for high-density tropical hardwood trees with large diameters. This system provides increment cores 15 mm in diameter and up to 1.35 m in length, allowing minimally invasive sampling of tropical hardwood tree species, which, up to the present, could not be collected by conventional 5 or 10 mm increment borers. This system provides a single core sample with ample amount of wood for multidisciplinary analyses, including ring width, stable isotope and wood anatomical measurements. The borer never gets stuck inside stems, even in hollowed trees, cores will never twist during coring, and the gasoline drill gives ample flexibility in the field. It is anticipated that the dendrochronological community will find our technique very useful in the pursuit of tropical tree ring research.     

Dendrochronological research of Scots pine (Pinus sylvestris L.) radial growth in vicinity of industrial pollution

The aim of this research is to investigate changes in the annual radial increment of Scots pine (Pinus sylvestris L.) in the vicinity of intensive (3–10 km) and moderate (11–20 km) industrial pollution during different growth periods (growth promotion, inhibition, and recovery). Low level of emission was beneficial for tree growth during the growth promotion period, and the annual radial increment in the zones of intensive and moderate pollution increased by approximately 15–25% and 10%, respectively. Severe loss was reported to forests during the growth inhibition period when nitrogen and sulphur dioxide emissions were 37–40 thousand (thou.) tons per year. About 40–45% tree radial increment loss was observed in the stands closest to the pollution source, and 15–20% loss was observed for the most distant stands. The stabilization of radial growth decrease and the beginning of recovery of damaged stands began in 1988–1992, when the annual amount of industrial emissions and environmental pollution were considerably reduced. The stabilization of radial growth and the initiation of recovery after pollution reduction were high for the most damaged stands. Their radial increment was stable and close to that of the control stands in 2000–2011. Stands with less damage growing further from the pollution source were recovered earlier, and their radial increment stabilized near the control increment in 1995–1999. The results of linear regression analysis demonstrated that the impact of pollution is different for stands growing at different distances from the plant, and the impact decreases with distance (R² = 0.78 and R² = 0.75, respectively; p < 0.05).     

Ecological Footprint: Refining the carbon Footprint calculation

Within the Ecological Footprint methodology, the carbon Footprint component is defined as the regenerative forest capacity required to sequester the anthropogenic carbon dioxide emissions that is not absorbed by oceans. A key parameter of the carbon Footprint is the Average Forest Carbon Sequestration (AFCS), which is calculated from the net carbon sequestration capacity of forests ecosystems.The aim of this paper is to increase the clarity and transparency of the Ecological Footprint by reviewing the rationale and methodology behind the carbon Footprint component, and updating a key factor in its calculation, the AFCS. Multiple calculation options have been set to capture different rates of carbon sequestration depending on the degree of human management of three types of forest considered (primary forests, other naturally regenerated forests and planted forests). Carbon emissions related to forest wildfires and soil as well as harvested wood product have been included for the first time in this update of the AFCS calculation. Overall, a AFCS value range of 0.73 ± 0.37 t C ha⁻¹ yr⁻¹ has been identified. The resulting carbon Footprint and Ecological Footprint values have then been evaluated based on this value range. Results confirm that human demand for ecosystem services is beyond the biosphere's natural capacity to provide them.     

Tree species diversity and community composition in a human-dominated tropical forest of Western Ghats biodiversity hotspot,India

The structure, function, and ecosystem services of tropical forest depend on its species richness, diversity, dominance, and the patterns of changes in the assemblages of tree populations over time. Long-term data from permanent vegetation plots have yielded a wealth of data on the species diversity and dynamics of tree populations, but such studies have only rarely been undertaken in tropical forest landscapes that support large human populations. Thus, anthropogenic drivers and their impacts on species diversity and community structure of tropical forests are not well understood. Here we present data on species diversity, community composition, and regeneration status of tropical forests in a human-dominated landscape in the Western Ghats of southern India. Enumeration of 40 plots (50 m × 20 m) results a total of 106 species of trees, 76 species of saplings and 79 species of seedlings. Detrended Correspondence Analysis ordination of the tree populations yielded five dominant groups, along disturbance and altitudinal gradients on the first and second axes respectively. Abundant species of the area such as Albizia amara, Nothopegia racemosa and Pleiospermum alatum had relatively few individuals in recruiting size classes. Our data indicate probable replacement of rare, localized, and old-growth ‘specialists’ by disturbance-adapted generalists, if the degradation is continuing at the present scale.     

High plant species richness indicates management-related disturbances rather than the conservation status of forests

There is a wealth of smaller-scale studies on the effects of forest management on plant diversity. However, studies comparing plant species diversity in forests with different management types and intensity, extending over different regions and forest stages, and including detailed information on site conditions are missing. We studied vascular plants on 1500 20 m × 20 m forest plots in three regions of Germany (Schwäbische Alb, Hainich-Dün, Schorfheide-Chorin). In all regions, our study plots comprised different management types (unmanaged, selection cutting, deciduous and coniferous age-class forests, which resulted from clear cutting or shelterwood logging), various stand ages, site conditions, and levels of management-related disturbances. We analyzed how overall richness and richness of different plant functional groups (trees, shrubs, herbs, herbaceous species typically growing in forests and herbaceous light-demanding species) responded to the different management types. On average, plant species richness was 13% higher in age-class than in unmanaged forests, and did not differ between deciduous age-class and selection forests. In age-class forests of the Schwäbische Alb and Hainich-Dün, coniferous stands had higher species richness than deciduous stands. Among age-class forests, older stands with large quantities of standing biomass were slightly poorer in shrub and light-demanding herb species than younger stands. Among deciduous forests, the richness of herbaceous forest species was generally lower in unmanaged than in managed forests, and it was even 20% lower in unmanaged than in selection forests in Hainich-Dün. Overall, these findings show that disturbances by management generally increase plant species richness. This suggests that total plant species richness is not suited as an indicator for the conservation status of forests, but rather indicates disturbances.     

Characteristics in coarse woody debris mediated by forest developmental stage and latest disturbances in a natural secondary forest of Pinus tabulaeformis

Coarse woody debris (CWD) characteristics are expected to reflect forest stand features. Few studies evaluated logging-induced stand characteristics of secondary coniferous forests by quantifying the quality and quantity in CWD. After selective logging, the form of secondary forest of Pinus tabulaeformis in the Qinling Mountains is inferior and the regeneration is poor. We measured the CWD characteristics of the forest which had an average CWD biomass amount of 12.56 t hm⁻², and was predominated by abundant logs (65.68%), followed by snags (33.13%). The CWD biomass of P. tabulaeformis and Toxicodendron vernicifluum was significantly higher than that of other species, which took up 85.51% of the total. Although there was no significant difference among different diameter sizes (P > 0.05), the CWD biomass of diameter 30–40 cm occupied 46.26% of the total (5.81 t hm⁻²). Similarly, the CWD biomass of decay class I and II accounted for 39.89% (5.01 t hm⁻²) and 33.04% (4.15 t hm⁻²) of the total CWD biomass respectively, despite no significant difference among those 5 decay classes (P > 0.05). The results indicated that the combination of young forest developmental stage caused by past selective logging and natural and anthropogenic disturbances such as strong wind, tapping lacquer, firewood collection, and illegal tree felling played a crucial role in distribution characteristics of CWD in this secondary forest of P. tabulaeformis.     

Tree allometry and improved estimation of carbon stocks and balance in tropical forests

Tropical forests hold large stores of carbon, yet uncertainty remains regarding their quantitative contribution to the global carbon cycle. One approach to quantifying carbon biomass stores consists in inferring changes from long-term forest inventory plots. Regression models are used to convert inventory data into an estimate of aboveground biomass (AGB). We provide a critical reassessment of the quality and the robustness of these models across tropical forest types, using a large dataset of 2,410 trees ≥ 5 cm diameter, directly harvested in 27 study sites across the tropics. Proportional relationships between aboveground biomass and the product of wood density, trunk cross-sectional area, and total height are constructed. We also develop a regression model involving wood density and stem diameter only. Our models were tested for secondary and old-growth forests, for dry, moist and wet forests, for lowland and montane forests, and for mangrove forests. The most important predictors of AGB of a tree were, in decreasing order of importance, its trunk diameter, wood specific gravity, total height, and forest type (dry, moist, or wet). Overestimates prevailed, giving a bias of 0.5–6.5% when errors were averaged across all stands. Our regression models can be used reliably to predict aboveground tree biomass across a broad range of tropical forests. Because they are based on an unprecedented dataset, these models should improve the quality of tropical biomass estimates, and bring consensus about the contribution of the tropical forest biome and tropical deforestation to the global carbon cycle. Electronic Supplementary Material Supplementary material is available for this article at     

Plant Biomass,Nutrient Concentration and Nutrient Storage in a Tropical Dry Forest in the South–west of Madagascar

Plant biomass, mineral composition and the amounts of nutrients in the different fractions of the vegetation were determined for a dense dry deciduous forest growing on light red sands in south-western Madagascar. Complete harvesting and soil coring were used to determine the above- and below-ground biomass respectively. The above-ground biomass, weighing 118 t ha⁻¹ (dry matter), was mostly (96%) made up of phanerophytes (woody trees and shrubs >25 cm tall). Dead material (litter and dead wood on the soil surface) represented 13.8 t ha⁻¹. These results fit well into the range of values reported for other tropical ecosystems. The below-ground biomass was 17.8 t ha⁻¹ giving a root/shoot ratio of 0.15. Rooting is superficial. The nutrient concentration in this dry forest on light reddish-brown sands is, as in other dry forests, considerably higher than that usually found for humid forests. Calcium is the most abundant element. The plant biomass Ca/K ratio is much higher than that of humid tropical forests. In spite of its high originality, this Madagascan dry forest has the same behaviour as other dry forests of the world.     

Biometric-based estimation of net ecosystem production in a mature Japanese cedar (<Emphasis Type="Italic">Cryptomeria japonica</Emphasis>) plantation beneath a flux tower

Quantification of carbon budgets and cycling in Japanese cedar (Cryptomeria japonica D. Don) plantations is essential for understanding forest functions in Japan because these plantations occupy about 20% of the total forested area. We conducted a biometric estimate of net ecosystem production (NEP) in a mature Japanese cedar plantation beneath a flux tower over a 4-year period. Net primary production (NPP) was 7.9 Mg C ha⁻¹ year⁻¹ and consisted mainly of tree biomass increment and aboveground litter production. Respiration was calculated as 6.8 (soil) and 3.3 (root) Mg C ha⁻¹ year⁻¹. Thus, NEP in the plantation was 4.3 Mg C ha⁻¹ year⁻¹. In agreement with the tower-based flux findings, this result suggests that the Japanese cedar plantation was a strong carbon sink. The biometric-based NEP was higher among most other types of Japanese forests studied. Carbon sequestration in the mature plantation was characterized by a larger increment in tree biomass and lower mortality than in natural forests. Land-use change from natural forest to Japanese cedar plantation might, therefore, stimulate carbon sequestration and change the carbon allocation of NPP from an increment in coarse woody debris to an increase in tree biomass.     

Diversity and spatial clustering of shade trees affect cacao yield and pathogen pressure in Costa Rican agroforests

The importance of the spatial organisation of individuals in explaining species coexistence within a community is widely recognised. However, few analyses of spatial structure have been performed on tropical agroforests.The main objective of this study was to highlight the links between spatial organisation of shade trees on the one hand, and shade tree species richness and cacao yield on the other, using data from 29 cacao agroforests in Costa Rica.A method of spatial statistics, Ripley's K-function, was used to analyse the spatial organisation of shade and cacao trees in the study plots. For each stand, the X and Y coordinates of ≥2.5-m-tall trees were recorded. In each plot we also assessed shade tree species richness and cacao yield (with total number of pods = number of pods damaged by frosty pod rot + number of healthy pods).Three types of stands were identified: the first was characterised by significant clustering of shade trees, the highest shade tree species richness (S = 6), and the highest number of damaged pods (139 pods ha^?1 year^?1). The second type was characterised by random spatial organisation of shade trees. The third type showed a trend towards regular organisation. Species richness of shade trees did not differ significantly between the last two types (S = 4 for both), nor did the number of damaged pods (56 pods ha^?1 year^?1 and 67 pods ha^?1 year^?1 respectively).Although the trends were not statistically significant for all the variables in our data set, the clustered spatial structure appears to favour a synergy between environmental (tree species richness), and provisioning (cacao production) services.     

Potential of dendrochronology in assessing carbon sequestration rates of Vitellaria paradoxa in southern Mali,West Africa

This study aimed to investigate the applicability of dendrochronology for assessing the growth dynamics and response to climate variability and to estimate the aboveground carbon stock and carbon sequestration potential of Vitellaria paradoxa in southern Mali. Twenty stem disks were collected from three land-use types (parklands, fallows and protected areas) in Koutiala and Yanfolila districts. We combined a standard dendrochronological approach with biomass allometric equations to estimate the growth and carbon stocks. The results showed that V. paradoxa forms distinct growth ring boundaries but most of the disks from parklands did not successfully cross-date due to management operations like pruning. The tree-ring width showed a significant standardized coefficient of regression with rainfall (r² = 0.66, p < 0.001) but insignificant correlation with temperature. One-way analysis of variance showed no significant difference (p > 0.05) for C-sequestration as well as for carbon stocks in aboveground biomass for both land-use types and sites. Mean values of the amount of C-sequestered in Yanfolila were 0.112 ± 0.0.065 Mg C ha⁻¹ yr⁻¹ in parklands, 0.075 ± 0.018 Mg C ha⁻¹ yr⁻¹ in fallows and 0.064 ± 0.028 Mg C ha⁻¹ yr⁻¹ in protected areas. In Koutiala, the values were 0.068±0.020 Mg C ha⁻¹ yr⁻¹ in the parklands and 0.053 ± 0.017 Mg C ha⁻¹ yr⁻¹ in the fallows. These results clearly indicate that dendrochronology can be applied to assess growth and carbon sequestration potential of V. paradoxa. These results also suggest that climate change could affect the growth and carbon sequestration potential of V. paradoxa. Given the limited size of our sample, figures on the amount of carbon are indicative calling for applying the tested approaches to larger samples and also to other tree species in West Africa.     

Soil and tree ring chemistry of Pinus banksiana and Populus tremuloides stands as indicators of changes in atmospheric environments in the oil sands region of Alberta,Canada

The impact of chronic air pollution such as increased CO₂ and NO_x emissions on forest ecosystems in the Athabasca oil sands region in Alberta, Canada, was investigated in Pinus banksiana (jack pine) and Populus tremuloides (trembling aspen, aspen) stands in two watersheds (NE7 and SM8) located at different distances from the main emission sources of oil sands mining and upgrading facilities, using δ¹³C, δ¹⁵N, and Ca/Al of soil and tree ring samples as indicators. Watershed NE7 was exposed to greater amounts of acid deposition due to its closeness to the mining and upgrading area. The δ¹⁵N in the forest floor was lower (p < 0.05) in NE7 (ranged from −1.42 to −0.87‰) than in SM8 (−0.54 to 1.43‰), implying a greater amount of recent deposition of ¹⁵N-depleted N in NE7. Tree ring δ¹³C gradually decreased over time for both tree species/watersheds, indicating the influence of ¹³C-depleted CO₂ emitted from industrial sources. Tree ring N concentration and δ¹⁵N were not different between watersheds and did not significantly change with time. Interestingly, however, the difference between watersheds (NE7–SM8) that is expressed as Diff_N (for N) increased with concomitant decreases in Diff_δ¹⁵N over time, implying greater increases in ¹⁵N-depleted N input in NE7 than in SM8. Such trends were stronger in aspen stands (R² = 0.64 and p < 0.001 for Diff_N and R² = 0.44 and p < 0.01 for Diff_δ¹⁵N between 1964 and 2009) than in jack pine stands. We conclude that δ¹⁵N in the forest floor and differences in N and δ¹⁵N of tree rings between watersheds are useful indicators reflecting the impact of spatial variations of air pollution on forest stands in the Athabasca oil sands region in western Canada.