Post-Fire Forest Pulse Recovery: Superiority of Generalized Additive Models (GAM) in Long-Term Landsat Time-Series Analysis

Wildfires are increasing globally and pose major challenges for assessing post-fire vegetation recovery and ecosystem resilience. We analyzed long-term Landsat time series in two contrasting fire-prone ecosystems in the United States and Australia. Vegetation area was extracted using the Enhanced Vegetation Index (EVI) with Otsu thresholding. Recovery to pre-fire baseline levels was modeled using linear, logistic, locally estimated scatterplot smoothing (LOESS), and generalized additive models (GAM), and their performance was compared using multiple metrics. The results indicated rapid recovery of Australian forests to baseline levels, whereas this was not the case for forests in the United States. Among climatic factors, temperature was the dominant parameter in Australia (Spearman ρ = 0.513, <i>p</i> < 10^-8), while no climatic variable significantly influenced recovery in California. Methodologically, GAM consistently performed best in both regions due to its success in capturing multiphase and heterogeneous recovery patterns, yielding the lowest values of AIC (United States: 142.89; Australia: 46.70) and RMSE_cv (United States: 112.86; Australia: 2.26). Linear and logistic models failed to capture complex recovery dynamics, whereas LOESS was highly sensitive to noise and unstable for long-term prediction. These findings indicate that post-fire recovery is inherently nonlinear and ecosystem-specific and that simple models are insufficient for accurate estimation, with GAM emerging as an appropriate method for assessing vegetation recovery using remote sensing data. This study provides a transferable approach using remote sensing and GAM to monitor forest resilience under accelerating global fire regimes.