Review of "An ensemble Kalman filter system with the Stony Brook Parallel Ocean Model v1.0" by Shun Ohishi et al.

Abstract

This study develops an ensemble Kalman filter (EnKF)-based regional ocean data assimilation system, in which the local ensemble transform Kalman filter (LETKF) is implemented with the Stony Brook Parallel Ocean Model (sbPOM) version 1.0 to assimilate satellite and in-situ observations at a daily frequency. A series of sensitivity experiments are performed with various settings of the incremental analysis update (IAU) and covariance inflation methods, for which the relaxation-to-prior perturbations and spread (RTPP and RTPS, respectively) and multiplicative inflation (MULT) are considered. We evaluate the geostrophic balance and the analysis accuracy compared with the control experiment in which the IAU and covariance inflation are not applied. The results show that the IAU improves the geostrophic balance, degrades the accuracy, and reduces the ensemble spread, and that the RTPP and RTPS have the opposite effect. The experiment using the combination of the IAU and RTPP results in significant improvement for both balance and accuracy when the RTPP parameter is 0.8–0.9. The combination of the IAU and RTPS improves the balance when the RTPS parameter is ≤ 0.8 and increases the accuracy for the parameter values between 1.0 and 1.1, but the balance and accuracy are not improved significantly at the same time. The experiments with MULT do not demonstrate sufficient skill in maintaining the balance and reproducing the surface flow field regardless of whether the IAU is applied or not. Therefore, the combination of the IAU and RTPP with the parameter of 0.8–0.9 is found to be the best setting for the EnKF-based ocean data assimilation system.