Research

Methods for combining models and observations, and the applications that test them.

Data assimilation combines imperfect models with incomplete, noisy observations to estimate the state of the atmosphere, ocean, and sea ice, and it starts every weather forecast. Our group develops ensemble, particle filter, and machine learning methods for data assimilation, especially where the usual assumptions of linear dynamics and Gaussian errors break down.

A twisted blue surface in three dimensions on which all prior states lie; red ensemble Kalman filter members sit well off the surface, while green local particle filter members stay close to it
Winds at three points around a Rankine vortex with an uncertain position. All prior states lie on a twisted surface; after assimilating radar winds, ensemble Kalman filter members (red) fall well off it, while local particle filter members (green) stay close (Poterjoy 2022a).

New data assimilation methodology

The local particle filter (LPF), introduced by Poterjoy (2016) and developed further by our group, drops the Gaussian assumptions of operational data assimilation while remaining practical for high-dimensional models (Poterjoy et al. 2019; Poterjoy 2022b). Where the states a model allows lie on a compact, curved manifold, ensemble Kalman filters pull members off it, while the LPF keeps them close (Poterjoy 2022a; figure); the two filters can also be blended adaptively (Kurosawa and Poterjoy 2023). We also estimate observation errors online, which works even with non-Gaussian errors but is sensitive to misspecified background uncertainty (Santer et al. 2025; Walsworth et al. 2023).

Current work Learning likelihood functions from data with kernel and machine learning methods, including approaches suited to small ensembles.

Side-by-side maps of upper-level vorticity over North America and the Atlantic in one ensemble member from each filter: the EnKF member is speckled with small-scale noise, while the LPF member shows smooth, continuous streams of vorticity wrapping around an extratropical cyclone
Top: 500-hPa cyclonic vorticity in the hourly LPF posterior mean, global Unified Forecast System. Bottom: 250-hPa vorticity in one EnKF and one LPF member, regional HAFS experiment, 19 August 2020.

Particle filters for global weather prediction

In weather models, ensemble Kalman filter members carry spurious small-scale corrections that LPF members largely avoid (Poterjoy 2022a; bottom panel).

Current work With NOAA, testing the LPF in the Unified Forecast System behind the Global Forecast System, where it matches the accuracy of the ensemble Kalman filter while estimating uncertainty more reliably, and bringing it into JEDI, the data assimilation framework used by NOAA and NASA.

Color scale for reflectivity, from 5 to 75 dBZ
Composite reflectivity and winds from an LPF member in a WRF-based Warn-on-Forecast System, assimilating radar with a Gaussian likelihood and a 25-dBZ threshold (left) or a likelihood learned from data (right). Red contours enclose winds above 30 m s−1 near the level of strongest rear inflow.

Data assimilation for mesoscale weather prediction

Thunderstorms evolve within minutes and are observed mostly by radar. Particle filter analyses of convection stay more consistent with the model than ensemble Kalman filter analyses (Poterjoy et al. 2017; McCurry et al. 2023). We also diagnose convective ensemble uncertainty (Santer et al. 2023) and helped plan the Maryland Mesonet with observing system simulation experiments (McCurry and Poterjoy 2024).

Current work Assimilating radar reflectivity with likelihoods learned from data, so weak echoes can be used without a threshold (figure).

Simulated satellite imagery from basin-scale HAFS analyses of Hurricanes Laura and Marco, August 2020 (animation by Kenta Kurosawa).

Hurricanes and high-impact weather

We develop ensemble and hybrid data assimilation for NOAA's hurricane models. Continuous cycling in HWRF exposed errors hidden by restarting from global analyses (Poterjoy et al. 2021); in HAFS, estimating radiance bias corrections within the system improved forecasts (Knisely and Poterjoy 2023), and global ensemble perturbations mitigate sampling errors (Kurosawa and Poterjoy 2025).

Current work Basin-wide HAFS experiments reveal cloud-driven model biases and show that the LPF keeps ensemble spread closer to forecast errors than the EnKF.

Left, a map of Arctic sea ice concentration with three marked points; right, histograms of observed sea ice concentration where the analysis is 0.14, 0.56, and 0.86, which are skewed, bounded by 0 and 1, and in one case bimodal
(a) Sea ice concentration (SIC) analysis from NOAA's coupled ocean–ice experiments. (b) Distributions of observed SIC where the analysis SIC is 0.14, 0.56, and 0.86; red markers show example locations. Figure created by Ph.D. student Henry Santer for his research (adapted from Wang et al. 2026).

Sea ice, polar prediction, and observing system design

Sea ice concentration is bounded, and its observation errors depend on the ice state (Wang et al. 2026; figure). We model atmospheric forcing errors for sea ice ensembles (Jung et al. 2026) and use observing system simulation experiments to inform requirements for a future NASA satellite mission (Choi et al. 2025).

Current work Learning sea ice observation errors from data for the LPF in the CICE model and, with NASA's PREFIRE science team, developing methods that assimilate far-infrared radiances directly.

Current support

NASA PREFIRE Science Team

Advancing Polar Sea Ice Prediction With PREFIRE Radiance Observations and Next-Generation Data Assimilation Methods

Co-I (PI: Alek Petty) · 2026–2029
NOAA

A Multi-University Consortium for Advanced Data Assimilation Research and Education (CADRE)

Institutional PI · 2024–2027
NOAA

Non-Gaussian Ensemble Data Assimilation Within a Rapidly Updating Global Prediction System

PI · 2023–2027
NSF

Online Uncertainty Quantification for Novel Atmospheric Measurements

PI · 2022–2026

We also gratefully acknowledge computing support from the NSF NCAR Computational and Information Systems Laboratory.