Source author record

Nathan Glatt-Holtz

Nathan Glatt-Holtz 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

15works
8topics
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

15 published item(s)

preprint2026arXiv

Bayesian Semi-Blind Deconvolution at Scale

Blind image deconvolution refers to the problem of simultaneously estimating the blur kernel and the true image from a set of observations when both the blur kernel and the true image are unknown. Sometimes, additional image and/or blur information is available and the term semi-blind deconvolution (SBD) is used. We consider a recently introduced Bayesian conjugate hierarchical model for SBD, formulated on an extended cyclic lattice to allow a computationally scalable Gibbs sampler. In this article, we extend this model to the general SBD problem, rewrite the previously proposed Gibbs sampler so that operations are performed in the Fourier domain whenever possible, and introduce a new marginal Hamiltonian Monte Carlo (HMC) blur update, obtained by analytically integrating the blur-image joint conditional over the image. The cyclic formulation combined with non-trivial linear algebra manipulations allows a Fourier-based, scalable HMC update, otherwise complicated by the rigid constraints of the SBD problem. Having determined the padding size in the cyclic embedding through a numerical experiment, we compare the mixing and exploration behaviour of the Gibbs and HMC blur updates on simulated data and on a real geophysical seismic imaging problem where we invert a grid with $300\times50$ nodes, corresponding to a posterior with approximately $80,000$ parameters.

preprint2020arXiv

The Generalized Langevin Equation with a power-law memory in a nonlinear potential well

The generalized Langevin equation (GLE) is a stochastic integro-differential equation that has been used to describe the velocity of microparticles in viscoelastic fluids. In this work, we consider the large-time asymptotic properties of a Markovian approximation to the GLE in the presence of a wide class of external potential wells. The qualitative behavior of the GLE is largely determined by its memory kernel $K$, which summarizes the delayed response of the fluid medium on the particles past movement. When $K$ can be expressed as a finite sum of exponentials, it has been shown that long-term time-averaged properties of the position and velocity do not depend on $K$ at all. In certain applications, however, it is important to consider the GLE with a power law memory kernel. Using the fact that infinite sums of exponentials can have power law tails, we study the infinite-dimensional version of the Markovian GLE in a potential well. In the case where the memory kernel $K$ is integrable (i.e. in the asymptotically diffusive regime), we are able to extend previous results and show that there is a unique stationary distribution for the GLE system and that the long-term statistics of the position and velocity do not depend on $K$. However, when $K$ is not integrable (i.e. in the asymptotically subdiffusive regime), we are able to show the existence of an invariant probability measure but uniqueness remains an open question. In particular, the method of asymptotic coupling used in the integrable case to show uniqueness does not apply when $K$ fails to be integrable.

preprint2016arXiv

Asymptotic Analysis for Randomly Forced MHD

We consider the three-dimensional magnetohydrodynamics (MHD) equations in the presence of a spatially degenerate stochastic forcing as a model for magnetostrophic turbulence in the Earth's fluid core. We examine the multi-parameter singular limit of vanishing Rossby number $ε$ and magnetic Reynold's number $δ$, and establish that: (i) the limiting stochastically driven active scalar equation (with $ε=δ=0$) possesses a unique ergodic invariant measure, and (ii) any suitable sequence of statistically invariant states of the full MHD system converge weakly, as $ε,δ\rightarrow 0$, to the unique invariant measure of the limit equation. This latter convergence result does not require any conditions on the relative rates at which $\varepsilon, δ$ decay. Our analysis of the limit equation relies on a recently developed theory of hypo-ellipticity for infinite-dimensional stochastic dynamical systems. We carry out a detailed study of the interactions between the nonlinear and stochastic terms to demonstrate that a Hörmander bracket condition is satisfied, which yields a contraction property for the limit equation in a suitable Wasserstein metric. This contraction property reduces the convergence of invariant states in the multi-parameter limit to the convergence of solutions at finite times. However, in view of the phase space mismatch between the small parameter system and the limit equation, and due to the multi-parameter nature of the problem, further analysis is required to establish the singular limit. In particular, we develop methods to lift the contraction for the limit equation to the extended phase space, including the velocity and magnetic fields. Moreover, for the convergence of solutions at finite times we make use of a probabilistic modification of the Grönwall inequality, relying on a delicate stopping time argument.

preprint2016arXiv

Invariant measures for passive scalars in the small noise inviscid limit

We consider a class of invariant measures for a passive scalar $f$ driven by an incompressible velocity field $\boldsymbol{u}$, on a $d$-dimensional periodic domain, satisfying $$ \partial_t f + \boldsymbol{u} \cdot \nabla f = 0, \qquad f(0)=f_0. $$ The measures are obtained as limits of stochastic viscous perturbations. We prove that the span of the $H^1$ eigenfunctions of the operator $\boldsymbol{u}\cdot\nabla$ contains the support of these measures. We also analyze several explicit examples: when $\boldsymbol{u}$ is a shear flow or a relaxation enhancing flow (a generalization of weakly mixing), we can characterize the limiting measure uniquely and compute its covariance structure. We also consider the case of two-dimensional cellular flows, for which further regularity properties of the functions in the support of the measure can be deduced. The main results are proved with the use of spectral theory results, in particular the RAGE theorem, which are used to characterize large classes of orbits of the inviscid problem that are growing in $H^1$.

preprint2014arXiv

Inviscid Limits for a Stochastically Forced Shell Model of Turbulent Flow

We establish the anomalous mean dissipation rate of energy in the inviscid limit for a stochastic shell model of turbulent fluid flow. The proof relies on viscosity independent bounds for stationary solutions and on establishing ergodic and mixing properties for the viscous model. The shell model is subject to a degenerate stochastic forcing in the sense that noise acts directly only through one wavenumber. We show that it is hypo-elliptic (in the sense of Hormander) and use this property to prove a gradient bound on the Markov semigroup.

preprint2014arXiv

Time Discrete Approximation of Weak Solutions for Stochastic Equations of Geophysical Fluid Dynamics and Applications

As a first step towards the numerical analysis of the stochastic primitive equations of the atmosphere and oceans, we study their time discretization by an implicit Euler scheme. From deterministic viewpoint the 3D Primitive Equations are studied with physically realistic boundary conditions. From probabilistic viewpoint we consider a wide class of nonlinear, state dependent, white noise forcings. The proof of convergence of the Euler scheme covers the equations for the oceans, atmosphere, coupled oceanic-atmospheric system and other geophysical equations. We obtain the existence of solutions weak in PDE and probabilistic sense, a result which is new by itself to the best of our knowledge.

preprint2013arXiv

Ergodic and Mixing Properties of the Boussinesq Equations with a Degenerate Random Forcing

We establish the existence, uniqueness and attraction properties of an ergodic invariant measure for the Boussinesq Equations in the presence of a degenerate stochastic forcing acting only in the temperature equation and only at the largest spatial scales. The central challenge is to establish time asymptotic smoothing properties of the Markovian dynamics corresponding to this system. Towards this aim we encounter a Lie bracket structure in the associated vector fields with a complicated dependence on solutions. This leads us to develop a novel Hörmander-type condition for infinite-dimensional systems. Demonstrating the sufficiency of this condition requires new techniques for the spectral analysis of the Malliavin covariance matrix.

preprint2013arXiv

Existence and Regularity of Invariant Measures for the Three Dimensional Stochastic Primitive Equations

We establish the continuity of the Markovian semigroup associated with strong solutions of the stochastic 3D Primitive Equations, and prove the existence of an invariant measure. The proof is based on new moment bounds for strong solutions. The invariant measure is supported on strong solutions, but is furthermore shown to have higher regularity properties.

preprint2013arXiv

Martingale and Pathwise Solutions to the Stochastic Zakharov-Kuznetsov Equation with Multiplicative Noise

We study in this article the stochastic Zakharov-Kuznetsov equation driven by a multiplicative noise. We establish, in space dimensions two and three the global existence of martingale solutions, and in space dimension two the global pathwise uniqueness and the existence of pathwise solutions. New methods are employed in the passage to the limit on a special type of boundary conditions and in the verification of the pathwise uniqueness of martingale solutions with a lack of regularity, where both difficulties arise due to the partly hyperbolic feature of the model.

preprint2013arXiv

On Inviscid Limits for the Stochastic Navier-Stokes Equations and Related Models

We study inviscid limits of invariant measures for the 2D Stochastic Navier-Stokes equations. As shown in \cite{Kuksin2004} the noise scaling $\sqrt{ν}$ is the only one which leads to non-trivial limiting measures, which are invariant for the 2D Euler equations. We show that any limiting measure $μ_{0}$ is in fact supported on bounded vorticities. Relationships of $μ_{0}$ to the long term dynamics of Euler in the $L^{\infty}$ with the weak$^{*}$ topology are discussed. In view of the Batchelor-Krainchnan 2D turbulence theory, we also consider inviscid limits for the weakly damped stochastic Navier-Stokes equation. In this setting we show that only an order zero noise (i.e. the noise scaling $ν^0$) leads to a nontrivial limiting measure in the inviscid limit.

preprint2012arXiv

Global Existence and Regularity for the 3D Stochastic Primitive Equations of the Ocean and Atmosphere with Multiplicative White Noise

The Primitive Equations are a basic model in the study of large scale Oceanic and Atmospheric dynamics. These systems form the analytical core of the most advanced General Circulation Models. For this reason and due to their challenging nonlinear and anisotropic structure the Primitive Equations have recently received considerable attention from the mathematical community. In view of the complex multi-scale nature of the earth's climate system, many uncertainties appear that should be accounted for in the basic dynamical models of atmospheric and oceanic processes. In the climate community stochastic methods have come into extensive use in this connection. For this reason there has appeared a need to further develop the foundations of nonlinear stochastic partial differential equations in connection with the Primitive Equations and more generally. In this work we study a stochastic version of the Primitive Equations. We establish the global existence of strong, pathwise solutions for these equations in dimension 3 for the case of a nonlinear multiplicative noise. The proof makes use of anisotropic estimates, $L^{p}_{t}L^{q}_{x}$ estimates on the pressure and stopping time arguments.

preprint2011arXiv

Parameter Estimation for the Stochastically Perturbed Navier-Stokes Equations

We consider a parameter estimation problem to determine the viscosity $ν$ of a stochastically perturbed 2D Navier-Stokes system. We derive several different classes of estimators based on the first $N$ Fourier modes of a single sample path observed on a finite time interval. We study the consistency and asymptotic normality of these estimators. Our analysis treats strong, pathwise solutions for both the periodic and bounded domain cases in the presence of an additive white (in time) noise.

preprint2010arXiv

Pathwise Solutions of the 2D Stochastic Primitive Equations

In this work we consider a stochastic version of the Primitive Equations (PEs) of the ocean and the atmosphere and establish the existence and uniqueness of pathwise, strong solutions. The analysis employs novel techniques in contrast to previous works in order to handle a general class of nonlinear noise structures and to allow for physically relevant boundary conditions. The proof relies on Cauchy estimates, stopping time arguments and anisotropic estimates.