Simulating the temporal and spatial dynamics of the North Sea using the new model GETM (general estuarine transport model)

Here we present results of a 1-year realistic North Sea simulation from the new model GETM (general estuarine transport model) and assess the capabilities of this model by comparing them to model results from the well-known HAMSOM (Hamburg shelf sea and ocean model) model, in situ data from the North Sea project and satellite-derived sea-surface temperature data. The annual cycle and the spatial variability of stratification and mixing in the North Sea is simulated. It is shown that the new model is successful in reproducing the general temporal and spatial dynamics of the North Sea. The major advantages of GETM for achieving improved results in this simulation are the implementation of general vertical coordinates, of a state-of-the-art turbulence model and of higher-order advection schemes. By exploiting the full capabilities of these features a more realistic simulation could be achieved. We found that the greatest differences in the model results are produced by applying advection schemes of different complexity. Here we are able to demonstrate that better advection schemes lead to stronger horizontal gradients and stronger vertical stratification during summer. When comparing these results to measurements from the North Sea project and to satellite data, we find that these stronger gradients are more realistic. Therefore, we consider it as essential to use such high-order advection schemes if the spatial variability of estuarine or shelf seas like the North Sea is to be resolved adequately. The advanced turbulence closure scheme also contributed to more realistic simulation of the vertical stratification. Finally, general vertical coordinates better resolve the shallow regions, but are also useful for the deeper regions, as they allow a better estimation of sea-surface temperature compared to traditional coordinates.Responsible Editor: Phil Dyke     

Sediment transport and resuspension due to combined motion of wave and current in the northern Adriatic Sea during a Bora event in January 2001: A numerical modelling study

The Adriatic Sea general circulation model coupled to a third generation wave model SWAN and a sediment transport model was implemented in the Adriatic Sea to study the dynamics of the sediment transport and resuspension in the northern Adriatic Sea (NAS) during the Bora event in January 2001. The bottom boundary layer (BBL) was resolved by the coupled model with high vertical resolution, and the mechanism of the wave–current interaction in the BBL was also represented in the model. The study found that, during the Bora event of 13–17 January 2001, large waves with significant wave height 2 m and period of 5 s were generated by strong winds in the northwestern shelf of the Adriatic where the direction of wave propagation was orthogonal to the current. The combined motion of the wave and current in the BBL increased the bottom stress over the western Adriatic shelf, resulting in stronger sediment resuspension there. Combining stronger bottom resuspension and strong upward vertical flux of resuspended sediments due to turbulent mixing, the model predicted that sediment concentration near the Po River was much higher than that predicted by the model run without wave forcing. The study also shows that wave–current interaction in the BBL reduced the western Adriatic Coastal Currents (WACCs) in the shallower north. It is concluded that wave forcing significantly changed the sediment distributions and increased the total horizontal fluxes over the western shelf. These results signified wave effect on sediment flux and distribution in the NAS, and suggested that waves cannot be neglected in the study of dynamics of sediment transport and resuspension in the shallow coastal seas. By including the tidal forcing in the coupled model, we also examined the effect of tides on the sediment transport dynamics in the NAS.     

High-resolution atmospheric forcing for regional oceanic model: the Iroise Sea

This study was aimed at modeling, as realistically as possible, the dynamics and thermodynamics of the Iroise Sea by using the Model for Applications at Regional Scale (MARS), a regional ocean 3D model. The horizontal resolution of the configuration in use is 2 km with 30 vertical levels. The 3D model of the Iroise Sea is embedded in a larger model providing open boundary conditions. As regards the atmospheric forcing, the originality of this study is to force the regional ocean model with the high-resolution (6 km) regional meteorological model, Weather Research and Forecasting (WRF). In addition, as the air surface temperature is highly sensitive to the sea surface temperature (SST), this regional meteorological model is improved by taking into account a regional climatologic SST to compute meteorological parameters. By allowing a better coherence between the SST and the temperature of the atmospheric boundary layer while giving a more realistic representation of heat fluxes exchanged at the air/sea interface, this forcing constitutes a noticeable improvement of the Iroise Sea modeling. The different sensitivity tests discussed here pinpoint the importance of entering, in WRF, SST data of sufficiently high quality before the computation of meteorological forcing when the aim is a study of dynamics and thermodynamics far away from the coast. On the other hand, when the target is the reproduction of coastal small-scale features in Iroise Sea modeling, the resolution of the meteorological forcing and the quality of SST are both paramount. The simulation of reference was carried out throughout the Summer and Autumn of year 2005 to allow comparisons with a campaign of surface current measurements by high-frequency radars conducted at the same period.     

Simulation of the 1953 storm surge in the North Sea

The 1953 North Sea floods, the Big Flood, was one of the worst natural disasters in Europe in modern times and is probably one of the most studied severe coastal floods. Several factors led to the devastating storm surge along the southern North Sea coast in combination of strong and sustained northerly winds, invert barometric effect, high spring tide, and an accumulation of the large surge in the Strait of Dover. However, the storm waves and their roles during the 1953 North Sea storm surge are not well investigated. Therefore, the effect of wave setup due to breaking waves in the storm surge processes is investigated through numerical experiments. A coupled process-based tide-wave-surge model was used to investigate and simulate the storm surge in the North Sea during January 31–February 1, 1953 and validated by comparing with historical water level records at tide gauges and wave observations at light vessels in the North Sea. Meteorological forcing inputs for the period, January 27–February 3, 1953 are reproduced from ERA-20C reanalysis data with a constant correction factor for winds. From the simulation results, it is found that, in addition to the high water due to wind setup, wave setup due to breaking waves nearshore play a role of approximately 10% of the storm surge peaks with approximately 0.2 m. The resulting modeling system can be used extensively for the preparedness of the storm surge and wave of extreme condition, and usual barotropic forecast.     

Synergy of satellite remote sensing and numerical modeling for monitoring of suspended particulate matter

Monitoring and modeling of the distribution of suspended particulate matter (SPM) is an important task, especially in coastal environments. Several SPM models have been developed for the North Sea. However, due to waves in shallow water and strong tidal currents in the southern part of the North Sea, this is still a challenging task. In general there is a lack of measurements to determine initial distributions of SPM in the bottom sediment and essential model parameters, e.g., appropriate exchange coefficients. In many satellite-borne ocean color images of the North Sea a plume is visible, which is caused by the scattering of light at SPM in the upper ocean layer. The intensity and length of the plume depends on the wave and current climate. It is well known that the SPM plume is especially obvious shortly after strong storm events. In this paper a quasi-3-D and a 3-D SPM transport model are presented. Utilizing the synergy of satellite-borne ocean color data with numerical models, the vertical exchange coefficients due to currents and waves are derived. This results in models that for the first time are able to reproduce the temporal and spatial evolution of the plume intensity. The SPM models consist of several modules to compute ocean dynamics, the vertical and horizontal exchange of SPM in the water column, and exchange processes with the seabed such as erosion, sedimentation, and resuspension. In the bottom layer, bioturbation via benthos and diffusion processes is taken into account.Responsible Editor: Jörg-Olaf Wolff     

Modeling studies of the upper ocean response to a tropical cyclone

A coupled ocean and boundary layer flux numerical modeling system is used to study the upper ocean response to surface heat and momentum fluxes associated with a major hurricane, namely, Hurricane Dennis (July 2005) in the Gulf of Mexico. A suite of experiments is run using this modeling system, constructed by coupling a Navy Coastal Ocean Model simulation of the Gulf of Mexico to an atmospheric flux model. The modeling system is forced by wind fields produced from satellite scatterometer and atmospheric model wind data, and by numerical weather prediction air temperature data. The experiments are initialized from a data assimilative hindcast model run and then forced by surface fluxes with no assimilation for the time during which Hurricane Dennis impacted the region. Four experiments are run to aid in the analysis: one is forced by heat and momentum fluxes, one by only momentum fluxes, one by only heat fluxes, and one with no surface forcing. An equation describing the change in the upper ocean hurricane heat potential due to the storm is developed. Analysis of the model results show that surface heat fluxes are primarily responsible for widespread reduction (0.5°–1.5°C) of sea surface temperature over the inner West Florida Shelf 100–300 km away from the storm center. Momentum fluxes are responsible for stronger surface cooling (2°C) near the center of the storm. The upper ocean heat loss near the storm center of more than 200 MJ/m² is primarily due to the vertical flux of thermal energy between the surface layer and deep ocean. Heat loss to the atmosphere during the storm’s passage is approximately 100–150 MJ/m². The upper ocean cooling is enhanced where the preexisting mixed layer is shallow, e.g., within a cyclonic circulation feature, although the heat flux to the atmosphere in these locations is markedly reduced.     

Effects of wave-induced forcing on a circulation model of the North Sea

The effect of wind waves on water level and currents during two storms in the North Sea is investigated using a high-resolution Nucleus for European Modelling of the Ocean (NEMO) model forced with fluxes and fields from a high-resolution wave model. The additional terms accounting for wave-current interaction that are considered in this study are the Stokes-Coriolis force, the sea-state-dependent energy and momentum fluxes. The individual and collective role of these processes is quantified and the results are compared with a control run without wave effects as well as against current and water-level measurements from coastal stations. We find a better agreement with observations when the circulation model is forced by sea-state-dependent fluxes, especially in extreme events. The two extreme events, the storm Christian (25–27 October 2013), and about a month later, the storm Xaver (5–7 December 2013), induce different wave and surge conditions over the North Sea. Including the wave effects in the circulation model for the storm Xaver raises the modelled surge by more than 40 cm compared with the control run in the German Bight area. For the storm Christian, a difference of 20–30 cm in the surge level between the wave-forced and the stand-alone ocean model is found over the whole southern part of the North Sea. Moreover, the modelled vertical velocity profile fits the observations very well when the wave forcing is accounted for. The contribution of wave-induced forcing has been quantified indicating that this represents an important mechanism for improving water-level and current predictions.     

Evaluation of a Finite-Element Sea-Ice Ocean Model (FESOM) set-up to study the interannual to decadal variability in the deep-water formation rates

The characteristics of a global set-up of the Finite-Element Sea-Ice Ocean Model under forcing of the period 1958–2004 are presented. The model set-up is designed to study the variability in the deep-water mass formation areas and was therefore regionally better resolved in the deep-water formation areas in the Labrador Sea, Greenland Sea, Weddell Sea and Ross Sea. The sea-ice model reproduces realistic sea-ice distributions and variabilities in the sea-ice extent of both hemispheres as well as sea-ice transport that compares well with observational data. Based on a comparison between model and ocean weather ship data in the North Atlantic, we observe that the vertical structure is well captured in areas with a high resolution. In our model set-up, we are able to simulate decadal ocean variability including several salinity anomaly events and corresponding fingerprint in the vertical hydrography. The ocean state of the model set-up features pronounced variability in the Atlantic Meridional Overturning Circulation as well as the associated mixed layer depth pattern in the North Atlantic deep-water formation areas.     

Wind-driven diversion of summer river runoff preconditions the Laptev Sea coastal polynya hydrography: Evidence from summer-to-winter hydrographic records of 2007–2009

This paper examines the role of atmospheric forcing in modifying the pathways of riverine water on the Laptev Sea shelf, using summer-to-winter hydrographic surveys from 2007 to 2009. Over the two consecutive winter seasons of 2007–2008 and 2008–2009 in the area of the winter coastal polynya, our data clearly link winter surface salinity fields to the previous summer conditions, with substantially different winter salinity patterns preconditioned by summer atmospheric forcing. In the summer of 2007, dominant along-shore westerly winds in the cyclonic regime force the Lena River runoff to flow eastward. In contrast, in the summer of 2008, dominant along-shore easterly winds over the East Siberian Sea and on-shore northerly winds over the Laptev Sea in the anticyclonic regime lock the riverine water in the vicinity of the Lena Delta. Over the coastal polynya area in the southeastern Laptev Sea these patterns precondition a surface salinity difference of 8–16 psu between the winters of 2008 and 2009. Overall, this indicates a residence time of at least half a year for riverine water on the Laptev Sea shelf. Future climate change associated with an enhanced summer cyclonicity over the eastern Arctic may turn more riverine water eastward along the eastern Siberian coast, resulting in weaker vertical density stratification over the Laptev Sea shelf, with possible impact on the efficiency of vertical mixing and polynya dense water production.     

Extreme storm surge modelling in the North Sea

This study shows that storm surge model performance in the North Sea is mostly unaffected by the application of temporal variations of surface drag due to changes in sea state provided the choice of a suitable constant Charnock parameter in the sea-state-independent case. Including essential meteorological features on smaller scales and minimising interpolation errors by increasing forcing data resolution are shown to be more important for the improvement of model performance particularly at the high tail of the probability distribution. This is found in a modelling study using WAQUA/DCSMv5 by evaluating the influence of a realistic air-sea momentum transfer parameterization and comparing it to the influence of changes in the spatial and temporal resolution of the applied forcing fields in an effort to support the improvement of impact and climate analysis studies. Particular attention is given to the representation of extreme water levels over the past decades based on the example of the Netherlands. For this, WAQUA/DCSMv5 is forced with ERA-Interim reanalysis data. Model results are obtained from a set of different forcing fields, which either (i) include a wave-state-dependent Charnock parameter or (ii) apply a constant Charnock parameter (α _{C h} =?0.032) tuned for young sea states in the North Sea, but differ in their spatial and/or temporal resolution. Increasing forcing field resolution from roughly 79 to 12 km through dynamically downscaling can reduce the modelled low bias, depending on coastal station, by up to 0.25 m for the modelled extreme water levels with a 1-year return period and between 0.1 m and 0.5 m for extreme surge heights.     

Mechanisms controlling the intra-annual mesoscale variability of SST and SPM in the southern North Sea

Thermal and optical remote sensing data were used to investigate the spatial and temporal distribution of sea surface temperature (SST) and of suspended particulate matter (SPM) in the southern North Sea. Monthly SST composites showed pronounced seasonal warming of the southern North Sea and delineated the English coastal and continental coastal waters. The East-Anglia Plume is the dominant feature of the English coastal waters in the winter and autumn SPM composites, and the Rhine region of freshwater influence (ROFI), including the Flemish Banks, is the dominant feature of the continental waters. These mesoscale spatial structures are also influenced by the evolution of fronts, such as the seasonal front separating well-mixed water in the southern Bight, from the seasonally stratified central North Sea waters. A harmonic analysis of the SST and SPM images showed pronounced seasonal variability, as well as spring-neap variations in the level of tidal mixing in the East Anglia Plume, the Rhine ROFI and central North Sea. The harmonic analysis indicates the important role played by the local meteorology and tides in governing the SST and near-surface SPM concentrations in the southern North Sea. In the summer, thermal stratification affects the visibility of SPM to satellite sensors in the waters to the north of the Flamborough and Frisian Fronts. Haline stratification plays an important role in the visibility of SPM in the Rhine ROFI throughout the year. When stratified, both regions typically exhibit low surface SPM values. A numerical model study, together with the harmonic analysis, highlights the importance of tides and waves in controlling the stratification in the southern North Sea and hence the visibility of SPM.     

The Mutual Response of the Tropical Squall Line and the Ocean

—The Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS) is used to investigate the mutual response of a tropical squall line and the ocean. Simulated squall line compares well with the observations, and consists of counterrotating vortices, and has a bow shape bulge toward the leading edge. In addition to these features, which are also shown in the previous numerical simulations, the unique results from the coupled simulation indicate that the air–sea interaction processes within the squall line are important. They affect both the atmosphere and the ocean locally. Simulated upper ocean displays significant response to the squall line with upwelling and baroclinicity. Depth of the ocean mixed layer in the coupled simulation becomes modified due to feedback processes. Ocean temperature acts as a destabilizing factor, and the salinity as a stabilizing factor. Surface turbulent fluxes from the coupled simulation are about 10% less than that of the uncoupled simulation. The SST in the coupled simulation decreases by about 0.21°C. Predicted squall line in the coupled simulation is weaker as compared to the uncoupled simulation. This is reflected in terms of differences in surface fluxes, cloud water, rain water and vertical velocities between the two simulations.     

The vertical mixing induced by winds and tides over the Yellow Sea in summer: a numerical study in 2012

Ocean Dynamics - In the Yellow Sea (YS), besides the energetic tidal forcing, winds also induce strong vertical mixing, especially in the summer season when the thermocline is well developed. The...     

Effects of the air–sea coupling time frequency on the ocean response during Mediterranean intense events

The near-sea surface meteorological conditions associated with the Mediterranean heavy precipitation events constitute, on a short time scale, a strong forcing on the ocean mixed layer. This study addresses the question of the optimal time frequency of the atmospheric forcing to drive an ocean model in order to make it able to capture the fine scale ocean mixed layer response to severe meteorological conditions. The coupling time frequency should allow the ocean model to reproduce the formation of internal low-salty boundary layers due to sudden input of intense precipitation, as well as the cooling and deepening of the ocean mixed layer through large latent heat fluxes and stress under the intense low-level jet associated with these events. In this study, the one-dimensional ocean model is driven by 2.4-km atmospheric simulated fields on a case of Mediterranean heavy precipitation, varying the time resolution of the atmospheric forcing. The results show that using a finer temporal resolution than 1 h for the atmospheric forcing is not necessary, but a coarser temporal resolution (3 or 6 h) modifies the event course and intensity perceived by the ocean. Consequently, when using a too coarse temporal resolution forcing, typically 6 h, the ocean model fails to reproduce the ocean mixed layer fine scale response under the heavy rainfall pulses and the strong wind gusts.     

层结稳定性“豁口”与北太平洋副热带中部模态水形成机制

本文首先指出北太平洋副热带中部模态水（简称中部模态水）的形成具有显著的“局地”特征,其形成海区在（165°E～160°W,38°N～42°N）区间. 海气通量分析表明单纯的外部大气强迫场（太阳短波辐射、净热通量和风应力旋度）不能解释中部模态水形成海区的“局地”性;进一步对上层海洋层结季节变化特征的分析发现秋季（9～10月）在北太平洋中部上层海洋（<75 m）（165°E～160°W,38°N～42°N）区间存在特殊的浮力频率低值区——层结稳定性“豁口”. 该层结稳定性“豁口”作为“预条件（Precondition Mechanism）”机制对中部模态水形成的“局地”特征给出了合理的解释. 在上述研究的基础上,基于一个上层海洋混合层热平衡方程,通过诊断分析揭示该层结稳定性“豁口”是由海表热力强迫、垂向挟卷、Ekman平流和地转平流效应共同导致的,“豁口”东、西边界的确定直接或间接地取决于海表热力强迫、Ekman冷平流和地转暖平流的纬向分布差异.     

Composition and daytime vertical distribution of the ichthyoplankton assemblage in the Central Cantabrian Sea shelf, during summer: An Eulerian study

During summer, wind driven coastal upwelling dominates in the Central Cantabrian Sea (southern Bay of Biscay). Nevertheless, atmospheric forcing is highly variable and wind pulses may cause noticeable and fast hydrographic responses in the shelf region. In this paper, the composition and vertical distribution of the summer ichthyoplankton assemblage during the daytime at a fixed station, located on the Central Cantabrian Sea shelf, are documented. Also, the impact of a short-time scale hydrographic event on the abundance and structure of the larval fish assemblage is examined. Significant small-scale temporal hydrographic variability was observed. Currents showed changes in speed and direction and significant changes in thermocline depth were also observed. A total of 34 taxa of fish larvae were identified. Engraulis encrasicolus eggs and larvae of the shelf-dwelling species Trachurus trachurus, Capros aper and E. encrasicolus dominated the ichthyoplankton assemblage. The distribution of E. encrasicolus eggs and fish larvae was vertically structured. E. encrasicolus egg concentration increased exponentially towards the surface. Fish larvae showed a subsurface peak of concentration and their vertical distribution was not conditioned by thermocline depths. The short term hydrographic event did not affect the vertical distribution of fish larvae but it accounted for significant temporal changes in larval fish assemblage structure and abundance. Results suggest that temperature and light intensity are important factors in the vertical distribution of fish larvae. They also indicate that the temporal monitoring of the larval fish assemblage in this region requires multiple sampling sites.     

Ocean response to strong precipitation events in the Gulf of Lions (northwestern Mediterranean Sea): a sensitivity study

The Mediterranean Sea is a region of intense air–sea interactions, with in particular strong evaporation over sea which drives the thermohaline circulation. The Mediterranean region is also prone to strong precipitation events characterized by low spatial extent, short duration, and high temporal variability. The impacts of intense offshore precipitation over sea, in the Gulf of Lions which is a spot for winter deep convection, are investigated using four sensitivity simulations performed at mesoscale resolution with the eddy-resolving regional ocean model NEMO-MED12. We use various atmospheric fields to force NEMO-MED12, downscaled from reanalyses with the non-hydrostatic mesoscale Weather Research and Forecasting model but differing in space resolutions (20 and 6.7 km) or in time frequencies (daily and three-hourly). This numerical study evidences that immediate, intense, and rapid freshening occurs under strong precipitation events. The strong salinity anomaly induced extends horizontally (≃50 km) as vertically (down to 50 m) and persists several days after strong precipitation events. The change in the space resolution of the atmospheric forcing modifies the precipitating patterns and intensity, as well as the shape and the dynamics of the low-salinity layer formed are changed. With higher forcing frequency, shorter and heavier precipitation falls in the ocean in the center of the Gulf of Lions, and due to a stronger vertical shear and mixing, the low-salinity anomaly propagates deeper.     

Influence of high-resolution wind forcing on hydrodynamic modeling of the Gulf of Lions

The impact of the choice of high-resolution atmospheric forcing on ocean summertime circulation in the Gulf of Lions (GoL; Mediterranean Sea) is evaluated using three different datasets: AROME (2.5 km, 1 h), ALADIN (9.5 km, 3 h), and MM5 (9 km, 3 h). A short-term ocean simulation covering a 3-month summer period was performed on a 400-m configuration of the GoL. The main regional features of both wind and oceanic dynamics were well-reproduced by all three atmospheric models. Yet, at smaller scales and for specific hydrodynamic processes, some differences became apparent. Inertial oscillations and mesoscale variability were accentuated when high-resolution forcing was used. Sensitivity tests suggest a predominant role for spatial rather than temporal resolution of wind. The determinant influence of wind stress curl was evidenced, both in the representation of a mesoscale eddy structure and in the generation of a specific upwelling cell in the north-western part of the gulf.     

Mechanisms controlling the seasonal mixed-layer temperature and salinity of the Indonesian seas

We examine the seasonal mixed-layer temperature (MLT) and salinity (MLS) budgets in the Banda–Arafura Seas region (120–138° E, 8–3° S) using an ECCO ocean-state estimation product. MLT in these seas is relatively high during November–May (austral spring through fall) and relatively low during June–September (austral winter and the period associated with the Asian summer monsoon). Surface heat flux makes the largest contribution to the seasonal MLT tendency, with significant reinforcement by subsurface processes, especially turbulent vertical mixing. Temperature declines (the MLT tendency is negative) in May–August when seasonal insolation is smallest and local winds are strong due to the southeast monsoon, which causes surface heat loss and cooling by vertical processes. In particular, Ekman suction induced by local wind stress curl raises the thermocline in the Arafura Sea, bringing cooler subsurface water closer to the base of the mixed layer where it is subsequently incorporated into the mixed layer through turbulent vertical mixing; this has a cooling effect. The MLT budget also has a small, but non-negligible, semi-annual component since insolation increases and winds weaken during the spring and fall monsoon transitions near the equator. This causes warming via solar heating, reduced surface heat loss, and weakened turbulent mixing compared to austral winter and, to a lesser extent, compared to austral summer. Seasonal MLS is dominated by ocean processes rather than by local freshwater flux. The contributions by horizontal advection and subsurface processes have comparable magnitudes. The results suggest that ocean dynamics play a significant part in determining both seasonal MLT and MLS in the region, such that coupled model studies of the region should use a full ocean model rather than a slab ocean mixed-layer model.     

Intercomparison of oceanic and atmospheric forced and coupled mesoscale simulations. Part I: Surface fluxes

A mesoscale non-hydrostatic atmospheric model has been coupled with a mesoscale oceanic model. The case study is a four-day simulation of a strong storm event observed during the SEMAPHORE experiment over a 500 × 500 km² domain. This domain encompasses a thermohaline front associated with the Azores current. In order to analyze the effect of mesoscale coupling, three simulations are compared: the first one with the atmospheric model forced by realistic sea surface temperature analyses; the second one with the ocean model forced by atmospheric fields, derived from weather forecast re-analyses; the third one with the models being coupled. For these three simulations the surface fluxes were computed with the same bulk parametrization. All three simulations succeed well in representing the main oceanic or atmospheric features observed during the storm. Comparison of surface fields with in situ observations reveals that the winds of the fine mesh atmospheric model are more realistic than those of the weather forecast re-analyses. The low-level winds simulated with the atmospheric model in the forced and coupled simulations are appreciably stronger than the re-analyzed winds. They also generate stronger fluxes. The coupled simulation has the strongest surface heat fluxes: the difference in the net heat budget with the oceanic forced simulation reaches on average 50 Wm⁻² over the simulation period. Sea surface-temperature cooling is too weak in both simulations, but is improved in the coupled run and matches better the cooling observed with drifters. The spatial distributions of sea surface-temperature cooling and surface fluxes are strongly inhomogeneous over the simulation domain. The amplitude of the flux variation is maximum in the coupled run. Moreover the weak correlation between the cooling and heat flux patterns indicates that the surface fluxes are not responsible for the whole cooling and suggests that the response of the ocean mixed layer to the atmosphere is highly non-local and enhanced in the coupled simulation.