AbstractA three‐dimensional positivity‐preserving and conservative transport model on a cubed‐sphere grid is presented using a multimoment finite‐volume (MMFV) method. The three‐point fourth‐order MMFV method with boundary gradient switching projection is utilized to ensure nonoscillatory numerical solutions in the horizontal, while the simple and computationally efficient semi‐Lagrangian piecewise rational method is used to allow a large time step in the vertical. By means of horizontal and vertical separation, a Strang‐type splitting technique is adopted to advance the numerical solutions in time. Non‐negative corrections are made to achieve positive‐definite preservation exactly. The resulting transport model is inherently conservative, high‐accuracy, and positive‐definite, and allows a large Courant number no matter how dense the vertical grid is. Various benchmark test cases, including three representative advection tests and an idealized terminator chemistry test, are conducted to validate the performance of the model presented. The numerical results indicate that the model is quite competitive in comparison with the existing models in the literature and looks very promising in real atmospheric simulations.