Data from: Fire season and drought influence fire effects on invasive grasses: A meta-analysis
Data files
Mar 14, 2025 version files 322.42 KB
-
InvasiveBluestemMetaAnalysisData.csv
172.73 KB
-
NativeGrassMetaAnalysisData.csv
128.51 KB
-
README.md
21.19 KB
Abstract
Fire can shape plant communities when species respond differently to fire characteristics like season and intensity. If invasive plants are more vulnerable than native species to specific aspects of fire, managers could use prescribed fires to control non-native species. We conducted a meta-analysis of fire effects on six closely related Afro-Eurasian and Australian grasses (Bothriochloa bladhii, B. ischaemum, B. pertusa, Dichanthium annulatum, D. aristatum, and D. sericeum, collectively “invasive bluestems”) that have invaded grasslands worldwide. Using data from 31 studies (only 3 from their native range), we evaluated whether these grasses (275 effect sizes) responded differently than native grasses (184 effect sizes) to fire season, fuel load, and drought. Original data included 15 different response variables (e.g., biomass, survival) that were converted to standardized mean difference for analysis. Fires in summer, fall, and early winter had negative effects on invasive bluestems; no fire season had significant positive effects. Most data were for B. ischaemum, but the other bluestems may also be vulnerable to summer fire. Native grasses did not show significant negative responses in any month. Drought (Keetch-Byram Drought Index) in the month before fire increased the negative effects of fire on invasive bluestems but not native grasses. Drought after fire led to similar negative effects on both groups. Unexpectedly, fuel load (which influences fire intensity) did not significantly influence fire effects in any analysis. At the fuel loads examined (~600 – 10,000 kg/ha dried herbaceous biomass), fire intensity may have been too low to cause meristem mortality. Between-study heterogeneity was large in all analyses (I2>80%), suggesting that additional factors beyond those reported in the studies influence fire effects. These factors could include plant phenology, fire behavior, weather conditions during the fire, and soil characteristics.
Synthesis and applications:
Fires during summer and fall, especially during dry conditions, could harm invasive bluestems relative to native grasses, likely due to subtle differences in heat sensitivity, phenology, and drought resistance. Other invasive species may have similar vulnerabilities to specific fire seasons and rainfall conditions that allow the use of fire as a control method.
https://doi.org/10.5061/dryad.xpnvx0ks4
Description of the data and file structure
Files and variables
File: InvasiveBluestemMetaAnalysisData.csv
Description: File “InvasiveBluestemMetaAnalysisData.csv” contains the processed data used in the invasive bluestem meta-analyses, along with a few extra columns that are useful for describing the studies (noted as “not used in the meta-analysis”) or where used to create graphs. Missing values are indicated as NA. KBDI = Keetch-Byram Drought Index, which represents the amount of rainfall needed to return the soil to saturation and changes daily based on temperature and rainfall. Values range from 0 (no moisture deficit) to 800.
Variables
- PubID: Numerical identifier for the publication.
- Site.Name: General location of the study. Not used in the meta-analysis.
- Site.SubsiteName: Specific location of the study. Named subsites within one site were treated as separate sites because they often had different soils or management histories. Used as “site” in the meta-analysis.
- Species: Species of grass as listed by the authors of the study.
- InvasiveGroup: “BotIsch” for Bothriochloa ischaemum (the most common grass) and “all other IB” for all other invasive bluestems
- Species.Range: In this location, is the species “native” or “introduced”? Note that some of the invasive bluestem studies were conducted in the bluestems’ native ranges.
- Hemisphere: “north” or “south”. Not used in the meta-analysis.
- Location.Country: Country where study was conducted. Not used in the meta-analysis.
- Location.StateProvince: State/province where the study was conducted. Not used in the meta-analysis.
- FireMonthYear: The month and year in which the fire was conducted. Used to calculate monthly KBDI variables and time since fire.
- Fire.Date: The exact date of the fire. Not provided by all studies. Used to calculate weekly KBDI variables.
- AdjFireMonth: Numerical month. Because we found more sites in the northern hemisphere, fire months from the southern hemisphere were adjusted by adding 6 months, making them seasonally equivalent to northern hemisphere months (e.g., “July” is always summer).
- Study.Type: ”experiment” (plots randomly assigned to a treatment) or “observational” (e.g., comparing burned and unburned sites after a wildfire or prescribed fire; plots not randomly assigned to a treatment).
- Fire.Type: ”prescribed” fire, “wildfire” or “burn barrel” (metal box/barrel to contain a plot-level burn)
- Sample.Date: Exact date on which the vegetation data were collected. Not provided by all studies. Not used in the meta-analysis.
- SampleMonthYear: Month and year in which the vegetation data were collected. Used to calculate time since fire.
- KBDI2MonthBeforeFire: Average KBDI for the two months before the month of fire (e.g., for January and February before a March fire). Not used in the meta-analysis but used for graphing.
- KBDI2MonthBeforeFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- KBDIMonthBeforeFire: Average KBDI for the month before the month of fire (e.g., February for a March fire). Not used in the meta-analysis but used for graphing.
- KBDIMonthBeforeFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- KBDIMonthOfFire: Average KBDI in the month of the fire (e.g., the month of March for a March fire). Not used in the meta-analysis but used for graphing.
- KBDIMonthOfFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- KBDIMonthAfterFire: Average KBDI in the month after the month of the fire (e.g., April for a March fire). Not used in the meta-analysis but used for graphing.
- KBDIMonthAfterFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- KBDI2MonthAfterFire: Average KBDI for the two months after the month of the fire (e.g., April and May for a March fire). Not used in the meta-analysis but used for graphing.
- KBDI2MonthAfterFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- KBDI10WksPreFire: Average KBDI for the 10 weeks before the exact fire date. Not available for all studies. Not used in the meta-analysis but used for graphing.
- KBDI10WksPreFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- KBDI6WksPreFire: Average KBDI for the 6 weeks before the exact fire date. Not available for all studies. Not used in the meta-analysis but used for graphing.
- KBDI6WksPreFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- KBDI2WksPreFire: Average KBDI for the 2 weeks before the exact fire date. Not available for all studies. Not used in the meta-analysis but used for graphing.
- KBDI2WksPreFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- KBDI2WksPostFire: Average KBDI for the 2 weeks after the exact fire date. Not available for all studies. Not used in the meta-analysis but used for graphing.
- KBDI2WksPostFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- KBDI6WksPostFire: Average KBDI for the 6 weeks after the exact fire date. Not available for all studies. Not used in the meta-analysis but used for graphing.
- KBDI6WksPostFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- KBDI10WksPostFire: Average KBDI for the 10 weeks after the exact fire date. Not available for all studies. Not used in the meta-analysis but used for graphing.
- KBDI10WksPostFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- Sample.Grazed: ”grazed” (during study period), “ungrazed”, “both” (data averaged across grazed and ungrazed treatments)
- Sample.Seeded: Seeded (“y”) or unseeded (“n”) with native species. Studies seeded with invasive bluestem species were excluded. Only used in invasive bluestem analyses (i.e., studies seeded with native species were not used in the native grass analyses).
- Fuelkgha: Dried herbaceous biomass at the time of fire in kg/ha. Not all studies provided this information. Not used in the meta-analysis but used for graphing.
- Fuelkgha.cen: Centered and scaled biomass (mean is subtracted and the remainder is divided by the standard deviation).
- TSFyr: Time (in years) between the month of the fire and the month of vegetation data collection.
- Site.SoilDepthA: Depth of A horizon in cm. Not all studies provided this information. Not used in the meta-analysis but used for graphing.
- Soil.cen: Centered and scaled soil depth (mean is subtracted and the remainder is divided by the standard deviation).
- AbsoluteLatitude: The absolute value of the study site latitude. Latitudes were provided in the studies or were estimated based on the described study location (e.g., nearest town, center of the county, etc.). Not used in the meta-analysis but used for graphing.
- Latitude.cen: Centered and scaled latitude (mean is subtracted and the remainder is divided by the standard deviation).
- Data.VariableType: The variable type from the original study, which included basal area, biomass, change in cover, change in frequency, cover, crown area, dead crown density, density, frequency, number of plants, number of seed heads, number of tillers, survival, and stem count.
- EffectSize: Standardized mean difference, where negative values indicate a negative effect on the target grass.
- Variance: Variance of the effect size
- Publication: Author and year of the publication. See Data Sources for full citations.
File: NativeGrassMetaAnalysisData.csv
Description: File “NativeGrassMetaAnalysisData.csv” contains the processed data used in the native grass meta-analyses, along with a few extra columns that are useful for describing the studies (noted as “not used in the meta-analysis”) or where used to create graphs. Missing values are indicated as NA. KBDI = Keetch-Byram Drought Index, which represents the amount of rainfall needed to return the soil to saturation and changes daily based on temperature and rainfall. Values range from 0 (no moisture deficit) to 800.
Variables
- PubID: Numerical identifier for the publication.
- Site.Name: General location of the study. Not used in the meta-analysis.
- Site.SubsiteName: Specific location of the study. Named subsites within one site were treated as separate sites because they often had different soils or management histories. Used as “site” in the meta-analysis.
- Species: Species of grass as listed by the authors of the study. If listed as “native grasses” or “native perennial grasses”, data do not include invasive bluestems.
- Phototype: Photosynthesis type (“C3” or “C4”) as defined by Cerros-Tlatilpa et al. (2011) and Osborne et al. (2014).
- Hemisphere: “north” or “south”. Not used in the meta-analysis.
- Location.Country: Country where study was conducted. Not used in the meta-analysis.
- Location.StateProvince: State/province where the study was conducted. Not used in the meta-analysis.
- FireMonthYear: The month and year in which the fire was conducted. Used to calculate monthly KBDI variables and time since fire.
- Fire.Date: The exact date of the fire. Not provided by all studies. Used to calculate weekly KBDI variables.
- AdjFireMonth: Numerical month. Because we found more sites in the northern hemisphere, fire months from the southern hemisphere were adjusted by adding 6 months, making them seasonally equivalent to northern hemisphere months (e.g., “July” is always summer).
- Study.Type: ”experiment” (plots randomly assigned to a treatment) or “observational” (e.g., comparing burned and unburned sites after a wildfire or prescribed fire; plots not randomly assigned to a treatment).
- Fire.Type: ”prescribed” fire, “wildfire” or “burn barrel” (metal box/barrel to contain a plot-level burn)
- Sample.Date: Exact date on which the vegetation data were collected. Not provided by all studies. Not used in the meta-analysis.
- SampleMonthYear: Month and year in which the vegetation data were collected. Used to calculate time since fire.
- KBDI2MonthBeforeFire: Average KBDI for the two months before the month of fire (e.g., for January and February before a March fire). Not used in the meta-analysis but used for graphing.
- KBDI2MonthBeforeFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- KBDIMonthBeforeFire: Average KBDI for the month before the month of fire (e.g., February for a March fire). Not used in the meta-analysis but used for graphing.
- KBDIMonthBeforeFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- KBDIMonthOfFire: Average KBDI in the month of the fire (e.g., the month of March for a March fire). Not used in the meta-analysis but used for graphing.
- KBDIMonthOfFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- KBDIMonthAfterFire: Average KBDI in the month after the month of the fire (e.g., April for a March fire). Not used in the meta-analysis but used for graphing.
- KBDIMonthAfterFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- KBDI2MonthAfterFire: Average KBDI for the two months after the month of the fire (e.g., April and May for a March fire). Not used in the meta-analysis but used for graphing.
- KBDI2MonthAfterFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- KBDI10WksPreFire: Average KBDI for the 10 weeks before the exact fire date. Not available for all studies. Not used in the meta-analysis but used for graphing.
- KBDI10WksPreFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- KBDI6WksPreFire: Average KBDI for the 6 weeks before the exact fire date. Not available for all studies. Not used in the meta-analysis but used for graphing.
- KBDI6WksPreFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- KBDI2WksPreFire: Average KBDI for the 2 weeks before the exact fire date. Not available for all studies. Not used in the meta-analysis but used for graphing.
- KBDI2WksPreFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- KBDI2WksPostFire: Average KBDI for the 2 weeks after the exact fire date. Not available for all studies. Not used in the meta-analysis but used for graphing.
- KBDI2WksPostFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- KBDI6WksPostFire: Average KBDI for the 6 weeks after the exact fire date. Not available for all studies. Not used in the meta-analysis but used for graphing.
- KBDI6WksPostFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- KBDI10WksPostFire: Average KBDI for the 10 weeks after the exact fire date. Not available for all studies. Not used in the meta-analysis but used for graphing.
- KBDI10WksPostFire.cen: Same as above but centered on 400 (the midpoint of KBDI) and scaled from -1 to 1.
- Sample.Grazed: ”grazed” (during study period), “ungrazed”, “both” (data averaged across grazed and ungrazed treatments)
- Fuelkgha: Dried herbaceous biomass at the time of fire in kg/ha. Not all studies provided this information. Not used in the meta-analysis but used for graphing.
- Fuelkgha.cen: Centered and scaled biomass (mean is subtracted and the remainder is divided by the standard deviation).
- TSFyr: Time (in years) between the month of the fire and the month of vegetation data collection.
- Site.SoilDepthA: Depth of A horizon in cm. Not all studies provided this information. Not used in the meta-analysis but used for graphing.
- Soil.cen: Centered and scaled soil depth (mean is subtracted and the remainder is divided by the standard deviation).
- AbsoluteLatitude: the absolute value of the study site latitude. Latitudes were provided in the studies or were estimated based on the described study location (e.g., nearest town, center of the county, etc.). Not used in the meta-analysis but used for graphing.
- Latitude.cen: Centered and scaled latitude (mean is subtracted and the remainder is divided by the standard deviation).
- Data.VariableType: The variable type from the original study, which included biomass, change in frequency, cover, crown area, frequency, number of plants, number of seed heads, number of tillers, survival, and stem count.
- EffectSize: Standardized mean difference, where negative values indicate a negative effect on the target grass.
- Variance: Variance of the effect size
- Publication: Author and year of the publication. See Data Sources for full citations.
NA- Not applicable
Code/software
The comma-separated values (csv) file can be imported directly into R or opened in a text file editor. Note that opening the file in Excel can re-assign data in the “FireMonthYear” column in unexpected ways.
Access information
Data was derived from the following sources:
Andruk, C.M. (2014). Restoration of central Texas savanna and woodland: The effects of fire, deer, and invasive species on plant community trajectories [PhD]. University of Texas at Austin.
Behr, W.L. (2023). The effects of prescribed fire on grasslands of the southern Great Plains [PhD]. University of Texas at Austin.
Blair, J.M. 2024. BGPVC Plant species composition in the Belowground Plot Experiment at Konza Prairie. Environmental Data Initiative. http://dx.doi.org/10.6073/pasta/2199a82aedbd8e462ac9edafd915d0fc.
Clark, M.D. (2014). Evaluating the interactive effects of seasonal prescribed fire and grazing on invasive grass abundance and woody brush encroachment [M.S.]. Texas A&M University.
Davis, F.H. (2011). Effects of prescribed burning on King Ranch bluestem at vegetative regrowth and flowering stages [M.S.]. Texas State University.
Foster, J. L., Clayton, M. K., Lesak, M. M., McCuistion, K., & Teinert, T. (2025). Combining Management Techniques for Short-Term Reduction of Introduced Old World Bluestems in South Texas Rangelands. Rangeland Ecology & Management, 99, 50–57. https://doi.org/10.1016/j.rama.2024.12.001 (Data provided by authors before publication)
Giefer, H., Harmoney, K., Ramirez, M., Tajchman, A., Duncan, Z., Lemmon, J., & Olson, K. (2023). Effects of Late-Summer Prescribed Fire on Botanical Composition, Soil Cover, and Forage Production in Caucasian Bluestem-Infested Rangeland in the Kansas Smoky Hills: Year 4 of 5. Kansas Agricultural Experiment Station Research Reports, 9(1), 1.
Han, X., Dendy, S.P., Garrett, K.A., Fang, L., & Smith, M.D. (2008). Comparison of damage to native and exotic tallgrass prairie plants by natural enemies. Plant Ecology, 198(2), 197–210.
Havill, S., Schwinning, S., & Lyons, K.G. (2015). Fire effects on invasive and native warm-season grass species in a North American grassland at a time of extreme drought. Applied Vegetation Science, 18, 637–649.
Lesak, M.M. (2016). Management of invasive Old World bluestems to restore native grasslands [M.S.]. Texas A&M University.
Lewis, T. (2006). Management for conservation of plant diversity in native grasslands of the Moree Plains, NSW [PhD]. University of New England.
Lewis, T., Reid, N., Clarke, P.J., & Whalley, R.D.B. (2010). Resilience of a high-conservation-value, semi-arid grassland on fertile clay soils to burning, mowing and ploughing. Austral Ecology, 35(4), 464–481.
Minton, D. (2006). Fire, erosion, and sedimentation in the Asan-Piti watershed and War in the Pacific NHP, Guam (Technical Report 150). Pacific Cooperative Studies Unit, University of Hawaii at Manoa.
Orr, D. (2010). Managing the grass-legume balance in Stylosanthes scabra cv. Seca pastures in central Queensland. Tropical Grasslands, 44.
Ramirez, M., Tajchman, A., Duncan, Z., Lemmon, J., & Olson, K. (2021). Effects of Late Summer Prescribed Fire on Botanical Composition, Soil Cover, and Forage Production in Caucasian Bluestem-Infested Rangeland in the Kansas Smoky Hills: Year 2 of 4. Kansas Agricultural Experiment Station Research Reports, 7(1).
Ramirez, M., Tajchman, A., Duncan, Z., Lemmon, J., & Olson, K. (2022). Effects of Late Summer Prescribed Fire on Botanical Composition, Soil Cover, and Forage Production in Caucasian Bluestem-Infested Rangeland in the Kansas Smoky Hills: Year 3 of 4. Kansas Agricultural Experiment Station Research Reports, 8(1).
Reed, H.E., Seastedt, T.R., & Blair, J.M. (2005). Ecological consequences of C4 grass invasion of A C4 grassland: A dilemma for management. Ecological Applications, 15(5), 1560–1569.
Reemts, C.M., McCaw, W.M., Greene, T.A., & Simmons, M.T. (2019). Short-term control of an invasive C4 grass with late-summer fire. Rangeland Ecology and Management, 72(1), 182–188.
Reemts, C.M., Picinich, C., & Greene, T.A. (2021). Late-summer fire provides long-term control of the invasive Old World Bluestem (Bothriochloa ischaemum). Southeastern Naturalist, 20(4), 589–601.
Ruckman, E.M., Schwinning, S., & Lyons, K.G. (2012). Effects of Phenology at Burn Time on Post-Fire Recovery in an Invasive C4 Grass. Restoration Ecology, 20(6), 756–763.
Silcock, R.G., Jones, P., Hall, T.J., & Waters, D.K. (2005). Enhancing pasture stability and profitability for producers in Poplar Box and Silver-leaved Ironbark woodlands. (NAP3.208). Meat & Livestock Australia Limited.
Simmons, M.T., Windhager, S., Power, P., Lott, J., Lyons, R.K., & Schwope, C. (2007). Selective and non-selective control of invasive plants: The short-term effects of growing-season prescribed fire, herbicide, and mowing in two Texas prairies. Restoration Ecology, 15(4), 662–669.
St. Clair, B. (2012). Growing season versus dormant season fire as a control method for an invasive species, King Ranch bluestem (Bothriochloa ischaemum) [M.S.]. Texas A&M University-Commerce.
Toomey, A.E. (2015). Effects of seasonal prescribed burning on Kleberg bluestem (Dichanthium annulatum) in south Texas [M.S.]. Texas A&M University - Kingsville.
Twidwell, D., Rogers, W., McMahon, E., Thomas, B., Kreuter, U., & Blankenship, T. (2012). Prescribed extreme fire effects on richness and invasion in coastal prairie. Invasive Plant Science and Management, 5, 330–340.
Whiting, C. (2022). Ecological relationships between invasive grasses, native grasses, and wildfire [PhD]. University of Texas at Austin.
Windhager, S.W. (1999). An assessment of the use of seeding, mowing, and burning in the restoration of an oldfield to tallgrass prairie in Lewisville, Texas [Ph.D.]. University of North Texas.
We searched Web of Science, Agricola, Proquest, and GoogleScholar on July 15, 2023 (Table S1). A separate search was conducted in each database for each of the six focal species (Bothriochloa bladhii, Bothriochloa ischaemum, Bothriochloa pertusa, Dichanthium annulatum, Dichanthium aristatum, and Dichanthium sericeum). The search term included scientific and common names as well as “fire OR burn*” (Table S2); search terms were the same for each database. After de-duplication, these initial searches identified 272 unique documents (including theses, dissertations, and unpublished reports).
Table S1: Details of databases searched
|
Database |
Details |
|
Web of Science |
Searched ‘topic’. Web of Science Core Collection included Science Citation Index Expanded (SCI-EXPANDED—1900-present); Social Sciences Citation Index (SSCI)—1900-present; Arts & Humanities Citation Index (AHCI)—1975-present; Conference Proceedings Citation Index—Science (CPCI-S)—1900-present; Conference Proceedings Citation Index—Social Science & Humanities (CPCI-SSH)—1990-present; Book Citation Index—Science (BKCI-S)—2005-present; Book Citation Index—Social Sciences & Humanities (BKCI-SSH)—2005-present; Emerging Sources Citation Index (SCI)—2005-present; Current Chemical Reactions (CCR-EXPANDED)—1985-present; Index Chemicus (IC)—1993-present. |
|
Agricola |
1967 to present; searched ‘TX All Text Fields’ |
|
ProQuest Dissertations and Theses Global |
1974 to present; searched ‘anywhere but full text-NOFT’ |
|
GoogleScholar |
No date restrictions; used search strings with only scientific names (otherwise the search term was too long). Stopping point was chosen when there were 10 non-relevant search results on a page. |
Table S2: Exact search terms used when searching the databases listed in Table 1.
|
(“Bothriochloa bladhii” OR “Bothriochloa caucasica” OR “Andropogon bladhii” OR “Bothriochloa intermedia” OR “Andropogon intermedius” OR “Andropogon caucasicus” OR “Caucasian bluestem” OR “Australian beard grass” OR “forest-bluegrass” OR “plains bluestem” OR “purple plume grass”) AND (fire OR burn*) |
|
(“Bothriochloa ischaemum” OR “Amphilophis ischaemum” OR “Andropogon ischaemum” OR “Dichanthium ischaemum” OR “King Ranch bluestem” OR “KR bluestem” OR “yellow bluestem” OR “plains bluestem” OR “bearded finger grass” OR “dogstooth grass” OR “Turkestan bluestem”) AND (fire OR burn*) |
|
(“Bothriochloa pertusa” OR “Andropogon pertusus” OR “Holcus pertusus” OR “Indian-bluegrass” OR “pitted beardgrass” OR “hurricane grass” OR “Indian couch grass” OR “Seymour grass” OR “Barbados sourgrass” OR “Antigua hay” OR “sweet pitted grass” OR “silver grass”) AND (fire OR burn*) |
|
(“Dichanthium annulatum” OR “Andropogon annulatus” OR “Andropogon nodosus” OR “Kleberg bluestem” OR “marvel grass” OR “Diaz bluestem” OR “Hindi grass” OR “ringed dichanthium” OR “sheda grass” OR “medio bluestem” OR “jargu grass” OR “Delhi grass” OR “vuda bluegrass” OR “two-flowered golden-beard” OR “Santa Barbara grass”) AND (fire OR burn*) |
|
(“Dichanthium aristatum” OR “Andropogon aristatus” OR “Angleton bluestem” OR “Angelton bluestem” OR “wildergrass”) AND (fire OR burn*) |
|
(“Dichanthium sericeum” OR “Andropogon sericeus” OR “Silky bluestem” OR “Queensland blue grass” OR “silky bluegrass” OR “slender bluegrass” OR “tassel bluegrass”) AND (fire OR burn*) |
Search results were stored in Rayyan (Ouzzani et al., 2016) and were screened manually by one of the authors. Papers were retained for further screening if they included prescribed fire or wildfire, at least one of the focal grass species (in the native or introduced range), and a measure of fire effects on the focal grass. Title and abstract screening eliminated 168 records.
Remaining records underwent full-text screening with additional criteria: studies included a comparison between burned and unburned treatments; no additional treatments (e.g., fertilizer, mowing) were applied to the measured populations; studies provided information needed for a quantitative meta-analysis (response means, sample size, and measures of variation); and studies provided the month or exact dates of the fires. Where the necessary data were not available in the papers, we attempted to contact the authors. Full-text screening eliminated 76 records.
We then conducted a forward/backward citation search based on the papers included after full-text screening, as well as articles that were themselves eliminated because data were combined across treatments (e.g., averaged across fertilizer levels, n = 20), no variance was reported (n = 8), or fire history was not available (n = 2). These searches were conducted in August 2023 and added three additional papers for a final total of 31 papers that met all search criteria.
We extracted data from tables, figures (using WebPlotDigitizer 4.6, Rohatgi 2022), published supplemental datasets, or data provided by the authors. We extracted means for burn and control (unburned) treatments along with sample sizes and standard deviation or standard error. Original data types included basal area, biomass, change in cover, change in frequency, cover, crown area, dead crown density, density, frequency, number of plants, number of seed heads, number of tillers, survival, and stem count. We extracted data for the focal grass species as well as any native grass species presented in the same papers. When multiple papers were published about the same study, we used the data from the most recent publication, but included additional data from earlier papers if they presented different information (e.g., cover vs biomass; more detailed treatment groupings), taking care to not duplicate data.
We compiled additional information to serve as moderators (defined in Table S3) including site, latitude, species name, species range of bluestems (native vs introduced), photosynthesis type of native grasses (C3 and C4 species according to Cerros-Tlatilpa et al. [2011] and Osborne et al. [2014]), seeding with native species (excluding focal bluestems), current grazing, study type (experimental, observational), fire type (prescribed fire, wildfire, burn box/burn barrel), time since fire, date/month of fire, and response type (e.g., frequency, cover). Because we found more sites in the northern hemisphere, fire months from the southern hemisphere were adjusted by adding 6 months, making them seasonally equivalent to northern hemisphere months (e.g., “July” is always summer). We also recorded fuel load (dried herbaceous biomass) and soil depth when available.
To examine the influence of drought on fire effects, we calculated the Keetch-Byram Drought Index (KBDI, Keetch, and Byram 1968; Alexander 1990) before and after each fire. This drought index represents the amount of rainfall needed to return the soil to saturation and changes daily based on temperature and rainfall. Values range from 0 (no moisture deficit) to 800. We used this drought index instead of rainfall because the index takes mean annual rainfall into account, allowing comparisons among regions with different climates. To calculate the index, we downloaded temperature and precipitation data from the closest weather station(s) to each study site from National Oceanic and Atmospheric Administration’s National Centers for Environmental Information (for the United States and Guam, National Centers for Environmental Information, 2023) or the Australian Bureau of Meteorology (for Australia, Australian Bureau of Meteorology 2023). For studies where the exact fire date was reported, we calculated the average drought index for 2, 6, and 10 weeks before and after the fire (as 6 separate variables). For studies where only the month of fire was reported, we calculated the average drought index for the month in which the fire happened, as well as for 1 and 2 months before and after the fire (creating 5 drought variables).
Continuous moderators (fuel load, soil depth, and the absolute value of latitude) were centered on their mean and scaled by their standard deviation across the whole dataset. The drought variables were centered on 400 (the mid-point of the Keetch-Byram Drought Index) and scaled from -1 to 1, keeping all drought variables on the same scale.
The adjusted standardized mean difference (SMD, Hedges’ g) was calculated for each pair of burned and unburned data as the effect size for meta-analysis (Hedges & Olkin, 1985). Negative SMD values indicate a negative effect of fire (e.g., lower cover in burned vs unburned plots). SMD was calculated using the “esc” package as SMD = (mean1 – mean2)/√(((N1-1)*sd12 + (N2-1)*sd12)/(N1+N2-2)), where N = sample size and sd = standard deviation (Ludecke, 2019). In a few cases, study data were proportions or counts so we first calculated odds ratio or binary proportion effect sizes and then converted those into SMD using the “esc” package.