Changes in body size spectra of benthic caridean shrimps ( Decapoda : Caridea ) and snails ( Gastropoda ) as response to seasonal variability

Ecologists have been largely interested in the description and understanding of the power scaling relationships between body size and abundance of organisms. Many studies have focused on estimating the exponents of these functions across taxonomic groups and spatial scales, to draw inferences about the processes underlying this pattern. The exponents of these functions usually approximate -3/4 at geographical scales, but they deviate from this value when smaller spatial extensions are considered. This has led to propose that body size-abundance relationships at small spatial scales may reflect the impact of environmental changes. This study tests this hypothesis by examining body size spectra of benthic shrimps (Decapoda: Caridea) and snails (Gastropoda) in the Tamiahua lagoon, a brackish body water located in the Eastern coast of Mexico. We measured water quality parameters (dissolved oxygen, salinity, pH, water temperature, sediment organic matter and chemical oxygen demand) and sampled benthic macrofauna during three different climatic conditions of the year (cold, dry and rainy season). Given the small size of most individuals in the benthic macrofaunal samples, we used body volume, instead of weight, to estimate their body size. Body size-abundance relationships of both taxonomic groups were described by tabulating data from each season into base-2 logarithmic body size bins. In both taxonomic groups, observed frequencies per body size class in each season were standardized to yield densities (i.e., individuals/m3). Nonlinear regression analyses were separately performed for each taxonomic group at each season to assess whether body size spectra followed power scaling functions. Additionally, for each taxonomic group, multiple regression analyses were used to determine whether these relationships varied among seasons. Our results indicated that, while body size-abundance relationships in both taxonomic groups followed power functions, the parameters defining the shape of these relationships varied among seasons. These variations in the parameters of the body size-abundance relationships seems to be related to changes in the abundance of individuals within the different body size classes, which seems to follow the seasonal changes that occur in the environmental conditions of the lagoon. Thus, we propose that these body size-abundance relationships are influenced by the frequency and intensity of environmental changes affecting this ecosystem. Rev. Biol. Trop. 64 (1): 33-44. Epub 2016 March 01.

The most widely studied of these allometric relationships is that between body size and the abundance of individuals (D), often measured as density (Basset, Sangiorgio, & Pinna, 2004;Marquet, Navarrete, & Castilla, 1990;Nee, Read, & Harvey, 1991;White et al., 2007).For most taxonomic groups of animals, the scaling exponents of these relationships usually approximates to -3/4 and, hence, D = a BS -3/4 (Brown et al., 2000;Cyr, 2000;Peters, 1986).This has led to propose the energetic equivalence rule, which states that the total use of energy by populations is independent of the body size of their individuals (Nee et al., 1991).This energetic equivalence rule assumes that species abundances in natural communities are constrained by the total amount of energy, and resources that are available to individual organisms, and scaling exponents close to -3/4, are expected when the energetic balance of the ecosystems is at equilibrium (Marquet et al., 1990(Marquet et al., , 2005)).On the other hand, the normalization constant summarizes the effects of resource supply rates, space availability, and other sources of variability affecting body size and species densities within ecosystems (Brown et al., 2004;Marquet et al., 2005).
The scaling relationship described above performs well at large geographical scales, such as continents, but its shape is still being matter of discussion on smaller spatial scales (White et al., 2007).Aquatic ecologists, for example, have usually analyzed the frequency distribution of body sizes regardless the species to which organisms belong, or body size spectra, and studies conducted in small water bodies have reported scaling exponents that largely deviate from the expected value of -3/4 (Boix, Sala, Quintana, & Moreno-Amich, 2004;Schmid, Tokeshi, & Schmid-Araya, 2000, 2002).Indeed, recent studies on benthic macroinvertebrates from Mediterranean lagoons have shown that local body size-abundance relationships can be polygonal in shape, with the highest population density per size class increasing with average individual body size of 0.25 mg and decreasing beyond this threshold (Barbone, Rosati, Pinna, & Basset, 2007;Basset et al., 2008).These differences between observed and expected body size-abundance relationships have been attributed to external disturbances that affect aquatic ecosystems (e.g., seasonal changes in water discharge regimens, anthropogenic disturbances or variations in resource supply rates from terrestrial ecosystems) instead of energetic or resources constraints (Barbone et al., 2007;Basset et al., 2008;Schmid et al., 2000).Therefore, the scaling relationships expected on geographical scales may not be attained at small spatial scales when local species assemblages are affected by external disturbances or cyclical environmental changes.Indeed, at small spatial scales, scaling exponents of body size spectra may vary idiosyncratically across time because environmental variability may prevent the community from reaching energetic equilibrium states.
To test these hypotheses we focused in a coastal brackish lagoon of Mexico that is recurrently affected by seasonal changes in water quality conditions.In this ecosystem, we assessed the body size spectra of two groups of benthic macroinvertebrates including infaunal and epifaunal caridean shrimps (Crustacea) and snails (Gastropoda).Body size spectra of these two groups were expected to fit power law functions, but we also expected changes in the parameters describing these relationships (a and b) after changes in water quality conditions.

Study site:
This study was conducted in the Tamiahua lagoon located in the North of the State of Veracruz, Mexico (21º15' N -97º19' W).It covers 88 000 ha and has an average depth of 2 m (Contreras, 1988).Four small rivers with irregular flow (Tampache, Carbajal, Cucharas and Tancochín) discharge their waters in this lagoon and a sand barrier separates it from the Gulf of Mexico (Cruz, 1968).The lagoon is connected to the Gulf of Mexico through a 15 m wide natural channel located in its Southern section, but its Northern section has been artificially connected to the sea through a 20 m wide channel.Local climate is tropical-type, but there are three clearly distinct seasons throughout the year.The coldest period occurs between October and February, when air temperature may fall below 15 °C.This cold season is characterized by frequent invasions of polar air masses from the North, which cause violent hurricane systems and some precipitation events (less than 400 mm across this season) (Contreras, 1988).It is followed by the dry season, between March and May, with air temperatures between 22 °C and 30 ºC and scarce precipitation events (less than 100 mm across this season) (Contreras, 1988).Finally, the rainy season occurs between June and September, concurring with the warmest period of the year.Air temperatures in this season surpass 40 °C and it concentrates up to 85 % of annual rainfalls (above 1 000 mm across the season) (Contreras, 1988).During the rainy season, there are strong increases in the flows of neighboring rivers, which increase the load of sediments and organic matter in the lagoon (Cruz, 1968).
Environmental variables: Water quality conditions in the lagoon were assessed in eight fixed sampling stations located 3-5 km away each other on the Western coast of the lagoon (Table 1).During 2008, dissolved oxygen, salinity, pH, and water temperature were measured at each sampling station at the middle of the three seasons (January 11 th , April 18 th and July 25 th ).For this we used a multiparameter probe (6600V2, YSI Incorporated, USA) that was immersed 20 cm deep in the water column.A one-liter core sample was also taken from the bottom of the lagoon at each sampling station.These samples were placed at 5 °C immediately after collection to reduce metabolic activity and taken to the Laboratory of Physicochemical Processes of University of the Americas Puebla (Mexico) to determine the amount of organic matter in sediments and the chemical oxygen demand (COD).These analyses were conducted within three days after collecting the sediment samples.Repeated measures ANOVA were used to compare these variables among seasons.In these analyses, sampling stations were considered as replicates and seasons were the repeated measures (Zar, 2010).When significant differences were detected, the Tukey's post hoc test was used to assess differences between seasons.

Macroinvertebrate sampling and body size measures:
Benthic macroinvertebrates were sampled in the stations mentioned above (Table 1) on the same dates in which the environmental variables were measured (January 11 th , April 18 th and July 25 th ).At each station we drew a 30 m line-transect located 20 m away from the coast line.On this transect, we collected one-liter sediment samples every 10m with a cylindrical hand corer (20 cm diameter x 32 cm depth).This resulted in four samples per transect, which were pooled to obtain a single composite sediment sample at each sampling station.Composite samples were sieved through a 0.5 mm mesh to remove sediments and macroinvertebrates retained in the sieve were preserved in formalin solution (7 %) buffered with sodium borate (28 g/L) (Marini, Pinna, Basset & Mancinelli, 2013;Sangiorgio et. al, 2014).Macroinvertebrates from each sample were separately stored and taken to the Laboratory of Entomology of the University of the Americas Puebla (Mexico), where they were washed with distilled water and preserved in ethanol solution (70 %).
Because caridean shrimps and snails were the most abundant organisms in the samples, these two taxonomic groups were selected to perform the body size spectra analyses.These analyses considered all species of caridean shrimps and snails retained in the samples after the sieving (i.e., all individuals larger than 0.5 mm).We were unable to identify collected specimens because detailed taxonomic descriptions of local fauna are still lacking.Nevertheless, previous studies reported that caridean shrimps of this lagoon are subsurface deposit feeders (algae grazers and scavengers) that belong to the families Palaemonidae and Hippolytidae (Álvarez, Villalobos, Rojas, & Robles, 1999).Snails of the lagoon belong to the families Acteonidae, Acteocinidae, Cerithiidae, Nassariidae, Turriidae and Vitrinellidae, which are subsurface carnivores and occasional scavengers (Reguero & García-Cuba, 1989;Reguero, García-Cuba, & Zúñiga, 1991).
Collected specimens were extremely small to be weighed and, for this reason, we used body volume instead of body mass as a surrogate of body size (Calder, 2000;McClain, 2004;Siemann, Tilman, & Haarstad, 1996).For this, each specimen was placed under a stereoscopic dissection microscope with scaled lens (STEMI SV6, Zeiss, Germany).In the case of shrimps we measured their length (l), from the rostrum to the last segment of the body, and width (w); antenna, flagellum, pereopods, pleopod and telson were not included in these measures because these structures are highly variable among species.Thus, the body volume (V shrimp ) of shrimps was estimated by following the mathematical formula of a cylinder: V shrimp = l w 2 .Because the body of snails is almost conical, to estimate their volume we measured the length of each shell (h) and the width of its base (b), and the body volume (V snail ) was estimated as: V snail = [p ( 1 / 2 b) 2 h] / 3.In both cases, decaying organisms and empty snail shells were not measured to ensure that only individuals that were alive at the time of sampling were included in the analyses described below.
Body size spectra: Body size spectra describe the frequency distribution of individuals according to body size classes, irrespective of the species they belong to (Sheldon, Prakash, & Sutcliffe, 1972;White et al., 2007).This method is commonly employed to assess body size abundance relationships when sampled individuals are too numerous and the species are hard to recognize, as occur in benthic macroinvertebrate assemblages (Dinmore & Jennings, 2004;Drgas, Radziejewska, & Warzocha 1998;Kaariainen & Bett, 2006;Mancinelli, Pinna, & Basset, 2008;Mouillot et al., 2006).
The main drawback of this method is defining objective criteria for establishing the width of the body size classes to which individuals are assigned.In most studies, the researchers arbitrarily define these body size classes, but this procedure can bias the shape of these relationships according to the width and number of body size classes that are established (Mancinelli et al., 2008).To avoid such arbitrariness we used the logarithmic binning method suggested by White, Enquist, & Green (2008).This method consists in constructing as many body size classes as necessary to cover the entire range of body sizes, but the width of the size classes increases by following a logarithmic scale.Once all these body size classes are defined, the individuals are assigned to the size class they belong to and counted.In our case, the width of the size classes was established by using a log 2 -scale, as recommended by White et al. (2008) for assessing body size spectra in small organisms.
Given that the number of individuals within each body size class was obtained from 32 L of sediment at each season (four 1 L samples x eight sampling stations), the density of individuals at each size class was standardized to a volume of 1 m 3 for both taxonomic groups (i.e., individuals/m 3 ).After that, nonlinear regression analyses were separately performed for each taxonomic group at each season to assess whether body size spectra fitted to power functions and hence estimate the normalization constants (a) and scaling exponents (b) of these relationships.For these analyses, we used nonlinear regression models with the following function: , where A is density of individuals and V is the median value of the different body volume classes (Peters, 1986).
The analyses described above allowed us to determine the values of a and b for the body size spectra of shrimps and snails at each season, but using nonlinear methods to statistically compare these values among seasons is extremely complex.For this reason, these relationships were linearized by applying the log 2transformation to the median value of body size classes and their respective abundances of individuals.In this way, the values of a and b were compared among seasons by using multiple linear regression analysis with categorical variables (Kutner, Nachtsheim, Christopher, Neter, & Li, 2005).Therefore, the log-transformed density of individuals within each size class was the continuous response variable, the log-transformed median value of the different size classes was the continuous predictive variable, and the sampling seasons (cold, dry and rainy season) were the three levels of the categorical predictive variable.Regression models used in these analyses also included all possible interaction terms between predictive variables, allowing us to assess differences among estimated parameters (Kutner et al., 2005).Differences in intercepts (log 2 a) and slopes (b) of regression functions obtained for the different seasons were assessed with t-tests (Kutner et al., 2005).

Physical-chemical variables:
All physical-chemical variables measured in the water column and sediments differed among seasons (repeated measures ANOVA p < 0.01 in all cases).Post hoc comparisons indicated that dissolved oxygen, salinity and COD significantly differed among all three seasons.The highest average value of dissolved oxygen was found in the dry season, while its lowest value was recorded in the rainy season (Fig. 1A).Both salinity (Fig. 1B) and COD (Fig. 1C) showed their higher values in the dry season, while their lower values were observed in the cold season.Water temperature (Fig. 1D) showed similar values in the dry and rainy seasons but, in both cases, these values were higher than in the cold season.Organic matter (Fig. 1E) and pH (Fig. 1F) were higher in the rainy season than in the other two seasons.
Body size spectra: Densities of shrimps and snails significantly decreased with increasing body size in the three seasons.Simple nonlinear regression analyses (Table 2) indicated that these body size spectra are well described by power functions (Fig. 2).Nevertheless, normalization constants and scaling exponents of these relationships changed across seasons in both taxonomic groups.For caridean shrimps, the power function with the highest normalization constant (a) and the lowest scaling exponent (b) was obtained in the dry season (Fig. 2).These values of a and b coincided with the period of the year in which higher abundances of individuals smaller than 1mm 3 in body size were recorded (up to 40 875 ind/m 3 ).The rainy and cold seasons showed lower densities of shrimps smaller than 1 mm 3 (below 20 000 ind/ m 3 ), which lead to lower values of a (Fig. 2).Conversely, the values of b in these two seasons were higher than those found in the dry season (Fig. 2).For snails, the highest value of a and the lowest value of b were obtained in the rainy season (Fig. 2).This coincided with the period of the year in which the largest density of snails with body sizes smaller than 1 mm 3 was recorded (up to 15 906 ind / m 3 ).In both the cold and dry seasons, densities of snails smaller than 1mm 3 were below 7000 ind / m 3 and body size spectra then showed lower values of a, but higher values of b, than those estimated in the rainy season (Fig. 2).

DISCUSSION
Overall, our results suggest that body size spectra of both taxonomic groups analyzed in this study can be described through power functions.Nevertheless, these results also indicated that the normalization constants (a) and the scaling exponents (b) describing these relationships varied across seasons of the year.Indeed, the scaling exponents rarely reached the widely expected value of -3/4 (Cyr, 2000).These deviations from the scaling exponent that is expected under energetic balance conditions suggest that the loads of matter and energy in this water body may be changing across seasons, which in turn impact on the local species assemblages and force them to depart from their energetic equivalence equilibrium states.On this issue, it is important to note that the body size spectra of both shrimps and snails attained greater values of a and lower values of b in those seasons when the abundance of individuals in the smaller body size classes increased (dry season for shrimps and rainy season for snails).This suggests that body size spectra of these organisms may be partially regulated by reproductive pulses that occur with changing environmental signals.
Previous studies have shown that reproduction of shrimps in brackish waters is triggered by increases in dissolved oxygen (Attrill & Power, 2000;Viegas, Martinho, Neto, & Pardal, 2007).In our case these conditions occurred in the dry season, coinciding with an increased abundance of shrimps smaller than 1 mm 3 in body volume.This suggests that the recruitment of caridean shrimps occurs in this season, when these organisms might take advantage of the increased inputs of oxygen to reach the energetic levels required to produce offspring.As a consequence of this increased external input of available energy, the value of b dropped below -8/5 (-1.16) in the body size spectrum obtained in dry season.
Reproductive events of snails in brackish wetlands, on the other hand, usually concur with greater loads of organic matter (Antonio et al., 2010;Naranjo-García, 2003;Naranjo-García & Meza-Meneses, 2000).In our study system this occurred during the rainy season, when the rivers that discharge their waters in the lagoon increase their flows and, consequently, the inputs of organic matter that comes from the surrounding terrestrial ecosystems.In this season, the value b in the body size spectrum was below -5/4 (-1.28), which could also  be due to the increased density of small newborn individuals that occurred because a higher availability of food resources.
For both taxonomic groups, the comparisons of body size spectra among seasons also indicated that the values of a declined, while values of b rose, as the number of individuals in the smaller body size classes decreased.
For shrimps this occurred during the cold and rainy seasons, while snails showed decreases in the abundance of smaller individuals during the cold and dry seasons.Such decreases in the abundance of smaller individuals may be due to increased predation intensity of newborn individuals by fish and other organisms, or because environmental conditions after breeding periods are adverse for the survival and development of juvenile stages (Boix et al., 2004;Gosselin & Qian, 1997;Kneib, 1987;Pechenik, 1999;Sola, 1996).Irrespective of the factors that caused these decreases in the abundance of smaller individuals, this may relax the contribution of smaller body size classes to the shape of body size spectra, leading to lower values of a and higher values of b as observed in this study.
Interestingly, the values of b in the body size spectra of shrimps were close to the expected value of -3/4 after the dry season (b = -0.71 in the rainy season and b = -0.75 in the cold season).This suggests that external climatic forces (e.g., changes in environmental conditions) or internal biotic pressures (e.g., increased predation) may pull down the excess of energy within this species assemblage after the reproductive season, hence approaching it to energetic equilibrium states.Scaling exponents of snails, on the other hand, were never close to this expected value of -3/4.Indeed, the body size spectra this group of organisms showed steeper slopes which approximated to -4/4 during both the dry (-0.98) and cold seasons (-1.07).These deviations from the expected -3/4 scaling exponent clearly suggest that this species assemblage is constantly departed from their energetic equilibrium state.However, elucidating what factors caused these deviations in the dry and cold seasons require detailed studies about the responses these species to different environmental factors, which are still lacking for most snails from the Gulf of Mexico.
As final remark it is important to highlight that, although scaling exponents of -3/4 seem to be the rule for body size-abundance relationships at geographical scales (Peters, 1986;Brown et al., 2000;Cyr, 2000), our results suggest that these values may vary idiosyncratically at local smaller scales.Indeed, our study suggests that species assemblages inhabiting ecosystems subjected to cyclical changes in environmental conditions would rarely reach their energy balance because they would not necessarily respond to local energetic constraints.Instead, since the inputs of matter and energy in these ecosystems seem to be controlled by external factors (e.g., seasonality), these species assemblages would be constantly restructured after these changes occur.

Fig. 1 .
Fig. 1.Physicochemical variables (averages ± 2 S.E.) measured at the cold, dry and rainy seasons in the Tamiahua Lagoon, Mexico, including dissolved oxygen (A), salinity (B), chemical oxygen demand (C), temperature of water (D), organic matter in sediments (E) and pH (F).Significant differences between seasons are indicated with different letters (Tukey test critical a = 0.05).Cold season Dry season Rainy Season Cold season Dry season Rainy Season Simple nonlinear regression analyses conducted to determine whether body size spectra of benthic shrimps and snails follow power functions at each season of the obtained from the regression ANOVA performed for each body size spectra, the degrees of freedom of the error (d.f.) in each analysis, the probability value (p) associated to the F-values (critical probability value = 0.05 in all cases), and the amount of variance in data explained by the power function (r 2).The table also indicates the estimated values of normalization constants (a) and scaling exponents (b) of body size spectra functions; values in brackets on the side of a and b indicate the amplitude of the respective 95% confidence intervals estimated by simple nonlinear regression analyses.

Fig. 2 .
Fig. 2. Body size spectra of shrimps (left column) and snails (right column) estimated at the cold (A-B), dry (C-D) and rainy (E-F) seasons.The power functions describing these relationships between the density of individuals and the body volume are provided in the respective panel.

Fig. 3 .
Fig. 3. Linearized body size spectra of shrimps (A) and snails (B) resulting from the multiple linear regression analyses conducted to determine differences in normalization constants and scaling exponents among the cold (solid circles-solid line), dry (empty circles-dashed line) and rainy (solid triangles-dotted line) seasons.

TABLE 1
Coordinates of the eight sites where water quality conditions were assessed and benthic macroinvertebrates were sampled at each season.