Predictive vegetation mapping in the Mediterranean context: Considerations and methodological issues

The need to map vegetation communities over large areas for nature conservation and to predict the impact of environmental change on vegetation distributions, has stimulated the development of techniques for predictive vegetation mapping. Predictive vegetation studies start with the development of a model relating vegetation units and mapped physical data, followed by the application of that model to a geographic database and over a wide range of spatial scales. This field is particularly important for identifying sites for rare and endangered species and locations of high biodiversity such as many areas of the Mediterranean Basin. The potential of the approach is illustrated with a mapping exercise in the alti-meditterranean zone of Lefka Ori in Crete. The study established the nature of the relationship between vegetation communities and physical data including altitude, slope and geomorphology. In this way the knowledge of community distribution was improved enabling a GIS-based model capable of predicting community distribution to be constructed. The paper describes the development of the spatial model and the methodological problems of predictive mapping for monitoring Mediterranean ecosystems. The paper concludes with a discussion of the role of predictive vegetation mapping and other spatial techniques, such as fuzzy mapping and geostatistics, for improving our understanding of the dynamics of Mediterranean ecosystems and for practical management in a region that is under increasing pressure from human impact.

A new Snow-SVAT to simulate the accumulation and ablation of seasonal snow cover beneath a forest canopy

A new snow-soil-vegetation-atmosphere transfer (Snow-SVAT) scheme, which simulates the accumulation and ablation of the snow cover beneath a forest canopy, is presented. The model was formulated by coupling a canopy optical and thermal radiation model to a physically-based multi-layer snow model. This canopy radiation model is physically-based yet requires few parameters, so can be used when extensive in-situ field measurements are not available. Other forest effects such as the reduction of wind speed, interception of snow on the canopy and the deposition of litter were incorporated within this combined model, SNOWCAN, which was tested with data taken as part of the Boreal Ecosystem-Atmosphere Study (BOREAS) international collaborative experiment. Snow depths beneath four different canopy types and at an open site were simulated. Agreement between observed and simulated snow depths was generally good, with correlation coefficients ranging between r^2=0.94 and r^2=0.98 for all sites where automatic measurements were available. However, the simulated date of total snowpack ablation generally occurred later than the observed date. A comparison between simulated solar radiation and limited measurements of sub-canopy radiation at one site indicates that the model simulates the sub-canopy downwelling solar radiation early in the season to within measurement uncertainty.

Factors affecting distribution of vegetation types on abandoned cropland in the hilly-gullied Loess Plateau region of China

A study was conducted in the forest-steppe region of the Loess Plateau to provide insight into the factors affecting the process of vegetation establishment, and to provide recommendations for the selection of indigenous species in order to speed up the succession process and to allow the establishment of vegetation more resistant to soil erosion. Four distinctive vegetation types were identified, and their distribution was affected not only by the time since abandonment but also by other environmental factors, mainly soil water and total P in the upper soil layers. One of the vegetation types, dominated by Artemisia scoparia, formed the early successional stage after abandonment while the other three types formed later successional stages with their distribution determined by the soil water content and total P. It can be concluded that the selection of appropriate species for introduction to accelerate succession should be determined by the local conditions and especially the total P concentration and soil water content.

An assessment of the MODIS collection 5 leaf area index product for a region of mixed coniferous forest

Canopy leaf area index (LAI), defined as the single-sided leaf area per unit ground area, is a quantitative measure of canopy foliar area. LAI is a controlling biophysical property of vegetation function, and quantifying LAI is thus vital for understanding energy, carbon and water fluxes between the land surface and the atmosphere. LAI is routinely available from Earth Observation (EO) instruments such as MODIS. However EO-derived estimates of LAI require validation before they are utilised by the ecosystem modelling community. Previous validation work on the MODIS collection 4 (c4) product suggested considerable error especially in forested biomes, and as a result significant modification of the MODIS LAI algorithm has been made for the most recent collection 5 (c5). As a result of these changes the current MODIS LAI product has not been widely validated. We present a validation of the MODIS c5 LAI product over a 121 km2 area of mixed coniferous forest in Oregon, USA, based on detailed ground measurements which we have upscaled using high resolution EO data. Our analysis suggests that c5 shows a much more realistic temporal LAI dynamic over c4 values for the site we examined. We find improved spatial consistency between the MODIS c5 LAI product and upscaled in situ measurements. However results also suggest that the c5 LAI product underestimates the upper range of upscaled in situ LAI measurements.

Vegetation and environmental dynamics in the southern Black Sea region since 18 kyr BP derived from the marine core 22-GC3

Sediments from the Black Sea, a region historically dominated by forests and steppe landscapes, are a valuable source of detailed information on the changes in regional terrestrial and aquatic environments at decadal to millennial scales. Here we present multi-proxy environmental records (pollen, dinoflagellate cysts, Ca, Ti and oxygen isotope data) from the uppermost 305 cm of the core 22-GC3 (42°13.53′N, 36°29.55′E) collected from a water depth of 838 m in the southern part of the Black Sea in 2007. The records span the last ~ 18 kyr (all ages are given in cal kyr BP). The pollen data reveal the dominance of the Artemisia-steppe in the region, suggesting rather dry/cold environments ~ 18–14.5 kyr BP. Warming/humidity increase during melt-water pulses (~ 16.1–14.5 kyr BP), indicated by δ18O records from the 22-GC3 core sediment and from the Sofular Cave stalagmite, is expressed in more negative δ13C values from the Sofular Cave, usually interpreted as the spreading of C3 plants. The records representing the interstadial complex (~ 14.5–12.9 kyr BP) show an increase in temperature and moisture, indicated by forest development, increased primary productivity and reduced surface run-off, whereas the switch from primary terrigenous to primary authigenic Ca origin occurs ~ 500 yr later. The Younger Dryas cooling is clearly demonstrated by more negative δ13C values from the Sofular Cave and a reduction of pines. The early Holocene (11.7–8.5 kyr BP) interval reveals relatively dry conditions compared to the mostly moist and warm middle Holocene (8.5–5 kyr BP), which is characterized by the establishment of the species-rich warm mixed and temperate deciduous forests in the low elevation belt, temperate deciduous beech-hornbeam forests in the middle and cool conifer forest in upper mountain belt. The border between the early and middle Holocene in the vegetation records coincides with the opening of the Mediterranean corridor at ~ 8.3 kyr BP, as indicated by a marked change in the dinocyst assemblages and in the sediment lithology. Changes in the pollen assemblages indicate a reduction in forest cover after ~ 5 kyr BP, which was likely caused by increased anthropogenic pressure on the regional vegetation.

Land-atmosphere energy exchange in Arctic tundra and boreal forest: available data and feedbacks to climate

This paper summarizes and analyses available data on the surface energy balance of Arctic tundra and boreal forest. The complex interactions between ecosystems and their surface energy balance are also examined, including climatically induced shifts in ecosystem type that might amplify or reduce the effects of potential climatic change. High latitudes are characterized by large annual changes in solar input. Albedo decreases strongly from winter, when the surface is snow-covered, to summer, especially in nonforested regions such as Arctic tundra and boreal wetlands. Evapotranspiration (QE) of high-latitude ecosystems is less than from a freely evaporating surface and decreases late in the season, when soil moisture declines, indicating stomatal control over QE, particularly in evergreen forests. Evergreen conifer forests have a canopy conductance half that of deciduous forests and consequently lower QE and higher sensible heat flux (QH). There is a broad overlap in energy partitioning between Arctic and boreal ecosystems, although Arctic ecosystems and light taiga generally have higher ground heat flux because there is less leaf and stem area to shade the ground surface, and the thermal gradient from the surface to permafrost is steeper. Permafrost creates a strong heat sink in summer that reduces surface temperature and therefore heat flux to the atmosphere. Loss of permafrost would therefore amplify climatic warming. If warming caused an increase in productivity and leaf area, or fire caused a shift from evergreen to deciduous forest, this would increase QE and reduce QH. Potential future shifts in vegetation would have varying climate feedbacks, with largest effects caused by shifts from boreal conifer to shrubland or deciduous forest (or vice versa) and from Arctic coastal to wet tundra. An increase of logging activity in the boreal forests appears to reduce QE by roughly 50% with little change in QH, while the ground heat flux is strongly enhanced.

Sensitivity of Bolivian seasonally-dry tropical forest to precipitation and temperature changes over glacial-interglacial timescales

We used fossil pollen to investigate the response of the eastern Chiquitano seasonally-dry tropical forest (SDTF), lowland Bolivia, to high-amplitude climate change associated with glacial–interglacial cycles. Changes in the structure, composition and diversity of the past vegetation are compared with palaeoclimate data previously reconstructed from the same record, and these results shed light on the biogeographic history of today’s highly disjunct blocks of SDTF across South America. We demonstrate that lower glacial temperatures limited tropical forest in the Chiquitanía region, and suggest that SDTF was absent or restricted at latitudes below 17°S, the proposed location of the majority of the hypothesized ‘Pleistocene dry forest arc’ (PDFA). At 19500 yrs b.p., warming supported the establishment of a floristically-distinct SDTF, which showed little change throughout the glacial–Holocene transition, despite a shift to significantly wetter conditions beginning ca. 12500–12200 yrs b.p. Anadenanthera colubrina, a key SDTF taxon, arrived at 10000 yrs b.p., which coincides with the onset of drought associated with an extended dry season. Lasting until 3000 yrs b.p., Holocene drought caused a floristic shift to more drought-tolerant taxa and a reduction in α-diversity (shown by declining palynological richness), but closed-canopy forest was maintained throughout. In contrast to the PDFA, the modern distribution of SDTF most likely represents the greatest spatial coverage of these forests in southern South America since glacial times. We find that temperature is a key climatic control upon the distribution of lowland South American SDTF over glacial-interglacial timescales, and seasonality of rainfall exerts a strong control on their floristic composition.

Differentiation between Neotropical rainforest, dry forest and savannah ecosystems by their modern pollen spectra and implications for the fossil pollen record

Accurate differentiation between tropical forest and savannah ecosystems in the fossil pollen record is hampered by the combination of: i) poor taxonomic resolution in pollen identification, and ii) the high species diversity of many lowland tropical families, i.e. with many different growth forms living in numerous environmental settings. These barriers to interpreting the fossil record hinder our understanding of the past distributions of different Neotropical ecosystems and consequently cloud our knowledge of past climatic, biodiversity and carbon storage patterns. Modern pollen studies facilitate an improved understanding of how ecosystems are represented by the pollen their plants produce and therefore aid interpretation of fossil pollen records. To understand how to differentiate ecosystems palynologically, it is essential that a consistent sampling method is used across ecosystems. However, to date, modern pollen studies from tropical South America have employed a variety of methodologies (e.g. pollen traps, moss polsters, soil samples). In this paper, we present the first modern pollen study from the Neotropics to examine the modern pollen rain from moist evergreen tropical forest (METF), semi-deciduous dry tropical forest (SDTF) and wooded savannah (cerradão) using a consistent sampling methodology (pollen traps). Pollen rain was sampled annually in September for the years 1999–2001 from within permanent vegetation study plots in, or near, the Noel Kempff Mercado National Park (NKMNP), Bolivia. Comparison of the modern pollen rain within these plots with detailed floristic inventories allowed estimates of the relative pollen productivity and dispersal for individual taxa to be made (% pollen/% vegetation or ‘p/v’). The applicability of these data to interpreting fossil records from lake sediments was then explored by comparison with pollen assemblages obtained from five lake surface samples.

Global vegetation and terrestrial carbon cycle changes after the last ice age

•In current models, the ecophysiological effects of CO2 create both woody thickening and terrestrial carbon uptake, as observed now, and forest cover and terrestrial carbon storage increases that took place after the last glacial maximum (LGM). Here, we aimed to assess the realism of modelled vegetation and carbon storage changes between LGM and the pre-industrial Holocene (PIH). •We applied Land Processes and eXchanges (LPX), a dynamic global vegetation model (DGVM), with lowered CO2 and LGM climate anomalies from the Palaeoclimate Modelling Intercomparison Project (PMIP II), and compared the model results with palaeodata. •Modelled global gross primary production was reduced by 27–36% and carbon storage by 550–694 Pg C compared with PIH. Comparable reductions have been estimated from stable isotopes. The modelled areal reduction of forests is broadly consistent with pollen records. Despite reduced productivity and biomass, tropical forests accounted for a greater proportion of modelled land carbon storage at LGM (28–32%) than at PIH (25%). •The agreement between palaeodata and model results for LGM is consistent with the hypothesis that the ecophysiological effects of CO2 influence tree–grass competition and vegetation productivity, and suggests that these effects are also at work today.

The influence of vegetation, fire spread and fire behaviour on biomass burning and trace gas emissions: results from a process-based model

A process-based fire regime model (SPITFIRE) has been developed, coupled with ecosystem dynamics in the LPJ Dynamic Global Vegetation Model, and used to explore fire regimes and the current impact of fire on the terrestrial carbon cycle and associated emissions of trace atmospheric constituents. The model estimates an average release of 2.24 Pg C yr−1 as CO2 from biomass burning during the 1980s and 1990s. Comparison with observed active fire counts shows that the model reproduces where fire occurs and can mimic broad geographic patterns in the peak fire season, although the predicted peak is 1–2 months late in some regions. Modelled fire season length is generally overestimated by about one month, but shows a realistic pattern of differences among biomes. Comparisons with remotely sensed burnt-area products indicate that the model reproduces broad geographic patterns of annual fractional burnt area over most regions, including the boreal forest, although interannual variability in the boreal zone is underestimated.

Climate and CO2 controls on global vegetation distribution at the last glacial maximum: analysis based on palaeovegetation data, biome modelling and palaeoclimate simulations.

The global vegetation response to climate and atmospheric CO2 changes between the last glacial maximum and recent times is examined using an equilibrium vegetation model (BIOME4), driven by output from 17 climate simulations from the Palaeoclimate Modelling Intercomparison Project. Features common to all of the simulations include expansion of treeless vegetation in high northern latitudes; southward displacement and fragmentation of boreal and temperate forests; and expansion of drought-tolerant biomes in the tropics. These features are broadly consistent with pollen-based reconstructions of vegetation distribution at the last glacial maximum. Glacial vegetation in high latitudes reflects cold and dry conditions due to the low CO2 concentration and the presence of large continental ice sheets. The extent of drought-tolerant vegetation in tropical and subtropical latitudes reflects a generally drier low-latitude climate. Comparisons of the observations with BIOME4 simulations, with and without consideration of the direct physiological effect of CO2 concentration on C3 photosynthesis, suggest an important additional role of low CO2 concentration in restricting the extent of forests, especially in the tropics. Global forest cover was overestimated by all models when climate change alone was used to drive BIOME4, and estimated more accurately when physiological effects of CO2 concentration were included. This result suggests that both CO2 effects and climate effects were important in determining glacial-interglacial changes in vegetation. More realistic simulations of glacial vegetation and climate will need to take into account the feedback effects of these structural and physiological changes on the climate.

Mid-Holocene and glacial-maximum vegetation geography of the northern continents and Africa.

BIOME 6000 is an international project to map vegetation globally at mid-Holocene (6000 14C yr bp) and last glacial maximum (LGM, 18,000 14C yr bp), with a view to evaluating coupled climate-biosphere model results. Primary palaeoecological data are assigned to biomes using an explicit algorithm based on plant functional types. This paper introduces the second Special Feature on BIOME 6000. Site-based global biome maps are shown with data from North America, Eurasia (except South and Southeast Asia) and Africa at both time periods. A map based on surface samples shows the method’s skill in reconstructing present-day biomes. Cold and dry conditions at LGM favoured extensive tundra and steppe. These biomes intergraded in northern Eurasia. Northern hemisphere forest biomes were displaced southward. Boreal evergreen forests (taiga) and temperate deciduous forests were fragmented, while European and East Asian steppes were greatly extended. Tropical moist forests (i.e. tropical rain forest and tropical seasonal forest) in Africa were reduced. In south-western North America, desert and steppe were replaced by open conifer woodland, opposite to the general arid trend but consistent with modelled southward displacement of the jet stream. The Arctic forest limit was shifted slighly north at 6000 14C yr bp in some sectors, but not in all. Northern temperate forest zones were generally shifted greater distances north. Warmer winters as well as summers in several regions are required to explain these shifts. Temperate deciduous forests in Europe were greatly extended, into the Mediterranean region as well as to the north. Steppe encroached on forest biomes in interior North America, but not in central Asia. Enhanced monsoons extended forest biomes in China inland and Sahelian vegetation into the Sahara while the African tropical rain forest was also reduced, consistent with a modelled northward shift of the ITCZ and a more seasonal climate in the equatorial zone. Palaeobiome maps show the outcome of separate, independent migrations of plant taxa in response to climate change. The average composition of biomes at LGM was often markedly different from today. Refugia for the temperate deciduous and tropical rain forest biomes may have existed offshore at LGM, but their characteristic taxa also persisted as components of other biomes. Examples include temperate deciduous trees that survived in cool mixed forest in eastern Europe, and tropical evergreen trees that survived in tropical seasonal forest in Africa. The sequence of biome shifts during a glacial-interglacial cycle may help account for some disjunct distributions of plant taxa. For example, the now-arid Saharan mountains may have linked Mediterranean and African tropical montane floras during enhanced monsoon regimes. Major changes in physical land-surface conditions, shown by the palaeobiome data, have implications for the global climate. The data can be used directly to evaluate the output of coupled atmosphere-biosphere models. The data could also be objectively generalized to yield realistic gridded land-surface maps, for use in sensitivity experiments with atmospheric models. Recent analyses of vegetation-climate feedbacks have focused on the hypothesized positive feedback effects of climate-induced vegetation changes in the Sahara/Sahel region and the Arctic during the mid-Holocene. However, a far wider spectrum of interactions potentially exists and could be investigated, using these data, both for 6000 14C yr bp and for the LGM.

Intercomparison of simulated global vegetation distribution in response to 6 kyr BP orbital forcing

The response of ten atmospheric general circulation models to orbital forcing at 6 kyr BP has been investigated using the BIOME model, which predicts equilibrium vegetation distribution, as a diagnostic. Several common features emerge: (a) reduced tropical rain forest as a consequence of increased aridity in the equatorial zone, (b) expansion of moisture-demanding vegetation in the Old World subtropics as a consequence of the expansion of the Afro–Asian monsoon, (c) an increase in warm grass/shrub in the Northern Hemisphere continental interiors in response to warming and enhanced aridity, and (d) a northward shift in the tundra–forest boundary in response to a warmer growing season at high northern latitudes. These broadscale features are consistent from model to model, but there are differences in their expression at a regional scale. Vegetation changes associated with monsoon enhancement and high-latitude summer warming are consistent with palaeoenvironmental observations, but the simulated shifts in vegetation belts are too small in both cases. Vegetation changes due to warmer and more arid conditions in the midcontinents of the Northern Hemisphere are consistent with palaeoenvironmental data from North America, but data from Eurasia suggests conditions were wetter at 6 kyr BP than today. The models show quantitatively similar vegetation changes in the intertropical zone, and in the northern and southern extratropics. The small differences among models in the magnitude of the global vegetation response are not related to differences in global or zonal climate averages, but reflect differences in simulated regional features. Regional-scale analyses will therefore be necessary to identify the underlying causes of such differences among models.

Postglacial vegetation, fire and climate dynamics at Central Chilean Patagonia (Lake Shaman, 44ºS)

A 19 cal ka BP pollen and charcoal record from Lake Shaman (44°S; 71°W, Chile) was analyzed to establish vegetation, fire and climate dynamics of the forest-steppe ecotone in Central Chilean Patagonia. Lake Shaman record indicates that the upper Río Cisnes valley was free of ice at around 19 cal ka BP. From this date and until 14.8 cal ka BP, a grass steppe with high proportions of shrubs associated to colder and drier conditions than present developed in this area. A continuous increase of Nothofagus accompanied by a decline in the steppe shrubs and sudden dominance of paludal over aquatic plants from 11 cal ka BP was associated to effective moisture increase but still under modern values. The replacement of the cold-dry grass-shrub steppe by a similar-than-present forest-steppe ecotone suggests an increase in temperature indicating the onset of the Holocene. At the same time, moderate fire activity suggested by the charcoal record could be related to major fuel availability as consequence of Nothofagus forest expansion. Between 8 and 3 cal ka BP, the record indicates the easternmost position of the forest-steppe ecotone suggesting the highest effective moisture with the establishment of seasonality between 5 and 3 cal ka BP. From 3 cal ka BP, the record indicates a retraction of the forest-steppe ecotone accompanied by a high pollen record variability and an increased fire activity. These late changes suggest decreased effective moisture associated with a high climatic variability. At regional and extra-regional scale, climatic changes at Lake Shaman's record are mostly associated to changes (latitudinal shifts and/or strengthening/weakening) of past Southern Westerlies that were previously recorded along Patagonia from the Lateglacial to the mid-Holocene. During the Late Holocene, a regional pattern characterized by high record variability emerges throughout Central Chilean Patagonia. This variability would be related to (1) low magnitude Southern Westerlies changes probably associated to ENSO and/or SAM or (2) the complex relationships between vegetation, fire and human occupations during the last 3 cal ka.

The response of vegetation on the Andean flank in western Amazonia to Pleistocene climate change

A reconstruction of past environmental change from Ecuador reveals the response of lower montane forest on the Andean flank in western Amazonia to glacial-interglacial global climate change. Radiometric dating of volcanic ash indicates that deposition occurred ~324,000 to 193,000 years ago during parts of Marine Isotope Stages 9, 7, and 6. Fossil pollen and wood preserved within organic sediments suggest that the composition of the forest altered radically in response to glacial-interglacial climate change. The presence of Podocarpus macrofossils ~1000 meters below the lower limit of their modern distribution indicates a relative cooling of at least 5°C during glacials and persistence of wet conditions. Interglacial deposits contain thermophilic palms suggesting warm and wet climates. Hence, global temperature change can radically alter vegetation communities and biodiversity in this region.