Insights from very‐large‐ensemble data assimilation experiments with a high‐resolution general circulation model of the Red Sea

• Published in: Quarterly journal of the Royal Meteorological Society 2024-10, Vol.150 (764), p.4235-4251
• Main Authors: Sanikommu, Sivareddy, Raboudi, Naila, El Gharamti, Mohamad, Zhan, Peng, Hadri, Bilel, Hoteit, Ibrahim
• Format: Article
• Language: English
• Subjects: Atmospheric forcing; Atmospheric models; Data assimilation; Data collection; ensembles; general circulation model experiments; General circulation models; geophysical sphere; Kalman filters; Localization; Medium-range forecasting; ocean; Ocean models; Oceans; Parameterization; Physics; regional and mesoscale modeling; Satellite observation; tools and methods; Weather forecasting
• ISSN: 0035-9009; 1477-870X
• DOI: 10.1002/qj.4813

Ensemble Kalman Filters (EnKFs), which assimilate observations based on statistics derived from an ensemble of samples of ocean states, have become the norm for ocean data assimilation (DA) and forecasting. These schemes are commonly implemented with inflation and localization techniques to increase their ensemble spread and to filter out spurious long‐range correlations resulting from the limited‐size ensembles imposed by computational burden constraints. Such ad‐hoc methods were found to be not necessary in ensemble DA experiments with simplified ocean/atmospheric models and large ensembles. Here, we conduct a series of one‐year‐long ensemble experiments with a fully realistic EnKF‐DA system in the Red Sea using tens ‐to thousands of ensemble members. The system assimilates satellite and in‐situ observations and accounts for model uncertainties by integrating a 4‐km‐resolution ocean model with European Center for Medium Range Weather Forecast (ECMWF) atmospheric ensemble fields, perturbed internal physics and initial conditions for forecasting. OceanOur results indicate that accounting for model uncertainties is more beneficial than simply increasing the ensemble size, with the improvements due to large ensembles leveling off at about 250 members. Besides, and in contrast to what is commonly observed with simplified models, the investigated ensemble DA system still required localization even when implemented with thousands of members. These findings are explained by: (i) amplified spurious long‐range correlations produced by the low‐rank nature of the ECMWF atmospheric forcing ensemble; and (ii) non‐Gaussianity generated by the perturbed internal physical parameterization schemes. Large‐ensemble forcing fields and non‐Gaussian DA methods might be needed to get full benefits from large ensembles in ocean DA. Spatial (a) and temporal (b) coverage of different observations. Spatial coverage for sea surface temperature (SST) (red) and sea surface height (SSH) (gray) is shown on a day and during the month of June, 2011, respectively. For temperature (blue) and salinity (green) profiles, the spatial coverage displayed corresponds to the whole year 2011. Panel b plots the number of temperature (blue solid line) and salinity (green solid line) profiles available during each month of the year 2011.