A first overview of the 53 year past hydrodynamical variability in the Bay of Biscay from a regional simulaon Guillaume Charria1, Frédéric Vandermeirsch1, Sébasen Theeen1, Charefeddine Assassi1, Raphaël Dussin2 1IFREMER; DYNECO/PHYSED, Brest, France; 2IMCS, Rutgers, New Brunswick, USA Objecves Contact: [email protected]

In the frame of the ENIGME project, we aim to describe and to analyze the past Main findings interannual variability over decadal to mul-decadal periods in the Bay of Biscay This first numerical experiment over the last 53 years highlights promising results about the mul-decadal hydrodynamical and the English Channel (41°N – 52.5°N / -15°W – 4°E). These invesgaons will variability in the Bay of Biscay and the English Channel. Indeed, even if the system is driven by atmospheric condions and the allow idenfying the capabilies and the limitaon of our numerical approaches fluxes at the open boundary condions, our results show a signature of a regional dynamics, which can partly explain for an applicaon in future scenarii. the differences in temperature between global simulaon and observaons. These results will lead improved esmaon of Three main groups of processes are considered the mul-decadal trends in the region including budget with new physical parameters (salinity and currents) not described in previous studies (Michel et al., 2009b). • Interannual evoluons (haline and thermal budgets and currents), Figure 1: Map of the North Atlanc with a • Shelf and slope current system and (sub)mesoscale dynamics, zoom on the studied region including the Following this simulaon, further experiments (at higher spaal resoluon - 2.5Km and from other numerical models) will be • in regional models. Bay of Biscay and the English Channel. performed in the frame of the ENIGME project, to deeply explore main physical processes in the Bay of Biscay and the English In the present study, first results from a long-term simulaon are presented. The Channel. hydrological content and the circulaon are described regards with in situ and remotely sensed data. The interannual evoluon of the temperature and the Temperature and Salinity over 53 years C06023 Water masses in the Bay of Biscay FRAILE-NUEZ ET AL.: MASS TRANSPORT IN THE BAY OF BISCAY C06023 Porcupine Bank (latitude 51°N–53°Nandlongitude salinity are also described. a) b) 12°W–14°W) due to the influence of the cyclonic subarctic (a) gyre. These authors suggest that the interannual variability The first overview of the may be due to either atmospheric forcing in winter or (b) advection of Central Water from further west. On the other Figure 5: T/S hand, the ENACW suffers a salinification due to winter interannual variability reproduced cooling of saline surface water near the European ocean diagrams (a) from margin, west and northwest of the Spanish and Portuguese by our regional simulaon gives coasts [Pollard and Pu, 1985]. Subsequently, a cool, saline Mode Water formed by winter convection subducts south- Fraile-Nuez et al. Model experiment ward into the permanent [van Aken,2001].A promising results. Indeed, in cloud of wide-ranging data points from above the thermo- cline reflects the effects of seasonal stratification and the (2008) during interaction with fresher water near the shelf edge [van Aken, temperature, at 5m depth (Figure 2001]. In the coastal domain of the Bay of Biscay, shallower VACLAN cruises in Numerical simulaon are performed using the Primive Equaon Ocean Model, MARS3D waters (<150 m) are affected by the freshwater; consequently, 3a) and 200m depth (Figure 3c), the interannual, seasonal and mesoscale variability in salinity (Duhaut et al., 2008; Lazure and Dumas, 2008) in the BACH4000 configuraon (without de patterns is affected by the river runoff. The Loire and Garrone August 2005, (b) rivers reduce surface salinity locally on the French continen- we observe a good agreement tal shelf [Puillat et al., 2004]. Low-salinity waters from the CORA-IBI from the CORA-IBI dynamics, Theeen et al., 2014). The modelled spaal resoluon is 4Km for 40 sigma vercal Adour and Garonne estuaries may cover a great part of between the climatology and our the eastern Cantabrian shelf in winter and spring time levels. The atmospheric forcings are based on the ECMWF reanalysis (ERA40 for the [Vincent and Kurc,1969;Vincent,1973;Botas et al.,1989; database, (c) from c) d) Figure 2. Potential temperature-salinity diagram of the Koutsikopoulos and Le Cann,1996]. simulaon. Furthermore, the VACLAN0905 stations. Approximate positions of the water [14] The intermediate waters, located in the density range (c) 3 (d) the regional 1958-1978 period; ERAINTERIM for the 1979-2010 period). The open boundary condions are types discussed in this paper are indicated as: STW, Surface sq 27.3–27.9 kg mÀ ,areSub-ArcticIntermediateWater Thermocline Water; ENACW, Eastern North-Atlantic (SAIW), Mediterranean Sea Water (MSW) and Labrador Sea regional simulaon is closer to the Central Water; MSW, Mediterranean Sea Water; SAIW, Water (LSW). The density ranges for these water masses simulaon based on the ORCA024-GRD100 DRAKKAR simulaon (global, ¼° spaal resoluon). Sub-Arctic Intermediate Water; LSW, Labrador Sea Water; show considerable overlap, which allows both diapycnal and observaons than the global NEADW, North–East Atlantic Deep Water; and LDW, isopycnal mixing [van Aken,2000b].Thepresenceof The simulaon in the Bay of Biscay/English Channel (Figure 1) extends from 1rst January 1958 Labrador Deep Water. the edge enhances the diapycnal mixing (BACH4000), and [Garrett, 1991] and forces the circulation to be aligned with results, suggesng a contribuon salinity and neutral densities (gn)illustratethepresenceof the topography as a subsurface eastern at (d) from the global to the 31rst December 2010 with daily outputs (Theeen et al., 2014 – poster n°B975). various water masses and their associated spatial and intermediate levels [Swallow et al.,1977;Ambar and Howe, density/depth distributions (Figures 3–5) for the Bay of 1979; Reid,1994;Paillet and Mercier,1997;Maze´ et al., of the regional dynamics resolved Biscay. In Figure 2 we present q/S diagram, on which we 1997]. The SAIW is a cold and fresh water mass that simulaon have included several selected constant potential density originates from the northwestern Atlantic. One branch of with higher spaal resoluon. At (isopycnal) contours (sq). Note that Figure 2 is the only one the SAIW subducts and recirculates in the northern cyclonic in this paper where we refer to potential density (at surface gyre, while the remainder subducts below the North Atlantic (GRD100) during pressure), while neutral density (gn) is used everywhere else Current and flows toward the eastern Atlantic margin [Arhan, the opposite, at 400m depth in this paper. The reason for this is that potential density is a 1990]. It is characterized by an intermediate salinity mini- August 2005. BACH4000 3 GRD100 CORA-IBI – an in situ database in the Ireland-Biscay-Iberia region direct function of q and S, so it is a straightforward quantity mum between 27.3–27.9 kg mÀ .TheMSWisknownto sq (Figure 3e), temperature are less to plot in q/S diagram, while neutral density is also a originate from the Mediterranean Sea and circulates through function of depth [Jackett and McDougall,1997]. the Strait of Gibraltar as an outflow. A branch of the MSW e) f) [12]Thermoclinewatersarelocatedbetweenthesurface spreads from Cape St. Vincent northward to at least as far as 3 A dataset has been aggregated over the Ireland-Biscay-Iberia (IBI) region (20°N – 60°N / coherent with observaons due and the isopycnal of sq =27.3kgmÀ (approximately Porcupine Bank at 53°Nandischaracterizedbyasalinity The main water masses in this region (Figure 5a3 ) are represented in both global 600 m). At this level, q/S curves have an inverse ‘‘S’’ shape maximum (S 36–36.2 and sq 27.7 kg mÀ ) in the Bay of 35°W – 12°W) for the period 1958-2012 (Szekely, 2014, pers. communicaon). This dataset is with a salinity minimum (Figure 2). This subsurface salinity Biscay [van Aken and Becker,1996; van Aken,2000b].The regional simulaon too diffusive minimum(Figure prohibits direct 5d interaction) and and mixing regional with water formation simulaons of LSW takes place in ( theFigures Labrador Sea and 5c is ) in agreement with vercal masses above [Pollard and Pu,1985;van Aken,2000b; characterized by a relatively deep salinity minimum centered Gonza´lez-Pola et al.,2005].Justabovethisminimumisthe an extracon from CORA4.0 (Cabanes et al., 2013) global dataset combined with the (solving of this issue in progress). at sq = 27.8 in the q/S diagram, overlying the deeper levels Aggregation of CORA and BOBY permanentprofiles, thermocline, from followed by a narrow band of floats, data where North–East in CORA-IBI Atlantic Deep Water dataset (NEADW) and (Figure 5b). At this stage, the points with the q/S characteristics of Eastern North-Atlantic Labrador Deep Water (LDW) are also observed. The abyssal BOBYCLIM dataset (Vandermeirsch et al., 2010). For the salinity (Figures 3b, d, f), Central Water (ENACW), approximately between the iso- water mass, LDW, is characterized by a low and near-bottom Profiles type comparison is qualitave and more quantave scores are under producon but 3 pycnals of 27.1 and 27.3 kg mÀ . salinity minimum whilst NEADW is characterized by a deep [13]Seasonalthermoclinewaters(STW)arelocated salinity maximum overlaying the LDW salinity minimum The Figure 2 displays the number of vercal despite the lack of observaons, 3 betweenwe can already noce that the MSW (Mediterranean Sea Water) is eroded in the isopycnals 26.4 and 27.1 kg mÀ (approximately [van Aken,2000a]. 150 m). Both temperature and salinity increase upward to the profiles during the whole period with the we observe a divergence of the bottom of the seasonal thermocline, where a salinity maxi- mumregional simulaon ( is observed. Huthnance et al. [2002] observedFigure 5c an 4. Inverse) with salinies lower than 36 for waters between Box Model opening of the internaonal databases interannual variation in ENACW salinity near the northwest [15] An inverse box model has been applied to Box A in soluons between the regional Iberian8°C margin and and a relatively 10°C. large-scale This variability nearresult order to obtain confirms the geostrophic flow the that obeys conservation improvements needed (e.g. open (SeaDataNet, since 1990), the advent of the and the global simulaons from boundary condions, posion 3of10 of vercal levels) in the resoluon of deeper ARGO profilers in 2000’s. We can also noce Figure 3: Yearly averaged temperature (a, c, e) and salinity (b, d, f) over a sub-region (11W-1W / 1979. This divergence can be due structures in the regional simulaon. the impact of measurements of opportunity 43N-50N) at 5m depth (a, b), 200m depth (c, d), and 400m depth (e, f) for the regional model to the switch in atmospheric with the RECOPESCA project since 2007. (BACH4000), the global simulaon (GRD100) and the regional interannual climatology forcings (ERA40 to ERAINTERIM) The modelled circulaon (BOBYCLIM) for the temperature only. and is under invesgaon. BOBYCLIM – a regional climatology CTD XBT PROFILERS OCL data RECOPESCA (a) The BOBYCLIM interannual climatology (hp://www.ifremer.fr/climatologie- – Model and Observaons gascogne/) and the associated dataset Figure 2: Number of vercal profiles a) (Michel et al., 2009a; Vandermeirsch et al., by sources in the CORA-IBI database. Main circulaon paerns are reproduced by Aer more than 50 years of integraon, model results show a good the regional simulaon. In