Comparative Analysis of Summer Upwelling and Downwelling Events in NW Spain: A Model-Observations Approach

: Upwelling and downwelling processes play a critical role in the connectivity between o ﬀ shore waters and coastal ecosystems, having relevant implications in terms of intense biogeochemical activity and global ﬁsheries production. A variety of in situ and remote-sensing networks were used in concert with the Iberia–Biscay–Ireland (IBI) circulation forecast system, in order to investigate two persistent upwelling and downwelling events that occurred in the Northwestern (NW) Iberian coastal system during summer 2014. Special emphasis was placed on quality-controlled surface currents provided by a high-frequency radar (HFR), since this land-based technology can e ﬀ ectively monitor the upper layer ﬂow over broad coastal areas in near-real time. The low-frequency spatiotemporal response of the ocean was explored in terms of wind-induced currents’ structures and immediacy of reaction. Mean kinetic energy, divergence and vorticity maps were also calculated for upwelling and downwelling favorable events, in order to verify HFR and IBI capabilities, to accurately resolve the prevailing surface circulation features, such as the locus of a persistent upwelling maximum in the vicinity of Cape Finisterre. This integrated approach proved to be well-founded to e ﬃ ciently portray the three-dimensional characteristics of the NW Iberian coastal upwelling system regardless of few shortcomings detected in IBI performance, such as the misrepresentation of the most energetic surface dynamics or the overestimation of the cooling and warming associated with upwelling and downwelling conditions, respectively. Finally, the variability of the NW Iberian upwelling system was characterized by means of the development of a novel ocean-based coastal upwelling index (UI), constructed from HFR-derived hourly surface current observations (UI HFR ). The proposed UI HFR was validated against two traditional UIs for 2014, to assess its credibility. Results suggest that UI HFR was able to adequately categorize and characterize a wealth of summer upwelling and downwelling events of diverse length and strength, paving the way for future investigations of the subsequent biophysical implications.


Introduction
Upwelling (UPW) and downwelling (DOW) phenomena represent a key role in the strong physical connectivity between offshore waters and coastal ecosystems. Wind-driven coastal UPW has been extensively studied along the eastern edges of the world's major ocean basins, as it has relevant (IBI) ocean forecast system [34] (Figure 1a) to adequately reproduce the three-dimensional features of two persistent UPW and DOW events during summer 2014 in the NW Iberian Peninsula (Figure 1b). To the extent that UPW and DOW circulation patterns have been previously reported to be recurring [7], their presence also becomes predictable, and thus anticipatory strategies can be efficiently prompted to face environmental affairs such as larval transport pathways or the destination of accidental oil spills [29]. To this end, quality-controlled HFR hourly current estimations, satellite products and in situ observations from buoys and a Conductivity-Temperature-Depth (CTD) device were used in concert with IBI hydrodynamic model. A thorough characterization of the ocean signatures associated with UPW and DOW activities can be achieved thanks to their complementary and interdependent nature. While traditional instrumented platforms provide a close approximation of "ground truth" for remotely sensed estimations, HFR observations enhance numerical simulations by resolving fine-scale processes in intricate regions with complex-geometry configurations. In turn, hydrodynamic models can reciprocally serve as integrative connectors of sparse in situ observations and gappy HFR surface current maps by offering a seamless predictive picture of the three-dimensional ocean state. Secondly and most importantly, to characterize the variability of the NW Iberian UPW system through the development of a novel basic ocean-based coastal UPW index, generated from HFR-derived hourly surface current estimations (UI HFR ). Assuming the immediate oceanic response to wind forcing, the proposed high-frequency index can act as a proxy of UPW and DOW conditions. In this context, UI HFR was compared with two preexisting six-hourly UPW indexes for the entire 2014, with the purpose of evaluating its consistency and variability. Likewise, the ability of UI HFR to categorize UPW and DOW events during summer 2014 was also qualitatively assessed.
This work is structured as follows: Section 2 outlines the observational and modelled data sources, along with the methodology adopted. Results are presented and interpreted in Section 3. A comprehensive discussion and future work are addressed in Section 4. Finally, principal conclusions are drawn in Section 5.

The Galician HFR System
A four-site long-range CODAR SeaSonde network, deployed along the Galician Coast (Figure 1b), was used in this work. The first two sites (owned by Puertos del Estado) were deployed in 2004: Silleiro (SILL) and Finisterre (FINI). Afterward, the network was extended northward, in December 2011, by installing two additional sites on phased approach: Vilán (VILA) and Prior (PRIO), owned by INTECMAR-Xunta de Galicia. Each single radar site operates at a central frequency of 4.86 MHz, with a 29.41 KHz bandwidth, measuring the following: (i) Hourly radial currents, moving toward or away from the site, that are representative of the upper 2 m of the water column. The maximum current speed, horizontal range and angular resolution are 100 cm·s −1 , 200 km and 5 • , respectively. All of those radial current vectors (from two or several sites) within a predefined search radius of 25 km are geometrically combined by applying an unweighted least squares fitting algorithm [35] to estimate hourly total current vectors on a Cartesian regular mesh of 6 × 6 km horizontal resolution. (ii) Thirty-minute wave estimations for five range cells, regularly spaced every 5.1 km, which extend radially from the site. For further details about this dataset, the reader is referred to Reference [22], as the present work is mainly focused on surface circulation.
The hourly data availability during 2014 was significantly high in the center of the domain (above 90%), decreasing to 60-70% in the westernmost borders ( Figure 1b). The spatially averaged data availability was almost 80% for the entire year, revealing that the HFR system performed within acceptable ranges despite the severe FINI breakdown since November 2014.
The specific geometry of the HFR domain and, hence, the intersection angles of radial vectors handicap the accuracy of the total current vectors resolved at each grid point. Such a source of Remote Sens. 2020, 12, 2762 5 of 32 uncertainty is quantified by a dimensionless parameter denominated Geometrical Dilution of Precision (GDOP) [36], which typically increases with the distance from the HFR sites. In this work, a cutoff filter of 2.3 was imposed for the GDOP, to get rid of those estimations affected by higher uncertainties. Consequently, GDOP values remain below 2 in the core of the spatial domain, reaching maximum values of 2.29 at the boundaries of the HFR areal coverage (Figure 1c).
Although this four-site HFR system has been working operationally since 2012, the continued data provision was dramatically impacted by a number of serious breakdowns, mainly due to severe weather episodes that typically battered the Galician Coast during wintertime. Additionally, the global economic recession also challenged the financial support to preserve the infrastructure core already implemented. Since the dataset corresponding to 2014 presented both the highest temporal availability and maximum areal coverage (Figure 1d), this work was focused on UPW and DOW events occurred during this specific year.
Finally, note that the present system is part of the cross-border HFR network in the NW of the Iberian Peninsula, reinforced in the framework of the RADAR ON RAIA project (Interreg V, a Spain-Portugal cooperative program). This initiative aims at capitalizing oceanographic knowledge through the extension and consolidation of this HFR network and the development of user-oriented products.

In Situ Buoys
Basic features of two deep ocean buoys, deployed within the HFR footprint (Figure 1b), are gathered in Table 1. Both in situ devices collect quality-controlled estimations of sea surface temperature (SST), among other physical parameters. Furthermore, since B1 buoy is equipped with an acoustic current meter and a wind sensor, it also provides hourly averaged current and wind vectors at 3.5 m depth and 3 m height, respectively.

Upwelling Indexes
Annual (2014) time series of two coastal upwelling indexes (UIs), based on 6-hourly data of sea level pressure (UI BAIXAS ) and wind (UI B1 ), were used as benchmark in this work to assess the credibility of the proposed ocean-based index. Both are freely available at the IEO (Instituto Español de Oceanografía) website: http://www.indicedeafloramiento.ieo.es/index_en.html. UI BAIXAS , which was calculated at 42 • N, 10 • W, was derived from the FNMOC (Fleet Numerical Meteorology and Oceanography Center, US Navy) database. UI B1 was computed by using wind data obtained from B1 buoy. For further information about these two UIs, the reader is referred to Reference [15].

CTD Device
Water temperature and salinity profiles were measured on a weekly basis, at the inner part of the Rías Baixas, using a CTD device (denoted by a purple box and a blue dot, respectively, in Figure 1b), which belongs to a network of 43 oceanographic stations operated by INTECMAR. Since the IBI ocean forecast system is not able to adequately resolve the jagged intricate coastline deep inside the Rías Baixas, only one of the outermost CTD stations was selected to conduct the comparisons: V5 (located at 42.18 • N-8.87 • W).

Satellite-Derived Products
Two daily satellite-derived products, gap-filled and interpolated on a regular grid, were used in this study, to characterize the UPW and DOW episodes: the first one for the SST (OSTIA: Operational Remote Sens. 2020, 12, 2762 6 of 32 Sea Surface Temperature and Ice Analysis) and the second one for the chlorophyll (CHL) concentration. Essential features of both products are summarized in Table 2. Additional information can be obtained in their respective Product User Manuals (PUMs), freely accessible through the Copernicus Marine Environment Monitoring Service (CMEMS) catalogue (https://marine.copernicus.eu/)

IBI Ocean Forecast System
The CMEMS IBI operational suite consists of a NEMO model v3.6 [37] application, run daily over a regional grid ( Figure 1a) with a horizontal resolution of ∼3 km and 50 unevenly distributed vertical levels. A best estimate of the ocean state (also called "hindcast") and a five-day-ahead forecast are routinely produced for a number of hydrodynamic variables: temperature, salinity, mixed layer depth, zonal and meridional velocity currents and sea surface height, among others. While hourly averaged estimations are provided at the sea surface, daily averaged values are computed for the rest of the three-dimensional water column. The system is driven every 3 h by high-resolution (1/8 • ) meteorological forcings provided by the European Center for Medium-Range Weather Forecast (ECMWF). Initial and lateral open-boundary conditions, imposed from daily 3D outputs from the parent system (CMEMS GLOBAL), are complemented with 11 tidal harmonics. The freshwater discharges are prescribed through 33 points corresponding to the main rivers present in the IBI area. This version of IBI ocean forecast system does not include an assimilation scheme. For further technical specifications, the reader is referred to Reference [38].

Methods
A battery of basic processing steps was applied to the sea surface temperature, salinity, current velocity and wind observations. Given the focus on the low-frequency characteristics of the flow field, any short gap of 6 h or less was filled by applying linear interpolation.
Remote Sens. 2020, 12, 2762 7 of 32 HFR estimations were bilinearly interpolated on the finer resolution IBI hindcast mesh, in order to adequately compare both datasets in a common regular domain. The prevailing wind-driven surface circulation during UPW and DOW episodes was examined from a Eulerian perspective, by means of 10-day averaged patterns of surface currents. Complementarily, spatial maps of the CC between observed and modelled current velocity vectors were also calculated, to study their concordance and infer those areas where the model performance might be more consistent. With the purpose of better understanding the flow dynamics, we estimated a wealth of ancillary diagnostics, such as the horizontal divergence (DIV), the vorticity (VOR) and the mean kinetic energy (MKE).
where N, u and v represent the total number of data, the zonal and meridional velocity fields, respectively. Another method to investigate the flow dynamics and seawater dispersion is the Instantaneous Rate of Separation (IROS), which is a Eulerian metric that determines how an infinitesimally small particle will be moved by an instantaneous velocity field, and is equal to the finite-time Lyapunov exponent (FTLE) at time t = 0 [40]. This diagnostic, successfully used in previous studies with HFR data [41,42], can be calculated from the sum of divergence and total strain: While the FTLE unveils a variety of features that dominate over longer time periods, IROS indicates how the particles react in the selected moment, estimating particle separation from a snapshot without integrating the flow over time. Since IROS does not require time integration, it is a rather simple and useful diagnostic to examine mean characteristics of the flow in the NW upwelling system. High values of IROS indicate potential regions of elevated particle dispersion, which has relevant implications for cross-shelf exchange of passive tracers between offshore and coastal waters.
Finally, a novel ocean-based coastal upwelling index (UI HFR ) was developed for the NW Iberian upwelling system. Assuming a direct relationship between atmospheric forcing and the prompt response of the upper ocean layer, a rather simple index was constructed from HFR-derived hourly estimations of sea surface currents, spatially averaged over the Galician continental shelf. Similar to the classical Ekman upwelling index [43], the UI HFR is defined as follows: where ρ w is the seawater density (1025 kg·m −3 ), C d is a dimensionless empirical drag coefficient (1.4 × 10 −3 ), ρ A is the air density in normal conditions (1.22 kg·m −3 ) and f is the Coriolis parameter. In this case, u and v denote the detided hourly time series of HFR zonal and meridional current velocities (m·s −1 ), respectively. The sign is changed to define positive (negative) magnitudes of UI HFR as response of the predominant equatorward (poleward) surface flow over the Galician shelf.

Skill Assessment of HFR Estimations
The credibility of HFR-derived current data has been extensively assessed worldwide [44][45][46] by conducting quantitative comparisons with independent observations provided by in situ platforms such as drifters, acoustic Doppler current profilers (ADCPs) or current meters (CMs). Routine validation exercises are pertinent since the accuracy of remotely sensed estimations might be negatively affected by inherent problems of radar technology such as radio frequency interferences, reflections from moving ships, antenna pattern distortion, adverse environmental conditions or hardware failures, among others [47].
In this context, a preliminary skill assessment of the Galician HFR system was performed for 2014. Hourly HFR current estimations at the grid point closest to B1 location were compared with concurrent observations from a CM installed in B1 (Figure 1b). Comparisons were attempted for both the zonal and meridional sub-inertial currents, obtained after applying a 10th-order digital low-pass Butterworth filter with a cutoff period of 30 h [48]. This scheme is adequate, as the study is mainly concerned with the low-frequency characteristics of the surface flow.
The visual resemblance between the HFR and CM time series is significantly high for both current components (Figure 2a (Table 3). According to the best linear fit of scatterplots, remote-sensed estimations seem to underestimate, to a small extent, the speed of the surface currents observed at B1 (Figure 2c,d).
Remote Sens. 2020, 12, x FOR PEER REVIEW 8 of 32 where ρw is the seawater density (1025 kg·m −3 ), Cd is a dimensionless empirical drag coefficient (1.4 × 10 −3 ), ρA is the air density in normal conditions (1.22 kg·m −3 ) and f is the Coriolis parameter. In this case, u and v denote the detided hourly time series of HFR zonal and meridional current velocities (m·s −1 ), respectively. The sign is changed to define positive (negative) magnitudes of UIHFR as response of the predominant equatorward (poleward) surface flow over the Galician shelf.

Skill Assessment of HFR Estimations
The credibility of HFR-derived current data has been extensively assessed worldwide [44][45][46] by conducting quantitative comparisons with independent observations provided by in situ platforms such as drifters, acoustic Doppler current profilers (ADCPs) or current meters (CMs). Routine validation exercises are pertinent since the accuracy of remotely sensed estimations might be negatively affected by inherent problems of radar technology such as radio frequency interferences, reflections from moving ships, antenna pattern distortion, adverse environmental conditions or hardware failures, among others [47].
In this context, a preliminary skill assessment of the Galician HFR system was performed for 2014. Hourly HFR current estimations at the grid point closest to B1 location were compared with concurrent observations from a CM installed in B1 (Figure 1b). Comparisons were attempted for both the zonal and meridional sub-inertial currents, obtained after applying a 10th-order digital low-pass Butterworth filter with a cutoff period of 30 hours [48]. This scheme is adequate, as the study is mainly concerned with the low-frequency characteristics of the surface flow.
The visual resemblance between the HFR and CM time series is significantly high for both current components (Figure 2a (Table 3). According to the best linear fit of scatterplots, remote-sensed estimations seem to underestimate, to a small extent, the speed of the surface currents observed at B1 (Figure 2c,d).   Although observations at a single point are unlikely to be fully representative of the oceanographic conditions found over the entire study area, hourly wind measured at B1 was low-pass filtered, depicted every 3 h and later analyzed as a proxy for the local open sea wind regime. An abrupt shift in local wind directions, from July to September 2014, is evidenced in monthly wind roses ( Figure 3). Northerly and northeasterly winds clearly predominated in July ( Figure 3a). By contrast, persistent southerly and southeasterly winds prevailed during September 2014, with strong gusts up to 12.7 m·s −1 (Figure 3b). Following References [54,55], we identified two anomalously long-lasting UPW and DOW events during summer 2014, by looking at persistently dominant wind directions and sustained high wind speeds above a predefined threshold of 3 m·s −1 . The meridional wind regime was permanently directed southward and northward, respectively, during 10 consecutive days when no relevant relaxation (below the imposed cutoff) in between was exhibited (Figure 3c,d). Although Reference [56] already documented the UPW conditions and the resulting coastal dynamics as response to two unusually persistent 10-day wind episodes, transient UPW conditions tend to last shorter time [3,4], with an average length of~3 days on the Galician shelf [19].
Remote Sens. 2020, 12, x FOR PEER REVIEW 9 of 32 Although observations at a single point are unlikely to be fully representative of the oceanographic conditions found over the entire study area, hourly wind measured at B1 was lowpass filtered, depicted every 3 h and later analyzed as a proxy for the local open sea wind regime. An abrupt shift in local wind directions, from July to September 2014, is evidenced in monthly wind roses ( Figure 3). Northerly and northeasterly winds clearly predominated in July ( Figure 3a). By contrast, persistent southerly and southeasterly winds prevailed during September 2014, with strong gusts up to 12.7 m·s −1 (Figure 3b). Following References [54,55], we identified two anomalously longlasting UPW and DOW events during summer 2014, by looking at persistently dominant wind directions and sustained high wind speeds above a predefined threshold of 3 m·s −1 . The meridional wind regime was permanently directed southward and northward, respectively, during 10 consecutive days when no relevant relaxation (below the imposed cutoff) in between was exhibited (Figure 3c,d). Although Reference [56] already documented the UPW conditions and the resulting coastal dynamics as response to two unusually persistent 10-day wind episodes, transient UPW conditions tend to last shorter time [3,4], with an average length of ~3 days on the Galician shelf [19].

Analysis of the Surface Circulation
According to the significant resemblance between observed and modelled mean circulation patterns, it can be stated that IBI seems to properly resolve the large-scale surface dynamics in the study area for both UPW and DOW events ( Figure 4). The UPW episode is characterized by a wind-induced southwestward flow, with an offshore advection of coastal waters to the open ocean, in accordance with Ekman´s theory that postulates a net movement to the right of wind direction in the northern hemisphere (Figure 4a,b). An advective acceleration of surface currents speed over the continental shelf is revealed in both datasets, especially in the periphery of CF.

Analysis of the Surface Circulation
According to the significant resemblance between observed and modelled mean circulation patterns, it can be stated that IBI seems to properly resolve the large-scale surface dynamics in the study area for both UPW and DOW events ( Figure 4). The UPW episode is characterized by a windinduced southwestward flow, with an offshore advection of coastal waters to the open ocean, in accordance with Ekman´s theory that postulates a net movement to the right of wind direction in the northern hemisphere (Figure 4a,b). An advective acceleration of surface currents speed over the continental shelf is revealed in both datasets, especially in the periphery of CF.
By contrast, a rather uniform poleward surface circulation is evidenced in response to prevailing southerly winds during the DOW episode (Figure 4d,e). On the southern inner shelf, the interplay between the existing topographic barriers and the cross-shore transport led to a marked directional change of the coastal flow to the left, intensified in form of a narrow jet, the well-documented IPC, which circulated along the NW Iberian shelf edge [8,57]. Since river plumes are mainly confined landward during DOW conditions [58], this northward coastal jet apparently was not significantly perturbated by impulsive-type freshwater outflows at this stage of the year.  By contrast, a rather uniform poleward surface circulation is evidenced in response to prevailing southerly winds during the DOW episode (Figure 4d,e). On the southern inner shelf, the interplay between the existing topographic barriers and the cross-shore transport led to a marked directional change of the coastal flow to the left, intensified in form of a narrow jet, the well-documented IPC, which circulated along the NW Iberian shelf edge [8,57]. Since river plumes are mainly confined landward during DOW conditions [58], this northward coastal jet apparently was not significantly perturbated by impulsive-type freshwater outflows at this stage of the year.
The primary IBI-HFR disagreement arises from the presence, detached from the coast, of a southern counter-clockwise eddy-like circulation structure ( Figure 4b) and an elongated meander ( Figure 4e) in the modelled currents under UPW and DOW conditions, respectively. Conversely, the HFR-derived circulation patterns seem to be more homogeneous, with no proof of the existence of any coherent vortex (Figure 4a,d).
Maps of CC reflect an overall higher IBI-HFR correspondence in regions close to the shoreline, with the CC index and the related phase lying in the ranges (0.6-0.9) and (0-10 • ), respectively, for both the UPW ( Figure 4c) and DOW (Figure 4f) events. These outcomes are in line with earlier works dealing with the same topic [59,60], where the CC coefficients and the associated (absolute) phases were reported to lie between 0.2 and 0.8 and between 0 • and 20 • , respectively. The degree of accordance tends to decline in areas near the limits of the HFR footprint, where higher GDOP values are indeed encountered, especially within the northernmost sectors ( Figure 1c): the CC coefficient drops to (0.2-0.4), and the associated phase increases up to 30 • (in absolute value) in nearby edges of the HFR domain ( Figure 4c,f). Accordingly, Reference [59] also obtained lower CC coefficients (∼0.4) in remote areas, distant from the coastline. Spatially averaged skill metrics that were derived from the IBI-HFR comparison of surface currents for the two selected events are gathered in Table 4. IBI performance was similar for both episodes, with comparable levels of accuracy for the zonal and meridional current components. Better statistics were obtained over the continental shelf than in open waters. Horizontal divergence at the sea surface was calculated to discriminate zones of contraction and expansion of the flow where vertical flux might be significant in terms of proximal origins and destinations of water particles [27]. Positive surface divergence is generally associated with a convergence at deeper levels and the subsequent vertical uplift of waters (i.e., UPW). By contrast, regions of negative surface divergence (i.e., convergence) are related to downward motions to lower levels (i.e., DOW). Therefore, in order to unveil localized areas of UPW and DOW conditions, maps of horizontal divergence were computed by using both the HFR and IBI datasets ( Figure 5). Under UPW-favorable winds, a common peak of positive divergence is exposed in the central portion of the domain (43 • N of latitude) and also in the periphery of CF, indicating accumulated upward vertical motions and strong UPW. A large belt of positive coastal divergence, extended from FINI to PRIO radar sites, is exhibited in the HFR dataset ( Figure 5a). On the contrary, in the case of IBI, the divergence over the continental shelf seems to be weaker and alternated with areas of horizontal convergence ( Figure 5b). Under dominant southerly winds, strong coastal convergence is evidenced north of CF, indicating that the surface flow is presumably sinking downward (Figure 5c,d). Both the HFR and IBI maps agree to determine an UPW core around CF (larger and stronger for the latter), confined between two bands of horizontal convergence. This is in agreement with previous historical works where CF was identified as a locus of permanent offshore advection [10,12,56]. vorticity maps derived from HFR and IBI data looked rather alike, with similar recirculation structures (Figure 5g,h). A predominantly CCW flow moved northeastward, until reaching the periphery of CF, where the surface currents tended to circulate in an anticyclonic way. As the poleward flow separated from the headland, the vorticity increased offshore of CF, inducing the positive divergence already observed in Figure 5c,d and, hence, the recurrent UPW locus extensively described in this manuscript. While the 10-day averaged surface circulation patterns represented in Figure 4 allow the characterization of the permanent structures and the mean flows, the MKE provides a complementary diagnostic accounting for the temporal variability along the analyzed UPW and DOW periods [61]. The comparison of daily MKE values for both episodes, spatially averaged on and off shelf, shows a good IBI-HFR correspondence, despite the fact that HFR represents a more energetic surface dynamics than IBI model ( Figure 6). Furthermore, the MKE is generally higher in coastal shallower areas (solid lines in Figure 6), while lower MKE values are observed at deep waters (dashed lines).
During the UPW event, the upper ocean over the shelf takes only three days for its MKE to reach a prominent peak around 400 (150) cm 2 ·s −2 in the case of HFR (IBI), followed by a steady weakening of the current field (Figure 6a). The evolution of daily MKE in open waters is similar, albeit rather attenuated. This is especially true in the case of IBI, with MKE (red dashed line) smoothly fluctuating in the range (80-100) cm 2 ·s −2 . Under DOW conditions, there is a significant correspondence over the shelf between HFR observations and IBI model, with the latter resolving properly both the abrupt rise of MKE during the second half of the event and the peak reached by the seventh day (Figure 6b). The horizontal vorticity was also examined as descriptor of the rotational flow (Figure 5e-h). Ten-day averaged maps of vorticity were computed from HFR data and IBI model outcomes to track recirculation patterns in time: Positive and negative values represent counterclockwise (CCW) and clockwise (CW) rotation, respectively. Under UPW conditions, the qualitative resemblance between observed and modelled maps was noticeable, with a strong cyclonic circulation (positive vorticity) over the northern shelf (Figure 5e,f) that transitioned into an offshore-directed anticyclonic flow (negative vorticity) after surrounding CF promontory. During DOW-favorable conditions, the vorticity maps derived from HFR and IBI data looked rather alike, with similar recirculation structures (Figure 5g,h). A predominantly CCW flow moved northeastward, until reaching the periphery of CF, where the surface currents tended to circulate in an anticyclonic way. As the poleward flow separated from the headland, the vorticity increased offshore of CF, inducing the positive divergence already observed in Figure 5c,d and, hence, the recurrent UPW locus extensively described in this manuscript.
While the 10-day averaged surface circulation patterns represented in Figure 4 allow the characterization of the permanent structures and the mean flows, the MKE provides a complementary diagnostic accounting for the temporal variability along the analyzed UPW and DOW periods [61]. The comparison of daily MKE values for both episodes, spatially averaged on and off shelf, shows a good IBI-HFR correspondence, despite the fact that HFR represents a more energetic surface dynamics than IBI model ( Figure 6). Furthermore, the MKE is generally higher in coastal shallower areas (solid lines in Figure 6), while lower MKE values are observed at deep waters (dashed lines).
Changes in daily averaged wind speed (i.e., module), derived from hourly estimations registered at B1 buoy, are in close qualitative agreement with the temporal evolution of MKE values under UPW and DOW states, with no apparent delay in the wind-driven response (Figure 6c,d). Therefore, it seems convenient to further analyze the role of the atmospheric forcing in the modulation of surface circulation.

Wind Forcing
In order to comprehend, to a greater extent, the correspondence between hourly subtidal circulation and prevalent wind conditions, we calculated CC maps for both UPW and DOW events (Figure 7). At first glance, observed and modelled surface currents seem to be readily induced by the forcing of strong winds, with a spatially averaged CC coefficient of 0.64 and 0.74, respectively, under During the UPW event, the upper ocean over the shelf takes only three days for its MKE to reach a prominent peak around 400 (150) cm 2 ·s −2 in the case of HFR (IBI), followed by a steady weakening of the current field ( Figure 6a). The evolution of daily MKE in open waters is similar, albeit rather attenuated. This is especially true in the case of IBI, with MKE (red dashed line) smoothly fluctuating in the range (80-100) cm 2 ·s −2 . Under DOW conditions, there is a significant correspondence over the shelf between HFR observations and IBI model, with the latter resolving properly both the abrupt rise of MKE during the second half of the event and the peak reached by the seventh day (Figure 6b). Such a noticeable increase, higher in the case of HFR, is presumably due to the reported migration of the IPC core from the shelf-break depths to the surface, becoming an intensified northward jet as response to the interplay between stronger DOW-favorable southerly wind regime and the local Galician topography [57]. Again, lower MKE values were observed off shelf during DOW-favorable wind regime (dashed lines).
Changes in daily averaged wind speed (i.e., module), derived from hourly estimations registered at B1 buoy, are in close qualitative agreement with the temporal evolution of MKE values under UPW and DOW states, with no apparent delay in the wind-driven response (Figure 6c,d). Therefore, it seems convenient to further analyze the role of the atmospheric forcing in the modulation of surface circulation.

Wind Forcing
In order to comprehend, to a greater extent, the correspondence between hourly subtidal circulation and prevalent wind conditions, we calculated CC maps for both UPW and DOW events (Figure 7). At first glance, observed and modelled surface currents seem to be readily induced by the forcing of strong winds, with a spatially averaged CC coefficient of 0.64 and 0.74, respectively, under UPW conditions (Figure 7a,b). By contrast, lower correlation coefficients were identified in the southernmost sector of the spatial coverage.
Remote Sens. 2020, 12, x FOR PEER REVIEW 14 of 32 UPW conditions (Figure 7a,b). By contrast, lower correlation coefficients were identified in the southernmost sector of the spatial coverage. During the DOW event, higher correlations (above 0.8) were detected in open waters ( Figure  7c,d), with a relevant drop (below 0.4) along the Western Iberian continental shelf (Figure 7d). Albeit not as noticeable, a similar tendency was encountered in HFR observations (Figure 7c). Both the HFR and IBI agree to determine a locus of minimum CC in the periphery of CF, indicating that the surface circulation around this coastal promontory might not be primarily driven by wind forcing but likely by a blend of factors, including topo-bathymetric constraints, local wind and tidal effects. The spatially averaged CC coefficients were 0.57 and 0.61 for the observed and modelled flow, respectively. Complementarily, the time-lagged CC investigation between subtidal currents and low-pass filtered wind time series (not shown) revealed the timely response of the upper layer flow, observed and modelled, to prevalent wind forcing. As expected, a peak of correlation (above 0.8) was generally detected during the first four hours and later accompanied by a gradual decrease in CC coefficient, During the DOW event, higher correlations (above 0.8) were detected in open waters (Figure 7c,d), with a relevant drop (below 0.4) along the Western Iberian continental shelf (Figure 7d). Albeit not as noticeable, a similar tendency was encountered in HFR observations (Figure 7c). Both the HFR and IBI agree to determine a locus of minimum CC in the periphery of CF, indicating that the surface circulation around this coastal promontory might not be primarily driven by wind forcing but likely by a blend of factors, including topo-bathymetric constraints, local wind and tidal effects. The spatially averaged CC coefficients were 0.57 and 0.61 for the observed and modelled flow, respectively.
Complementarily, the time-lagged CC investigation between subtidal currents and low-pass filtered wind time series (not shown) revealed the timely response of the upper layer flow, observed and modelled, to prevalent wind forcing. As expected, a peak of correlation (above 0.8) was generally detected during the first four hours and later accompanied by a gradual decrease in CC coefficient, in accordance with previous works [53,62,63] which reported correlation peaks in the range (0.42-0.64) at zero lag.

Sea Surface Temperature Analysis
The 10-day averaged map of satellite-derived SST and CHL during the UPW episode highlights the importance of the coastline orientation (Figure 8a). A marked cooling (below 17 • C) along the Western Galician Coast is evidenced, with surrounding warmer waters reaching up to 19 • C. Two main cores of lower SST can be observed: one located to the south of the Rias Baixas and the second one in the vicinity of CF. Isolines of constant CHL exhibit a spatial distribution in concordance with thermal fronts location. Peaks of CHL (above 7 mg·m −3 ) are detected close to the western shoreline, thus suggesting accrued injection of subsurface nutrients. By contrast, the sea surface in the Northern Galician shelf is warmer (around 17.5 • C) and not so nutrient-rich, with lower CHL concentrations, lying between 0.2 and 0.5 mg·m −3 . These results underline the importance of the coastal orientation and confirm that stronger UPW conditions are developed in the western coast of the NW Iberian system, in concordance with References [9,10]. Under DOW conditions, the core of lowest SST (below 18 °C) was confined close to the Northern Galician Coast, wrapped in a strip of warmer waters (around 19.5 °C) that expanded southwestward over the shelf, far beyond B1 location (Figure 8f). The concentration of CHL revealed a significantly lower nutrient fertilization of surface waters, with isolines closer to the western shoreline but, by contrast, more distantly spaced from the northern coast, where a peak of chlorophyll (1 mg·m −3 ) was detected. By comparing the 10-day averaged maps obtained for UPW ( Figure 8a) and DOW ( Figure  8f), it can be readily deduced that thermal fronts (i.e., zones with a pronounced horizontal SST gradient) are often associated with higher biological productivity in coastal regions. Such a relationship between front structures and CHL concentration has been previously documented In order to analyze the SST evolution during the 10-day UPW event, the SST difference between the end and the beginning of such episode was explored by means of OSTIA satellite estimations ( Figure 8b) and IBI model outputs (Figure 8c). Both share similarities like the overall summer warming in open-waters or the signal of upwelled waters (cooler than the original surface water), represented by a characteristic band of low SST close to the western coast. However, discrepancies arise in the intensity of the wind-driven cooling, which is broader and more intense in the case of IBI (|∆T| > 2 • C), extended along the northern shelf in the form of a narrow strip of cool water. On the contrary, remote-sensed estimations indicate a general and relevant warming north of VILA radar site (Figure 8b). An additional difference in the modelled SST is the offshore cooling (Figure 8c), probably related to southern cyclonic recirculation cell (Figure 4b), that is absent from OSTIA map ( Figure 8b).
As OSTIA is a satellite-derived product with intrinsic uncertainties, we also used in situ SST observations from B1 and B2 buoys as a reference benchmark to elucidate unequivocally the magnitude of SST variations during the two wind-induced events here analyzed. Since OSTIA is a daily satellite product, hourly estimations from IBI and both buoys were daily averaged, to gain consistency in the data sampling.
As shown in Figure 8d, in situ observations at B1 (represented by a blue line) exhibit a steep temperature decrease of 2 • C, followed by a slight recovery by the end of the UPW event, with a related net cooling of 1.5 • C. The modelled SST at the closer grid point (red line), albeit exposing a positive daily bias, is in clear qualitative accordance with in situ observations in terms of temporal evolution and net decrease of SST under UPW conditions. However, OSTIA estimations at the B1 location (green line) show a moderate drop of 1 • C and a marked recovery until reaching back the pre-UPW conditions (18 • C), in accordance with the whitish color around B1 in Figure 8b. In situ observations provided by B2 reveal a slight cooling of 0.5 • C and a gradual warming during the rest of the episode that led to a positive gradient in temperature (blue line in Figure 8e). Nevertheless, IBI outputs show a sharper decrease of 1.3 • C, followed by a water heating of 1 • C and an eventual small SST drop, yielding a final negative balance of 0.5 • C. This is already shown in Figure 8c, where the high resolution of IBI can effectively resolve sub-mesoscale dynamics like meanders and the intrusion of UPW filaments off shelf, like the one detached from the coast and presumably wrongly extended northward beyond B2 location. OSTIA estimations (green line) barely captured the SST decrease at the first stage of the UPW event, experiencing later a significant warming until 18.5 • C. Despite such a positive bias, satellite-derived estimations could replicate (and slightly overestimate) the overall warming at B2 location.
Under DOW conditions, the core of lowest SST (below 18 • C) was confined close to the Northern Galician Coast, wrapped in a strip of warmer waters (around 19.5 • C) that expanded southwestward over the shelf, far beyond B1 location (Figure 8f). The concentration of CHL revealed a significantly lower nutrient fertilization of surface waters, with isolines closer to the western shoreline but, by contrast, more distantly spaced from the northern coast, where a peak of chlorophyll (1 mg·m −3 ) was detected. By comparing the 10-day averaged maps obtained for UPW ( Figure 8a) and DOW (Figure 8f), it can be readily deduced that thermal fronts (i.e., zones with a pronounced horizontal SST gradient) are often associated with higher biological productivity in coastal regions. Such a relationship between front structures and CHL concentration has been previously documented [64,65].
The SST evolution during the DOW episode was also investigated with OSTIA satellite estimations ( Figure 8g) and IBI model outputs (Figure 8h). Again, the resemblance between both maps is relevant, with an overall SST increment in coastal areas, ranging from 0.5 to 2 • C. In the case of OSTIA, the warming is smoothly extended over the entire continental shelf and expose a vast peak of 2 • C in the northernmost region, surrounding B2 buoy location (Figure 8g). By contrast, IBI exhibits a narrow belt of warmer waters, confined in shallower coastal areas. A significant increase of SST (above 2 • C) is detected over the edge of the continental break, just in the periphery of B2 buoy location and detached from the coast (Figure 8h). Moreover, note that IBI seems to indicate a cooling in the close vicinity of B1, whereas OSTIA estimations reveal no substantial SST trend in that place (Figure 8g). The daily SST evolution registered at B1 buoy exhibits a net decrease of temperature, with an abrupt drop of 2 • C during the first days and a moderate recovery afterward (blue line, Figure 8i). The IBI model (red line), despite the observed bias (1 • C) in the daily outputs at the grid point closer to B1 emplacement, seems to properly resolve the overall cooling of 1 • C during the DOW event. On the contrary, OSTIA estimations scarcely captured the SST drop and the subsequent warming, reflecting a negligible SST variation between the first and the last day of the analyzed DOW episode. The situation at B2 location was, conversely, rather different (Figure 8j). The SST collected at B2 (blue line) experienced an abrupt increase, followed by an equivalent sharp drop and again by an eventual increase of almost 2 • C. The IBI model (red line) failed to capture those oscillations but agreed to reproduce adequately the net warming, with rather similar initial and final SST values. This time, OSTIA remote estimations properly reproduced the ocean water heating, despite the water temperature overestimation by the end of the event.

Vertical Analysis
In this section, the IBI model outputs were compared against independent (i.e., not assimilated) profiles of temperature and salinity, collected during the entire summer 2014, at one CTD station deployed inside of Rias Baixas (V5, denoted in Figure 1b), with the aim of portraying the vertical structure and gaining an improved insight into the three-dimensional flow geometry under UPW and DOW conditions. To this aim, IBI daily estimations at the closest grid point were used.
As shown in Figure 9a,b, the IBI-CTD accordance in the temporal evolution of the temperature field is relevant. The IBI performance appears to be sound in terms of timing, intensity and vertical extension of successive water-cooling and -warming episodes, likely related to the alternation of wind-driven UPW and DOW conditions, respectively. Especially relevant is the increase of temperature by the second half of September, where the subduction of warmer waters (20 • C) reached 30 m in depth. The agreement between CTD observations and IBI estimations, in terms of onset, strength and duration of this DOW event, is significantly high, as reflected by the moderate differences, evidenced in Figure 9c. By contrast, higher discrepancies can be observed in July 2014 when the main UPW event analyzed here took place. Negative differences seem to indicate that IBI overestimated the vertical uplift of colder waters (Figure 9c).
With regards to the salinity field (Figure 10), the model appeared to overrate the noticeable intrusion of fresher (33.5 PSU) waters into deeper levels (up to 25 m depth), related to downward motions under dominant southerly wind conditions (Figure 10b). However, IBI properly captured the upward flux of deep, dense and saltier (35.5 PSU) waters into the sea surface for several dates (14th and 28th of July, and 19th of August) coincident with UPW-favorable events, as indicated by the negligible differences observed in Figure 10c.
In order to get a more realistic diagnosis of IBI three-dimensional performance, instead of computing the bias (where positive and negative differences might partially cancel each other, leading to a more benevolent skill score), the mean absolute difference (MAD), vertically averaged over the entire study period, was used. A MAD of 0.79 • C (0.48 PSU) was derived from the IBI-CTD comparison for the temperature and salinity fields, respectively. Similar qualitative results were recently reported by Reference [66] during a validation experiment of a MOHID modelling suite against V5 CTD station. Remote Sens. 2020, 12, x FOR PEER REVIEW 18 of 32

New HFR-Derived Upwelling Index: UIHFR
In order to assess the consistency of the proposed UIHFR, an annual validation against two preexisting coastal UIs was conducted for the entire 2014. According to the best linear fit of scatterplot between UIHFR y UIBAIXAS, the agreement between this novel ocean-based index and UIBAIXAS seemed to be significantly high (Figure 11a). The slope and the correlation coefficient were close to 1 and 0.8, respectively. The comparison against UIB1 (not shown) revealed a moderate agreement, with a lower correlation coefficient (0.49).  Figure 1b) by a CDT device (a) and IBI model (b) at the closest grid point. (c) Differences of SAL between the modelled and the observed profiles. Note that this figure is not a Hovmöller diagram but a temporal concatenation of specific dates: Only those CTD in situ observations that successfully fulfilled the quality control have been depicted.

New HFR-Derived Upwelling Index: UI HFR
In order to assess the consistency of the proposed UI HFR , an annual validation against two preexisting coastal UIs was conducted for the entire 2014. According to the best linear fit of scatterplot between UI HFR y UI BAIXAS , the agreement between this novel ocean-based index and UI BAIXAS seemed to be significantly high (Figure 11a). The slope and the correlation coefficient were close to 1 and 0.8, respectively. The comparison against UI B1 (not shown) revealed a moderate agreement, with a lower correlation coefficient (0.49). m 3 ·s −1 ·km −1 in the vicinity of SILL site. A secondary minimum, less intense (about −2000 m 3 ·s −1 ·km −1 ), was detected in northern areas, around the VILA site (Figure 11f). By contrast, a prominent peak of UIHFR (~3000 m 3 ·s −1 ·km −1 ) was confined in the southernmost region under UPW-favorable conditions in May (Figure 11d). Later, a dipole-like structure was observed in July (Figure 11e), with two main cores of maximum UIHFR (close to 2000 m 3 ·s −1 ·km −1 ). Since both cores were coincident with each center of cooler and fertilized water previously evidenced in the satellite-derived SST and CHL map ( Figure  8a), these results seem to confirm the viability of the proposed HFR-derived index.  The temporal evolution of monthly averaged UIs showed that UI HFR could appropriately reproduce the annual cycle along 2014, with a predominant UPW phase during the central part of the year (May-July) and prevailing DOW-favorable conditions during the rest of the months (Figure 11b). The visual resemblance with UI BAIXAS is noticeable, although UI HFR appears to generally overestimate the intensity of DOW scenarios. In contrast, the accordance with UI B1 was substantially lower despite the similar evolution.
The spatial distribution of the monthly UI HFR values exhibited recurrent features over the Galician continental shelf (Figure 11c-f). During DOW-favorable months, this coastal index persistently reached two localized minimums in the periphery of SILL and VILA radar sites. In February, UI HFR dropped to −4800 m 3 ·s −1 ·km −1 (Figure 11c). In September, it only decreased until −3900 m 3 ·s −1 ·km −1 in the vicinity of SILL site. A secondary minimum, less intense (about −2000 m 3 ·s −1 ·km −1 ), was detected in northern areas, around the VILA site (Figure 11f). By contrast, a prominent peak of UI HFR (~3000 m 3 ·s −1 ·km −1 ) was confined in the southernmost region under UPW-favorable conditions in May (Figure 11d). Later, a dipole-like structure was observed in July (Figure 11e), with two main cores of maximum UI HFR (close to 2000 m 3 ·s −1 ·km −1 ). Since both cores were coincident with each center of cooler and fertilized water previously evidenced in the satellite-derived SST and CHL map (Figure 8a), these results seem to confirm the viability of the proposed HFR-derived index.
With the purpose of gaining further insight into the capabilities of the new UI HFR , we selected a variety of UPW and DOW events of diverse duration and strength, during summer 2014, and later characterized them by means of IROS metric. During July, the evolution of the three UIs shared some similarities (Figure 12a): three UPW episodes were clearly evidenced (persistently positive UI magnitudes), while a short and more elusive DOW event was categorized according to UI BAIXAS and UI B1 negative values. Both the time-averaged surface circulation and the spatial distribution of IROS under UPW conditions presented common peculiarities (Figure 12b,c,e). A prominent peak of IROS was clearly evidenced along the northern shallower coastal waters, expanding westward offshore in the form of a longitudinal belt (Figure 12b,c). As a result of a temporal outage that affected the northern HFR sites, the spatial coverage was substantially reduced during the last part of the month (Figure 12e). However, the strip of higher IROS values was still partially observable. By contrast, the surface circulation and IROS patterns related to the short DOW episode were kind of noisy (Figure 12d). With the purpose of gaining further insight into the capabilities of the new UIHFR, we selected a variety of UPW and DOW events of diverse duration and strength, during summer 2014, and later characterized them by means of IROS metric. During July, the evolution of the three UIs shared some similarities (Figure 12a): three UPW episodes were clearly evidenced (persistently positive UI magnitudes), while a short and more elusive DOW event was categorized according to UIBAIXAS and UIB1 negative values. Both the time-averaged surface circulation and the spatial distribution of IROS under UPW conditions presented common peculiarities (Figure 12b,c,e). A prominent peak of IROS was clearly evidenced along the northern shallower coastal waters, expanding westward offshore in the form of a longitudinal belt (Figure 12b,c). As a result of a temporal outage that affected the northern HFR sites, the spatial coverage was substantially reduced during the last part of the month (Figure 12e). However, the strip of higher IROS values was still partially observable. By contrast, the surface circulation and IROS patterns related to the short DOW episode were kind of noisy (Figure 12d).  During August 2014, the resemblance between the three different Uis was again noticeable (Figure 13a). A relevant peak of UPW (above 4000 m 3 ·s −1 ·km −1 ) was captured by each UI, although UI HFR indicated not only more intense UPW conditions but also anticipated in the timeline. Four different episodes were eventually categorized. The first one, comprised between the 11th and 14th of August, was chosen according to high UI HFR values, while the other remaining indexes indicated weak UPW or even DOW conditions. However, the time-averaged map of surface currents confirmed the typical offshore deflection of the flow, associated with UPW-favorable conditions (Figure 13b). Equally, we could get a glimpse of the aforementioned belt of high IROS magnitudes between FINI and VILA sites, despite the limited areal coverage. Two additional UPW events were selected, and the resulting time-averaged maps revealed similar characteristics in terms of surface circulation and IROS distribution (Figure 13c,e), in spite of the different time span (nine days versus two days). Again, a rather noisy circulation map was obtained under DOW conditions (Figure 13d), although a longitudinal belt of high IROS values was still present between FINI and VILA sites. These peculiarities were in clear accordance with the previous DOW event analyzed in July (Figure 12d).
Remote Sens. 2020, 12, x FOR PEER REVIEW 22 of 32 During August 2014, the resemblance between the three different Uis was again noticeable (Figure 13a). A relevant peak of UPW (above 4000 m 3 ·s −1 ·km −1 ) was captured by each UI, although UIHFR indicated not only more intense UPW conditions but also anticipated in the timeline. Four different episodes were eventually categorized. The first one, comprised between the 11th and 14th of August, was chosen according to high UIHFR values, while the other remaining indexes indicated weak UPW or even DOW conditions. However, the time-averaged map of surface currents confirmed the typical offshore deflection of the flow, associated with UPW-favorable conditions (Figure 13b). Equally, we could get a glimpse of the aforementioned belt of high IROS magnitudes between FINI and VILA sites, despite the limited areal coverage. Two additional UPW events were selected, and the resulting time-averaged maps revealed similar characteristics in terms of surface circulation and IROS distribution (Figure 13c,e), in spite of the different time span (nine days versus two days). Again, a rather noisy circulation map was obtained under DOW conditions (Figure 13d), although a longitudinal belt of high IROS values was still present between FINI and VILA sites. These peculiarities were in clear accordance with the previous DOW event analyzed in July (Figure 12d).  During September 2014, the three indexes generally looked alike (Figure 14a) and jointly identified a relatively weak two-day DOW episode (days five and six). The associated circulation pattern depicted the well-documented poleward flow (Figure 14b). A delimited peak of IROS was observed, confined in coastal areas between FIST and VILA sites. A strong DOW event was classified according to the extremely persistent negative values of UI HFR and UI BAIXAS (around −6000 m 3 ·s −1 ·km −1 ). Again, the surface waters flowed northward with higher IROS values detected over the shelf break and the NW sector (Figure 14c). Between the 22nd and 25th of September, UI BAIXAS and UI B1 indicated the existence of predominant UPW conditions, in opposition to UI HFR , which simply denoted a steady relaxation of DOW until reaching a calm state (Figure 14a). The related circulation pattern indeed exhibited a prevailing poleward flow over the shelf, whereas a cyclonic recirculation could be observed in the southernmost sector of the HFR domain, likely induced by northerly winds that started to constantly blow (Figure 14d). Finally, a brief DOW event (27-28 of September) was selected according to UI HFR negative values, although the other two indexes reflected calm conditions. The associated map of surface currents confirmed the validity of the proposed ocean-based index: The flow was effectively oriented to the north, and an elongated strip of high IROS values was encountered along the shelf break (Figure 14e).
Remote Sens. 2020, 12, x FOR PEER REVIEW 23 of 32 During September 2014, the three indexes generally looked alike (Figure 14a) and jointly identified a relatively weak two-day DOW episode (days five and six). The associated circulation pattern depicted the well-documented poleward flow (Figure 14b). A delimited peak of IROS was observed, confined in coastal areas between FIST and VILA sites. A strong DOW event was classified according to the extremely persistent negative values of UIHFR and UIBAIXAS (around −6000 m 3 ·s −1 ·km −1 ). Again, the surface waters flowed northward with higher IROS values detected over the shelf break and the NW sector (Figure 14c). Between the 22nd and 25th of September, UIBAIXAS and UIB1 indicated the existence of predominant UPW conditions, in opposition to UIHFR, which simply denoted a steady relaxation of DOW until reaching a calm state (Figure 14a). The related circulation pattern indeed exhibited a prevailing poleward flow over the shelf, whereas a cyclonic recirculation could be observed in the southernmost sector of the HFR domain, likely induced by northerly winds that started to constantly blow (Figure 14d). Finally, a brief DOW event (27-28 of September) was selected according to UIHFR negative values, although the other two indexes reflected calm conditions. The associated map of surface currents confirmed the validity of the proposed ocean-based index: The flow was effectively oriented to the north, and an elongated strip of high IROS values was encountered along the shelf break (Figure 14e).

Discussion
According to the results derived from the annual comparison of filtered HFR hourly time series against independent in situ current observations (Figure 2), it can de concluded that the HFR system performance was reliable for 2014, operating within suitable ranges and in line with the typical values previously reported in the literature (Table 3). RMSE (correlation) values were lower (higher) for the zonal velocity component than for the meridional: 5.72 cm·s −1 vs. 8.59 cm·s −1 (0.74 vs. 0.53). To some extent, such disparity might be due to the geometry of the HFR system where the site closest to B1 buoy, SILL, is installed at almost the same latitude (42 • N in Figure 1b). Thus, SILL's electromagnetic signal was the strongest over the B1 location and provided reliable radial currents, moving toward or away from the site in the latitudinal west-east axis.
Detected differences between HFR and CM current time series could be, in part, attributed a number of relevant elements, as indicated by References [67,68]: (i) sensors' limitations (and the related instrumental noise); (ii) mismatch in the horizontal sampling (HFR provided current estimations spatially averaged over extensive areas, whereas PC gave point measurements); (iii) disparity in the depth sampling (HFR estimations were vertically integrated and representative of the upper 2 m of the water column, while CM collected at a specific 3.5 m deeper level); (iv) variations in the mixed layer and thermohaline stratification (that might impact on the vertical shear of ocean currents, especially in summertime when the two selected UPW and DOW events occurred).
The inspection of the time-averaged surface circulation patterns unveiled that the HFR adequately resolved well-known oceanographic aspects of UPW and DOW processes in the Galician Coast, such as the prevailing southwestward flow or the settlement of the poleward IPC along the NW Iberian shelf edge, respectively [7][8][9][10][11][12][13][14][15][16][17]. Equally, modelled currents appeared to reproduce satisfactorily the large-scale surface dynamics for both episodes, in close qualitative agreement with HFR-derived maps (Figure 4b,e). The correspondence was higher in regions near the shoreline, with CC coefficients generally exceeding 0.7 (Figure 4c,f). In terms of overall skill metrics ( Table 4) that were spatially averaged over the common domain, the CC, zonal and meridional correlation coefficients were 0.50, 0.51 and 0.44 (0.50, 0.51 and 0.42) under UPW (DOW) conditions, respectively. These magnitudes are similar to those documented in earlier works focused on HFR-model comparisons [59,60,69] and provide added confidence in both the approach adopted and the quantitative consistency of the IBI model. Over the Galician continental shelf, model performance was even better for both the zonal and the meridional current velocities: Higher correlations and lower RMSE emerged in the ranges (0.58-0.66) and (8.26-11.87) cm·s −1 , respectively. In summary, the accuracy of IBI to portray the surface dynamics was not only acceptable but also similar for the UPW and DOW episodes.
The principal IBI-HFR differences stemmed from the existence of two southern cyclonic recirculation cells in the modelled currents under UPW and DOW conditions (Figure 4b,e). This is likely attributable to the significantly higher horizontal resolution of IBI, which is based on an eddy-resolving NEMO model application (∼3 km versus 6 km). Moreover, Reference [20] previously opened the debate over the convenience or not of imposing a search radius of 25 km (already mentioned in Section 2.1) to estimate HFR total current vectors that presumably might lead to an oversmoothness of the surface circulation field. This is particularly true in the NW Iberian waters, where coastal meanders and migratory eddies have been documented to exist and have ecological implications as nutrient conveyors [70,71]. In this context, previous studies demonstrated that small-scale coastal processes are poorly reproduced (or even misrepresented) in coarser grid meshes, suggesting the strong necessity of finer horizontal resolution in both hydrodynamic models [72] and HFR systems [46] to properly capture such littoral phenomena. Since this topic remains open and disputable, further investigations should be conducted to shed light on it.
Analysis of surface divergence and vorticity fields from both IBI and HFR datasets corroborated not only the key role of the Galician shoreline orientation in modulating UPW conditions but also the importance of CF promontory and its ambient waters as a locus of permanent positive vorticity and surface divergence, independently of the dominant along-shore wind regime ( Figure 5). IBI was apparently successful at reproducing broad-scale features of the HFR-derived patterns of divergence and vorticity. To the authors' knowledge, this is the first attempt to characterize UPW and DOW conditions in the Galician waters by means of divergence and vorticity fields derived from both hourly modelled and remote-sensed surface currents. Similar initiatives were effectively addressed in the west coast of the USA [26,27,29], also identifying localized areas of positive divergence at the tip of headlands, such as Point Loma [26] or Point Reyes [27] in California, that ultimately induced topographic upwelling. With regards to the main discrepancies between the modelled and observed vorticity, they were identified in the southern sector of the common spatial domain (Figure 5f,h). Two broad cores of positive vorticity were found in IBI outcomes, presumably associated with the coherent and elongated vortexes in the along-shore direction previously described (Figure 4b,e). Over the rest of the domain, the IBI-HFR correspondence was fairly good.
Spatially averaged MKE was computed as an ancillary indicator of the temporal variability of the flow during the selected UPW and DOW periods ( Figure 6). Although IBI statistically reproduced the daily evolution of MKE, the overall kinetic energy of the modelled currents was lower, indicating that a fraction of the total flow variance was not entirely replicated. Similar conclusions in relation to the limited kinematic content of the modelled flow were previously drawn by References [59,60]. The moderate but systematic underestimation of IBI currents energy at the uppermost level might be related to a limited representation of vertical stratification during the warm summer period. Indeed, results exposed in Section 3.2.5 of this paper seem to confirm this conjecture, as IBI tended to overestimate (in intensity and spatial extension) both the cooling and warming during the UPW and DOW episodes, respectively (Figure 9), thus driving an enhanced vertical mixing that might have partially slowed down the wind-induced currents speed.
The relationship of atmospheric forcing and flow variability was analyzed by means of CC maps between the wind collected at B1 buoy and observed and modelled surface currents (Figure 7). Results appeared to indicate a relevant sensitivity of IBI ocean forecast system to wind stresses in the Galician tidally dominated environment. A time-lagged vector correlation analysis was applied to investigate the reaction of the upper-layer circulation for a 24-hour interval (not shown). Maximum correlations (above 0.7) were found during the first four hours, illustrating that only a brief time is required for the subtidal current to establish the equilibrium with local wind, thus staying in conformity with Reference [56,73]. Once again, HFR and IBI exhibited a similar qualitative performance. Notwithstanding, results should be interpreted with caution since the spatial inhomogeneity of the wind field in Galicia has been previously suggested [10,56]. Although Reference [74] acknowledged that B1 buoy was free of possible errors in the correction of the terrain roughness and was, hence, the better station to represent ocean-atmosphere interactions at both Rias Baixas and the adjacent western shelf, local wind in the complex Northern Galician Coast may substantially differ and cause a decrease of maximum correlation between wind and surface currents.
The 10-day averaged map of satellite-derived SST and CHL during the UPW and DOW events confirmed the relevance of the Galician coastline's shape (Figure 8). Abrupt changes in coastal orientation can induce noticeable wind stress fluctuations and, hence, different UPW conditions with subsequent biophysical implications, as here observed and previously documented by References [9,10,12,56]. Stronger upwelling, with cooler and more enriched waters, was evidenced in the Western Galician shelf [9,10].
The comparison of the SST field observed by OSTIA mission and modelled by IBI during the UPW event revealed interesting differences. Despite sharing similar features, the intensity of the wind-driven cooling was broader and more intense in the case of IBI. These results are in line with a recent work published by Reference [75] where a variety of L4 satellite-derived SST products (including OSTIA) were reported to underestimate the thermal fingerprint of strong coastal UPW, not only in the NW Iberian system but also in other regions worldwide. While performing satisfactorily over broad areas for most of the year, such products seemed to misrepresent UPW-induced abrupt temperature transitions in coastal waters. Notwithstanding, it is worth noting that both in situ and remote-sensing devices are affected by technical limitations and intrinsic uncertainties so SST observations from B1 and B2 buoys should only be considered as a close approximation of "ground truth". Additionally, IBI model daily outputs were compared against temperature and salinity profiles collected at one CTD station, with the aim of gaining further insight into the three-dimensional flow geometry under UPW and DOW conditions (Figures 9 and 10). Therefore, the skill metrics here presented appear to indicate that IBI represented appropriately the main signatures of the vertical thermohaline structure along summer 2014 in the intricate Rias Baixas, where the irregular coastal morphology and abrupt gradients in local bathymetry might handicap precise hydrodynamic models forecast. Focusing on the UPW event here studied (7th-16th July), we see that the results seem to suggest that IBI misrepresented, to some extent, the associated cooling, but accurately resolved the associated salinity variability. On the contrary, the model slightly overestimated the temperature of sinking water masses under persistent DOW conditions (12th-21st September), but failed to capture their salinity. Nonetheless, since this is a single isolated example, ancillary validation exercises along the Galician Coast should be performed in the future. The detected discrepancies within the Rias Baixas might be in part attributed to discordant time sampling: While CTD devices capture instantaneous pictures of the water column state, the IBI model provides daily averaged estimations that could hypothetically include high-frequency effects of a variety of coastal processes, including riverine freshwater discharges, tidal mixing or land-sea breezes.
A variety of constraint factors inherently restrict the IBI-HFR comparison and must be taken into account when interpreting the results previously exposed. Such elements encompass remote-sensor limitations and the aforementioned GDOP, discrepancies in both the depth sampling and the horizontal grid mesh, or any potential limitation in the meteorological forcing used to run IBI ocean forecasting system (which is out of the scope of the present paper but probably deserves a detailed examination elsewhere). In the same line, model comparisons against in situ buoys, satellite-derived SST products or CTD profiles are handicapped by similar idiosyncratic aspects. As data assimilation provides the integrative framework for maximizing the joint utility of observations and ocean forecasting systems, a SAM-2 assimilation scheme [76] was later incorporated in the following version of IBI system (released in April 2018). Since then, altimetry data (along-track sea-level anomalies), OSTIA SST product and CORA (Coriolis Ocean database ReAnalysis) datasets (including ARGO, gliders, drifting buoys, etc.) are routinely assimilated into IBI operational chain in order to enhance its predictive skills. Additional efforts should be also devoted to optimizing the observations-model strategy by combining HFR synoptic circulation maps and high-resolution coastal simulations. In this context, classic observing system experiments (OSEs) could be useful to evaluate the impact of assimilating HFR data on model prognostic capabilities.
With the aim of characterizing the variability of the NW Iberian UPW system, a novel and rather simple coastal UPW index was generated from HFR-derived hourly surface current estimations, the so-called UI HFR . Assuming the prompt and direct reaction of the upper ocean layer to intense and prolonged wind forcing, we think it seemed reasonable to develop an ocean-based proxy of UPW and DOW conditions. The proposed high-frequency index was compared with two preexisting six-hourly UPW indexes for the entire 2014. Results confirmed the consistency of UI HFR in terms of mean values and temporal variability (Figure 11a,b). The spatial distribution of UI HFR along the Galician continental shelf was also sound (Figure 11c-f). The map of UI HFR corresponding to July (Figure 11e) matched well with the remotely sensed SST and CHL fields associated with the 10-day UPW event (Figure 8a), confirming the existence of relevant biophysical interactions in the NW Iberian coastal UPW system.
The skill of UI HFR to categorize and describe UPW and DOW events during summer 2014 was also qualitatively evaluated. The most intense episodes were properly captured by UI HFR , in agreement with the other two indexes (Figures 12a, 13a and 14a). Discrepancies arise from relatively short events of moderate intensity: While UI BAIXAS and UI B1 indicated predominant DOW (Figures 12a and 13a) and UPW (Figure 14a) conditions, UI HFR appeared to contradict the formers, revealing by contrast relaxation episodes (i.e., UI HFR values close to zero). The related time-averaged circulation maps exhibited noisy patterns with recirculation structures: While some sectors depicted a northward circulation, other areas presented a southward flow (Figures 12d and 13d). This was probably due to the transient response of the sea surface dynamics to an abrupt change in the dominant wind regime. These results support the potential of the proposed UI HFR , albeit more exhaustive investigations should be conducted to determine not only the required thresholds to accurately discriminate UPW/DOW events but also its applicability to other regions worldwide.
The IROS diagnostic was calculated for hourly HFR surface current maps and later time-averaged to investigate some aspects of the flow dynamics under specific UPW/DOW events. High values of IROS indicate potential areas for strong dispersion of passive tracers. Under UPW conditions, maps of IROS revealed a significant offshore-directed advection in the form of a longitudinal band of high values, spatially comprised between FINI and VILA sites (Figures 12 and 13) where the coastal orientation abruptly changes. Under predominant northerly winds, an elevated westward dispersion of particles in this subregion is suggested, regardless of the length and intensity of the UPW episode. The analyzed events ranged from two consecutive days of moderate UI HFR (below 500 m 3 ·s −1 ·km −1 ) to 10 days of high UI HFR (above 2000 m 3 ·s −1 ·km −1 ), leading to similar IROS magnitudes and spatial distribution, as exposed in Figures 13e and 12c, respectively.
Conversely, DOW events were characterized by an acceleration of the northward jet, depicting a strong poleward advection associated with high IROS values along the shelf-break, likely due to the interplay between southerly winds and local topography (Figure 14b-e). Another peak of IROS was recurrently observed in the NW sector of the HFR domain, Once again, the ocean response was similar independently of the strength and life span of the chosen DOW episodes. The selected events ranged from two consecutive days of moderately negative UI HFR (about −700 m 3 ·s −1 ·km −1 ) to 10 days of extremely negative UI HFR (far below −5000 m 3 ·s −1 ·km −1 ), yielding similar IROS magnitudes, as shown in Figure 14e,c, respectively.
In order to expand and complete the present work, future efforts should be devoted to analyze a broader number of UPW and DOW events from a Lagrangian perspective, with a focus on particle tracking, transport and dispersion processes, residence times and water renewal mechanisms. Equally, the skill assessment of HFR Lagrangian trajectories against drifting buoys should remain as a preliminary priority before addressing this topic [46,67]. An integrated HFR-model strategy could be effective to better understand the connectivity between offshore and littoral areas along with the biophysical interactions at coastal scales, especially in the vulnerable Rias Baixas ecosystem.
Finally, a question that has emerged in recent years is the role of global warming in potential variations of UPW strength and timing that could consequently affect both the delivery of nutrients into the euphotic zone and the marine fish stocks. Former researches have investigated changes on UPW trends worldwide, reporting contradictory results that are highly dependent on the region and season selected and the characteristics of the dataset used [77][78][79]. With a focus on the NW Iberian system, recent studies have hypothesized that coastal UPW will be not as effective in the future, regardless of the tendency for an intensification of both the Azores High and UPW-favorable winds [80]. The UPW imprint is expected to be less intense during the ongoing century due to the sea surface warming and the enhanced thermal stratification of the water column that might inhibit nutrients exchange through vertical mixing.

Conclusions
In the present contribution, the first focus was placed on the synergistic blend of multi-platform observations and the CMEMS IBI ocean forecast system, working in concert to adequately depict the main oceanographic features of two persistent UPW and DOW events in the NW Iberian coastal system during summer 2014.
Results seem to suggest that the HFR performance was sound and credible for 2014, operating within suitable ranges and providing reliable surface current estimations in near-real time that could be effectively used for coastal monitoring and the characterization of recurring UPW and DOW episodes. Furthermore, the examination of HFR-derived surface circulation patterns confirms that this land-based technology properly resolved well-known oceanographic aspects of the Galician UPW and DOW processes, in accordance with the historical literature.
The significant resemblance between observed and modelled wind-driven currents revealed the strong connectivity between offshore and coastal ecosystems and also confirmed the importance of CF promontory as a locus of stationary peak of upwelling, in agreement with previous studies. Skill metrics showed that IBI performance was quantitatively similar for both events, with a higher degree of accuracy over the continental shelf. The primary role of the predominant wind regime in the prompt modulation of surface circulation, nutrient supply and productivity in Galicia was also evidenced. Again, a considerable accordance between IBI and the HFR was revealed in terms of surface patterns of divergence and vorticity. Equally, a good agreement in the MKE evolution was also observed, despite the reduced kinematic energy of the modelled flow.
Complementarily, the skill assessment of IBI against temperature and salinity profiles collected at one CTD station inside of Rias Baixas demonstrated that the principal signatures of the vertical thermohaline structure along summer 2014 were satisfactorily captured by the model. In conclusion, IBI may be deemed as a consistent regional ocean forecast system, able to properly resolve the large-scale surface dynamics in the NW Iberian coastal system, regardless of few shortcomings, such as the overestimation of the cooling (warming) associated with UPW (DOW) conditions. The adopted approach appears to be well-founded to efficiently characterize the hydrodynamic features of the coastal circulation but also to identify ocean modelling aspects where IBI system is still susceptible to improvement in future operational versions.
The second goal of this contribution was to develop a purely ocean-based UPW index for the NW Iberian coastal system, constructed from HFR-derived hourly surface current estimations (UI HFR ). The proposed index was validated for 2014 against traditional indexes, based on six-hourly observations of sea level pressure (UI BAIXAS ) and wind (UI B1 ). Results obtained were promising, revealing the following: (i) a noticeable concordance between UI HFR and UB AIXAS at different timescales; (ii) a proven ability of UI HFR to categorize properly a wealth of relevant UPW and DOW episodes that took place during summer 2014, regardless of their length and intensity; (iii) the existing biophysical interactions in the NW Iberian coastal UPW system, hence, underlying the strong link between HFR-derived wind-induced circulation and the modulation of SST and CHL fields at the uppermost layer.
Additionally, a Eulerian diagnostic denominated IROS was computed to investigate additional features of the surface flow under UPW/DOW favorable conditions. IROS emerged as a simple useful, easily computable metric that is able to determine areas of high dispersive nature. This has direct implications for the connectivity between open ocean and littoral waters though the cross-shelf transport of pollutants and larvae communities. Equally, IROS may deliver essential information for policy-making and the timely management of search-and-rescue (SAR) operations or accidental oil spills.