Fully Lagrangian numerical solutions of unbalanced frontogenesis and frontal collapse

Abstract

Numerical simulation has failed to answer some fundamental questions about atmospheric frontogenesis because of the artificial minimum resolved scale in grid point and spectral models alike. A fully Lagrangian primitive-equation numerical model is developed for nonturbulent, slab-symmetric flow on an f-plane. With physical position treated as an explicit function of particle label and time, and model grid deforms to follow natural changes in disturbance length scales. Exact conservation of volume and potential vorticity, as well as of basic tracer variables, is demonstrated, and details of the truncation error for energy conservation are obtained for the case of second-order central differencing in label space. The Lagrangian model is used to simulate frontogenesis by horizontal wind deformation in a dry, Boussinesq atmosphere, with no prior assumption of hydrostatic or geostrophic balance. For realistic choices of the parameters governing the rate of frontogenesis, imbalances alone are found to be insufficient to prevent frontal collapse. For small values of the normalized potential vorticity, the ageostrophic secondary circulation is weaker than in the corresponding balanced solutions, and frontal collapse is accordingly delayed. -from Author