Previous estimates of land-atmosphere interaction (the impact of soil moisture on precipitation) have been limited by a lack of observational data and by the model dependence of computational estimates. To counter the second limitation, a dozen climate-modeling groups have recently performed the same highly controlled numerical experiment as part of a coordinated comparison project. This allows a multimodel estimation of the regions on Earth where precipitation is affected by soil moisture anomalies during Northern Hemisphere summer. Potential benefits of this estimation may include improved seasonal rainfall forecasts.

The Global Land–Atmosphere Coupling Experiment (GLACE) is a model intercomparison study focusing on a typically neglected yet critical element of numerical weather and climate modeling: land–atmosphere coupling strength, or the degree to which anomalies in land surface state (e.g., soil moisture) can affect rainfall generation and other atmospheric processes. The 12 AGCM groups participating in GLACE performed a series of simple numerical experiments that allow the objective quantification of this element for boreal summer. The derived coupling strengths vary widely. Some similarity, however, is found in the spatial patterns generated by the models, with enough similarity to pinpoint multimodel “hot spots” of land–atmosphere coupling. For boreal summer, such hot spots for precipitation and temperature are found over large regions of Africa, central North America, and India; a hot spot for temperature is also found over eastern China. The design of the GLACE simulations are described in full detail so that any interested modeling group can repeat them easily and thereby place their model’s coupling strength within the broad range of those documented here.

The 12 weather and climate models participating in the Global Land–Atmosphere Coupling Experiment (GLACE) show both a wide variation in the strength of land–atmosphere coupling and some intriguing commonalities. In this paper, the causes of variations in coupling strength—both the geographic variations within a given model and the model-to-model differences—are addressed. The ability of soil moisture to affect precipitation is examined in two stages, namely, the ability of the soil moisture to affect evaporation, and the ability of evaporation to affect precipitation. Most of the differences between the models and within a given model are found to be associated with the first stage—an evaporation rate that varies strongly and consistently with soil moisture tends to lead to a higher coupling strength. The first-stage differences reflect identifiable differences in model parameterization and model climate. Intermodel differences in the evaporation–precipitation connection, however, also play a key role.

Peatlands store large amounts of soil carbon and freshwater, constituting an important component of the global carbon and hydrologic cycles. Accurate information on the global extent and distribution of peatlands is presently lacking but is needed by Earth system models (ESMs) to simulate the effects of climate change on the global carbon and hydrologic balance. Here, we present Peat-ML, a spatially continuous global map of peatland fractional coverage generated using machine learning (ML) techniques suitable for use as a prescribed geophysical field in an ESM. Inputs to our statistical model follow drivers of peatland formation and include spatially distributed climate, geomorphological and soil data, and remotely sensed vegetation indices. Available maps of peatland fractional coverage for 14 relatively extensive regions were used along with mapped ecoregions of non-peatland areas to train the statistical model. In addition to qualitative comparisons to other maps in the literature, we estimated model error in two ways. The first estimate used the training data in a blocked leave-one-out cross-validation strategy designed to minimize the influence of spatial autocorrelation. That approach yielded an average r² of 0.73 with a root-mean-square error and mean bias error of 9.11 % and −0.36 %, respectively. Our second error estimate was generated by comparing Peat-ML against a high-quality, extensively ground-truthed map generated by Ducks Unlimited Canada for the Canadian Boreal Plains region. This comparison suggests our map to be of comparable quality to mapping products generated through more traditional approaches, at least for boreal peatlands.

Plant functional types (PFTs) are used to represent vegetation distribution in land surface models (LSMs). Previous studies have shown large differences in the geographical distribution of PFTs currently used in various LSMs, which may arise from the differences in the underlying land cover products but also the methods used to map or reclassify land cover data to the PFTs that a given LSM represents. There are large uncertainties associated with existing PFT mapping methods since they are largely based on expert judgement and therefore are subjective. In this study, we propose a new approach to inform the mapping or the cross-walking process using analyses from sub-pixel fractional error matrices, which allows for a quantitative assessment of the fractional composition of the land cover categories in a dataset. We use the Climate Change Initiative (CCI) land cover product produced by the European Space Agency (ESA). Previous work has shown that compared to fine-resolution maps over Canada, the ESA-CCI product provides an improved land cover distribution compared to that from the GLC2000 dataset currently used in the CLASSIC (Canadian Land Surface Scheme Including Biogeochemical Cycles) model. A tree cover fraction dataset and a fine-resolution land cover map over Canada are used to compute the sub-pixel fractional composition of the land cover classes in ESA-CCI, which is then used to create a cross-walking table for mapping the ESA-CCI land cover categories to nine PFTs represented in the CLASSIC model. There are large differences between the new PFT distributions and those currently used in the model. Offline simulations performed with the CLASSIC model using the ESA-CCI-based PFTs show improved winter albedo compared to that based on the GLC2000 dataset. This emphasizes the importance of accurate representation of vegetation distribution for realistic simulation of surface albedo in LSMs. Results in this study suggest that the sub-pixel fractional composition analyses are an effective way to reduce uncertainties in the PFT mapping process and therefore, to some extent, objectify the otherwise subjective process.

Canada's boreal forests and tundra ecosystems are responding to unprecedented climate change with implications for the global carbon (C) cycle and global climate. However, our ability to model the response of Canada's terrestrial ecosystems to climate change is limited and there has been no comprehensive, process‐based assessment of Canada's terrestrial C cycle. We tailor the Canadian Land Surface Scheme Including Biogeochemical Cycles (CLASSIC) to Canada and evaluate its C cycling performance against independent reference data. We utilize skill scores to assess model performance against reference data alongside benchmark scores that quantify the level of agreement between the reference data sets to aid in interpretation. Our results demonstrate CLASSIC's sensitivity to prescribed vegetation cover. They also show that the addition of five region‐specific Plant functional types (PFTs) improves CLASSIC's skill at simulating the Canadian C cycle. CLASSIC performs well when tailored to Canada, falls within the range of the reference data sets, and meets or exceeds the benchmark scores for most C cycling processes. New region‐specific land cover products, well‐informed PFT parameterizations, and more detailed reference data sets will facilitate improvements to the representation of the terrestrial C cycle in regional and global land surface models. Incorporating a parameterization for boreal disturbance processes and explicitly representing peatlands and permafrost soils will improve CLASSIC's future performance in Canada and other boreal regions. This is an important step toward a comprehensive process‐based assessment of Canada's terrestrial C cycle and evaluating Canada's net C balance under climate change.