Intercomparison of the Charnock and COARE bulk wind stress formulations for coastal ocean modelling

Abstract. The accurate parameterisation of momentum and heat transfer across the air–sea interface is vital for realistic simulation of the atmosphere–ocean system. In most modelling applications accurate representation of the wind stress is required to numerically reproduce surge, coastal ocean circulation, surface waves, turbulence and mixing. Different formulations can be implemented and impact the accuracy of the instantaneous and long-term residual circulation, the surface mixed layer, and the generation of wave-surge conditions. This, in turn, affects predictions of storm impact, sediment pathways, and coastal resilience to climate change. The specific numerical formulation needs careful selection to ensure the accuracy of the simulation. Two wind stress parameterisations widely used in the ocean circulation and the storm surge communities respectively are studied with focus on an application to the NW region of the UK. Model–observation validation is performed at two nearshore and one estuarine ADCP (acoustic Doppler current profiler) stations in Liverpool Bay, a hypertidal region of freshwater influence (ROFI) with vast intertidal areas. The period of study covers both calm and extreme conditions to test the robustness of the 10 m wind stress component of the Coupled Ocean–Atmosphere Response Experiment (COARE) bulk formulae and the standard Charnock relation. In this coastal application a realistic barotropic–baroclinic simulation of the circulation and surge elevation is set-up, demonstrating greater accuracy occurs when using the Charnock relation, with a constant Charnock coefficient of 0.0185, for surface wind stress during this one month period.