How does 4DVar data assimilation affect the vertical representation of mesoscale eddies? A case study with observing system simulation experiments (OSSEs) using ROMS v3.9

Accurate estimates and forecasts of ocean eddies in key regions such as western boundary currents are important for weather and climate, biology, navigation, and search and rescue. The dynamic nature of mesoscale eddies requires data assimilation to produce accurate eddy timings and locations in ocean model simulations. However, data assimilating models are rarely assessed below the surface due to a paucity of observations; hence it is not clear how data assimilation impacts the subsurface eddy structure. Here, we use a suite of observing system simulation experiments to show how the subsurface representation of eddies is changed within data assimilating simulations even when assimilating nearby observations. We examine in detail two possible manifestations of how the data assimilation process impacts three-dimensional eddy structure, namely, by producing overly active baroclinic instability and through inaccurate vertical mode structure. Therefore, in DA simulations, subsurface temperature structures can be too deep and too warm, particularly in dynamic eddy features. Our analyses demonstrate the need for further basic research in ocean data assimilation methodologies to improve the representation of the subsurface ocean structure.