Mapping post-fire forest regeneration and vegetation recovery using a combination of very high spatial resolution and hyperspectral satellite imagery

Careful evaluation of forest regeneration and vegetation recovery after a fire event provides vital information useful in land management. The use of remotely sensed data is considered to be especially suitable for monitoring ecosystem dynamics after fire. The aim of this work was to map post-fire forest regeneration and vegetation recovery on the Mediterranean island of Thasos by using a combination of very high spatial (VHS) resolution (QuickBird) and hyperspectral (EO-1 Hyperion) imagery and by employing object-based image analysis. More specifically, the work focused on (1) the separation and mapping of three major post-fire classes (forest regeneration, other vegetation recovery, unburned vegetation) existing within the fire perimeter, and (2) the differentiation and mapping of the two main forest regeneration classes, namely, Pinus brutia regeneration, and Pinus nigra regeneration. The data used in this study consisted of satellite images and field observations of homogeneous regenerated and revegetated areas. The methodology followed two main steps: a three-level image segmentation, and, a classification of the segmented images. The process resulted in the separation of classes related to the aforementioned objectives. The overall accuracy assessment revealed very promising results (approximately 83.7% overall accuracy, with a Kappa Index of Agreement of 0.79). The achieved accuracy was 8% higher when compared to the results reported in a previous work in which only the EO-1 Hyperion image was employed in order to map the same classes. Some classification confusions involving the classes of P. brutia regeneration and P. nigra regeneration were observed. This could be attributed to the absence of large and dense homogeneous areas of regenerated pine trees in the study area.     

Assessing Forest Fire and Vegetation Recovery in the Black Hills,South Dakota

After 110 years of sustained fire suppression, the 2000 Jasper fire consumed about 33,785 ha (83,500 acres), or 12% of the Black Hills National Forest. We mapped the severity of the Jasper fire using Landsat imagery and then investigated post-fire vegetation regeneration conditions using field data, Quickbird imagery, and regression modeling. We found that fire scar and severity could be delineated and mapped accurately based on remotely sensed and field-acquired data. Results also revealed that vegetative recovery relative to burn severity, topography, and soil factors could be modeled effectively using local geographically weighted regression (GWR). Further regeneration assessment revealed that severely or heavily burned areas were more rapidly re-vegetated with grasses, forbs, and woody shrubs in the short term. The field survey showed that prescribed burns retard crown fires and that after eight years ponderosa pine seedlings have not re-established.     

Semi-automated mapping of burned areas in semi-arid ecosystems using MODIS time-series imagery

Understanding spatial and temporal patterns of burned areas at regional scales, provides a long-term perspective of fire processes and its effects on ecosystems and vegetation recovery patterns, and it is a key factor to design prevention and post-fire restoration plans and strategies. Remote sensing has become the most widely used tool to detect fire affected areas over large tracts of land (e.g., ecosystem, regional and global levels). Standard satellite burned area and active fire products derived from the 500-m Moderate Resolution Imaging Spectroradiometer (MODIS) and the Satellite Pour l’Observation de la Terre (SPOT) are available to this end. However, prior research caution on the use of these global-scale products for regional and sub-regional applications. Consequently, we propose a novel semi-automated algorithm for identification and mapping of burned areas at regional scale. The semi-arid Monte shrublands, a biome covering 240,000 km² in the western part of Argentina, and exposed to seasonal bushfires was selected as the test area. The algorithm uses a set of the normalized burned ratio index products derived from MODIS time series; using a two-phased cycle, it firstly detects potentially burned pixels while keeping a low commission error (false detection of burned areas), and subsequently labels them as seed patches. Region growing image segmentation algorithms are applied to the seed patches in the second-phase, to define the perimeter of fire affected areas while decreasing omission errors (missing real burned areas). Independently-derived Landsat ETM+ burned-area reference data was used for validation purposes. Additionally, the performance of the adaptive algorithm was assessed against standard global fire products derived from MODIS Aqua and Terra satellites, total burned area (MCD45A1), the active fire algorithm (MOD14); and the L3JRC SPOT VEGETATION 1 km GLOBCARBON products. The correlation between the size of burned areas detected by the global fire products and independently-derived Landsat reference data ranged from R² = 0.01–0.28, while our algorithm performed showed a stronger correlation coefficient (R² = 0.96). Our findings confirm prior research calling for caution when using the global fire products locally or regionally.     

Forest structural diversity characterization in Mediterranean landscapes affected by fires using Airborne Laser Scanning data

ABSTRACT

Forest fires can change forest structure and composition, and low-density Airborne Laser Scanning (ALS) can be a valuable tool for evaluating post-fire vegetation response. The aim of this study is to analyze the structural diversity differences in Mediterranean Pinus halepensis Mill. forests affected by wildfires on different dates from 1986 to 2009. Several types of ALS metrics, such as the Light Detection and Ranging (LiDAR) Height Diversity Index (LHDI), the LiDAR Height Evenness Index (LHEI), and vertical and horizontal continuity of vegetation, as well as topographic metrics, were obtained in raster format from low point density data. In order to map burned and unburned areas, differentiate fire occurrence dates, and distinguish between old and more recent fires, a sample of pixels was previously selected to assess the existence of differences in forest structure using the Kruskal–Wallis test. Then, k-nearest neighbors algorithm (k-NN), support vector machine (SVM) and random forest (RF) classifiers were compared to select the most accurate technique. The results showed that, in more recent fires, around 70% of the laser returns came from grass and shrub layers, yielding low LHDI and LHEI values (0.37–0.65 and 0.28–0.46, respectively). In contrast, the areas burned more than 20 years ago had higher LHDI and LHEI values due to the growth of the shrub and tree strata. The classification of burned and unburned areas yielded an overall accuracy of 89.64% using the RF method. SVM was the best classifier for identifying the structural differences between fires occurring on different dates, with an overall accuracy of 68.79%. Furthermore, SVM yielded an overall accuracy of 75.49% for the classification between old and more recent fires.     

Prototyping an artificial neural network for burned area mapping on a regional scale in Mediterranean areas using MODIS images

Each year thousands of ha of forest land are affected by forest fires in Southern European countries such as Spain. Burned area maps are a valuable instrument for designing prevention and recovery policies. Remote sensing has increasingly become the most widely used tool for this purpose on regional and global scales, where a large variety of techniques and data has been applied. This paper proposes a semiautomatic method for burned area mapping on a regional scale in Mediterranean areas (the Iberian Peninsula has been used as a study case). A Multi-layer Perceptron Network (MLPN) has been designed and applied to MODIS/Terra Surface Reflectance Daily L2G Global 500m SIN Grid multitemporal composite monthly images. The compositing criterion was based on maximum surface temperature. The research covered a six year period (2001–2006) from June to September, when most of the forest fires occur. The resulting burned area maps have been validated using official fire perimeters and compared with MODIS Collection 5 Burned Area Product (MCD45A1). The MLPN shown as an effective method, with a commission error of 29.1%, in the classification of the burned areas, while the omission error was of 14.9%. The results were compared with the MCD45A1 product, which had a slightly higher commission error (30.2%) and a considerably higher omission error (26.2%), indicating a high underestimation of the burned area.     

测绘地理信息智能化监测在森林防灭火工作中的应用

针对目前森林防灭火事前、事中和事后的监测需求,本文综合利用测绘地理信息技术,构建了测绘地理信息服务森林防灭火的综合性智能化监测技术体系。通过获取多源地理信息数据,构建森林防灭火数据库,面向日常防灭火能力评估、灾前常态化监测、灾中动态监测、灾后评估及灾后长期植被恢复监测,提供了一系列的智能化监测服务。通过在小珠山森林火灾中的实践检验了该体系的实用性,结果显示该监测体系可以为森林防灭火工作提供可靠、精准的智能化地理信息解决方案。     

Combining AVHRR imagery with CORINE Land Cover data to observe forest fires and to assess their consequences

This paper describes an operational application of AVHRR satellite imagery in combination with the satellite-based land cover database CORINE Land Cover (CLC) for the comprehensive observation and follow-up of 10000 fire outbreaks and of their consequences in Greece during summer 2000. In the first stage, we acquired and processed satellite images on a daily basis with the aim of smoke-plume tracking and fire-core detection at the national level. Imagery was acquired eight times per day and derived from all AVHRR spectral channels. In the second stage, we assessed the consequences of fire on vegetation by producing a burnt-area map on the basis of multi-annual normalised vegetation indices using AVHRR data in combination with CLC. In the third stage we used again CLC to assess the land cover of burnt areas in the entire country. The results compared successfully to available inventories for that year. Burnt area was estimated with an accuracy ranging from 88% to 95%, depending on the predominant land cover type. These results, along with the low cost and high temporal resolution of AVHRR satellite imagery, suggest that the combination of low-resolution satellite data with harmonised CLC data can be applied operationally for forest fire and post-fire assessments at the national and at a pan-European level.     

中国林火卫星遥感预警监测技术研究进展

中国目前已形成了地面巡护、近地面监测、航空巡护和卫星监测等4级立体林火监测层次,但森林火灾仍是造成中国森林资源损失、森林生态环境安全和人身伤害的主要林业灾害。为对林火预警监测技术研究提供技术借鉴参考,本文从可燃物参数估测、烟区识别、着火点检测、森林大火燃烧动态监测、森林火烧迹地制图、森林火灾受害程度评价、森林燃烧生物量估算和火后植被恢复监测等8个方面,对中国近二十多年来开展的林火卫星遥感预警监测应用技术的研究进展、存在的技术问题和发展趋势进行了分析,并对构建服务于生态文明建设的天—空—地一体化的林火预警监测技术体系建设进行了展望。     

Variation in tree mortality and regeneration affect forest carbon recovery following fuel treatments and wildfire in the Lake Tahoe Basin, California, USA

Background

Forest fuel treatments have been proposed as tools to stabilize carbon stocks in fire-prone forests in the Western U.S.A. Although fuel treatments such as thinning and burning are known to immediately reduce forest carbon stocks, there are suggestions that these losses may be paid back over the long-term if treatments sufficiently reduce future wildfire severity, or prevent deforestation. Although fire severity and post-fire tree regeneration have been indicated as important influences on long-term carbon dynamics, it remains unclear how natural variability in these processes might affect the ability of fuel treatments to protect forest carbon resources. We surveyed a wildfire where fuel treatments were put in place before fire and estimated the short-term impact of treatment and wildfire on aboveground carbon stocks at our study site. We then used a common vegetation growth simulator in conjunction with sensitivity analysis techniques to assess how predicted timescales of carbon recovery after fire are sensitive to variation in rates of fire-related tree mortality, and post-fire tree regeneration.

Results

We found that fuel reduction treatments were successful at ameliorating fire severity at our study site by removing an estimated 36% of aboveground biomass. Treated and untreated stands stored similar amounts of carbon three years after wildfire, but differences in fire severity were such that untreated stands maintained only 7% of aboveground carbon as live trees, versus 51% in treated stands. Over the long-term, our simulations suggest that treated stands in our study area will recover baseline carbon storage 10?C35?years more quickly than untreated stands. Our sensitivity analysis found that rates of fire-related tree mortality strongly influence estimates of post-fire carbon recovery. Rates of regeneration were less influential on recovery timing, except when fire severity was high.

Conclusions

Our ability to predict the response of forest carbon resources to anthropogenic and natural disturbances requires models that incorporate uncertainty in processes important to long-term forest carbon dynamics. To the extent that fuel treatments are able to ameliorate tree mortality rates or prevent deforestation resulting from wildfire, our results suggest that treatments may be a viable strategy to stabilize existing forest carbon stocks.     

Multi‐temporal Mapping of Burned Forest over Canada Using Satellite‐based Change Metrics

Abstract

A procedure for continental‐scale mapping of burned boreal forest at 10‐day intervals was developed for application to coarse resolution satellite imagery. The basis of the technique is a multiple logistic regression model parameterized using 1998 SPOT‐4 VEGETATION clear‐sky composites and training sites selected across Canada. Predictor features consisted of multi‐temporal change metrics based on reflectance and two vegetation indices, which were normalized to the trajectory of background vegetation to account for phenological variation. Spatial‐contextual tests applied to the logistic model output were developed to remove noise and increase the sensitivity of detection. The procedure was applied over Canada for the 1998‐2000 fire seasons and validated using fire surveys and burned area statistics from forest fire management agencies. The area of falsely mapped burns was found to be small (3.5% commission error over Canada), and most burns larger than 10 km2 were accurately detected and mapped (R² = 0.90, P<0.005, n = 91 for burns in two provinces). Canada‐wide satellite burned area was similar, but consistently smaller by comparison to statistics compiled by the Canadian Interagency Forest Fire Centre (by 17% in 1998, 16% in 1999, and 3% in 2000).     

Assessment of Post-fire Soil Erosion Risk in Fire-Affected Watersheds Using Remote Sensing and GIS

Soil erosion is a prominent cause of land degradation and desertification in Mediterranean countries. The detrimental effects of soil erosion are exemplified in climate (in particular climate change), topography, human activities, and natural disasters. Forest fires, which are an integral part of Mediterranean ecosystems, are responsible for the destruction of above-and below-ground vegetation that protects against soil erosion. Under this perspective, the estimation of potential soil erosion, especially after fire events, is critical for identifying watersheds that require management to prevent sediment loss, flooding, and increased ecosystem degradation. The objective of this study was to model the potential post-fire soil erosion risk following a large and intensive wildland fire, in order to prioritize protection and management actions at the watershed level in a Mediterranean landscape. Burn severity and preand post-fire land cover/uses were mapped using an ASTER image acquired two years before the fire, air photos acquired shortly after the fire, and a Landsat TM image acquired within one month after-fire. We estimated pre-and post-fire sediment loss using an integrated GIS-based approach, and additionally we analyzed landscape erosion patterns. The overall accuracy of the severity map reached 83%. Severe and heavy potential erosion classes covered approximately 90% of the total area following the fire, compared to 55% before. The fire had a profound effect on the spatial erosion pattern by altering the distribution of the potential erosion classes in 21 out of 24 watersheds, and seven watersheds were identified as being the most vulnerable to post-fire soil erosion. The spatial pattern of the erosion process is important because landscape cover heterogeneity induced especially by fire is a dominant factor controlling runoff generation and erosion rate, and should be considered in post-fire erosion risk assessment.     

Modeling of multi-strata forest fire severity using Landsat TM Data

Most of fire severity studies use field measures of composite burn index (CBI) to represent forest fire severity and fit the relationships between CBI and Landsat imagery derived differenced normalized burn ratio (dNBR) to predict and map fire severity at unsampled locations. However, less attention has been paid on the multi-strata forest fire severity, which represents fire activities and ecological responses at different forest layers. In this study, using field measured fire severity across five forest strata of dominant tree, intermediate-sized tree, shrub, herb, substrate layers, and the aggregated measure of CBI as response variables, we fit statistical models with predictors of Landsat TM bands, Landsat derived NBR or dNBR, image differencing, and image ratioing data. We model multi-strata forest fire in the historical recorded largest wildfire in California, the Big Sur Basin Complex fire. We explore the potential contributions of the post-fire Landsat bands, image differencing, image ratioing to fire severity modeling and compare with the widely used NBR and dNBR. Models using combinations of post-fire Landsat bands perform much better than NBR, dNBR, image differencing, and image ratioing. We predict and map multi-strata forest fire severity across the whole Big Sur fire areas, and find that the overall measure CBI is not optimal to represent multi-strata forest fire severity.     

Integrating modelling and remote sensing to identify ecosystem performance anomalies in the boreal forest,Yukon River Basin,Alaska

Abstract

High-latitude ecosystems are exposed to more pronounced warming effects than other parts of the globe. We develop a technique to monitor ecological changes in a way that distinguishes climate influences from disturbances. In this study, we account for climatic influences on Alaskan boreal forest performance with a data-driven model. We defined ecosystem performance anomalies (EPA) using the residuals of the model and made annual maps of EPA. Most areas (88%) did not have anomalous ecosystem performance for at least 6 of 8 years between 1996 and 2004. Areas with underperforming EPA (10%) often indicate areas associated with recent fires and areas of possible insect infestation or drying soil related to permafrost degradation. Overperforming areas (2%) occurred in older fire recovery areas where increased deciduous vegetation components are expected. The EPA measure was validated with composite burn index data and Landsat vegetation indices near and within burned areas.     

Burnt area mapping of Bandipur National Park,India using IRS 1C/1D LISS III data

The Bandipur National Park situated in the Western Ghats of Karnataka State, is one of the biodiversity hotspots of the world. During recent years, this park has witnessed repeated fires, affecting considerable areas under vegetation. The temporal satellite data from 1997 to 2006 have been analyzed to map the burnt areas using Remote Sensing (RS) and Geographic Information System (GIS) techniques. The vegetation cover is moist deciduous, dry deciduous, scrub forests and teak plantation. Information on extent of the burnt areas and the type of vegetation affected were derived forest range-wise. The fire prone regions have been identified by integrating vegetation type/density, road and settlement network and past history of forest fire occurrence, by assigning subjective weightage according to their fire-inducing capability or their sensitivity to fire. Comparison between each temporal dataset in terms of the extent of burnt area was also carried out to interpret fire incidence pattern. Three categories of fire risk regions such as Low, Moderate and High fire intensity zones were identified and it was found that almost 40% of the study area falls under low risk zone. An evaluation of the existing fire management systems and the implication of fire prevention programmes has been discussed, besides an assessment of causal factors for fire incidence in the park.     

Mapping burns and natural reforestation using thematic Mapper data

Remote sensing techniques are specially suitable to detect and to map areas affected by forest fires. In this work, Landsat 5 Thematic Mapper (TM) data has been used to study a number of forest fires that occurred in the province of Valencia (Spain) and to monitor the vegetation regeneration over burnt areas.

A reference area (non‐burnt forest) was established to assess the change produced by fire. The radiance in the thermal band (10.4–12.5 μm) and the normalized difference in reflectance between near 1R (0.76–0.90 μm) and middle IR (2.08–2.35 μm) were the most suitable parameters to map burnt areas. This index can also be used for monitoring vegetation regeneration in burnt areas. About a month after the fire, the burns show temperatures of 5–6 °C higher than those found in the reference area, and the vegetation index shows negative values whereas the reference area values remain positive. The differences between the burns and the reference area for the vegetation index decrease with time as vegetation regenerates.     

Carbon fluxes in ecosystems of Yellowstone National Park predicted from remote sensing data and simulation modeling

Background

A simulation model based on remote sensing data for spatial vegetation properties has been used to estimate ecosystem carbon fluxes across Yellowstone National Park (YNP). The CASA (Carnegie Ames Stanford Approach) model was applied at a regional scale to estimate seasonal and annual carbon fluxes as net primary production (NPP) and soil respiration components. Predicted net ecosystem production (NEP) flux of CO₂ is estimated from the model for carbon sinks and sources over multi-year periods that varied in climate and (wildfire) disturbance histories. Monthly Enhanced Vegetation Index (EVI) image coverages from the NASA Moderate Resolution Imaging Spectroradiometer (MODIS) instrument (from 2000 to 2006) were direct inputs to the model. New map products have been added to CASA from airborne remote sensing of coarse woody debris (CWD) in areas burned by wildfires over the past two decades.

Results

Model results indicated that relatively cooler and wetter summer growing seasons were the most favorable for annual plant production and net ecosystem carbon gains in representative landscapes of YNP. When summed across vegetation class areas, the predominance of evergreen forest and shrubland (sagebrush) cover was evident, with these two classes together accounting for 88% of the total annual NPP flux of 2.5 Tg C yr^-1 (1 Tg = 10¹² g) for the entire Yellowstone study area from 2000-2006. Most vegetation classes were estimated as net ecosystem sinks of atmospheric CO₂ on annual basis, making the entire study area a moderate net sink of about +0.13 Tg C yr^-1. This average sink value for forested lands nonetheless masks the contribution of areas burned during the 1988 wildfires, which were estimated as net sources of CO₂ to the atmosphere, totaling to a NEP flux of -0.04 Tg C yr^-1 for the entire burned area. Several areas burned in the 1988 wildfires were estimated to be among the lowest in overall yearly NPP, namely the Hellroaring Fire, Mink Fire, and Falls Fire areas.

Conclusions

Rates of recovery for burned forest areas to pre-1988 biomass levels were estimated from a unique combination of remote sensing and CASA model predictions. Ecosystem production and carbon fluxes in the Greater Yellowstone Ecosystem (GYE) result from complex interactions between climate, forest age structure, and disturbance-recovery patterns of the landscape.     

Mapping Wildfire Burn Severity in Southern California Forests and Shrublands Using Enhanced Thematic Mapper Imagery

Abstract

Wildfire is a major disturbance agent in Mediterranean Type Ecosystems (MTEs). Providing reliable, quantitative information on the area of burns and the level of damage caused is therefore important both for guiding resource management and global change monitoring. Previous studies have successfully mapped burn severity using remote sensing, but reliable accuracy has yet to be gained using standard methods over different vegetation types. The objective of this research was to classify burn severity across several vegetation types using Landsat ETM imagery in two areas affected by wildfire in southern California in June 1999. Spectral mixture analysis (SMA) using four reference endmembers (vegetation, soil, shade, non‐photosynthetic vegetation) and a single (charcoal‐ash) image endmember were used to enhance imagery prior to burn severity classification using decision trees. SMA provided a robust technique for enhancing fire‐affected areas due to its ability to extract sub‐pixel information and minimize the effects of topography on single date satellite data. Overall kappa classification accuracy results were high (0.71 and 0.85, respectively) for the burned areas, using five canopy consumption classes. Individual severity class accuracies ranged from 0.5 to 0.94.     

应用迁移学习的林火烈度初始评估研究

林火发生后,开展森林生态系统烈度信息的初始评估,能够为灾后生态修复管理措施的快速实施提供定量依据。为了改善传统林火烈度评估模型的时效性,本研究利用历史过火区域的实地调查数据,构建基于迁移学习的烈度评估模型,并将其应用于2020年3月30日发生的西昌泸山森林大火烈度初始评估研究中。研究结果表明：迁移学习算法能够将源区域和目标区域的遥感影像光谱转换为多个新的特征变量,在这些新特征变量构成的投影空间中,源区域和目标区域样本具有相似的特征分布。在此基础上,基于源区域历史实地调查数据构建的烈度评估模型,能够迁移应用于目标区域的烈度评估。在本研究林火烈度的初始评估中,基于迁移学习的烈度评估模型精度较高,总体精度为71.20%,Kappa系数为0.64。与该模型对比,未进行迁移学习的支持向量回归模型精度较低,其总体精度为58.00%,Kappa系数为0.48。同时,基于dNDVI、dLST和dNBR指数的经验回归模型精度最低,其总体精度分别为：20.80%、34.8%和24.80%,Kappa系数分别为：0.01、0.19和0.06。本研究可为林火灾后管理措施的快速响应,提供一种新的思路和参考。     

Forest Fire Burnt Area Assessment in the Biodiversity Rich Regions Using Geospatial Technology: Uttarakhand Forest Fire Event 2016

The hills of Uttarakhand witness forest fire every year during the summer season and the number of these fire events is reported to have increased due to increased anthropogenic disturbances as well as changes in climate. These fires cause significant damage to the natural resources which can be mapped and monitored using satellite images by virtue of its synoptic coverage of the landscape and near real time monitoring. This study presents burnt area assessment caused by the fire episode of April 2016 to the forest vegetation. Digital classification of satellite images was done to extract the burnt area which was found to be 3774.14 km², representing 15.28% of the total forest area of the state. It also gives an account of cumulative progression of forest fire in Uttarakhand using satellite images of three dates viz. 23rd, 27th May and 2nd June, 2016. Results were analyzed at district, administrative and forest division level using overlay analysis. Separate area statistics were given for different categories of biological richness, forest types and protected areas affected by forest fire. The burnt area assessment can be used in mitigation planning to prevent drastic ecological impacts of the forest fire on the landscape.     

A satellite data driven approach to monitoring and reporting fire disturbance and recovery across boreal and temperate forests

The regular and consistent measurements provided by Earth observation satellites can support the monitoring and reporting of forest indicators. Although substantial scientific literature espouses the capabilities of satellites in this area, the techniques are under-utilised in national reporting, where there is a preference for aggregating ad hoc data. In this paper, we posit that satellite information, while perhaps of low accuracy at single time steps or across small areas, can produce trends and patterns which are, in fact, more meaningful at regional and national scales. This is primarily due to data consistency over time and space. To investigate this, we use MODIS and Landsat data to explore trends associated with fire disturbance and recovery across boreal and temperate forests worldwide. Our results found that 181 million ha (9 %) of the study area (2 billion ha of forests) was burned between 2001 and 2018, as detected by MODIS satellites. World Wildlife Fund biomes were used for a detailed analysis across several countries. A significant increasing trend in area burned was observed in Mediterranean forests in Chile (8.9 % yr⁻¹), while a significant decreasing trend was found in temperate mixed forests in China (-2.2 % yr⁻¹). To explore trends and patterns in fire severity and forest recovery, we used Google Earth Engine to efficiently sample thousands of Landsat images from 1991 onwards. Fire severity, as measured by the change in the normalized burn ratio (NBR), was found to be generally stable over time; however, a slight increasing trend was observed in the Russian taiga. Our analysis of spectral recovery following wildfire indicated that it was largely dependent on location, with some biomes (particularly in the USA) showing signs that spectral recovery rates have shortened over time. This study demonstrates how satellite data and cloud-computing can be harnessed to establish baselines and reveal trends and patterns, and improve monitoring and reporting of forest indicators at national and global scales.