Forest disturbance regimes across biomes are being altered by interactive effects of global change. Establishing baselines for assessing change requires detailed quantitative data on past disturbance events, but such data are scarce and difficult to obtain over large spatial and temporal scales. The integration of remote sensing with dense time series analysis and cloud computing platforms is enhancing the ability to monitor historical disturbances, and especially non-stand replacing events along climatic gradients. Since the integration of such tools is still scarce in Mediterranean regions, here, we combine dense Landsat time series and the Continuous Change Detection and Classification - Spectral Mixture Analysis (CCDC-SMA) method to monitor forest disturbance in continental Spain from 1985 to 2023. We adapted the CCDC-SMA method for improved disturbance detection creating new spectral libraries representative of the study region, and quantified the year, month, severity, return interval, and type of disturbance (stand replacing, non-stand replacing) at a 30 m resolution. In addition, we characterised forest disturbance regimes and trends (patch size and severity, and frequency of events) of events larger than 0.5 ha at the national scale by biome (Mediterranean and temperate) and forest type (broadleaf, needleleaf and mixed). We quantified more than 2.9 million patches of disturbed forest, covering 4.6 Mha over the region and period studied. Forest disturbances were on average larger but less severe in the Mediterranean than in the temperate biome, and significantly larger and more severe in needleleaf than in mixed and broadleaf forests. Since the late 1980s, forest disturbances have decreased in size and severity while increasing in frequency across all biomes and forest types. These results have important implications as they confirm that disturbance regimes in continental Spain are changing and should therefore be considered in forest strategic planning for policy development and implementation.

Forests provide habitat for diverse species and play a key role in mitigating climate change. Remote sensing enables efficient monitoring of many forest attributes across vast areas, thus supporting effective and efficient management strategies. This study aimed to identify an effective combination of remote sensing sensors for estimating biophysical variables in Mediterranean and temperate forests that can be easily translated into an operational context. Aboveground biomass (AGB), canopy height (CH), and forest canopy cover (FCC) were estimated using a combination of optical (Sentinel-2, Landsat) and radar sensors (Sentinel-1 and TerraSAR-X/TanDEM-X), along with records of past forest disturbances and topography-related variables. As a reference, lidar-derived AGB, CH, and FCC were used. Model performance was assessed not only with standard approaches such as out-of-bag sampling but also with completely independent lidar-derived reference datasets, thus enabling evaluation of the model’s temporal inference capacity. In Mediterranean forests, models based on optical imagery outperformed the radar-enhanced models when estimating FCC and CH, with elevation and spectral indices being key predictors of forest structure. In contrast, in denser temperate forests, radar data (especially X-band relative heights) were crucial for estimating CH and AGB. Incorporating past disturbance data further improved model accuracy in these denser ecosystems. Overall, this study underscores the value of integrating multi-source remote sensing data while highlighting the limitations of temporal extrapolation. The presented methodology can be adapted to enhance forest variable estimation across many forest ecosystems.