Understanding and mapping forest fire vulnerability is essential for informed landscape management and disaster risk reduction, especially in the context of increasing anthropogenic and climatic pressures. This study aims to model and spatially predict forest fire vulnerability across Romania using two machine learning algorithms: MaxEnt and XGBoost. We integrated forest fire occurrence data from 2006 to 2024 with a suite of climatic, topographic, ecological, and anthropogenic predictors at a 250 m spatial resolution. MaxEnt, based on presence-only data, achieved moderate predictive performance (AUC = 0.758), while XGBoost, trained on presence–absence data, delivered higher classification accuracy (AUC = 0.988). Both models revealed that the impact of environmental variables on forest fire occurrence is complex and heterogeneous, with the most influential predictors being the Fire Weather Index, forest fuel type, elevation, and distance to human proximity features. The resulting vulnerability and uncertainty maps revealed hotspots in Sub-Carpathian and lowland regions, especially in Mehedinți, Gorj, Dolj, and Olt counties. These patterns reflect historical fire data and highlight the role of transitional agro-forested landscapes. This study delivers a replicable, data-driven approach to wildfire risk modelling, supporting proactive management and emphasising the importance of integrating vulnerability assessments into planning and climate adaptation strategies.

The objective of this study was to analyze the efficiency of individual tree identification and stand volume estimation from LiDAR data. The study was located in Norway spruce [Picea abies (L.) Karst.] stands in southwestern Romania and linked airborne laser scanning (ALS) with terrestrial measurements through empirical modelling. The proposed method uses the Canopy Maxima algorithm for individual tree detection together with biometric field measurements and individual trees positioning. Field data was collected using Field-Map real-time GIS-laser equipment, a high-accuracy GNSS receiver and a Vertex IV ultrasound inclinometer. ALS data were collected using a Riegl LMS-Q560 instrument and processed using LP360 and Fusion software to extract digital terrain, surface and canopy height models. For the estimation of tree heights, number of trees and tree crown widths from the ALS data, the Canopy Maxima algorithm was used together with local regression equations relating field-measured tree heights and crown widths at each plot. When compared to LiDAR detected trees, about 40-61% of the field-measured trees were correctly identified. Such trees represented, in general, predominant, dominant and co-dominant trees from the upper canopy. However, it should be noted that the volume of the correctly identified trees represented 60-78% of the total plot volume. The estimation of stand volume using the LiDAR data was achieved by empirical modelling, taking into account the individual tree heights (as identified from the ALS data) and the corresponding ground reference stem volume. The root mean square error (RMSE) between the individual tree heights measured in the field and the corresponding heights identified in the ALS data was 1.7-2.2 meters. Comparing the ground reference estimated stem volume (at trees level) with the corresponding ALS estimated tree stem volume, an RMSE of 0.5-0.7 m³ was achieved. The RMSE was slightly lower when comparing the ground reference stem volume at plot level with the ALS-estimated one, taking into account both the identified and unidentified trees in the LiDAR data (0.4-0.6 m³).