Abstract. High-Frequency (HF) radars measure the ocean currents at various spatial and temporal scales. These include tidal currents, wind-driven circulation, density-driven circulation and Stokes drift. Sequential assimilation methods updating the model state have been proven successful to correct the density-driven currents by assimilation of observations such as sea surface height, sea surface temperature and in-situ profiles. However, the situation is different for tides in coastal models since these are not generated within the domain, but are rather propagated inside the domain through the boundary conditions. For improving the modeled tidal variability it is therefore not sufficient to update the model state via data assimilation without updating the boundary conditions. The optimization of boundary conditions to match observations inside the domain is traditionally achieved through variational assimilation methods. In this work we present an ensemble smoother to improve the tidal boundary values so that the model represents more closely the observed currents. To create an ensemble of dynamically realistic boundary conditions, a cost function is formulated which is directly related to the probability of each perturbation. This cost function ensures that the perturbations are spatially smooth and that the structure of the perturbations satisfies approximately the harmonic linearized shallow water equations. Based on those perturbations an ensemble simulation is carried out using the full three-dimension General Estuarine Ocean Model (GETM). Optimized boundary values are obtained using all observations within the assimilation period using the covariances of the ensemble simulation.
An ensemble smoother scheme is presented to assimilate high-frequency (HF) radar surface currents to improve tidal boundary conditions and wind forcings of a circulation model of the German Bight. To create an ensemble of dynamically realistic tidal boundary conditions, a cost function is formulated which is directly related to the probability of each perturbation. This cost function ensures that the perturbations are spatially smooth and that the structure of the perturbations satisfies approximately the harmonic linearized shallow water equations. Based on those perturbations an ensemble simulation is carried out using the full three-dimensional General Estuarine Ocean Model (GETM). Optimized boundary values are obtained using all observations within the assimilation period using the covariances of the ensemble simulation. The approach acts like a smoother scheme since past and future observations are taken into account. The final analysis is obtained by rerunning the model using the optimal perturbation of the boundary conditions. The analyzed model solution satisfies thus the model equations exactly and does not suffer from spurious adjustments often observed with sequential assimilation schemes. Model results are also compared to independent tide gage data. The assimilation also reduces the model error compared to those sea level observations. The same scheme is also used to correct surface winds. Surface winds are crucial for accurately modeling the marine circulation in coastal waters. The method is validated directly by comparing the analyzed wind speed to in situ measurements and indirectly by assessing the impact of the corrected winds on sea surface temperature (SST) relative to satellite SST.
Abstract. High-Frequency (HF) radars measure the ocean surface currents at various spatial and temporal scales. These include tidal currents, wind-driven circulation, density-driven circulation and Stokes drift. Sequential assimilation methods updating the model state have been proven successful to correct the density-driven currents by assimilation of observations such as sea surface height, sea surface temperature and in-situ profiles. However, the situation is different for tides in coastal models since these are not generated within the domain, but are rather propagated inside the domain through the boundary conditions. For improving the modeled tidal variability it is therefore not sufficient to update the model state via data assimilation without updating the boundary conditions. The optimization of boundary conditions to match observations inside the domain is traditionally achieved through variational assimilation methods. In this work we present an ensemble smoother to improve the tidal boundary values so that the model represents more closely the observed currents. To create an ensemble of dynamically realistic boundary conditions, a cost function is formulated which is directly related to the probability of each boundary condition perturbation. This cost function ensures that the boundary condition perturbations are spatially smooth and that the structure of the perturbations satisfies approximately the harmonic linearized shallow water equations. Based on those perturbations an ensemble simulation is carried out using the full three-dimensional General Estuarine Ocean Model (GETM). Optimized boundary values are obtained by assimilating all observations using the covariances of the ensemble simulation.
'%3e%3cpath fill='%23012169' d='M0 0v30h60V0z'/%3e%3cpath stroke='%23fff' stroke-width='6' d='m0 0 60 30m0-30L0 30'/%3e%3cpath stroke='%23C8102E' stroke-width='4' d='m0 0 60 30m0-30L0 30' clip-path='url(%23b)'/%3e%3cpath stroke='%23fff' stroke-width='10' d='M30 0v30M0 15h60'/%3e%3cpath stroke='%23C8102E' stroke-width='6' d='M30 0v30M0 15h60'/%3e%3c/g%3e%3c/svg%3e)
