Extratropical cyclone induced sea surface temperature anomalies in the 2013/14 winter

. The 2013–2014 winter averaged sea surface temperature (SST) was anomalously cool in the mid-North Atlantic region. This season was also unusually stormy, with extratropical cyclones passing over the mid-North Atlantic every 3 d. However, the processes by which cyclones contribute towards seasonal SST anomalies are not fully quan-tiﬁed. In this paper a cyclone identiﬁcation and tracking method is combined with European Centre for Medium-Range Weather Forecasts (ECMWF) atmosphere and ocean reanalysis ﬁelds to calculate cyclone-relative net surface heat ﬂux anomalies and resulting SST changes. Anomalously large negative heat ﬂux is located behind the cyclones’ cold front, resulting in anomalous cooling up to 0.2 K d − 1 when the cyclones are at maximum intensity. This extratropical-cyclone-induced “cold wake” extends along the cyclones’ cold front but is small compared to climatological variability in the SSTs. To investigate the potential cumulative effect of the passage of multiple cyclone-induced SST cooling in the same location, we calculate Earth-relative net surface heat ﬂux anomalies and resulting SST changes for the 2013– 2014 winter period. Anomalously large winter averaged negative heat ﬂux occurs in a zonally orientated band extending across the North Atlantic


Introduction
The interaction of the ocean and atmosphere has long been recognized as an important element of oceanic cyclogenesis.In the tropics, sub-saturation of air above the sea surface in the vicinity of tropical cyclones results in strong surface heat flux which cools the upper ocean in the wake of tropical cyclones (Kleinschmidt, 1951;Fisher, 1958).In addition, strong winds enhance mixing in the upper ocean, which can result in the transport of cool water to the surface.Similarly, in the midlatitudes it was observed by Pettersen et al. (1962) that surface heat fluxes are largest in the advancing cold air mass behind an extratropical cyclone's cold front.Although fluxes are typically an order of magnitude smaller in extratropical cyclones than tropical cyclones, they still have the potential to influence the underlying ocean.Alexander and Scott (1997) analysed the association of ocean heat fluxes with propagating extratropical cyclones on synoptic timescales.They found increased fluxes directed from the ocean to the atmosphere in the western parts of the Published by Copernicus Publications on behalf of the European Geosciences Union.
cyclones, and the eastern parts of the cyclones were associated with decreased fluxes directed from the ocean to the atmosphere.These results have been confirmed by many subsequent studies (Persson et al., 2005;Nelson et al., 2014;Schemm and Sprenger, 2015;Dacre et al., 2019) and suggest a close association between the cyclones and surface turbulent fluxes in the midlatitudes.However, a cyclone compositing study by Rudeva and Gulev (2011) found that although composites of fluxes show locally very strong positive fluxes in the rear part of the cyclone, the total air-sea turbulent fluxes provided by cyclones were not significantly different from the averaged background fluxes in the North Atlantic, suggesting at least partial cancellation of the flux anomalies associated with cyclones.
The anomalous surface heat fluxes generated by cyclones can create sea surface temperature (SST) anomalies known as the cold wake effect.Case studies of winter cyclones in the north-western Atlantic have found SST cooling in the rear part of cyclones of between 0.4 and 2 K (Ren et al., 2004;Nelson et al., 2014;Kobashi et al., 2019).This is largely due to enhanced turbulent fluxes behind the cold front; however Kobashi et al. (2019) also attribute part of the cooling to cloud shielding of incoming solar radiation, although this is possibly due to the more southerly latitude of the cyclone in their study.Cooling may also result from an episodic wind effect, known as resonant wind-driven mixing, which occurs when the rotation rate of the winds at a fixed point matches the rotation of the wind-driven currents (Crawford and Large, 1996).The magnitude of the cooling has been shown to depend on the cyclone's intensity (Yao et al., 2008;Rudeva and Gulev, 2011), with stronger cyclones creating enhanced surface fluxes and hence increased cooling.In addition, there is some seasonality in the cooling magnitude, with largest cooling occurring in late summer and autumn, when the ocean surface mixed layer in the Northern Hemisphere is shallower (Ren et al., 2004;Kawai and Wada, 2011).Finally, the relationship between enhanced turbulent fluxes is regionally dependent.For example, Tanimoto et al. (2003) showed that in the central North Pacific, enhanced turbulent fluxes can generate local SST variations, but in regions where ocean dynamics are important, such as the western Pacific, the SST anomalies formed in the early winter determine the mid-and late-winter turbulent heat flux anomalies rather than the passage of cyclones.Similarly, Buckley et al. (2015) find that in the Gulf Stream region, ocean dynamics are important in setting the upper-ocean heat content anomalies on interannual timescales and that air-sea heat fluxes damp anomalies created by the ocean.
The effect of extratropical-cyclone-induced fluxes on longer timescale variability has been investigated in both the Atlantic and Pacific.Several studies have shown that wintertime fluxes in the Gulf Stream are characterized by episodic high flux events due to cold-air outbreaks from North America associated with the passage of extratropical cyclones (Zolina and Gulev, 2003;Shaman et al., 2010;Parfitt et al., 2016Parfitt et al., , 2017;;Ogawa and Spengler, 2019).Similar relationships have been found in the high-latitude South Pacific by Papritz et al. (2015).The influence of these enhanced fluxes is to restore the low-level atmospheric baroclinicity destroyed by the passage of the cyclones (Papritz and Spengler, 2015;Vannière et al., 2017), but the impact on the underlying SSTs has not been quantified.
While the role of individual cyclones in local SSTs has been studied, the cumulative effect of cyclone-induced SST changes over individual seasons has not received much attention in the literature.During the winter of 2013-2014 a cold anomaly developed in the SST in the mid-Atlantic.Grist et al. (2016) found that during this winter, enhanced sensible and latent heat fluxes occurred in the North Atlantic, with latent heat fluxes being largest in the eastern North Atlantic and sensible heat flux anomalies stronger in the western North Atlantic, resulting in a reduction in ocean heat content of the subpolar gyre.The extent to which extratropical cyclones were responsible for these enhanced fluxes and associated cooling is not well understood.
In this paper we investigate both the SST cooling associated with individual cyclones and the SST cooling associated with the passage of multiple cyclones over the same location in the 2013-2014 season to determine how significant cyclones were in contributing to the observed cooling.
2 Data and analysis methods

ERA-Interim data
ERA-Interim is a global atmospheric reanalysis dataset (Dee et al., 2011).The data assimilation system used to produce ERA-Interim is based on the Integrated Forecasting System (Cy31r2).The system includes a four-dimensional variational analysis with a 12 h analysis window.The spatial resolution of the dataset is approximately 80 km (T255 spectral) on 60 vertical levels from the surface up to 0.1 hPa; 6 hourly ERA-Interim data have been used in this study to determine extratropical-cyclone-related SST changes.We analyse several reanalysis fields from ERA-Interim which are described in this section.
The net surface thermal radiation (Q LW ) is the thermal radiation emitted by the atmosphere and clouds reaching the surface minus the amount emitted by the surface.Surface solar radiation (Q SW ) is the amount of solar radiation reaching the surface (both direct and diffuse) minus the amount reflected by the surface.Surface latent heat flux (Q E ) is the exchange of latent heat with the surface through turbulent diffusion, and the surface sensible heat flux (Q H ) is the exchange of sensible heat with the surface.The magnitudes of Q E and Q H depend on the wind speed and moisture and temperature differences between the surface and the lower atmosphere.
Weather Clim.Dynam., 1, 27-44, 2020 www.weather-clim-dynam.net/1/27/2020/ The net surface heat flux (Q N ) is given by the sum of Q SW , Q LW , Q H and Q E .The European Centre for Medium-Range Weather Forecasts (ECMWF) convention for vertical fluxes is positive downwards.We also analyse ERA-Interim SSTs, which are the temperatures of sea water near the surface.ERA-Interim SSTs are taken from different sources, depending on the dates of the following reanalyses: NCEP 2D-Var (two-dimensional variational analysis) SST (January 1989-June 2001), NOAA Optimum Interpolation SST V2 (July-December 2001), NCEP Real-Time Global SST (January 2002-January 2009) and Met Office Operational SST (February 2009(February -2015)).

ECMWF Ocean Reanalysis System (ORAS5)
The ECMWF Ocean Reanalysis System 5 (ORAS5) is a global eddy-permitting ocean-sea ice reanalysis system.It provides an estimate of the historical ocean state from 1979 to present.The ocean model resolution in ORAS5 is 0.25 • in the horizontal (approximately 25 km in the tropics and increasing to 9 km in the Arctic) and 75 levels in the vertical.ORAS5 uses the Nucleus for European Modelling of the Ocean (NEMO v3.4.1) ocean coupled to the Louvainla-Neuve Sea Ice Model (LIM2).It includes a prognostic thermodynamic-dynamic sea-ice model with assimilation of sea-ice concentration data (Zuo et al., 2017(Zuo et al., , 2019)).The ocean mixed layer is the layer immediately below the ocean air-sea interface and is typically tens of metres deep in summer, while values of several hundred metres may be reached in winter.In this paper the interannually and monthly varying mixed layer depth (MLD) from ORAS5 is used to calculate the SST tendencies ( SST) due to surface heat flux between 1989 and 2016.The ORAS5 MLD is the first depth at which the density difference, compared to density at 10 m depth, reaches 0.01 kg m −3 .The ECMWF Ocean Reanalysis MLD agrees well with the observationally based MLD estimates in the mid-Atlantic (Toyoda et al., 2017).However, in deep convective regions, such as the Labrador Sea, the density difference MLD definition can overestimate MLD (Courtois et al., 2017).

Cyclone identification
Following Dacre et al. (2012) we identify and track the position of the 200 most intense cyclones in 20 years of the ERA-Interim dataset  during wintertime only (DJF) using the tracking algorithm of Hodges (1995).Tracks are identified using 6-hourly 850 hPa relative vorticity, truncated to T42 resolution to emphasize the synoptic scales.The 850 hPa relative vorticity features are filtered to remove stationary or short-lived features that are not associated with extratropical cyclones.The intensity of the cyclones is measured by the maximum T42 vorticity.The 200 most intense DJF cyclone tracks with maximum intensity in the North Atlantic (30-90 • N, 70-10 • W) are used in this study, and they account for 19 % of the total number of identified North Atlantic cyclones during this period.These tracks are shown in Fig. 1a.The cyclones generally propagate in a northeasterly direction, from the East Coast of the US towards Iceland.The position of the cyclones 24 h before maximum intensity (max − 24) are predominantly over the Gulf Stream (Fig. 1b).By maximum intensity (max) the majority of cyclones are located east of Newfoundland (Fig. 1c).During the decaying stage of the cyclones' evolution (max + 24) the cyclones are more uniformly distributed across the North Atlantic (Fig. 1d).

Cyclone-relative composites
The fields, described in Sect.2.1, are extracted from ERA-Interim at each of the 6 hourly locations of the cyclone within a 30 • radius surrounding the cyclone centre.Cyclone-relative composites are produced by averaging over all cyclones.Following Catto et al. (2010), before compositing, the fields are rotated according to the direction of travel of each cyclone such that the direction of travel in the composite becomes the same for all cyclones.Since the cyclones have quite different propagation directions, performing the rotation ensures that mesoscale features such as warm and cold fronts are approximately aligned and are not smoothed out by the compositing.As this method assumes that the cyclones all intensify and decay at the same rate, only the 200 most intense cyclones are included in the composite.The 850 hPa relative vorticity mean and standard deviation of the of these cyclone are shown in (Dacre et al., 2012;their Fig. 2c).Limiting the number of cyclones produces a more homogeneous group in terms of their evolution but will bias the mean fields towards being typical of the most intense cyclones.Data are only included in the composite where grid points lie over the ocean surface.

Cyclone and environmental flow partition
To understand the part played by cyclones in the development of the cold anomaly during the winter of 2013-2014 the net surface heat flux is partitioned into environmental and cyclone components.To partition Q N into a part due to environmental flow and a part associated with cyclones we combine the cyclone tracks for that season with a masking method.A cyclone mask is calculated for each time step where the regions influenced by a cyclone are given a value of 1 (i.e. they are inside the cyclone mask), and regions that are not influenced are given a value of zero (i.e. they are outside the cyclone mask; Hawcroft et al., 2012;Sinclair and Dacre, 2019).The heat flux associated with cyclones is asymmetric around the cyclone.Large negative flux is found close to the cyclone centre and extending along the trailing cold front, which lies to the west of the cyclone centre.To account for this, a cyclone mask at a given time, t, is created by identifying the position of a cyclone and also its position during the prewww.weather-clim-dynam.net/1/27/2020/ Weather Clim.Dynam., 1, 27-44, 2020 vious 30 h representing the cyclone's cold wake.A mask is created which extends 14 • in a radial circle from each track point, creating an elongated oval shape which encompasses both the cyclone centre and the cold wake (Fig. 2).This is equivalent to taking a swath with circles around all cyclone positions in the previous 30 h.The Q N anomaly when cyclones are not present is associated with the environmental flow.The Q N anomaly within the cyclone masks is associated with cyclones embedded within the environmental flow.

Sea surface temperature tendency
In order to calculate the SST tendency ( SST) we must take into account the MLD through which the heating or cooling due to the surface fluxes is mixed (h).Mixing in the ocean is assumed to occur between the surface and the MLD obtained from ORAS5.Assuming a well-mixed layer, the SST tendency due to Q N and MLD ( SST Q N ) is given by where ρ is the density of sea water, 1024 kg m −3 , and c p is the specific heat capacity of sea water, 4000 J kg −1 K −1 .
The SST tendency anomaly due to air-sea interactions (ASIs), SST ASI , for the 2013-2014 season is determined by subtracting the climatological SST tendency from the 2013-2014 SST tendency and is given by where the overbar represents the 1989-2015 climatological value.The SST tendency anomaly can be separated into the anomaly associated with (i) anomalous Q N (term 1 in Eq. 2), Weather Clim.Dynam., 1, 27-44, 2020 www.weather-clim-dynam.net/1/27/2020/ (ii) anomalous MLD (term 2 in Eq. 2) and (iii) anomalous entrainment through the base of the mixed layer, Q ENT (term 3 in Eq. 2).Since we have no measurements of the entrainment flux anomaly across the ocean boundary layer, it is estimated to be 20 % of the surface Q N anomaly (Stull, 1988).This estimate for the entrainment flux neglects contributions made by wind-driven turbulence.

North Atlantic heat flux and SST tendency climatologies
Figure 3 shows the average heat flux for the period DJF 1989-2015 over the North Atlantic.The net surface solar radiation is positive, with a meridional gradient of 50 W m −2 1000 km −1 , (Fig. 3a) and the net thermal radiation (Fig. 3b) is negative, with a magnitude between −50 and −100 W m −2 .The sensible heat flux (Fig. 3c) is generally positive over land and negative over the ocean, with negative flux between −50 and −150 W m −2 in the Gulf Stream and Davis Strait regions, caused by the advection of cold air from the land over relatively warm oceans in DJF.The latent heat flux (Fig. 3d) is generally negative, with a band of enhanced negative flux extending in a north-eastwards direction from the East Coast of the US towards Iceland, with the values > 200 W m −2 found in the west of the North Atlantic.The net heat flux, Q N (Fig. 3e), is negative over the majority of the domain.A combination of Q H and Q E results in maximum negative heat flux > 300 W m −2 in the Gulf Stream re-gion.The largest Q N values are co-located with the position of the North Atlantic storm track (Fig. 3f), which also exhibits a pronounced south-west-north-east tilt similar to the storm track in Fig. 1a. Figure 4a and b show the climatological SST and the average SST change over DJF ( SST TOT ) for the period 1989-2015.During DJF the North Atlantic cools by an average of 2 K (Fig. 4b).The cooling is greatest where the SST gradient is largest in the Gulf Stream region (dotted box in Fig. 4a), with cooling of over 6 K over the winter.Unlike the climatological Q N shown in Fig. 3e the region of highest SST TOT does not extend over the North Atlantic towards Iceland but remains close to the East Coast of the US.
In order to calculate the SST tendency due to Q N we must take into account the MLD through which heating or cooling at the surface is mixed using Eq. ( 1). Figure 4c shows the climatological MLD.The MLD ranges from a few tens of metres close to Newfoundland to over 1000 m in the Labrador and Norwegian seas, where deep convection occurs.On average the MLD deepens by 50 % between December and February outside the deep convection regions (not shown).Taking the MLD into account restricts the largest SST Q N to the western North Atlantic region (Fig. 4d).The difference between SST TOT and SST Q N is shown in Fig. 4e.Differences > 6 K are found close to the coast, where the western boundary currents transport warm waters north, resulting in reduced cooling over the DJF period than would be expected due to Q N alone.

Cyclone-relative heat flux composites
The north-south SST gradient in Fig. 4a is important for creating negative heat flux when cold dry air, from higher latitudes, is advected over relatively warm ocean surfaces.The largest heat flux occurs when large differences in temperature or moisture between the surface and overlying atmosphere are co-located with enhanced wind speeds.Given the spatial distribution of SSTs, atmospheric temperatures and moisture content, this is likely to occur when there are anomalously strong meridional winds such as those found ahead of and behind extratropical cyclones.Thus, whilst the maximum climatological SST gradient controls where locally high surface flux occurs, the presence of extratropical cyclones is likely to control when the largest heat flux occurs and its magnitude (Rudeva and Gulev, 2011;Ogawa and Spengler, 2019).In addition, the similarities between the spatial pattern of climatological Q N (Fig. 3e) and cyclone numbers per season (Fig. 3f) suggest that cyclones play a role in generating Q N .In this section we investigate Q N in a cyclone-relative frame of reference using the methodology described in Sect.2.3 and 2.4.
Figure 5 shows composite cyclone-centred heat flux for the 200 most intense cyclones occurring between 1989 and 2009.The cyclones are at maximum intensity (maximum 850 hPa relative vorticity) when they are located at the centre www.weather-clim-dynam.net/1/27/2020/ Weather Clim.Dynam., 1, 27-44, 2020  of the domain and were rotated so that they all travel from left to right.Q SW (Fig. 5a) contribution to Q N is positive, and, like the Earth-relative perspective (Fig. 3a), there is a weak gradient.However, because the cyclones typically travel in a north-eastwards direction and they have been rotated, this gradient is also rotated.Q LW (Fig. 5b) is negative everywhere and small, similar to the Earth-relative perspective.Q LW is slightly enhanced in the region behind the cold front, potentially due to a reduction of cloud, which reduces the downwelling radiation.Q H (Fig. 5c) shows a dipole structure at www.weather-clim-dynam.net/1/27/2020/ Weather Clim.Dynam., 1, 27-44, 2020 this stage in the cyclones' evolution, with negative flux behind the cyclone centre and positive flux in the cyclones' warm sector.Q E (Fig. 5d) is negative everywhere, has its minimum in an extended region behind the cyclone and is reduced ahead of the cyclone.Q N (Fig. 5e) is therefore negative surrounding the cyclone centre, with a minimum behind the cyclone.Negative flux > 200 W m −2 occurs within 1000 km of the cyclone centre but extends almost 2000 km behind the cyclone due to a combination of Q H and Q E occurring behind the cold front.
In order to illustrate the processes leading to negative Q H behind the cyclone cold front and positive Q H in the warm sector it is necessary to examine the temperature characteristics of the different air masses in these regions.Fig-ure 6a shows the composite 10 m air temperature overlaid with 925 hPa winds, and Fig. 6b shows the SST.The 10 m air temperature exhibits a wave-like structure, whilst the SST gradient is more linear.This results in a large negative near surface temperature difference (SST > 10 m temperature) behind the cold front.Cyclonic winds advect relatively cold air over a warm ocean surface behind the cyclone, resulting in negative Q H . Ahead of the cyclone cyclonic winds advect warm air over the ocean surface; 6 h before maximum intensity, SST < 10 m air temperature ahead of the cyclone (not shown), but at maximum intensity the difference is close to zero (white contour in Fig. 6c shows −0.4 K).Since Q H is a 6 h average, this results in the positive Q H observed in Fig. 5c. Figure 6d  the saturation specific humidity at the SST respectively.The 10 m specific humidity also has a pronounced wave structure, with drier air behind the cold front and moister air in the warm sector.Behind the cold front the saturation deficit is > 4 g kg −1 (Fig. 6f), causing evaporation of moisture from the surface and large negative Q E observed in Fig. 5d.Ahead of the cyclone the moisture deficit is < 1 g kg −1 , significantly reducing the magnitude of Q E .Figure 7a-c show composite Q N centred on cyclones at different stages in their evolution.As the cyclones start to intensify negative Q N behind the cold front strengthens (Fig. 7b).During the mature stages of the cyclone evolution (Fig. 7c) Q N begins to decrease and to wrap cyclonically around the cyclone centre.This is due to the fact that the cold front typically rotates cyclonically towards the warm front as the cyclone reaches maturity.These results are consistent with the findings of Rudeva and Gulev (2011), who found that turbulent heat flux increases with cyclone intensity.
Interestingly, throughout the cyclone life cycle there is a secondary minimum in Q N occurring approximately 2000 km ahead of the cyclone location.This secondary minimum does not change magnitude, so it is unlikely to be affected by the cyclones at the centre of the composite.It is possible that this second minimum is due the presence of a downstream cyclone, indicated by the composite mean sealevel pressure (mslp) contours, which extends towards the upper-right quadrant of the domain.
Since many of the 200 cyclones contributing to the composite Q N (Fig. 5e) are generated over the Gulf Stream region it is possible that large negative Q N occurring behind of the cyclone centre could be an artefact of their preferential cyclogenesis over a region of climatologically large negwww.weather-clim-dynam.net/1/27/2020/ Weather Clim.Dynam., 1, 27-44, 2020 At all stages of the cyclone evolution negative Q N behind the cyclone centre is more than 0.5 SD greater than the mean for strong cyclones and is more than 1 SD greater than the mean at maximum intensity (Fig. 7e).Ahead of the cyclone Q N is 0.4-0.8SD greater than the mean at maximum intensity (Fig. 7e) due to warm air advection in the warm sector of the cyclone.Note that towards the edges of the domain many of the grid points are over land, so they have been excluded; therefore the sample size contributing to the composites is small, resulting in a noisy field.

Cyclone-relative SST tendency
Using Eq. ( 1) the evolution of the cyclone-relative SST tendencies due to Q N can be estimated.Figure 8 shows Weather Clim.Dynam., 1, 27-44, 2020 www.weather-clim-dynam.net/1/27/2020/ SST Q N per day for cyclones at different stages in their evolution.The patterns of SST Q N are similar to the patterns of Q N (Fig. 7a-c), showing that the surface flux is the dominant variable in the cooling and that variations in mixed layer depth are less important.As for Q N , SST Q N increases as the cyclones reach maximum intensity, with maximum SST tendencies of 0.2 K d −1 occurring in the cold sector behind the cold front.The normalized SST Q N anomalies are calculated by subtracting the climatological SST TOT and dividing by the standard deviation of SST TOT (Fig. 8d-f).
SST Q N is larger than the climatological mean behind the cold front.However in this case the intense cyclones only reduce SST by up to 0.25 times the standard deviation.SST Q N cooling associated with each individual extratropical cyclone is of the order of 0.1-0.2K d −1 (Fig. 8a-c); therefore if many cyclones track over the same location we might expect to see a signature of the storm track in the seasonal Q N and SST anomaly patterns.Figure 9a shows the tracks of cyclones in the North Atlantic region.Applying the cyclone-masking methodology described in Sect.2.6 to the 2013-2014 winter cyclones, we see that the cyclones track in a south-west-north-east direction in a narrow band that extends from the East Coast of the US towards the UK (Fig. 9b).This season was unusually stormy in the UK, with cyclones passing over the UK every 3 d (Priestley et al., www.weather-clim-dynam.net/1/27/2020/ Weather Clim.Dynam., 1, 27-44, 2020 Weather Clim.Dynam., 1, 27-44, 2020 www.weather-clim-dynam.net/1/27/2020/ 2017).The effect of propagation speed is taken into account in the masking methodology, since multiple time steps for a single cyclone contribute to the seasonal mask fraction.
We also apply the cyclone-masking methodology to winter cyclones between 1989 and 2015 (Fig. 9c) and find that, in comparison to the 1989-2015 average, in 2013-2014 the storm track was more active, with a higher number of cyclones, and also more zonal than usual (Fig. 9d).This was associated with an anomalously strong and zonally elongated upper-level westerly jet in 2013-2014 (Kendon and McCarthy, 2015).More cyclones travelled towards western Europe and Iceland than usual.It was shown in Fig. 7 that intensifying cyclones create negative Q N behind the cold front, where cold dry air is advected over a warm ocean surface.Therefore we hypothesize that an anomalously strong and zonal storm track will result in an anomalously strong and zonally orientated seasonal Q N anomaly.The 2013-2014 DJF season Q N is shown in Fig. 9e, and the 2013-2014 Q N anomaly is shown in Fig. 9f.The seasonal Q N anomaly has a tripole pattern, with stronger negative Q N in the mid-North Atlantic and weaker negative Q N over the Gulf Stream and in the Norwegian and Greenland seas.In the mid-Atlantic this is consistent with the shift in the storm track, with cyclones travelling more zonally from the US towards western Europe rather than northeastwards towards Iceland.
Figure 10a and b show UK Met Office synoptic analysis charts at 00:00 UTC on 24 and 20 December 2013 respectively.At both times there is a low-pressure centre situated to the west of Scotland and long trailing cold fronts extending across the North Atlantic.Figure 10c and d show Q N at the corresponding times, and Fig. 10e and f show the cyclone masks.The red lines show the full track of the cyclones, and the elongated oval masks the location of the cyclones at the analysis time and 30 h earlier.The mask captures Q N surrounding the cyclone centres and the cyclone's cold wake.Sensitivity tests have been performed using a radial circle ranging from 12 to 16 • and masks extending 24 to 36 h prior to the analysis time.

2013-2014 SST anomalies
Figure 11a shows the 2013-2014 DJF SST tendency anomaly ( SST ) that is associated with anomalous Q N .As expected SST due to anomalous Q N closely resembles the Q N anomaly (shown in Fig. 9f), with anomalous cooling in the mid-North Atlantic, where the flux anomaly is negative and anomalous warming (less cooling than climatology) in the Gulf Stream and Norwegian Sea regions.Small differences are due to spatial inhomogeneity in the North Atlantic climatological MLD. Figure 11b shows the 2013-2014 DJF SST that is associated with anomalous MLD.The 2013-2014 MLD is shallower than the climatological average over much of the domain, particularly near the Gulf Stream region, and deeper than climatology in the mid-Atlantic region.
In the mid-North Atlantic the enhanced negative Q N results in negative buoyancy and thus mixing, deepening the MLD.Therefore, the surface flux decreases the temperature over a deeper layer of the ocean than usual, which reduces the direct SST cooling due to Q N .At the same time, the increased MLD mixes colder water from below into the mixed layer, which cools the surface indirectly.Figure 11c shows the sum of the SST due to anomalous Q N , anomalous MLD and anomalous entrainment ( SST ASI ).It shows the same tripole pattern as the SST TOT anomaly (Fig. 11d), which has an average SST cooling anomaly of −1.0 K in the mid-North Atlantic region (black box in Fig. 11d).In the mid-North Atlantic region, the SST ASI accounts for 68 % of the total anomalous SST cooling in the mid-North Atlantic.The largest discrepancies occur along the eastern coast of North America, suggesting that ocean dynamics are responsible for transporting warmer and colder water into these regions via the western boundary currents.For example, off the coast of Nova Scotia the SSTs are very cold, since the Labrador current flows south from the Arctic Ocean.Advection of relatively warm air over these cold SSTs in a region where the mixed layer depth is shallow results in a positive SST tendency anomaly.

2013-2014 cyclone and environmental flow SST anomalies
In this section we partition the Q N anomaly into a part associated with the environmental flow (i.e.outside the cyclone masks described in Sect.2.6) and a part associated with the presence of cyclones (inside the cyclone masks described in Sect.2.6). Figure 12a shows the contribution to the DJF Q N anomaly due to both cyclones and the environmental flow, and Fig. 12b shows the contribution that is due to the environmental flow only.Both Fig. 12a and b show a tripole pattern, with anomalously negative heat flux between 40 and 60 • N and anomalously positive heat flux over the Gulf Stream and Norwegian seas.This suggests that the overall pattern is controlled by the environmental flow.
The anomalously positive heat flux in the Norwegian Sea and Gulf Stream region may be due to a reduced number of cold-air outbreaks (Papritz and Spengler, 2017;Papritz and Grams, 2018;Parfitt et al., 2016).In the mid-North Atlantic the negative Q N anomaly is enhanced, compared to contribution made by the environment (doubled from 32 % to 68 % of the total seasonal Q N anomaly), when cyclones are present.
In the Norwegian Sea and Mediterranean Sea the positive Q N anomaly due to the environmental flow is suppressed by the presence of cyclones.Thus cyclones embedded within the environmental flow pattern tend to increase the negative surface heat flux, consistent with the cyclone-relative results in Sect.3.2.Varying the size of the cyclone mask radius from 12 to 16 • results in a contribution to the total Q N anomaly when cyclones are present in the mid-North Atlantic that ranges from 62 % to 72 %, and varying the length of the cyclone www.weather-clim-dynam.net/1/27/2020/ Weather Clim.Dynam., 1, 27-44, 2020 mask from 24 to 36 h results in a cyclone contribution from 62 % to 71 %.This suggests that the conclusion that cyclones enhance the Q N anomaly does not depend strongly on the choices used to define the cyclone mask.
Figure 12c shows the contribution to the 2013-2014 SST ASI due to both cyclones and the environmental flow, and Fig. 12d shows the contribution that is due to the environmental flow only.As for the Q N anomaly, both show a tripole pattern, suggesting that the environmental flow controls the anomalously large cooling in the mid-North Atlantic and the anomalously small cooling in the Norwegian Sea and Gulf Stream regions.In the mid-North Atlantic region, the negative SST ASI is also enhanced (from 28 % to 41 % of the total seasonal SST ) when cyclones are present.This is less than their added contribution to the 2013-2014 Q N anomaly because the enhanced negative Q N cools a deeper layer of the ocean than usual, which reduces the direct SST cooling.

Conclusions
In this paper we investigate both the SST cooling associated with individual cyclones and the SST cooling associated with the passage of multiple cyclones over the same location in the 2013-2014 season to determine how significant cyclones were in contributing to the anomalously large cooling that occurred during the 2013-2014 winter.We find that enhanced air-sea exchange of heat and moisture in the cold sector behind the cold front of an extratropical cyclone can lead to cooling of up to 0.2 K d −1 for the strongest cyclones.This cooling is relatively small compared to the variability in SST tendency in the North Atlantic; thus the cold wake associated with the passage of an individual extratropical cyclone is difficult to observe in the instantaneous SST field.
During the 2013-2014 DJF season there was a zonal band of anomalously large negative heat fluxes extending from the East Coast of the US towards Europe.The mixed layer depth was also anomalously deep in the mid-North Atlantic due to enhanced mixing and entrainment of water into the mixed layer from below.The combination of these air-sea interacwww.weather-clim-dynam.net/1/27/2020/ Weather Clim.Dynam., 1, 27-44, 2020 tions accounts for 68 % of the total SST anomaly.Thus air-sea interactions were very important in determining the anomalous SST cooling between December 2013 and February 2014.
In 2013-2014 the environmental flow pattern was anomalously zonal compared to climatology over the North Atlantic.This resulted in anomalously negative heat flux between 40 and 60 • N and anomalously positive heat flux over the Norwegian and Gulf Stream regions.When cyclones were present, heat flux from the ocean to the atmosphere was doubled in the mid-North Atlantic region.The enhancement of the negative heat flux caused a direct cooling of the ocean but also led to increased entrainment and thus a deeper mixed layer.This deepening reduces the overall effect of the cooling, as the heat loss acts on a greater volume of water than normal.Consequently, while the SST tendency anomaly in the mid-North Atlantic was enhanced by the presence of cyclones, it was by a smaller amount than might be expected due to a doubling of the heat flux.We conclude that both the environmental flow and extratropical cyclones embedded within this flow played important roles in determining the extreme 2013-2014 winter season SST cooling.

Figure 1 .
Figure 1.(a) Tracks of the 200 most intense DJF North Atlantic storms between December 1989 and February 2009 (black).Position of the cyclones (b) 24 h before maximum intensity (max − 24), (c) at maximum intensity (max) and (d) 24 h after maximum intensity (max + 24; red crosses).Overlaid on the DJF North Atlantic 1989-2015 net surface heat flux climatology, Q N (W m −2 ).

Figure 2 .
Figure 2. Schematic of cyclone-masking method.Red and blue lines show the approximate position of the cyclone warm and cold front respectively.Red crosses show the position of a given cyclone at time (t) and at 6-hourly intervals up until 30 h previously (−30).The outer circle shows the 30 • radius circle used to produce the composites in Figs. 5 to 8. The smaller dashed circles show example locations of a 14 • radius circle centred on the cyclone position at t = 0 and t = 30.The grey shading shows the extent of the cyclone mask, encompassing the cyclone's centre and cold wake.

Figure 3 .
Figure 3. North Atlantic 1989-2015 heat flux climatologies (W m −2 ).(a) Surface solar radiation Q SW , (b) surface thermal radiation Q LW , (c) surface sensible heat flux Q H , (d) surface latent heat flux Q E and (e) net surface heat flux Q N .Note change in scale in (e).Positive flux is into the surface, and negative flux is into the atmosphere.(f) 1989-2015 climatological number of cyclones per DJF season within 12 • of grid point.

Figure 4 .
Figure 4. North Atlantic DJF 1989-2015 (a) sea surface temperature (SST; K), black dotted box outlines region of the Gulf Stream, (b) SST TOT (K per DJF), (c) mixed layer depth (m), (d) SST Q N (K per DJF), and (e) difference between SST TOT and SST Q N .

Figure 5 .
Figure 5. Composite cyclone-centred heat fluxes (W m −2 ) for cyclones at maximum intensity.(a) Surface solar radiation Q SW , (b) surface thermal radiation Q LW , (c) surface sensible heat flux Q H , (d) surface latent heat flux Q E and (e) net surface heat flux Q N overlaid with mslp.Note the change in colour scale in (e).Positive fluxes are into the surface, and negative fluxes are into the atmosphere.The radius of the circles is 3000 km, and the cyclones travel from left to right.The blue and red lines in (e) represent the approximate positions of the cold and warm fronts at max intensity.

Figure 7 .
Figure 7. Evolution of cyclone-centred (a-c) net surface heat flux Q N (filled contours, W m −2 ), (d-f) normalized net heat flux anomaly (filled contours) and mslp (black contours; hPa).Cyclones reach the centre of the domain (a, d) 24 h before maximum intensity, (b, e) at maximum intensity and (c, f) 24 h after maximum intensity.In (d)-(f) negative normalized heat flux anomalies indicate anomalously large heat flux into the atmosphere, and positive anomalies indicate anomalously small heat flux into the atmosphere, compared to climatology.The blue and red lines in (e) represent the approximate positions of the cold and warm fronts at max intensity.

Figure 8 .
Figure 8. Evolution of cyclone-centred (a-c) SST Q N (filled contours; K 12 h −1 ), (d-f) normalized SST Q N anomaly (filled contours) and mslp (black contours; hPa).Cyclones reach the centre of the domain (a, d) 24 h before maximum intensity, (b, e) at maximum intensity and (c, f) 24 h after maximum intensity.In (d)-(f) negative normalized SST Q N anomalies indicate anomalously large SST cooling due to Q N , and positive anomalies indicate anomalously small SST cooling, compared to climatology.

Figure 10 .
Figure 10.(a, b) UK Met Office synoptic analysis charts, (c, d) Q N (W m −2 ), and (e, f) cyclone mask overlaid with tracks for (a, c, e) 00:00 UTC 24 December 2013 and (b, d, f) 00:00 UTC 20 December 2013.Red crosses show the position of cyclones at the analysis time and 30 h previously.

Figure 12 .
Figure 12. 2013-2014 anomalous Q N associated with (a) both cyclones and the environmental flow and (b) the environmental flow.2013-2014 anomalous SST due to anomalous Q N , MLD and entrainment associated with (c) both cyclones and the environmental flow and (d) the environmental flow.