Our novel scheme enhances large-eddy simulations (LESs) for atmosphere–land interactions. It couples LES subgrid closure with Monin–Obukhov similarity theory (MOST), overcoming MOST's limitations. Validated over diverse land surfaces, our approach outperforms existing methods, aligning well with field measurements. Robustness is demonstrated across varying model resolutions. MOST's influence strengthens with decreasing grid spacing, particularly for sensible heat flux.