Source author record

Andrew J. Majda

Andrew J. Majda 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.data-an physics.flu-dyn Machine Learning math.NA math.PR math.ST nlin.CD Statistics Theory

Catalog footprint

What is connected

6works

9topics

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.

preprint2020arXiv

Anomalous waves triggered by abrupt depth changes: laboratory experiments and truncated KdV statistical mechanics

Recent laboratory experiments of Bolles et al. (2019) demonstrate that an abrupt change in bottom topography can trigger anomalous statistics in randomized surface waves. Motivated by these observations, Majda et al. (2019) developed a theoretical framework, based on deterministic and statistical analysis of the truncated Korteweg-de Vries (TKdV) system, that successfully captures key qualitative features of the experiments, including the robust emergence of anomalous statistics and heightened skewness in the outgoing wavefield. Here, we extend these parallel experimental and modeling efforts with several new findings that have resulted from a synergetic interaction between the two. By precisely relating model parameters to physical ones, we calibrate the model inverse temperature to the specific conditions present in the experiments, thereby permitting a quantitative comparison. We find theoretically predicted distributions of surface displacement to match the experimental measurements with surprising detail. Prompted by the presence of surface slope in the TKdV Hamiltonian, we present new experimental measurements on surface slope statistics and compare them to model predictions. Analysis of some deterministic trajectories of TKdV elucidates the experimental length and time scales required for the statistical transition to a skewed state. Finally, the theory predicts a peculiar relationship between the outgoing displacement skewness and the change in slope variance, specifically how their ratio depends on the wave amplitude and depth ratio. New experimental measurements confirm this prediction in spectacular fashion.

preprint2016arXiv

Rigorous accuracy and robustness analysis for two-scale reduced random Kalman filters in high dimensions

Contemporary data assimilation often involves millions of prediction variables. The classical Kalman filter is no longer computationally feasible in such a high dimensional context. This problem can often be resolved by exploiting the underlying multiscale structure, applying the full Kalman filtering procedures only to the large scale vari- ables, and estimating the small scale variables with proper statistical strategies, including multiplicative inflation, representation model error in the observations, and crude localization. The resulting two-scale reduced filters can have close to optimal numerical filtering skill based on previous numerical evidence. Yet, no rigorous explanation exists for this success, because these modifications create unavoidable bias and model error. This paper contributes to this issue by establishing a new error analysis framework for two different reduced random Kalman filters, valid independent of the large dimension. The first part of our results examines the fidelity of the covariance estimators, which is essential for accurate uncertainty quantification. In a simplified setting, this is demonstrated by showing the true error covariance is dominated by its estimators. In general settings, the Mahalanobis error and its intrinsic dissipation can indicate covariance fidelity. The second part develops upper bounds for the covariance estimators by comparing with proper Kalman filters. Combining both results, the classical tools for Kalman filters can be used as a-priori performance criteria for the reduced filters. In applications, these criteria guarantee the reduced filters are robust, and accurate for small noise systems. They also shed light on how to tune the reduced filters for stochastic turbulence.

preprint2015arXiv

Data-driven prediction strategies for low-frequency patterns of North Pacific climate variability

The North Pacific exhibits patterns of low-frequency variability on the intra-annual to decadal time scales, which manifest themselves in both model data and the observational record, and prediction of such low-frequency modes of variability is of great interest to the community. While parametric models, such as stationary and non-stationary autoregressive models, possibly including external factors, may perform well in a data-fitting setting, they may perform poorly in a prediction setting. Ensemble analog forecasting, which relies on the historical record to provide estimates of the future based on past trajectories of those states similar to the initial state of interest, provides a promising, nonparametric approach to forecasting that makes no assumptions on the underlying dynamics or its statistics. We apply such forecasting to low-frequency modes of variability for the North Pacific sea surface temperature and sea ice concentration fields extracted through Nonlinear Laplacian Spectral Analysis. We find such methods may outperform parametric methods and simple persistence with increased predictive skill.

preprint2015arXiv

Multiscale Asymptotics for the Skeleton of the Madden-Julian Oscillation and Tropical-Extratropical Interactions

A new model is derived and analyzed for tropical-extratropical interactions involving the Madden-Julian oscillation (MJO). The model combines (i) the tropical dynamics of the MJO and equatorial baroclinic waves and (ii) the dynamics of barotropic Rossby waves with significant extratropical structure, and the combined system has a conserved energy. The method of multiscale asymptotics is applied to systematically derive a system of ordinary differential equations (ODEs) for three-wave resonant interactions. Two novel features are (i) a degenerate auxiliary problem with overdetermined equations due to a compatibility condition (meridional geostrophic balance) and (ii) cubic self-interaction terms that are not typically found in three-wave resonance ODEs. Several examples illustrate applications to MJO initiation and termination, including cases of (i) the MJO, equatorial baroclinic Rossby waves, and barotropic Rossby waves interacting, and (ii) the MJO, baroclinic Kelvin waves, and barotropic Rossby waves interacting. Resonance with the Kelvin wave is not possible here if only dry variables are considered, but it occurs in the moist model here through interactions with water vapor and convective activity.

preprint2012arXiv

Nonlinear Laplacian spectral analysis: Capturing intermittent and low-frequency spatiotemporal patterns in high-dimensional data

We present a technique for spatiotemporal data analysis called nonlinear Laplacian spectral analysis (NLSA), which generalizes singular spectrum analysis (SSA) to take into account the nonlinear manifold structure of complex data sets. The key principle underlying NLSA is that the functions used to represent temporal patterns should exhibit a degree of smoothness on the nonlinear data manifold M; a constraint absent from classical SSA. NLSA enforces such a notion of smoothness by requiring that temporal patterns belong in low-dimensional Hilbert spaces V_l spanned by the leading l Laplace-Beltrami eigenfunctions on M. These eigenfunctions can be evaluated efficiently in high ambient-space dimensions using sparse graph-theoretic algorithms. Moreover, they provide orthonormal bases to expand a family of linear maps, whose singular value decomposition leads to sets of spatiotemporal patterns at progressively finer resolution on the data manifold. The Riemannian measure of M and an adaptive graph kernel width enhances the capability of NLSA to detect important nonlinear processes, including intermittency and rare events. The minimum dimension of V_l required to capture these features while avoiding overfitting is estimated here using spectral entropy criteria.

preprint2011arXiv

Two-dimensional Moist Stratified Turbulence and the Emergence of Vertically Sheared Horizontal Flows

Moist stratified turbulence is studied in a two-dimensional Boussinesq system influenced by condensation and evaporation. The problem is set in a periodic domain and employs simple evaporation and condensation schemes, wherein both the processes push parcels towards saturation. Numerical simulations demonstrate the emergence of a moist turbulent state consisting of ordered structures with a clear power-law type spectral scaling from initially spatially uncorrelated conditions. An asymptotic analysis in the limit of rapid condensation and strong stratification shows that, for initial conditions with enough water substance to saturate the domain, the equations support a straightforward state of moist balance characterized by a hydrostatic, saturated, vertically sheared horizontal flow (VSHF). For such initial conditions, by means of long time numerical simulations, the emergence of moist balance is verified. Specifically, starting from uncorrelated data, subsequent to the development of a moist turbulent state, the system experiences a rather abrupt transition to a regime which is close to saturation and dominated by a strong VSHF. On the other hand, initial conditions which do not have enough water substance to saturate the domain, do not attain moist balance. Rather, the system remains in a turbulent state and oscillates about moist balance. Even though balance is not achieved with these general initial conditions, the time scale of oscillation about moist balance is much larger than the imposed time scale of condensation and evaporation, thus indicating a distinct dominant slow component in the moist stratified two-dimensional turbulent system.