Toward a PV-based algorithm for the dynamical core of hydrostatic global models

The diabatic contour-advective semi-Lagrangian (DCASL) algorithms previously constructed for the shallow-water and multilayer Boussinesq primitive equations are extended to multilayer non-Boussinesq equations on the sphere using a hybrid terrain-following-isentropic (sigma-) vertical coordinate. It is shown that the DCASL algorithms face challenges beyond more conventional algorithms in that various types of damping, filtering, and regularization are required for computational stability, and the nonlinearity of the hydrostatic equation in the sigma- coordinate causes convergence problems with setting up a semi-implicit time-stepping scheme to reduce computational cost. The prognostic variables are an approximation to the Rossby-Ertel potential vorticity Q, a scaled pressure thickness, the horizontal divergence, and the surface potential temperature. Results from the DCASL algorithm in two formulations of the sigma- coordinate, differing only in the rate at which the vertical coordinate tends to with increasing height, are assessed using the baroclinic instability test case introduced by Jablonowski and Williamson in 2006. The assessment is based on comparisons with available reference solutions as well as results from two other algorithms derived from the DCASL algorithm: one with a semi-Lagrangian solution for Q and another with an Eulerian grid-based solution procedure with relative vorticity replacing Q as the prognostic variable. It is shown that at intermediate resolutions, results comparable to the reference solutions can be obtained.Publisher PDFPeer reviewe

The combined Lagrangian advection method

We present and test a new hybrid numerical method for simulating layerwise-two-dimensional geophysical flows. The method radically extends the original Contour-Advective Semi-Lagrangian (CASL) algorithm by combining three computational elements for the advection of general tracers (e.g. potential vorticity, water vapor, etc.): (1) a pseudospectral method for large scales, (2) Lagrangian contours for intermediate to small scales, and (3) Lagrangian particles for the representation of general forcing and dissipation. The pseudo-spectral method is both efficient and highly accurate at large scales, while contour advection is efficient and accurate at small scales, allowing one to simulate extremely finescale structure well below the basic grid scale used to represent the velocity field. The particles allow one to efficiently incorporate general forcing and dissipation

Remote sensing of mesospheric dust layers using active modulation of PMWE by high-power radio waves

This paper presents the first study of the modulation of polar mesospheric winter echoes (PMWE) by artificial radio wave heating using computational modeling and experimental observation in different radar frequency bands. The temporal behavior of PMWE response to HF pump heating can be employed to diagnose the charged dust layer associated with mesospheric smoke particles. Specifically, the rise and fall time of radar echo strength as well as relaxation and recovery time after heater turn-on and turnoff are distinct parameters that are a function of radar frequency. The variation of PMWE strength with PMWE source region parameters such as electron-neutral collision frequency, photodetachment current, electron temperature enhancement ratio, dust density, and radius is considered. The comparison of recent PMWE measurements at 56 MHz and 224 MHz with computational results is discussed, and dust parameters in the PMWE generation regime are estimated. Predictions for HF PMWE modification and its connection to the dust charging process by free electrons is investigated. The possibility for remote sensing of dust and plasma parameters in artificially modified PMWE regions using simultaneous measurements in multiple frequency bands are discussed. © 2016. American Geophysical Union. All Rights Reserved

Geostrophic balance preserving interpolation in mesh adaptive shallow-water ocean modelling

The accurate representation of geostrophic balance is an essential requirement for numerical modelling of geophysical flows. Significant effort is often put into the selection of accurate or optimal balance representation by the discretisation of the fundamental equations. The issue of accurate balance representation is particularly challenging when applying dynamic mesh adaptivity, where there is potential for additional imbalance injection when interpolating to new, optimised meshes. In the context of shallow-water modelling, we present a new method for preservation of geostrophic balance when applying dynamic mesh adaptivity. This approach is based upon interpolation of the Helmholtz decomposition of the Coriolis acceleration. We apply this in combination with a discretisation for which states in geostrophic balance are exactly steady solutions of the linearised equations on an f-plane; this method guarantees that a balanced and steady flow on a donor mesh remains balanced and steady after interpolation onto an arbitrary target mesh, to within machine precision. We further demonstrate the utility of this interpolant for states close to geostrophic balance, and show that it prevents pollution of the resulting solutions by imbalanced perturbations introduced by the interpolation

Correlations of control variables in variational data assimilation

Variational data assimilation systems for numerical weather prediction rely on a transformation of model variables to a set of control variables that are assumed to be uncorrelated. Most implementations of this transformation are based on the assumption that the balanced part of the flow can be represented by the vorticity. However, this assumption is likely to break down in dynamical regimes characterized by low Burger number. It has recently been proposed that a variable transformation based on potential vorticity should lead to control variables that are uncorrelated over a wider range of regimes. In this paper we test the assumption that a transform based on vorticity and one based on potential vorticity produce an uncorrelated set of control variables. Using a shallow-water model we calculate the correlations between the transformed variables in the different methods. We show that the control variables resulting from a vorticity-based transformation may retain large correlations in some dynamical regimes, whereas a potential vorticity based transformation successfully produces a set of uncorrelated control variables. Calculations of spatial correlations show that the benefit of the potential vorticity transformation is linked to its ability to capture more accurately the balanced component of the flow

CABARET in the ocean gyres

Author Posting. © The Author(s), 2009. This is the author's version of the work. It is posted here by permission of Elsevier B.V. for personal use, not for redistribution. The definitive version was published in Ocean Modelling 30 (2009): 155-168, doi:10.1016/j.ocemod.2009.06.009.A new high-resolution Eulerian numerical method is proposed for modelling quasigeostrophic ocean dynamics in eddying regimes. The method is based on a novel, second-order non-dissipative and lowdispersive conservative advection scheme called CABARET. The properties of the new method are compared with those of several high-resolution Eulerian methods for linear advection and gas dynamics. Then, the CABARET method is applied to the classical model of the double-gyre ocean circulation and its performance is contrasted against that of the common vorticity-preserving Arakawa method. In turbulent regimes, the new method permits credible numerical simulations on much coarser computational grids.Supports from the Royal Society of London and from the Mary Sears Visitor Grant are acknowledged by SK with gratitude. The work of VG was supported by the Russian Foundation for Basic Research (RFBR), grant 06-01-00819a. Funding for PB was provided by the NSF grant 0725796