Fixed-point time series, repeat survey and high-resolution modeling reveal event scale responses of the Northwestern Iberian upwelling
Introduction
The Northwestern Iberian Margin (NWIM) forms the northernmost part of the North Atlantic Eastern Boundary Upwelling System (EBUS), which like other EBUS is governed by the seasonally varying location of a mid-latitude atmospheric high-pressure system (Relvas et al., 2007). In the Eastern North Atlantic, the anticyclonic Azores High migrates between offshore the Iberian Peninsula during summer and south of the Azores during winter. The intensification of the Iceland Low in winter permits the eastward passage of low pressure systems that enhances variability of atmospheric circulation and causes significant precipitation over the NWIM (Vitorino et al., 2002). On interannual scales the strength of the winter variability is related to the North Atlantic Oscillation (NAO) (Miranda et al., 2002).
As with other EBUS, the large scale oceanic circulation along the roughly meridional coast is dominated by a slow equatorward flow, in this case the Portugal Current (Saunders, 1982), on the eastern limb of the subtropical gyre. The region of the NWIM concentrated on here stretches from Cape Mondego to Cape Finisterre (40°N to 43°N in Fig. 1). North of Cape Silleiro (42°06′N), the coastline is made up of a series of embayments, the Rías Baixas; to the south, it is relatively smooth and oriented approximately NNW as far as Porto (41°06′N). The shelf deepens gently to the shelf edge at the 200-m isobath, and is bounded by a steep slope that plunges to 2000–4000 m. To the south of Porto Canyon (41°20′N) the shelf is around 50–60 km wide, while to the north it narrows. Inshore of the 100-m isobath, the inner shelf generally widens in the southward direction.
The upwelling season starts in late spring in response to the northerly winds resulting from the build-up of the Azores high (Wooster et al., 1976). The continental shelf circulation responds on the scale of an inertial period to the atmospheric forcing and, under constant winds, reaches equilibrium after 3 or more days (McClain et al., 1986). Offshore Ekman transport in the surface layer produces upwelling of subsurface cold, nutrient-rich central waters at the coast (Torres et al., 2003, Rossi et al., 2013). As a consequence, the sea level falls nearshore and an equatorward geostrophic jet develops along the front between coastal cold upwelled waters and the offshore warmer waters (Ambar and Fiuza, 1994). The front typically moves offshore with prevailing upwelling-favorable winds (Haynes et al., 1993) as far as the shelf edge (Rossi et al., 2010, Rollo et al., 2020). As the upwelling season progresses, the front becomes unstable and gives rise to filaments of cold, offshore flowing water at preferred locations, including one south of Cabo Silleiro at 42°N (Haynes et al., 1993, Cordeiro et al., 2015). This filament may be associated with the region of intensified upwelling or colder upwelled waters between Cape Silleiro and Porto (Fig. 1 Relvas et al., 2007, Cordeiro et al., 2018). However, Álvarez-Salgado et al. (2003) found that >70 percent of the variability of the coastal alongshore wind stress that forces the system occurs at time scales shorter than 10 days. The seasonal cycle of the NWIM, therefore, is partly masked by mesoscale phenomena occurring on time scales of days to weeks.
The mesoscale variability manifests as intermittent occurrence of upwelling, jets and counter-flows, meanders, eddy development, and upwelling filaments, superimposed on the seasonal changes (Peliz et al., 2002, Peliz et al., 2005, Serra et al., 2002). While extended periods of equatorward wind appear to result in exclusively equatorward currents over outer shelf and upper slope (Rollo et al., 2020), even brief periods of calm or northward wind can allow coastal downwelling and counterflows to develop (Torres et al., 2003) and to extend into the Rias Baixas (Barton et al., 2016).
In winter, wind has a persistent northward component and the Iberian Poleward Current (IPC) advects warm waters from the south along the continental slope and shelf edge (Frouin et al., 1990, Haynes and Barton, 1990, Peliz et al., 2005). The relatively warm, and also saline, surface water of the IPC can be tracked into the Bay of Biscay and further north (Pingree and Le Cann, 1990). This flow has been attributed to the interaction of the meridional density gradient and intermittent southerly winds with the continental slope and shelf (Huthnance, 1984, Peliz et al., 2003). The IPC is thought to be present all year, although it weakens and spreads offshore in spring to become confined to the subsurface during summer, below the upwelling jet (Huthnance et al., 2002). At the start of the upwelling season in May and June, the offshore IPC has been observed to coexist with shelf and slope equatorward flows (Torres and Barton, 2007). Following the end of the upwelling season, from September through January, as southerly winds are increasingly prevalent, the IPC becomes more surface intensified and jet-like (Teles-Machado et al., 2016). The poleward flow is frequently associated with eddies and smaller scale instabilities spun up either side of the current core (Peliz et al., 2003). The year-round variability of forcing results in occurrence of sporadic winter upwelling events (Vitorino et al., 2002), which can have significant biological repercussions (Santos et al., 2004, Ribeiro et al., 2005).
The passage of atmospheric low pressure frontal systems is responsible for precipitation events and increased river outflow on the Iberian coast. The sporadic injections of buoyant river outflow modify the typical annual cycle of stratification of the nearshore surface layers associated with local heat and momentum fluxes between atmosphere and ocean. Winter mixing generally homogenizes temperature and salinity down to 100 m, occasionally 200 m (Reboreda et al., 2014) but shallow buoyant plumes form because of continental runoff of fresh waters and local precipitation. Summer stratified conditions are modified over the shelf by upwelling and also by river outflow. In the absence of other forcing, waters discharged from western Iberian rivers form a plume deflected adjacent to the coast in the direction of coastal trapped waves, as expected (Chao and Boicourt, 1986, Horner-Devine et al., 2015). Peliz et al. (2002) reported that the individual river plumes can merge into a single Western Iberian Buoyant Plume (WIBP). Under northerly winds, westward Ekman transport spreads the plume seawards across the shelf to disappear by mixing with oceanic waters, while under southerly winds, the plume is constrained against the coast and the halocline below it deepens nearshore, even reaching the bottom and becoming vertical with strong winds (Otero et al., 2008, Mendes et al., 2016). During periods of high continental runoff and persistent downwelling, the plume waters from the rivers to the south may reach and enter the Rías Baixas (Sordo et al., 2001). There the surface fresh waters accumulate and warm through shortwave radiation (Barton et al., 2015) so that in subsequent upwelling events, these waters are exported to the shelf, with lower salinity and higher temperature than the locally upwelled waters (Álvarez-Salgado et al., 2000, Torres and Barton, 2007).
Basic understanding of the shelf and slope conditions has been obtained from extensive observational efforts over the years, although these have often been sporadic, with different goals and restricted to sub-regions. However, understanding is limited because slow, often low-resolution ship surveys made in different locations during different years under rapidly changing conditions complicate the interpretation because the NWIM is driven by processes on scales from years to days acting together in a complex manner.
The mean circulation and its inter-annual and seasonal changes in the NWIM have been studied through numerical modeling (Nolasco et al., 2013, Teles-Machado et al., 2016). Idealized studies have helped comprehend the development of the coastal upwelling and the formation of filaments (Røed and Shi, 1999, Meunier et al., 2010, Rossi et al., 2010), while the poleward slope currents have been modeled by Peliz et al. (2003). Other model configurations, with enough spatial resolution to resolve processes on scales down to kilometers, have addressed the winter dynamics of plumes and their interaction with the Rías Baixas (Otero et al., 2008, Otero et al., 2013, Sousa et al., 2014, Mendes et al., 2016).
For the first time, we integrate high-resolution modeling with systematic spatial and Eulerian observations on different scales over an annual cycle to explore details of the shelf/slope hydrography and circulation of the NWIM and its relation to the open ocean (Section 2). A year of cross-shelf transects and moored observations, complemented by a model able to resolve the transient processes, has produced new insights into details of shelf conditions (Section 3). The results are placed in the context of experience in this and other upwelling regions (Section 4) before remaining gaps and recommendations for further study are discussed.
Section snippets
Observational data
From November 2008 to December 2009, repeated surveys were made of a zonal transect from nearshore (8.93°W) to the shelf edge (9.44°W) along 42.1°N off Cape Silleiro (Fig. 1). Eleven sections were made at roughly monthly intervals aboard the Research Vessel (RV) Mytilus comprised of 7 conductivity-temperature-depth (CTD) stations 6.9 km apart. Underway current measurements were performed with an RDI 300 kHz vessel mounted Acoustic Doppler Current Profiler (vmADCP).
The transect was sampled three
Time series of ocean and atmospheric data
Time series of hourly alongshore winds and near-surface and near-bottom temperature, salinity and alongshore component of currents from both model and mADCP observations were compared (Fig. 2). CTD casts sampled at the site were used to verify the instantaneous temperature and salinity, while the near bottom temperature variability was examined with the continuous record from the mADCP.
Time series of temperature and currents at 3 m depth from the Silleiro buoy 55 km offshore (green star in Fig.
Modeling considerations
Regional models, such as ROMS-AGRIF, need boundary and initial conditions from observations (typically oceanic climatologies such as World Ocean Atlas) or from other models with larger domains. In the case of the configuration developed here (NWIc), the boundary was taken from the IBI-MFC reanalysis by CMEMS with 1/12° horizontal resolution. The reanalysis is reinitialized every week using assimilation of Argo floats and satellite products (Aznar et al., 2016), hence providing realistic
Conclusions
The purpose of this study was to combine observation and high-resolution numerical modeling of the NW Iberian upwelling system, taking into account atmospheric forcing, large scale circulation, river inflow and topography, to investigate the complexity of the circulation at different scales. Through a validation procedure that compared modeled and observed time series of temperature, salinity and current velocity, and vertical across-shelf transects, it was shown that the model provided
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgments
This work was carried out under the framework of the multidisciplinary project “Canaries-Iberian Marine Ecosystem Exchanges (CAIBEX)” (CTM2007-66408-C02-01/MAR)(Spanish Ministry of Education and Science). NGF Cordeiro was supported by the Portuguese Science and Technology Foundation (FCT) through PhD fellowship SFRH/BD/101070/2014. Thanks are due for the financial support to CESAM (UID/AMB/50017/2013), to FCT/MEC through national funds, and the co-funding by the FEDER, within the PT2020
References (77)
- et al.
New production of the NW Iberian shelf during the upwelling season over the period 1982–1999
Deep-Sea Res. Part I-Oceanogr. Res. Pap.
(2002) - et al.
The Portugal coastal counter current off NW Spain: new insights on its biogeochemical variability
Prog. Oceanogr.
(2003) - et al.
Surface waters of the nw iberian margin: upwelling on the shelf versus outwelling of upwelled waters from the Rías Baixas
Estuar. Coast. Shelf Sci.
(2000) - et al.
Strengths and weaknesses of the CMEMS forecasted and reanalyzed solutions for the Iberia-Biscay-Ireland (IBI) waters
J. Mar. Syst.
(2016) - et al.
Vertical structure, turbulent mixing and fluxes during Lagrangian observations of an upwelling filament system off Northwest Iberia
Prog. Oceanogr.
(2001) - et al.
Coastal upwelling and downwelling forcing of circulation in a semi-enclosed bay: Ria de Vigo
Prog. Oceanogr.
(2015) - et al.
Ocean surface wind simulation forced by different reanalyses: comparison with observed data along the Iberian Peninsula coast
Ocean Model.
(2012) - et al.
Two-way nesting in split-explicit ocean models: algorithms, implementation and validation
Ocean Model.
(2012) - et al.
Variability of the western Galician upwelling system (NW Spain) during an intensively sampled annual cycle. An EOF analysis approach
J. Mar. Syst.
(2008) - et al.
Ocean margin exchange–water flux estimates
J. Mar. Syst.
(2002)