Eloisa Raluy-López, Domingo Muñoz-Esparza, Juan Pedro Montávez
Abstract. High-resolution simulations at both mesoscale and microscale increasingly rely on detailed terrain datasets, but terrain-following coordinate models can suffer from numerical instabilities in steep-slope regions. To address this issue, terrain smoothing is typically applied in numerical weather prediction models, though conventional global smoothing unnecessarily reduces resolution across the entire domain. This study presents a localized terrain smoothing approach designed to prevent numerical instabilities while preserving terrain details. Different smoothing strategies were tested for efficiency, computational cost, and terrain preservation. The final approach applies a Gaussian filter with adaptive standard deviation within a localized 3 × 3 grid, with a blending factor of 0.2, and treating all the steep-slope points simultaneously. Integrated into the NCAR's FastEddy® LES and WRF mesoscale community models, this technique effectively prevents terrain-driven instabilities in high-resolution simulations over complex terrain. The proposed local filtering method helps minimizing loss of terrain detail and avoiding the need for excessively strong numerical filtering during run time to stabilize the simulations. This method is computationally efficient, easy to implement, and adaptable to other models, providing a robust solution to improve numerical stability while maintaining high-resolution terrain features.