Environmental Factors Affecting Spatial Dinoflagellate Cyst Distribution in Surface Sediments Off Aveiro-Figueira da Foz (Atlantic Iberian Margin)

Resting cysts of planktonic dinoflagellates, once produced, sink to the seabed where they can remain viable for a long time. These cysts have important ecological roles, such as acting as the inoculum for the development of planktonic populations. Moreover, dinoflagellate cyst records from depth sediment cores are broadly used as a proxy to infer past environmental conditions. In this study, the main objective was to obtain information on the relationships between the spatial distribution of modern dinoflagellate cysts and present-day hydrography in the NW Iberian shelf. Cyst assemblages were analyzed in 51 surface sediment samples with varying grain sizes, collected at different water depths, following nine transects perpendicular to the coast, between Aveiro and Figueira da Foz (Atlantic Iberian margin). Multivariate statistical analyses revealed marked land-sea and latitudinal gradients in the distribution of cysts, and helped investigate how environmental factors [water depth, grain size, sea-surface temperature (SST), sea-surface salinity (SSS), bottom temperature (BTT) and surface chlorophyll-a concentration (CHL)] influence modern dinoflagellate cyst composition and abundances. Three main ecological signals were identified in the modern dinoflagellate cyst assemblages: (1) the heterotroph signal as the main upwelling signal; (2) the dominance of P. reticulatum and L. polyedra signal, indicative of warm stratified conditions, possibly reflecting transitional environments between more active inshore upwelling and warmer offshore waters; and (3) the G. catenatum signal for the presence of mid-shelf upwelling fronts. The almost absence of viable cysts of the toxic and potentially toxic species G. catenatum L. polyedra and P. reticulatum suggests that in the study area, for these species, there is no build-up of significant cyst beds and thus planktonic populations must depend on other seeding processes. These results are the first detailed modern distribution of dinoflagellate cysts in the NW Iberian Atlantic margin (off Portugal), and show a good correspondence with hydrographic features of summer upwelling season in the study area, meaning that they are reflecting water column characteristics and therefore may be used as supporting evidence for the interpretation of stratigraphic cyst records and reconstruction of past marine ecosystems in W Iberia.

Resting cysts of planktonic dinoflagellates, once produced, sink to the seabed where they can remain viable for a long time. These cysts have important ecological roles, such as acting as the inoculum for the development of planktonic populations. Moreover, dinoflagellate cyst records from depth sediment cores are broadly used as a proxy to infer past environmental conditions. In this study, the main objective was to obtain information on the relationships between the spatial distribution of modern dinoflagellate cysts and present-day hydrography in the NW Iberian shelf. Cyst assemblages were analyzed in 51 surface sediment samples with varying grain sizes, collected at different water depths, following nine transects perpendicular to the coast, between Aveiro and Figueira da Foz (Atlantic Iberian margin). Multivariate statistical analyses revealed marked land-sea and latitudinal gradients in the distribution of cysts, and helped investigate how environmental factors [water depth, grain size, sea-surface temperature (SST), sea-surface salinity (SSS), bottom temperature (BTT) and surface chlorophyll-a concentration (CHL)] influence modern dinoflagellate cyst composition and abundances. Three main ecological signals were identified in the modern dinoflagellate cyst assemblages: (1) the heterotroph signal as the main upwelling signal; (2) the dominance of P. reticulatum and L. polyedra signal, indicative of warm stratified conditions, possibly reflecting transitional environments between more active inshore upwelling and warmer offshore waters; and (3) the G. catenatum signal for the presence of mid-shelf upwelling fronts. The almost absence of viable cysts of the toxic and potentially toxic species G. catenatum L. polyedra and P. reticulatum suggests that in the study area, for these species, there is no build-up of significant cyst beds and thus planktonic populations must depend on other seeding processes. These results are the INTRODUCTION Dinoflagellates are a highly diverse group of protists that along with diatoms and coccolithophores are one of the most abundant groups of coastal marine phytoplankton, significantly contributing to primary productivity (Falkowski and Knoll, 2007). As a group, dinoflagellates present different trophic strategies: some are autotrophs (photosynthetic), others heterotrophs, but most of them are mixotrophic (able to combine both trophic modes) (Schnepf and Elbrächter, 1992;Hansen, 2011;Stoecker et al., 2017). The trophic strategy and life cycle traits together with different environmental preferences determine the ecology of the different species (Hansen, 2011).
Around 13-16% of living dinoflagellate species are known to form benthic resting stages (cysts) as part of their lifecycle (Head, 1996). Since then, several new species have been described and their cyst-theca relationship established (e.g., Luo et al., 2019). Once formed, cysts are accumulated in the bottom sediments where they can remain viable for a long time, playing important ecological roles. Benthic resting cysts are key for the survival of the species, acting as overwintering stages, inoculum for the formation of planktonic blooms or as a genetic reservoir (Anderson et al., 2005;Lundholm et al., 2011;Bravo and Figueroa, 2014;Ellegaard and Ribeiro, 2018). In the present-day warming scenario, recent work has highlighted the importance of considering cyst physiology when investigating the response of cyst forming dinoflagellates to environmental change (Brosnahan et al., 2020).
Given their capacity to fossilize, dinoflagellate cyst records from sediment cores are broadly used as a proxy to infer past environmental conditions (see Ellegaard et al., 2017 for a review). Stratigraphic cyst records, and more recently ancient DNA (De Schepper et al., 2019), represent a very valuable tool to study past biodiversity and reconstruct past marine ecosystems. However, interpretations of the cyst record can be biased by several factors related to the variable representation of the planktonic populations in the sediments. Importantly, post-depositional processes such as sediment reworking, horizontal transport and selective degradation of cysts (e.g., Zonneveld and Brummer, 2000;Zonneveld et al., 2018) can complicate the interpretation of the environmental signals (Ellegaard et al., 2017). Moreover, when inferring palaeoenvironmental data from cyst records we must consider that (1) not all dinoflagellate species produce resting cysts -i.e., they will not be represented in the cyst records -and (2) cyst productivity rates can vary between species -meaning some species can be over-or underrepresented in the sediment cyst record. Certain studies provide data which may help in overcoming these limitations, such as studies on cyst-theca relationships and cyst production dynamics (e.g., Susek et al., 2005;Bringué et al., 2013;Matsuoka and Head, 2013), and particularly those oriented to explore the relationships between modern cyst distributions in surface sediments and environmental gradients (e.g., Wall et al., 1977;Dale, 1996;Pospelova et al., 2008;Zonneveld et al., 2009;de Vernal et al., 2020;Van Nieuwenhove et al., 2020). Although a huge effort has been made to document in the Northern hemisphere the distribution of dinoflagellate cysts in surface sediments and their relation to major environmental variables [see Van Nieuwenhove et al. (2020) for an historical review] there are not many studies that investigate at a regional scale dinoflagellate cyst assemblages and their relation to present day hydrography, particularly coastal ecosystems of the Atlantic Iberian margin are still largely underrepresented in reference datasets .
Several dinoflagellates are responsible for harmful algal blooms (HAB) events, many of which produce resting cysts. These events can be associated with various threats in coastal systems, such as seafood poisoning or water discoloration (Hallegraeff et al., 2003). The W Iberian coast is frequently affected by these blooms which raise important economic and health safety concerns (Vale et al., 2008). In this region, within cyst producing HAB species, Gymnodinium catenatum is of particular concern given its association with the human syndrome Paralytic Shellfish Poisoning (PSP). Blooms of G. catenatum have been recorded in the West Iberian margin since late 1970s of the last century (Fraga et al., 1988), and since then have been thoroughly studied. These blooms are characterized by a high decadal and annual unpredictability. The conditions leading to bloom seeding and initiation are still not fully understood, although seeding from planktonic populations is considered the most plausible hypothesis (Fraga et al., 1993;Moita et al., 1998;Sordo et al., 2001;Amorim et al., 2004;Bravo et al., 2010;Smayda and Trainer, 2010). Blooms generally occur during the summer upwelling season in relation to upwelling-relaxation-downwelling cycles and the autumn upwelling-downwelling transition (Fraga et al., 1993;Moita et al., 1998;Pitcher et al., 2010). Upwelling plumes are important features on the dynamics of G. catenatum blooms. Moita et al. (2003) related the presence of G. catenatum in the inshore side of the upwelling plume front with the asymmetric hydrodynamic conditions of upwelling plumes rooted at major capes in the W Iberian coast. The presence of G. catenatum in frontal zones and downwelling conditions has been attributed to its high swimming velocities as a result of the long chain life-form strategy that characterizes this species (Fraga et al., 1988;Moita et al., 2003;Pitcher et al., 2010). Evidence from the sub-fossil cyst record suggests this species recently colonized the NE Atlantic possibly by geographical range expansion from NW Africa (Ribeiro et al., 2012).
Lingulodinium polyedra (=L. machaerophorum) and Protoceratium reticulatum [=Operculodinium centrocarpum sensu Wall and Dale (1966)] are two other cyst-forming species of concern because of the production of a group of potentially nuisance toxins, the yessotoxins. Lingulodinium polyedra is also known to cause water discoloration in W and SW Iberia (Amorim et al., , 2004, with high negative impacts on tourism. Blooms have been reported in summer and early autumn associated with warm stratified conditions (∼17 • C) in the inshore side of nutrient rich upwelling filaments (Amorim et al., , 2004. On a global scale, blooms of L. polyedra have been reported from temperate regions since the beginning of last century, occurring in summer-early fall (Lewis and Hallett, 1997). In the NE Atlantic, the distribution of L. polyedra, based on surface sediment cyst distributions, spans from the equatorial zone to the temperate/sub-polar boundary with a distribution center in NW Africa and the Gulf of Cadiz (Dale, 1996;Lewis and Hallett, 1997;Zonneveld et al., 2013;de Vernal et al., 2020).
Regarding Protoceratium reticulatum, very little is known on the ecology of the planktonic stage in this area, with only few records reported in phytoplankton studies off Portugal (Moita et al., 1998;Loureiro et al., 2011;Silva et al., 2015). In contrast, the cyst record in modern and subfossil sediments from the Atlantic Iberian margin, particularly in the NW Iberian coast, is characterized by the common occurrence of cysts of P. reticulatum -and also of G. catenatum and L. polyedra -(e.g., Amorim, 2001;Amorim et al., 2004;Sprangers et al., 2004;Ribeiro and Amorim, 2008;Ribeiro et al., 2016;García-Moreiras et al., 2015. Based on surface sediment cyst distributions, P. reticulatum is one of the most cosmopolitan species (e.g., Wall et al., 1977;Dale, 1996;Zonneveld et al., 2013). The proportional increase of this species (=Operculodinium centrocarpum) in cyst assemblages seems to be associated with unstable transitional conditions such as the neritic-oceanic boundary or the influence of river plumes (Dale, 1996;. Investigating the oceanographic mechanisms that control bloom dynamics of these and other HAB species is essential for the development of HAB forecasting tools and coastal monitoring programs (Wells et al., 2015;Berdalet et al., 2017). In this context, the study of dinoflagellate cyst distribution patterns in recent sediments and of the overlaying abiotic conditions is particularly interesting to understand how cysts can work as environmental proxies, identify the presence of potential cyst beds and detect rare species (Blanco, 1995;Dale, 1996;Orlova et al., 2004;Anderson et al., 2005;Ellegaard et al., 2017).
In this study, modern dinoflagellate cyst assemblages (concentration and relative abundance) and grain-size were analyzed in 51 surface sediment samples off Aveiro-Figueira da Foz (Atlantic Iberian margin). The well documented hydrography of the study area (Relvas et al., 2007;Oliveira et al., 2019) allowed the investigation of how well the cyst assemblages reflected the water column characteristics. The relationships between environmental variables -grain-size, water depth, sea-surface temperatures (SST), sea-surface salinities (SSS), bottom temperatures (BTT) and chlorophyll-a concentration (CHL) -and community composition were investigated by multivariate statistics to study how environmental gradients affect the present-day spatial distribution of dinoflagellate cysts. This work represents the first detailed dataset on modern dinoflagellate cyst distributions on the NW Portuguese shelf strongly influenced by seasonal upwelling, and may provide reference data for the interpretation of environmental signals from stratigraphic dinoflagellate cyst records to reconstruct past marine ecosystems in coastal environments from related or similar areas.

Study Area
The sampling area is located between the latitudes 40 • 7 and 40 • 49 N and the longitudes −8 • 44 and −9 • 18 W, offshore Figueira da Foz-Aveiro (Iberian Atlantic margin) ( Figure 1A). From south to north, the main freshwater sources in the study area are the Mondego river and the Ria de Aveiro ( Figure 1B). The Mondego river has a catchment basin of 6670 km 2 and an annual average runoff of 500 m 3 s −1 (Cunha and Dinis, 2002;Marques et al., 2002). Several rivers drain to the Ria de Aveiro -an estuary with coastal lagoon characteristics -with a total catchment area of ∼3600 m 2 , Vouga being the largest river, with an average runoff of 29 m 3 s −1 (Dias et al., 1999;Da Silva and Oliveira, 2007). Moreover, during periods of upwelling favorable winds, the study area may also be influenced by the Douro river plume, whose outlet is situated ∼35 km north of the study site (Fernández-Nóvoa et al., 2017). The Douro river has a larger catchment basin (98 073 km 2 ) and a mean discharge of 700 m 3 s −1 , with a significant sediment load to the NW Iberian shelf Oliveira et al., 2002;Fernández-Nóvoa et al., 2017).
Seasonal upwelling (spring-early autumn) and the variation of low-salinity buoyant plumes are among the main oceanographic features affecting the phytoplankton dynamics within the study area (Fiúza et al., 1982;Peliz et al., 2002;Moita et al., 2003;Oliveira et al., 2019). Latitudinal and onshore-offshore environmental gradients (e.g., SST, SSS) have been described during upwelling events, and sub-mesoscale oceanographic processes were suggested as major factors affecting the horizontal distributions of chlorophyll-a and biota (Moita, 2001;Oliveira et al., 2019).
In summer, the mean Chl-a and SST distributions have a typical coastal upwelling pattern: a coastal band of high Chl-a coinciding with low SST, resulting from the upwelling of cold, nutrient-rich, sub-surface waters. The resulting gradient is not uniform and usually presents a maximum at some distance from the coast along the so-called upwelling fronts. The distribution of the SST fronts in the study area have been studied by Relvas et al. (2007) who showed that these fronts are most frequently located around the 50 isobath, and identified a low frontal probability zone at mid-shelf ( Figure 1C).  In the same area, Oliveira et al. (2019) described the bifurcation of the equatorward alongshore flow of upwelled waters originating further north (from 41-42 • N) and the presence of a cyclonic area over the mid-shelf (∼100 m isobath) (see their Figure 10) with a location that coincides with the low frontal probability zone of Relvas et al. (2007). SSS distribution is influenced by the low-saline waters of the Western Iberia Buoyant Plume (WIBP), composed mainly by freshwater loads from the Douro River ( Figure 1A). During strong upwelling events, the WIBP is advected equatorward and covers most of the northern portion of the study area (Oliveira et al., 2019), with SSS increasing southwards and offshore.
During winter the predominant SW winds induce poleward currents along the shelf. However, winter upwelling events may be recorded generating equatorward currents (Relvas et al., 2007). Moita (2001) recorded short pulses of upwelling in the study area in winter, but these were not persistent enough to break the halo stratification at the surface induced by the well-developed buoyant low salinity plume. It is within these low salinity plumes that the low winter Chl-a maxima were recorded (Moita, 2001).
In the study area, the spatial distribution of the bottom sediments is influenced by bottom topography, river inputs and coastal hydrodynamics. The topography of the study region is characterized by a wide continental shelf with isobaths almost parallel to the coast (Figures 1A,B). The coast has no major protrusions except for Cape Mondego, where a persistent upwelling filament has been reported to occur during the upwelling season (Sousa and Bricaud, 1992). Fine sedimentary deposits are associated with rocky outcrops, which act as a barrier to the cross-shelf movement of sediment and also to the general circulation (Oliveira et al., 2001). Latitudinally, two main circulation patterns affecting sediment distribution may be identified. In summer, higher intensity and frequency of upwelling promotes sediment transport southwards and offshore. In winter, when precipitation and storms are more intense and frequent, increased riverine supplies promote the transport of resuspended fine materials northwards (Fiúza et al., 1982;Fiúza, 1983;Vitorino, 1989).

Field Sampling
During the Hydrographic Institute of Portugal (IH)/AQUIMAR Cruise (March 2019), 51 surface sediment samples were collected with a Smith-McIntyre grab in coastal environments between the Ria de Aveiro and the Mondego (Figueira da Foz) outlets (Figure 1). Sampling stations followed nine land-sea transects perpendicular to the coast, corresponding to different grainsizes and water depths ( Table 1). For dinoflagellate cyst analyses, plexiglass tubes (3.6 cm internal diameter) were inserted in the sediment and the top 1-cm layer was collected and stored at 4 • C in the dark for further analysis in the laboratory. For sediment characterization, the bulk grab sample, representative of the top 20 cm of surface sediment coverage, was collected and frozen (−18 • C) until further analysis.

Sediment Processing for Dinoflagellate Cyst Analyses
A replicate of each sample (2 to 30 cm 3 , depending on cyst concentration) was taken for dinoflagellate cyst analyses, which included wet-sieving with distilled water through a 150 µm-Nylon mesh and a 20 µm-stainless steel mesh (Endeccots) after sonication (60 s) (Elmasonic S50R). The retained fraction (between 150 and 20 µm) was centrifuged at 3600 rpm (∼2510 × g) (Eppendorf 5804 R), and the supernatant was then removed, and dinoflagellate cysts were further concentrated by centrifuging the homogenized pellet at 1000 rpm (∼190 × g) in a high-density solution of sodium polytungstate (∼2. 016 g cm −3 ) (Bolch, 1997;Amorim et al., 2001). The floating organic fraction was collected, rinsed twice by centrifugation at 3600 rpm (∼2510 × g) (Eppendorf 5804 R) and recovered with filtered seawater in a final volume of 1 to 10 ml, depending on cyst concentration. Replicated sediment samples were used for dry weight determination (drying at 60 • C until constant weight) and % moisture. The sample dry weight was determined using the previously calculated % moisture.
An inverted light microscope (Leica DMi1) was used for cyst counting and identification under 200x and 400x magnifications. One or more Sedgewick-Rafter chambers (1 ml) (Graticules Optics, United Kingdom) were counted to obtain a minimum of 100 cysts when possible, including both empty and cysts with cell contents (referred as full cysts onward), with an average of 264 cysts per sample. This method was used to obtain relative (percentage values of the total cyst assemblage) and absolute abundances. The latter were determined relative to the volume (cysts . ml −1 ) and dry weight (cysts . g −1 ) ( Supplementary  Tables 1, 2). Phase-contrast morphological examination (and photography) of some specimens were performed on a Zeiss Axiovert 200 microscope. All photographs (both bright-field and phase-contrast) were taken with a Zeiss Axiocam HRc camera. Finally, to keep a permanent collection as a backup, aliquots of all samples were mounted on slides using glycerine jelly and sealed with wax.

Notes on the Identification and Nomenclature of Dinoflagellate Cysts
The identification and nomenclature of dinoflagellate cysts follow Zonneveld and Pospelova (2015), Gurdebeke et al. (2019), Mertens et al. (2020), andVan Nieuwenhove et al. (2020). When species-level identification was not possible, the identification was done at the genus or higher level. The biological names were used preferentially, however, in those cases where the cysttheca relationship is not clear, the paleontological name was used instead (e.g., Quinquecuspis concreta; see Table 2). This was also applied in the cases when the name refers to one of various cyst (morpho)types that are currently linked to only one species (e.g., Votadinium calvum, which illustrates one cyst type of Protoperidinium oblongum).
Unidentifiable spiny brown cysts were grouped as "spiny brown cysts" or SBC ( Table 2) for data analysis purposes and inter-sample cyst-record comparison, due to some common difficulties with their identification (see Radi et al., 2013). Unidentifiable Round Brown cysts (RBC) included folded and broken round brown cysts that were not assigned to any specific genera (probably Brigantedinium sp., Dubridinium sp., etc.), and also specimens of Diplopsalis-type ( Table 2). The latter includes light brown to grayish round cysts, with smooth surface and theropylic archeopyle, which correspond to cysts of Diplopsalis spp., and may also include other round brown cysts with problematic morphology such as Diplopelta sp. and Diplopsalopsis sp. and others described in Mertens et al. (2020).
For the same reason, Echinidinium delicatum and E. granulatum, both characterized by having hollow processes (Zonneveld and Pospelova, 2015), were grouped in E. delicatum/granulatum. On the other hand, taxa that could be misidentified during routine counts were included in the same group, namely cyst types for which Because of their presumably similar ecological affinities, for clustering and RDA analyses, the round brown cysts Brigantedinium spp., Diplopsalis-type, Dubridinium spp., and unidentifiable round (smooth) brown cysts (which include broken or folded brownish cysts with round outline that probably correspond to one of the previous genera) -were grouped in RBC. Additionally, unidentified cysts were excluded and those groups for which identification was doubtful in some samples were grouped to the genus level. A total of 23 cyst types and square-root transformed percentage data (%) were used for RDA analyses (taxa that never contributed >1% were excluded). Note that very consistent results were obtained with or without data transformation.
Six variables were investigated in this study, namely: water depth (WD), sediment mean diameter (Gran), bottom temperature (BTT), sea-surface salinity (SSS), satellite-derived sea-surface temperature (SST), and sea-surface chlorophyll-a concentration (CHL). Initially, six parameters for each variable (except for WD and Gran) were explored by RDA to check their influence on the ordination of samples, namely: the global average, the global median value (50th percentile or p50), the value for the month with the lowest mean value, the value for the month with the highest mean value, and the averages of the trimester with the highest and lowest mean values. These parameters were calculated from daily values from 2003 to 2019; except for BTT and SSS, which were calculated from a data series corresponding to the period 2018-2020. These time intervals were selected as they correspond to the highest resolution of environmental data to which access was available. In the final RDA diagrams (Table 1), only the parameters with the highest weight in the ordination of samples were included. Our choices were based on statistical scores R 2 , the variance explained and the p-value from random permutation tests (999 permutations at a significance level of P ≤ 0.05). Additionally, redundant variables that strongly correlated with more significant variables were discarded (see Oksanen et al., 2015).
Oceanographic variables were obtained from satellite remote sensing and hydrodynamic numerical models. The daily numerical model solutions for bottom temperatures (BTT) and sea-surface salinities (SSS) were gathered from the E.U. Copernicus Marine Monitoring Service (CMEMS), namely from the Operational IBI (Iberian Biscay Irish) Ocean Analysis and Forecasting system, which is based on a (eddyresolving) NEMO model application run at 1/36 • horizontal resolution (Sotillo et al., 2015). The two satellite-derived variables, available as daily, 1 km resolution maps, were the sea-surface temperature (SST), from the Multi-scale Ultra-high Resolution (MUR) SST analysis (JPL, 2015;Chin et al., 2017), and chlorophyll-a (CHL) from CMEMS based on the Copernicus-GlobColour processor, namely from the North Atlantic Chlorophyll product (OCEANCOLOUR_ATL_CHL_L4_REP_ OBSERVATIONS_009_098).

Estimation of Sediment Accumulation Rates
No box-corers were collected during the cruise, so establishing accurate sedimentation rates for sediment samples investigated in this study was not possible. However, in nearby areas (north of the Douro mud patch), Jouanneau et al. (2002) performed Pb analyses in 31 box-core samples collected at depths ranging from −39 to −213 m, and found that accumulation rates ranged between 0.05-0.40 cm year −1 . Considering the results by Jouanneau et al. (2002), we estimated that the top 1 cm of surface sediment samples collected off Aveiro-Figueira da Foz recorded between 2.5 and 20 years of cyst deposition. Moreover, the top 20 cm of surface sediment samples used for grain-size analyses would represent <400 years of deposition. We believe that this temporal mismatch between the two type of samples is not a major problem since the obtained grain-size distribution accurately represents the main (modern) distribution patterns observed in the study area, being consistent with available sediment cartography, such as the Seabed Substrate maps from EMODnet Geology (available at 1 ) and the sedimentological charts from the Hydrographic Institute of Portugal 2 ,used as base to compile the EMODnet information. Therefore, we considered that the obtained grain-size distribution was well representative of modern sediment distribution and represents the best available data to be compared with the modern distribution of dinoflagellate cyst assemblages.

Grain Size and Dinoflagellate Cyst Distribution
Mean grain diameter ranged from 0.06 to 2.20 mm. Sands dominated all samples with percentages varying between 51 FIGURE 2 | Spatial distribution of cysts (empty cysts and cysts with cell contents) (cysts . ml −1 ) (A) and sediment grain size (B) along the nine transects studied. Table 3). Coarsest samples usually corresponded with stations at mid-depths. Grain size also increased in northern transects, in front of Ria de Aveiro, with gravel percentage ranging 40-49% in samples B3, B14, B15, B16, B17, and B28. In contrast, proportions of finer sediments (silt and clay) increased southwards and in deeper samples ( Figure 2B).
A total of 45 dinoflagellate cyst types were identified, excluding unidentifiable types round brown cysts (RBC), spiny brown cysts (SBC), and brown peridinoids ( Table 2). Morphotype richness or number of cyst types per sample is shown in Supplementary Table 4 (Unidentifiable types are counted as one type), with a maximum of 33 and a minimum of 10 cyst types. Morphotype richness generally increased northwards and at mid depths.
Full cysts of the HAB species G. catenatum and P. reticulatum occurred in negligible amounts while no full cysts of L. polyedra were recorded. In contrast, many full cysts of calcareous taxa (Scrippsiella trochoidea, spherical-type cf. Scrippsiella sp., cf. Ensiculifera tyrrhenica) and RBC were identified (Supplementary Table 4).
The cysts of cf. E. tyrrhenica have a calcareous wall, yellowish under bright field microscopy. It is characterized by a flat side with pentagonal outline and a triangular side suggesting a "Napoleon hat" (Montresor et al., 1993) (Images 2 and 4, Supplementary Plate I). Cysts are 30-37 µm in diameter and the cyst content has a clear red accumulation body (Images 1-4, Supplementary Plate I). These morphological features coincide with those described for cysts of Pentapharsodinium tyrrhenicum (Balech) Montresor, Zingone & D.Marino ex Head, described in Montresor et al. (1993), which has been recently transferred to the genus Ensiculifera (Li et al., 2020). To date, on the Atlantic Iberian margin, E. tyrrhenica has only been reported in this study and in Amorim (2001). In the latter study, a dinocyst survey was done along the whole Portuguese coast but E. tyrrhenica was only recorded in Ria de Aveiro, located in the vicinity of the present study area.
The highest proportions of autotrophic cysts (38.2-78.4%) were recorded in the B82-B66 transect, located in front of FIGURE 4 | Spatial trends in relative abundances (%) of heterotrophic (orange) and autotrophic (blue) taxa along the nine land-sea transects studied.

Multivariate Analyses and Environmental Influences
According to cluster analysis constrained by latitude, cyst assemblages could be classified in three groups (southern, intermediate and northern samples). Clustering also indicated that the major change in cyst composition was between northern and the remaining (southern + intermediate) samples, cyst assemblages from southern and intermediate samples being more similar (Figures 8A,B, 9A,B).
Redundant discriminant analysis (RDA) helped in investigating the main environmental factors influencing dinoflagellate cyst distribution. The final RDA model (Figures 9A,B) significantly ordinated the samples according to cyst composition and the influence of environmental variables (p-value < 0.05). The variance explained by the two main axes (RDA1 e RDA2) was 30.4%. Mean sea-surface temperature (SST) of the warmest trimester (July to September) and water depth (WD) were the most relevant environmental variables (i.e., those contributing the most to explain the variability of cyst assemblage distribution), followed by grain size (Gran) and chlorophyll-a concentration (CHL). Despite the small amplitude of variation of the SST values used for the RDA (see SD, Table 1), it represents an average of daily values for a period of 16 years (2003 to 2019), and therefore, small changes in SST could actually represent relevant changes in the environment.
The main variability gradient was defined by RDA1 (explaining 25.1% of the total variance), which was positively correlated with SST, SSS and WD and negatively correlated with CHL, BTT and Gran. It also reflected relevant differences between northern samples (generally negative RDA1 values), and southern/intermediate samples (generally positive values on RDA1), in good agreement with the cluster classification results (Figure 8). Northern samples positively correlated with CHL, BTT and Gran, whereas southern samples positively correlated with SST, SSS and WD ( Figure 9A).
Heterotrophic cyst taxa ordinated to the left side of the biplot, with CHL, BTT and Gran. The cyst taxa with the most negative scores on RDA1 were Q. concreta and RBC ( Figure 9B). In contrast, autotrophic taxa generally presented positive values on RDA1, with SST, SSS and WD. G. catenatum and P. reticulatum had the highest values on RDA1. G. catenatum, G. microreticulatum, and Scripsiella sp., plotted on the upper right quadrant with WD. In contrast, P. reticulatum, L. polyedra, S. mirabilis/hyperacanthus, and Spiniferites sp. plotted on the lower right quadrant, representing warmer conditions ( Figure 9B). The distribution of the heterotrophic species P. americanum did not follow the general pattern detected for other heterotrophic species. P. americanum ordinated to the right, with most of the autotrophic species.
An inshore-offshore environmental gradient was also reflected in RDA. Deeper samples were ordinated to the upper right quadrant of the biplot (Figure 9A), with higher abundances of G. catenatum (Figure 9B). On the other hand, shallower samples FIGURE 6 | Dinoflagellate cyst percentages obtained for the studied 51 surface samples. Samples are arranged by nine transects perpendicular to the coast, from north to south (top to bottom). Those cyst types for which identification was doubtful were summed to a higher category, i.e., Brigantedinium spp. includes Brigantedinium sp. and cysts of Protoperidinium avellana and P. cariacoense; Dubridinium spp. includes Dubridinium sp. and D. caperatum; Impagidinium spp. includes Impagidinium sp., I. paradoxum, and I. aculeatum; Spiniferites spp. includes Spiniferites sp., S. bentorii, and S. delicatus; RBC includes unidentifiable round browns and cf. Diplopsalis-type. ordinated to the lower left quadrant of the biplot (Figure 9A), with higher proportions of Protoperidinium species and most of heterotrophs ( Figure 9B).

DISCUSSION
This study investigated the distribution of benthic dinoflagellate resting cysts in recent sediments from the shelf off Aveiro-Figueira da Foz (Atlantic Iberian margin) in relation with present day environmental drivers. These studies are particularly relevant since cysts are the only fossilizable stage of dinoflagellates and understanding how they reflect present day environmental conditions will contribute to a better interpretation of environmental signals from stratigraphic cyst records.

Distribution of Sediments and Cyst Concentrations
Coastal hydrodynamics and topography play an important role in the distribution of sediments and dinoflagellate cysts. The main topographic and morphological features in the shelf area covered by this study that may affect the distribution patterns of sediments and dinoflagellate cysts are the Cape Mondego and sheltered rocky outcrops ( Figure 1B). Coarse sediments are common in the shelf of the Iberian Atlantic coast, while sheltered deposits of mud (silt and clay) are strongly influenced by the vicinity of rocky outcrops that act as traps favoring the accumulation of fine sediments (Figure 1B; Dias et al., 2002;Oliveira et al., 2007). Additionally, in coastal areas, sediments can be resuspended by waves and tidal currents and while sands are deposited rapidly once conditions calm down, finer sediments (cysts included) are more susceptible to be transported offshore . This could partially explain the decrease in sediment grain size with water depth observed in the study area ( Figure 2B). Changes in total cyst concentrations along the 9 transects perpendicular to the coast (Figure 2A) were also highly related to sediment grain size. Both cyst concentrations and fine-grained sediments (silt and clay) increased seawards and southwards (Figures 2A,B, 3A). This relationship, recognized since early works on dinoflagellate cyst assemblages from recent sediments, supports the interpretation of cysts behaving as fine silt particles in the sedimentary regime (Dale, 1976).

Environmental Influences on the Dinoflagellate Cyst Assemblages
The dinoflagellate cyst assemblages (concentrations and relative abundances) revealed marked inshore-offshore and latitudinal gradients that may reflect different environmental factors. It should be noted that shelf surface sediments are the result of several years of deposition, reworking and bioturbation. Likewise, the associated dinoflagellate cyst assemblages represent an integration in time and space of the species encystment dynamics and of sedimentation processes. They reflect the dominant ecological processes over several years rather than a seasonal signature (Dale, 1976).
Along the W Iberian shelf, the oceanographic process regarded as the major source of seasonal and spatial phytoplankton variability is coastal upwelling, which usually lasts from spring to early autumn, with a peak in July (Moita, 2001;Pitcher et al., 2010). Phytoplankton succession is characterized by a strong link to the mixing-stratification patterns both at the seasonal and the upwelling event scale. In general, spring and summer upwelling events are characterized by the dominance of chain-forming diatoms, and the summer-autumn stratified conditions by dinoflagellates. In some years, the autumn upwelling-downwelling transition, as well as summer upwelling relaxation, are dominated by chain forming dinoflagellates such as the HAB species G. catenatum (Moita et al., 1998;Pitcher et al., 2010). In this study, Redundant discriminant analysis (RDA) helped to explore the relationships between environmental parameters -grain-size (Gran), water depth (WD), seasurface temperatures (SST), sea-surface salinities (SSS), bottom temperatures (BTT) and chlorophyll-a concentration (CHL) ( Table 1) -and the dinoflagellate cyst assemblages in surface samples (Figures 9A,B). It should be noted that the spatial resolution of the environmental data is a bit lower than the dinoflagellate cyst records. Despite this, the available SST data seems to well represent the main environmental distribution gradients that predominated in the study area for the last decades (Peliz et al., 2002;Oliveira et al., 2019).
Redundant Discriminant Analysis suggests the existence of two main environmental regimes (Figures 9A,B). One (positive values on RDA1), which included the southern and offshore (deeper) sites, was characterized by higher SST, lower BTT and lower primary productivity (lower CHL), suggesting warm stratified conditions. The other (negative values on RDA1), which included the northern sector and inshore (shallower) sites, was characterized by lower SST, higher BTT and enhanced primary productivity (higher CHL), suggesting upwelling influence.
These two main environmental regimes were reflected on the cyst assemblages by a clear trophic segregation. In the first case, autotrophic taxa, such as P. reticulatum, L. polyedra, S. mirabilis/hyperacanthus, G. catenatum, and G. microreticulatum, dominated the assemblage, while in the second case heterotrophs were the dominant group, particularly RBC, including Brigantedinium spp., and Q. concreta (Figure 9B).
Heterotroph dominance has been previously identified as the primary cyst signal for upwelling (Dale, 1996;Bringué et al., 2014). This interpretation has been supported by a 2-year sediment trap study in Santa Barbara Basin, North-East Pacific, influenced by seasonal upwelling (Bringué et al., 2013). The study covered several cycles of upwelling and relaxation and indicated that active upwelling was characterized by an increase in the fluxes of Brigantedinium spp. (Protoperidinium spp.) and other heterotrophic species such as cysts of Protoperidinium conicum. The authors also reported that biogenic silica, used as a proxy for diatom production, showed a strong positive and significant correlation with heterotrophic cyst types, suggesting a causal predator-prey relationship between the two groups. Very little is known on the effect of turbulence in heterotrophic dinoflagellates. Smayda (2002) reviewed the importance of swimming speed in the capacity of dinoflagellates to overcome vertical and horizontal velocity rates. Heterotrophs, namely protoperidinioids, were among the species that presented swimming speeds compatible with surviving upwelling velocity rates. This would allow them to prey on species characteristic of active upwelling, such as diatoms. Evidence from other studies investigating dinoflagellate cysts in relation to hydrological conditions and marine productivity also related the distribution of heterotrophic dinoflagellate cysts to the distribution of their prey, especially diatoms (e.g., Price and Pospelova, 2011;Elshanawany and Zonneveld, 2016). In the study area, the upwelling pattern is characterized by an onshore narrow band of cold water with the isotherms almost following the bottom contours (Fiúza, 1983;Relvas et al., 2007), with a zone of higher upwelling probability running parallel to the coast ( Figure 1C). It is noteworthy that this pattern can also be seen in the distribution of heterotroph dominance (particularly RBC) in inshore samples and the northern sector ( Figure 5B). In W Iberia, as in other upwelling areas, diatoms are the dominant group during active upwelling when vertical stratification is weakened (Moita, 2001;Moita et al., 2003;Pitcher et al., 2010;Villamaña et al., 2017) suggesting the observed distribution of heterotrophic dinoflagellates is reflecting the abundance of their prey as observed in similar systems. These results are consistent with the use of the "heterotroph signal" as an upwelling productivity signal in W Iberia (Dale, 1996).
The interpretation of the heterotroph dominance as an upwelling signal is preferred to other nutrient enrichment processes, such as nutrient loads from river input, since this would also promote increased stratification which would favor autotrophic dinoflagellates relative to other non-flagellated groups such as diatoms. The latitudinal gradient on heterotrophic dominance (higher abundance in northern sites) can at least partially be explained by a recurrent upwelling front observed on the shelf north of the study site, offshore the Douro estuary (e.g., Haynes et al., 1993;Cordeiro et al., 2015), and the coincidence with a zone of higher upwelling probability ( Figure 1C; Relvas et al., 2007), i.e., the greater abundances of heterotrophic cysts in northern samples responding to increased productivity ( Figure 9B).
The southern and offshore sites, characterized by warm stratified conditions (positive side of RDA1), were dominated by autotrophic cyst taxa (Figures 9A,B). However, the statistical analysis suggested two groups that may represent distinct ecological signals. One group, characterized by P. reticulatum, L. polyedra and S. mirabilis/hyperacanthus, and a second group characterized by G. catenatum, G. microreticulatum and the heterotrophic species P. americanum.
The first group showed a strong (positive) correlation with SST ( Figure 9B). Spiniferites mirabilis has a wide geographical distribution in temperate to equatorial regions and occurs in all major upwelling areas (Zonneveld et al., 2013). In the Northeastern Pacific, influenced by upwelling, Pospelova et al. (2008) described S. mirabilis associated with a warm water assemblage which also included L. polyedra. In the Po river estuary (Eastern Mediterranean), Zonneveld et al. (2009) also found an association between S. mirabilis and warm upper waters. Regarding L. polyedra, globally, modern distribution of cysts has been related to coastal environments of low and middle latitudes, with a distribution center south of the Iberian Peninsula, in the Gulf of Cadiz, and the northwest African coast (Lewis and Hallett, 1997;Marret et al., 2020;de Vernal et al., 2020). At regional scales, many studies related high abundances of this cyst to warm and stratified environments (e.g., Leroy et al., 2013;Ribeiro et al., 2016). The flagellate stage of L. polyedra is known to perform daily vertical migrations (Lewis and Hallett, 1997 for a review) which may represent an important competitive advantage in stratified environments. Along the S and W Iberian shelf, planktonic blooms of L. polyedra have been reported in late summer-early autumn associated with warm stratified waters (SST of ∼17 • C) that were adjoined by upwelling plumes (Amorim et al., , 2004. In agreement with these observations, results from the above-mentioned sediment trap study by Bringué et al. (2013) in the seasonal upwelling system off Southern California, showed that higher fluxes of theca and cysts of L. polyedra coincided with periods of higher SST, stronger stratification and reduced primary productivity. It was even suggested that 15 • C could represent a lower temperature threshold for increased L. polyedra fluxes. These observations support the interpretation that in seasonal upwelling systems, as suggested herein, the increased abundance of L. polyedra is the main signal for an increase in SST and water stratification. This is in good agreement with previous reports of cysts of L. polyedra (=Lingulodinium machaerophorum) from areas adjoining upwelling systems off South California and Northwest Africa (Dale, 1996).
Protoceratium reticulatum (=Operculodinium centrocarpum) is considered one of the most cosmopolitan species (Wall et al., 1977;Dale, 1996;Marret et al., 2020;de Vernal et al., 2020). Recent work has shown the existence of at least three genotypes with different ecological preferences (Wang et al., 2019). Despite this, cysts of P. reticulatum have been successfully used as part of ecological signals related with environmental instability in different marine ecosystems. This is the case for the coastal/oceanic transition, frontal zones or the instability associated with the propagation of river plumes (Dale, 1996;. Studies on the fossil and sub-fossil record of dinoflagellate cysts in the NW Iberian shelf have found increased abundances of P. reticulatum associated with major environmental shifts (Amorim and Dale, 2006;Ribeiro et al., 2016;García-Moreiras et al., 2018. Relvas et al. (2007) investigated the average effect of upwelling structures along the Portuguese coast and found a relatively weak upwelling signal and low frontal probability zone over the mid-shelf (between ∼100 and 200 m depth) of the present study area (Figure 1C), coincident with observations by Oliveira et al. (2019) of a cyclonic area limited inshore and offshore by upwelled waters (their Figure 10). The high relative abundance of P. reticulatum and the presence of L. polyedra in the outer stations of the southern transects is here interpreted as characterizing the location of the transitional zone between more active inshore upwelling and warmer offshore waters (compare Figures 1C, 7). This would be in agreement with the concept of upwelling relaxation taxa proposed by Reynolds (2001, 2003) for L. polyedra and P. reticulatum. The robustness of this signal could be further tested in the future extending the sampling area offshore well into the low frontal probability zone. Ribeiro et al. (2016) investigated the main changes in cyst records in the NW Iberian shelf during the last ∼150 years. They reported a many fold increase in total cyst abundance that was concomitant with the increase in L. polyedra and P. reticulatum. This was interpreted as reflecting the combined effect of warming and increased water stability. In the Ria de Vigo, further north in the west Iberian margin (RdV, Figure 1A), García-Moreiras et al. (2018) also interpreted L. polyedra as an indicator of warm stratified environments, related to river inputs. Their interpretations are reinforced by a good correspondence between the highest abundances of this species in the dinoflagellate cyst record and a period of increased annual precipitations, as indicated by climate data directly measured and also inferred from independent proxies (tree ring reconstruction). Data presented in the present study based on recent cyst distributions (Figures 7, 9B) supports these interpretations. Besides, they are in line with various studies worldwide (e.g., Amorim and Dale, 2006;Leroy et al., 2013;Ribeiro et al., 2016;García-Moreiras et al., 2018;de Vernal et al., 2020) that indicated a strong correlation between SST and L. polyedra abundances, supporting the utility of this cyst species as an indicator of warm and stratified environments.
The second group, characterized by high abundances of G. catenatum, the single species with the highest relative abundances in the study area, seems to reflect a different ecological niche from other autotrophs. It plotted the most positive on RDA1 with a strong association with deeper waters (Figures 7, 9B). Latitudinally, it was particularly abundant (percentages) in the intermediate transects (i.e., red dots in Figure 9A).
The upwelling pattern in the study region is characterized by two upwelling frontal zones that tend to follow the depth contours: a narrow coastal band and further offshore close to the shelf break (Relvas et al., 2007), which coincides with the zones of high frontal probability represented in Figure 1C. The high relative abundances of G. catenatum in deeper samples following the 100 m isobath may be linked to these upwelling fronts. Indeed, Moita et al. (1998) reported high abundances of the vegetative stage at mid-shelf, related to the offshore displacement of blooms. More recent observations recorded the presence of G. catenatum in the inshore side of upwelling plumes where currents may be reduced and G. catenatum may develop close to the core of the upwelling plume without being advected away (Moita et al., 2003). This is consistent with the fact that coastal outbreaks of this species typically occur during the first late-summer/autumn downwelling events, when these offshore populations may be advected to the coast (e.g., Bravo et al., 2010;Díaz et al., 2019). High abundances of G. catenatum cysts may thus reflect the association of this species with upwelling fronts in W Iberia. Several authors have proposed the importance of chain formation as an adaptive strategy in upwelling areas (Fraga et al., 1988;Smayda and Reynolds, 2003). More recently, Smayda (2010) emphasized that although chain formation may be an advantage in turbulent waters associated with upwelling, chain-forming dinoflagellates are the exception and not the rule in upwelling systems and probably other factors such as temperature tolerance play an important role in species selection.
The previous interpretations of the distribution of dinoflagellate cyst assemblages in surface sediments and its comparison with present-day environmental data suggest that, in the study area, cysts are mainly reflecting the hydrographic conditions observed during the summer upwelling season. This season corresponds to the period of highest primary productivity (Moita, 2001). Interestingly, Abrantes and Moita (1999) investigated the planktonic communities of diatoms and coccolithophores along the Portuguese coast under winter (non-upwelling) and summer (upwelling) conditions and how they compared with the sediment record. Their results indicated that the sediment assemblages reflected mainly the water column distribution observed during the upwelling season.
Regarding the distribution of P. americanum, the weight in the ordination analysis is low and interpretations should be considered with caution. However, our results suggest that similarly to G. catenatum this species may prefer transitional waters between less productive and stratified environments (lower right quadrant, Figure 9B), and more productive environments associated to upwelling fronts (left side of RDA, Figure 9B). Existing studies on the distribution of this cyst and environmental gradients do not reflect a distinct distribution from other heterotrophs, indicating that its distribution is associated -but not restricted to -productive waters near upwelling cells (Dale, 1996;Marret and Zonneveld, 2003;Ribeiro and Amorim, 2008;Zonneveld et al., 2013). Zonneveld et al. (2013) also refer the preference for fully marine environments. Concordantly, the RDA indicates the distribution of P. americanum associated with more saline (and deeper) environments in the study area (Figure 9B), if compared to other heterotrophs.
Finally, a few considerations on the distribution of HAB species and the role of benthic cysts in HAB formation. In the present study, we recorded cysts of three dinoflagellate species that may be potentially toxic and produce Harmful Algal Blooms (HAB) in the Atlantic Iberian margin, namely G. catenatum, L. polyedra and P. reticulatum. Very high percentages of G. catenatum cysts (Figure 7) reflect a long history of blooms along the NW coast of Portugal. However, the almost absence of viable cysts in surface sediments is additional evidence supporting the hypothesis that this species does not rely on benthic cyst beds for seeding planktonic blooms (Fraga et al., 1993;Moita et al., 1998;Bravo et al., 2010). Although blooms of G. catenatum have been studied in W Iberia for more than 40 years, it is still not fully understood which factors other than life cycle traits (physical or biological) may be involved in the initiation and successful development of blooms.
The yessotoxin producers L. polyedra and P. reticulatum reached high relative abundance in the southern transects (∼20 and 40% respectively). However, as for G. catenatum, the percentage of cysts with cell contents was very low suggesting that in the study area there is not a build-up of significant cyst beds of these two species. This interpretation is supported by a previous cyst survey along the coast of Portugal, where the main distribution center for L. polyedra cysts occurred in the south coast (Amorim et al., 2004). Regarding P. reticulatum, the absence of a good correspondence between reports of vegetative planktonic populations and the cyst record has been previously recognized (Dale, 1976). Several hypotheses were put forward which included the mismatch between phytoplankton studies and the occurrence of planktonic populations of P. reticulatum, and/or a species specific high cyst:motile cell ratio during the encystment process. So far, no harmful events have been associated with P. reticulatum in W Iberia. However, the high abundance of cysts in sediments prompts investigation on the ecology and toxicity of regional strains.

CONCLUSION
The study of the distribution of dinoflagellate cyst assemblages in 51 surface sediment samples off Aveiro-Figueira da Foz revealed marked land-sea and latitudinal gradients. Summer coastal upwelling was identified as the main ecological gradient driving dinoflagellate cyst assemblages in the study area.
Redundant analysis (RDA) on relative abundances revealed the existence of two main environmental regimes. One included the southern and offshore sites and was characterized by higher SST, lower BTT and lower primary productivity (lower CHL), suggesting warm, stratified and less productive environments. The other regime, which included the northern sector and inshore sites, was characterized by lower SST, higher BTT and enhanced primary productivity (CHL), suggesting the influence of upwelling.
Three main ecological signals were identified in the dinoflagellate cyst assemblages: i) The heterotroph signal as the main upwelling signal; ii) The dominance of P. reticulatum and L. polyedra signal, indicative of warm stratified conditions, possibly reflecting transitional environments between more active inshore upwelling and warmer offshore waters; iii) The G. catenatum signal, the main signal for the presence of mid-shelf upwelling fronts.
The almost absence of viable cysts of the toxic and potentially toxic species G. catenatum L. polyedra and P. reticulatum in surface sediments, suggests that in the study area, for these three species, there is no build-up of significant cyst beds and thus planktonic populations must depend on other seeding processes.
These results are the first detailed distributions of modern dinoflagellate cysts in the NW Iberian Atlantic margin (off Portugal). Despite the uncertainties related to sediment and cyst transport, post-depositional processes and the possible disparity between the time scale of the cyst record and the time scale of the environmental data used, cyst distributions show a fairly compelling coincidence with hydrographic features during the summer upwelling season in the study area. This means that modern dinoflagellate cyst assemblages are reflecting water column characteristics and may be used as supporting evidence for the interpretation of stratigraphic cyst records and reconstruction of past marine ecosystems in W Iberia.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

AUTHOR CONTRIBUTIONS
AA participated in the oceanographic campaign to obtain surface sediment samples and water depth data. AA, PO, and AO planned the research. Sediment texture analyses and statistics were performed by AO and AS. Oceanographic data collection and analysis was performed by PO. Dinoflagellate cyst analyses and microscopic observation and photography were conducted by IG-M and AA. Multivariate statistical analyses were performed by IG-M. All authors participated in data discussion, figure elaboration, and manuscript writing.

FUNDING
This work was a contribution to HABWAVE project LISBOA-01-0145-FEDER-031265, co-funded by EU ERDF funds, within the PT2020 Partnership Agreement and Compete 2020, and national funds through Fundação para a Ciência e Tecnologia, I.P.(FCT, I.P.) and FCT, I.P. under the project UIDB/04292/2020, and also to AQUIMAR project MAR2020 No. MAR-02.01.01-FEAMP-017. This work was also supported by funding from the European Union's Horizon 2020 Research and Innovation Programme under grant agreement N 810139: Project Portugal Twinning for Innovation and Excellence in Marine Science and Earth Observation -PORTWIMS. IG-M was supported by a postdoctoral fellowship from Xunta de Galicia, Spain (ref. ED481B-2019-074, 2019).

ACKNOWLEDGMENTS
Teresa Moita (CCMAR -Centro de Ciências do Mar, Universidade do Algarve, Portugal) is greatly acknowledged for her invaluable help in the preparation of the distribution maps. We also thank Melissa Hatherly for her great assintance in the processing of sediment samples for cyst analyses. We are grateful to all the researchers, students and the crew of "N. O. Auriga", who were most helpful in ensuring the success of the cruise, and especially to Raquel Melo for her technical assistance on-board and geographical data base management. We also thank the reviewers for their comments which greatly improved the manuscript. This study has been conducted using E.U. Copernicus Marine Service Information and data obtained from the Physical Oceanography Distributed Active Archive Centre (PO.DAAC) at JPL (Jet Propulsion Laboratory).