Source author record

Ian Grooms

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

8 published item(s)

preprint2022arXiv

A comparison of nonlinear extensions to the ensemble Kalman filter: Gaussian Anamorphosis and Two-Step Ensemble Filters

Ensemble Kalman filters are based on a Gaussian assumption, which can limit their performance in some non-Gaussian settings. This paper reviews two nonlinear, non-Gaussian extensions of the Ensemble Kalman Filter: Gaussian anamorphosis (GA) methods and two-step updates, of which the rank histogram filter (RHF) is a prototypical example. GA-EnKF methods apply univariate transforms to the state and observation variables to make their distribution more Gaussian before applying an EnKF. The two-step methods use a scalar Bayesian update for the first step, followed by linear regression for the second step. The connection of the two-step framework to the full Bayesian problem is made, which opens the door to more advanced two-step methods in the full Bayesian setting. A new method for the first part of the two-step framework is proposed, with a similar form to the RHF but a different motivation, called the `improved RHF' (iRHF). A suite of experiments with the Lorenz-`96 model demonstrate situations where the GA-EnKF methods are similar to EnKF, and where they outperform EnKF. The experiments also strongly support the accuracy of the RHF and iRHF filters for nonlinear and non-Gaussian observations; these methods uniformly beat the EnKF and GA-EnKF methods in the experiments reported here. The new iRHF method is only more accurate than RHF at small ensemble sizes in the experiments reported here.

preprint2020arXiv

A fast tunable blurring algorithm for scattered data

A blurring algorithm with linear time complexity can reduce the small-scale content of data observed at scattered locations in a spatially extended domain of arbitrary dimension. The method works by forming a Gaussian interpolant of the input data, and then convolving the interpolant with a multiresolution Gaussian approximation of the Green's function to a differential operator whose spectrum can be tuned for problem-specific considerations. Like conventional blurring algorithms, which the new algorithm generalizes to data measured at locations other than a uniform grid, applications include deblurring and separation of spatial scales. An example illustrates a possible application toward enabling importance sampling approaches to data assimilation of geophysical observations, which are often scattered over a spatial domain, since blurring observations can make particle filters more effective at state estimation of large scales. Another example, motivated by data analysis of dynamics like ocean eddies that have strong separation of spatial scales, uses the algorithm to decompose scattered oceanographic float measurements into large-scale and small-scale components.

preprint2016arXiv

Investigations of non-hydrostatic, stably stratified and rapidly rotating flows

We present an investigation of rapidly rotating (small Rossby number $Ro\ll 1$) and stratified turbulence where the stratification strength is varied from weak (large Froude number $Fr\gg1$) to strong ($Fr\ll1$). The investigation is set in the context of a reduced model derived from the Boussinesq equations that efficiently retains anisotropic inertia-gravity waves with order-one frequencies and highlights a regime of wave-eddy interactions. Numerical simulations of the reduced model are performed where energy is injected by a stochastic forcing of vertical velocity, which forces wave modes only. The simulations reveal two regimes characterized by the presence of well-formed, persistent and thin turbulent layers of locally-weakened stratification at small Froude numbers, and by the absence of layers at large Froude numbers. Both regimes are characterized by a large-scale barotropic dipole enclosed by small-scale turbulence. When the Reynolds number is not too large a direct cascade of barotropic kinetic energy is observed, leading to total energy equilibration. We examine net energy exchanges that occur through vortex stretching and vertical buoyancy flux and diagnose the horizontal scales active in these exchanges. We find that the baroclinic motions inject energy directly to the largest scales of the barotropic mode, implying that the large-scale barotropic dipole is not the end result of an inverse cascade within the barotropic mode.

preprint2015arXiv

On Gakerkin approximations for the surface-active quasigeostrophic equations

We study the representation of solutions of the three-dimensional quasigeostrophic (QG) equations using Galerkin series with standard vertical modes, with particular attention to the incorporation of active surface buoyancy dynamics. We extend two existing Galerkin approaches (A and B) and develop a new Galerkin approximation (C). Approximation A, due to \cite{flierl1978}, represents the streamfunction as a truncated Galerkin series and defines the potential vorticity (PV) that satisfies the inversion problem exactly. Approximation B, due to \cite{tulloch_smith2009b}, represents the PV as a truncated Galerkin series and calculates the streamfunction that satisfies the inversion problem exactly. Approximation C, the true Galerkin approximation for the QG equations, represents both streamfunction and PV as truncated Galerkin series, but does not satisfy the inversion equation exactly. The three approximations are fundamentally different unless the boundaries are isopycnal surfaces. We discuss the advantages and limitations of approximations A, B, and C in terms of mathematical rigor and conservation laws, and illustrate their relative efficiency by solving linear stability problems with nonzero surface buoyancy. With moderate number of modes, B and C have have superior accuracy than A at high wavenumbers. Because B lacks conservation of energy, we recommend approximation C for constructing solutions to the surface-active QG equations using Galerkin series with standard vertical modes.

preprint2014arXiv

Asymptotic Behavior of Heat Transport for a Class of Exact Solutions in Rotating Rayleigh-Bénard Convection

The non-hydrostatic, quasigeostrophic approximation for rapidly rotating Rayleigh-Bénard convection admits a class of exact `single mode' solutions. These solutions correspond to steady laminar convection with a separable structure consisting of a horizontal planform characterized by a single wavenumber multiplied by a vertical amplitude profile, with the latter given as the solution of a nonlinear boundary value problem. The heat transport associated with these solutions is studied in the regime of strong thermal forcing (large reduced Rayleigh number $\widetilde{Ra}$). It is shown that the Nusselt number $Nu$, a nondimensional measure of the efficiency of heat transport by convection, for this class of solutions is bounded below by $Nu\gtrsim \widetilde{Ra}^{3/2}$, independent of the Prandtl number, in the limit of large reduced Rayleigh number. Matching upper bounds include only logarithmic corrections, showing the accuracy of the estimate. Numerical solutions of the nonlinear boundary value problem for the vertical structure are consistent with the analytical bounds.

preprint2014arXiv

Bounds on Heat Transport in Rapidly Rotating Rayleigh-Bénard Convection

The heat transport in rotating Rayleigh-Bénard convection is considered in the limit of rapid rotation (small Ekman number $E$) and strong thermal forcing (large Rayleigh number $Ra$). The analysis proceeds from a set of asymptotically reduced equations appropriate for rotationally constrained dynamics; the conjectured range of validity for these equations is $Ra \lesssim E^{-8/5}$. A rigorous bound on heat transport of $Nu \le 20.56Ra^3E^4$ is derived in the limit of infinite Prandtl number using the background method. We demonstrate that the exponent in this bound cannot be improved on using a piece-wise monotonic background temperature profile like the one used here. This is true for finite Prandtl numbers as well, i.e. $Nu \lesssim Ra^3$ is the best upper bound for this particular setup of the background method. The feature that obstructs the availability of a better bound in this case is the appearance of small-scale thermal plumes emanating from (or entering) the thermal boundary layer.

preprint2013arXiv

Stochastic superparameterization in quasigeostrophic turbulence

In this article we expand and develop the authors' recent proposed methodology for efficient stochastic superparameterization (SP) algorithms for geophysical turbulence. Geophysical turbulence is characterized by significant intermittent cascades of energy from the unresolved to the resolved scales resulting in complex patterns of waves, jets, and vortices. Conventional SP simulates large scale dynamics on a coarse grid in a physical domain, and couples these dynamics to high-resolution simulations on periodic domains embedded in the coarse grid. Stochastic SP replaces the nonlinear, deterministic eddy equations on periodic embedded domains by quasilinear stochastic approximations on formally infinite embedded domains. The result is a seamless algorithm which never uses a small scale grid and is far cheaper than conventional SP, but with significant success in difficult test problems. Various design choices in the algorithm are investigated in detail here, including decoupling the timescale of evolution on the embedded domains from the length of the time step used on the coarse grid, and sensitivity to certain assumed properties of the eddies (e.g. the shape of the assumed eddy energy spectrum). The different design choices are compared and contrasted on a stringent test suite for quasigeostrophic turbulence involving two-layer dynamics on a beta-plane forced by an imposed background shear. The algorithms developed are expected to be particularly useful in providing accurate and efficient stochastic parameterizations for use in ensemble-based state estimation and prediction.