more than 90% of the increasing radiative forcing imbalance is reflected as global ocean warming [ alexander et al . ,
the deceleration of the upper ocean heat storage since the mid2000s [ levitus et al .
cheng and zhu , 2015 ] has instigated an active search for the  missing heat  in the deep ocean [ meehl et al .
, 2011 ; trenberth and fasullo , 2010 ; katsman and van oldenborgh , 2011 ; balmaseda et al . , 2013 ; kosaka and xie , 2013 ] . the enhanced downward heat transfer necessary for this vertical heat redistribution has been proposed to be related to the hiatus in global warming , although the hiatus started earlier and lasted for the 15  years period 19982013 [ trenberth et al .
two natural mechanisms have been proposed to explain the deceleration of the rise in global mean surface temperature and the efficient transfer of heat to the deep ocean : ( i ) the intensification of zonal winds causing la nialike conditions in the tropical pacific [ meehl et al . , 2011 ;
trenberth and fasullo , 2010 ; katsman and van oldenborgh , 2011 ; balmaseda et al . , 2013 ; kosaka and xie , 2013 ] and ( ii ) decadal variations in the strength of deep convection and currents in the north atlantic [ meehl et al . , 2011 ;
katsman and van oldenborgh , 2011 ; chen and tung , 2014 ; drijfhout et al . ,
modeling simulations show that la nia surface cooling in the eastern and central tropical pacific can decrease the rate of increase in atmospheric temperature ( although the similarity of observed and simulated deceleration is still debated [ trenberth , 2015 ; karl et al . ,
2015 ] ) and induce an anomalous heat flux into the ocean [ kosaka and xie , 2013 ; england et al . , 2014 ; balmaseda et al . ,
associated to this anomalous heat flux into the ocean , subsurface heat uptake in the western tropical pacific and indian oceans occurs [ goddard , 2014 ; england et al . , 2014 ; meehl et al . , 2013 ; lee et al . ,
2015 ] mostly above 300  m ( see figure s1 in the supporting information ) [ chen and tung , 2014 , figures 1 and 3 ; nieves et al . ,
2015 ] , thus not largely involving transfer of heat into the deep ocean . regarding the second mechanism , vertical heat inventories of the north atlantic show notable heat transfer from the upper ( 0300  m ) to intermediate ( 300700  m ) and deep ( > 700  m ) layers since the mid2000s ( see figure s1 and corresponding comments ) [ alexander et al . , 2013 , figure 3.1 ; chen and tung , 2014 , figure 6 ; williams et al .
such recent downward heat transfer follows decades of sustained surface warming and increasing ocean heat content in the basin [ drijfhout et al . , 2014 ;
however , concomitant with the downward heat export after the mid2000s , the increase in ocean heat content in the upper 2000  m of the north atlantic paused [ chen and tung , 2014 , figures 1 and 3 ; hkkinen et al .
, 2015 , figure 1 ] , making the bulk rate of change of ocean heat content per unit area in the basin almost zero .
bulk rates of change in ocean heat content per unit area are considered an estimation of the ocean heat uptake ( q , w  m ) at ocean basin scales [ e.g. , mauritzen et al .
, 2012 ; guemas et al . , 2013 ; chen and tung , 2014 ] .
changes in ocean circulation can generate ocean heat content anomalies of opposite signs between different ocean layers and regions ( e.g. , subtropical versus subpolar gyre ) , but the anomalies compensate when integrated at ocean basin scales and decadal timescales [ lozier et al . , 2008 ; mauritzen et al . ,
therefore , one can conclude that the heat redistribution in the north atlantic since the mid2000s as a whole has apparently occurred without extra heat uptake from the atmosphere .
overall , during the past decade , different basins have played contrasting roles on the earth 's atmosphereocean heat budget distribution . in particular , the north atlantic would not have contributed notably as in previous decades to atmospheric heat uptake , in contrast to that observed in the tropical pacific and southern oceans [ roemmich et al . ,
2015 ] , but to the quick transfer of heat accumulated during the previous decades in the upper layers to the deep ocean . what are the causes for this shift ?
deep convection in the subpolar north atlantic under highsalinity conditions might account for the quick transfer of heat to the deep ocean [ chen and tung , 2014 ] .
according to these authors , recent highsalinity conditions may be caused by a strengthening of the atlantic meridional overturning circulation ( amoc ) .
agreeing that high salinity in deep convection regions favors heat transfer to the deep ocean [ mauritzen et al . ,
first , available observational records do not show a strengthening of the amoc [ kanzow et al .
second , deep convection in the subpolar north atlantic is restricted to the labrador and irminger seas ( see figure  1 ) , but far from these prominent deep ocean ventilation sites , ocean observations seem to show the largest deep ocean heat uptake at midlatitudes in the eastern north atlantic ( figure  1 ) [ desbruyres et al . ,
warming pattern in the north atlantic at middepth . temperature differences ( t in color , c ) in the 300700  m layer between the 20072012 and the 19972002 periods ( i.e. , average values are estimated for each time range ) .
figure  1 is the zoom into the north atlantic in figure s1 ocean reanalysis system 4 ( oras4 ) ( d ) ( see figure s1 procedure description for details ) .
arrows represent ( in red ) warm and salty upper water currents , ( in dark blue ) cold and fresh deep water currents , ( in black , dotted ) the subpolar gyre , and ( in black , dashed ) the subtropical gyre .
deep convection sites in the irminger and greenland seas as well as the position of the santander standard section in the bay of biscay are marked by green squares .
midlatitudes in the eastern north atlantic ( mlena ) , considered as the oceanic area including the mode water formation area of the eastern north atlantic central water ( enacw ) and the oceanic waters of the bay of biscay and surrounding the iberian peninsula to north africa as in somavilla et al .  
we reexamine the mechanism responsible for the quick transfer of ocean heat from upper to deeper layers in the north atlantic focusing on the role of ena modal waters .
highquality hydrographic data with global coverage are required for providing observational evidence to the mid2000s shift in the north atlantic .
these data are obtained from regularly observed hydrographic sections , which are affected by gaps in coverage both in time and space . to overcome this issue , we combined standard hydrographic observations ( salinity and temperature profiles ) at oceanographic sections regularly maintained in the bay of biscay and the greenland sea ( see figure  1 ) with other data sets that offer more global coverage such as argo floats , globalgridded ocean reanalysis ( european centre for mediumrange weather forecasts ocean reanalysis system 4 , oras4 ) and objective analysis ( met office hadley centre monthly objective analysis , en4 , and world ocean database 2001 , woa01 ) . along the work , temperature and salinity time series from oceanographic sections are systematically presented together with those constructed from argo floats and reanalysis data . in this manner
, we prove that the changes detected at oceanographic sections are representative of larger oceanic areas ( e.g. , figures  2 and 5 ) , and it serves as an important validation for both argo floats and gridded data sets affected by factors such as changes in observation coverage around 2005 [ roemmich et al . , 2015 ] , contamination from model errors in ocean reanalysis data [ balmaseda et al . ,
2013 ] , or the gapfilling strategy in climatological data [ cheng and zhu , 2015 ] . with a similar purpose of increasing the reliability of individual data sets , warming / cooling and saltening / freshening maps obtained from different gridded data are compared ( figures s1 and s2 )
thus , we increase the confidence of the enhanced warming at midlatitudes in the eastern north atlantic as shown in figure  1 , since similar warming distributions showing enhanced heat uptake at midlatitudes in the ena and surrounding the subpolar north atlantic margins are observed in other gridded products ( figure s1 based on en4 ; figure 3c in chen and tung   based on ishii data ; figure 4 to 6 in hkkinen et al .  
based on noaa / national oceanographic data center database ; figure 3.1 on the intergovernmental panel on climate change report [ alexander et al . ,
2013 ] constructed from an update of the annual analysis of levitus et al . ) . a complete description of the data sets , methods , and further details are included in the supporting information . strong mixinginduced deep ocean heat injection at mlena .
( a ) mean temperature ( white line ) between the shallower winter mixed layer depth achieved in warm winters and the enacw upper bound ( 100300  m ) and ( black line ) in the depth range of the enacw core ( 300550  m ) in the enacw formation area ( 4348n , 1020w ) obtained from oras4 reanalysis data set .
( b and c ) same as figure 2a but for temperature and salinity at  
 =27.227.3  kg  m ( corresponding to the + 0.1  kg  m 2005 denser enacw core ) .
the blue circles indicate temperature and salinity at  
 =27.227.3  kg  m as recorded at station 7 of the santander standard section ( figure  1 ) .
the magenta line in figure 2b shows the background average warming trend ( + 0.1c  decade ) at the intermediate depth ( 300700  m ) for the global ocean ( table s1 ) .
the winter of 2005 was exceptionally cold and dry in southwestern europe [ shein , 2006 ; somavilla et al . ,
, the ocean surface at midlatitudes in the eastern north atlantic ( mlena , see figure  1 for geographical reference ) experienced the highest densification ( buoyancy loss ) since the 1960s , causing the winter mixed layer to reach very deep levels likely unprecedented for decades [ somavilla et al .
the induced deep convective mixing resulted in homogeneization of properties of the water column from the surface to density levels of  
 =27.227.3  kg  m , at depth ranges of 350500  m , well beyond the normal airsea interaction extent at the core of enacw (  
 =27.127.2  kg  m ) .
this mixing of ( warmer and saltier ) upper water and ( colder and fresher ) modal waters made the upper waters colder ( 0.4c ) , creating a denser ( + 0.1  kg  m ) , warmer ( + 0.23c ) , and more saline ( + 0.015 ) type of enacw ( figures  2 and s1 ) [ somavilla et al . , 2009 ] . such temperature increase at the enacw core ( + 0.23c ) resulted in a sudden heat content increase of 2.410  j  m for the 350500  m layer ( 3.9  0.610  j  m for the 300700  m layer , table s1 ) .
this suggests that 70% of the estimated 5.810  j  m of heat gained slowly but persistently between 1994 and 2005 in the shallower levels ( 0.7c in the 100300  m layer , table s1 and figure  2a ) was transferred into a denser variety of mode waters through the extreme winter mixing .
in addition , considering that such heat injection took place during winter 2005 , the heat uptake ( q ) at intermediate depths would be 30  w  m ( 5  w  m for the 300700  m layer ) , between 65 and 100 times faster than the average ocean heat uptake at such intermediate depths ( 300700  m ) between 2000 and 2010 ( see table s1 ) .
note that the heat uptake at intermediate depths originates from the heat transfer from the upper layers and not from the atmosphere .
forced by atmospheric cooling , the convective mixing would have released heat to the atmosphere ( see comments in table s1 for additional information ) .
most of the enhanced warming observed at intermediate depths at mlena shown in figure  1 ( + 0.33c gained between 1997 and 2012 ) is explained by the increase in the modal waters during the winter of 2005 ( + 0.23c ) .
the estimated overall heat content increase in the region is 1.2    10  j and accounts for 44% of the total heat content increase at intermediate depths in the north atlantic from 36 to 65n from 1997 to 2012 ( table s2 ) .
evidence suggests that this extreme winter mixing episode may have contributed to the deep ocean heat injection at higher latitudes by inducing ocean circulation changes that enabled northward propagation of the relatively salty enacw .
mlena climatological circulation charts and specific studies indicate slow southward flow west of the bay of biscay and iberia [ pollard et al . , 1996 ; pingree , 1993 ; reid , 1994 ] ( figure  3a ) .
the 27.15 isopycnal became shallower from the midatlantic eastward and peaked at about 18w , beyond which it deepened slowly .
after winter 2005 , the isopycnals became shallower all the way against the european margin ( figures  3b3d ) .
the associated zonal dynamic height sections were thus altered , and the zonal dynamic height slope eastward of 22w changed from almost flat or negative to positive ( figures  3e3 g ) .
hence , the associated meridional flow shifts from southward to northward . estimated as the geostrophic meridional current between 22 and 12w ( v 
 g=/f 

l , where  is the difference between the dynamic height at 12w and 22w , f the coriolis parameter , and l the distance between the two longitudes ) , the meridional current after 2005 estimated from the mean filtered dynamic height from argo floats is + 0.7  cm  s ( northward ) with respect to 0.4  cm  s ( southward ) estimated from woa01 climatological dynamic height ( with reference level at 1500  m both for argo floats and woa01 data ) .
oras4 meridional current data show a comparable northward reversal of the flow centered at 18w from 2005 to 2010 ( mean value =  + 1  cm  s , maximum value  =  + 3  cm  s ) with respect to its southward direction ( mean value  =  0.5  cm  s , maximum value  =  3  cm  s ) from 1990 to 2005 ( see figure s3 ) .
the inversion of the zonal dynamic height slope is a robust feature ( e.g. , with reference level at 1000  m , the meridional current after 2005 estimated from argo floats is + 0.6  cm  s ( northward ) with respect to 0.35  cm  s ( southward ) estimated from woa01 climatological dynamic height ) as also observed in oras4 meridional current ( figure s3 ) .
a similar flow reversal was predicted by a sensitivity experiment in a global ocean circulation model , and anomalous northward advection of salt was obtained as a result of surface cooling at 45n in the ena [ svellec et al . ,
( a ) map of the northeast atlantic showing the location of the nwse , we , and swne sections across the enacw formation area for which the data are shown in figures 3b3 g .
( b  d ) location of the  
  = 27.15 at the nwse , ew , and swne sections marked in figure 3a based on woa01 climatological data ( grey ) and available argo floats in the corresponding section between 2005 and 2010 ( cyan ) .
the blue line indicates the mean filtered value of all the argo floats each 2 of longitude .
( e  g ) idem for the dynamic height between 300 and 1500  m ( see section 2.3 in the supporting information for details ) .
black arrows in figure 3a indicate the southward flow west of the bay of biscay based on pollard et al .
, inferred in figures 3e3 g from the climatological data ( grey dots ) ; and the red arrow indicates the reversed flow after 2005 , shown in figures 3e3 g by argo floats ( blue line ) .
it is worth noting that the maximum water exchange between the subtropical and the subpolar gyres occurs at a density level of  
 =27.227.4  kg  m [ brambilla and talley , 2006 ] .
thus , the + 0.1  kg  m denser core of enacw (  
 =27.227.3  kg  m after 2005 ) might easily get incorporated to the major branch of the north atlantic current flowing northward .
the anomalous northward flow of salty enacw is also inferred directly from largescale hydrography ( see figure  4b ) as will be described in detail in the next section . increasing salinities at middepths after 2005 .
salinity difference ( s in color ) in the 300700  m layer between ( a ) the 20022005 and the 19972002 periods and ( b ) the 20072012 and the 19972002 periods based on en4 data set ( i.e. , average values are estimated for each time range ) .
figure  4 is the zoom into the north atlantic in figure s2 en4 ( c and d ) ( see figure s2 procedure description for details ) .
due to the weak circulation and dominance of mesoscale activity in the mlena region , finding direct evidence of a flow reversal as the one inferred from dynamic height and reanalysis current data is difficult .
argobased deep displacements [ ollitrault and rannou , 2013 ] are an available source , and luckily , a number of floats drifting at enacw levels were active in the region ( see figure s4 ) . the outcome from the analysis of their lagrangian displacements
upper water temperature and salinity in the subpolar gyre and the nordic seas increased since early the 2000s ; this tendency is attributed to the arrival of warmer and saltier western subtropical waters caused by the contraction of the subpolar gyre [ htn et al . , 2005 ; woollings , 2011 ] .
though the higher salinities peaking between 2007 and 2010 in these areas ( see figure  4b ) [ holliday et al . , 2008 ; beszczynskamller and dye , 2013 ; http://ocean.ices.dk/iroc/ ] would be consistent with a higher transport of western subtropical waters , that higher transport does not account for the drop in nutrient concentration in the eastern subpolar north atlantic ( rockall trough , figure  1 ) observed from 2007 onward [ johnson et al . , 2013 ] .
this drop in nutrient concentration is explained by a higher contribution of enacw in the area [ johnson et al . ,
2013 ] , thus adding further evidence in support of the circulation reversal previously described .
the enhanced northward flow in the ena would have reinforced the northwestward retreat of the subpolar front [ hkkinen et al .
, 2013 ] , broadening the passage for both enacw and western subtropical waters toward the north ( see figure  4b ) .
both the eastern and western north atlantic central water types have the potential to increase temperature and salinity at higher latitudes .
however , enacw are saltier than western waters arriving at the eastern subpolar north atlantic [ johnson et al . , 2013 ; hkkinen et al . ,
2013 ] , providing an alternative explanation for the highestsalinity anomalies observed between 2007 and 2010 in the subpolar north atlantic and nordic seas .
assuming mean anomaly propagation speeds of 2 to 3  cm  s [ larsen et al .
2008 ] , a lag of 2 to 5  years is expected for the salinity anomaly ( from the mlena ) to reach first the eastern subpolar north atlantic and then the western subpolar gyre and northern nordic seas . in these areas ,
the arrival of salty southern waters due to the circulation shift in the mlena would have contributed to the highsalinity conditions required for deep ocean heat uptake through density compensated anomalies [ mauritzen et al . , 2012 ; chen and tung , 2014 ] without requiring an increase of the amoc as suggested by chen and tung . such circulation shift and enhanced northward flow in the mlena do not seem to contradict the observed weakening of amoc [ srokosz and bryden , 2015 ] . indeed , a strong northward flow east of 22w and south of 48n in the a priori recirculation area at the ' ovide ' section crossing the north atlantic from portugal to greenland was observed in 2006 , coincident with the lowest amoc measured through the section between 1997 and 2010 [ mercier et al . ,
consistent with the aforementioned spreading timescale and presence of the saltiest upper waters after 2007 in these areas , the heat and salt uptake below 700  m in the subpolar gyre and the nordic seas increased after 2007 ( figures s1 and s2 ) , coinciding with more intensive deep convection activity in both the subpolar gyre ( irminger and labrador seas , [ yashayaev and loder , 2009 ; kieke and yashayaev , 2015 ; de jong et al . , 2012 ; piron et al . ,
2015 ] ) and nordic seas ( greenland sea , [ beszczynskamller and dye , 2013 ] ) .
after years of sustained warming and increasing salinity in the upper layers , the first deep convection events at each of these sites resemble the enacw transformation in 2005 . in the deep water formation region of the greenland sea ,
a strong mixinginduced deep ocean heat injection event occurred through density compensated anomalies in 2008 ( figure  5 ) .
upper ocean salinity in this region had been increasing since the mid1990s , and the 2008 winter cooling made the surface waters dense enough to trigger deep convection with temperature and salinity higher than typical values ( figure s5 ) . as a consequence of the overturning ,
the heat and salt gained slowly over time by the upper layers were directly injected to great depths , with a resulting sudden temperature and salinity increase of 0.16c and 0.01 , respectively ( figures  5c and 5d ) . with a less stratified water column ,
convection occurred at greater depths in the subsequent years as compared to previous decades in the greenland sea [ beszczynskamller and dye , 2013 ] .
coincident with the arrival of highly saline waters in the upper layers of the subpolar gyre , a strong winter mixing occurred in the irminger sea in 2007 [ piron et al . ,
figures  5a and 5b show sudden increases in temperature ( 0.15c ) and salinity ( 0.01 ) at depth , related to this event . as in the greenland sea , subsequent deep convection events followed in later years ( 2008 , 2009 , and 2012 ) [ de jong et al . , 2012 ; piron et al . ,
2015 ] . finally , coinciding with the arrival of highly saline waters in the upper levels , deep convection has reached levels deeper than 1700  m in the labrador sea in 2008 after more than a decade of weak convection , being repeated in 2012 and 2014 [ kieke and yashayaev , 2015 ] . as in the greenland and irminger seas , the first deep convection events appear to be associated with downward transfer of heat and salt [ see yashayaev and loder , 2009 , figure 2a ; kieke and yashayaev , 2015 , figure 4b ] . in general , the steplike temperature increases at depth at deep convection sites caused by the strong mixinginduced heat injection events are equivalent to heat uptake pulses between 16 and 25  w  m , which is between 50 and 100 times faster than the average heat uptake below 700  m in the global ocean ( see table s1 ) . as in the mlena ,
concurrent with the deep ocean heat injection , the deep convective mixing events would have released heat to the atmosphere , favoring deceleration of the north atlantic heat uptake from the atmosphere , as apparently observed since the mid2000s [ chen and tung , 2014 ; hkkinen et al . ,
deep ( a ) temperature and ( b ) salinity ( 14001600  dbar ,  
 =27.7827.8  kg  m ) in the irminger sea ( 6062n ; 3540w , green square in figure 
grey dots indicate vertically averaged argo temperature and salinity values in the pressure range of interest from individual profiles .
deep ( c ) temperature and ( d ) salinity ( 15001800  dbar ,  
 =28.05828.062  kg  m ) in the greenland sea ( 7476n ; 06w , see figure  1 ) .
blue circles indicate vertically and spatially averaged conductivitytemperaturedepth ( ctd ) casts in the area and pressure range of interest and grey dots vertically averaged temperature and salinity of individual argo floats .
the magenta line in figures 5a and 5c shows the background global ocean average warming trend below 700  m ( + 0.015c decade ) ( table s1 ) .
blue lines highlight the years 2007 ( in figures 5a and 5b ) and 2008 ( in figures 5c and 5d ) .
the primary cause of the extremely cold winter in 2005 in the mlena was the occurrence of atmospheric blocking anomalies in the ena .
in addition to its effects on ocean circulation described in section  2 ( i.e. , regional enhanced mixing altering the largescale density structure ) , blocking anomalies in the ena generate a dipole of temperature anomalies between europe ( cold ) and greenland ( warm ) associated with a strong northward meridional atmospheric flow [ barriopedro et al . , 2008 ] with a potential to reinforce largescale ocean circulation changes . blocking anomalies over the ena occurred previously in the late 1950s and 1960s and extreme winter mixing events ( similar to that in 2005 )
enhanced northward flow of enacw may have occurred and contributed ( along with western subtropical waters ) to the intrusion of warmer and saltier waters into the eastern subpolar north atlantic and nordic seas during the 1960s [ holliday et al . ,
2005 ] , as well as enhanced heat uptake by the subpolar gyre through densitycompensated temperature anomalies [ mauritzen et al . , 2012 ; chen and tung , 2014 ; hkkinen et al . , 2013 ; williams et al . ,
the heat uptake in the 1960s had smaller amplitude compared to the 2005 event , and this may be related to the sustained warming trend in the global upper ocean in the few decades prior to 2005 .
the ultimate driver for the occurrence of blocking anomalies at midlatitudes is under debate [ kintisch , 2014 ] .
both warmer arctic and tropical pacific have been associated with a wavier jet stream responsible for an increase in blocking events in the north atlantic and higher frequency of extreme weather events in general [ francis and vavrus , 2012 ; trenberth et al . ,
the blocking pattern developed in winter 2005 in the ena region is not inconsistent with natural climate variability [ barnes et al .
however , overall global warming is expected to affect the development of such atmospheric features [ alexander et al . ,
2013 ] as well as the subsequent ocean response . estimating the global significance of ocean heat transfer from upper to deeper levels by the mechanism described here
requires a better model representation of the frequency and duration of wintertime blocking events , which have been underestimated in most global climate models thus far [ hov et al . ,
since the mid2000s , in contrast to the enhanced ocean heat uptake observed in other ocean basins , the north atlantic deep ocean heat content has significantly increased without apparently intervening ocean heat uptake from the atmosphere . about 40% of the increase in deep ocean heat content increase below 300  m
is centered at midlatitudes in the eastern north atlantic and occurred relatively quickly around the mid2000s far from prominent deep ocean ventilation sites .
we suggest that a sequence of connected strong mixinginduced heat injection events originating from this region triggered large transfer of heat from upper to deeper ocean layers , implying a deceleration or halt of the north atlantic heat uptake from the atmosphere .
a large transformation of modal water masses in the eastern north atlantic in the winter of 2005 played a major role in the observed quick transfer of heat from upper to deep ocean layers but also aided a reverse of the southward circulation of saltier waters from the eastern side of the temperate atlantic with far reaching consequences ; the northward flow of saltier waters into the subpolar north atlantic and nordic seas favored additional heat injection through deep convection events in later years .
thus , we agree that existing highsalinity conditions in the late 2000s in areas of deep convection have aided heat transfer to the deep ocean , but the suggestion that an increase of the amoc is behind such salinity increase is not supported by observations
. a greater influence of the saltier enacw and broadened passage for both enacw and western subtropical waters toward the north appears to be a more plausible explanation .
anomalous atmospheric patterns such as the one behind this shift are not unique to the last decade , although they may have been exacerbated under global warming .
 abstractprior to the 2000s , the north atlantic was the basin showing the greatest warming . however , since the mid2000s during the socalled global warming hiatus , large amounts of heat were transferred in this basin from upper to deeper levels while the dominance in terms of atmospheric heat capture moved into the indopacific . 
 here we show that a large transformation of modal waters in the eastern north atlantic ( ena ) played a crucial role in such contrasting behavior . 
 first , strong winter mixing in 2005 transformed ena modal waters into a much saltier , warmer , and denser variety , transferring upper ocean heat and salt gained slowly over time to deeper layers . 
 the new denser waters also altered the zonal dynamic height gradient reversing the southward regional flow and enhancing the access of saltier southern waters to higher latitudes . 
 then , the excess salinity in northern regions favored additional heat injection through deep convection events in later years . 
Introduction
Data Sets
Transformation of NorthEast Atlantic Modal Waters in 2005 and Implications for BasinScale Circulation
Increasing Salinities and Deep Ocean Heat Injection at Higher Latitudes
Further Remarks
Conclusions
Supporting information