Scale-Specific Patterns of the Production of the Charophyte Chara aspera in the Brackish Baltic Sea: Linking Individual and Community Production and Biomass Growth

An understanding of the nature of scale-dependence in macroalgal production and to quantify how primary production is incorporated into autotrophic biomass requires an assessment of primary production at multiple scales. In this study we experimentally estimated seasonal variability in individual photosynthesis, community production and growth in biomass of the charophyte Chara aspera in the Baltic Sea together with the key environmental variables known to regulate the production of macroalgae. The experiments suggested that the production and growth of C. aspera was defined largely by light and temperature at all studied scales. However, the algal production at the same light levels was systematically lower at the community than an individual level, suggesting the importance of shelf-shading and below-ground processes in natural plant assemblages. Moreover, the observation scale also defined response types between the environment, plant production and growth. This implies that the patterns of variability of the production and growth of macrophytes should always be interpreted in the context of scale and any multiscale model development should involve experimental validation at all important scales.


INTRODUCTION
Primary production is important to all ecosystems as it supplies energy to other living organisms, thereby forming the foundation of food webs. Macroalgae are the dominant primary producers in coastal temperate marine ecosystems, representing some of the highest recorded per-area basis production values (Field et al., 1998;Gattuso et al., 2006). As macroalgae contribute significantly to total global system carbon production (Hill et al., 2015) it is important to understand better how macroalgal primary production links to growth in biomass and how different environmental factors affect these patterns.
Traditionally, the primary production of macroalgae has been quantified using detached pieces of algal thally or whole plants Gómez et al., 1997) and more recently, the production of macroalgae has been measured directly in natural assemblages (Binzer and Middelboe, 2005;Tait and Schiel, 2010). While the first approach enables quantification of accurate response functions between different environmental variables and primary production (e.g., irradiance-net photosynthesis curves), the latter technique allows the scaling up of primary production in natural conditions. The estimates of primary production are highly dependent on the scale of observation with photosynthetic production being higher at individual than at community levels (Binzer et al., 2006;Pärnoja et al., 2014). Likewise, scale affects how macroalgae respond to changes in ambient conditions. The photosynthesis of individual plants often shows a typical light saturation curve (Lobban et al., 1985), whereas macroalgal photosynthesis in assemblages increases in a roughly linear fashion up to maximum levels of natural irradiance (Middleboe and Binzer, 2004). Changes in biomass integrate the net photosynthetic fixation of inorganic carbon into autotrophic biomass as well as the processes of natural loss of biomass and thereby permits estimation of large-scale variation in realized production over time. Expectedly, when compared to the patterns of photosynthesis, another set of combinations of environmental factors can influence change in macroalgal biomass or the same relationships can have very different functional shapes. In this respect, laboratory experiments have shown species-specific responses of macroalgae to environmental change in which macroalgal growth may increase, decrease or remain unchanged with increasing rates of plant photosynthesis (de los Santos et al., 2000;Israel and Hophy, 2002;Kim et al., 2016). Thus, understanding the nature of scale-dependence in macroalgal production can help to elucidate some of the possible controls by which primary production is incorporated into autotrophic biomass. Despite this importance, knowledge of such controls is limited as few studies have been undertaken to date that quantify the patterns of macroalgal production at different scales with most being performed under laboratory conditions (e.g., de los Santos et al., 2000;Israel and Hophy, 2002;Binzer and Middelboe, 2005;Pärnoja et al., 2014;Kim et al., 2016). Due to artificially induced conditions, as well as the use of separate phytoelements (thallus pieces, leaves), however, laboratory incubations often fail to characterize primary production to that which occurs in nature (Binzer and Middelboe, 2005;Binzer et al., 2006). Primary production is controlled by both abiotic and biotic factors (Field et al., 1998;Hauxwell et al., 2003) but abiotic constraints are often dominant in natural environments (Hill et al., 2015). The availability of light and nutrients plays a crucial role in regulating primary production of macroalgae (Field et al., 1998;Hauxwell et al., 2003;Binzer et al., 2006). Light limitation may arise from the amount of irradiance arriving at the sea surface, the optical characteristics of water and self-shading within algal assemblages. The first two processes vary over multiple temporal and spatial scales and define the availability of light field above underwater canopies and the maximum photosynthetic rates of macroalgae (Kirk, 1994;Anthony et al., 2004). The third process, shelf-shading, is one of the most critical biological limitations because it sets an actual threshold for realized photosynthesis along vertical light gradient within macroalgal assemblages (Binzer et al., 2006;Tait and Schiel, 2010). The effect is stronger with increasing algal densities (Kotta et al., 2008;Paalme et al., 2013;Pärnoja et al., 2014) and in extreme cases, e.g., within giant kelp forests, only a fraction of light above algal canopy reaches the basal subcanopy species (Reed and Foster, 1984). This has large but mostly unknown affects on photosynthesis and growth of macroalgal assemblages. Even though algal photosynthesis is saturated in the upper layers of algal canopies with increasing irradiance, community photosynthesis will remain unsaturated because photosynthetic tissue in the lower layers of algal canopy has an unused photosynthetic potential (Binzer and Middelboe, 2005;Tait and Schiel, 2010).
Nutrient availability strongly affects the production of macroalgae but the responses are mostly taxon-specific and highly influenced by the local light climate (Raven and Hurd, 2012). Importantly, macroalgae can store nutrients in their tissue in order to overcome temporal disparity between nutrient availability and optimal growth season. For example, the marine Laminariales grow fastest during summer when irradiance is the highest, but when nutrients are depleted above the thermocline. During the peak growth the algae use their nitrogen and phosphorus reserves that were stored during winter when the nutrient availability is greater but light levels are lower (Chapman and Craigie, 1977;Lüning, 1990). In addition, the Charales can effectively transport nutrients from rhizoids (in sediment with high N and P but no light) to shoots (in water with low N and P but high light irradiance) (Vermeer et al., 2003). However, when nutrients are plentiful in the water column and light levels are sufficient to support photosynthesis, direct nutrient uptake from the water through the outer cells of the branchlets and stems is expected (Raven, 2003).
Temperature controls primary productivity of macroalgae to a lesser extent than light and nutrients unless the temperature is above or below the tolerances of the species. Within tolerance limits, macroalgae exhibit relative uniform responses to temperature change (Raven and Smith, 1978;Wiencke and tom Dieck, 1990). The species-specific temperature demands for survival and production relate to different types of temperature stress the macroalgal species have locally adapted (Breeman and Pakker, 1994). Controlled experimental trials further demonstrate the high plasticity of response of macroalgae to specific limiting temperature conditions in which a long-term acclimation of macroalgae at extreme temperatures significantly relaxes this temperature limitation (Nejrup et al., 2013). Thus, estimations of primary production in natural assemblages are increasingly relevant given the growing need for obtaining realistic response patterns to a complex of environmental factors.
Charophytes are found in diverse aquatic habitats, so this group consists of both early colonizers of newly-formed water bodies as well as species preferring more stable environmental conditions and habitats (Blindow and van de Weyer, 2016). The lush charophyte meadows play a underpinning role by providing food and/or refuge to zooplankton (Puche et al., 2020), benthic invertebrates (Kotta et al., 2004), fish and fish larvae (Hargeby et al., 2005), microbes (Kataržytė et al., 2017), and waterfowl (Schmieder et al., 2006). They also reduce water turbidity by stabilizing the sediment with rhizoids (Nõges et al., 2003), immobilizing nutrients (Kufel and Kufel, 2002), and inhibiting phytoplankton biomass (Hilt, 2015). Charophytes are most common in freshwater habitats but are also found in brackish water. The Baltic Sea is one of the best-studied brackish regions with respect to charophytes. In this marginal environment, the shallow low salinity embayments are often dominated by charophytes . The most common and widespread species in such habitats is Chara aspera Willd., a species characterized by a broad environmental niche space (Herkül et al., 2018).
To address the importance of scale on the factors regulating macroalgal production, we quantified experimentally the seasonal variability in individual photosynthesis, community production and biomass growth of the charophyte C. aspera in the brackish Baltic Sea. We then analyzed how seasonal change in the assemblage of environmental variables manifested in the production of macroalgae at these different scales.
We hypothesized that C. aspera with its well-known broad environmental tolerances can maintain high photosynthetic and growth potential across wide gradients in light, temperature, and water nutrients. We also expected that the relationships between the environment and photosynthesis when assessed on individual plants are strong and simple. At the community level, however, charophytes have an ability to adjust photosynthesis and growth to a different set of environmental conditions. Moreover, below-ground processes may occasionally contribute to the oxygen fluxes. Therefore relationships between light, nutrient concentration, temperature, community photosynthesis and growth are weaker and more complex compared to individual plants.

Description of Study Area
The study was carried out in Rame Bay from March 2014 to March 2015. Rame Bay is a shallow and semi-enclosed bay in the northeastern Baltic Sea. Depending on the prevailing wind direction, the bay receives eutrophicated waters from the Gulf of Riga or less eutrophied waters from the northern Baltic Proper (Suursaar et al., 1998). The maximum depth of the bay is 1.5 m. The bottom is composed of fine sand and a thick layer of fine silt (Torn and Martin, 2003). Salinity varies between 2 and 5.6 psu with lower values during heavy rainfall. Rapidly increasing water temperatures during spring, high temperatures in summer, and a long period of subzero temperatures in winter are typical. The ice-cover of shallow bays in the region lasts 4-5 months on average (Jevrejeva and Leppäranta, 2002).
Most of Rame Bay is covered by a rich charophyte community. Five species of charophytes are found in the bay with C. aspera being the most widespread species. The seasonal dynamic of charophytes follows general growth patterns of the northern part of the temperate zone, which is characterized by the highest growth rates at the beginning of summer and degradation of the community in autumn (Torn et al., 2006(Torn et al., , 2010. C. aspera hibernate both as detached mats and bulbils in the study area (current study, Torn et al., 2010).

Sampling of Environmental Variables and Communities of C. aspera
The study was carried out at depths of 0.9-1.3 m. The study design combined a traditional sampling of water quality variables and benthic communities with the results from the latest oceanographic technology to achieve an accurate environmental assessment of the Rame Bay area. In addition to these observations, the production of C. aspera was measured experimentally at three organizational scales (for more details see below). C. aspera can be visually determined from other charophytes in the study area as it has clearly distinguishable globular white bulbils on the rhizoids. This feature enabled to run the production of C. aspera without a need for prior determination in the laboratory.
The oceanographic instruments were deployed 20 cm above the seabed to record short-term variability of temperature, light availability, and water velocity in the study site from March 2014 to March 2015. These results provided information on the environment around the experimental site during and in between samplings. Temperature was measured with a General Oceanics thermologger. Light availability was recorded by an Odyssey photosynthetically active radiation (PAR) logger. Temperature and PAR was logged at 60 min intervals. Water velocity was recorded by a JFE Advantech Logger version of the INFINITY 2-D electro-magnetic current meter AEMD-USB. During deployment the instrument recorded current speed once per minute. Water level was measured using a JFE Advantech Logger version of the INFINITY-WH AWH-USB wave height meter. During deployment the instrument recorded pressure once per second. Nutrient concentrations in water were measured triplicate once a month. Immediately after collection, samples were preserved and brought to the laboratory in an insulated ice box. We used standard ISO procedures for measuring total N and total P (ISO 11905 and ISO 15681). We did not use neither filtering nor fractioning for our samples. The analyses were made using a Skalar San + + autoanalyzer.
The charophyte community was sampled for coverage and biomass once in a month. Total coverage of all macrophytes and percentage cover of all visually distinguishable species were recorded by SCUBA diver during each sampling. Biomass samples were collected in three replicates by SCUBA diving using a 20 × 20 cm frame. Samples were stored deep frozen (-18 • C) until analysis. In the laboratory all species were identified to the species level. The dry weight of species was obtained after drying the individuals at 60 • C for 2 weeks (i.e., the time needed to achieve the same weight for two consecutive measurements). Growth in biomass was measured as the difference in total biomass of the plant community between two consecutive sampling months. The chlorophyll-a concentration was measured spectrophotometrically by Libra S32 after extraction of chlorophyll from algal tissue in ethanol. To determine the phosphorus content of the algae, the material was mineralized by heating with HNO 3 and HClO 4 (Miller, 1998).
The content of the elements in the resulting solution was determined by MP-AES spectrometry. The Dumas method was used to measure the nitrogen and carbon content of the algae (Buckee, 1994).

Photosynthetic Production of Individual Plants
In each study visit the net photosynthetic rate of C. aspera was measured in triplicate from 600 ml glass bottles filled with seawater from the study site. Bottles with about 0.5 g dry weight of C. aspera and without algae (controls) were incubated in situ on special frames at a depth of 1 meter around midday. The hourly net production rates (mg O 2 g dry weight −1 h −1 ) were calculated from the differences in oxygen concentrations, measured over the incubation period (1-3 h). Oxygen concentration and water temperature in incubation bottles was measured with OptiOx optical dissolved oxygen sensor connected to SevenGo DO meter (Mettler Toledo). The dry weight of C. aspera was determined after drying at 60 • C for 48 h. During the experiment the photosynthetically active radiation (PAR; µmol m −2 s −1 ) at the incubation depth was measured using Li 188B (Li-Cor Inc.) quantum meter.

Production in Algal Communities
The production of C. aspera communities was measured in triplicate from transparent plexiglas chambers holding 29 l of water with a surface area 850 cm 2 . The chambers were placed randomly on the seafloor over the homogeneous monospecific stand of C. aspera. A rubber seal assured an airtight seal between the chamber and the surrounding environment.
Oxygen concentration in the chamber was measured every minute using a calibrated Optode-type oxygen sensor (Aanderaa Instruments) connected to a data logger (data recorder by JFE Advantech Co., Ltd.). This instrument also provides data on water temperature. Changes in dissolved oxygen averaged over minute intervals were used as a proxy of community (algal assemblage + sediment) net production. During deployment irradiance above the charophytes was measured every minute using a calibrated spherical quantum sensor connected to a data logger (ultra-miniature logger for light intensity by JFE Advantech Co., Ltd.). Each production experiment lasted 3 h on average.
After deployment, macroalgae within incubation chambers were harvested and stored in a deep freezer at −20 • C. The subsequent sorting and determination of species were performed in the laboratory using a stereomicroscope. The dry weight of each macroalgal community was obtained after drying the individuals at 60 • C for 2 weeks. Combining this information with the oxygen flux measurements described above enabled us to express net primary production in mg O 2 g dry weight per macrophyte min −1 .

Statistical Analyses
The contribution of different environmental variables on the individual, community photosynthesis and growth in biomass of C. aspera was explored using the Boosted Regression Trees technique (BRT). BRT models are capable of handling different types of predictor variables and their predictive performance is superior to most traditional modeling methods [see e.g., comparisons with GLM, GAM and multivariate adaptive regression splines, Elith et al. (2006) and Leathwick et al. (2006)]. Overfitting is often regarded as a problem in statistical modeling, but can be overcome by using independent data sets. The BRT modeling iteratively develops a large ensemble of small regression trees constructed from random subsets of the data. Each successive tree predicts the residuals from the previous tree to gradually boost the predictive performance of the overall model (Elith et al., 2008). Important parameters in building BRT models are the learning rate and tree complexity. The learning rate determines the contribution of each tree to the growing model and tree complexity defines the depth of interactions allowed in a model. A tree complexity of 1 assesses only main effects; A tree complexity >1 includes interactions. Different combinations of these parameters may yield variable predictive performance but generally a lower learning rate and inclusion of interactions gives better results (Elith et al., 2008). In the current study, the model learning rate was kept at 0.001 and tree complexity at 5. In order to avoid potential problems of overfitting, unimportant variables were removed using a simplify tool. In order to eliminate non-informative variables, the simplify tool progressively simplifies the model, then re-fits the model and sequentially repeats the process until a stopping criterion is reached. Such simplification is most useful for small data sets where redundant predictors may degrade performance by increasing variance. Model performance was evaluated using the cross-validation statistics calculated during model fitting (Hastie et al., 2009). The BRT modeling was done in R using the gbm package (Elith et al., 2008). Standard errors for the predictions and pointwise standard errors for the partial dependence curves, produced by R package "pdp" (Greenwell, 2017), were estimated using bootstrap (100 replications). Multicollinearity can be an issue with BRT modeling when assessing if and when environmental variables are of ecological interest. Thus, prior to modeling, the Pearson correlation analysis between all environmental variables was calculated in order to avoid including highly correlated variables into the model. The correlation analysis showed that most variables were only weakly intercorrelated (r < 0.5).

Seasonal Changes in Abiotic Environmental Conditions
Water temperature and light availability above the seabed followed the same seasonal pattern with higher values measured in July and lower values in winter (Figure 1).   and December. August stood out among other summer months with the highest nutrient concentration and variability of water movement.

Seasonal Changes in the Community of C. aspera
Charophytes had high biomass in summer and autumn, whereas the algal cover was stably high throughout the year including winter (Figure 2). C. aspera overwintered as a sparse and flat biomass residue. Small fresh shoots were observed only in April and the full community established by July. The rapid decay in biomass occurred during the formation of ice cover. Nutrient content (nitrogen and phosphorus) of C. aspera was high in decomposing algal material in winter and low in summer. The chlorophyll a content of C. aspera plants varied between months with no clear seasonal patterns (somewhat lower values in winter and higher values in summer).

Seasonal Changes in the Individual and Community Production of C. aspera
Net production rates of C. aspera measured at an individual and community level followed similar seasonal patterns over the whole study period with maximum values observed in June and July and lowest values during winter months. The only exceptions were abnormally low community production values in May and August related to the presence of high amounts of organic debris on the sediment. Similarly, the biomass growth of C. aspera was high from June to September and low in other months (Figure 3).

Relationships Between Environment, Individual, Community Production and Biomass Growth of C. aspera
Boosted regression trees models on the production of C. aspera accounted for a significant proportion of the variability at the studied scales with r 2 values as follows: individual photosynthesis 0.855, community production 0.854 and change in algal biomass 0.573 (Table 1). Temperature was the best overall predictor in the models of individual photosynthesis and change in algal biomass and the second-best predictor in the model of community production. However, the thresholds of algal production to increasing temperature varied among the studied scales. Light availability was another generic predictor of the production of FIGURE 3 | Seasonal dynamics in individual and community net production rates and biomass growth of C. aspera in the study area (means ± Standard Error). Extreme events (negative values) of community net production associated with the presence of high amounts of organic debris on the sediment are shown in red. A moving average smoother was used to visualize the trend in biomass growth.
C. aspera in all models. Like temperature, the complexity of (functional form) relationships between light and production did not change among models and the differences were due mostly to different light thresholds in different models. Among the studied environmental variables, nutrient concentration in water and chloroplast content in algal tissue were not significant and the variables were not retained in the models.
Photosynthesis in individual plants was related to temperature and light intensity with these variables accounting for nearly 100% of the predicted variability in the production of C. aspera at this scale. Increasing temperature increased the production of individual plants up to a threshold of 20 • C. Above this tipping point, variability in plant production was independent of the ambient water temperature. The production of individual plants had two clear phases: low production at light levels <500 µmol quanta m −2 s −1 and high production at light levels >600 µmol quanta m −2 s −1 with light saturation remaining at 550-600 µmol quanta m −2 s −1 (Figure 4).
Like photosynthesis in individual plants, the production of benthic communities correlated positively with light intensity. However, the relationship was characterized by a smoother response and a higher light saturation point at 800 µmol quanta m −2 s −1 . Importantly, community production decreased when the coverage of C. aspera attained >95%, indicating greater light limitation in denser algal communities. The relationship between temperature and community production exhibited two phases: low production at temperatures <13 • C and high production at temperatures >14 • C. BRT analysis excluded abnormally low community production values in May and August related to the presence of high amounts of organic debris on the sediment and high levels of belowground respiration. This suggests that when temperature exceeds 21 • C, the community net production is expected to decline sharply ( Figure 5).
Temperature, light climate and underwater currents were the key predictors of change in algal biomass accounting for >70% of the predicted variability in the biomass change. Importantly, high levels of growth were sustained at much lower light intensities, already above 100 µmol quanta m −2 s −1 compared to thresholds of 600 and 800 µmol quanta m −2 s −1 in the models of individual and community production. Similarly, high growth values were predicted only at velocities <1 cm s −1 . Moreover, increasing Coverage of Chara aspera 13 10 Water velocity 23 Water level 17 The contributions of environmental variables are relative contributions to the models' total variable explained (%). water level was associated with reduced biomass of macroalgal communities (Figure 6).

DISCUSSION
In this study we experimentally quantified seasonal variability in individual photosynthesis, community production and growth in biomass of the charophyte C. aspera in the Baltic Sea along with the key environmental variables that are known to determine the production of macroalgae. We found that temperature and light availability were the best overall predictors of the production of C. aspera, followed by macroalgal coverage, flow velocity and water level. Importantly, the observed scale notably regulated the importance of different predictors, as well as the relationships between the environment and plant responses in the models. This implies that although the physiological responses of C. aspera should always be interpreted in the context of scale, temperature and light availability predict a significant proportion in the variability of charophyte production. Low light can limit the production potential of temperate macroalgae, particularly during low-light seasons or in turbid waters (Wallentinus, 1978;Littler, 1980); many studies have suggested a direct dependence of charophytes on light availability (Libbert and Walter, 1985;Blindow, 2000;Schneider et al., 2006;Schubert et al., 2018). However, due to high tolerances of C. aspera to low light environments (Blindow, 1992;Sand-Jensen et al., 2000), we hypothesized that the species could maintain positive photosynthesis and growth potential in a wide range of light conditions. In general, this held true for individual photosynthesis, community production and growth in biomass. However, individual photosynthesis was typically higher than community production under the same light conditions. This difference in photosynthesis can be partly attributed to lower light availability (light limitation) in a three-dimensional algal canopy than to individual algae exposed to ambient light levels. For the same reason, the photosynthesis-irradiance curve at community level was smoother than that of individual plants. Specifically, different parts of the algal canopy are expected to receive different amounts of light and therefore are expected to respond differentially. For example, the upper parts may often be oversaturated with light, whereas lower parts are undersaturated. This results in different light adaptation (different light saturation points) along the spatial light gradient and a smooth response of the algal community to light. This also applies for C. aspera, which is known to increase the light saturation point of photosynthesis with increasing irradiance from 100 to 400 µmol quanta m −2 s −1 . Similar responses of net photosynthesis along light gradients have been observed for C. aspera in other areas (Kovtun-Kante et al., 2014) and for other charophyte species, such as C. baltica and C. canescens (Küster et al., 2004). Beside shelfshading, sediment respiration is an additional explanation for the differences between individual and community photosynthesis. At high load of organic debris and/or elevated temperatures algal respiration and sediment decomposition may dominate over benthic production (Carvalho et al., 2005).
Differential responses of photosynthesis and growth to light irradiance is also related to the different approaches used to calculate production at different scales. Specifically, the photosynthesis of individual plants and communities is expressed as relative change per plant biomass, whereas biomass growth is reported on an absolute scale. In spring and early summer with abundant light, photosynthesis of individual plants and communities is enhanced. However, the absolute change in biomass is actually low because the plants themselves are small. Later, when the plants are taller, biomass increases are greater despite moderate light conditions and low photosynthesis of individual plants and communities. This is happening because some inorganic carbon is fixed into autotrophic biomass even at low photosynthesis levels between 100 and 500 µmol quanta m −2 s −1 .
In addition to light irradiance reaching the algal canopy, water turbidity affects light availability and the production of charophytes in many habitats (Hellström, 1991;Schneider et al., 2006). Likewise, the sediments in Rame Bay are composed of fine fractions (very fine sand, silt, and clay) and even slow underwater currents (1 cm s −1 ) can bring silt and mud particles into the water column and thereby modulate the light environment as well as deposit sediment particles on plant thalli. Consequently, less light reaches plant photopigments and the rate of community photosynthesis and biomass growth decreases. This process is supported by experimental evidence in which long exposure to high levels of suspended sediment leads to cessation of growth and ultimately death of charophytes (Henricson et al., 2006). On the other hand, an experimental reduction of light availability to 25% of the natural irradiance does not induce a chronic compensatory reduction of the photosynthesis of charophyte communities (Kovtun-Kante et al., 2014). This suggests that charophytes recover easier from longer-term light reduction than sedimentation and may explain why water velocity (i.e., the intensity of resuspension) was the key predictor in the models of biomass growth. In addition to underwater currents, water level explained a significant part of charophyte growth in biomass. The current charophyte meadows are located in very shallow water (1 m) and a decline in water level of only a few tens of centimeters is expected to improve light conditions in the algal canopy.
It has been suggested that an increased production of photosynthetic pigments and chloroplasts supports more efficient production of charophytes and helps plants to recover when light becomes limiting (Andrews et al., 1984;Küster et al., 2000;Marquardt and Schubert, 2009). Our study did not support this, as variability in photosynthesis and growth along light gradient was independent of the chloroplast content of charophyte tissue. The chloroplast content was measured during fieldwork but due to lack of significance, the variable was not retained in the models and therefore not reported. It has been also suggested that charophytes elongate quickly when light levels are low. This adaptation enables charophytes to reach the water surface where light is more plentiful (Andrews et al., 1984;Henricson et al., 2006). However, the shoot elongation in C. aspera is salinity-dependent and brackish water plants fail to elongate under impoverished light conditions . Our study corroborated this as underwater light conditions were not inversely related to plant production throughout the light gradient.
We also hypothesized that C. aspera could maintain positive production and growth potential in a wide range of temperature conditions. This hypothesis held true for individual and community photosynthesis but not for growth in biomass. Interestingly, photosynthesis in individual plants increased with increasing temperature throughout the gradient, whereas this temperature dependence translated into relatively stable positive community photosynthesis with somewhat higher photosynthesis at temperatures >14 • C. Nevertheless, at temperatures >22 • C benthic production declined substantially due to the dominance of below-ground processes (e.g., sediment decomposition) and algal respiration over algal production (Carvalho et al., 2005;Torn et al., 2006). Thus, C. aspera should sustain positive growth in its Baltic Sea habitats even under relatively extreme temperature conditions, which substantiates earlier studies showing that temperature does not significantly control photosynthesis and growth of charophytes (Raven et al., 1979;Libbert and Walter, 1985;Torn et al., 2006). Due to such broad temperature tolerance, charophytes are expected to benefit from contemporary climate change and broaden their distribution range in the Baltic Sea region (Torn et al., 2020).
Finally, we expected that C. aspera could maintain high photosynthesis and growth potential across wide gradients in water nutrients. This hypothesis was fully supported as nutrient concentrations in water were not significant in all models. Although charophytes are able to acquire nutrients from the water column (Raven, 2003), nutrient concentration in water has been shown to play a marginal effect on the physiological state of charophytes (Blindow, 1988;Munsterhjelm, 2005). Charophytes in our study area seem to rely on nutrients stored in their tissues as shown by the inverse relationship between macroalgal nutrient content and the photosynthesis of charophyte individuals (R 2 = 0.15). Such strategy is widely used among long-living macroalgae in order to overcome temporal and/or spatial mismatch between nutrient availability and optimal growth season (Chapman and Craigie, 1977;Lüning, 1990;Raven, 2003).
Prior to the study, we anticipated that models of photosynthesis in individual plants would involve fewer important environmental predictors than the models of photosynthesis and growth of charophyte communities. This prediction also held true. Although water temperature and light irradiance account for most of variability in the photosynthesis of individual plants, this share is less than in models of community photosynthesis and growth. The existence of a weak negative relationship between algal cover and community photosynthesis at the community level suggests that density dependence does play some role in photosynthesis regulation (Binzer and Middelboe, 2005;Tait and Schiel, 2010). Moreover, water properties (water level and flow velocity) were important in regulating charophyte growth. However, as judged by the importance of different predictors in the models, the effects of macroalgal density and water properties were much weaker than light dependence of production and growth of C. aspera.
The availability of different forms of inorganic carbon may limit the photosynthesis of submerged vegetation including C. aspera (van den Berg et al., 2002). Shallow water habitats of our study area are also characterized by a large amplitude of natural variability in pH and pCO2 with a daily range of pH recorded between 8 and 9 (Pajusalu et al., 2013). Among different charophyte communities inhabiting our study area C. aspera is way less sensitive to inorganic carbon limitation compared to C. horrida and C. tomentosa (Pajusalu et al., 2015). Thus, as compared to light and temperature, the limitation of inorganic carbon is likely not a prominent stress to the studied community.
Our study did not focus on interspecific interactions that potentially shape the photosynthesis and growth of C. aspera. A wealth of experimental studies has demonstrated that marine herbivores play key roles in regulating macroalgal communities (Hayward, 1988;Sala and Graham, 2002;Kotta and Witman, 2009). However, previous experiments suggest that herbivory or plant fragmentation induced by herbivores is not crucial in the studied charophyte communities (Kotta et al., 2004). Moreover, C. aspera largely dominates the studied charophyte community and interspecific competition among algae and submerged vegetation is slight. It is plausible that spatial heterogeneity or differential utilization of abiotic resources are behind the pattern in distribution of species (Chambers and Prepas, 1990).

CONCLUSION
A large set of experimental data is a prerequisite to understand and accurately model the patterns of production of any macroalgal assemblages. In this study we experimentally quantified seasonal variability in individual photosynthesis, community production and growth in biomass of the charophyte C. aspera in the Baltic Sea along with the key environmental variables known to determine the production of macroalgae. The experiments suggested that the production and growth of C. aspera was defined largely by light and partly by temperature at all studied scales. However, algal production was systematically lower at the community than at an individual level under the same light conditions, suggesting the importance of shelfshading and below-ground processes in natural plant assemblages (e.g., Carvalho et al., 2005;Binzer et al., 2006). Moreover, the observation scale also defined response types between the environment, plant production and growth. This implies that the patterns of variability of the production and growth of macrophytes should always be interpreted in the context of scale and any multiscale model development should involve experimental validation at all important scales.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.