Yamaguchi and Feingold (2012) note that the cloud fields in their Weather Research and Forecasting (WRF) large-eddy simulations (LESs) of marine stratocumulus exhibit a strong sensitivity to time stepping choices. In this study, we reproduce and analyze this sensitivity issue using two stratocumulus cases, one marine and one continental. Results show that (1) the sensitivity is associated with spurious motions near the moisture jump between the boundary layer and the free atmosphere, and (2) these spurious motions appear to arise from neglecting small variations in water vapor mixing ratio (qv) in the pressure gradient calculation in the acoustic substepping portion of the integration procedure. We show that this issue is remedied in the WRF dynamical core by replacing the prognostic equation for the potential temperature ? with one for the moist potential temperature ?m=?(1+1.61qv), which allows consistent treatment of moisture in the calculation of pressure during the acoustic substeps. With this modification, the spurious motions and the sensitivity to the time stepping settings (i.e., the dynamic time step length and number of acoustic substeps) are eliminated in both of the example stratocumulus cases. This modification improves the applicability of WRF for LES applications, and possibly other models using similar dynamical core formulations, and also permits the use of longer time steps than in the original code.
Revised: March 18, 2016 |
Published: October 29, 2015
Citation
Xiao H., S. Endo, M. Wong, W. Skamarock, J.B. Klemp, J.D. Fast, and W.I. Gustafson, et al. 2015.Modifications to WRF’s dynamical core to improve the treatment of moisture for large-eddy simulations.Journal of Advances in Modeling Earth Systems 7, no. 4:1627-1642.PNNL-SA-112067.doi:10.1002/2015MS000532