结合中尺度大气模式MM5和Shao的风蚀起沙方案,对敦煌莫高窟顶灌木林带和草方格沙障的固沙效应进行数值模拟研究。结果表明:灌木林带的建立,使其地表风蚀起沙的临界摩擦速度(ut)达到1.15m/s,是林带建立前的4.1倍,地表的顺风向沙粒通量(Q)和垂直尘粒通量(F)均减小为0,完全抑制了地表风蚀起沙;草方格沙障的设置,使ut达到0.53m/s,为设置前的1.9倍,Q和F分别减少约3/5和7/8,固沙效果非常明显。一般情况下,ut不会大于1.0m/s,充分说明灌木林带比草方格具有更强的固沙效果。MM5模拟的摩擦速度和观测的摩擦速度在变化趋势上呈较好的线性相关,相关系数达0.8以上,说明MM5对该区域摩擦速度的模拟结果可靠。