1 1 Does long-term fire exclusion in an Australian tropical savanna 2 result in a biome shift? A test using the re-introduction of fire. 3 4 5 Kenneth ScottA, Samantha A. SetterfieldA*, Michael M. DouglasA,B, Catherine L. ParrC, Jon SchatzD, and Alan N. AndersenD. 6 7 A 8 Darwin, NT, 0909, Australia. 9 B Research Institute for Environment and Livelihoods, Charles Darwin University, Tropical Rivers and Coastal Knowledge Research Hub, Charles Darwin University, 10 Darwin, NT, 0909, Australia. 11 C 12 Environment, Oxford OX1 3QY, England 13 D 14 0822, Australia. University of Oxford, Environmental Change Institute, School of Geography & CSIRO Ecosystem Sciences, Tropical Ecosystems Research Centre, Winnellie, NT, 15 16 * Corresponding author: 17 Dr Samantha Setterfield 18 Research Institute for the Environmental and Livelihoods 19 T. +61 8 8946 6756 | F. +61 8 8946 6847 | M. 0427 456 031 20 [email protected] 21 22 23 Running Title: Re-introducing fire to an Australian savanna 2 1 Abstract 2 The structure of tropical savanna ecosystems is influenced by fire frequency and 3 intensity. There is particular interest in the extent to which long-term fire exclusion can 4 result in a shift from savanna to forest vegetation that is not easily reversed by the re- 5 introduction of fire. This study examined changes in the structure and composition of a 6 long-unburnt site within the northern Australian savannas following a long period of 7 active fire exclusion (>20 years), and the effect of the re-introducing fire through 8 experimental fire regimes, including fires in the early and late dry season at a range of 9 frequencies. After the long period of fire exclusion, the vegetation community was 10 characterised by a well developed mid-story and canopy layer, low grass cover, 11 substantially higher densities of woody sprouts and saplings than frequently burnt 12 savanna. The community composition included a high proportion of rainforest-affiliated 13 species. Three years of experimental fires had no detectable effect on the overall 14 composition of grass-layer and woody plants but had an effect on woody vegetation 15 structure. Continued fire exclusion further increased the density of woody stems, 16 particularly in the mid-storey (2.0-4.99 m), whereas moderate-intensity fires (> 800 kW 17 m-1) significantly reduced the density of mid-storey stems. The re-introduction of higher 18 moderate intensity fire events resulted in the vegetation in some compartments reverting 19 to the open savanna structure typical of frequently burnt sites. Such rapid reversibility 20 suggests that in general, the woody thickening resulting from long-term fire exclusion 21 did not represent a biome shift to a non-savanna state. However, there was a small 22 proportion of the site that could not sustain the fires applied to them because grass cover 23 was very low and patchy and therefore appeared to have crossed an ecological threshold 24 towards closed forest. 25 Keywords. savanna, fire, biome shift, woody thickening, tropical woodland 3 1 2 INTRODUCTION State-and-transition models are particularly useful for describing the range of 3 complex ecosystem processes and multiple vegetation states observed in tree-grass 4 ecosystems (Westoby et al. 1989; Angassa & Oba 2008). According to state-and- 5 transition theory, vegetation can shift from one stable state to another, and this shift can 6 be irreversible without dramatic intervention (Briske et al. 2003; Stringham et al. 2003; 7 Bestelmeyer et al. 2003; Briske et al. 2006). Substantial changes in composition and 8 structure may occur within a vegetation state, but negative feedbacks reinforce 9 ecosystem stability, and changes are easily reversed by removing the stressor/s that 10 11 drove the initial change (Bestelmeyer et al. 2003; Briske et al. 2006). In extreme cases a transition from one state to another can involve such a 12 transformation in ecosystem structure and function that it represents a shift from one 13 biome to anther. This is the case throughout the seasonal tropics, where savanna and 14 forest represent contrasting biomes in regions of high rainfall (Boulière 1983; Bond 15 2008; Bond & Parr 2010). Savanna vegetation is dominated by fire-resistant, shade- 16 intolerant, sclerophyll tree species that form an open canopy over a continuous grass 17 layer, whereas forest is dominated by fire-sensitive, shade-tolerant mesophyll tree 18 species that form a closed canopy with little or no grass cover (Ratnam et al. 2011; 19 Thompson et al. 2001; Wilson et al. 1990). Long-term fire exclusion sees a transition 20 from savanna to forest, as juvenile woody plants to recruit into the midstorey layer, and 21 the increased shade reduces grass biomass (Brookman-Amissah et al. 1980; San Jose & 22 Farinas 1991; Silva et al. 2001; Russell-Smith et al. 2003; Woinarski et al. 2004; 23 Williams RJ et al. 2003). The reduction in grass biomass results in reduced fire 4 1 frequency, which reinforces the transition from savanna to forest (Ratnam et al. 2011; 2 Scholes & Archer 1997; Russell-Smith et al. 2003; Sharp & Whittaker 2003). 3 Alternatively, frequent burning has been shown to reverse woody thickening, 4 particularly in combination with drought and browsing by large herbivores 5 (Freudenberger et al. 1997; Roques et al. 2001; Augustine & McNaughton 2004; Holdo 6 2007, Midgley et al. 2010). Most savanna woody species are capable of re-sprouting 7 after fire, although high mortality occurs in response to higher-intensity fires, 8 particularly of the smallest and largest individuals (Midgley et al. 2010; Setterfield 9 1997; Williams RJ et al. 2003; Werner 2005). Conversely, dry season burning 10 generally has little effect on grass-layer plants, with compositional change attributed to 11 inter-annual rainfall variability (O’Connor 1985; O’Connor & Everson 1998; Williams 12 RJ et al. 2003; Williams PR et al. 2003b). Site and stand characteristics, and the 13 interaction between these, may alter these responses because they determine fire 14 behaviour and the response of the vegetation community to fire behaviour (Williams et 15 al. 2003; Bowman & Fensham 1991). Fensham et al. (2003) demonstrated a major (78- 16 89%) reduction in tree basal area following burning of a fire-protected savanna, where 17 rainforest-associated species and obligate seeders were particularly susceptible to 18 repeated burning. 19 Open savanna communities characterised by a continuous understorey of C4 20 grasses and a scattered overstory of Eucalyptus species dominate northern Australia 21 (Setterfield and Hutley 2008). These communities are frequently burnt every 1-3 yrs 22 (Andersen et al. 2005). Within the savanna matrix are small (typically <5ha) pockets of 23 closed monsoon rainforest (Russell-Smith 1993; Russell-Smith and Setterfield 2006). A 24 range of studies in the savannas of northern Australia have documented transition 5 1 towards closed forest in the long-term absence of fire, characterised by increases in 2 stem density, development of the midstory or subcanopy, decrease in grass cover and 3 increase in rainforest-associated species (Russell-Smith et al. 2003; Williams PR et al. 4 2003a; Woinarski et al. 2004). However, it is unclear at what point along this transition 5 the shift from savanna to forest can be considered to have occurred. One approach to 6 addressing this is to focus on the ease with which vegetation change is reversible. If the 7 re-introduction of fire results in a rapid reversion to open savanna then it can be argued 8 that a biome shift has not yet occurred. 9 Here we use this approach to assess if >20 years of active fire exclusion at a 10 savanna site in northern Australia has resulted in a biome shift. First, we describe the 11 state of the vegetation in the long-term absence of fire, with the expectation that it will 12 have characteristics of forest structure, composition and function, such as a substantial 13 proportion of forest species, a relatively closed canopy, and patchy grass cover (and 14 therefore a low capacity to carry fire). We then assess the vegetation response to the re- 15 introduction of fire focussing on the extent to which these forest characteristics are 16 reversed. 17 METHODS 18 Study area 19 The study was undertaken within the Territory Wildlife Park (TWP) at Berry 20 Springs, 30 km southeast of Darwin, Northern Territory, Australia (12º 41′ 42.25″ S, 21 130º 58′ 50.36″ E). The region experiences a monsoonal tropical climate, with a hot 22 summer wet season between October and March, and a mild winter dry season between 23 April and September. Mean annual rainfall at the Middle Point Automatic Weather 6 1 Station (36 km to the north-east) is 1399 mm, of which 90% falls in the wet season 2 (Bureau of Meteorology, unpublished data, 2007). Rainfall recorded at the Park during 3 the study (2004-2007) was 17% higher than the long-term average at Middle Point, 4 principally as a result of a cyclone-related rainfall event in April 2006 (Bureau of 5 Meteorology and L Smith, unpublished data, 2007). 6 Fire has been actively excluded from the Park since at least 1981. Only two 7 patchy fires had occurred in the study site between 1981 and the application of 8 experimental fire regimes in 2004. In 1980, vegetation structure at the site was typical 9 of the broader region; a sparse midstorey (approximately 2-5 m high) over a grass-layer 10 of tall perennial (Heteropogon triticeus) and annual (Sorghum intrans) grasses 11 (Sivertsen et al. 1980). Eucalyptus tetrodonta F.Muell. and Eucalyptus miniata 12 A.Cunn. ex Shauer were dominant overstorey trees (Sivertsen et al. 1980; Fig. 1). At 13 the beginning of this study (2004), a moderately- to well-developed midstorey was 14 present throughout the site (Fig. 1), particularly in association with the Quaternary 15 colluvial sediments of the permanent lagoon in the south (Scott et al. 2009). Common 16 midstorey shrubs included Calytrix exstipulata DC., Planchonia careya (F.Muell.) 17 Kunth, Petalostigma pubescens Domin and Exocarpus latifolius R.Br. Common 18 herbaceous plants in the grass-layer were low-biomass species, such as the annual grass 19 Pseudopogonatherum contortum (Brongn.) A.Camus, sedge Fimbristylis sp. Charles 20 Darwin (J.L.Egan 5300) and perennial grass Eriachne triseta Nees ex Steud. (Scott et 21 al. 2009). 22 Experimental Design 23 The fire experiment at the TWP consisted of a line of 18 1-ha compartments 24 located in the north-western corner of TWP, arranged in three contiguous blocks (North, 7 1 Centre and South) to accommodate the environmental variability at the site. A detailed 2 description of the variation is provided in Scott et al. (2009). In summary, soil moisture 3 content of the surface soil was significantly higher in the North and Central blocks 4 where the soils are gravelly clay loams on a small ridge, and exhibit an horizon of 5 massive laterite of Tertiary origin, 10–50 cm below the surface (Scott 2008). In the 6 South block, soils were sandy loams associated with Quaternary colluvial sediments of 7 a permanent lagoon, significantly less gravelly (25% gravel >2mm in the upper 5 cm 8 soil surface), positioned at the foot of the ridge, and did not contain massive laterite in 9 the upper 1m of the soil profile (Scott 2008). Vegetation structure varied along the 10 edaphic gradient, with significantly higher density of woody stems in the midstorey in 11 the South block (~2000 stems ha-1) compared with the North and Central Blocks (1000- 12 1100 stems ha-1). 13 Fire regimes were allocated as treatments to three compartments each, one per 14 block, according to a randomised complete block design. This allowed variation in 15 edaphic characteristics to be assessed in the experiment. Fire regimes were: unburnt 16 (UN), annual early fires (Annual Early), a single early season fire in 2004 (Early 04), a 17 single early season fire in 2005 (Early 05), and biennial late (Biennial Late). All fire 18 regimes were first applied in 2004, with the exception of the Early 05 regime which was 19 first burnt in 2005 (Table 1). Experimental fires were lit with drip torches along the 20 leeward sides of each compartment to be burnt (back fire). After the fire had 21 established within the compartment, fire lines were lit on the windward side unless the 22 back fire was intense and moving rapidly to the windward side. This paper reports on 23 changes in vegetation in the first three years following the reintroduction of fire. The 24 early season fires were typically lower intensity, and therefore lower flame height, than 8 1 late season fires which reflects the fire behaviour in this region (Table 1; Andersen et al. 2 2005; Setterfield et al. 2010). The exception was the Early 05 treatment which had a 3 moderate intensity fire, and contrasted with the early fires in other treatments. There 4 was variability within fire regimes due to the edaphic and vegetation differences 5 between compartments and blocks. For example, very low fire intensities occurred in 6 some compartments, particularly in the South block, where dense woody vegetation and 7 high canopy cover resulted in sparse grass cover (Fig 1b, Scott et al. 2009). 8 Plant composition and abundance 9 To determine if the reintroduction of fire caused a substantial change in structure 10 and composition of woody plants, surveys were completed in May 2004, before fire 11 regimes were first implemented, and again in March 2007. Sampling was undertaken 12 within two parallel belt transects (80 × 6 m, approximately 30 m apart) within each 13 compartment, where the abundance of woody plant species (including standing dead 14 stems) was recorded. For stems < 2 m, abundance was categorised according to two 15 height strata: sprouts (0-0.49 m) and small saplings (0.5-1.99 m) (collectively known as 16 the grass-layer stratum). For stems > 2 m tall (and > 3 cm dbh), the height and dbh of 17 each individual was recorded, which enabled categorisation into the following height 18 strata: lower midstorey (2.0-4.99 m), upper midstorey (5.0-9.99 m) and overstorey (10.0 19 m+). Stems > 2m, but < 3 cm dbh, were regarded as large saplings. The precise height 20 of large saplings was not measured, although they were typically < 5 m (i.e. in the lower 21 midstorey), and contributed only a very small proportion (2.4%) of the total woody stem 22 density (and therefore basal area). Basal area (m2 ha-1) was calculated by multiplying 23 stem area (obtained from dbh measurements) by stem density, considering each woody 9 1 stem > 2 m height and > 3 cm dbh. Eight uncommon vine species were present but 2 were not included in the site assessment. 3 Most savanna woody plants can regenerate vegetatively following fire. To 4 document the extent of structural change associated with stem mortality and subsequent 5 basal sprouting (where no mortality takes place), stems in the 2007 survey were noted if 6 this response had occurred (individuals are termed ‘basal sprouters’). 7 To assess the effect of fire regimes on the grass-layer species, the abundances of 8 grass, sedge and forb species were sampled annually in the late wet season from 2004 to 9 2007. The density of each species was counted within 18 2-m2 permanent quadrats per 10 compartment; nine quadrats were spaced, every 10 m, along the centre line of each 11 swathe established to sample the woody vegetation. To aid the interpretation of floristic 12 responses to fire, tree litter and canopy cover were also assessed annually at these 13 quadrats. The percent cover of grass and leaf litter was estimated visually (to the nearest 14 5%), and percent canopy cover recorded using the mean of four cardinal points on a 15 forestry densitometer immediately outside each quadrat (Lemmon 1957). 16 Sampling occurred in the late wet season (March-May) within a narrow window 17 of opportunity when the ability to detect and identify grass-layer plants was maximised, 18 given the presence of mature, flowering individuals. A slightly later sampling period in 19 2004 made for unreliable identification of annual species (grasses, forbs and sedges), 20 because they had already set seed and begun to cure. With the exception of perennial 21 grasses and selected annual grass species, most analyses therefore consider the surveys 22 of 2005 to 2007 inclusive. 10 1 For some other plants, identification to species level was not possible given the 2 absence of flowering material during the sampling period. Such species were included 3 as aggregates of two or more species with a similar vegetative appearance: Eriachne 4 agrostidea F.Muell included E. ciliata R.Br.; Panicum mindanaense Merr. included 5 Yakirra nulla Lazarides & R.D. Webster and Urochloa polyphylla (R.Br.) R.D. 6 Webster; Eriachne triseta Nees ex Steud. included Aristida holathera Domain and E. 7 burkittii Jansen; and Eriachne avenacea R.Br. included E. major (Ewart & O.B.Davies 8 Lazarides). 9 10 Statistical analysis Analysis of variance was used to assess changes in the abundance of woody 11 vegetation within each height stratum, species richness, and basal area. Data from 12 transects were averaged and the variables were analysed with multivariate ANOVA 13 with the following factors: Year (fixed), Regime (fixed), Block (random), with Block 14 nested in Regime. Variables were square-root or log-transformed if necessary, in this 15 and all other analyses, to satisfy the normality and variance homogeneity assumptions 16 of ANOVA. 17 Analyses of woody stem abundance were also conducted separately (in each 18 stratum) for the rainforest affiliated species (defined by Liddle et al. 1994 as species 19 that are “either exclusive to, or typically components of, perennially or seasonally 20 closed canopy vine-forest or vine-thicket vegetation”), given expected differences in 21 their responses to fire compared with savanna species (Fensham et al. 2003). Analyses 22 of woody stem abundance was also conducted for species based on the height of those 23 species when they are mature individuals (grass-layer species (<2m), 6 spp.; shrub 24 species (2 -10m) , 58 spp.; tree species (>10m), 12 spp.). 11 1 Two-way ANOVA was used to determine whether fire regime influenced the 2 incidence of basal sprouting. The number of woody stems exhibiting the basal 3 sprouting response (as a proportion of the total live stem density) was compared 4 between fire regime treatments in a design consisting of Regime (fixed), and Block 5 (random), with block nested in fire regime. 6 To make generalisations on grass-layer plant responses, species were 7 categorised according to life form: annual or perennial grasses (Poaceae), non- 8 leguminous forbs, sedges (Cyperaceae) and legumes (Fabaceae). Grass cover, species 9 richness (2 m2 scale), total density, and the density of each life form group, were 10 analysed with a repeated measures ANOVA to assess change over time within fire 11 regime treatments. Factors comprised Fire Regime (fixed factor), Block (random 12 factor), Transect (random factor) and Year (repeated measure, fixed factor), with 13 transect nested in fire regime and block. 14 Patterns of woody and grass-layer species composition at the site were analysed 15 separately with ordination and ANOSIM (analysis of similarity) using Primer 6 (Clarke 16 & Gorley 2001). The number of individuals within each sampling unit was converted to 17 a density value, and then averaged to the compartment level. Data were log (x+1) or 18 square root transformed, and resemblance matrices constructed using the Bray-Curtis 19 measure of similarity. One-way ANOSIM was used to determine whether species 20 composition differed between fire regime treatments and blocks at the end (2007) of the 21 sampling period. 22 Canopy, grass and litter cover prior to fire were compared between blocks using 23 one-way ANOVA. Tukey’s HSD posthoc tests were used to identify which blocks were 24 significantly different from each other. Litter and Canopy cover was compared between 12 1 fire regime treatments in each survey year (2004-2007 inclusive) using a non-parametric 2 Kruskal-Wallis Ranks test (as a result of positively skewed distributions with zero 3 counts). 4 5 RESULTS 6 The two woody vegetation surveys recorded 84 plant species at a mean density 7 of 9945 stems ha-1. Approximately half of the species (34 spp.), and a third of the stems 8 overall (36.1%), were rainforest-associated species. Rainforest species had a greater 9 prominence in the grass-layer strata, where they comprised 38.5% of the stems, 10 compared to 19.1% in the midstorey and 15.9% in the overstorey. The overwhelming 11 majority of woody stems were sprouts (73.6%), followed by small saplings (12.6%), 12 stems in the lower (7.1%) and upper (3.5%) midstorey, large saplings (2.4%) and 13 overstorey trees (0.7%). 14 Site characteristics varied along the edaphic gradient prior to the reintroduction 15 of fire. Canopy cover in 2004 was significantly lower in the North (58.6 ± 3.8%) and 16 Central (57 ± 2.1 %) blocks compared to the South block (80.6 ± 5.3; F2,15 =9.5, 17 P<0.01). Similarly litter cover in 2004 was significantly lower in the North (22.8 ± 18 2.5%) and Central (23.6 ± 4.5%) blocks compared to the South block (43.4 ± 3.9; F2,15 19 =9.62, P<0.01). 20 Edaphic conditions were clearly an important influence on vegetation with a 21 significant effect of block overall evident for abundance of sprouts, small saplings, and 22 stems in the upper midstorey and overstorey (Table 2). In all cases, abundance of these 23 groups was higher in the South block (sandy loams soils) than the North or Central 13 1 block (lateritic soils; Fig. 2). This trend was particularly evident for rainforest species 2 in these strata (Table 2). 3 Leaf litter cover increased at a similar rate between 2003 and 2007 in the Annual 4 Early, Early 04 and Unburnt regimes, demonstrating that lower fire intensities had little 5 effect on litter accumulation (Fig. 4). By contrast, the higher intensity fire in the Early 6 05 regime in 2006, and late dry season fires in the Biennial Late regime, markedly 7 reduced litter cover (Fig. 4). As a consequence, there was a significant difference in 8 litter cover between regimes in 2005 (H(4, n=30) = 15.25, P<0.01) and 2007 (H(4, n=30) = 17.34, 9 P<0.01), after some of the higher intensity fires. 10 Canopy cover remained relatively constant during the study period in the 11 Unburnt regime and in regimes with low-intensity or infrequent fires (i.e. Annual Early 12 and Early 04 regimes, Fig. 4). In contrast, canopy cover decreased (by approximately 13 15%) in the first year after the moderate-intensity Early 05 fire, and decreased (by 14 approximately 10%) in response to the second fire of the Biennial Late regime, the most 15 intense fire of the experiment, resulting in a significant difference in canopy cover 16 between regimes in 2007 (H(4, n=30) = 10.92, P<0.05). 17 Effect of fire on woody plant species 18 The reintroduction of fire altered stand structure through reduced abundance of 19 the large saplings and the lower midstorey stems, although the magnitude of this change 20 was determined by fire intensity resulting in Year by Regime interactions (Table 2). 21 The abundance of large sapling stems increased or remained relatively stable within 22 most plots, but decreased with the repeated moderate-intensity fires of the Biennial Late 23 regime (Fig. 5). The abundance of lower midstorey stems decreased within the fires of 14 1 the Early 05 and Biennial Late regimes but was maintained or increased within the UN, 2 and lower intensity Annual Early and Early 04 regimes (Fig. 5). A similar pattern 3 occurred for the abundance of the tree species, with the latter also demonstrating a 4 significant year by regime interactions (Table 2). The mean abundance of lower 5 midstorey stems overall (and tree species in that group) showed a significant effect of 6 year overall (Table 2). Stems in the upper midstorey strata and overstorey strata were 7 more resilient to burning than those in the lower midstorey, showing no effect of Year, 8 nor a Year x Regime interaction (Table 2; Fig. 5). 9 The fire regimes in this experiment had little effect on the sprout layer with 10 increases in mean sprout abundance evident in all fire regime treatments over time (Fig. 11 6); resulting in a significant Year effect. The significant increase was evident for the 12 rainforest, midstorey and overstorey groups in the sprout layer (Table 2). The 13 abundance of overstorey species in the sprout layer decreased in the higher-intensity 14 Early 05 regime in contrast to the increases observed within all other fire regimes, 15 resulting in a significant year by regime interaction (Table 2). The abundance of small 16 saplings also showed a significant increase over time, primarily due to the increase in 17 overstorey species (Table 2). 18 Despite the reduction in large saplings and midstorey stems, basal area was not 19 significantly affected by year, fire regime or block. However, there were significant 20 inter-annual changes in species richness and stem abundance. Species richness 21 increased significantly over time within all but the Biennial Late regime (Fig. 7), from a 22 mean of 24.7 species/transect in 2004, to 27.1 species/transect in 2007. There was no 23 significant effect of Fire regime or Block, or interaction between factors (Table 2). 15 1 Rainforest species showed a significant effect of block in all strata (Table 2), 2 reflecting the general pattern observed for the community (Fig. 3). Rainforest species in 3 the sprout and lower midstorey strata showed a significant effect of Year, whereby 4 sprouts increased considerably over time (particularly in the Annual Early and UN 5 regimes) and lower midstorey stems decreased, particularly in the higher intensity 6 regimes (Early 05 and Biennial Late; Fig. 8). 7 The abundance of standing dead stems tripled during the experiment, resulting 8 in a significant effect of Year and Block but no difference was detected between fire 9 regimes (Table 2). The mean number of dead stems was 38.5 stems ha-1 in 2004 10 (ranging from 82.3 ha-1 in the South block to 10.4 stems ha-1 in the North block), to a 11 mean of 160.8 dead stems ha-1 in 2007 (ranging from 259.4 stems ha-1 in the South 12 block to 59.4 stems ha-1 in the North block). Despite this increase, standing dead stems 13 (represented by 59 species) still only comprised a small proportion (1%) of the total 14 number of stems in 2007. The species with the highest proportion of standing dead 15 stems (with n > 10 individuals in total) were Calytrix exstipulata (Myrtaceae; 8.0%), 16 Acacia holosericea (Mimosaceae; 7.6%) and Acacia oncinocarpa (Mimosaceae; 7.3%) 17 in the grass-layer, and Grevillea pteridifolia (Proteaceae; 44.4%), Alphitonia excelsa 18 (Rhamnaceae; 42.2%) and Acacia latescens (Mimosaceae; 23.1%) in the midstorey / 19 overstorey. 20 In 2007, basal sprouting was evident on 59 species. Basal sprouting occurred on 21 7.4% of the woody stems, and the vast majority of basal sprouters (live tissue 22 component, rather than the often taller dead portion) were present in the grass-layer 23 (87.0%), some were present in the midstorey layer (12.8%), and only two individuals 24 were present in the overstorey layer (0.1%). There was a significant difference between 16 1 fire regime treatments in the proportion of individuals that exhibited basal sprouting (F 2 4, 10 3 the UN and the Early 04 regime (approximately 4 % of total stems in the regimes), 4 whereas basal sprouters comprised about 10% of stems in the Annual Early and Early 5 05 regimes, and 20% in the higher-intensity Biennial Late regime. The species for 6 which basal sprouters were most common (with n > 10 individuals in total) included 7 Lophostemon lactifluus (Myrtaceae; 43.8%), Syzygium eucalyptoides subsp. bleeseri 8 (Myrtaceae; 32.6%) and Corymbia bleeseri (Myrtaceae; 31.4%) for stems in the grass- 9 layer, and Syzygium eucalyptoides subsp. bleeseri (Myrtaceae; 40.0%), Acacia = 6.98; P < 0.01; Fig. 9). The lowest proportion of basal sprouting was observed in 10 lamprocarpa (Mimosaceae; 28.1%) and Acacia mimula (Mimosaceae; 19.6%) in the 11 midstorey / overstorey. 12 Effects on grass and legumes species 13 In total, 116 herbaceous grass-layer plant species were identified during the 14 annual surveys, at a mean density of 70 individuals m-2 (range: 0-1802 m-2). The grass- 15 layer was dominated by annual grasses (36.9% of total density, 17 spp.) and annual 16 forbs (28.9%, 56 spp.). Sedges (14.9%, 11 spp.), perennial grasses (13.4%, 18 spp.) and 17 legumes (5.9%, 14 spp.) together comprised the remaining one-third of the grass-layer 18 plants. The most common species (according to density) were the annual grass 19 Pseudopogonatherum contortum (Poaceae; 11.7 ± 1.1 SE individuals per m-2), sedge 20 Fimbristylis sp. Charles Darwin (Cyperaceae; 7.0 ± 0.7 per m-2) and perennial grass 21 Eriachne triseta (Poaceae; 6.1 ± 0.3 per m-2). 22 Grass cover was low at the commencement of the fire experiment (12.9 ± 2.2 23 %), and with a small overall increase to 18.0 ± 2.3 after the three years for experimental 24 fires). There was no significant increase in grass cover was in the higher intensity 17 1 Biennial Late regime, where cover increased from 11.7 ± 7.1 % in 2004 to 27.9 ± 3.1 2 %). However, the total density of individuals in the grass layer, and the density of 3 annual grasses and legumes specifically, showed a significant effect of year overall, 4 with density increasing over time (Table 3). The mean density of all life forms, except 5 sedges, progressively increased in each year of the experiment (Table 4). For both 6 annual grasses and total density, a significant year by regime interaction was not evident 7 (Table 3), indicating that these increases occurred in spite of the different fire regime 8 treatments applied. Total grass-layer density and annual grass density increased by ~ 9 45% and 100% respectively over the sampling period. The increase in density during 10 the experiment coincided with an increase in the total rainfall and number of rain days 11 during the wet season (November to February inclusive) prior to sampling (from 2005 12 onwards; Table 4). The increase in grass-layer plant density also coincided with a 13 decrease in the total rainfall and number of rain days in the late dry season (August to 14 October inclusive) prior to sampling. Species richness and the density of all life form 15 groups showed a significant effect of block overall (Table 3). 16 Perennial grasses and legumes showed a significant year by regime interaction, 17 with changes in density over time determined by the fire regime treatments imposed 18 (Table 3). Perennial grass density increased gradually in the UN regime, to double 19 within three years. A similar pattern occurred within the Early 04 regime which was 20 burnt early in the experiment but protected from fire thereafter (Fig. 10 a). Regular 21 burning and regimes with higher-intensity fires tended to limit population increases; 22 perennial grass density in such regimes (e.g. Biennial Late) remained relatively stable 23 during the experiment. 18 1 The significant regime × year interaction of legume density was the result of a 2 substantial increase in density in the first year after fire, particularly after the higher- 3 intensity fires of the Early 05 and Biennial Late regimes (Fig. 10 b). The post-fire 4 response was smaller after lower-intensity fires of the Annual Early regime, and no 5 increases were evident in the UN regime. Indeed, density tended to decrease after the 6 initial post-fire increases (e.g. Fig. 10 b, Early-05, 2006-2007). 7 Density of the annual grass Sorghum intrans was low (≤ 0.06 individuals per 8 m2) before burning commenced and remained so three years after the application of fire 9 regime treatments. There was, however, a notable increase in abundance towards the 10 end of the experiment, apparently unrelated to fire regime (Fig. 10 c). 11 Floristic composition 12 After three years, there was no significant difference in woody (R = -0.079; P = 13 0.721) or grass-layer (R = -0.119, P = 0.756) species composition between fire regime 14 treatments. At this time, however, there was a strong, significant effect of Block 15 evident for both components, with the open forest community in the North and Central 16 Block separating from the South Block (woody plants, R = 0.369, P = 0.002; grass-layer 17 plants, R = 0.606, P = 0.002). This dissimilarity of the woody component was primarily 18 due to the differences in midstorey taxa (Table 5), with five of the six species that 19 contributed more than 3% towards the dissimilarity being midstorey species, and the 20 other being a canopy tree. Four of these species (including three rainforest species) were 21 most common in the ecotone community of the South Block. Most grass-layer plant 22 species that contributed to the dissimilarity in composition between the woodland 23 (North and Central Blocks) and ecotone (South Block) were more abundant in the 19 1 woodland, with the notable exception of E. triseta, which was much more abundant in 2 the ecotone (Table 5). 3 DISCUSSION 4 Vegetation structure after long-term fire exclusion 5 The high density of woody plants at the site at the commencement of the study, 6 particularly in the lower midstorey, suggests that a significant change in vegetation 7 structure occurred during the period of fire exclusion. The density of woody sprouts 8 and small saplings at the site (500 – 1000 stems/ha) was considerably higher than 9 typically occurs in Australian savannas (Bowman & Panton 1993; Bowman & Panton 10 1995; Vigilante & Bowman 2004). In contrast, woody stem density was similar to that 11 at another Australian savanna site (Munmarlary; 1400 midstorey stems/ha) where fire 12 had been excluded for >20 years, and where midstorey stems were previously absent 13 (Russell-Smith et al. 2003). At our site, edaphic characteristics strongly affected the 14 development of the midstorey, with significantly higher midstorey stem density in the 15 South block, where the sandy loam soils with low gravel content are likely to result in 16 higher plant available moisture compared with the North and Central blocks where the 17 soils are clay loams with high gravel content and a shallow lateritic layer (Scott et al. 18 2009). The South block is located close to a permanent freshwater lagoon, and therefore 19 the watertable is likely to be higher throughout the dry season than on the ridge in the 20 north (Scott et al. 2009). This process of woody thickening is also supported by the 21 trends documented in the unburnt regime during the three years of this study, when stem 22 density increased throughout most strata. An increase in the midstorey during this 23 period is consistent with previous fire experiments in the savannas of northern 24 Australia, Africa and South America, where fire exclusion has enabled the recruitment 20 1 of woody sprouts into the midstorey above the flame zone and away from the 2 competition of grass-layer plants (Brookman-Amissah et al. 1980; San Jose & Farinas 3 1991; Silva et al. 2001; Russell-Smith et al. 2003; Woinarski et al. 2004). 4 A large proportion of individuals at the site, particularly in the grass-layer, were 5 rainforest species; however, it is unclear the extent to which this reflects species 6 composition prior to fire exclusion, or recruitment after. Many of the most common 7 rainforest-associated species (e.g. Alphitonia excelsa and Acacia lamprocarpa) also 8 occur in frequently burnt savannas, but show a strong association with fire protection 9 over time (Bowman & Panton 1995; Woinarski et al. 2004). Edaphic characteristics 10 were critically important in the composition of the midstorey, with substantially higher 11 representation of rainforest species in the South block. These edaphic characteristics 12 clearly influenced the structure and composition of the midstorey within this study site 13 irrespective of fire regime. 14 Grass-layer plant composition has changed during the period of long term fire exclusion 15 since 1981. Vegetation surveys at the time showed that the grass-layer was dominated 16 by the tall annual Sarga intrans (Andersen and Hoffmann 2010), and this species 17 remains dominant immediately outside the Park (Scott 2008). Therefore, in the long- 18 term absence of fire, the composition of the grass understory changed to dominance by 19 the smaller perennial Eriachne triseta and low-growing annuals (e.g. 20 Pseudopogonatherum contortum and Eriachne agrostidea). This change has been 21 documented elsewhere in northern Australia (Russell-Smith et al. 2003; Williams PR et 22 al. 2003a; Woinarski et al. 2004). Mechanistic experiments have independently shown 23 the negative influence of canopy shading and litter cover on the establishment and seed 24 production of grasses, including annual Sorghum (Cook et al. 1998; Scott et al. 2010a). 21 1 Effects of re-introducing fire 2 Reintroducing moderate intensity fire to the long-unburnt savanna caused the 3 vegetation to transition back toward open savanna structure, principally through a 4 reduction in lower midstorey stem density. However, fire intensity was the important 5 factor determining this change. The abundance of midstorey stems increased or 6 remained similar in the UN regime and with repeated, low-intensity fire (Annual Early, 7 mean fire intensity = 116 kW m-1), demonstrating woody plant density had not peaked 8 during the 20 years of fire exclusion. By contrast, the abundance of large saplings 9 decreased in the Late regime (808 and 2121 kW m-1), and the number of lower 10 midstorey stems decreased considerably in both the Biennial Late and Early 05 regime 11 (1110 kW m-1). At least part of this reduction in midstorey stem density could be 12 attributed to top-kill, where individuals remained alive, but were redistributed into 13 smaller size classes (e.g. sprouts). Minimal changes to basal area were recorded (Table 14 2), indicating that most structural change involved individuals with thin stems (i.e. in 15 the midstorey). 16 In general, the resilience of sprouts and small saplings, response of midstorey 17 stems and resistance of overstorey trees to fire in this study is well known in northern 18 Australia. The susceptibility of midstorey stems to burning is consistent with the notion 19 that smaller (lower dbh) plants experience greater rates of stem mortality (top-kill) and 20 have a lower probability of surviving fire (Williams et al. 1999; Werner et al. 2006; 21 Liedloff & Cook 2007; Prior et al. 2009). The fire intensities in this study were low. 22 Had fire intensities been higher (e.g. > 5000 kW m-1), mortality and top-kill may have 23 been much higher and basal area significantly reduced, as a consequence of mortality in 22 1 large overstorey trees, as evidenced by the landscape-scale Kapalga fire experiment 2 (Williams et al. 1999). 3 The reintroduction of fire had a limited effect on grass-layer species 4 composition. Populations of grass-layer plants are highly resilient to fire as a 5 consequence of vegetative regeneration and the annual input of seeds into fire- 6 protected, dry season seed banks (O’Connor 1985; Belsky 1992; O’Connor & Pickett 7 1992; Scott et al. 2010b, Williams PR et al. 2003a). The increase in perennial grass 8 density in compartments remaining unburnt for several years appears to reflect 9 increased seedling recruitment in the absence of fire because there was a uniformly low 10 mortality rate of common, adult perennial grasses recorded at the site in the unburnt 11 regime and the burnt regimes (Scott et al. 2010c). Enhanced seedling recruitment in the 12 unburnt areas could have resulted from an increase in seed production or seed viability, 13 or lower rates of seed or seedling mortality, in those compartments. The marked 14 increase in legume density after fire (Fig. 10 b) was a direct result of elevated 15 temperatures in the soil overcoming physical dormancy (Auld & O’Connell 1991, 16 Morrison et al. 1998, Bell & Williams 1998). The response of seeds to heating has been 17 documented for several legume species in north-eastern Australia, including those at the 18 study site (Scott et al. 2010b; Williams PR et al. 2003b). However, even though 19 perennial grasses and legumes showed significant changes in abundance over time 20 according to fire regime, these life forms were only minor components of the vegetation 21 and so did not significantly influence composition overall. 22 Inter-annual rainfall variability was strongly related to the density of grass-layer 23 plants (Table 4). Rainfall variability has been demonstrated to affect grass-layer species 24 composition in Australia, Africa and South America, particularly in response to long 23 1 periods of drought (O’Connor 1994; O’Connor & Everson 1998; Williams RJ et al. 2 2003; Bisigato & Bertiller 2004). Seedlings of grass-layer plants are vulnerable to 3 desiccation, so mortality would be exacerbated in drier wet seasons with prolonged 4 periods of low rainfall intensity or frequency (Andrew & Mott 1983; O’Connor 1994, 5 1996; Veenendaal et al. 1996). A later start to the wet season, as occurred in 2006- 6 2007, could conversely enhance seedling survival by allowing most seed germination to 7 occur closer to the mid-wet season period of reliable rainfall, thereby minimising post- 8 germination mortality (Crowley & Garnett 1999). 9 10 Reversibility of change following long-term fire exclusion State-and-transition conceptual models accommodate the possibility that 11 different stable vegetation states are possible within a location, and vegetation state 12 persists unless there is a significant change in an ecosystem driver (Bestelmeyer et al. 13 2003; Briske et al. 2006). An open savanna state is maintained by burning. Fire results 14 in some mortality of midstorey and overstorey stems, and the suppression of woody 15 sprouts in the grass-layer. The open canopy supports a continuous grass-layer, which, 16 due to an annual cycle of fuel accumulation and curing, reinforces the regime of 17 frequent burning (Williams et al. 2003). At the beginning of this study, the savanna 18 structure reflected the effects of a long period of fire exclusion with a well-developed 19 midstorey (2-10m) represented by ~500 – 1000 stems/ha. In our study, fires of 20 moderate intensity (1000 – 2000 kWm-1) reduced the density of the midstorey, with a 21 decrease from 800 to 350 stems/ha in the Biennial Late regime, and from ~1100 to 650 22 stems/ha in the Early 05 regime. This study showed, however, that fire intensity and 23 frequency are critical to these processes, with a maintenance (at ~1000 stems / ha) of 24 the midstorey in the UN regime and increased development of the midstorey (from 24 1 ~500 to 750 stems/ha) in the low intensity Annual Early regime. Moderate or high 2 intensity fires over several years may be required to achieve the sparse midstorey 3 characteristic of frequently burnt savanna in this region, and fully restore the structure 4 and / or function of frequently burnt communities typical of the region’s Eucalyptus 5 dominated savanna. 6 The reduction in lower midstorey stems in higher-intensity regimes was driven 7 by the responses of canopy and rainforest species in particular (Table 2). It would 8 seem, therefore, that both canopy and rainforest species experience a recruitment 9 bottleneck in the lower midstorey in these regimes. For the short-term at least, this has 10 little consequence for canopy species because of the resilience and longevity of trees 11 already in the overstorey. Over the longer-term however, fire-free periods will be 12 required for recruitment into the overstorey (Setterfield 2002). Similarly, fire-free 13 periods would be required for rainforest species to dominate the midstorey and 14 overstorey vegetation layers. 15 Even if the flammable annual grass Sorghum intrans was replaced by less 16 flammable grasses such as the perennial Eriachne triseta, the changes still indicate the 17 presence of a pre-threshold state across most of the experimental catchments because 18 the resulting composition carried fires that modified vegetation structure (positive 19 feedbacks). This change in composition and others, such as the increase in perennial 20 grass density with fire exclusion (Fig. 10 a) and increase in legume density after fire 21 (Fig. 10 b), did not cause any structural or functional degradation of the ecosystem. In 22 the time-frame examined, the loss of S. intrans was unidirectional; the reintroduction of 23 fire was unable to reinstate its former dominance. This suggests some longer-term 24 implications, given poor seed dispersal and a lack of long-term seed dormancy within 25 1 the seed bank (Andrew & Mott 1983). The immediate (1 year) post-fire increase in the 2 abundance of Sorghum ecarinatum Lazarides reported on the Arnhem Plateau was 3 probably a result of its occurrence in over half of the quadrats prior to the fire (Russell- 4 Smith et al. 2002). 5 A small proportion of the site (perhaps 10%) appears to have already crossed an 6 ecological threshold, characterised by the predominance of positive feedbacks which 7 reinforce a change between open woodland and closed forest states. There were patches 8 in the South Block (approximately 1000 m2) where woody plant density was very high 9 (~ 10 500 woody stems m-1), presumably a result of more favourable edaphic conditions 10 supporting this midstory development in the absence of fire (Bowman & Minchin 1987; 11 Medina & Silva 1990; Scott et al. 2009). As a consequence, the higher canopy shading 12 and higher amounts of litter (see Fig. 8) resulted in a fuel layer dominated by litter with 13 significantly fewer herbaceous grass-layer plant species (Scott et al. 2009). These litter 14 patches acted as a network of fire breaks, and despite efforts to re-ignite them, were 15 unable to sustain fire. A dominance of litter in the fuel layer therefore reinforced 16 woody thickening on a transition towards closed forest. 17 18 19 Acknowledgements 20 Funding for this project, including a PhD scholarship for K. Scott and a postdoctoral 21 fellowship for C. Parr, was provided by the Bushfire Cooperative Research Centre. A 22 grant from the Bushfire CRC allowed the K. Scott to write this manuscript as a 23 Research Associate at Charles Darwin University. Many people assisted with the 26 1 conduct of experimental fires, vegetation sampling, plant identification and data 2 analysis. We particularly thank Bushfires NT, Jenni Low Choy, Ian Cowie, the late 3 Jack Cusack, Robert Eager, Jon Edgar, Andrew Edwards, Garry Fox, John Henderson, 4 Anthony Kerr, Wawan Sastrawan Manullang, Keith McGuiness, Sally Mitchell, Ken 5 Nicholl, Dave Pearson, Liz Smith and Dick Williams. The Territory Wildlife Park 6 shared our vision towards creating the site as an educational and research facility, and 7 allowed us to access it regularly. 27 1 REFERENCES 2 Andersen, A.N., Cook G.D, Corbett L.K., Douglas M.M., Eager R.W., Russell-Smith J., 3 Setterfield S.A., Williams R.J., & Woinarski J.C.Z. (2005) Fire frequency and 4 biodiversity conservation in Australian tropical savannas: implications from the 5 Kapalga fire experiment. Austral Ecol. 30, 155–167. 6 Andrew M.H. & Mott J.J. (1983) Annuals with transient seed banks: the population 7 biology of indigenous Sorghum species of tropical north-west Australia. Aust. J. 8 Ecol. 8, 265-276. 9 10 Angassa A. & Oba G. (2008). Effects of management and time on mechanisms of bush encroachment in southern Ethiopia. African J. Ecol. 46, 186-196. 11 Augustine D.J. & McNaughton S.J. (2004) Regulation of shrub dynamics by native 12 browsing ungulates on East African rangeland. J. App. Ecol. 41, 45-58. 13 14 15 16 Auld T.D. & O’Connell M.A. (1991) Predicting patterns of post-fire germination in 35 eastern Australian Fabaceae. Aust. J. Ecol. 16, 53-70. Bell D.T. & Williams D.S. (1998) Tolerance of thermal shock in seeds. Aust. J. Bot. 46, 221-233. 17 Belsky A.J. (1992) Effects of grazing, competition, disturbance and fire on species 18 composition and diversity in grassland communities. J. Veg. Sci. 3, 187-200. 19 Bestelmeyer B.T., Brown J.R., Havstad K.M., Alexander R., Chavez G. & Herrick J.E. 20 (2003) Development and Use of State-and-Transition Models for Rangelands. J. 21 Range Mgt. 56, 114-126. 22 23 Bisigato A.J. & Bertiller M.B. (2004) Seedling recruitment of perennial grasses in degraded areas of Patagonian Monte. Rang. Ecol. Mgt. 57, 191-196. 28 1 2 3 4 5 6 7 8 9 10 11 12 13 Bond W.J. (2008). What limits trees in C4 grasslands and savannas? Ann. Rev. Ecol, Evol. Syst. 39, 641-659. Bond W.J. & Parr C.L. (in press) Beyond the forest edge: ecology, diversity and conservation of the grass biomes. Biol. Cons. Boone R.B. & Wang G. (2007) Cattle dynamics in African grazing systems under variable climates. J. Arid Env. 70, 495-513. Boulière F. (1983) Tropical Savannas. Ecosystems of the World. Volume 13. Elsevier Publishing, Amsterdam. Bowman D.M.J.S. & Fensham R.J. (1991) Response of a monsoon-forest savanna boundary to fire protection, Weipa, northern Australia. Aust. J. Ecol. 16, 111-118. Bowman D.M.J.S. & Minchin P.R. (1987) Environmental relationships of woody vegetation patterns in the Australian monsoon tropics. Aust. J. Bot. 35, 151-169. Bowman D.M.J.S. & Panton W.J. (1993) Factors that control monsoon-rainforest 14 seedling establishment and growth in north Australian Eucalyptus savanna. J. 15 Ecol. 81, 297-304. 16 Bowman D.M.J.S. & Panton W.J. (1995) Munmarlary revisited: Response of a north 17 Australian Eucalyptus tetrodonta savanna protected from fire for 20 years. Aust. J. 18 Ecol. 20, 526-531. 19 20 21 Briske D.D., Fuhlendorf S.D. & Smeins F.E. (2003). Vegetation dynamics on rangelands: a critique of the current paradigms. J. App. Ecol. 40, 601-614. Briske D.D., Fuhlendorf S.D. & Smeins F.E. (2006). A Unified Framework for 22 Assessment and Application of Ecological Thresholds. Range. Ecol Mgt. 59, 225- 23 236. 29 1 Brookman-Amissah J., Hall J.B., Swaine M.D. & Attakorah J.Y. (1980) A re- 2 assessment of a fire protection experiment in north-eastern Ghana savanna. J. 3 App. Ecol. 17, 85-99. 4 5 6 Clarke K.R. & Gorley R.N. (2001) Primer v5: User manual/tutorial. Primer-E Ltd., Plymouth. Cook G.D., Ronce O. & Williams R.J. (1998) Suppress thy neighbour: The paradox of a 7 successful grass that is highly flammable, but harmed by fire. In: Proceedings 8 13th Conference on Fire and Forest Meteorology, Lorne, October 1996, pp. 71-4 9 International Association of Wildland Fire, Hot Springs, South Dakota. 10 11 12 Crowley C.M. & Garnett S.T. (1999) Seeds of the annual grasses Schizachyrium spp. as a food resource for tropical granivorous birds. Aust. J. Ecol. 24, 208-220. Fensham R.J., Fairfax A.J., Butler D.W. & Bowman D.M.J.S. (2003) Effects of fire and 13 drought in a tropical eucalypt savanna colonized by rain forest. J. Biogeog. 30, 14 1405-1414. 15 16 Holdo R.M. (2007) Elephants, fire and frost can determine community structure and composition in Kalahari woodlands. Ecol. App. 17, 558-568. 17 Hutley, L.B. and Setterfield, S.A. (2008). Savannas. In Jorgensen, SE and Faith, BD 18 (Ed.), Encyclopedia of Ecology. Oxford, U.K.: Elsevier B.V. (pp. 3143-3154). 19 20 21 Lemmon, P.E. (1956). A spherical densiometer for estimating forest overstory density. For. Sci. 2: 314-320. Liddle D.T., Russell-Smith J., Brock J., Leach G.J. & Connors G.T. (1994) Atlas of the 22 vascular rainforest plants of the Northern Territory. Flora of Australia 23 Supplementary Series No. 3. ABRS, Canberra. 30 1 Liedloff A.C. & Cook G.D. (2007) Modelling the effects of rainfall variability and fire 2 on tree populations in an Australian tropical savanna with the FLAMES 3 simulation model. Ecol. Mod. 201, 269-282. 4 Medina E. & Silva J.F. (1990) Savannas of northern South America: a steady state 5 regulated by water-fire interactions on a background of low nutrient availability. 6 J. Biogeog. 17, 403-413. 7 Midgley, J.J., Lawes, M.J. & Chamaille-Jammes, S. (2010). Savanna woody plant 8 dynamics: the role of fire and herbivory, separately and synergystically. Aust. J. 9 Bot. 58, 1-11. 10 Morrison D.A., McClay K., Porter C. & Rish S. (1998) The role of the lens in 11 controlling heat-induced breakdown of testa-imposed dormancy in native 12 Australian legumes. Annals Bot. 82, 35-40. 13 O’Connor T.G. (1985) A synthesis of field experiments concerning the grass-layer in 14 the savanna regions of southern Africa. South African National Scientific 15 Programmes Report 114. CSIR - Commonwealth Scientific and Industrial 16 Research, Pretoria, South Africa. 17 18 O’Connor T.G. (1994) Composition and population responses of an African savanna grassland to rainfall and grazing. J. App. Ecol. 31, 155-171. 19 O’Connor T.G. (1996) Hierarchical control over seedling recruitment of the bunch- 20 grass Themeda triandra in a semi-arid savanna. J. App. Ecol. 33, 1094-1106. 21 O’Connor T.G. & Everson T.M. (1998) Population dynamics of perennial grasses in 22 African savanna and grassland. In: Population Biology of Grasses (ed G.P. 23 Cheplick) pp. 333-365. Cambridge University Press: Cambridge, UK. 31 1 O’Connor T.G. & Pickett G.A. (1992) The influence of grazing on seed production and 2 seed banks of some African savanna grasslands. J. App. Ecol. 29, 247-260. 3 Parr, C.L. & Andersen, A.N. (2008). Fire resilience of ant assemblages in long unburnt 4 5 savanna of northern Australia. Austral Ecol. 33: 830-838. Prior L.D., Murphy B.P. & Russell-Smith J. (2009). Environmental and demographic 6 correlates of tree recruitment and mortality in north Australia savannas. Forest 7 Ecol. Mgt. 257, 66-74. 8 9 10 11 Ratnam, J., Bond, W.J., Fensham, R.J., Hoffmann, W.A., Archibold, S. Lehmann, C.E.R., Anderson, M.T., Higgins, S.I., Sankaran, M. (2011). When is a ‘forest’ a savanna, and why does it matter? Global Ecol. and Biogeog. 20, 653-660. Roques K.G., O’Connor T.G. & Watkinson A.R. (2001) Dynamics of shrub 12 encroachment in an African savanna: Relative influences of fire, herbivory, 13 rainfall and density dependence. J. App. Ecol. 38, 268-280. 14 Russell-Smith J., Ryan P.G. & Cheal D.C. (2002) Fire regimes and the conservation of 15 sandstone heath in monsoonal northern Australia: frequency, interval, patchiness. 16 Biol. Cons. 104, 91-106. 17 Russell-Smith, J., Stanton, P.J., Whitehead, P.J. & Edwards, A. (2004) Rain forest 18 invasion of eucalypt-dominated woodland savanna, Iron Range, north-eastern 19 Australia: I. Successional processes. Journal of Biogeography, 31, 1293–1303. 20 Russell-Smith J., Whitehead P.J., Cook G.D. & Hoare J.L. (2003) Response of 21 Eucalyptus-dominated savanna to frequent fires: Lessons from Munmarlary, 22 1973-1996. Ecol. Mon. 73, 349-375. 32 1 2 3 4 5 6 7 8 9 10 11 San José J.J. & Farinas M.R. (1991) Temporal changes in the structure of a Trachypogon savanna protected for 25 years. Acta Oecol. 12, 237-247. Scheffer, M. & Carpenter, S.R. (2003) Catastrophic regime shifts in ecosystems: Linking theory to observation. Trends Ecol. Evol. 18, 648-656. Scholes R.J. & Archer S.R. (1997) Tree-grass interactions in savannas. Ann. Rev. Ecol. Syst. 28, 517-544. Scott K.A. (2008) The distribution and dynamics of grass-layer plants in a tropical savanna. PhD thesis. Charles Darwin University, Darwin. Scott K.A., Setterfield S.A., Andersen A.N. & Douglas M.M. (2009) Correlates of grass species composition in an Australian savanna. Aust. J. Bot. 57, 10-17. Scott K.A., Setterfield S.A., Douglas M.M. & Andersen A.N. (2010a) Environmental 12 factors influencing the establishment, height and fecundity or the annual grass 13 Sorghum intrans in an Australian tropical savanna. J. Trop. Ecol. 26: 313-322 14 Scott K.A., Setterfield S.A., Douglas M.M. & Andersen A.N. (2010b) Soil seed banks 15 confer resilience to savanna grass-layer plants during seasonal disturbance. Acta 16 Oecol. 36: 202-210 17 Scott K.A., Setterfield S.A., Douglas M.M. & Andersen A.N. (2010c) Fire tolerance of 18 perennial grass tussocks in a savanna woodland. Austral Ecol. 35: 858-861. 19 Setterfield, S.A. (1997) The impact of experimental fire regimes on seed production in 20 21 22 two tropical eucalypt species in northern Australia. Aust. J. Ecol. 22, 279-287. Setterfield, S. A. (2002). Seedling establishment in an Australian tropical savanna: effects of seed supply, soil disturbance and fire. J. of Appl. Ecol. 39:949–959. 33 1 Setterfield S.A., Rossiter-Rachor N.A., Hutley L.B., Douglas M.M. & Williams R.J. 2 (2010) Turning up the heat: the impacts of Andropogon gayanus (gamba grass) 3 invasion on fire behaviour in northern Australian savannas. Diversity and 4 Distributions 16, 854–61. 5 Sharp B.R. & Whittaker R.J. (2003) The irreversible cattle-driven formation 6 transformation of a seasonally flooded Australian savanna. J. Biogeog. 30, 783- 7 802. 8 Silva J.F., Zambrano A. & Farinas M.R. (2001) Increase in the woody component of 9 seasonal savannas under different fire regimes in Calabozo, Venezuela. J. 10 11 Biogeog. 28, 977-983. Sivertsen D.P., McLeod P.J. & Henderson R.L. (1980) Land Resources of the Berry 12 Springs Nature Park and Proposed Development Area. Land Conservation Unit, 13 Northern Territory Government, Darwin. 14 Specht R.L. (1981) Foliage projective cover and standing biomass. Vegetation 15 Classification in Australia (Eds A. N. Gillison & D. J. Anderson), pp. 10–21. 16 CSIRO Publishing, Canberra. 17 18 19 Stringham T.K., Krueger W.C. & Shaver P.L. (2003) State and transition modelling: An ecological process approach. J. Range Mgt. 56, 106-113. Veenendaal E.M., Ernst W.H.O. & Modise G.S. (1996) Effect of seasonal rainfall 20 pattern on seedling emergence and establishment of grasses in a savanna in south- 21 eastern Botswana. J. Arid Env. 32, 305-317. 34 1 Vigilante T. & Bowman D.M.J.S. (2004) Effects of fire history on the structure and 2 floristic composition of woody vegetation around Kalumburu, North Kimberley, 3 Australia: a landscape-scale natural experiment. Aust. J. Bot. 42, 381-404. 4 Walker B.H., Langridge J.L. & McFarlane F. (1997) Resilience of an Australian 5 savanna grassland to selective and non-selective perturbations. Aust. J. Ecol. 22, 6 125-135. 7 Werner P.A. (2005) Impact of feral water buffalo and fire on growth and survival of 8 mature savanna trees: An experimental field study in Kakadu National Park, 9 northern Australia. Austral Ecol. 30, 625-647. 10 Werner P.A., Cowie I.D. & Cusack J.S. (2006) Juvenile tree growth and demography in 11 response to feral water buffalo in savannas of northern Australia: An experimental 12 field study in Kakadu National Park. Aust. J. Bot. 54, 1-14. 13 14 15 Westoby M., Walker B. & Noy-Meir I. (1989) Opportunistic management for rangelands not at equilibrium. J. Range Mgt. 42, 266-274. Williams P.R., Congdon R.A., Grice A.C. & Clarke P.J. (2003a) Effect of fire regime 16 on plant abundance in a tropical eucalypt savanna of north-eastern Australia. 17 Austral Ecol. 28, 327-338. 18 Williams P.R., Congdon R.A., Grice A.C. & Clarke P.J. (2003b) Fire-related cues break 19 seed dormancy of six legumes of tropical eucalypt savannas in north-eastern 20 Australia. Austral Ecol. 28, 507-514. 21 Williams R.J., Cook G.D., Gill A.M. & Moore P.H.R. (1999) Fire regime, fire intensity 22 and tree survival in a tropical savanna in northern Australia. Aust. J. Ecol. 24, 50- 23 59. 35 1 Williams R.J., Müller W.J., Wahren C., Setterfield S.A. & Cusack J. (2003) Vegetation. 2 In: Fire in tropical savannas: the Kapalga experiment (Eds. A.N. Andersen, G.D. 3 Cook & R.J. Williams) pp.79-106. Springer, New York. 4 Woinarski J.C.Z., Risler J. & Kean L. (2004) Response of vegetation and vertebrate 5 fauna to 23 years of fire exclusion in a tropical Eucalyptus open forest, Northern 6 Territory, Australia. Austral Ecol. 29, 156-176 36 Table 1. Burning schedule for experimental fire regimes at the Territory Wildlife Park. Values are mean intensity (kW m-1) with 1 SE in parenthesis (data unavailable for Early 04). – indicates no fire event. Fire regime Year 2004 2005 2006 Annual Early 144.0 (18.1) 77.7 (36.7) 125.3 (54.1) Early 04 n/a - - Early 05 - 1110.3 (569.0) - Biennial Late 808.4 (625.0) - 2121.3 (975.3) 37 Table 2. Summary ANOVA results for variation in species richness, and the abundance of woody plants in each stratum between 2004 and 2007. Basal area, species richness and standing dead stems consider all strata together, whereas abundance variables are categorised by stratum for all species combined (e.g. sprouts, 0-0.49 m), and then by species groups based on the maximum height that they grow too (i.e. groundlayer, shrub, tree species) within each stratum. Abundance was also analysed for rainforest species). Significant F-ratios are indicated by asterisks (*P < 0.05, **P < 0.01, ***P < 0.001). Grass-layer species in the small sapling stratum, overstorey species in the large sapling stratum, and monsoon forest species in the overstorey stratum, were too uncommon for separate analysis. Response variable F-ratios Year (Y) Fire regime (R) Block (B) Y×R Y×B F (1, 10) F (4, 10) F (10, 15) F (4, 10) F (10, 15) Species richness 12.98** 1.49 2.24 2.37 1.46 Standing dead stems 23.91*** 0.66 4.53** 1.27 2.41 11.28** 0.76 5.19** 1.58 1.56 10.52** 0.73 8.13*** 1.81 1.14 1.10 0.31 1.67 2.14 0.44 Shrub 10.62** 0.73 3.41* 1.34 2.40 Tree 9.31* 0.74 4.24** 4.98* 0.07 5.39* 1.04 12.35*** 1.24 3.76* Live stem abundance Grass-layer Sprouts (0.0-0.49 m) Rainforest Ground-layer Small saplings (0.5-1.99 38 m) Rainforest 1.25 0.04 11.84*** 1.58 0.55 Shrub 3.68 0.76 11.16*** 1.22 2.76* 19.27** 1.82 2.62* 2.41 0.50 0.32 0.54 1.96 3.89* 1.70 0.49 0.45 4.57** 2.31 1.08 < 0.01 0.42 2.01 3.42 1.46 7.07* 3.44 2.13 5.19* 5.42** Rainforest 5.27* 0.44 6.32*** 1.30 7.96*** Shrub 4.78 3.56* 1.80 4.45* 3.47* Tree 17.09** 1.45 2.68* 12.64*** 0.37 2.46 0.37 3.38* 2.18 3.44* Rainforest 0.75 0.19 7.03*** 1.77 3.37* Shrub 1.03 0.42 2.01 1.86 5.62** Tree 4.39 1.29 1.13 1.37 0.54 2.50 0.79 3.58* 2.91 1.42 Tree Midstorey Large saplings (> 2.0 m & < 3 cm dbh) Rainforest Shrub Lower Midstorey (2.0-4.99 m) Upper Midstorey (5.0-9.99 m) Overstorey (10.0 m +) 39 Table 3. Summary for the repeated measures ANOVA analysis of grass-layer species richness, total density, and the density of individual life form groups. Statistically significant test statistics are indicated by asterisks (*P < 0.05, **P < 0.01, ***P < 0.001). †Data examined annual surveys between 2004 and 2007 inclusive; all other variables examined data between 2005 and 2007 inclusive; late sampling in 2004 precluded the sampling of some life form group. Response variable Source of variation and F-value Fire regime (R) Block (B) Year (Y) R×B R×Y B×Y R×B×Y F (4, 8) F (2, 15) F (2, 4) F (8, 15) F (8, 16) F (4, 30) F (16, 30) 2.4 74.0*** 4.0 5.9** 2.1 7.9*** 3.0** 2.0 18.1*** 16.1* 1.9 1.2 4.0* 2.2* Annual grasses 0.9 28.3*** 27.9** 4.1** 1.8 6.4*** 3.6** Perennial grasses† 2.2 77.9*** F (3, 6) = 4.1 4.4** F(12,24) = 2.2* F(6,45) = 3.2* F24,45=1.5 Forbs 2.3 27.8*** 1.1 2.5 2.2 2.6 0.7 Sedges 1.7 18.3*** 5.7 1.9 0.9 2.6 2.4* Legumes 3.7 34.9*** 8.0* 4.1** 3.9** 4.8** 3.7** Species richness Abundance Total 40 Table 4. Inter-annual variation in the abundance of grass-layer plants (categorised by life form; fire regimes pooled) and rainfall in the late dry season (August to October inclusive) and wet season (November to February inclusive). Abundance data for most life-form groups are unavailable for 2004, as many annual species had already cured by the time of sampling. Year 2004 2005 2006 2007 Total - 60.1 (7.0) 62.2 (7.0) 86.7 (11.7) Annual grasses - 16.0 (3.0) 23.8 (4.7) 38.3 (6.0) 7.0 (0.6) 7.7 (0.6) 9.1 (0.8) 11.3 (1.1) Forbs - 17.6 (2.9) 18.5 (3.0) 24.1 (7.0) Sedges - 15.8 (2.8) 6.8 (1.4) 8.5 (1.9) Legumes - 3.0 (0.5) 3.9 (0.7) 5.4 (0.7) 64 101 33 2 1050 581 789 1115 Late dry season prior 7 6 5 1 Wet season prior 80 62 67 78 Abundance (individuals m-2 ± 1 SE) Perennial grasses Rainfall (mm) Late dry season prior† Wet season prior‡ Rainfall (no. rain days) †Rainfall collected within 2 km of the study site. ‡Rainfall at the site was not consistently recorded in one wet season (2007) and so data presented are from Middle Point AWS 36 km to the north-east (Bureau of Meteorology, unpublished data). Mean rainfall at Middle Point AWS (for the period 19572004) is 71 mm on 9 rain days in the late dry season period, and 965 mm on 72 rain days for the wet season period (Bureau of Meteorology, unpubl. Data) 41 1 Table 5. Mean density (untransformed) of woody and grass-layer plant species in 2007, in each block at the Territory Wildlife Park study site, which 2 contributed at least 3% of the dissimilarity in species composition between blocks. Species in bold font were most abundant in the ecotone community 3 (Southern block). Note separate SIMPER analyses were conducted for woody plants and grass-layer plants. † Monsoon forest species (Liddle et al. 1999) Species Structural / lifeform Mean density by block Dissimilarity Cumulative group (stems ha-1, woody plants; (%) dissimilarity individuals ha-1, grass-layer plants) North Centre South (woodland) (woodland) (ecotone) (%) Woody plants Exocarpus latifolius† Shrub 319.4 4,586.8 13,506.9 10.05 10.05 Petalostigma pubescens Shrub 125.0 17.4 3,208.3 5.82 15.87 Calytrix exstipulata Shrub 1,975.7 2,885.4 2,006.9 4.93 20.8 Alphitonia excelsa† Shrub 152.8 159.7 2,534.7 4.77 25.57 Petalostigma quadriloculare Shrub 1,034.7 1,590.3 1,423.6 3.76 29.33 Syzygium eucalyptoides† Tree 194.4 27.8 1,197.9 3.09 32.43 42 Grass-layer plants Pseudopogonatherum Annual grass 28,2824.1 15,9861.1 0.0 contortum Thaumastochloa major Eriachne triseta Schizachyrium fragile Fimbristylis sp. Charles 10.5 Annual grass 9,8611.1 64,907.4 3,379.6 5.64 16.15 Perennial grass 20,277.8 5,7731.5 139,675.9 5.45 21.6 Annual grass 98,518.5 6,296.3 3,796.3 5.18 26.77 Sedge 65,972.2 60,000.0 1,574.1 3.72 30.49 3.65 34.14 Darwin (JL Egan 5300) Spermacoce exserta 10.5 Herb 9,351.9 106,064.8 79,398.2 43 1 Figure captions 2 Fig. 1. The Territory Wildlife Park study site with savanna with a Eucalyptus dominated 3 overstory and a (a) moderate to (b) well-developed midstorey. 4 Fig. 2. Total abundance of woody stems within each height stratum (sprouts, < 0.5 m; 5 small saplings, 0.5-1.99 m; large saplings, > 2 m but < 3 cm dbh; lower midstorey, 2.0- 6 4.99 m; upper midstorey, 5.0-9.99 m; overstorey, 10.0+ m), in 2004 before regime 7 treatments were first implemented (white bars), and three years later in 2007 (black 8 bars). Data are presented by block: a) North (woodland), b) Centre (woodland), and c) 9 South (ecotone). Data are mean stem density ± SE, where transect (n = 10) is the unit of 10 replication. Note differences in scale. 11 Fig. 3. Total abundance of woody stems of monsoon forest species (sensu Liddle et al. 12 1994) within each height stratum (sprouts, < 0.5 m; small saplings, 0.5-1.99 m; large 13 saplings, > 2 m but < 3 cm dbh; lower midstorey, 2.0-4.99 m; upper midstorey, 5.0-9.99 14 m), in 2004 before regime treatments were first implemented (white bars), and three 15 years later in 2007 (black bars). Data are presented by block: a) North (woodland), b) 16 Centre (woodland), and c) South (ecotone). Data are mean stem density ± SE, where 17 transect (n = 10) is the unit of replication. Rainforest species were very uncommon in 18 the overstorey stratum and data are note presented. Note differences in scale. 19 Fig. 4. Mean (+ 1 SE) values of a) litter cover and b) canopy cover, during the four 20 years (2004-2007) of the manipulative fire experiment. Arrows and associated numbers 21 indicate the timing of fires for each treatment and their intensity, and ‘n/a’ indicates that 22 fire intensity was not assessed (kW m-1 ± 1 SE; see Table 1). 23 Fig. 5. Total abundance of stems within each height stratum above the sprout layer (small 24 saplings, 0.5-1.99 m; large saplings, > 2 m but < 3 cm dbh; lower midstorey, 2.0-4.99 25 m; upper midstorey, 5.0-9.99 m; overstorey, 10.0+ m), in 2004 before regime treatments 26 were first implemented (black bars), and three years later in 2007 (white bars). Data are 27 presented by fire regime treatment: a) Annual Early, b) Early 05, c) Biennial Late, d) 28 Early 04, and e) UN (see Table 1 for details of intensity and frequency). Data are mean 29 stem density ± SE, where transect (n = 6) is the unit of replication. Note differences in 30 scale. 44 1 Fig. 6. Total abundance of sprouts (0-0.49 m) within each of the fire regime treatments 2 (see Table 1 for details of intensity and frequency) of the experiment, in 2004 before 3 regimes were first implemented (white bars), and three years later in 2007 (black bars). 4 Fig. 7. Species richness of all woody species within each of the fire regime treatments 5 (see Table 1 for details of intensity and frequency) of the experiment, in 2004 before 6 regimes were first implemented (white bars), and three years later in 2007 (black bars). 7 Fig. 8. Total abundance of woody stems of monsoon forest species (sensu Liddle et al. 8 1994) within each height stratum (sprouts, < 0.5 m; small saplings, 0.5-1.99 m; large 9 saplings, > 2 m but < 3 cm dbh; lower midstorey, 2.0-4.99 m; upper midstorey, 5.0-9.99 10 m; overstorey, 10.0+ m), in 2004 before regime treatments were first implemented 11 (white bars), and three years later in 2007 (black bars). Data are presented by fire 12 regime treatment: a) Annual Early, b) Early 05, c) Biennial Late, d) Early 04, and e) UN 13 (see Table 1 for details of intensity and frequency). Data are mean stem density ± SE, 14 where transect (n = 6) is the unit of replication. Note differences in scale. 15 Fig. 9. The abundance of basal resprouts in 2007 as a proportion of live stems within 16 each of the fire regime treatments (Annual Early, Early 05, Biennial Late, Early 04 and 17 UN). 18 Fig. 10. Mean (+ 1 SE) density of (a) perennial grasses, (b) legumes, and (c) Sorghum 19 intrans, categorised by fire treatment and year. Perennial grasses and legumes showed a 20 significant regime × year interaction in the experiment (see Table 1), whereas S. intrans 21 did not. Arrows and associated numbers indicate the timing of fires for each treatment 22 and their intensity (kW m-1 ± 1 SE; R. Williams, unpublished data; see Table 1). 45 Figure 1 a) b) 46 Figure 2 a) Sprouts North Centre South Centre South b) Small saplings North c) Large saplings North Centre South 47 d) Lower midstorey North Centre South e) Upper midstorey North Centre South Centre South f) Overstorey North 48 Figure 3 a) Sprouts North Centre South b) Small saplings North c) Large saplings Centre South 49 North Centre South d) Lower midstorey North Centre South e) Upper midstorey North Centre South 50 Figure 4 a) 100 90 125 ± 54 78 ± 37 80 70 60 1110 ± 569 144 ± 18 2121 ± 975 808 ± 625 50 40 30 20 Litter cover (%) 10 0 04 05 06 07 04 05 Annual early 06 07 04 Early 05 05 06 07 Biennial late 100 90 80 70 n/a 60 50 40 30 20 10 0 04 05 06 07 04 05 Early 04 b) 85 06 07 Unburnt 144 ± 18 1110 ± 569 78 ± 37 80 125 ± 54 75 808 ± 625 70 2121 ± 975 65 60 55 50 Canopy cover (%) 45 40 04 05 06 07 04 05 Annual early 06 07 Early 05 80 n/a 70 65 60 55 50 45 40 04 05 Early 04 06 07 04 05 06 Unburnt 05 06 Biennial late 85 75 04 07 07 51 Figure 5 a) Annual Early b) Early 05 52 c) Biennial Late d) Early 04 e) UN 53 Figure 6 54 Figure 7 55 Figure 8 a) Sprouts b) Small saplings c) Large saplings 56 d) Lower midstorey e) Upper midstorey 57 Figure 9 25 Proportion of live stems (%) 20 15 10 5 0 Annual Early Early 05 Biennial Late Early 04 Unburnt 58 Figure 10 20 125 ± 54 Perennial grass density (individuals m -2 ) a) 15 78 ± 37 144 ± 18 1110 ± 569 808 ± 625 2121 ± 975 04 06 10 5 04 05 06 07 04 Annual Early 05 06 07 Early 05 Biennial Late 20 15 10 5 04 05 Early 04 06 07 04 05 05 06 Unburnt 07 07 59 b) 14 2121 ± 975 12 10 125 ± 54 808 ± 625 Legume density (individuals m -2) 8 144 ± 18 78 ± 37 6 1110 ± 569 4 2 0 05 06 07 05 06 Annual Early 07 05 06 Early 05 07 Biennial Late 14 12 10 8 6 4 2 0 05 06 07 05 06 Early 04 c) 07 Unburnt 0.6 0.5 Sorghum intrans density (individuals m -2) 0.4 0.3 0.2 144 ± 18 125 ± 54 78 ± 37 1110 ± 569 808 ± 625 2121 ± 975 0.1 0.0 04 05 06 07 04 Annual Early 05 06 07 Early 05 0.5 0.4 0.3 0.2 0.1 0.0 05 06 Early 04 07 04 05 06 Unburnt 05 06 Biennial Late 0.6 04 04 07 07
© Copyright 2026 Paperzz