The Mediterranean forests, particularly Cork oak (Quercus suber L., 1927), make a major contribution to the fight against climate change through Carbon sequestration. Hence, there is a great interest in the accurate quantification of biomass and carbon stock. In this context, this study aims at assessing the performance of a new approach, based on the combination of Unmanned aerial vehicle airborne Aerial laser scanning (ALS-UAV) and Terrestrial laser scanning (TLS) data, in the determination of dendrometric parameters (Circumference at 1.30 m and Tree Height), and consequently the estimation of biomass and carbon stock, considering field data as reference. This study takes the Maamora forest in Morocco as an example of a Mediterranean Cork oak forest. The methodology consists of collecting data at three levels: the entire area level for an ALS-UAV scan, the plot and tree levels for TLS surveys, as well as field data collection. Afterwards, dendrometric parameters (Circumference at 1.30 m and the Tree height) were estimated using individual tree segmentation and biomass; the carbon stock (aboveground, belowground, and total) was estimated using allometric equations. The comparison of the estimated dendrometric parameters with those measured in the field shows a strong relationship, with a Pearson coefficient of 0.86 and 0.83, a correlation coefficient (R²) of 0.81 and 0.71, and a Root mean square error (RMSE) of 1.84 cm and 0.47 m, respectively. Concerning the biomass and carbon stock estimation, the proposed approach gives a satisfactory accuracy, with a Pearson coefficient of 0.77, an R2 of 0.83, and an RMSE of 36.40 kg for biomass and 20.24 kg for carbon stock.

Spatialization of biomass and carbon stocks is essential for a good understanding of the forest stand and its characteristics, especially in degraded Mediterranean cork oak forests. Furthermore, the analysis of biomass and carbon stock changes and dynamics is essential for understanding the carbon cycle, in particular carbon emissions and stocks, in order to make projections, especially in the context of climate change. In this research, we use a multidimensional framework integrating forest survey data, LiDAR UAV data, and extracted vegetation indices from Landsat imagery (NDVI, ARVI, CIG, etc.) to model and spatialize cork oak biomass and carbon stocks on a large scale. For this purpose, we explore the use of univariate and multivariate regression modeling and examine several types of regression, namely, multiple linear regression, stepwise linear regression, random forest regression, simple linear regression, logarithmic regression, and quadratic and cubic regression. The results show that for multivariate regression, stepwise regression gives good results, with R² equal to 80% and 65% and RMSE equal to 2.59 and 1.52 Mg/ha for biomass and carbon stock, respectively. Random forest regression, chosen as the ML algorithm, gives acceptable results, explaining 80% and 60% of the variation in biomass and carbon stock, respectively, and an RMSE of 2.74 and 1.72 Mg/ha for biomass and carbon stock, respectively. For the univariate regression, the simple linear regression is chosen because it gives satisfactory results, close to those of the quadratic and cubic regressions, but with a simpler equation. The vegetation index chosen is ARVI, which shows good performance indices, close to those of the NDVI and CIG. The assessment of biomass and carbon stock changes in the study area over 35 years (1985–2020) showed a slight increase of less than 10 Mg/ha and a decrease in biomass and carbon stock over a large area.