Annual spatial distributions of carbon sources and sinks in Canada's forests at 1 km resolution are computed for the period from 1901 to 1998 using ecosystem models that integrate remote sensing images, gridded climate, soils and forest inventory data. GIS‐based fire scar maps for most regions of Canada are used to develop a remote sensing algorithm for mapping and dating forest burned areas in the 25 yr prior to 1998. These mapped and dated burned areas are used in combination with inventory data to produce a complete image of forest stand age in 1998. Empirical NPP age relationships were used to simulate the annual variations of forest growth and carbon balance in 1 km pixels, each treated as a homogeneous forest stand. Annual CO2 flux data from four sites were used for model validation. Averaged over the period 1990–1998, the carbon source and sink map for Canada's forests show the following features: (i) large spatial variations corresponding to the patchiness of recent fire scars and productive forests and (ii) a general south‐to‐north gradient of decreasing carbon sink strength and increasing source strength. This gradient results mostly from differential effects of temperature increase on growing season length, nutrient mineralization and heterotrophic respiration at different latitudes as well as from uneven nitrogen deposition. The results from the present study are compared with those of two previous studies. The comparison suggests that the overall positive effects of non‐disturbance factors (climate, CO2 and nitrogen) outweighed the effects of increased disturbances in the last two decades, making Canada's forests a carbon sink in the 1980s and 1990s. Comparisons of the modeled results with tower‐based eddy covariance measurements of net ecosystem exchange at four forest stands indicate that the sink values from the present study may be underestimated.

[1999] for the United States suggests that North America (> 15øN), was probably a C sink of 0.2-0.5 Pg C yr '• during 1980s, much less than that of 1.7 Pg C yr" estimated by Fan et al.[1998] using an atmospheric inversion method.

[1] Most climate records and climate change scenarios projected by general circulation models are for atmospheric conditions. However, permafrost distribution as well as ecological and biogeochemical processes at high latitudes is mainly controlled by soil thermal conditions, which may be affected by atmospheric climate change. In this paper, the changes in soil temperature during the twentieth century in Canada were simulated at 0.5°latitude/longitude spatial resolution using a process-based model. The results show that the mean annual soil temperature differed from the mean annual air temperature by À2°to 7°C, with a national average of 2.5°C. Soil temperature generally responded to the forcing of air temperature but in complex ways. The changes in annual mean soil temperature during the twentieth century differed from that of air temperature by À3°to 3°C from place to place, and the difference was more significant in winter and spring. On average, for the whole of Canada the annual mean soil temperature at 20 cm depth increased by 0.6°C, while the annual mean air temperature increased by 1.0°C. Three mechanisms were investigated to explain this differentiation: air temperature change, which altered the thickness and duration of snow cover, thereby altering the response of soil temperature; seasonal differences in changes of air temperature; and changes in precipitation. The first two mechanisms generally buffer the response of soil temperature to changes in air temperature, while the effect of precipitation is significant and varies with time and space. This complex response of soil temperature to changes in air temperature and precipitation would have significant implications for the impacts of climate change.

No abstract

Annual spatial distributions of carbon sources and sinks in Canada's forests at 1 km resolution are computed for the period from 1901 to 1998 using ecosystem models that integrate remote sensing images, gridded climate, soils and forest inventory data. GIS-based fire scar maps for most regions of Canada are used to develop a remote sensing algorithm for mapping and dating forest burned areas in the 25 yr prior to 1998. These mapped and dated burned areas are used in combination with inventory data to produce a complete image of forest stand age in 1998. Empirical NPP age relationships were used to simulate the annual variations of forest growth and carbon balance in 1 km pixels, each treated as a homogeneous forest stand. Annual CO 2 flux data from four sites were used for model validation. Averaged over the period 1990-1998, the carbon source and sink map for Canada's forests show the following features: (i) large spatial variations corresponding to the patchiness of recent fire scars and productive forests and (ii) a general south-to-north gradient of decreasing carbon sink strength and increasing source strength. This gradient results mostly from differential effects of temperature increase on growing season length, nutrient mineralization and heterotrophic respiration at different latitudes as well as from uneven nitrogen deposition. The results from the present study are compared with those of two previous studies. The comparison suggests that the overall positive effects of nondisturbance factors (climate, CO 2 and nitrogen) outweighed the effects of increased disturbances in the last two decades, making Canada's forests a carbon sink in the 1980s and 1990s. Comparisons of the modeled results with tower-based eddy covariance measurements of net ecosystem exchange at four forest stands indicate that the sink values from the present study may be underestimated.