Forest carbon stock is an important indicator reflecting a forest ecosystem’s structures and functions. Its spatial distribution is valuable for managing natural resources, protecting ecosystems and biodiversity, and further promoting sustainability, but accurately mapping the forest carbon stock distribution in a large area is a challenging task. This study selected Changting County, Fujian Province, as a case study to explore a method to map the forest carbon stock distribution using the integration of airborne Lidar, Sentinel-2, and ancillary data in 2022. The Bayesian hierarchical modeling approach was used to estimate the local forest carbon stock based on airborne Lidar data and field measurements, and then the random forest approach was used to develop a regional forest carbon stock estimation model based on the Sentinel-2 and ancillary data. The results indicated that the Lidar-based carbon stock distribution effectively provided sample plots with good spatial representativeness for modeling regional carbon stock with a coefficient of determination (R2) of 0.7 and root mean square error (RMSE) of 12.94 t/ha. The average carbon stocks were 48.55 t/ha, 55.51 t/ha, and 57.04 t/ha for Masson pine, Chinese fir, and broadleaf forests, respectively. The carbon stock in non-conservation regions was 15.2–16.1 t/ha higher than that in conservation regions. This study provides a promising method through the use of airborne Lidar data as a linkage between sample plots and Sentinel-2 data to map the regional carbon stock distribution in those subtropical regions where serious soil erosion has led to a relatively sparse forest canopy density. The results are valuable for local government to make scientific decisions for promoting ecosystem restoration due to water and soil erosion.