Tree diameter distributions are essential for the calculation of stem volume and biomass, as well as simulation of growth and yield and to understand timber assortments. Accurate and reliable prediction of tree diameter distributions is critical for optimizing forest structure compositions, scheduling silvicultural operations and promoting sustainable management. In this study, we investigated the potential of airborne Light Detection and Ranging (LiDAR) data for predicting tree diameter distributions using a bimodal finite mixture model (FMM) and a multimodal k-nearest neighbor (KNN) model (compared to the unimodal Weibull model (UWM)) over a subtropical planted forest in southern China. To do so, we first evaluated the capability of various LiDAR predictions (i.e., the bimodality coefficient (BC) and Lorenz-based indicators) to stratify forest structural types into unimodal and multimodal stands. Once the best LiDAR prediction for the differentiation was determined, the parameters of UWM (in non-specific and species-specific models) and FMM (in structure-specific models) were estimated by LiDAR-derived metrics and the tree diameter distributions of stands were generated by the estimated LiDAR parameters. When KNN was applied for constructing diameter distributions, optimal KNN strategies, including number of neighbors k, response configurations and imputation methods (i.e., Most Similar Neighbor (MSN) and Random Forest (RF)) for different species were heuristically determined. Finally, the predictive performance of estimated LiDAR the parameters of UWM, FMM and KNN for predicting diameter distributions were assessed. The results showed that LiDAR-predicted Lorenz-based indicators performed best for differentiation. Parameters of UWM and FMM were predicted well and the species-specific models had higher accuracies than the non-specific models. Overall, RF imputation from KNN with an optimal response set (i.e., DBH) were was stable than MSN imputation when k = 5 neighbors. In addition, the inclusion of bimodal FMM for differentiated all plots generally produced a more accurate result (Mean eR = 40.85, Mean eP = 0.20) than multimodal KNN (Mean eR = 52.19, Mean eP = 0.26), whereas the UWM produced the lowest performance (Mean eR = 52.31, Mean eP = 0.26). This study demonstrated the benefits of multimodal models with LiDAR for estimating diameter distributions for supporting forest inventory and sustainable forest management in subtropical planted forests.