Abstract Fourth-order-centered operators on regular hexagonal grids with the ZM-grid arrangement are described. The finite-volume method is used and operators are defined at hexagonal cell centers. The gradient operator is calculated from 12 surrounding cell center scalars. The divergence operator is defined from 12 surrounding cell corner vectors. A linear combination of local or interpolated values generates cell corner values used to calculate the operators. The flux-divergence operator applies the same cell corner values as those used in the gradient and divergence operators. The fourth-order convergence of the gradient and divergence operators is obtained in numerical tests using sufficiently smooth and differentiable test functions. The flux-divergence operator is formally second-order accurate. However, the results from a cone advection test show that the flux-divergence operator performs better than a commonly used second-order flux-divergence operator. Numerical dispersion and phase error are small because mean wind advection is computed with fourth-order accuracy.