Data from: Management inputs, site conditions, and fire history shape outcomes of invasive plant control and native recovery
Data files
Dec 06, 2025 version files 62.86 KB
-
EAP24-0881_data.csv
55.01 KB
-
README.md
7.85 KB
Abstract
Millions of dollars and countless hours are spent on invasive plant management, and the field of invasion ecology has gained increasing attention in recent decades. Yet, despite these efforts to control and understand plant invasions, successful management is often elusive. Budgetary constraints are a common factor limiting invasive plant management programs, and therefore, optimizing control strategies are essential. However, such optimization requires data on management inputs and outcomes, and these data are often missing, lacking, or underutilized. To address this knowledge gap and identify predictors of successful weed control in natural areas, we examined nearly 20 years of invasive plant treatment data in the world’s largest urban national park – Santa Monica Mountains National Recreation Area of southern California. We resurveyed 279 sites that had undergone control in the past two decades, collecting data on the abundance of native and invasive plant species to evaluate long-term management outcomes. We used multiple statistical approaches to identify management inputs and site characteristics that are predictors of eradication, invasive plant cover, and native species recovery. We found that the greater the initial size or percent cover of an infestation, the lower the probability of eradication. We also found that weed infestations on steeper slopes and in areas that have burned more frequently are less likely to be eradicated. Promisingly, our results also showed that greater reductions in invasives generally benefited native plant communities, though not in all cases. These analyses also highlighted that persistence is key; more frequent treatments (both chemical and nonchemical) and greater investment of labor resulted in larger reductions in invasive plants. Our results highlight how site characteristics and limited resources can complicate invasive plant management, while demonstrating the value of analyzing treatment and monitoring data to identify effective control strategies and guide adaptive management decisions.
Dataset DOI: 10.5061/dryad.fj6q5747x
Description of the data and file structure
This dataset contains site-level monitoring data for multiple invasive plant species managed within the Santa Monica Mountains National Recreation Area (SMMNRA). Each row represents a single site or infestation polygon with associated measurements of infestation size, percent cover, target species identity, slope, elevation, fire history, and treatment outcomes.
The dataset supports analyses of:
- invasive species abundance and change through time
- site conditions associated with invasion severity
- treatment effectiveness
- environmental drivers (topography, fire regime, vegetation type)
Data were compiled from field observations, GIS analyses, and land-management records.
Files and variables
File: EAP24-0881_data.csv
Description: Primary dataset containing 44 columns and 279 rows (one per site).
Variables:
| Variable | Description |
| site_ID | Unique identifier for each site or infestation polygon. |
| peak_initial_cover | Maximum percent cover of the target invasive species at the initial assessment. |
| site_ha | Total site area (hectares). |
| inf_size_ha | Initial infestation area of the target invasive species (hectares). |
| target | Four-letter invasive species code for the focal species being monitored or treated. |
| target_cover_23 | Percent cover of the target invasive species in 2023. |
| eradicated | Whether the target species was eradicated from the site (Yes/No). |
| inf_size_23 | Infestation area of the target invasive species in 2023 (hectares). |
| cover_change | Change in percent cover from baseline to 2023 (percentage points). |
| change_infst_ha | Change in infestation area from baseline to 2023 (hectares). |
| trt_start | Year initial treatment began. |
| trt_end | Year most recent treatment occurred. |
| yrs_trtd | Total number of years during which treatment occurred. |
| yrs_since_trtd | Number of years since the last treatment. |
| treatment_type | Treatment category: chemical, nonchemical, or both. |
| foliar_spot | Number of times foliar spot-spray treatment was applied. |
| foliar_broadcast | Number of times foliar broadcast-spray treatment was applied. |
| mow | Number of times mowing treatment was conducted. |
| brushcut | Number of times brush-cutting or brush-removal treatment occurred. |
| manual | Number of times manual hand-pulling or hand-removal treatments were conducted. |
| cut_apply | Number of times cut-stump herbicide treatment was applied. |
| basal_bark | Number of times basal-bark herbicide treatment was applied. |
| times_treated | Total number of treatment actions at the site (sum of all methods). |
| labor_hours | Total labor hours spent treating the site. |
| hours_per_ha | Labor hours standardized by site area (hours per hectare). |
| survey_date | Date the site was surveyed in 2023. |
| native_cover | Percent cover of native plant species in 2023. |
| prop_native | Proportion of total vegetation cover composed of native species (0–1). |
| nonnative_cover_23 | Percent cover of nonnative species in 2023. |
| native_rich | Number of native plant species observed (species richness). |
| natives_m2 | Native plant density standardized per square meter. |
| nonnative_rich | Number of nonnative plant species observed (species richness). |
| aspect_azm | Aspect in azimuth degrees (0–360°). |
| min_slope | Minimum slope at the site (degrees). |
| max_slope | Maximum slope at the site (degrees). |
| mean_slope | Mean slope across the site (degrees). |
| major_veg | Dominant vegetation type at the site (e.g., chaparral, grassland). |
| fire_freq | Number of recorded fires at the site. |
| fire_return_interval | Average interval between fires (years). |
| last_fire | Year of the most recent fire. |
| yrs_postfire | Number of years since the most recent fire. |
| elev | Elevation (meters above sea level). |
| restored | Whether ecological restoration activities were conducted (Yes/No). |
| fuel_mod | Whether fuel-modification work occurred (Yes/No). |
Target species (column 'target'):
Code/software
All statistical analyses and data processing were conducted using R version 4.4.1 (R Core Team 2024), an open-source programming language and environment for statistical computing. Analyses relied on base R functions as well as several freely available R packages. No proprietary software was used. All workflows can be executed on Windows, macOS, or Linux systems capable of running R 4.4.1 or later.
The following packages were used in the analyses:
- tidyverse (v2.0.0) for data manipulation, cleaning, and visualization
- lme4 (v1.1–35.1) for fitting linear and generalized linear mixed-effect models
- car (v3.1–2) for ANOVA and regression diagnostics
- MASS (v7.3–60.0) for additional statistical tools
- randomForest (v4.7–1.1) for random forest models
- rfPermute (v2.5.2) for permutation-based importance testing in random forests
- ggplot2 (v3.4.4) for data visualization
- stats (base R) for parametric and nonparametric hypothesis testing
- vegan (v2.6–4) for diversity calculations (when evaluating richness metrics)
All code used to run the analyses can be executed using these packages with no additional dependencies.
