Source author record

K. Shafer Smith

K. Shafer Smith appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.ao-ph physics.comp-ph

Catalog footprint

What is connected

2works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2019arXiv

Fast and Accurate Computation of Vertical Modes

The vertical modes of linearized equations of motion are widely used by the oceanographic community in numerous theoretical and observational contexts. However, the standard approach for solving the generalized eigenvalue problem using second-order finite difference matrices produces $O(1)$ errors for all but the few lowest modes, and increasing resolution quickly becomes too slow as the computational complexity of eigenvalue algorithms increase as $O(n^3)$. Existing methods are therefore inadequate for computing a full spectrum of internal waves, such as needed for initializing a numerical model with a full internal wave spectrum. Here we show that rewriting the eigenvalue problem in stretched coordinates and projecting onto Chebyshev polynomials results in substantially more accurate modes than finite-differencing at a fraction of the computational cost. We also compute the surface quasigeostrophic modes using the same methods. All spectral and finite difference algorithms are made available in a suite of Matlab classes that have been validated against known analytical solutions in constant and exponential stratification.

preprint2014arXiv

The effect of surface buoyancy gradients on oceanic Rossby wave propagation

Motivated by the discrepancy between satellite observations of coherent westward propagating surface features and Rossby wave theory, this paper revisits the planetary wave propagation problem, taking into account the effects of lateral buoyancy gradients at the ocean's surface. The standard theory for long baroclinic Rossby waves is based on an expansion of the quasigeostrophic stretching operator in normal modes, $ϕ_n(z)$, satisfying a Neumann boundary condition at the surface, $ϕ_n'(0) = 0$. Buoyancy gradients are, by thermal wind balance, proportional to the vertical derivative of the streamfunction, thus such modes are unable to represent ubiquitous lateral buoyancy gradients in the ocean's mixed layer. Here, we re-derive the wave propagation problem in terms of an expansion in a recently-developed "surface-aware" (SA) basis that can account for buoyancy anomalies at the ocean's surface. The problem is studied in the context of an idealized Charney-like baroclinic wave problem set in an oceanic context, where a surface mean buoyancy gradient interacts with a constant interior potential vorticity gradient that results from both $β$ and the curvature of the mean shear. The wave frequencies, growth rates and phases are systematically compared to those computed from a two-layer model, a truncated expansion in standard baroclinic modes and to a high-vertical resolution calculation that represents the true solution. The full solution generally shows faster wave propagation when lateral surface gradients are present. Moreover, the wave problem in the SA basis best captures the full solution, even with just a two or three modes.