Allochthonous matter: an important factor shaping the phytoplankton community in the Baltic Sea

Abstract It is well-known that nutrients shape phytoplankton communities in marine systems, but in coastal waters allochthonous dissolved organic matter (ADOM) may also be of central importance. We studied how humic substances (proxy of ADOM) and other variables influenced the nutritional strategies, size structure and pigment content of the phytoplankton community along a south–north gradient in the Baltic Sea. During the summer, the proportion of mixotrophs increased gradually from the phosphorus-rich south to the ADOM-rich north, probably due to ADOM-fueled microbes. The opposite trend was observed for autotrophs. The chlorophyll a (Chl a): carbon (C) ratio increased while the levels of photoprotective pigments decreased from south to north, indicating adaptation to the darker humic-rich water in the north. Picocyanobacteria dominated in phosphorus-rich areas while nanoplankton increased in ADOM-rich areas. During the winter–spring the phytoplankton biomass and concentrations of photoprotective pigments were low, and no trends with respect to autotrophs and mixotrophs were observed. Microplankton was the dominant size group in the entire study area. We conclude that changes in the size structure of the phytoplankton community, the Chl a:C ratio and the concentrations of photoprotective pigments are indicative of changes in ADOM, a factor of particular importance in a changing climate.


I N T R O D U C T I O N
Phytoplankton communities are shaped by bottom-up and top-down factors (Chisholm, 1992;Kiørboe, 1993;Mousing et al., 2014), and their community properties are in turn important for aquatic food web functioning (Havens, 1998;Dahlgren et al., 2010). Small phytoplankton cells possess advantages over larger cells for resource acquisition, growth rate and photosynthetic rate under low-nutrient and light-limiting conditions (Grover, 1989;Raven, 1998). Dominance of small phytoplankton results in more internal trophic levels in the food web and thus reduced food web efficiency (Legendre and Rassoulzadegan, 1995). Larger phytoplankton cells dominate in turbulent and nutrient-rich environments, inducing shorter food chains with higher food web efficiency (Legendre and Rassoulzadegan, 1995). Because of the interrelationship between nutrients and temperature in aquatic systems, changes in the size structure of the phytoplankton community are often caused by concurrent variations of these factors (Agawin et al., 2000;Mousing et al., 2014).
Due to the interaction with the heterotrophic microbial food web, the size distribution and function of the phytoplankton community is not only influenced by factors like nutrient availability and temperature. In coastal areas and semi-enclosed seas like the Baltic Sea, the influence of allochthonous dissolved organic matter (ADOM) can be strong, fueling the heterotrophic microbial food web (Andersson et al., 2015;Figueroa et al., 2016). Under such conditions mixotrophic phytoplankton may also be enhanced, since they feed, e.g. on heterotrophic bacteria (Andersson et al., 1989). Changes in ADOM input can be influenced both by climate-related factors, such as temperature, precipitation and hydrology or by alterations in acidity or land-use activity (Evans et al., 2006;Erlandsson et al., 2008). Climate change projections indicate that both the temperature and precipitation will increase in northern Europe. The river inflow of freshwater to the Baltic Sea will thus increase, causing decreased salinity and increased concentrations of ADOM (Meier et al., 2012;Andersson et al., 2015). How an increase in ADOM will affect phytoplankton structure in the Baltic Sea is poorly understood.
The brownish color of ADOM attenuates light in the water and hampers primary production (Thrane et al., 2014;Seekell et al., 2015). It also leads to changes in the cellular concentration of pigments in phytoplankton (Falkowski and Raven, 2007). Phytoplankton respond to decreasing light by increasing their concentrations of chlorophyll a (Chl a) and accessory photosynthetic pigment content in order to harvest as much light as possible (MacIntyre et al., 2002;Behrenfeld et al., 2005).
Under poorer light conditions smaller cell sizes with large S/V ratios are also beneficial, because the lightharvesting pigments are positioned closer to the cell surface and are more evenly distributed in the cells (Raven, 1998;Kirk, 2011). In high light conditions the photoprotective pigments in phytoplankton will increase and larger cells become promoted. Furthermore, the taxonomic composition and nutrient availability influence pigment concentrations, e.g. as nitrogen is essential in the synthesis of pigment-proteins complexes (Dubinsky and Stambler, 2009;Edwards et al., 2015;Spilling et al., 2015). Even though such physiological acclimatization's are known from laboratory studies, it is not known whether similar responses occur in natural systems.
The objective of this study was to elucidate how bottom-up factors, such as ADOM, nitrogen, phosphorus and temperature, influence properties in the phytoplankton community in the semi-enclosed Baltic Sea. Samples were collected along a south-north gradient during a summer and a winter-spring campaign, to elucidate factors governing the general size structure, nutritional strategy and pigment content of the phytoplankton community. The following hypotheses were tested: (i) high phosphorus and nitrogen concentrations lead to large phytoplankton cell sizes, (ii) high ADOM concentrations promote mixotrophs and drive phytoplankton to increase their chlorophyll a to carbon ratio (Chl a:C) and (iii) low ADOM drives phytoplankton to increase their levels of photoprotective pigments. This study contributes to the general understanding of marine ecosystem system function as well as consequences of climate change, since in the Baltic Sea the ADOM and nutrient concentrations are expected to change during the course of the next 100 years (Andersson et al., 2015).

Field sampling
Spatial and temporal variation in the phytoplankton community and physicochemical variables were investigated during a late summer (23-25 August 2011) and a winter-spring (16-19 March 2012) sampling campaign. Fourteen stations were sampled along a south-north gradient in the Baltic Sea (Fig. 1, Table S1, Supplementary material online). Ice cover was present in the Bothnian Bay in March, and reached 50 cm in the northernmost part of the basin. Samples were collected from a Ferry Box system, installed on a cargo ship (TransPaper) traveling between Gothenburg (Sweden) and Kemi (Finland). Temperature and salinity were measured using an SBE Temp 38/SBE TSG 45 sensor. Water samples were handled and preserved on board within 2 h of collection.

Chemical analyses
For the analysis of total phosphorus (Tot P) and total nitrogen (Tot N), 50 ml of unfiltered water were frozen (-20°C) in Falcon tubes and stored until analysis. Measurements were carried out using a Braan and Luebbe TRAACS 800 autoanalyzer, following standard analytical methods (Grasshoff et al., 1983). Tot N and Tot P were considered to indicate the carrying capacity of the system, and also to reflect the nutrients available to the phytoplankton. This assumption is based on results of other studies performed in the Baltic Sea, where Tot N and Tot P concentrations were shown to be positively correlated with inorganic nutrients (data not shown). Furthermore, nutrient turnover has been shown to be relatively rapid in aquatic systems: Kress et al. (2005) showed that in the Mediterranean Sea the nutrient turnover time was <1 h when nutrient limitation prevailed, while when nutrients were replete the nutrient turnover time increased to 3-4 days. In the Baltic Proper, Tot P turnover related to the phytoplankton production was estimated to be 3.2 days (Håkansson and Bryhn, 2008).
Humic substances (HSs) were analyzed and used as a proxy for ADOM. Samples of 100 ml were kept in amber-colored glass bottles in the dark at 4°C and measured with a Perkin Elmer LS 30 fluorometer at 350/450 excitation/emission wavelengths. Calibration standards  were prepared from quinine dihydrogen sulfate dihydrate in 0.05 M sulfuric acid (Hoge et al., 1993;Wedborg et al., 1994). Sulfuric acid (0.05 M) was used as a blank.

Pigments
Samples (200-800 ml) were filtered onto 25 mm GF/F filters, frozen in liquid nitrogen and stored at -80°C. To extract the pigments, filters were placed in 3 ml 95% acetone containing vitamin E as an internal standard, sonicated in an ice bath for 15 min and extracted at 4°C overnight in the dark. After extraction, the samples were filtered through a 0.2 µm Teflon syringe filter into vials, injected into a Schimadzu LC-10ADVP (High Performance Liquid Chromatography: HPLC) and analyzed according to van Heukelem and Thomas (van Heukelem and Thomas, 2001), with slight modifications for the local conditions (Schlüter et al., 2014) to identify and quantify the phytoplankton pigments. The photoprotective pigment index (PI), a measure of the physiological state of the phytoplankton community, was calculated according to Moreno (Moreno et al., 2012): Phytoplankton For the analysis of nanoplankton and microplankton, two 50 ml samples were preserved with 2% acidic Lugol's solution, settled in sedimentation chambers for 24-48 h and counted at ×100 (microplankton) and ×400 (nanoplankton) magnification using a Leica DM IRB inverted microscope (Utermöhl, 1958). The nutritional strategy of the phytoplankton (autotrophic, heterotrophic or mixotrophic) was classified according to Olenina (Olenina et al., 2006). To study the size structure, the phytoplankton cells were divided into four size categories: <2 µm (picoplankton), 2-10 µm (ultraplankton), 10-20 µm (nanoplankton) and >20 µm (microplankton), based on measurements of the longest cell axis. Tightly connected cells of filamentous cyanobacteria were grouped into the >20 µm fraction.
Samples for the analysis of picocyanobacteria were fixed with glutaraldehyde (2% final concentration) and filtered (2-5 ml) onto 0.6 µm black polycarbonate filters. The samples were counted using an epifluorescence microscope (Nikon Eclipse TE 2000-U) at ×1000 magnification, using green excitation light (510-560 nm, emission wavelength > 590 nm). At least 200 cells per sample, in 20 randomly positioned view fields, were counted.
The phytoplankton biovolume was calculated according to the size and geometry of the cells (Olenina et al., 2006). Cell carbon was calculated from the biovolume (Menden-Deuer and Lessard, 2000), and carbon biomass concentrations from cell abundance and cell carbon. The Chl a:C ratio indicated the chlorophyll a content in the phytoplankton cells, and was calculated by dividing the Chl a concentrations with the phytoplankton carbon biomass concentrations.

Statistical analyses
Physicochemical and biological variables were tested for normality using the Shapiro-Wilk test and homogeneity by Levene's test. Variables that did not fulfill for a normal distribution were either logarithmically transformed or arcsin-square root transformed (contribution of mixotrophs, autotrophs to the total biomass and Chl a:C ratio, photoprotective PI). The Student's t-test was used to test differences in phytoplankton total biomass, phytoplankton size structure, chlorophyll a and physicochemical variables between the summer and winterspring sampling seasons. To explore relationships between physicochemical and biological variables, Pearson's rank correlation was performed. Principal component analysis (PCA) was used to show the main patterns of phytoplankton size structure, nutritional strategies and pigment content in relation to potentially explanatory variables during the summer and winterspring samplings. Non-metric multidimensional scaling (NMDS) plot based on Bray-Curtis similarity matrix was performed to visualize differences in phytoplankton composition between stations and seasons. To test similarities in phytoplankton composition between sampling periods, ANOSIM analysis was conducted. Station 1 was excluded from analyses because the values differed substantially from the other stations as a result of the strong influence of the North Sea. Data analyses were performed in SPSS Statistics 23, Primer 6 and Canoco 5.

R E S U L T S
The physicochemical variables in general showed spatial and temporal variations (Table S1, supplementary material online). The salinity was 18-20 in the Kattegat and decreased to 2-3 in the Bothnian Bay. The average temperature was 17 and 2°C in the summer and winter-spring, respectively. Tot P showed a decreasing trend from south to north during both summer and winter-spring, while Tot N showed a similar spatial pattern only during the summer ( Fig. 2A and B). The Tot P concentrations were relatively similar during both periods (Student's t-test: t = 0.056, df = 24, P = 0.956), while the Tot N concentrations were higher during summer (Student's t-test: t = 6.786, df = 24, P < 0.001). The concentrations of HSs increased from south to north with highest values in the Bothnian Bay (Table S, supplementary material online, Fig. 2C). Overall HSs did not differ between seasons (Student's t-test: t = 0.645, df = 24, P = 0.525). Pearson correlation analysis showed that many of the physicochemical variables were correlated to each other (Table S2). For example, HSs were negatively correlated to Tot P, salinity and temperature during both the summer and the winterspring samplings.

Distribution of Chl a
The Chl a concentrations were higher in summer than in winter-spring (Fig. 2D, Student's t-test: t = 6.668, df = 24, P < 0.001). During summer the concentrations were~1 µg L −1 at most of the stations. In winter-spring there was an overall decreasing trend from the southern Baltic Proper (0.8 µg L −1 ) to the northern Bothnian Bay (0.1 µg L −1 ). However, it is notable that during both seasons highest concentrations were observed at the entrance of the Gulf of Bothnia, while the lowest occurred in the Kattegat.

The Chl a:C ratio and PI
The Chl a:C ratio increased from south to north during the summer, while the opposite trend was observed during the winter-spring (Fig. 2E). However, the lowest Chl a:C ratio was observed at station 1 (Kattegat), during both seasons. In the summer, the Chl a:C ratio showed a positive correlation with HSs and a negative correlation with salinity, temperature, Tot N and Tot P (Table I). In  the winter-spring, the Chl a:C ratio showed a negative correlation with HSs and a positive correlation with temperature, salinity and Tot P (Table I). Overall, our data showed a 10-fold difference in the Chl a:C ratio.
The PI was ten times higher in the summer than in the winter-spring (Fig. 2F). In the summer, the PI showed a hump-shaped pattern from south to north, with the highest values in the northern Baltic Proper (stations 4-7). During the winter-spring relatively low values were observed. The PI correlated negatively with HSs and positively with salinity, temperature, Tot N and Tot P during the summer (Table I).

Phytoplankton biomass and nutritional strategy
During the summer the total phytoplankton biomass increased from the Kattegat station to the station located between the islands of Öland and Gotland (Fig. 3A). Northwards the phytoplankton biomass gradually decreased. The contribution of autotrophs was highest in the south and decreased gradually towards the north, while the proportion of mixotrophs increased northwards ( Fig. 3C and E). Heterotrophs constituted an insignificant part of the biomass, and therefore they were not taken into account in the analysis. The total phytoplankton biomass was approximately five times higher in summer than in winter-spring (Fig. 3B) (Student's t-test: t = 7.634, df = 24, P < 0.001). In winter-spring, peaks of total pytoplankton biomass were observed in Kattegat, stations located close to islands (4, 7 and 11) and at the northernmost station in the Bothnian Bay (Fig. 3B). No specific trend in relation to nutritional strategy was observed during the winterspring ( Fig. 3D and F).
In the summer the total biomass was positively correlated with salinity, temperature, Tot N and Tot P, and negatively correlated with HSs (Table I). In winterspring a positive correlation with Tot N was found (Table I).

Size distribution and composition of phytoplankton
In the summer picoplankton (picocyanobacteria) was the dominant size group in all basins, while their relative contribution decreased gradually towards the north (Fig. 3A). Microplankton showed a similar spatial pattern, constituting~25% of the biomass in the Baltic Proper and Bothnian Sea but only 4% in the Bothnian Bay. Ultraplankton and nanoplankton increased in importance towards the north and formed~50-60% of the biomass in the northernmost stations in the Bothnian Bay (Fig. 3A). Picoplankton and microplankton were positively correlated with salinity, temperature, Tot N and Tot P and negatively with HSs (Table I). In contrast, the ultraplankton and nanoplankton showed a positive correlation with HSs and a negative correlation with other physicochemical variables. The second largest groups were Dinophyceae in the Kattegat, Cyanophyceae (colony-forming and/or filamentous cyanobacteria) in the Baltic Proper and Bothnian Sea, and Chlorophyceae in the Bothnian Bay (Fig. 4A).
During the winter-spring period, microplankton (Dinophyceae) dominated the phytoplankton community at all stations, forming on average~50% of the biomass (Figs 3B and 4B). Picoplankton (picocyanobacteria) was the second most dominant fraction, and lowest concentrations were observed in the ice-covered Bothnian Bay. Ultraplankton and nanoplankton generally constituted <10% of the phytoplankton biomass. Positive correlations between picoplankton and Tot N and Tot P were observed, while microplankton only showed a positive correlation with Tot N (Table I). Picoplankton Table I: Pearson correlation coefficients between total phytoplankton biomass, different size groups (<2, 2-10, 10-20 and >20 µm), autotrophs (AU), mixotrophs (MX), chlorophyll a (Chl a), chlorophyll a: carbon ratio (Chl a:C), photoprotective pigment index (PI) and physicochemical variables [salinity, temperature (Temp), total nitrogen (Tot N), total phosphorus (Tot P), humic substances (HSs)]. (Studen's t-test: t = 10.944, df = 24, P < 0.001) and ultraplankton biomass (Student's t-test: t = 4.065, df = 24, P < 0.001) showed significant differences between the seasons. NMDS ordination analysis showed that the phytoplankton community structure differed in the summer and the winter-spring (Fig. 5, ANOSIM global R = 0.746, P < 0.01). Furthermore, the phytoplankton community composition showed more differences in the south-north gradient during the winter-spring than during the summer, as indicated by the NMDS plot.

Factors shaping the phytoplankton community
In the PCA the first two axes explained 79.9% of the variance of the summer data, and 70.1% of the variance of the winter-spring data (Table II). In summer the highest principal component loadings were constituted by HSs (positive loading) and Tot P and salinity (negative loadings), while in winter-spring Tot N was the most important factor (positive loadings). The PCA indicated that in the summer HSs was a driver for the occurrence of ultra-and nanoplankton (partly consisting of mixotrophs), and high chlorophyll content in the phytoplankton cells, i.e. high Chl a:C ratio (Fig. 6A). This environmental condition was predominant at the most northerly stations. Autotrophic pico-and microplankton and the photoprotective pigments index were promoted by high Tot P and salinity, which were the prevalent environmental conditions at the more southerly stations. The PCA indicated that during winterspring high Tot N concentration promoted e.g. total phytoplankton biomass and Chl a (Fig. 6B). Fig. 3. Phytoplankton biomass divided into four size groups ( <2, 2-10, 10-20 and >20 µm) (A and B), relative proportion (%) of autotrophs (C and D) and mixotrophs (E and F) during summer and winter-spring sampling periods.

D I S C U S S I O N
Our study indicates that ADOM, in addition to the more well-known factors (nitrogen, phosphorus and salinity), governs the large-scale patterns of phytoplankton nutritional strategy, size structure and pigment content in the Baltic Sea. During the summer, the relative proportion of autotrophic picoplankton was positively correlated with Tot P. This was somewhat unexpected because, applying ecological assumptions (Legendre and Rassoulzadegan, 1995), low-nutrient concentrations should lead to a dominance of smaller cells. Smaller phytoplankton cells have higher affinity for nutrient uptake and higher photosynthetic efficiencies than larger cells (Raven, 1998), but we did not find such responses at the ecosystem scale. Large autotrophic phytoplankton cells, >20 µm, showed a similar response as picoplankton. The decrease of the relative proportion of autotrophic pico-and microphytoplankton towards the north during the summer may thus have been caused by low-nutrient availability, deeper mixing and higher concentrations of HSs, which reduce the light availability in the water column (Kuosa, 1991). It is well-known that the northernmost part of the Baltic Sea (Bothnian Bay) is strongly phosphorus limited, while the Baltic Proper and the Swedish west coast generally are nitrogen-limited systems (Graneli et al., 1990;Andersson et al., 1996;Tamminen and Andersen, 2007).
Ultraplankton and nanoplankton were positively influenced by HSs during the summer. This indicates that these size groups have a different functional ecology compared with picocyanobacteria and larger phytoplankton. In fact, many of the species belonging to these groups are potentially mixotrophs, such as chrysophyceans (e.g. Dinobryon spp. and unidentified pigmented flagellates), dinophyceans and prymnesiophyceans (Chrysochomulina spp.). However, some purely autotrophic chlorophyceans (e.g. Oocystis spp.) also occurred in these size groups. Mixotrophs are likely to be promoted by HSs, as some of them feed on bacteria, which are capable of utilizing part of the ADOM (Andersson et al.,   2015; Figueroa et al., 2016). These potential mixotrophs mainly occurred in the north, where the concentrations of HSs were higher and conditions more similar to humic lakes than the open sea. Our results are in agreement with studies in freshwater systems, where brownwater lakes have been shown to have higher abundances of mixotrophs compared with clear-water lakes (Bergström et al., 2003;Saad et al., 2013). This finding is less commonly recognized in marine systems, because usually only the coastal zone is markedly affected by ADOM. However, because of the high influence of river inflow in the north of the Baltic Sea, a strong structuring role of HSs was detected even in the offshore waters. Previous studies in the Baltic Sea have shown that mixotrophs are most abundant during the summer because of low-nutrient levels, which should promote phytoplankton with a diverse (i.e. mixotrophic) feeding mode and low salinity (Andersson et al., 1996;Hajdu et al., 1996;Dahl et al., 2005). However, our study indicates that low phosphorus and salinity and high humic concentrations selectively promoted this group in the summer. This study shows that unicellular picocyanobacteria can be a dominant component of the phytoplankton community, forming 40-90% of the total phytoplankton biomass in the summer. The observed abundances of picocyanobacteria are comparable to previous studies in the Baltic Sea (Andersson et al., 1996;Hajdu et al., 2007). Stal et al. (1999) reported that 65% of the phytoplankton biomass in the Baltic Proper during late summer was composed of picoplankton, while the second most dominant group was nitrogen-fixing cyanobacteria (Aphanizomenon spp. and Nodularia spp.). Hajdu et al. (2007) showed that during the decline phase of cyanobacterial blooms in late summer, small diatoms, nanoflagellates, unicellular and colony-forming picocyanobacteria increase in abundance (Hajdu et al., 2007). From satellite image analysis we know that our summer sampling campaign was performed 1-2 weeks after the decline of extensive cyanobacterial blooms in the Baltic Proper (Hansson and Öberg, 2011), proving that the sampling was performed during the post-bloom period. In agreement with some earlier studies conducted during late phases of the summer bloom (Albertano et al., 1997;Stal et al., 1999), we found that unicellular picocyanobacteria and filamentous nitrogen-fixing cyanobacteria were the dominant phytoplankton groups. Thus, the general ecological relationship between nutrient availability and plankton size structure (Legendre and Rassoulzadegan, 1995) can be confounded by complex species interactions and successional timing during the year.
During winter-spring, autotrophic dinoflagellates dominated the phytoplankton community. This group has been reported to become more dominating during the spring bloom, possibly due to stronger vertical stratification induced by warmer seawater temperatures, which favor motile dinoflagellates over immotile diatoms (Wasmund and Uhlig, 2003). However, the present data-set does not directly support this theory, since  dinoflagellates were the dominant group also in the coldest ice-covered basin, the Bothnian Bay.
Our study showed that the Chl a:C ratio and PI in the phytoplankton community was governed by a combination of different physicochemical variables. Thus phytoplankton inhabiting dynamic and changing environments, e.g. estuaries, may have a more varied pigment content than phytoplankton in the open sea. Chl a is commonly used as a proxy for total phytoplankton biomass in both monitoring and research programs (Andersson et al., 1996;HELCOM, 2013). However, our study shows that Chl a is a poor proxy for total phytoplankton biomass in the Baltic Sea because the Chl a:C ratio can vary by a factor of 10. We interpret the higher Chl a:C ratios in the northern Baltic during the summer as an effect of low light as a result of increased water color, which forced the phytoplankton to increase their chlorophyll a content to be able to capture light and sustain photosynthesis. Moreover, the increase in the Chl a:C ratios from south to north may also be caused by decreasing nitrogen limitation, as nitrogen is a major component of Chl a. On the other hand, it may appear strange that the Chl a:C ratio was highest where mixotrophic phytoplankton were most common. In some mixotrophs the chloroplasts are rudimentary and photosynthesis is used just as a survival mechanism when particulate food is scarce (e.g. Andersson et al., 1989). However, a large range of mixotrophic ecotypes occurs in aquatic systems, from almost purely autotrophic to almost purely heterotrophic (Jones, 2000). Furthermore, as the mixotrophs in our study formed at most only 10% of the total biomass, it is unlikely that they have a strong influence on the Chl a: C of the entire phytoplankton community.
The PI in the phytoplankton community was relatively high in the southern Baltic where concentrations of HSs were low. We interpret this to be a direct effect of a higher light intensity in the areas where the Tot P was also relatively high. In the northern region the concentrations of HSs were high, causing significant attenuation of the photosynthetically active radiation in the water column (Hoikkala et al., 2015). Previous studies have shown that the PI is governed by light conditions and that adaptation can occur in a very short period of time (Moisan et al., 1998), and that the photoprotective pigments increase in non-turbid clear waters and in surface water (Brunet et al., 1993;Riegman and Kraay, 2001). Certain cyanobacteria, e.g. Synechococcus spp., have some of the highest variations in the zeaxanthin:Chl a ratio under changing environmental conditions (Veldhuis and Kraay, 1990;Schlüter et al., 2000). Since picocyanobacteria dominated the phytoplankton community, they probably contributed to the observed geographical PI changes.

C O N C L U S I O N S
In conclusion, we found that in the summer ADOM was one of the major factors governing the phytoplankton community, while Tot N played a major role in the winter-spring. The strong influence of ADOM in the northern region of the Baltic Sea (Bothnian Bay) favored ultraplankton and nanoplankton, which partly is formed by mixotrophic phytoplankton, capable of feeding on bacteria and other particulate matter. Furthermore, the brownish color of the ADOM causes shading in the seawater, which seemed to drive phytoplankton to increase their chlorophyll a content in the cells and to reduce the photoprotective pigments. Regional climate change projections indicate that within a 100-year period precipitation will increase in northern Europe, which will lead to increased river inflow of ADOM and a freshening of the Baltic Sea. We suggest that the observed changes in the phytoplankton community from north to south partly mirror future changes. The present northerly phytoplankton community structure in the future will transfer somewhat further south in the Baltic due to climate change.

S U P P L E M E N T A R Y D A T A
Supplementary data can be found online at Journal of Plankton Research online

A C K N O W L E D G M E N T
We thank the captain and the crew of the merchant vessel TransPaper, and TransAtlantic AB for hosting us during the sampling seasons. Dr Iveta Jurgensone, Latvian Institute of Aquatic Ecology, is gratefully acknowledged for phytoplankton analysis, and Merete Allerup, DHI Water & Environment, Denmark, for assistance with the analysis and interpretation of HPLC samples.