Abstract
The present day distribution of Miombo savanna-woodland in Mozambique has been attributed to an expansion due to the clearing of original coastal forests through agriculture and use of fire. Here, we test this hypothesis using palaeoecological data from Lake Nhauhache, situated in the Vilankulo region. Our analysis shows that Brachystegia, one of the main constituents of the Miombo, has varied over time, and its variability seems to be driven by hydrological changes related to climatic variability rather than by land-use changes. The analyses show that Brachystegia was most common during
Introduction
The coastal lowlands of southern Mozambique and the Vilankulo area are commonly reported as a forest or forest transition region (White, 1983) or as a Miombo savanna-woodland with a coastal ticket or forest along the sea stretch (Wild and Grandvaux Barbosa, 1967). The forest or forest transition region, the Indian Ocean coastal belt, is described as a mosaic of semi-deciduous forest, thickets, woodland, savannas and edaphic grasslands (Moll and White, 1978; Werger, 1978). Despite representation in vegetation maps, islands of semi-deciduous forests, are, however, only present in the Inhambane region and by the Save River (Wild and Fernandes, 1967). Miombo is a colloquial term of savanna-woodland dominated by the genera Brachystegia, Julbernardia and/or Isoberlinia (Frost, 1996). In the phytogeographies of Werger (1978), the Miombo savanna-woodland was suggested to be invasive on the coast with its original centre of endemism in the south-central Africa (hence the name Zambesian phytochorion). The expansion of the Miombo savanna-woodland has been attributed to the use of fire and farming in the coastal regions. In previous discussions on landscape dynamics (e.g. Hall, 1981; Morais, 1988), it has been assumed that the coastal regions were forested prior to the arrival of farmers in southern Africa c.
Earlier palaeoenvironmental investigations in the Vilankulo region have shown that the coastal area in
Here, we present a new vegetation record from the Vilankulo region to explore the Miombo/forest dynamics. The analysis is based on stratigraphic occurrences of pollen, spores and microfossil charcoal from Lake Nhauhache and combined with the diatom assemblage presented previously (Holmgren et al., 2012). The results are discussed in relation to changes in land use, and in regional vegetation and hydrological changes, with a specific focus on the long-term ecology of the Miombo.
There are several processes that can be responsible for the hypothesised transition from forests to a Miombo savanna-woodland. One process is land-use change and swidden agriculture, whereby farmers converted coastal forests into agricultural lands and fallow periods were too short to allow forests to return. Another process that would have promoted Miombo constituents on the cost of forest species is a possible shift in fire regimes (e.g. seasonality, intensity and frequency). As the Miombo constituents are well adjusted to the influence of fire, they sprout from the root, and Miombo woodlands in the African interior have also been reported to burn approximately in four cycles (Chidumayo, 1997). These two processes are likely to be interlinked as forest fires usually tend to spread mainly on the forests floor, for example, not affecting trees once they are matured. A third possibility that will be tested here is that the presence of Miombo constituents can be explained as a long-term constituent of the coastal landscape and as unrelated or weakly related to land use. Indirect evidence of farming communities in southern Mozambique is recorded from
Vilankulo region and research background
The Vilankulo region is marked by the presence of several lake systems, most of which are closed lakes with a characteristic circular or semi-circular shape. The Pleistocene dune system of the region forms an aquifer independent from the interior that is fed mainly by rainfall (Coetsee and Hartley, 2001). Vilankulo receives most of its rainfall (832 mm/yr) between November and March, when temperatures are the highest. The region experiences a c. 20-year rainfall cyclicity (Tyson and Preston-Whyte, 2000), in addition to high inter-annual rainfall variability. The dune system of the coastal plain is characterised by weakly developed soils with a low organic content and moderate weathering (Fränzle, 1984; Wood, 2001). These soils support a savanna and woodland savanna where J. globiflora and S. spinosa are well represented, interspersed with grasses. B. spiciformis is rarely found in the Vilankulo region (Massinga in Ekblom, 2004; Telford and De Castro, 2001). In some parts, the coastal dunes support a dense shrub characterised by Commiphora zanzibarica, Phyllanthus reticulatus and Turrea nilotica. These occur with Grewia monticola, Deinbollia oblongifolia, Clerodendrum glabrum and Acalypha glabrata and also climbers such as Cocculus hirta (Massinga in Ekblom, 2004). Coastal forests, constituted by species such as Chlorophora excelsa, Ficus spp., Morus mesozygia, Celtis africana, Afzelia quanzensis, Dialium schlechteri and Brachyleana discolor, are absent in the Vilankulo region today (Ekblom, 2004; Wild and Fernandes, 1968). However, patches of closed canopies sometimes occur and are characterised by species, such as Millettia stuhlmanni, Alchornea laxiflora, Euclea natalensis and A. quanzensis. Along streams, species such as Ficus trichopoda and Myrica cf. pilulifera are recorded together with J. globiflora and, occasionally, B. spiciformis (Telford and De Castro, 2001).
The earliest date of farming communities in Mozambique comes from south of Maputo, where early farming community–style pottery (Matola) was found in layers dated to
Palaeoecological investigations have been carried out in Lake Nhaucati and Lake Xiroche in the Chibuene area (Ekblom, 2008). These records span from
Lake Nhauhache
Lake Nhauhache is situated at 21°58′50″S, 35°17′39″E (at 6 m a.s.l.), about 3 km northwest of Vilankulo town (Figure 1). Nhauhache is uppermost in a system of smaller coastal lakes linked to the Inhamalene River, which feeds into the Indian Ocean. The lake varies in size depending on regional rainfall. In July 2008, at the time of coring, it measured 150 m in diameter and was 2.2 m deep. The lake is characteristic of other lakes in this region, circular with steep surrounding slopes. In the southeast, the flatter topography allows for drainage of the lake during high lake levels.

Location of the Vilankulo region; (inset) Vilankulo town, Lake Nhauhache and other investigated lakes in the vicinity.
The surroundings of the lake are today densely populated and cultivated, and the biological production within the lake is high. Hydrophilous plants, such as Nymphaea lotus, Hyparrhenia sp. and Phragmites sp., are found in the limnic and telmatic zones. Local residents reported that the lake was dry from
Methodology
Lake Nhauhache was sampled in July 2008 using a dinghy and a Russian coring device. The bathymetry was examined using a handheld echo-sounder, and a sediment core was subsequently collected in the deepest part of the lake, at c. 2.2 m water depth. Lacustrine sediment with a thickness of 2.3 m was cored before it reached a layer of pure dune sand. All sediment depths in this paper are referred to in terms of centimetre from the bottom of the lake. The lithology of the lake sediment is homogeneous and composed of reddish to grey sandy gyttja. In the lower portion of the core, nodules of iron are observed, and a 1-cm-thick layer of rootlets appears at 355 cm depth. There are no clear signs of erosion horizons or palaeosoil surfaces. As the lake is small, we expect it to reflect local conditions of vegetation change (Jacobsen and Bradshaw, 1981; Prentice, 1985).
Pollen
Pollen and microscopic charcoal were extracted through digestion of cellulose with NaOH and concentration of pollen through acetolysis (Faegri and Iversen, 1992; Moore et al., 1991). A known amount of exotic Lycopodium marker spores to the sediment to calculate total concentration of pollen grains per square centimetre of sediment. Every 8 cm (or more) was sampled in the upper part of the core. As the original focus of this investigation was on the period
The pollen of indigenous African cereal grasses (Pennisetum americanum, Sorghum bicolor and Eleusine coracana) overlap in size and are morphologically similar to wild grass pollen (Tomlinson, 1973). It has been noted in pollen analyses from the Mozambique interior (Ekblom and Gillson, 2010) that larger grass pollen (e.g. <40 µm) tend to occur in greater numbers together with maize pollen; thus, it is likely that the grass pollen >40 µm represent indigenous cereal grasses. We have therefore used a tentative size separation of grasses of <40 µm here. Glycerol tends to swell pollen grains to some extent (Cushing, 1961). There may therefore be an overrepresentation of grains in the <40 µm size class due to the effect of glycerol. We have only included pollen grains that are decisively >40 µm in size, and we did not note any significant changes in sizes in other taxa.
Charcoal
The area of charcoal, that is, black opaque angular particles, was estimated in relation to sediment volume (cm2/cm3) and added Lycopodium spores, using the point count method, with counting minimum of in total of 50 items (Clark, 1982, 1988). This gives a rough estimate of the amount of charcoal in the sample, although a larger count, of 200 Lycopodium/charcoal, would be ideal for higher statistical reliability (Finsinger and Tinner, 2005). We also separated two main classes of charcoal: 10–100 µm interpreted as linked to local–regional fires (here referred to as ‘regional charcoal’) and above 100 µm (here referred to as ‘local charcoal’) interpreted as linked to local fires (see review in Pitkänen et al., 1999).
Spores
All spores were counted in relation to the added Lycopodium spores and presented as concentrations. Most types were identified or grouped together into main type groups (Graf and Chmura, 2006; Jarzen and Elsik, 1986; Prager et al., 2006; Van Geel and Aptroot, 2006; Van Geel et al., 1983, 1986, 2003). Here, we will focus on the coprophilous fungi, mainly Sordariaceae, but also Coniochaeta lignaria (Marinova and Atanassova, 2006; Van Geel and Aptroot, 2006; Van Geel et al., 2003). Algal spores from Spirogyra spores and Glomus (only found in very low numbers and not shown here) were also distinguished. Botryococcus, most probably Botryococcus braunii, occurred in large numbers throughout the core but was not quantified as they occur as mass material.
Diatoms
The diatom assemblage has been presented in detail elsewhere (Holmgren et al., 2012). Only the ecological groups are presented here, that is, planktonic, benthic, halophilous, aerophilous and ungrouped taxa. Every 3–4 cm intervals were sampled for diatoms, and the samples were oxidised with 17% H2O2. Clay particles were removed by decanting from 100 mL beakers in 2-h intervals (Battarbee, 1986). The remaining samples were mounted on glass slides using Naphrax and studied under a light microscope using oil immersion. Approximately 500 diatom frustules were counted per level using standard references (Cleve-Euler, 1953, 1955; Gasse, 1986; Krammer and Lange-Bertalot, 1986, 1988, 1991a, 1991b; Tynni, 1979).
Dating
Four samples were accelerator mass spectroscopy (AMS) 14C dated at the Radiocarbon Laboratory in Poznan, Poland, to complement the 10 samples presented by Holmgren et al. (2012; Table 1). Due to the lack of macrofossils in the sediments, bulk sediments (1 cm thickness each) were used for dating. OxCal 4.0 (Bronk Ramsey, 1995) was used for calibration to calendar scale with the calibration curve for the southern hemisphere (McCormac et al., 2004; Figure 2). The new 14C dates strongly support the previously proposed age model (Holmgren et al., 2012). The record spans the period from c. 2300
Calibration dates for Lake Nhauhache. Probability at 1σ and 2σ. Calibration curve for the southern hemisphere used (McCormac et al., 2004). New dates are marked in bold.

Age–depth diagram with the range of dates published in Holmgren et al. (2012, grey) and additional dates (black). The probability distribution (2σ) of each date is shown in the respective graph. The pollen and diatom zones are shown to the right.
Results
The pollen assemblage has been divided into three zones through cluster analysis. The analysis is based on unweighted paired group average and Euclidean similarity index using the software PAST (Hammer et al., 2001). This method is unconstrained by depth and more suitable for cores with possible temporal gaps (Figure 3). The most commonly represented species is presented here (Figures 4 and 5). The tree/shrub group has been divided into two main ecological groupings to assist interpretation, a riverine forest/forest group and a savanna/generalist group (Table 2).

Cluster analyses of analysed levels, the analysis that is unconstrained by depth is based on unweighted paired group average and Euclidean similarity index using the software PAST (Hammer et al., 2001).

Pollen diagram, with terrestrial pollen percentages grouped into arboreal types and savanna and generalist types. The lithology consists of sandy gyttja. Herbs and aquatics are also shown. The pollen types with a low representation have been exaggerated (hatched silhouette) with a factor of 5.

Summary pollen diagrams with the main ecological groups shown in percentages. The lithology consists of sandy gyttja. Distribution of spores (number/cm2), regional and local charcoal (number/cm2) and total pollen and spore concentrations (left). The pollen diagram has been divided into three zones based on the pollen assemblage. A summary diatom diagram is shown to the right with the four zones based on the diatom assemblage indicated for comparison. The sparse representation of charcoal in the lower part of the diagram is due to the low sampling resolution here.
Represented pollen taxa and ecological groupings in the Lake Nhauhache sediment core.
Nhauhache 1: 210–135 cm depth (ad 200–700)
Total pollen concentrations are generally low in this zone. Trees and shrubs are well represented, particularly Stereospermum, Combretaceae, Myrsine Africana, Alchornea (well represented throughout the core), Moraceae and Trema, are here combined into one curve as their representation is low, and they have similar ecological requirements. The identification of Stereospermum remains provisional (see Table 2). Brachystegia shows its highest values in this zone compared with other parts of the core. There are high values of Nymphaea throughout the zone together with Laurembergia and Hydrocharitaceae. Poaceae and Cyperaceae show moderate levels. The herbs (Chenopodiaceae, Asteraceae and Acanthaceae) are relatively well represented. Charcoal particles and Poaceae >40 µm are few in this zone. Spores are represented in low numbers and represented mainly by fern spores.
Nhauhache 2: 135–42.5 cm depth (ad 700–1500)
This zone could possibly be divided into two subzones 135–90 cm and 90–42.5 cm (
Nhauhache 3: 42.5–10 cm depth (ad 1500–1800)
The uppermost zone is of short duration and is distinguished by an increase in trees/shrubs and in Cyperaceae and the aquatics Nymphaea and most markedly of Typha. Among the trees/shrubs, Alchornea remains well represented together with Myrtaceae and Phyllanthus. Stereospermum and Brachystegia occur only in low numbers, and the latter is absent in the upper part of the zone. Dodonaea is common (as also in the upper part of the previous zone). The herbs show a significant increase, particularly the Lamiaceae group. Fern and fern/moss spores occur in moderate numbers, and there is an overall increase in spores, particularly spores associated with coprophilous fungi C. lignaria and Sordariaceae that are very common in the uppermost part of the core. Local and regionally transported charcoal remains well represented together with grasses >40 µm. Pollen concentrations in this zone are higher than in previous zone, and the amount of degraded grains remains high.
Discussion
Summarising the changes in vegetation communities
The results of the Lake Nhauhache pollen analysis show that Brachystegia was most common in the beginning of the sequence. At
In the period
From
The period from
Testing the hypothesis of Miombo expansion
Brachystegia is the only species of the Miombo that frequently occurs in the pollen diagram, Julbernardia type has been identified but very rarely, and this is a pattern also in other pollen analyses as discussed in section ‘Introduction’. Therefore, Brachystegia is here treated as a representative of the Miombo. Brachystegia is underrepresented in the pollen diagram, and even a low presence is therefore interpreted to indicate considerable coverage. It is unlikely that the low representation of forest taxa throughout the diagram can be explained as the result of land use and/or shift in fire regime as was posited in section ‘Introduction’. It cannot be excluded that farmers were active in the area from
Our results indicate that Brachystegia has varied over time, with the highest representation at
Grass pollen <40 µm that may represent cereal grasses is represented throughout the core, and in the base of the diagram, this is most likely wild grasses. An increase in grass pollen <40 µm can be seen in the transition between zones 2a and 2b at
Brachystegia disappears at the top of the record, from
Conclusion
The Lake Nhauhache sequence shows that Brachystegia was most common at
Footnotes
Funding
The work was, in part, carried out within the framework of Karin Holmgrens VR-funded project Holocene climate variability in southern Africa.
This research received no specific grant from any funding agency in the public, commercial or not-for-profit sectors.
