Generation of state‐dependent ensemble perturbations based on time‐varying seawater density for GloSea5 initialization

In this study, we developed a flow‐dependent oceanic initialization system for initializing the oceanic temperature and salinity in the Global Seasonal forecast system version 5 (GloSea5). Our algorithm overcomes the limitation of stationary perturbations for Ensemble Optimal Interpolation (EnOI) by spreading observed information along isopycnal lines to create three‐dimensional snapshot density states. The proposed algorithm, which we call state‐dependent ensemble‐based EnOI (SD‐EnOI), takes into account changes in the background error covariance over time without relying on ensemble model simulations. To evaluate the quality of the oceanic initial conditions (ICs) produced by SD‐EnOI, we compared them with those generated by the Global Ocean Data Assimilation and Prediction System version 1 (GODAPS1) operated by the Korea Meteorological Agency (KMA) throughout January 2017 to December 2017. Our findings show that the thermal construction of the SD‐EnOI ICs is more realistic than that of GODAPS1, particularly in the tropical Pacific region. The strong warm bias in sea surface temperature (SST) and the shallow mixed‐layer depth bias observed in the GODAPS1 ICs are not shown in SD‐EnOI. Due to the more realistic oceanic thermal structure present in the SD‐EnOI ICs, their use in retrospective forecast experiments resulted in a systematic reduction in climatological SST drift in the central‐eastern Pacific for forecasts up to four lead months compared to using GODAPS1 ICs. This demonstrates the significant impact of the initialization process on the quality of dynamical seasonal forecasts.