A numerical phytoclimatic classification of China

The relationship is examined between vegetation and climate using climatic variables collected from 644 meteorological stations located throughout China. Multivariate methods are applied directly to the raw climatic data in order to define climatic clusters and to examine the relationship between the clusters and vegetation types. This approach is based on the concept of multidimensional climatic space defined by the combination of climatic variables. Phytoclimatic classes are defined on the basis of the distribution of vegetation types in climatic clusters and a new phytoclimatic classification of China is proposed. Patterns of climatic changes between neighbouring phytoclimatic classes are described. Two indexes of the influence of climate on vegetation are proposed based on discriminant analysis.     

Climate and vegetation at the Eurosiberian-Mediterranean boundary in the Iberian Peninsula

Abstract. The Northern Iberian Peninsula is dominated by various types of vegetation from deciduous oak and ash to evergreen oak woodlands. A recent vegetation map of Spain portrays vegetation series which are characterized in terms of their phytogeographic region or bioclimatic (altitudinal) belt. The aim of this paper is to determine whether the areas comprised by both phytogeographic regions (Eurosiberian and Mediterranean) in the study area, as established from the phytogeographic characterization of the vegetation, can be discriminated by climatic variables using multivariate methods, and to compare these with other conventional approaches. In addition, bioclimatic (altitudinal) belts and the main vegetation types were tested for discrimination by climatic variables. Conventional climatic criteria as well as discriminant and principal component analysis were applied to climatic data from 205 meteorological stations for which vegetation information had been taken from the vegetation map. Conventional criteria are good predictors of the phytogeographic division (Mediterranean and Eurosiberian regions) in the study area. Results were improved by multiple discriminant analysis based on climatic data of the dry period of the year (June to September). Both regions in the study area can be predicted with over 95 % accuracy. Using the same multivariate procedure and temperature data the bioclimatic (altitudinal) belts of the study area can be predicted with over 90 % accuracy. The main vegetation groups of the study area can also be predicted with over 80 % accuracy. Ordination analysis supported the results of the discriminant analysis. Empirical models have been generated to predict the phytogeographic- and belt character of any station in the area. The significance of the various periods of the year for discriminating regions and belts is evaluated. The responsiveness to climatic events during the year may be region specific. This study confirms the strong relationship between climate and vegetation in the Northern Iberian Peninsula, particularly regarding the Eurosiberian-Mediterranean boundary.     

Ecological response surfaces for North American boreal tree species and their use in forest classification

Abstract. Empirical ecological response surfaces were derived for eight dominant tree species in the boreal forest region of Canada. Stepwise logistic regression was used to model species dominance as a response to five climatic predictor variables. The predictor variables (annual snowfall, degree-days, absolute minimum temperature, annual soil moisture deficit, and actual evapotranspiration summed over the summer months) influence the response of plants more directly than the annual or monthly measures of temperature and precipitation commonly used in response surface modeling. The response surfaces provided estimates of the probability of species dominance across the spatial extent of North America with a high degree of success. Much of the variation in the probability of dominance is apparently related to the species' individualistic response to climatic constraints within different airmass regions. A forest type classification for the Canadian boreal forest region was derived by a cluster analysis based on the probability estimates. Five major forest types were distinguished by the application of a stopping rule. The predicted forest types showed a high degree of geographic correspondence with the distribution of forest types in the actual vegetation mosaic. The distribution of the predicted types also bears a direct relationship to seasonal airmass dynamics in the boreal forest region.     

中国不同植被类型归一化植被指数对气候变化和人类活动的响应

不同植被类型对外界干扰和环境变化的敏感性不同。为厘清中国不同类型植被的动态变化特征及其对外界环境变化的响应,综合利用趋势分析、残差分析和情景模拟方法,在明确2000-2015年间我国不同植被类型归一化植被指数(NDVI)时空变化基础上,对气候变化和人类活动两大驱动要素在不同植被类型NDVI变化中的相对贡献进行了定量评估和归因。研究结果表明:(1)2000-2015年,我国植被NDVI整体呈增加趋势,且其空间占比高达84.1%。其中,森林植被的改善状况最佳,显著增加的面积占到了森林总面积的82.4%;而荒漠植被的改善状况相对较差,仅有22.3%的区域呈显著增加趋势。(2)人类活动在我国植被变化中占主导地位。植被改善区和植被退化区人类活动的相对贡献分别为76.4%和60.0%,且人类活动对植被的影响更多与管理方式而非土地利用类型转变有关。(3)不同类型植被对气候变化和人类活动的响应差异显著。对于植被改善区,除沼泽外,人类活动对各类型植被NDVI变化的贡献率均在70%以上,尤其是对农作物的贡献率最高,达到80.7%;对于植被退化区,人类活动影响较大的植被类型为沼泽和农作物,表明2000-2015年间我国沼泽受到了更强烈人类活动的负面影响。研究有助于增强对不同植被类型对全球变化响应机制的理解,并为促进生态建设和植被恢复工作的有效实施提供科学参考。     

Spatial information analysis of chemotactic trajectories

During bacterial chemotaxis, a cell acquires information about its environment by sampling changes in the local concentration of a chemoattractant, and then uses that information to bias its motion relative to the source of the chemoattractant. The trajectory of a chemotaxing bacteria is thus a spatial manifestation of the information gathered by the cell. Here we show that a recently developed approach for computing spatial information using Fourier coefficient probabilities, the k-space information (kSI), can be used to quantify the information in such trajectories. The kSI is shown to capture expected responses to gradients of a chemoattractant. We then extend the k-space approach by developing an experimental probability distribution (EPD) that is computed from chemotactic trajectories collected under a reference condition. The EPD accounts for connectivity and other constraints that the nature of the trajectories imposes on the k-space computation. The EPD is used to compute the spatial information from any trajectory of interest, relative to the reference condition. The EPD-based spatial information also captures the expected responses to gradients of a chemoattractant, although the results differ in significant ways from the original kSI computation. In addition, the entropy calculated from the EPD provides a useful measure of trajectory space. The methods developed are highly general, and can be applied to a wide range of other trajectory types as well as non-trajectory data.     

Predicting physiognomic vegetation types with climate variables

A quantitative terrestrial vegetation model was produced which consists of:

1. A world classification of important terrestrial plant growth forms (life forms);
2. A set of predictive variables representing the main climatic correlates of these forms; and
3. Empirically obtained hypothetical limiting values defining an ecoclimatic envelope for each plant form (relative to the climatic variables).

The model was applied to a world climatic data-base (1 225 sites) in order to substantiate the hypothesized life-form status of the plant types by accurately predicting their actual world distributions. Particular combinations of forms are interpreted as vegetation formation types by reference to growth-form dominance considerations. Model validation was attempted by comparing predicted and actually occurring vegetation at independent sites on all continents. Prediction accuracy of 85% for individual plant types and 50% for vegetation structure (exact combination of actually occurring dominant forms) suggests that general macroclimatic conditions are much more important than any other factors (such as complex specific interactions) in determining general ecological structure on most sites.     

Zonal and azonal vegetation revisited: How is wetland vegetation distributed across different zonobiomes

In ‘zonal’ vegetation, climatic factors are the main influence on growth and performance and the climate determines the vegetation type completely, which makes this vegetation dominant in the landscape. If vegetation is ‘azonal’ however, local stresses are assumed to have an overwhelming influence on plant performance and climatic influences will be minimal; typically, this vegetation occurs only in small patches in the landscape. In this study I ask whether wetland plant communities, as they are described for South Africa, are evenly distributed among different terrestrial vegetation types, to test whether they are zonal or azonal. Three contingency tables were construed based on the counts of wetland vegetation records, defined on three hierarchical levels (Main Clusters, Community Groups and Community) and their occurrence in the country (at the level of Biome, Bioregion and terrestrial vegetation type). An ‘azonality index’ was calculated as the sum of all Chi‐square values for each wetland vegetation type divided by the total number of records. The overall correlation between hydroperiod and the azonality index was very weak. At the finest level, terrestrial vegetation types were clustered on the basis of having similar combinations of wetland community types. Eighteen different ‘wetland ecoregions’ have been defined, on the basis of wetland vegetation types occurring within them. Instead of regarding wetland vegetation as azonal, it should rather be regarded as ‘intrazonal’, meaning that climate does have an impact but many vegetation types are widespread across climatic regions. The reason why community types in wetlands are widespread is due to the monodominance of a single widespread, often clonal, species. The different wetland ecoregions do not correspond to terrestrial biomes, so it is expected that wetland vegetation responds differently to climate than terrestrial vegetation.     

Climate and vegetation in China II. Distribution of main vegetation types and thermal climate

The authors examined relationships between Kira's warmth index (WI) and four other important thermal indices: the sums of daily mean temperatures above 5°C and 10°C, Thornthwaite's potential evapotranspiration (PE) and Holdridge's annual biotemperature. The thermal records of 671 meteorological stations evenly located all over China were used to make these comparisons. Close correlations were found within the four relationships, and accordingly WI was used to analyse the thermal distributions of the main vegetation types. Vegetation types around the 671 stations were read from a vegetation map with a scale of 1/4000000. Vegetation types at 269 stations corresponded to the natural or seminatural vegetation, and 29 vegetation types were distinguished by arranging the 269 data into the same or similar types. The geographical distribution of these 29 types and the corresponding main climatic features were described. The relations between WI and distribution of these vegetation types were discussed in detail. As a result, WI values (°C month) corresponding to the vegetation zones could be summarized as follows: (1) arctic or alpine vegetation zone: 0–15; (2) boreal or subalpine vegetation zone: 15-(50–55); (3) cool-temperate vegetation zone: (50–55)–(80–90); (4) warm-temperate vegetation zone: (80–90)–(170–180). These values almost coincided with Kira's values. Chinese postgraduate student in Japan sent by the Chinese Government.     

A vegetation-ecological approach to the classification and evaluation of potential natural vegetation of the Fagetea crenatae region in Tohoku (northern Honshu), Japan

A new vegetation-ecological approach is proposed for classification and evaluation of vegetation zones by means of phytosociological landscape analysis, based on the potential natural vegetation. The study area is the “Fagetea crenatae region” of the cool-temperate zone of Tohoku (northern Honshu) and the northern parts of Kanto. The area was divided into 953 geographic quadrats on a base map at a scale of 1 ∶ 500000. Based on climax complexes of vegetation in each quadrat, 55 community sub-groups were distinguished as basic units of community complex and vegetation landscapes. The community sub-groups were then grouped into 17 larger community groups by the phytosociological table method. As a result, three phytogeographic vegetation zones (Japan Sea side, inland areas and Pacific side) were classified. For each of these community sub-groups, five geographical and climatic variables (average altitude, mean annual temperature, Kira's warmth index, annual precipitation and mean annual maximum snow depth) were averaged, and the community sub-groups in the same community group, which resembled each other ecologically, were assembled into 28 clusters. The clusters were combined into 11 ecological groups by means of Pearson's similarity ratio of geographical and climatic characteristics. By comparing these ecological groups as a vegetation complex, four phytogeographic vegetation zones (Japan Sea side, inland areas, Pacific side and northern Honshu) corresponding to each potential natural vegetation region with distinct environmental characteristics, were newly classified.     

用于全球变化研究的中国植物功能型划分

 植物功能型(Plant functional types, PFTs)作为沟通植物的结构和功能与生态系统属性的桥梁,随着全球变化与植被的关系研究的深入而受到广泛重视。植物功能型的划分依赖于研究的背景、尺度和要解决的问题。为了区域尺度全球变化研究的需要,该文提出了一个基于植物关键特征的植物功能型划分方法。该方法首先选择了6项植物特征,包括3项冠层特征：木本-草本、常绿-落叶和针叶-阔叶,以及3项生理特征：光合途径(C3 / C4)、植物的水分需求和热量需求,作为划分植物功能型的关键特征;然后,先根据植物冠层特征划分得到5个基本类型,再根据水分和热量条件进行详细划分,得到29种备选类型;需要时,再根据研究目的从这29种备选类型中选择所需类型。根据这个方法,在充分考虑了我国季风气候条件下特有的水热配置和高海拔环境对植物的形态和功能特征影响的基础上,从备选类型中选择了一套适合中国气候和植被特征的植物功能型体系。这套体系包括18类植物功能型,其中含7类‘树’功能型、6类‘灌木’功能型和5类‘草’功能型,另根据需要设置2类‘裸地’功能型。并且根据植物的生理生态特征和中国植被的地理分布确定了用于限制植物功能型分布的气候因子,这些气候因子包括绝对最低温度、最暖月平均温度、有效积温、年最热月平均温和最冷月平均温之差、湿润指数、年均降水量。应用表明,这套植物功能型可用于模拟我国植被在当前气候条件下的分布。该研究为发展适于我国的植被模型和区域气候模型、评估全球变化对我国植被的影响及植被变化对气候的反馈作用提供依据与参数。     

Interactive effects of vegetation type and elevation on aboveground and belowground properties in a subarctic tundra

An improved knowledge of how contrasting types of plant communities and their associated soil biota differ in their responses to climatic variables is important for better understanding the future impacts of climate change on terrestrial ecosystems. Elevational gradients serve as powerful study systems for answering questions on how ecological processes can be affected by changes in temperature and associated climatic variables. In this study, we evaluated how plant and soil microbial communities, and abiotic soil properties, change with increasing elevation in subarctic tundra in northern Sweden, for each of two dominant but highly contrasting vegetation types, namely heath (dominated by woody dwarf shrubs) and meadow (dominated by herbaceous species). To achieve this, we measured plant community characteristics, microbial community properties and several soil abiotic properties for both vegetation types across an elevation gradient of 500 to 1000 m. We found that the two vegetation types differed not only in several above‐ and belowground properties, but also in how these properties responded to elevation, pointing to important interactive effects between vegetation type and elevation. Specifically, for the heath, available soil nitrogen and phosphorus decreased with elevation whereas fungal dominance increased, while for the meadow, idiosyncratic responses to elevation for these variables were found. These differences in belowground responses to elevation among vegetation types were linked to shifts in the species and functional group composition of the vegetation. Our results highlight that these two dominant vegetation types in subarctic tundra differ greatly not only in fundamental aboveground and belowground properties, but also in how these properties respond to elevation and are therefore likely to be influenced by temperature. As such they highlight that vegetation type, and the soil abiotic properties that determine this, may serve as powerful determinants of how both aboveground and belowground properties respond to strong environmental gradients.     

Representation of tundra vegetation by pollen in lake sediments of northern Alaska

Aim To understand better the representation of arctic tundra vegetation by pollen data, we analysed pollen assemblages and pollen accumulation rates (PARs) in the surface sediments of lakes. Location Modern sediment samples were collected from seventy‐eight lakes located in the Arctic Foothills and Arctic Coastal Plain regions of northern Alaska. Methods For seventy of the lakes, we analysed pollen and spores in the upper 2 cm of the sediment and calculated the relative abundance of each taxon (pollen percentages). For eleven of the lakes, we used ²¹⁰Pb analysis to determine sediment accumulation rates, and analysed pollen in the upper 10–15 cm of the sediment to estimate modern PARs. Using a detailed land‐cover map of northern Alaska, we assigned each study site to one of five tundra types: moist dwarf‐shrub tussock‐graminoid tundra (DST), moist graminoid prostrate‐shrub tundra (PST) (coastal and inland types), low‐shrub tundra (LST) and wet graminoid tundra (WGT). Results Mapped pollen percentages and multivariate comparison of the pollen data using discriminant analysis show that pollen assemblages vary along the main north–south vegetational and climatic gradients. On the Arctic Coastal Plain where climate is cold and dry, graminoid‐dominated PST and WGT sites were characterized by high percentages of Cyperaceae and Poaceae pollen. In the Arctic Foothills where climate is warmer and wetter, shrub‐dominated DST, PST and LST were characterized by high percentages of Alnus and Betula pollen. Small‐scale variations in tundra vegetation related to edaphic variability are also represented by the pollen data. Discriminant analysis demonstrated that DST sites could be distinguished from foothills PST sites based on their higher percentages of Ericales and Rubus chamaemorus pollen, and coastal PST sites could be distinguished from WGT sites based on their higher percentages of Artemisia. PARs appear to reflect variations in overall vegetation cover, although the small number of samples limits our understanding of these patterns. For coastal sites, PARs were higher for PST than WGT, whereas in the Arctic Foothills, PARs were highest in LST, intermediate in DST, and lowest in PST. Main conclusion Modern pollen data from northern Alaska reflect patterns of tundra vegetation related to both regional‐scale climatic gradients and landscape‐scale edaphic heterogeneity.     

Arable weed vegetation of the northeastern part of the Czech Republic: effects of environmental factors on species composition

Factors determining changes in species composition of arable field weed vegetation in the northeastern part of the Czech Republic were studied. Gradsect sampling, i.e. a priori stratified selection of sampling sites, was used for the field research. Using this method, a data set of 174 vegetation plots, covering a whole range of basic environmental characteristics in the study area, was compiled in 2001–2003. A set of environmental variables (altitude, annual precipitation, mean annual temperature, soil type, soil pH and crop type) together with date of sampling was obtained for each plot. Ordination methods were used to determine the effects of variables on arable weed composition. For each variable, the gross and net effect on weed species composition were calculated. All variables considered in this study had a significant effect on weed species composition and explained 7.25% of the total variation in species data. Major changes in weed species composition in the study area were associated with different crop types. The second most important gradient in the variability of weed vegetation in the study area was associated with altitudinal and climatic changes followed by seasonal changes and different soil types and pH. Our results show that on a regional scale, the relative importance of different crop types and their associated management on changes in arable weed species composition is higher than the relative importance of climatic variables. The relative importance of climatic variables decreases with their decreasing length of gradient.     

The potential contribution of vegetation ecology to biodiversity research

The contribution of vegetation ecology to the study of biodiversity depends on better communication between the different research paradigms in ecology. Recent developments in vegetation theory and associated statistical modelling techniques are reviewed for their relevance to biodiversity. Species composition and collective properties such as species richness vary as a continuum in a multi-dimensional environmental space; a concept which needs to be incorporated into biodiversity studies. Different kinds of environmental gradients can be recognised and species responses to them vary. Species response curves of eucalypts to an environmental gradient of mean annual temperature have been shown to exhibit a particular pattern of skewed response curves. Generalised linear modelling (GLM) and generalised additive modelling (GAM) techniques are important tools for biodiversity studies. They have successfully distinguished the contribution of environmental (climatic) and spatial (history and species dispersal ability) variables in determining forest tree composition in New Zealand. Species richness studies are examined at global, regional and local scales. At all scales, direct and resource environmental gradients need to be incorporated into the analysis rather than indirect gradients e.g. latitude which have no direct physiological influence on biota. Evidence indicates that species richness at the regional scale is sensitive to environment, confounding current studies on local/regional species richness relationships. Plant community experiments require designs based on environmental gradients rather than dependent biological properties such as productivity or species richness to avoid confounding the biotic components. Neglect of climatic and other environmental gradients and the concentration on the collective properties of species assemblages has limited recent biodiversity studies. Conservation evaluation could benefit from greater use of the continuum concepts and statistical modelling techniques of vegetation ecology. The future development of ecology will depend on testing the different assumptions of competing research paradigms and a more inclusive synthesis of ecological theory.     

Spontaneous vegetation succession in disused gravel‐sand pits: Role of local site and landscape factors

Questions : What is the variability of succession over a large geographical area? What is the relative importance of (1) local site factors and (2) landscape factors in determining spontaneous vegetation succession? Location : Various regions of the Czech Republic, Central Europe. The regions represent two categories characterized by agrarian lowlands, with a relatively warm and dry climate, and predominant woodland uplands with a relatively cold and wet climate. Methods : Gravel‐sand pits ranged in age from 1–75 years since abandonment. Three types of sites were distinguished: dry, wet and hydric in shallow flooded sites. Vegetation relevés were recorded with species cover (%) visually estimated using the space‐for‐time substitution approach. Local site factors, such as water table and soil characteristics, and landscape characteristics, namely climatic parameters, presence of nearby (semi‐) natural plant communities and main land cover categories in the wider surroundings, were evaluated. Results : Ordination analyses showed that water table was the most important local site factor influencing the course of spontaneous vegetation succession. Succession was further significantly influenced by soil texture, pH, macroclimate, the presence of some nearby (semi‐) natural communities and some land cover categories in the wider surroundings. Spontaneous vegetation succession led to the formation of either shrubby grassland, deciduous woodland, alder and willow carrs, and tall sedge or reed and Typha beds in later stages depending predominantly on the site moisture conditions. Conclusions: Although the water table was the most influential on the course of vegetation succession, the landscape factors together explained more vegetation variability (44%) than local site factors (23%).     

Determinants of bird species richness: role of climate and vegetation structure at a regional scale

We examined the respective roles of climate and vegetation structure on geographical variation in bird species richness. The Province of Buenos Aires (central-eastern Argentina) was divided into 146 squares of 50 km on a side. For each square we evaluated the number of bird species, the value of thirteen climatic variables, and the value of a vegetation strata index. The climatic matrix was analyzed by Principal Component Analysis (PCA), and the first factors resulting from PCA were considered as multifactorial climatic gradients. Simple and Partial Correlation Analysis among bird species richness, vegetation strata, and the first two factors derived from PCA (65% of total variation) indicated that bird richness distribution was determined by the availability of vegetation strata, associated with different vegetation types that, at the same time, were influenced by the climatic conditions summarized in the first climatic factor (a gradient of precipitation, relative humidity, annual termical amplitude, and frost occurrence). This relationships reflect the complexity of factors that can act directly as well as indirectly on the geographical patterns in species richness. Also, we evaluated the importance of study scale comparing our results with previous studies at macrogeographic and local scales, found out that the vegetation structure was the principal determinant of bird species richness at this three geographical scales.     

中国热带森林植被类型研究历史和划分探讨

热带森林是我国森林植被的重要组成部分,明确其森林植被类型分类对于生物多样性维持机制研究和保护管理等都具有重要意义。该文以中国热带森林植被分类研究中存在的问题为出发点,通过阐述我国各省区植被类型分类的研究历史,在综合考虑多种生物和非生物影响因素的基础上,提出一个新的热带森林植被类型分类框架以供探讨。结果表明:(1)尽管针对我国热带森林的分布范围和群落特征等都已开展了诸多研究,但对我国热带森林植被类型的划分依据和分类体系仍存在争议。(2)尽管我国的热带森林都处于季风气候区带内,但许多热带地区的森林植被类型并不只是受季风影响,而是气候带、关键气候因子、地形、土壤反馈和物种适应等多种因素共同作用的结果。(3)我国的热带森林植被包括5个植被型,即非典型性热带雨林、热带季雨林、热带山地雨林、热带山顶苔藓矮林(热带云雾林)和热带针叶林,其中热带季雨林植被型包含4种植被亚型[热带落叶季雨林、热带半落叶(半常绿)季雨林、热带常绿季雨林和热带石灰岩(石山)季雨林]。(4)阐明了上述热带森林植被型和植被亚型在我国各省区的分布情况,并提出未来有必要对人工恢复后的热带森林进行评估和植被类型划分。综上所述,该文提出一个新的热带森林植被类型分类框架,以期为今后基于不同地区开展热带森林比较研究提供参考。     

Modern vegetation–pollen–climate relationships for the Pampa grasslands of Argentina

Aim To analyse the relationships between potential natural vegetation, pollen and climate in order to improve the interpretation of fossil pollen records and provide the background for future quantitative palaeoclimatic reconstructions. Location Pampa grasslands of Argentina, between 33–41° S and 56–67° W. Methods Modern pollen data were obtained from a pollen data base developed by the Grupo de Investigación de Paleoecología y Palinología, Universidad Nacional de Mar del Plata, Argentina (143 surface samples and 17 pollen types). Analysis of pollen and climate data involved multivariate statistics (cluster analysis and principal components analysis), scatter diagrams, Pearson’s correlation and isopoll mapping. Results Vegetation patterns at regional scales (grasslands and xerophytic woodlands) and local scales (edaphic communities) were identified by cluster analysis of pollen surface samples. The main climatic variables that appear to constrain the vegetation distribution and abundance of taxa are mean annual precipitation, annual effective precipitation and summer temperature. Individual pollen types such as Chenopodiaceae, Apiaceae, Cyperaceae, Prosopis, Schinus, Condalia microphylla and other xerophytic taxa are good indicators of moisture regime. Many pollen types are significantly correlated with summer temperature. The modern vegetation–pollen–climate relationships vary in a broadly predictable manner, supporting the contention that fossil pollen assemblages can be related to particular climatic characteristics. Main conclusions An expanded suite of modern analogues facilitated new insights into vegetation–pollen–climate relationships at the regional scale in Pampa grasslands. Relationships between individual pollen types and climate are appraised at a regional scale and new modern analogues are presented. The results provide the basis for improved vegetation and climate reconstruction from fossil records of the study area.     

Mapping forest vegetation patterns in an Atlantic-Mediterranean transitional area by integration of ordination and geostatistical techniques

Aims Forest vegetation variability may be explained by the complex interplay among several spatial structuring factors, including climate and topography. We modelled the spatial variability of forest vegetation assemblages and significant environmental variables along a complex environmental gradient or coenocline to produce a detailed cartographic database portraying the distribution of forests along it.Methods We combined an analysis of ordination coenoclines with kriging over 772 field data plots from the third Spanish National Forest Inventory in an Atlantic–Mediterranean transitional area (northern Spain).Important findings The best fitted empirical semivariogram revealed a strong spatial structure of forest species composition along the complex environmental gradient considered (the climatic–topographic gradient from north to south). The steady and gradual increase of semivariance with a marked lag distance indicates a gradual turnover of forest assemblages according to the climatic–topographic variations (regional or local). Two changes in the slope of the semivariogram suggest the existence of two different scales of spatial variation. The interpolation map by Kriging of forest vegetation assemblages along the main coenocline shows a clear spatial distribution pattern of trees and shrubs in accordance with the spatial variation of significant environmental variables. We concluded that the multivariate geostatistical approach is a suitable technique for spatial analysis of forest systems employing data from national forest inventories based on a regular network of field plots. The development of an assortment of maps describing changes in vegetation assemblages and variation in environmental variables is expected to be a suitable tool for an integrated forest management and planning.