Wetland maps including open water extent dynamics based on ENVISAT ASAR WS for Siberia, 2007 and 2008, links to GeoTIFFs

Wetlands store large amounts of carbon, and depending on their status and type, they release specific amounts of methane gas to the atmosphere. The connection between wetland type and methane emission has been investigated in various studies and utilized in climate change monitoring and modelling. For improved estimation of methane emissions, land surface models require information such as the wetland fraction and its dynamics over large areas. Existing datasets of wetland dynamics present the total amount of wetland (fraction) for each model grid cell, but do not discriminate the different wetland types like permanent lakes, periodically inundated areas or peatlands. Wetland types differently influence methane fluxes and thus their contribution to the total wetland fraction should be quantified. Especially wetlands of permafrost regions are expected to have a strong impact on future climate due to soil thawing. In this study ENIVSAT ASAR Wide Swath data was tested for operational monitoring of the distribution of areas with a long-term SW near 1 (hSW) in northern Russia (SW = degree of saturation with water, 1 = saturated), which is a specific characteristic of peatlands. For the whole northern Russia, areas with hSW were delineated and discriminated from dynamic and open water bodies for the years 2007 and 2008. The area identified with this method amounts to approximately 300,000 km**2 in northern Siberia in 2007. It overlaps with zones of high carbon storage. Comparison with a range of related datasets (static and dynamic) showed that hSW represents not only peatlands but also temporary wetlands associated with post-forest fire conditions in permafrost regions. Annual long-term monitoring of change in boreal and tundra environments is possible with the presented approach. Sentinel-1, the successor of ENVISAT ASAR, will provide data that may allow continuous monitoring of these wetland dynamics in the future complementing global observations of wetland fraction.

ERA-interim reanalysis debiased at FLUXNET sites

(preliminary) Exchanges of carbon, water and energy between the land surface and the atmosphere are monitored by eddy covariance technique at the ecosystem level. Currently, the FLUXNET database contains more than 500 sites registered and up to 250 of them sharing data (Free Fair Use dataset). Many modelling groups use the FLUXNET dataset for evaluating ecosystem model's performances but it requires uninterrupted time series for the meteorological variables used as input. Because original in-situ data often contain gaps, from very short (few hours) up to relatively long (some months), we develop a new and robust method for filling the gaps in meteorological data measured at site level. Our approach has the benefit of making use of continuous data available globally (ERA-interim) and high temporal resolution spanning from 1989 to today. These data are however not measured at site level and for this reason a method to downscale and correct the ERA-interim data is needed. We apply this method on the level 4 data (L4) from the LaThuile collection, freely available after registration under a Fair-Use policy. The performances of the developed method vary across sites and are also function of the meteorological variable. On average overall sites, the bias correction leads to cancel from 10% to 36% of the initial mismatch between in-situ and ERA-interim data, depending of the meteorological variable considered. In comparison to the internal variability of the in-situ data, the root mean square error (RMSE) between the in-situ data and the un-biased ERA-I data remains relatively large (on average overall sites, from 27% to 76% of the standard deviation of in-situ data, depending of the meteorological variable considered). The performance of the method remains low for the Wind Speed field, in particular regarding its capacity to conserve a standard deviation similar to the one measured at FLUXNET stations.

Plant Functional Type (PFT) map over Siberia derived from GlobCover 2005 dataset, link to dataset in NetCDF format

High-latitude ecosystems play an important role in the global carbon cycle and in regulating the climate system and are presently undergoing rapid environmental change. Accurate land cover data sets are required to both document these changes as well as to provide land-surface information for benchmarking and initializing Earth system models. Earth system models also require specific land cover classification systems based on plant functional types (PFTs), rather than species or ecosystems, and so post-processing of existing land cover data is often required. This study compares over Siberia, multiple land cover data sets against one another and with auxiliary data to identify key uncertainties that contribute to variability in PFT classifications that would introduce errors in Earth system modeling. Land cover classification systems from GLC 2000, GlobCover 2005 and 2009, and MODIS collections 5 and 5.1 are first aggregated to a common legend, and then compared to high-resolution land cover classification systems, vegetation continuous fields (MODIS VCFs) and satellite-derived tree heights (to discriminate against sparse, shrub, and forest vegetation). The GlobCover data set, with a lower threshold for tree cover and taller tree heights and a better spatial resolution, tends to have better distributions of tree cover compared to high-resolution data. It has therefore been chosen to build new PFT maps for the ORCHIDEE land surface model at 1 km scale. Compared to the original PFT data set, the new PFT maps based on GlobCover 2005 and an updated cross-walking approach mainly differ in the characterization of forests and degree of tree cover. The partition of grasslands and bare soils now appears more realistic compared with ground truth data. This new vegetation map provides a framework for further development of new PFTs in the ORCHIDEE model like shrubs, lichens and mosses, to represent the water and carbon cycles in northern latitudes better. Updated land cover data sets are critical for improving and maintaining the relevance of Earth system models for assessing climate and human impacts on biogeochemistry and biophysics.

Establishment, expansion and regeneration of aspen (Populus tremula) in the Torneträsk area, subarctic Sweden

In subarctic Sweden, recent decadal colonization and expansion of aspen (Populus tremula L.) were recorded. Over the past 100 years, aspen became c. 16 times more abundant, mainly as a result of increased sexual regeneration. Moreover, aspen now reach tree-size (>2 m) at the alpine treeline, an ecotone that has been dominated by mountain birch (Betula pubescens ssp. czerepanovii) for at least the past 4000 years. We found that sexual regeneration in aspen probably occurred seven times or more within the last century. Whereas sexual regeneration occurred during moist years following a year with an exceptionally high June-July temperature, asexual regeneration was favored by warm and dry summers. Disturbance to the birch forest by cyclic moth population outbreaks was critical in aspen establishment in the subalpine area. At the treeline, aspen colonization was less determined by these moth outbreaks, and was mainly restricted by summer temperature. If summer warming persists, aspen spread may continue in subarctic Sweden, particularly at the treeline. However, changing disturbance regimes, future herbivore population dynamics and the responses of aspen's competitors birch and pine to a changing climate may result in different outcomes.

Maps of subsoil temperature and active layer depth of Yakutian ASSR (Autonomous Soviet Socialist Republic of the Soviet Union)

Based on the map of landscapes and permafrost conditions in Yakutia (Merzlotno-landshaftnaya karta Yakutskoi0 ASSR, Gosgeodeziya SSSR, 1991), rasterized maps of permafrost temperature and active-layer thickness of Yakutia, East Siberia were derived. The mean and standard deviation at 0.5-degree grid cell size are estimated by assigning a probability density function at 0.001-degree spatial resolution. Spatial pattern of both variables are dominated by a climatic gradient from north to south, and by mountains and the soil type distribution. Uncertainties are highest in mountains and in the sporadic permafrost zone in the south. The maps are best suited as a benchmark for land surface models which include a permafrost module.

High-resolution digital elevation model (DEM) of a thermokarst depression in Siberian ice-rich permafrost deposits

This data set provides a high-resolution digital elevation model (DEM) of a thermokarst depression (~7 km²) on ice-complex deposits in the Arctic Lena Delta, Siberia. The DEM based on a geodetic field survey and was used for quantitative land surface analyses and detailed description of the thermokarst depression morphology. Detailed morphometrical analyses, volume calculations, and solar radiation modeling were performed and statistically analyzed by Ulrich et al. (2010) to investigate the asymmetrical thermokarst depression development and directed lake migration previously proposed by Morgenstern et al. (2008). Furthermore, the high-resolution DEM in combination with satellite data allowed detailed analyses of spatial and temporal landscape changes due to thermokarst development (Günther, 2009).

Catchment characteristics and soil and water chemistry of the San Francisco basin and subbasins, Ecuador

In this study, the Mean Transit Time and Mixing Model Analysis methods are combined to unravel the runoff generation process of the San Francisco River basin (73.5 km**2) situated on the Amazonian side of the Cordillera Real in the southernmost Andes of Ecuador. The montane basin is covered with cloud forest, sub-páramo, pasture and ferns. Nested sampling was applied for the collection of streamwater samples and discharge measurements in the main tributaries and outlet of the basin, and for the collection of soil and rock water samples. Weekly to biweekly water grab samples were taken at all stations in the period April 2007-November 2008. Hydrometric data, Mean Transit Time and Mixing Model Analysis allowed preliminary evaluation of the processes controlling the runoff in the San Francisco River basin. Results suggest that flow during dry conditions mainly consists of lateral flow through the C-horizon and cracks in the top weathered bedrock layer, and that all subcatchments have an important contribution of this deep water to runoff, no matter whether pristine or deforested. During normal to low precipitation intensities, when antecedent soil moisture conditions favour water infiltration, vertical flow paths to deeper soil horizons with subsequent lateral subsurface flow contribute most to streamflow. Under wet conditions in forested catchments, streamflow is controlled by near surface lateral flow through the organic horizon. Exceptionally, saturation excess overland flow occurs. By absence of the litter layer in pasture, streamflow under wet conditions originates from the A horizon, and overland flow.

Tropical climate and vegetation simulations during the Heinrich event 1 using an Earth System Model of Intermediate Complexity (EMIC) - the University of Victoria Earth System-Climate Model (UVic ESCM)

Abrupt climate changes from 18 to 15 thousand years before present (kyr BP) associated with Heinrich Event 1 (HE1) had a strong impact on vegetation patterns not only at high latitudes of the Northern Hemisphere, but also in the tropical regions around the Atlantic Ocean. To gain a better understanding of the linkage between high and low latitudes, we used the University of Victoria (UVic) Earth System-Climate Model (ESCM) with dynamical vegetation and land surface components to simulate four scenarios of climate-vegetation interaction: the pre-industrial era, the Last Glacial Maximum (LGM), and a Heinrich-like event with two different climate backgrounds (interglacial and glacial). We calculated mega-biomes from the plant-functional types (PFTs) generated by the model to allow for a direct comparison between model results and palynological vegetation reconstructions. Our calculated mega-biomes for the pre-industrial period and the LGM corresponded well with biome reconstructions of the modern and LGM time slices, respectively, except that our pre-industrial simulation predicted the dominance of grassland in southern Europe and our LGM simulation resulted in more forest cover in tropical and sub-tropical South America. The HE1-like simulation with a glacial climate background produced sea-surface temperature patterns and enhanced inter-hemispheric thermal gradients in accordance with the "bipolar seesaw" hypothesis. We found that the cooling of the Northern Hemisphere caused a southward shift of those PFTs that are indicative of an increased desertification and a retreat of broadleaf forests in West Africa and northern South America. The mega-biomes from our HE1 simulation agreed well with paleovegetation data from tropical Africa and northern South America. Thus, according to our model-data comparison, the reconstructed vegetation changes for the tropical regions around the Atlantic Ocean were physically consistent with the remote effects of a Heinrich event under a glacial climate background.

Global 2009 1km MODIS (MOD35/MOD09) cloud frequency and MOD35 processing path

Identifying cloud interference in satellite-derived data is a critical step toward developing useful remotely sensed products. Most MODIS land products use a combination of the MODIS (MOD35) cloud mask and the 'internal' cloud mask of the surface reflectance product (MOD09) to mask clouds, but there has been little discussion of how these masks differ globally. We calculated global mean cloud frequency for both products, for 2009, and found that inflated proportions of observations were flagged as cloudy in the Collection 5 MOD35 product. These erroneously categorized areas were spatially and environmentally non-random and usually occurred over high-albedo land-cover types (such as grassland and savanna) in several regions around the world. Additionally, we found that spatial variability in the processing path applied in the Collection 5 MOD35 algorithm affects the likelihood of a cloudy observation by up to 20% in some areas. These factors result in abrupt transitions in recorded cloud frequency across landcover and processing-path boundaries impeding their use for fine-scale spatially contiguous modeling applications. We show that together, these artifacts have resulted in significantly decreased and spatially biased data availability for Collection 5 MOD35-derived composite MODIS land products such as land surface temperature (MOD11) and net primary productivity (MOD17). Finally, we compare our results to mean cloud frequency in the new Collection 6 MOD35 product, and find that landcover artifacts have been reduced but not eliminated. Collection 6 thus increases data availability for some regions and land cover types in MOD35-derived products but practitioners need to consider how the remaining artifacts might affect their analysis.