Evaluation of Scalar Advection Schemes in the Advanced Research WRF Model Using Large-Eddy Simulations of Aerosol–Cloud Interactions

Abstract In the Advanced Research Weather Research and Forecasting Model (ARW), versions 3.0 and earlier, advection of scalars was performed using the Runge–Kutta time-integration scheme with an option of using a positive-definite (PD) flux limiter. Large-eddy simulations of aerosol–cloud interactions using the ARW model are performed to evaluate the advection schemes. The basic Runge–Kutta scheme alone produces spurious oscillations and negative values in scalar mixing ratios because of numerical dispersion errors. The PD flux limiter assures positive definiteness but retains the oscillations with an amplification of local maxima by up to 20% in the tests. These numerical dispersion errors contaminate active scalars directly through the advection process and indirectly through physical and dynamical feedbacks, leading to a misrepresentation of cloud physical and dynamical processes. A monotonic flux limiter is introduced to correct the generally accurate but dispersive solutions given by high-order Runge–Kutta scheme. The monotonic limiter effectively minimizes the dispersion errors with little significant enhancement of numerical diffusion errors. The improvement in scalar advection using the monotonic limiter is discussed in the context of how the different advection schemes impact the quantification of aerosol–cloud interactions. The PD limiter results in 20% (10%) fewer cloud droplets and 22% (5%) smaller cloud albedo than the monotonic limiter under clean (polluted) conditions. Underprediction of cloud droplet number concentration by the PD limiter tends to trigger the early formation of precipitation in the clean case, leading to a potentially large impact on cloud albedo change.