# Impact of resolution on the atmosphere-ocean coupling along the Gulf Stream in global high resolution models

E.E. Tsartsali, R.J. Haarsma, P.J. Athanasiadis, Alessio Bellucci, H. De Vries, S. Drijfhout, I.E. De Vries, D. Putrahasan, M.J. Roberts, Emilia Sanchez-Gomez, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

E.E. Tsartsali, R.J. Haarsma, P.J. Athanasiadis, Alessio Bellucci, H. De Vries, et al., Impact of resolution on the atmosphere-ocean coupling along the Gulf Stream in global high resolution models. Climate Dynamics, 2022, 58 (11), pp.3317-3333. 10.1007/s00382-021-06098-9. hal-04740244

HAL Id: hal-04740244

https://cnrs.hal.science/hal-04740244v1

Submitted on 16 Oct 2024

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Impact of resolution on the atmosphere-ocean coupling along the Gulf Stream in global high resolution models

E. E. Tsartsali<fnref n="1" />,<fnref n="2" /> ID · R. J. Haarsma<fnref n="2" /> · P. J. Athanasiadis<fnref n="1" /> · A. Bellucci<fnref n="1" />,<fnref n="3" /> · H. de Vries<fnref n="2" /> · S. Drijfhout<fnref n="2" /> · I. E. de Vries<fnref n="2" />,<fnref n="4" /> · D. Putrahasan<fnref n="5" /> · M. J. Roberts<fnref n="6" /> · E. Sanchez-Gomez<fnref n="7" /> · C. D. Roberts<fnref n="8" />

Received: 19 February 2021 / Accepted: 5 December 2021 / Published online: 19 February 2022
© The Author(s) 2022

## Abstract

We have investigated the horizontal resolution dependence of the ocean–atmosphere coupling along the Gulf Stream, of simulations made by six Global Climate Models according to the HighResMIP protocol, and compared it with reanalysis and remote sensing observations. Two ocean–atmosphere interaction mechanisms are explored in detail: The Vertical Mixing Mechanism (VMM) associated with the intensification of downward momentum transfer, and the Pressure Adjustment Mechanism (PAM) associated with secondary circulations driven by pressure gradients. Both VMM and PAM are found to be active even in the eddy-parameterized models. However, increasing ocean and/or atmosphere resolution leads to enhanced ocean–atmosphere coupling and improved agreement with reanalysis and observations. Our results indicate that while one part of the stronger air–sea coupling is attributable to the refinement of the oceanic component to eddy-permitting, optimal results are obtained only by further increase of the atmosphere resolution too. The use of the eddy-resolving model show weaker or same coupling strength over the eddy-permitting resolution. We conclude that at least eddy-permitting ocean resolution and comparable atmosphere resolution are required for a reliable ocean–atmosphere coupling along the Gulf Stream.

**Keywords** Air–sea interaction · Mesoscale coupling · Gulf Stream · Vertical mixing mechanism (VMM) · Pressure adjustment mechanism (PAM)

## 1 Introduction

There is growing evidence that strong sea surface temperature (SST) gradients, oceanic fronts and eddies can force the atmosphere by affecting the near surface wind. At basin scales there is a negative correlation between SST and surface wind speed. Higher wind speeds are located over cooler SSTs, as the atmosphere is forcing the ocean through evaporative cooling and entrainment of colder thermocline waters to the surface (Cayan 1992; Liu et al. 1994; Alexander et al. 2002; Deser et al. 2010). In contrast, at smaller scales, in regions of large mesoscale oceanic activity, higher wind speeds are associated with warm SST anomalies (Small et al. 2008; O’Neill et al. 2005; Xie 2004). This reversed relation indicates that at mesoscale, the ocean is forcing the atmosphere. The surface wind response to the spatial variability of SST generates divergence and curl of the near surface wind (Chelton et al. 2004). The proposed potential mechanisms, which are responsible for the SST influence on surface winds, are the VMM and the PAM. Additionally, there are recent studies suggesting that atmospheric fronts

<figure><img src="image_3.png" /></figure>

<footnote>E. E. Tsartsali
eirini.tsartsali@cmcc.it</footnote>

<footnote n="1">¹ Fondazione Centro Euro-Mediterraneo sui Cambiamenti Climatici (CMCC), Bologna, Italy</footnote>

<footnote n="2">² Koninklijk Nederlands Meteorologisch Instituut (KNMI), De Bilt, The Netherlands</footnote>

<footnote n="3">³ Present Address: CNR-ISAC, Bologna, Italy</footnote>

<footnote n="4">⁴ Eidgenössische Technische Hochschule Zürich (ETHZ), Zurich, Switzerland</footnote>

<footnote n="5">⁵ Max Planck Gesellschaft zur Foerderung der Wissenschaften E.V. (MPI-M), Hamburg, Germany</footnote>

<footnote n="6">⁶ Met Office, Exeter EX1 3PB, UK</footnote>

<footnote n="7">⁷ CECI, Université de Toulouse, CERFACS/CNRS, Toulouse, France</footnote>

<footnote n="8">⁸ European Centre for Medium Range Weather Forecasting (ECMWF), Reading, UK</footnote>
and storm systems can also affect and contribute to the time-
mean near-surface wind convergence (Parfitt and Seo 2018;
O'Neill et al. 2017; Parfitt and Czaja 2016), but this contri-
bution is beyond the scope of this study.

The VMM describes that over warm SST anomalies,
there is an intensification of the downward momentum
transfer due to the destabilization of the overlying marine
atmospheric boundary layer (MABL), which results in accel-
eration of surface winds (Wallace et al. 1989; Sweet et al.
1981). Acceleration of the wind over warm waters, where
it is blowing parallel to the SST front, results in horizon-
tal wind shear, leading to low-level vorticity. In contrast,
acceleration where the wind is blowing across the SST front,
results in wind stress divergence (convergence) when the
wind flow is from a cold to a warm (warm to cold) area
(Chelton and Xie 2010; Chelton et al. 2004; O'Neill et al.
2003). According to this, the downwind (crosswind) SST
gradients are proportional to the wind stress divergence
(curl) (Chelton et al. 2001, 2004; O'Neill et al. 2003).

On the other hand, PAM describes the response of sec-
ondary circulations in the MABL which are driven by pres-
sure gradients due to the differential heating of the atmos-
phere on each side of the front, leading in turn to surface
wind stress convergence (divergence) over warm (cold) SST
patterns (Hsu 1984; Lindzen and Nigam 1987; Warner et al.
1990; Wai and Stage 1989). We can visualize this mecha-
nism by the analogy of the sea-breeze circulation proposed
by Hsu (1984), according to which the warm ocean along
the Gulf Stream (GS) has the role of the land. Minobe et al.
(2008) noticed that the wind divergence pattern co-locates
with the Laplacian of sea-level pressure (∇²SLP), and that
both have many similarities with the pattern of the sign-
reversed Laplacian of SST (−∇²SST).

The importance of these mechanisms is not limited to
the fundamental understanding of the Earth system. The
SST is not only modifying the near surface wind fields but
also can affect the free troposphere influencing that way the
Northern Hemisphere from weather to climate time-scales
(Minobe et al. 2008, 2010; Ma et al. 2015; Kobashi et al.
2008; Lee et al. 2018, Roberts et al. in press). This is through
the wind convergence in the MABL that anchors a band
of precipitation and deep upward motion extending to the
free troposphere, with a remote impact by forcing of atmos-
pheric planetary waves (Minobe et al. 2008). Moreover, SST
induced modifications of the wind speeds can affect cyclo-
genesis, and influence the storm tracks and storm intensity
(Xie et al. 2002; Nakamura et al. 2004; Small et al. 2014;
Woollings et al. 2010). In addition, the near surface wind
stress response to the SSTs can possibly feed back on the
ocean by altering, for example, Ekman pumping (Chelton
et al. 2007; Spall 2007).

Due to the coarse resolution of climate models, we hypothesize that the impact of the ocean on the atmosphere at small scales is not correctly represented, with implications for the atmospheric dynamics such as position and strength of the storm tracks (Piazza et al. 2016; Parfitt et al. 2017; Small et al. 2014, 2019) and that higher resolution is needed. Recently, a new set of global high resolution climate simulations following the HighResMIP protocol (Haarsma et al. 2016) have become available. Using these simulations, we investigate the ocean–atmosphere coupling in the GS region that is characterized by sharp fronts and large eddy activity. This region has also been the focus of earlier studies of VMM and PAM (e.g. Wai and Stage 1989; Minobe et al. 2008; Chelton and Xie 2010). Air–sea interaction over the GS region is of great interest since it affects cyclogenesis and storm track characteristics (Small et al. 2014), induces deep tropospheric response and influences the large scale atmospheric circulation (Minobe et al. 2008). The High-ResMIP protocol, requiring all simulations to be performed with a low and high resolution model configuration, allows us to investigate the spatial resolution dependence of the ocean–atmosphere interaction. The impact of increased resolution in HighResMIP simulations on climate and variability has been documented in several studies, including Moreno-Chamarro et al. (2021) who discuss the impact on biases. We have focused here on the VMM and PAM mechanisms and compared the results with available observations.

In the following Sect. 2 we describe the models’ con-
figurations, the observation/reanalysis datasets and the
techniques used for this analysis. In Sect. 3 we present our
results, followed by a discussion in Sect 4. Finally, a sum-
mary of the conclusions is given in Sect. 5.

## 2 Data and methods

2.1 Data

In the current study, we analyse the “control-1950” simula-
tions of the High Resolution Model Intercomparison Project
(HighResMIP) (Haarsma et al. 2016), which is part of phase
6 of the Coupled Model Intercomparison Project (CMIP6)
(Eyring et al. 2016). This experiment consists of 100-year
coupled simulations of different horizontal resolutions,
with external forcings kept constant at 1950’s values. All
the simulations are forced with greenhouse gases (GHG)
and aerosols loadings using 1950s climatology (fixed val-
ues). The initialization of the coupled runs was provided
by ~50-year spin-up integrations starting from the observed
ocean and atmosphere state in 1950, with fixed forcing rep-
resentative of the 1950s. Details of the experimental set-up

<figure><img src="image_4.png" /></figure>
of the HighResMIP spin-up and control-1950 simulations are described in Haarsma et al. (2016).

The six GCMs examined in this study are from the EU funded Horizon-2020 PRIMAVERA project. Each of them consists of more than one configuration (15 in total). The differences between the configurations of each model are appertain to the horizontal resolution of the atmosphere, ocean or both components. Resolution-dependent param- eterizations might differ among the configurations. Particu- larly, the Gent and McWilliams (1990) scheme for the ocean eddy parameterization is deactivated in the HadGEM3- GC3.1, ECMWF-IFS and CMCC-CM2 configurations of oceanic resolution 25km or higher. For each model con- figuration, one member was analyzed. We argue that for the ocean-atmosphere coupling along the GS, natural variability is sufficiently sampled in a single 100-year integration. This is supported by the fact that applying our analysis to only the first 50 years and only the last 50 years of the simulations, yields similar results (Supplemental Fig. S1). The model configurations and their characteristics are summarized in Table 1. Detailed descriptions of the configurations can be found in Roberts et al. (2019) for HadGEM3-GC3.1, in Rob- erts et al. (2018b) for ECMWF-IFS, in Haarsma et al. (2020) for EC-Earth3P, in Voldoire et al. (2019) for CNRM-CM6, in Gutjahr et al. (2019) for MPI-ESM1-2 and in Cherchi et al. (2019) for CMCC-CM2. A short description of these models is provided in the Supplemental Material (text S1).

For the atmospheric components, the effective resolution in addition to the nominal resolution is listed in Table 1 as computed in Klaver et al. (2020). The effective resolution is based on the shape of the kinetic energy spectrum and

provides an estimate of the boundary of the simulated scales that are not affected by the resolution of the model.

Satellite observations combined with reanalysis data, for the period January 2007 to December 2018, are used for comparison with the model output. This period is cho- sen because of the availability of high resolution satellite observations. The 11-year observational period is, however, much smaller than the 100-year control simulation which can increase the uncertainties. Below we briefly describe the three different observational and reanalysis data-sets used in this study.

SST measurements are obtained from the National Oce- anic and Atmospheric Administration (NOAA) 0.25° daily Optimum Interpolation SST version 2 (OISST.v2) (Ban- zon et al. 2016). This data-set is constructed by combin- ing observations from the Advanced Very High Resolution Radiometer (AVHRR) and other platforms (ships, buoys) on a regular global grid. Optimum interpolation (OI) is used to produce a spatially complete map.

Wind stress divergence and curl are obtained from Advanced SCATterometer (ASCAT) on MetOp-A (Figa- Saldaña et al. 2002). It is a product of the European Organi- zation for the Exploitation of Meteorological Satellites (EUMETSAT) Ocean and Sea Ice Satellite Application Facility (OSI SAF) provided through the Royal Netherlands Meteorological Institute (KNMI). We use daily values from the L3 REP product, sampled on a 0.25° grid.

SLP is derived from the European Centre for Medium Range Weather Forecasts Reanalysis 5 atmospheric ReA- nalysis (ERA5) (Hersbach et al. 2020). Data processing for ERA5 is carried out by using ECMWF’s Earth System

<table><caption>**Table 1** Specifications of the different simulations</caption><thead><tr><th rowspan="2">Models</th><th rowspan="2">Configurations</th><th colspan="3">Atmosphere</th><th colspan="2">Ocean</th></tr><tr><th>Nominal Resolution (km)</th><th>Effective Resolution (km)</th><th>Vertical Levels</th><th>Nominal Resolution</th><th>Vertical Levels</th></tr></thead><tbody><tr><td rowspan="4">HadGEM3 GC3.1</td><td>HH</td><td>50</td><td>185</td><td rowspan="4">85</td><td>1/12°</td><td rowspan="4">75</td></tr><tr><td>HM</td><td>50</td><td>185</td><td>0.25°</td></tr><tr><td>MM</td><td>100</td><td>364</td><td>0.25°</td></tr><tr><td>LL</td><td>250</td><td>≥ 625</td><td>1°</td></tr><tr><td rowspan="3">ECMWF IFS</td><td>HR</td><td>25</td><td>185</td><td rowspan="3">91</td><td>0.25°</td><td rowspan="3">75</td></tr><tr><td>MR</td><td>50</td><td>253</td><td>0.25°</td></tr><tr><td>LR</td><td>50</td><td>253</td><td>1°</td></tr><tr><td rowspan="2">EC-Earth3P</td><td>HR</td><td>50</td><td>238</td><td rowspan="2">91</td><td>0.25°</td><td rowspan="2">75</td></tr><tr><td>LR</td><td>100</td><td>351</td><td>1°</td></tr><tr><td>CNRM CM6</td><td>HR</td><td>50</td><td>313</td><td rowspan="2">91</td><td>0.25°</td><td rowspan="2">75</td></tr><tr><td></td><td>LR</td><td>250</td><td>≥ 625</td><td>1°</td></tr><tr><td rowspan="2">MPIESM-1-2</td><td>XR</td><td>50</td><td>256</td><td rowspan="2">95</td><td>0.4°</td><td rowspan="2">40</td></tr><tr><td>HR</td><td>100</td><td>364</td><td>0.4°</td></tr><tr><td>CMCC CM2</td><td>VHR4</td><td>25</td><td>182</td><td rowspan="2">26</td><td>0.25°</td><td rowspan="2">50</td></tr><tr><td></td><td>HR4</td><td>100</td><td>571</td><td>0.25°</td></tr></tbody></table>
model IFS, cycle 41r2. ERA5 has spatial resolution of 31 km, hourly time resolution and 137 vertical levels. We use monthly values on a 0.25° grid.

## 2.2 Methods

We analyze the three-month boreal winter season from December to February (DJF). The rationale for this choice is that due to the increased temperature difference between the atmosphere and ocean, and the existence of stronger winds, the air-sea coupling is stronger in winter (Minobe et al. 2010; Putrasahan et al. 2013).

We use monthly values of the surface wind stress components, SST and SLP from the 100-year control simulations. We choose to work with monthly values in order to remove the effects of energetic synoptic weather variability (Song et al. 2009; Chelton et al. 2007; Chelton and Xie 2010). Since we are investigating two mechanisms which explain the atmospheric response to the ocean, we first bilinearly interpolate the SST onto the atmospheric grid of each model configuration. In that way, our computations are not affected by the oceanic variability at spatial scales that cannot be resolved by the atmosphere. Following O'Neill et al. (2003), we compute the downwind SST gradient as the projection of the SST gradient on the wind stress vector. The crosswind SST gradient is the component of the SST gradient that is perpendicular to the projection of the SST gradient on the wind stress vector. In addition, the wind stress divergence and curl, the $\nabla^2$SLP and the $-\nabla^2$SST are computed. For the calculation of the derivatives, a central difference scheme is used.

We compute the detrended anomalies of the fields. Next, following Roberts et al. (2016), we applied a high-pass spatial filter to these anomalies, by subtracting smoothed fields that are computed by the application of a 18° longitude × 6° latitude box car filter from the original fields. This filter is similar to a loess spatial high-pass filter with half-power points at 30° longitude and 10° latitude which was applied in previous studies to isolate the mesoscale signal (Bryan et al. 2010; Maloney and Chelton 2006).

The linear relationship, in terms of least squares regression, is investigated between the anomalies of the three aforementioned pairs of fields ([wind stress divergence, downwind SST gradient], [wind stress curl, crosswind SST gradient], [$\nabla^2$ SLP, $-\nabla^2$SST]). The slope of the linear fit is the so-called coupling coefficient, as its magnitude indicates the strength of the air-sea coupling. To compute the coupling coefficients, we create monthly binned histograms of the SST-related fields in the GS domain (30°–50° N, 80°–40° W, boxed domain in Fig. 1). Then, we create their monthly histograms by determining the averaged atmospheric fields for each bin of the associated SST-related fields. Next the monthly histograms over the analyzed time period are averaged and plotted in a scatter

<figure><img src="image_5.png" /><figcaption>**Fig. 1** Wintertime climatology of ASCAT near surface wind stress divergence (red) and convergence (blue), in Pa/10000km (colours) and NOAA OISST sea surface temperatures (contours). The sea surface temperature contour interval is 2° C. The thick black line denotes the 18° C isotherm. The red boxed region (40°–80° W, 30°–50° N) is used for the analysis on coupling coefficients</figcaption></figure>

diagram. The linear fit is made over the mean histogram, resulting in one coupling coefficient for each configuration.

For a fair comparison with the observations, we handle the latter in the same way as the model data. The only difference is that for the ASCAT products we use daily values of the wind stress divergence/curl and then we compute the monthly averages for the points with 20 or more available observations during a month. We proceed this way, since due to the inhomogeneous distribution of the observations over space and time, the fields appear noisy and this ‘noisy’ character is further amplified by the derivatives operation. The coupling coefficient though, is not affected significantly by first averaging the winds stress and then computing the divergence (Supplemental Fig. S2).

The use of the monthly mean values in our analysis could be considered a limitation. As we do not compute the variables of interest with higher time-resolved data (e.g. daily) before averaging, errors may occur in the non-linear terms of the downwind and crosswind SST gradients. For this reason, we also compute the coupling coefficient by first applying the operators on the daily values before averaging, for the configurations of the HadGEM3-GC3.1 model and the observations, in order to investigate if any significant differences are induced. Any regridding, required before the averaging, was made using a conservative interpolation method. Even though the magnitude of the coupling differs between the two analyses, the dependence of the coupling on the horizontal resolution is not sensitive to the order of the aforementioned steps (not shown here).

## 3 Results

The 11-year winter mean climatology of ASCAT near surface wind stress divergence and NOAA OISSTv.2 SSTs are shown in Fig. 1. As previous studies have shown (Minobe

<figure><img src="image_6.png" /></figure>
et al. 2010; Chelton et al. 2004; Takatama et al. 2012), observations reveal wind stress divergence (convergence) along the cold (warm) flank of the GS. Strong divergence (red in Fig. 1) is observed in the region of Grand Banks where the GS turns to the north. Westerlies cross the GS resulting in divergence due to the wind speed differences on each side of the strong SST gradients.

The ability of the models to reproduce these wind features and their dependence on spatial resolution is shown in Fig. 2, where the 100-year DJF climatology of wind stress divergence and SST contours are presented for the different configurations of the models. Generally, all models are capable of capturing the divergence and convergence zones on the cold and warm flank of the GS respectively. However, configurations with a low-resolution oceanic component of about 1°, where the ocean eddy dynamics are parameterized, display much lower absolute values of divergence. Especially east of Newfoundland, the divergence values are smaller by a factor of four or more compared to the observations. This is associated with the much weaker

SST gradients of the eddy-parameterized configurations (Supplemental Fig. S3). By increasing the oceanic resolution from 1° to 0.4°–0.25°, representing eddy-permitting dynamics, the strong SST gradients at about 50° W around the Grand Banks are better resolved, in line with Kirtman et al. (2012). However, a substantial underestimation of the divergence field is observed even in eddy-permitting models. Even though there is a large spread among the models, enhancing the atmospheric resolution improves the simulated divergence field. The HadGEM3-CG3.1-HH, HM, MM and CMCC-CM2-VHR4 are the configurations with the best performance, with divergence patterns and magnitudes comparable to observations.

In observations, the GS separates from the US coast near Cape Hatteras. Only HadGEM3-CG3.1-HH with an ocean eddy-resolving resolution of 1/12° and atmosphere resolution of 50 km correctly simulates this separation. Coarser resolutions (oceanic and atmospheric) result in separations further to the north along the eastern coast of North America. This is a well-known feature and has

<figure><img src="image_7.png" /><figcaption>**Fig. 2** Same as in Fig. 1 but for the model simulations. The configurations of each model are presented in each row, with higher resolution configurations on the left and lower on the right. From right to</figcaption></figure>

left: eddy-parameterized models, eddy-permitting models, eddy-permitting models with further increase of the atmospheric resolution, eddy-resolving models
been addressed in more detail in previous studies (Bryan et al. 2007; Chassignet and Xu 2017; Seidov et al. 2019). Consequently, the associated convergence zones of wind stress are also shifted northwards for coarser resolutions.

# 3.1 Coupling through vertical mixing mechanism

## 3.1.1 Spatial anomalies patterns and correlation

To illustrate the resolution dependence of downwind SST gradient and near surface wind stress divergence and their collocation, we present in Fig. 3 (left) the high-pass filtered downwind SST gradient and the near surface wind

stress divergence anomalies for the different configurations of the HadGEM3-GC3.1 model and observations (ASCAT-OISST). For all configurations, the relationship between the two fields is evident as wind stress divergence anomalies are associated with positive downwind SST gradient anomalies, while convergence anomalies are found over negative downwind SST gradient anomalies. The spatial covariability of the fields is large for all configurations. However, the ocean eddy-parameterized (1°) HadGEM3-GC3.1-LL results in lower absolute values of anomalies for both fields by a factor of five, compared to the eddy-permitting (0.25°) and eddy-resolving (1/12°) configurations. In addition, the eddy-parameterized

<figure><img src="image_8.png" /><figcaption>**Fig. 3** Right: DJF high-pass filtered anomalies of wind stress divergence (contours) and downwind SST gradients (°C/100km) (colours), for the HadGEM3-GC3.1 configurations and ASCAT/NOAA-OISST. Left: DJF high-pass filtered anomalies of wind stress curl (contours) and crosswind SST gradients (°C/100km) (colours), for the HadGEM3-GC3.1 configurations and ASCAT/NOAA-OISST.</figcaption></figure>

Means are over one representative December-February period as the simulations do not represent actual years. The contours intervals (c) differ for each resolution for better visualization and are noted (in Pa/10000km) at the left top corner in each case. The colorbar range differs for the LL configuration
configuration exhibits patterns of a much larger scale. Configurations with eddy-permitting oceanic resolution, result in increased magnitude of the downwind SST gradient anomalies and smaller spatial features. The collocation of the fields becomes more evident, by further increasing the atmospheric resolution, from 100 to 50 km. Further increase of the oceanic resolution to eddy-resolving, leads to enhanced values of downwind SST gradient anomalies without a substantial associated increase in the wind stress divergence. Similar dependence on resolution is observed for the crosswind SST gradient and wind stress curl variability as illustrated in Fig. 3 (right).

The other models reveal a similar resolution dependence of the SST influence on near surface winds as the HadGEM3-GC3.1 model (Supplemental Fig. S4–8). All the eddy-parameterized configurations are incapable of reproducing the small scale features. Due to their coarse resolution, it is not possible to resolve the mesoscale activity of the GS resulting in small anomalies of the aforementioned fields. Increasing oceanic and atmospheric resolution results

in higher amplitude of anomalies with smaller scale structure and better agreement with observations.

To further investigate the relation between the near surface winds and the SSTs, and its variability in space in the GS region, Fig. 4 shows the temporal correlation between the downwind SST gradient and the wind stress divergence for all models and observations/reanalysis. All models reveal strongest correlation along the strong SST gradient of the GS, as expected. Models with eddy-parameterized ocean components (1°) exhibit low correlations over a large area. The largest correlations with values larger than 0.5 are obtained at the end of the GS region east of Nova Scotia. Approximately 50% of the values of the correlation coefficient in eddy-parameterized models are lower than 0.2.

The eddy-permitting (0.4°–0.25°) models generally reveal higher correlations up to 0.8 over a much larger area than the eddy-parameterized models. The overall increase of the correlation in the ECMWF-IFS-MR, reveals that the oceanic resolution has a primary role in this improvement. However,

<figure><img src="image_9.png" /><figcaption>**Fig. 4** Correlation of timeseries of spatially high-pass filtered wind stress divergence and downwind SST gradient anomalies. Correlations are computed using monthly values for the winter period (DJF).</figcaption></figure>

The stippling shows statistically non-significant correlations at 95% significance level using two-sided student's t-test
even though there is spread between the models, the correlation appears quite sensitive to the atmospheric resolution too. We note that the band of positive correlation along the GS is further enhanced, when only the atmospheric resolution increases (HadGEM3-HM, CMCC-VHR4, MPIESM-XR and ECMWF-IFS-HR), displaying a 50–100% increase of the area with correlations larger than 0.6. The only eddy-resolving (1/12°) model HadGEM3-CG31-HH displays lower correlations than the eddy-permitting version with the same atmospheric resolution HadGEM3-CG31-HM. Possible reasons for this unexpected behavior will be discussed below.

In general, ASCAT data display lower correlations in comparison with the models, except for the eddy-parame-terized ones. Similar behavior was noticed by Roberts et al. (2016) where the correlation between the wind stress and the SSTs was examined. This can possibly be explained by the temporal and spatial inhomogeneity of the satellite observa-tions. In Fig. 1, we can notice the “noisy” character of the observed divergence field, which is indicative of a certain

degree of uncertainty in the observations due to the limited period covered by the data.

Figure 5 displays the correlation coefficient between the high-pass filtered crosswind SST gradient anomalies and the associated near surface wind stress curl (note that reversed colors are used compared to Fig. 4). The dependence of the amplitude and spatial structure of correlation on the spatial resolution of the model components is similar to the one found for the downwind component of SST gradient and the divergence of the wind stress. Negative (positive) cross-wind SST gradient anomalies are associated with cyclonic (anticyclonic) wind stress anomalies resulting in anti-correlation. The correlation is systematically lower for vorti-city than for divergence. The possible reasons for this will be discussed below. In addition, while the HadGEM3-HH shows lower correlations than the HadGEM3-HM for the relationship between the downwind SST gradient and the wind stress divergence, here the two configurations display similar results.

<figure><img src="image_10.png" /><figcaption>**Fig. 5** Correlation of timeseries of spatially high-pass filtered wind stress curl and crosswind SST gradient anomalies. Correlations are computed using monthly values for the winter period (DJF). The stip-</figcaption></figure>

pling shows statistically non-significant correlations at 95% significance level using two-sided student's t-test

<figure><img src="image_11.png" /></figure>
### 3.1.2 Coupling coefficient

To quantify the strength of air-sea coupling in the GS region, we focus on the red box area defined in Fig. 1 (30°–50° N, 80°–40° W). Figure 6 shows the binned scatter plots of spatially high-pass filtered downwind SST gradients and wind stress divergence for the HadGEM3-GC3.1 configurations. The strength of the VMM is determined by the slope of the linear regression. Slopes are defined to be statistically significant when the *p*-value <= 0.05, based on a Wald test with t-distribution. The HadGEM3-GC3.1-LL results in a four times weaker slope than observed, $s = 0.18 \times 10^{-2}$ Pa/°C. The other three configurations of the HadGEM3-GC3.1 model exhibit much higher coupling with $s = 0.52 \times 10^{-2}$ Pa/°C, $s = 0.7 \times 10^{-2}$ Pa/°C and $s = 0.6 \times 10^{-2}$ Pa/°C for the MM, HM and HH respectively. The lower coupling of HH configuration in comparison with the HM is in agreement with the results in Fig. 4.

The results for all the models and the observations are synthesized in Table 2, displaying not only the coupling coefficients between the wind stress divergence and the downwind SST gradient but also those between the wind stress curl and the crosswind SST gradient. In addition, the 95% confidence intervals of the slopes are presented, together with the correlation coefficients in parenthesis. All the values are statistically significant.

All the eddy-parameterized models (HadGEM3-LL, ECMWF-LR, EC-Earth-LR and CNRM-LR) exhibit linear relationships with coupling coefficients not exceeding $0.3 \times 10^{-2}$ Pa/°C between the two pairs of fields. Models of eddy-permitting resolutions of 0.4° and 0.25°, give coupling coefficients of larger magnitude. One part of this increase is attributable to the better representation of the mesoscale oceanic activity by the finer oceanic component. This is evident by the comparison of the LR and MR configurations of the ECMWF-IFS, where there is an increase of the coupling coefficient by only refining the oceanic resolution to eddy-permitting. However, the CMCC-CM2 HR4 produce coupling coefficients that are similar to the eddy-parameterized models especially in the case of the relationship between the downwind SST gradient and the wind stress divergence. On the other hand, the VMM is sensitive to the atmospheric resolution too. Increasing the atmospheric resolution results in enhanced coupling and leads to values similar in amplitude with the observed for the case of wind stress divergence. The HadGEM3-HH shows weaker or same coupling strength compared to the HadGEM3-HM.

The coupling coefficients between the wind stress curl and the crosswind SST gradient are for all models and their configurations lower than the ones between the downwind component of SST gradient and the divergence of the wind stress. This is in agreement with the observations, but the differences in the models are larger. Reduction is also

<figure><img src="image_12.png" /><figcaption>**Fig. 6** Binned scatter plots for the HadGEM3-GC3.1 configurations and observations of high-pass filtered wind stress divergence and downwind SST gradients over the region of Gulf Stream (red box in Fig. 1). The calculations were made by averaging the monthly binned histograms for Dec-Feb, over the years 1950-2050. The bin range is 0.1°C/100km. The error bars represent the ±1 standard deviation over the averaged samples in each bin. The red line is the linear fit in terms of least squares. Statistically significant slopes are represented with green color text. The criterion for a significant value is $p \le 0.05$ (based on a Wald test using t-distribution)</figcaption></figure>

noticed in the temporal correlations as discussed above. This result is consistent with other studies and has been attributed to many causes. Chelton et al. (2001) and Kilpatrick et al. (2014) attribute this to differences in the Marine Atmospheric Boundary Layer (MABL) response for cross-front and along-front winds. The finite timescale for adjustment of the MABL, allows the boundary layer to come into equilibrium with the SST for an along-front flow, in contrast to a cross-front flow where the downwind SST gradients are sharper (Chelton et al. 2001). That may enhance the wind stress divergence response on the downwind SST gradient
<caption>**Table 2** Wintertime (Dec-Feb) coupling coefficients and correlations (in parenthesis) for the models and observations. First column: near surface wind stress divergence ($\nabla \cdot \boldsymbol{\tau}$) and downwind component of SST gradient (Downwind $\nabla$SST). Second column: near surface wind stress curl ($\nabla \times \boldsymbol{\tau}$) and crosswind component of SST gradient (Crosswind $\nabla$SST). The models' results are based on 100-year simulations (1950–2050) while the observations (last row) consisted of 11-year records (2007–2018). The margin of errors indicates the 95% confidence interval. All the values are statistically significant. The criterion for a statistically significant value is $p \le 0.05$ (based on a Wald test using t-distribution)</caption>

<table><thead><tr><th>Models</th><th>Config.</th><th>Atm. Resol. (km)</th><th>Oceanic Resol.</th><th>∇·τ and Downwind ∇SST</th><th>∇×τ and Crosswind ∇SST</th></tr></thead><tbody><tr><td rowspan="4">HadGEM3<br>GC3.1</td><td>HH</td><td>50</td><td>1/12°</td><td>0.6±0.022 (0.99)</td><td>-0.4±0.019 (-0.99)</td></tr><tr><td>HM</td><td>50</td><td>0.25°</td><td>0.7±0.021 (1)</td><td>-0.4±0.022 (-0.99)</td></tr><tr><td>MM</td><td>100</td><td>0.25°</td><td>0.52±0.015 (1)</td><td>-0.32±0.014 (-0.99)</td></tr><tr><td>LL</td><td>250</td><td>1°</td><td>0.18±0.038 (0.84)</td><td>-0.11±0.046 (-0.61)</td></tr><tr><td rowspan="3">ECMWF<br>IFS</td><td>HR</td><td>25</td><td>0.25°</td><td>0.75±0.014 (1)</td><td>-0.38±0.016 (-0.99)</td></tr><tr><td>MR</td><td>50</td><td>0.25°</td><td>0.52±0.024 (0.99)</td><td>-0.27±0.019 (-0.98)</td></tr><tr><td>LR</td><td>50</td><td>1°</td><td>0.31±0.051 (0.9)</td><td>-0.19±0.046 (-0.81)</td></tr><tr><td rowspan="2">EC-Earth3P</td><td>HR</td><td>50</td><td>0.25°</td><td>0.58±0.012 (1)</td><td>-0.45±0.009 (-1)</td></tr><tr><td>LR</td><td>100</td><td>1°</td><td>0.26±0.035 (0.93)</td><td>-0.02±0.026 (-0.93)</td></tr><tr><td rowspan="2">CNRM<br>CM6</td><td>HR</td><td>50</td><td>0.25°</td><td>0.8±0.021 (1)</td><td>-0.52±0.01 (-1)</td></tr><tr><td>LR</td><td>250</td><td>1°</td><td>0.3±0.058 (0.87)</td><td>-0.12±0.042 (-0.7)</td></tr><tr><td rowspan="2">MPIESM-1-2</td><td>XR</td><td>50</td><td>0.4°</td><td>0.76±0.019 (1)</td><td>-0.54±0.012 (-1)</td></tr><tr><td>HR</td><td>100</td><td>0.4°</td><td>0.6±0.014 (1)</td><td>-0.48±0.011 (-1)</td></tr><tr><td>CMCC<br>CM2</td><td>VHR4</td><td>25</td><td>0.25°</td><td>0.68±0.009 (1)</td><td>-0.38±0.02 (-0.99)</td></tr><tr><td></td><td>HR4</td><td>100</td><td>0.25°</td><td>0.34±0.021 (0.98)</td><td>-0.22±0.013 (-0.98)</td></tr><tr><td>ASCAT and NOAA-OISSTs</td><td></td><td></td><td></td><td>0.75±0.023 (1)</td><td>-0.67±0.022 (-1)</td></tr></tbody></table>

due to the increased turbulent fluxes in the non-equilibrium state. Differences in the ocean current speed effects on the cross-isotherm and along-isotherm wind stress can also cause differences in the divergence and curl responses (e.g. Chelton et al. 2004). O'Neill et al. (2010) explain the weaker wind curl response to the SST gradients, by separating the effect of the SST-induced wind speed and direction gradients on the wind curl and divergence. They found that while the SST-induced wind speed anomalies contribute equally to the curl and divergence, the SST-forced wind-direction gradients diminish the wind curl response to crosswind SST gradients and enhances the wind divergence response to downwind SST gradients.

## 3.2 Coupling through pressure adjustment mechanism

Figure 7 shows the observed winter-averaged $\nabla^2$SLP (top) and the $-\nabla^2$SST (bottom), by the ERA5 and the NOAA-OISST products. Negative pressure anomalies are located at the warm flank of the GS while positive pressure anomalies are found over the cold flank. $\nabla^2$SLP exhibits many similarities with the $-\nabla^2$SST and the wind stress divergence (Fig. 1). These results are consistent with Minobe et al. (2008) and suggest the existence of PAM over the GS region, according to which warm (cold) anomalies induce low (high) pressure

anomalies and low level convergence (divergence). This is also evident in the models (Supplemental Fig. S9–10). However, the eddy-parameterized models exhibit much weaker $-\nabla^2$SST and a north-shifted band of positive $\nabla^2$SLP, due to the unrealistic separation of the GS. Higher resolution configurations are in better agreement with observations.

### 3.2.1 Correlation

Spatial patterns of the temporal correlation between the high-pass filtered $-\nabla^2$SST and the associated $\nabla^2$SLP (Fig. 8), highlight a similar dependence on resolution for PAM as for VMM. The eddy-parameterized models only display correlations larger than 0.5° east of Nova Scotia, due to the weaker SST gradients and the shifted, more zonally oriented path of the GS. The eddy-permitting models show a systematic increase of the correlations, exceeding the observed, with values above 0.5° starting from Cape Hatteras and continue all along the SST front. Enhancing atmosphere resolution leads to higher correlations and improves further the spatial structure. The eddy-resolving HadGEM3-CG3.1-HH shows no improvement compared to the eddy-permitting HadGEM3-CG3.1-HM. A second band of increased correlation is noticed south of 30° in the extratropics, both in observations and some models. However, the climatological $-\nabla^2$SST and $\nabla^2$SLP do not reveal a signal in that region
<figure><img src="image_13.png" /><figcaption>**Fig. 7** Wintertime climatology of ERA5 $\nabla^2$SLP (top) $(10^{-9} Pa/m^2)$ and NOAA OISSTv2 $\nabla^2$SST (bottom) $(10^{-10} C/m^2)$. The contours represent the NOAA OISST sea surface temperatures with contour interval 2° C. The thick black line denotes the 18° C isotherm</figcaption></figure>

(Fig. 7), suggesting that this band is not associated with the PAM mechanism over the GS. We speculate that this band of correlation is a reflection of the enhanced ocean–atmosphere coupling in the tropics (Kushnir et al. 2002).

### 3.2.2 Coupling coefficient

Since the Laplacian operator acts as a high-pass filter, it is sensitive to the grid on which it is computed. We argue that the linear relationship between the Laplacians is scale-dependent, and for this reason a fair comparison between the different resolutions is not possible. In order to demonstrate this we compute the $-\nabla^2$SST and the $\nabla^2$SLP both from ERA5, with the same numerical scheme, first using every grid point (effective resolution of 0.25 degrees, “HR”) and subsequently using every fourth grid point (effective resolution of 1 degree, “LR”). The coupling coefficient for the HR is 0.26 while for the LR is 0.41 (Supplemental Fig. S11). This increase of the coupling coefficient, by computing the Laplacians over larger distances using the exact same dataset, indicates that comparing the coupling coefficients between the Laplacians across different model resolutions,

is not simply a comparison between the realism of the models. On the other hand, the temporal correlations are not affected by the computation of the Laplacians on different grid spacings (Supplemental Fig. S12). For this reason, we do not proceed to a comparison of the coupling coefficients for the PAM mechanism and we base our conclusions on the temporal correlations, which represent a more reliable metric for this case.

## 4 Discussion

Generally, we see that horizontal model resolution plays an important role in the accurate representation of air–sea interactions. Even though, the eddy-parameterized configurations are able to capture the oceanic forcing, they underestimate the coupling strength and exhibit lower temporal/ spatial correlation and coupling coefficients. Increasing the resolution leads to simulations that are in better agreement with the observations. The oceanic resolution appears to be a key feature for more realistic simulations, which is in agreement with other studies (e.g. Bryan et al. (2010); Bellucci et al. (2021); Putrasahan et al. (2013); Parfitt et al. (2017)). However, only after a further increase of the atmospheric resolution, we obtain results comparable to the observations.

Bellucci et al. (2021) investigate the air–sea interaction over the GS region using the same simulations as here. Looking at the correlation and co-variance patterns between the SST and the turbulent heat fluxes (THF) they also find that enhancing the spatial resolution leads to better representation of the coupled ocean–atmosphere processes. However, they attribute the improvements on the finer oceanic resolution, since their results show less sensitivity to the resolution of the atmospheric component. Our results, on the other hand, indicate that for an optimal coupling, in addition to an eddy-permitting oceanic component, the atmospheric resolution should not be much coarser than the ocean resolution. The results from the eddy-permitting models for which are available two atmospheric resolutions (HadGEM3-GC3.1, ECMWF-IFS, MPIESM-1-2, CMCC-CM2), support this hypothesis. Experiments with eddy-parameterized oceanic components coupled with high-resolution atmospheric models could provide a better understanding of the relative importance of the atmospheric and oceanic resolutions, yet such simulations are not available in this set of experiments.

A peculiar result is that for the correlations as well as the coupling coefficients, the eddy-resolving HadGEM3-CG3.1-HH shows lower or similar values compared to the eddy-permitting HM version of this model with the same atmospheric component. This result is in contrast with the change from eddy-parameterized to eddy-permitting models, where there is a consistent increase in correlations and coupling
<figure><img src="image_14.png" /><figcaption>**Fig. 8** Correlation of timeseries of spatially high-pass filtered $\nabla^2$SLP and $-\nabla^2$SST anomalies. Correlations are computed using monthly</figcaption></figure>

values for the winter period (DJF). The stippling shows statistically non-significant correlations at 95% significance level using two-sided student's t-test

coefficients among all models. However, this is consistent with other studies. Roberts et al. (2016), comparing an eddy-permitting to an eddy-resolving version of the Met Office climate model, find no significant improvement of the SST and winds stress relationship. Bellucci et al. (2021), find little improvements on the relationship between the SST and THF when comparing an eddy-resolving configuration of the MPI-ESM model to the HighResMIP simulations of the same model.

Presently, we do not have a clear explanation for this behavior, but we note that the atmospheric resolution of 50 km is much coarser than the ocean resolution of 1/12° (~10km) and this ocean eddy variability at those small scales cannot be seen by the atmosphere. If in the change from ocean eddy-permitting to eddy-resolving models a significant fraction of the ocean variability is transferred to those small scales, this will not be resolved by the atmosphere, resulting in weaker or similar correlations and coupling coefficients. Moreton et al. (2021), using three different

configurations of the HadGEM3-GC3.1 model, they found evidence that increasing the oceanic resolution at a constant atmospheric resolution can deteriorate the representation of air-sea interaction. However, Bryan et al. (2010), comparing two eddy-resolving versions of the NCAR Community Climate System Model (CCSM3.5) with 0.5° and 0.25° atmospheric resolution, find no significant changes. They suggest that implementing more accurate parameterizations of the atmospheric boundary layer would be more beneficial than increasing the atmospheric resolution. To understand better these results, more in-depth analyses, in a multi-model approach, is required.

Despite the general tendency of the models to exhibit a better agreement with observations by enhancing the resolution, there is an inter-model spread which can be attributed to factors other than the spatial resolution. The coupling strength can be affected by the vertical resolution of the atmosphere and the number of vertical levels in the boundary layer (Chelton and Xie 2010). Moreover, different

<figure><img src="image_15.png" /></figure>
physical parameterizations in the convection and boundary layer schemes can significantly alter the surface wind response to the SST variability (Song et al. 2009; Perlin et al. 2014).

An important impact of the underestimation of the air-sea coupling that is not studied here is that on the storm track activity. Piazza et al. (2016) show that a simulation forced by smoothed SST fields over the GS region, leads to a northward shift of the storm track density towards the North America coast. Similar results are obtained by Woollings et al. (2010) and Small et al. (2014). This shift of the main storm track can also result in reduced storm activity over remote regions such as Greenland and the Mediterranean (Piazza et al. 2016). Haarsma et al. (2019) analysing coupled seasonal forecast experiments at different resolutions, find stronger injection of wave activity into the storm track with increasing oceanic resolution. In order to assess these impacts further, dedicated experiments have been conducted, and the results will be presented in future work.

Finally, even though we focus our analysis on the GS, the ocean forcing on the atmosphere has been also identified in other regions with strong ocean fronts and eddy activity, such as the Kuroshio and its Extention (KE) (e.g. Nonaka and Xie 2003; Putrasahan et al. 2013), the Antarctic Circumpolar Current (e.g. O'Neill et al. 2003), the Agulhas Return Current (e.g. O'Neill et al. 2005; Song et al. 2009) and the Brazil-Malvinas confluence (e.g. Tokinaga et al. 2005; Pezzi et al. 2009). Both observational and model analysis indicate that while the coupling coefficients are comparable for the KE and GS, higher values are retrieved in the Southern Ocean (Chelton et al. 2004; Maloney and Chelton 2006; Bryan et al. 2010). Despite the regional differences of the coupling magnitude, which might be related to differences in the MABL structure, we generally expect similar resolution dependence of the coupling strength in all the regions. Past studies have shown that in eddy-parameterized models the coupling is weak or even absent over all regions with intense mesoscale activity, indicating the oceanic resolution as a key feature for the air-sea coupling (e.g. Bryan et al. 2010; Small et al. 2019). However, we should be careful extrapolating our results to other regions since discrepancies may arise due to different intensities and orientations of the ocean fronts and/or air-sea temperature contrasts. For example, Bellucci et al. (2021) found that while eddy-parameterized models are able to reproduce the ocean forcing on the atmosphere over GS and KE, this is not the case for the Southern Ocean. This may be related to geographical dependencies of the transition length scales setting the transition from ocean-driven to atmosphere-driven regimes (Bishop et al. 2017; Small et al. 2019).

# 5 Conclusions

The GS region, with its sharp SST front and large eddy activity, is an area with strong ocean-atmosphere interactions, especially during winter when cold, dry continental air-masses are passing over. It is also a source region for atmospheric baroclinic instability, shaping the strength and position of the storm tracks that affect the weather and climate in Europe. Simulating accurately the ocean-atmosphere interaction is therefore key for correctly simulating the North Atlantic jet and storm track and European climate (Haarsma et al. 2019; Piazza et al. 2016; Small et al. 2014).

Using the simulations of six GCMs made according to the HighResMIP protocol, we have evaluated the ocean-atmosphere interaction and compared it with reanalyses and available observations. Two mechanisms have been investigated in detail: VMM and PAM. A clear dependence on resolution is found with increased ocean-atmosphere coupling strength and better agreement with reanalysis and observations, for both increasing ocean and atmosphere resolution, in agreement with Czaja et al. (2019). Based on the available simulations we conclude that eddy-permitting ocean resolution (0.4° to 0.25°) and comparable atmosphere resolution are required for a realistic simulation of the two mechanisms. Configurations of eddy-parameterized ocean resolution exhibit lower correlation and coupling coefficients, indicating that the ocean-atmosphere coupling is weaker in the GS region. An intriguing result is that the only model with an eddy-resolving ocean revealed weaker or similar coupling and correlation coefficients compared to its eddy-permitting version. Further analyses is required to better understand this result.

Since most GCMs used to inform the IPCC AR5 and AR6 are eddy-parameterized models, we argue that their representation of the ocean-atmosphere interaction along the GS and other western boundary currents is flawed, with implications for the dynamics of the storm tracks and the associated weather and climate effects. Based on our results, we conjecture that they underestimate both the PAM and the VMM. Because of the limited length of the simulations and the many other drivers that affect the storm tracks when the resolution is changed (drivers that originate in the tropics or in the Arctic), we could not investigate how incorrectly representing ocean-atmosphere interaction affects the storm track dynamics. New, designed on purpose experiments are presently being analyzed.

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1007/s00382-021-06098-9.

**Acknowledgements** We thank two anonymous reviewers for their comments/suggestions, which significantly helped to improve this manuscript. EET, RH, AB, PA, HdV, SD, IEdV, DP, MJR, ESG, CDR
acknowledge PRIMAVERA funding received from the European Commission under Grant Agreement 641727 of the Horizon 2020 research program.

**Funding** This work was supported by Horizon 2020 Framework Programme (641727).

**Availability of data and material (data transparency)** All data used in this study are publicly available on line. The model output is available on the Earth System Grid Federation (ESGF) under the references: HadGEM3-GC3.1 (Roberts 2017b, c, a, 2018), ECMWF-IFS (Roberts et al. 2017a, 2018a, 2017b), EC-Earth3P (EC-Earth 2019, 2018), CNRM-CM6 (Voldoire 2019b, a), MPI-ESM1-2 (von Storch et al. 2017a, b), CMCC-CM2 (Scoccimarro et al. 2017b, a). The observational/reanalysis data are available at: the Copernicus Climate Change Service Climate Data Store (CDS) (DOI:https://doi.org/10.24381/cds.68d2bb30) for ERA5; the Copernicus Marine Environment Monitoring Service (https://resources.marine.copernicus.eu/?option=com_csw&view=details&product_id=WIND_GLO_WIND_L3_REP OBSERVATIONS_012_005) for ASCAT satellite data; http://clime.xp.knmi.nl/select.cgi?sstoiv2_monthly_mean for the OISST.v2.

**Code availability** The computer code used in the present study is available from the corresponding author upon reasonable request.

## Declarations

**Conflicts of interest** The authors declare that they have no conflict of interest.

**Ethics approval** Not applicable.

**Consent to participate** Not applicable.

**Consent for publication** Not applicable.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

## References

Alexander MA, Bladé I, Newman M, Lanzante JR, Lau NC, Scott JD (2002) The atmospheric bridge: The influence of enso teleconnections on air-sea interaction over the global oceans. J Clim 15(16):2205-2231. 10.1175/1520-0442(2002)015%3C2205:TA BTIO%3E2.0.CO;2

Banzon V, Smith TM, Chin TM, Liu C, Hankins W (2016) A long-term record of blended satellite and in situ sea-surface temperature for climate monitoring, modeling and environmental

studies. Earth Syst Sci Data 8(1):165-176. https://doi.org/10.5194/essd-8-165-2016

Bellucci A, Athanasiadis PJ, Scoccimarro E, Ruggieri P, Gualdi S, Fedele G, Haarsma RJ, Garcia-Serrano J, Castrillo M, Putrahasan D, Sanchez-Gomez E, Moine MP, Roberts CD, Roberts MJ, Seddon J, Vidale PL (2021) Air-sea interaction over the gulf stream in an ensemble of highresmip present climate simulations. Clim Dyn. https://doi.org/10.1007/s00382-020-05573-z

Bishop SP, Small RJ, Bryan FO, Tomas RA (2017) Scale dependence of midlatitude air-sea interaction. J Clim 30(20):8207-8221. https://doi.org/10.1175/JCLI-D-17-0159.1

Bryan FO, Hecht MW, Smith RD (2007) Resolution convergence and sensitivity studies with north atlantic circulation models. Part I: the western boundary current system. Ocean Model 16(3-4):141-159. https://doi.org/10.1016/j.ocemod.2006.08.005

Bryan FO, Tomas R, Dennis JM, Chelton DB, Loeb NG, McClean JL (2010) Frontal scale air-sea interaction in high-resolution coupled climate models. J Clim 23(23):6277-6291. https://doi.org/10.1175/2010JCLI3665.1

Cayan DR (1992) Latent and sensible heat flux anomalies over the northern oceans: driving the sea surface temperature. J Phys Oceanogr 22(8):859-881. https://doi.org/10.1175/1520-0485(1992)022%3C0859:LASHFA%3E2.0.CO;2

Chassignenet EP, Xu X (2017) Impact of horizontal resolution (1/12 to 1/50) on gulf stream separation, penetration, and variability. J Phys Oceanogr 47(8):1999-2021. https://doi.org/10.1175/jpo-d-17-0031.1

Chelton DB, Xie SP (2010) Coupled ocean-atmosphere interaction at oceanic mesoscales. Oceanography 23(4):52-69. https://doi.org/10.5670/oceanog.2010.05

Chelton DB, Esbensen SK, Schlax MG, Thum N, Freilich MH, Wentz FJ, Gentemann CL, McPhaden MJ, Schopf PS (2001) Observations of coupling between surface wind stress and sea surface temperature in the eastern tropical pacific. J Clim 14(7):1479-1498. https://doi.org/10.1175/1520-0442(2001)014%3C1479:OOCBSW%3E2.0.CO;2

Chelton DB, Schlax MG, Freilich MH, Milliff RF (2004) Satellite measurements reveal persistent small-scale features in ocean winds. Science 303(5660):978-983. https://doi.org/10.1126/science.1091901

Chelton DB, Schlax MG, Samelson RM (2007) Summertime coupling between sea surface temperature and wind stress in the california current system. J Phys Oceanogr 37(3):495-517. https://doi.org/10.1175/JPO3025.1

Cherchi A, Fogli PG, Lovato T, Peano D, Iovino D, Gualdi S, Masina S, Scoccimarro E, Materia S, Bellucci A et al (2019) Global mean climate and main patterns of variability in the CMCC-cm2 coupled model. J Adv Model Earth Syst 11(1):185-209. https://doi.org/10.1029/2018MS001369

Czaja A, Frankignoul C, Minobe S, Vannière B (2019) Simulating the midlatitude atmospheric circulation: what might we gain from high-resolution modeling of air-sea interactions? Curr Clim Chang Rep 5(4):390-406. https://doi.org/10.1007/s40641-019-00148-5

Deser C, Alexander MA, Xie SP, Phillips AS (2010) Sea surface temperature variability: patterns and mechanisms. Ann Rev Mar Sci 2:115-143. https://doi.org/10.1146/annurev-marine-120408-151453

EC-Earth (2018) Ec-earth-consortium ec-earth3p-hr model output prepared for cmip6 highresmip. https://doi.org/10.22033/ESGF/CMIP6.2323

EC-Earth (2019) Ec-earth-consortium ec-earth3p model output prepared for cmip6 highresmip. https://doi.org/10.22033/ESGF/CMIP6.2322

Eyring V, Bony S, Meehl GA, Senior CA, Stevens B, Stouffer RJ, Taylor KE (2016) Overview of the coupled model intercomparison project phase 6 (cmip6) experimental design and organization.
Geosci Model Dev 9(5):1937–1958. https://doi.org/10.5194/gmd-9-1937-2016

Figa-Saldaña J, Wilson JJ, Attema E, Gelsthorpe R, Drinkwater MR, Stoffelen A (2002) The advanced scatterometer (ASCAT) on the meteorological operational (METOP) platform: a follow on for european wind scatterometers. Can J Remote Sens 28(3):404–412. https://doi.org/10.5589/m02-035

Gutjahr O, Putrasahan D, Lohmann K, Jungclaus JH, von Storch JS, Brüggemann N, Haak H, Stössel A (2019) Max planck insti- tute earth system model (mpi-esm1.2) for the high-resolution model intercomparison project (highresmip). Geosci Model Dev 12(7):3241–3281. https://doi.org/10.5194/gmd-12-3241-2019

Haarsma R, Acosta M, Bakhshi R, Bretonnière PAB, Caron LP, Cas- trillo M, Corti S, Davini P, Exarchou E, Fabiano F, Fladrich U, Fuentes Franco R, García-Serrano J, von Hardenberg J, Koenigk T, Levine X, Meccia V, van Noije T, van den Oord G, Palmeiro FM, Rodrigo M, Ruprich-Robert Y, Le Sager P, Tourigny E, Wang S, van Weele M, Wyser K (2020) Highresmip versions of ec-earth: Ec-earth3p and ec-earth3p-hr. description, model performance, data handling and validation. Geoscientific Model Development Discussions 2020:1–37 https://doi.org/10.5194/gmd-2019-350, https://www.geosci-model-dev-discuss.net/gmd-2019-350/

Haarsma RJ, Roberts MJ, Vidale PL, Senior CA, Bellucci A, Bao Q, Chang P, Corti S, Fučkar NS, Guemas V et al (2016) High resolu- tion model intercomparison project (highresmip v1.0) for cmip6. Geosci Model Dev 9(11):4185–4208. https://doi.org/10.5194/gmd-9-4185-2016

Haarsma RJ, García-Serrano J, Prodhomme C, Bellprat O, Davini P, Drijfhout S (2019) Sensitivity of winter north Atlantic-European climate to resolved atmosphere and ocean dynamics. Sci Rep 9(1):1–8. https://doi.org/10.1038/s41598-019-49865-9

Hersbach H, Bell B, Berrisford P, Hirahara S, Horányi A, Muñoz- Sabater J, Nicolas J, Peubey C, Radu R, Schepers D et al (2020) The era5 global reanalysis. Quart J R Meteorol Soc 146(730):1999–2049. https://doi.org/10.1002/qj.3803

Hsu S (1984) Sea-breeze-like winds across the north wall of the gulf stream: an analytical model. J Geophys Res Oceans 89(C2):2025– 2028. https://doi.org/10.1029/JC089iC02p02025

Kilpatrick T, Schneider N, Qiu B (2014) Boundary layer convergence induced by strong winds across a midlatitude SST front. J Clim 27(4):1698–1718. https://doi.org/10.1175/JCLI-D-13-00101.1

Kirtman BP, Bitz C, Bryan F, Collins W, Dennis J, Hearn N, Kinter JL, Loft R, Rousset C, Siqueira L et al (2012) Impact of ocean model resolution on CCSM climate simulations. Clim Dyn 39(6):1303– 1328. https://doi.org/10.1007/s00382-012-1500-3

Klaver R, Haarsma R, Vidale PL, Hazeleger W (2020) Effective reso- lution in high resolution global atmospheric models for climate studies. Atmos Sci Lett 21(4):e952. https://doi.org/10.1002/asl. 952

Kobashi F, Xie SP, Iwasaka N, Sakamoto TT (2008) Deep atmospheric response to the north pacific oceanic subtropical front in spring. J Clim 21(22):5960–5975. https://doi.org/10.1175/2008JCLI2311.1

Kushnir Y, Robinson W, Bladé I, Hall N, Peng S, Sutton R (2002) Atmospheric GCM response to extratropical SST anomalies: syn- thesis and evaluation. J Clim 15(16):2233–2256. https://doi.org/ 10.1175/1520-0442(2002)015<2233:AGRTES>2.0.CO;2

Lee RW, Woollings TJ, Hoskins BJ, Williams KD, O’Reilly CH, Masato G (2018) Impact of gulf stream SST biases on the global atmospheric circulation. Clim Dyn 51(9):3369–3387. https://doi. org/10.1007/s00382-018-4083-9

Lindzen RS, Nigam S (1987) On the role of sea surface temperature gradients in forcing low-level winds and convergence in the trop- ics. J Atmos Sci 44(17):2418–2436. https://doi.org/10.1175/1520- 0469(1987)044%3C2418:OTROSS%3E2.0.CO;2

Liu WT, Zhang A, Bishop JK (1994) Evaporation and solar irradiance as regulators of sea surface temperature in annual and interannual

changes. J Geophys Res Oceans 99(C6):12623–12637. https://doi. org/10.1029/94JC00604

Ma J, Xu H, Dong C, Lin P, Liu Y (2015) Atmospheric responses to oceanic eddies in the kuroshio extension region. J Geophys Res Atmos 120(13):6313–6330. https://doi.org/10.1002/2014JD0229 30

Maloney ED, Chelton DB (2006) An assessment of the sea surface temperature influence on surface wind stress in numerical weather prediction and climate models. J Clim 19(12):2743–2762. https:// doi.org/10.1175/JCLI3728.1

Minobe S, Kuwano-Yoshida A, Komori N, Xie SP, Small RJ (2008) Influence of the gulf stream on the troposphere. Nature 452(7184):206–209. https://doi.org/10.1038/nature06690

Minobe S, Miyashita M, Kuwano-Yoshida A, Tokinaga H, Xie SP (2010) Atmospheric response to the gulf stream: seasonal vari- ations. J Clim 23(13):3699–3719. https://doi.org/10.1175/2010J CLI3359.1

Moreno-Chamarro E, Caron LP, Loosveldt Tomas S, Gutjahr O, Moine MP, Putrasahan D, Roberts CD, Roberts MJ, Senan R, Terray L et al (2021) Impact of increased resolution on long- standing biases in highresmip-primavera climate models. Geo- sci Model Dev Dis. https://doi.org/10.5194/gmd-2021-209

Moreton SM, Ferreira D, Roberts MJ, Hewitt HT (2021) Air-sea turbulent heat flux feedback over mesoscale eddies. Earth and Space Science Open Archive ESSOAr

Nakamura H, Sampe T, Tanimoto Y, Shimpo A (2004) Observed associations among storm tracks, jet streams and midlatitude oceanic fronts. Earth’s Clim Ocean Atmos Interact Geophys Monogr 147:329–345. https://doi.org/10.1029/147GM18

Nonaka M, Xie SP (2003) Covariances of sea surface tempera- ture and wind over the kuroshio and its extension: evidence for ocean-to-atmosphere feedback. J Clim 16(9):1404–1413. https://doi.org/10.1175/1520-0442(2003)16<1404:COSST A>2.0.CO;2

O’Neill LW, Chelton DB, Esbensen SK (2003) Observations of SST- induced perturbations of the wind stress field over the southern ocean on seasonal timescales. J Clim 16(14):2340–2354. https:// doi.org/10.1175/2780.1

O’Neill LW, Chelton DB, Esbensen SK, Wentz FJ (2005) High-reso- lution satellite measurements of the atmospheric boundary layer response to SST variations along the agulhas return current. J Clim 18(14):2706–2723. https://doi.org/10.1175/JCLI3415.1

O’Neill LW, Chelton DB, Esbensen SK (2010) The effects of SST- induced surface wind speed and direction gradients on midlati- tude surface vorticity and divergence. J Clim 23(2):255–281. https://doi.org/10.1175/2009JCLI2613.1

O’Neill LW, Haack T, Chelton DB, Skyllingstad E (2017) The gulf stream convergence zone in the time-mean winds. J Atmos Sci 74(7):2383–2412. https://doi.org/10.1175/JAS-D-16-0213.1

Parfitt R, Czaja A (2016) On the contribution of synoptic transients to the mean atmospheric state in the gulf stream region. Quart J R Meteorol Soc 142(696):1554–1561. https://doi.org/10.1002/ qj.2689

Parfitt R, Seo H (2018) A new framework for near-surface wind convergence over the kuroshio extension and gulf stream in wintertime: the role of atmospheric fronts. Geophys Res Lett 45(18):9909–9918. https://doi.org/10.1029/2018GL080135

Parfitt R, Czaja A, Kwon YO (2017) The impact of SST resolution change in the era-interim reanalysis on wintertime gulf stream frontal air-sea interaction. Geophys Res Lett 44(7):3246–3254. https://doi.org/10.1002/2017GL073028

Perlin N, De Szoeke SP, Chelton DB, Samelson RM, Skyllingstad ED, O’Neill LW (2014) Modeling the atmospheric boundary layer wind response to mesoscale sea surface temperature per- turbations. Monthly Weather Rev 142(11):4284–4307. https:// doi.org/10.1175/MWR-D-13-00332.1
Pezzi LP, de Souza RB, Acevedo O, Wainer I, Mata MM, Garcia CA, de Camargo R (2009) Multiyear measurements of the oceanic and atmospheric boundary layers at the Brazil-Malvinas confluence region. J Geophys Res Atmos. https://doi.org/10.1029/2008JD011379

Piazza M, Terray L, Boé J, Maisonnave E, Sanchez-Gomez E (2016) Influence of small-scale north atlantic sea surface temperature patterns on the marine boundary layer and free troposphere: A study using the atmospheric arpege model. Clim Dyn 46(5-6):1699-1717. https://doi.org/10.1007/s00382-015-2669-z

Putrasahan DA, Miller AJ, Seo H (2013) Isolating mesoscale coupled ocean-atmosphere interactions in the kuroshio extension region. Dyn Atmos Oceans 63:60-78. https://doi.org/10.1016/j.dynatmoce.2013.04.001

Roberts C, Vitart F, Balmaseda M (in press) Hemispheric impact of north atlantic ssts in sub-seasonal forecasts. Geophys Res Lett https://doi.org/10.1029/2020GL091446

Roberts CD, Senan R, Molteni F, Boussetta S, Keeley S (2017a) Ecmwf ecmwf-ifs-hr model output prepared for cmip6 high-resmip. https://doi.org/10.22033/ESGF/CMIP6.2461

Roberts CD, Senan R, Molteni F, Boussetta S, Keeley S (2017b) Ecmwf ecmwf-ifs-lr model output prepared for cmip6 high-resmip. https://doi.org/10.22033/ESGF/CMIP6.2463

Roberts CD, Senan R, Molteni F, Boussetta S, Keeley S (2018a) Ecmwf ecmwf-ifs-mr model output prepared for cmip6 high-resmip. https://doi.org/10.22033/ESGF/CMIP6.2465

Roberts CD, Senan R, Molteni F, Boussetta S, Mayer M, Keeley SP (2018b) Climate model configurations of the ecmwf integrated forecasting system (ecmwf-ifs cycle 43r1) for highresmips. Geosci Model Dev. https://doi.org/10.5194/gmd-11-3681-2018

Roberts M (2017a) Mohc hadgem3-gc31-hm model output prepared for cmip6 highresmip. https://doi.org/10.22033/ESGF/CMIP6.446

Roberts M (2017b) Mohc hadgem3-gc31-ll model output prepared for cmip6 highresmip. https://doi.org/10.22033/ESGF/CMIP6.1901

Roberts M (2017c) Mohc hadgem3-gc31-mm model output prepared for cmip6 highresmip. https://doi.org/10.22033/ESGF/CMIP6.1902

Roberts M (2018) Mohc hadgem3-gc31-hh model output prepared for cmip6 highresmip. https://doi.org/10.22033/ESGF/CMIP6.445

Roberts MJ, Hewitt HT, Hyder P, Ferreira D, Josey SA, Mizielinski M, Shelly A (2016) Impact of ocean resolution on coupled air-sea fluxes and large-scale climate. Geophys Res Lett 43(19):10-430. https://doi.org/10.1002/2016GL070559

Roberts MJ, Baker A, Blockley EW, Calvert D, Coward A, Hewitt HT, Jackson LC, Kuhlbrodt T, Mathiot P, Roberts CD, Schiemann R, Seddon J, Vannière B, Vidale PL (2019) Description of the resolution hierarchy of the global coupled hadgem3-gc3.1 model as used in cmip6 highresmip experiments. Geosci Model Dev 12(12):4999-5028, 10.5194/gmd-12-4999-2019, https://www.geosci-model-dev.net/12/4999/2019/

Scoccimarro E, Bellucci A, Peano D (2017a) Cmcc cmcc-cm2-hr4 model output prepared for cmip6 highresmip. https://doi.org/10.22033/ESGF/CMIP6.1359

Scoccimarro E, Bellucci A, Peano D (2017b) Cmcc cmcc-cm2-vhr4 model output prepared for cmip6 highresmip. https://doi.org/10.22033/ESGF/CMIP6.1367

Seidov D, Mishonov A, Reagan J, Parsons R (2019) Eddy-resolving in situ ocean climatologies of temperature and salinity in the northwest Atlantic ocean. J Geophys Res Oceans 124(1):41-58. https://doi.org/10.1029/2018JC014548

Small R, Xie S, O'Neill L, Seo H, Song Q, Cornillon P, Spall M, Minobe S et al (2008) Air-sea interaction over ocean fronts and

eddies. Dyn Atmos Oceans 45(3-4):274-319. https://doi.org/10.1016/j.dynatmoce.2008.01.001

Small RJ, Tomas RA, Bryan FO (2014) Storm track response to ocean fronts in a global high-resolution climate model. Clim Dyn 43(3-4):805-828. https://doi.org/10.1007/s00382-013-1980-9

Small RJ, Msadek R, Kwon YO, Booth JF, Zarzycki C (2019) Atmosphere surface storm track response to resolved ocean mesoscale in two sets of global climate model experiments. Clim Dyn 52(3):2067-2089. https://doi.org/10.1007/s00382-018-4237-9

Song Q, Chelton DB, Esbensen SK, Thum N, O'Neill LW (2009) Coupling between sea surface temperature and low-level winds in mesoscale numerical models. J Clim 22(1):146-164. https://doi.org/10.1175/2008JCLI2488.1

Spall MA (2007) Effect of sea surface temperature-wind stress coupling on baroclinic instability in the ocean. J Phys Oceanogr 37(4):1092-1097. https://doi.org/10.1175/JPO3045.1

von Storch JS, Putrasahan D, Lohmann K, Gutjahr O, Jungclaus J, Bittner M, Haak H, Wieners KH, Giorgetta M, Reick C, Esch M, Gayler V, de Vrese P, Raddatz T, Mauritsen T, Behrens J, Brovkin V, Claussen M, Crueger T, Fast I, Fiedler S, Hagemann S, Hohenegger C, Jahns T, Kloster S, Kinne S, Lasslop G, Kornblueh L, Marotzke J, Matei D, Meraner K, Mikolajewicz U, Modali K, Müller W, Nabel J, Notz D, Peters K, Pincus R, Pohlmann H, Pongratz J, Rast S, Schmidt H, Schnur R, Schulzweida U, Six K, Stevens B, Voigt A, Roeckner E (2017a) Mpi-mmpiesm1.2-hr model output prepared for cmip6 highresmip. https://doi.org/10.22033/ESGF/CMIP6.762

von Storch JS, Putrasahan D, Lohmann K, Gutjahr O, Jungclaus J, Bittner M, Haak H, Wieners KH, Giorgetta M, Reick C, Esch M, Gayler V, de Vrese P, Raddatz T, Mauritsen T, Behrens J, Brovkin V, Claussen M, Crueger T, Fast I, Fiedler S, Hagemann S, Hohenegger C, Jahns T, Kloster S, Kinne S, Lasslop G, Kornblueh L, Marotzke J, Matei D, Meraner K, Mikolajewicz U, Modali K, Müller W, Nabel J, Notz D, Peters K, Pincus R, Pohlmann H, Pongratz J, Rast S, Schmidt H, Schnur R, Schulzweida U, Six K, Stevens B, Voigt A, Roeckner E (2017b) Mpi-mmpi-esm1.2-xr model output prepared for cmip6 highresmip. https://doi.org/10.22033/ESGF/CMIP6.10290

Sweet W, Fett R, Kerling J, La Violette P (1981) Air-sea interaction effects in the lower troposphere across the north wall of the gulf stream. Monthly Weather Rev 109(5):1042-1052. 10.1175/1520-0493(1981)109%3C1042:ASIEIT%3E2.0.CO;2

Takatama K, Minobe S, Inatsu M, Small RJ (2012) Diagnostics for near-surface wind convergence/divergence response to the gulf stream in a regional atmospheric model. Atmos Sci Lett 13(1):16-21. https://doi.org/10.1002/asl.355

Tokinaga H, Tanimoto Y, Xie SP (2005) SST-induced surface wind variations over the Brazil-Malvinas confluence: satellite and in situ observations. J Clim 18(17):3470-3482. https://doi.org/10.1175/JCLI3485.1

Voldoire A (2019a) Cnrm-cerfacs cnrm-cm6-1-hr model output prepared for cmip6 highresmip. https://doi.org/10.22033/ESGF/CMIP6.1387

Voldoire A (2019b) Cnrm-cerfacs cnrm-cm6-1 model output prepared for cmip6 highresmip. https://doi.org/10.22033/ESGF/CMIP6.1925

Voldoire A, Saint-Martin D, Sénési S, Decharme B, Alias A, Chevalier M, Colin J, Guérémy JF, Michou M, Moine MP et al (2019) Evaluation of cmip6 deck experiments with cnrm-cm6-1. J Adv in Model Earth Syst 11(7):2177-2213. https://doi.org/10.1029/2019MS001683

Wai MMK, Stage SA (1989) Dynamical analyses of marine atmospheric boundary layer structure near the gulf stream oceanic front. Quart J R Meteorol Soc 115(485):29-44. https://doi.org/10.1002/qj.49711548503
Wallace JM, Mitchell T, Deser C (1989) The influence of sea-surface temperature on surface wind in the eastern equatorial pacific: seasonal and interannual variability. J Clim 2(12):1492–1499. https://doi.org/10.1175/1520-0442(1989)002%3C1492: TIOSST%3E2.0.CO;2

Warner TT, Lakhtakia MN, Doyle JD, Pearson RA (1990) Marine atmospheric boundary layer circulations forced by gulf stream sea surface temperature gradients. Monthly Weather Rev 118(2):309–323. https://doi.org/10.1175/1520-0493(1990)118%3C0309:MABLCF%3E2.0.CO;2

Woollings T, Hoskins B, Blackburn M, Hassell D, Hodges K (2010) Storm track sensitivity to sea surface temperature resolution in a regional atmosphere model. Clim Dyn 35(2–3):341–353. https://doi.org/10.1007/s00382-009-0554-3

Xie SP (2004) Satellite observations of cool ocean–atmosphere inter-action. Bull Am Meteorol Soc 85(2):195–208. https://doi.org/10.1175/BAMS-85-2-195

Xie SP, Hafner J, Tanimoto Y, Liu WT, Tokinaga H, Xu H (2002) Bathymetric effect on the winter sea surface temperature and climate of the yellow and east China seas. Geophys Res Lett 29(24):81–1. https://doi.org/10.1029/2002GL015884

**Publisher's Note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.