Short-term vegetation responses to the first prescribed burn in an urban pine rockland preserve

Background: Preserving fire-dependent ecosystems can mitigate biodiversity loss from urbanization, but prescribing fire is challenging near human habitation. Consequently, dereliction of fire-dependent forests is widespread in urban fragments. Natural disturbance-based management, like prescribing fire, is gaining global acceptance, yet it is unclear what affects prolonged exclusion have on the initial regeneration of isolated plant communities immediately after fire is reintroduced. We took advantage of the first prescribed low-intensity burn on a university pine rockland nature preserve in South Florida, USA, to gain insight. We measured the changes in plant community composition and vegetation cover 1 week before the prescribed burn, and again 1, 2, and 14 weeks after to assess the early and short-term stages of recovery. Results: The fire consumed substantial leaf litter, surface fuels, and canopy leaves, increasing sunlight availability to the understory and exposing bare ground. Many woody plants perished within a week post-burn, particularly invasive shrubs; however, germinating and resprouting plant growth were rapid. By 14 weeks, vegetation covered more of the ground than before the burn, although the upper canopy remained relatively open. Rarefied species richness was recovered by 14 weeks but did not exceed pre-burn levels. Invasive species richness was also maintained post-burn. Despite no overall changes in the community structure, our correspondence analysis and analysis of similarity of the plant community suggest high species turnover from the pre-burn to the final community surveyed, with an intermediate turnover in between. Conclusion: The endangered pine rockland ecosystem, like many fire-dependent ecosystems, is threatened by habitat loss and fire suppression. Managing urban preserves with periodic burns is essential for supporting habitat for endemic species while decreasing demands for manual and time-intensive maintenance. Our study demonstrates that seedling recruitment from early plantings of native species can contribute significantly and immediately to restoration efforts in a fire-excluded urban preserve; however, many changes were ephemeral. Supplemental burns are likely necessary to further reduce vegetation density and sustain changes to the community composition.


Background
Restoring nature preserves in urban areas is a promising avenue for conservation in the modern era (Dunn et al. 2006). If managed with ecological processes as a guiding principle, urban preserves can create stepping stones and corridors for species and ecosystem functions to persist. Further, they connect residents to nature and provide an experiential educational tool for science and conservation awareness (Elmqvist et al. 2015).
An ongoing challenge in urban land management is maintaining ecosystem functioning with minimal human interference (Lundholm 2015). Conversely, the task of managing fire-dependent forests is laden with complexities (van Wilgen et al. 2012). Negative public perception and safety concerns further impede fire management (Ryan et al. 2013). As such, fire is commonly excluded in the urban landscape (Koptur 2006), despite the increasing need to conserve fire-dependent forest fragments (Jones and Koptur 2017). Instead, suppressing fire results in losses of plant functional groups (Folke et al. 2004;Johnstone et al. 2016), such as the replacement of fire-adapted plants by fire-intolerant hardwoods. Turnover of the plant community is often accompanied by changes in the fuel structure and increased vegetation density (Mitchell et al. 2009). Furthermore, long-term fire suppression can result in alternative stable states that are likely not reversible. For example, firedependent pine rockland ecosystems in South Florida, USA, can transition into late successional hardwood hammocks in as little as 20 to 30 years without fire (Ross et al. 2009). In urban pine rockland preserves, this transition by native hardwoods is often accompanied by encroaching invasive species (Possley et al. 2014).
Urban fire restoration is a relevant and growing field of study. Less than 2% of pine rockland remains intact in South Florida, and nearly 20% remaining occurs in small urban fragments (USFWS 1999). Expanding the footprint of this globally endangered ecosystem is gaining interest in the conservation community, but restoration outcomes of introduced fire are not well studied in constructed preserves. Furthermore, evidence of initial vegetation responses to prescribed burning following prolonged exclusion is scarce (Kobziar et al. 2015). Short-term restoration assessments generally include surveys few to many months post-burn (Abrahamson and Abrahamson 1996;Abrahamson et al. 2021). Although delayed surveys are informative for long-term restoration goals, immediate responses such as seedling recruitment and surface litter removal may remain undetected with increasing time between surveys. For example, fast-growing native and invasive species may respond rapidly to increased nutrients and sunlight when fire has long been excluded (Abrahamson and Abrahamson 1996). Under such circumstances, seedlings may be outcompeted, yet these seed bank contributions are significant for restoration potential (Cohen 2004).
We investigated the initial plant community responses to prescribed fire in a university pine rockland preserve with no previous known fire history. Namely, we were interested in immediate responses and those occurring within less than 4 months. Our specific goal was to understand how vegetation structure and plant community composition initially respond to fire, including seedling recruitment and resprouting of native and invasive species.

Site description
The prescribed burn was conducted at Florida International University's Nature Preserve (25.7548°N, 80.3794°W) in Miami, FL. The 0.4-ha pine rockland (Pineland) is centered within a 3.7-ha preserve, designed in 1978. Cleared trails surround the Pineland, separating it from a wetland ecosystem and hardwood hammocka South Florida late successional forest. The initial and subsequent plantings include species known to occur in pine rockland, and a subset of this accumulated planting list was present in our pre-burn survey. Besides supplemental plantings and occasional removals of invasive plants, the preserve was managed only by hand for the 40 years before our first prescribed burn. Consequently, hardwood trees and invasive shrubs and grasses were abundant.
Pine rockland is maintained by a fire frequency of 3-10 years (Snyder et al. 1990). Forests managed with frequent fire regimes are sparsely dominated by a Florida slash pine (Pinus elliottii Engelm) canopy with few midcanopy palms and shrubs (Snyder et al. 2005), like saw palmetto (Serenoa repens (W. Bartram) Small) and Jamaican nettle tree (Trema micrantha (L.) Blume). The herbaceous understory is species rich and includes several endemic and rare plants, such as Small's milkpea (Galactia smallii H.J. Rogers ex Herndon; Gann et al. 2020;O'Brien 1998). Soils are shallow, and the oolitic limestone is often visible between herbaceous plants. The prolonged absence of fire alters species composition as well as vegetation and fuel structure. Specifically, herbaceous plant diversity wanes, vegetation becomes shrubbier, and heavier fuel loads accumulate with time since burn (O'Brien 1998).

Prescribed burn
The prescribed burn was conducted by the Florida Department of Agriculture and Consumer Services Division of Forestry on 15 March 2016. The approximate temperature before fire initiation was 82°F with 7 mph average wind speeds and 60% relative humidity. The whole Pineland within the preserve was burned. Backing fires, strip heading fires, heading fires, flanking fires, and dot fires were the burn techniques prescribed. The maximum acceptable crown scorch was between 10 and 40%. Although the fire consumed much leaf litter, patches of unburned litter were present throughout the site, indicating a low-intensity fire.

Field design
At each survey, we established ten 1-m 2 quadrats from randomly selected points inside the Pineland, excluding a 5-m buffer from the perimeter. We measured canopy closure, ground cover composition, and plant composition within the quadrats for four survey periods: days before the burn (T0), 1 week (T1), 4 weeks (T4), and 14 weeks (T14) after the burn. Our surveys occurred from March to early June before the onset of the peak rainy season. Weekly average temperatures and rainfall were 76°F and 0 cm at T0 and 82°F and 1.17 cm at T14.
We measured the changes in environmental variables that are responsive to fire and important for firedependent understory plant communities, including light availability and leaf litter, bare soil, and vegetation cover. To measure light availability, we used a spherical densitometer at 1 m above each plot, and we averaged the proportion of the sky hemisphere obscured by canopy vegetation from four cardinal directions (Jennings 1999). We then estimated ground cover composition as the proportion of quadrats (1-100%) covered by leaf litter, by shrubby and woody vegetation, by exposed soils or bare limestone, and by herbaceous vegetation. To analyze the changes in plant community structure and composition, we approximated species abundances as the sum of their presences or absences (incidences; Colwell et al. 2004) within all quadrats (n = 10) at each survey (n = 4). Incidences are a proxy for abundances when individuals are difficult to ascertain, such as in resprouting and vegetatively clonal individuals (Colwell et al. 2004).

Analyses
All analyses were performed in R (ver. 3.4.0). We analyzed the changes in the proportions of canopy closure before and post-recovery using a beta-regression (Cribari-Neto and Zeileis 2010) with weeks since burn as the categorical response variable. To test for compositional changes in ground cover, we ran a compositional analysis of variance (van den Boogaart and Tolosana-Delgado 2013). We first reclassified the vegetation cover for each sample as the summed proportions of herbaceous and shrubby plant cover; proportions of bare ground and leaf litter were maintained. Where ground cover classes were observed to be less than 1%, we imputed the rounded zeroes using a robust linear model method (Templ et al. 2011) before taking the isometric log ratios between the cover types.
We rarefied species richness for each survey, extrapolated it to double the sample size (Hsieh et al. 2016), and calculated the confidence intervals for comparison of species richness at different stages of recovery (Colwell et al. 2004;Chao et al. 2014). We tested if the native status (native or invasive) composition was different between surveys by performing a chi-squared test of differences on counts of invasive and native species present at each survey. We then tested the probability that frequencies of plants present at each survey came from the same distributional species pool. We calculated a Jaccard binary distance matrix (Oksanen 2001) from species incidences across quadrats at each survey (n = 40), and we performed an analysis of similarity (ANO-SIM) on the distance matrix with weeks since prescribed burn as a categorical explanatory variable (permutations = 9999). To visualize these compositional changes, we performed a correspondence analysis (CA; Legendre and Legendre 2012) of species incidences before the prescribed burn and post-recovery. We displayed the first two axes of the CA, preserving the X 2 distance among survey periods ( Fig. 3; Additional file 1: Table S1).

Results
Much of the upper canopy leaves and branches burned in the fire, reducing closure from 75% (± 0.03) to 36% (± 0.02). New leaves filled in slightly to 46% closure (± 0.03) by T14, but the overall available sunlight to the understory increased (Φ = 29.4, F 3,36 = 13.5, P < 0.001) after burning. Similarly, fire consumed leaf litter on the ground and increased bare soil, but vegetation filled in all available soil space by T14 (Fig. 1). Woody plant cover was initially reduced post-burn but recovered to higher levels by T14. These changes to proportional ground cover compositions of vegetation, bare soil, and leaf litter were significantly different across survey periods (Pillai's trace 1.16, F 3, 36 = 16.7, P < 0.001). However, T1 and T4 shared similar canopy cover and ground compositional measurements.
There was a significant effect of time since burn on observed species incidences (F 3,192 = 9.9, P < 0.001). Of the 49 unique species identified across all surveys, 37 were native species and 12 were invasive. From this species pool, 27 were present before the prescribed burn (T0), compared to 38 species after 14 weeks (T14). Additionally, five species were observed within 4 weeks postburn that were not detected before the burn or 14 weeks after burning. Despite losses and gains of species over time, confidence intervals of T0 and T14 completely overlapped after extrapolating species accumulation curves to asymptotic coverage (Fig. 2). Similarly, invasive species richness did not significantly change across surveys (X 3 2 = 3.19, P > 0.05). While ten invasive species occurred before the burn, seven were detected in T14.
Despite the similarity in extrapolated species richness, ANOSIM results indicate that plant composition was more different across surveys than within surveys (permutations = 9999, Global R = 0.22, P < 0.001). For instance, 23 species occurring in the three surveys after the burn were not present in pre-burn surveys, including two invasive shrub species, one native shrub, and 20 native herbs (S1). We see 79% of this variability in species incidences separated along the first two axes of the CA (total inertia = 0.82; Fig. 3). Survey period scores from the CA for T1 and T2 separated from T14 along the first axis (0.34 Eigenvalues) and from T0 on the second axis (0.30 Eigenvalues). Further, the composition of our final survey (T14) partitioned from the pre-burn composition (T0) on both CA axes (Fig. 3).

Discussion
The objective of the prescribed burn on the university campus preserve was to restore the plant composition of the Pineland, reminiscent of a natural pine rockland ecosystem. Specific goals were to encourage the establishment of native herbaceous plants ( Barrios et al. 2011;O'Brien 1998) by reducing biomass of fire-sensitive hardwoods (Mitchell et al. 2009), decreasing canopy light interception (Cannon et al. 2017), consuming leaf litter (Ormeño et al. 2009), and providing germination opportunities from the seed bank (Liu et al. 2005). Our shortterm pre-and post-burn assessment permitted early observations of competition and rapid replacement between seedling recruitment and surviving plants. These immediate responses can inform both natural restoration potential and sources contributing to less successful restoration outcomes on a longer time frame.
Despite no former fire management in the Pineland, the low-intensity prescribed burn was effective at removing some canopy and surface fuels with little damage to fire-resistant species. Nonetheless, structural vegetation changes appeared to be short-lived, such as described in similar fire-suppressed pine forests (Abrahamson and Abrahamson 1996). For example, leaf litter was consumed but soon replenished with charred leaves and branches from the canopy. Additionally, exposed soil cover increased immediately post-burn but was replaced by a high density of germinating and resprouting plants within 14 weeks, including many shrub seedlings.
Plant compositional changes were also variable over a short time frame. Mortality was high for some species after the burn, but several new species germinated. Perennial fast-growing herbs like the invasive lantana (Lantana strigocamara R.W. Sanders) replaced large shrubs like lancewood (Damburneya coriacea (Sw.) Trofimov & Rohwer) from the pre-burn community. However, juveniles of shrubs, like Jamaican nettle tree and spine fiddlewood (Citharexlyum spinosum L.), that were not abundant before the burn were frequent in the last survey. The earliest post-burn surveys also revealed important seed bank contributions. For instance, short-lived herbs, like goldenrods (Solidago L. spp.) and groundcherries (Physalis walteri Nutt.), germinated and reproduced within weeks post-burn.
We had few expectations for managing invasive species with fire, since no consensus exists on using fire to reduce invasive plants (Alba et al. 2015), largely due to species-specific differences. Indeed, invasive species richness was unchanged, but species did vary in their responses. Some invasive plants became more common (e.g., lantana and coatbuttons (Tridax procumbens L.)), others were not affected (e.g., Chinese ladder brake (Pteris vittata L.)), and some became less common, particularly invasive woody plants like shoebutton (Ardisia elliptica Thunb.) and Brazilian pepper (Schinus terebinthifolius Raddi). Invasive shrubs like Brazilian pepper can alter fire behavior and reduce the effectiveness of prescribed burning for restoration. However, when the invasive shrub's abundance is low, fires can cause high Fig. 1 Compositions of ground cover in a managed pine rockland preserve in South Florida measured prior to the first prescribed burn (T0), 1 week post-burn (T1), 4 weeks (T4), and 14 weeks later (T14). The position of points in each ternary diagram illustrates the percent of ground surface within a quadrat covered by bare ground (soil), dead plant material (litter), and live plants (vegetation). Color depth of points separately indicates the percent of vegetation cover consisting of shrubs or trees. Ellipses within diagrams are produced from means and variances of compositions within quadrats mortality and be an effective restoration tool (Stevens and Beckage 2009). Low pre-burn abundances of Brazilian pepper and other invasive shrubs in the Pineland may explain the higher mortality of these species that we encountered.
Nevertheless, periodic fires can increase resiliency and inhibit the initial establishment of exotic plants in uninvaded communities (Ross et al. 2002;Debuse and Lewis 2014). Such factors are critical for small urban preserves, where invasive plants are more prevalent in the surrounding landscape (Koptur 2006), as is the situation for the Pineland preserve. Future restoration efforts could consider manual reductions of invasive plants before and after fire prescriptions to increase fire responses to undesirable shrubs, like Brazilian pepper (Stevens and Beckage 2009), and invasive herbs that benefit from fire, like lantana (Sharma et al. 2005).
Outcomes of the first prescribed fire in the Pineland are consistent with vegetation responses to fire reintroductions in analogous fire-suppressed pine forests (Abrahamson and Abrahamson 1996), reiterating that a single fire may not be sufficient for structural restoration when fire has long been excluded. Similarly, a single introduced fire and tree trimming improved plant community composition two decades after converting an old field to scrub habitat in Florida; however, structural and compositional similarity to a natural scrub ecosystem was not fully achieved with one fire (Schmalzer and Foster 2019). Although efforts aimed at introducing fire in urban areas are minimal (Jones and Koptur 2017), our results indicate that prescribed fire in constructed ecosystems is a positive step forward in maintaining these complex ecosystems.
Maintaining fire-dependent ecosystems at the wildland-urban interface could release some pressures of habitat loss and fragmentation caused by urbanization (Diamond and Heinen 2016). Periodic fire can stimulate germination, growth, and flowering in some plants ( Barrios et al. 2011;Pyke 2017), which can have a positive effect on associated animals (Dell et al. 2019). When these ecosystems are distributed across the landscape at different stages of burn recovery, diversity increases spatially (Turner 1989) and temporally. The Pineland in our study is a constructed ecosystem, but our post-burn seedling recruitment suggests that these planted native species made large contributions to the seed bank. These outcomes are promising for conserving rare and endangered ecosystems like pine rockland in urban fragments.

Conclusion
The early stages of recovery that we observed in this pine rockland forest offers insight into restoration Fig. 2 Estimated species accumulation curve before and at multiple periods after a prescribed fire in a managed pine rockland on the University Campus Nature Preserve in South Florida, USA. Solid lines are plant species richness interpolated by species incidences in each sample quadrat, points indicate observed richness, and dashes are extrapolated species richness to complete sample coverage. Quadrats were surveyed 1 week before the prescribed burn (T0; circle; n = 10), 1 week post-burn (T1; triangle), 4 weeks after (T4; square), and at 14 weeks (T14; rhombus). Confidence intervals of interpolated and extrapolated species richness are shaded with corresponding colors of sample times as indicated in the legend outcomes in previously unburned fragments. Despite the high density of herbs and shrubs 14 weeks postburn, temporal diversity was high. Many novel herbs germinated post-burn, and some short-lived annuals made important contributions back to the seed bank. The burn effectively consumed much biomass and fuel loads, while cycling nutrients (Freidenreich et al. 2020), with minimal human interference. Such feasibility of prescribed burns in unburned urban patches are critical for ecosystems like the globally endangered pine rockland, in which nearly 20% of the remaining habitat occurs within urban fragments averaging less than 5 ha (USFWS 1999).
Abbreviations ANOSIM: Analysis of similarities; CA: Correspondence analysis; Pineland: Pine rockland forest in the Florida International University Nature Preserve; T0: Days before the prescribed burn; T1: One week after the prescribed burn; T4: Four weeks after the prescribed burn; T14: Fourteen weeks after the prescribed burn; X 2 : Chi-squared
Additional file 1: Table S1. Plant species present in the fire-dependent pine rockland preserve at Florida International University, Miami, FL before (T0), and one (T1), four, (T4) and fourteen weeks (T14) after the first prescribed fire. Correspondence analysis axes scores represent scores for species incidences at each survey period. Vegetation type and regional native status in South Florida were provided by Institute for Regional Conservation's online database (Gann et al. 2020).

Fig. 3
Correspondence analysis (CA) of relative species frequencies following a prescribed burn within a previously unburned pine rockland fragment in South Florida, USA. Observations occurred 1 week before the prescribed burn (T0), 1 week post-burn (T1), at 4 weeks (T4), and at 14 weeks (T14). Species were scaled to preserve Χ 2 distances between survey times. ANOSIM results indicate that the means of ranked dissimilarities between the groups are more different than within the groups (R = 0.22, P < 0.001). The total contribution to inertia was 0.82, and this CA graph displays 79% of the variability. Only species contributing larger than 90% to the inertia are labeled (see Supplementary Table), and they are coded by vegetation type and regional native status in South Florida (Gann et al. 2020). Some random noise was added to species scores to visualize (± 0.07 points) overlapping labels