Ecología

Meghan I. Zolá-Rodríguez ^a, Mariana Cuautle ^b, *, Marco Daniel Rodríguez-Flores ^c, Citlalli Castillo-Guevara ^b

^a Akumal Monkey Sanctuary & Rescued Animals, Camino a Uxuxubi s/n Predio Santa Pilar Lote 16, 77776 Akumal, Quintana Roo, Mexico

^b Universidad Autónoma de Tlaxcala, Centro de Investigación en Ciencias Biológicas, Km 10.5 Carretera Tlaxcala-San Martín Texmelucan, 90120 San Felipe Ixtacuixtla, Tlaxcala, Mexico

^c Blue Marlin Conservation (Conservation Diver), Sunset Beach Gili Air, Gili Indah, Pemenang, North Lombok Regency, West Nusa Tenggara 83355, Indonesia

*Corresponding author: mcuautle2004@gmail.com (M. Cuautle)

Received: 26 June 2024; accepted: 27 August 2024

Abstract

This study examines the impact of disturbance on the food preferences and dominance of an ant community in a temperate ecosystem in Mexico. The study focused on 2 types of vegetation: native oak forest and induced grassland (disturbed vegetation). Observations were conducted to record the food elements carried by ants to their nests. These data were analyzed using x² tests. Tuna and honey baits were placed near the nests to record the presence of ants in 5-minute periods. We used a binomial model to determine whether the probability of finding an ant foraging at the baits was affected by vegetation type, bait type, and/or ant species. Additional baits were used to determine the ant dominance indices. T-tests and ANOVAs were used to compare dominance indices between vegetation types, baits, and ant species. No significant differences were observed in food preferences between vegetations. However, some species showed a preference for honey (i.e., carbohydrates), which could be limited in ground-level environments. Ants showed a submissive behavior in both vegetation types. This research shows that ants could optimize their nutrient intake, enabling them to survive efficiently even when facing disturbances, instead of increasing dominance. 

Keywords: Ant nest; Dominance index; Feeding habits; Compensation hypothesis; Carbohydrates; Proteins

© 2024 Universidad Nacional Autónoma de México, Instituto de Biología. Este es un artículo Open Access bajo la licencia CC BY-NC-ND

(http://creativecommons.org/licenses/by-nc-nd/4.0/).

Impacto del disturbio en las preferencias alimentarias y dominancia de las hormigas (Hymenoptera: Formicidae) en un bosque templado de México

Resumen

Este estudio examina el impacto del disturbio en las preferencias alimentarias y dominancia de una comunidad de hormigas en un ecosistema templado en México, en bosque de encino nativo y pastizal inducido (vegetación perturbada). Se registraron los alimentos transportados por las hormigas a sus nidos. Estos datos fueron analizados utilizando pruebas de x². Se colocaron cebos de atún y miel cerca de los nidos para registrar la presencia de hormigas. Utilizamos un modelo binomial para determinar si la probabilidad de encontrar una hormiga en los cebos se veía afectada por la vegetación, cebo o especie de hormiga. Los índices de dominancia se determinaron usando cebos. Se emplearon pruebas t y Anova para comparar los índices de dominancia entre tipos de vegetación, cebos y especies de hormigas. No hubo diferencias significativas en las preferencias alimentarias entre tipos de vegetación, pero algunas especies mostraron una preferencia por la miel (carbohidratos), que podría ser un recurso limitado a nivel del suelo. Las hormigas mostraron un comportamiento sumiso en ambos tipos de vegetación. Esta investigación muestra que las hormigas podrían optimizar su ingesta de nutrientes, permitiéndoles sobrevivir bajo condiciones de disturbio, en lugar de aumentar su dominancia.

Palabras clave:Nidos de hormigas; Índice de dominancia; Hábitos alimenticios; Hipótesis de compensación; Carbohidratos; Proteínas

Introduction

Land use change stands as one of the primary factors contributing to global change and its adverse effects on biodiversity (Ellis et al., 2010; Foley et al., 2005; Sala et al., 2000). Land use change has reduced biodiversity through the loss, modification, and fragmentation of habitats; degradation of soil and water; and overexploitation of native species (Foley et al., 2005), and its effects depend strongly on the type, severity, frequency and timing of disturbance (Foley et al., 2005; White & Jentsch, 2001). Ants are one of the most dominant insects both ecologically and numerically (Rico-Gray & Oliveira, 2007; Schultheiss et al., 2022; Toro et al., 2012) that participate in different ecological processes —e.g., nutrient recycling, soil formation, decomposition, seed dispersion (Toro et al., 2012). Ants are model organisms to study the effect of disturbance because they respond to environmental change (Agosti et al., 2000; Andersen, 2000). It is often observed that disturbances favor behaviorally dominant ants, including invasive species. Vonshak and Gordon (2015), observed that native ant richness was highest in natural habitats, and alien species richness was highest in urban habitats, along an urban-rural gradient in the San Francisco Bay area. Nevertheless, the specific effect will depend on factors like the type of disturbance, including fires, floodings, or even treefall gaps (Cerdá et al., 2013). Hoffman and Andersen (2003) studied the response of ant functional groups to disturbance finding that the dominant Dolichoderinae and “hot climate specialists” tend to be favoured by low levels of disturbance. The opportunist and generalized Myrmicinae have wide habitat tolerances but are sensitive to competitive interactions, while “cryptic species” and “specialist predators” are highly sensitive to disturbance.

Feeding habits and foraging strategies represent vital life history traits among ants, playing a significant role in their ecological dynamics. These traits are contingent upon the availability of resources within their respective habitats (Andersen, 2000; Davidson, 2005), and the ensuing competition for those resources (Arnan et al., 2012). The feeding habits of ants play a pivotal role in determining their ecological function within the ecosystem (Spotti et al., 2015). However, despite their significance, the food habits and preferences of ants remain relatively understudied and represent one of the least understood aspects of their biology (Houdria et al., 2015). 

The effects of disturbance (e.g., land use change) on ant feeding habits (especially foraging strategies) have not been studied in depth (Castillo-Guevara et al., 2019; Hernández-Flores et al., 2016; Radnan et al., 2018), and the way ants respond to habitat complexity can provide clues as to how ants could respond to disturbance. When the external disturbing factors are of high intensity (natural-fires or anthropogenic-deforestation), the disturbances can initiate a directional regression manifested as gradual or rapid simplification of the horizontal and vertical structure of a community, leading to the replacement of complex communities by a simpler one (Łaska, 2001). The level of habitat complexity can influence ant foraging patterns. For instance, higher complexity may decrease interspecific interactions and ant recruitment or minimize the trade-off between resource discovery and dominance (Parr & Gibb, 2012).

Examining the effect of soil surface complexity on food exploitation in the context of change from grassland to shrubland in Australia, Radnan et al. (2018) discovered that substrate complexity (wood debris, leaf litter, or no substrate) influenced the discovery time, ant size, and monopolization index of tuna and honey baits within testing arenas. However, at a larger scale of vegetation type, the effect was not observed. Similarly, Castillo-Guevara et al. (2019) found that the dominance level was similar in natural oak vegetation in comparison to an agricultural area in a temperate community; however, ant foraging strategies differed between the 2 communities.

Competition has been regarded as the primary factor influencing the structure of ant communities (Cerdá et al., 2013). Through the use of aggressive strategies, ants exert influence on the spatial distribution, abundance, and behavior of other ants. These strategies encompass a range of tactics, including the deployment of repellent chemicals and the establishment of territories (Cerdá et al., 2013). Moreover, ants can establish interspecific dominance hierarchies primarily based on variations in food collection behavior and aggressiveness. While there are various proposals for categorizing such hierarchies, the classification proposed by Vepsäläinen and Pisarski (1982) and Savolainen et al. (1989), as described by Cerdá et al. (2013), is considered the most well-defined hierarchical system from an ecological perspective. This classification system is founded on the aggressive behavior of ants and its impact on other ant species, and includes 3 categories: the dominant, subdominant and subordinate category. Dominant species exhibit highly aggressive behavior, exhibit numerical predominance over other species, fiercely defend their territories, and establish mutual exclusivity within their ecological communities. Subdominant species, while not actively defending territories, display a remarkable propensity for aggressively safeguarding their food resources (Cerdá et al., 2013). Lastly, subordinate species, characterized by small colonies devoid of recruitment systems, employ strategies to avoid physical confrontations with other colonies and species. Nonetheless, they exhibit a strong commitment to defending their nests against potential threats (Cerdá et al., 2013). Different foraging strategies, such as the subordinate species’ ability to discover resources before dominant species (dominance-discovery trade-off, Fellers, 1987) or their capacity to forage across a broader temperature range compared to dominant species —dominance-thermal tolerance trade-off (Fellers, 1989), can contribute to the coexistence of ants within a community.

The present study focuses on the analysis of ant feeding habits and dominance hierarchy in response to land use change within a temperate ecosystem located in central Mexico. A habitat with the native oak forest vegetation (complex habitat), was compared to a nearby area where the land use had been altered to induced grassland (simplified habitat). The study hypothesized: i) it is expected that in the oak forest, food items transported to the nest will be more varied than in the induced grassland. This variation at the species level is also anticipated due to the change in land use. The conversion from a complex habitat (oak forest) to a simplified habitat (induced grassland) will likely provide fewer resources variety for the ants; ii) it is expected that as a result of the disturbance in the induced grassland, ants will be more generalized in low heterogeneous habitats due to the dominance of generalist ant species. Additionally, the arrival times at the baits will be shorter due to the absence of leaf litter; iii) a lower dominance index is expected in the induced grassland than in the oak forest, due to the presence of subdominant species and the modification of the dominance hierarchy system by the land use change.

Materials and methods

The study was conducted within “Flor del Bosque” State Park, located in the municipality of Amozoc de Mota, in the state of Puebla. The coordinates of this protected reserve are 19°00’00”-19°01’50” N, 98°20’35”-98°20’53” W. The State Park encompasses an area of 664.03 hectares, characterized by altitudes ranging from 2,225 to 2,400 m asl. The annual average temperature in this region fluctuates between 14 °C and 16 °C, with the majority of rainfall occurring during the summer months, contributing to an average annual precipitation of 750 to 950 mm. It is important to note that the park experiences a distinct dry season lasting approximately 6 months, from November to April, as reported by Costes et al. (2006). The native vegetation of the reserve primarily consists of oak forest [Quercus castanea (Née), Q. laurina (Bonpl.), Q. laeta (Liebm.) (Fagaceae)]. However, human disturbances have led to the presence of induced grassland areas and, to a lesser extent, eucalyptus plantations [Eucalyptus spp. L´Hér) (Myrtaceae)] within the park, as documented by Costes et al. (2006). Given that the disturbed vegetation within the reserve primarily consisted of induced grassland, our research aimed to establish a meaningful comparison between native and altered vegetation.

Ant communities were surveyed once a month during specific periods (April, August and October 2015, and January to March 2016). The survey methodology involved the establishment of 6 transect plots (400 m × 20 m), with 3 plots located in the oak forest and 3 in the induced grassland. The spatial distribution of these transects can be observed in Figure 1. To locate ant nests, the transects were traversed, and various substrates such as leaf litter, stones, trunks, and branches were lifted and examined. Additionally, a total of 5 tuna baits and 5 honey baits were placed along the transects at 10 m intervals to attract ants and to follow them to locate their nest. The nests found during the survey were georeferenced for accurate spatial documentation. A 5-minute observation period was designated for each nest, wherein all food items transported by the ants to their nests were recorded. These observations were conducted between 9:00 a.m. and 3:00 p.m. on the sampling days. Each recorded item was subsequently categorized into one of the following groups: 1) plant elements, encompassing any component of plants except for seeds, 2) seeds, and 3) arthropods. After each observation period, 2 bait samples, enclosed within Petri dishes, were positioned in close proximity to the nests (approximately 10 cm), maintaining a distance of approximately 5 centimeters from each other. Various studies exploring the feeding habits, preferences, and foraging strategies of ants have employed diverse bait types. Typically, these baits consist of carbohydrates such as honey and other sources representing proteins such as tuna (Houdria et al., 2015; Lynch et al., 1980; Spotti et al., 2015; Trigos-Peral et al., 2016). The bait materials employed in the study consisted of honey and tuna, which were consistently provided in a standardized amount of one tablespoon, equivalent to approximately 5 grams. To quantify the number of individual ants on each type of food resource, we placed baits and observed them for 5 min. Ants responded very quickly to the baits; therefore, the observation time was limited to 5 min. Within this duration, the number of individuals was recorded for each bait type. After the observation periods, 1 to 3 ants were collected from each nest using a vacuum cleaner or tweezers and then preserved in Eppendorf tubes containing 70% ethanol. In the laboratory, the ants were separated, mounted, and identified to the genus level using the key by Mackay and Mackay (1989). In certain instances, the species identification was feasible by comparing the collected specimens with those present in the Entomological Collection of the Universidad de las Américas Puebla (UDLAP).

To analyze the food preferences of the ant community, we considered 3 factors. First, we determined whether the number of foragers carrying plants, seeds, or arthropods to the nest differed between the whole ant communities of each vegetation type or among ant species. We used the appropriate contingency tables and chi-square tests for such comparisons. Second, we determined whether the presence of foragers at each bait during each minute of the 5-minute period was influenced by vegetation type, bait type (honey or tuna), ant species or interactions among these factors. We used a binomial model to evaluate this response variable, with “ant presence” coded as 1 and “absence” as 0. To select the most relevant explanatory variables, we employed a stepwise forward approach. We began with the simplest model and successively added each response variable. The significance of adding each variable was assessed using an x² test, comparing the previous model with the new model. If a variable was found to be statistically significant, it was retained in the model. Otherwise, it was removed, and we proceeded with the next explanatory variable. This stepwise procedure allowed us to build the final model with the most significant factors explaining the presence of foragers at each bait. No overdispersion was detected in the selected model. Post-hoc analyses were conducted when significant differences were detected. The statistical analysis of data was performed using R software (R Core Team, 2022).

We evaluated the dominance hierarchy of the ants in each vegetation type (oak forest and induced grassland), ant species, and bait types using the same transects as for the detection of the nests (Fig. 1), during October 2015 and January, February and March (2016). At each transect, 9 sampling points were established on the ground, spaced 10 m apart. Each sampling point consisted of a pair of Petri dishes with baits: one with honey and one with tuna, placed less than 5 cm apart. This resulted in a total of 540 baits (5 replicates × 6 transects × 9 sampling points × 2 baits; 2 samplings in February). These types of baits have been widely used for ant dominance hierarchy studies (Dáttilo et al., 2014; Parr & Gibb, 2012; Trigos-Peral et al., 2016).

For the first 3 sets of bait, we recorded the arrival times of the ants during a one-hour observation period. Most of the ant species were identified in the field. For the ants that were not identified, a few individuals (2 to 3 specimens) were collected. In the case of the remaining 6 sets of bait, they were filled with water. Ants that fell onto the Petri dish after 2 h were collected. The presence of water did not deter the ants from visiting, but it allowed us to determine which ants had been attracted to the baits. The baits were placed and retrieved from different transects in the field between 9:00 a.m. and 5:00 p.m. The order of the transects was changed on different sampling days to mitigate the potential impact of the time of day. Ant specimens collected from the baits were preserved in 70% ethanol and transported to the laboratory for further identification.

Figure 1. Localization of the transects used in the study zone (mapped by Luna F.).

We chose numerical dominance to assess the dominance hierarchy of the ants —i.e., ordering of ant species based on their numerical or behavioral dominance by vegetation types (Andersen, 1992; Cerdá et al., 1997; Stuble et al., 2017), ant species and bait type. As numerical and behavioral dominance are highly correlated, this method has been well-established and documented in the ant literature (Dáttilo et al., 2014; Dejean & Corbara, 2003; Parr, 2008; Parr & Gibb, 2012; Santini et al., 2007). This method indicates which species are consistently present at the baits, and which ones dominate the baits numerically and thus monopolize them (Parr, 2008). We represented numerical dominance using the numerical dominance index (DI) for each morphospecies calculated by the formula: DI = (Di)/(Di + Si), where, Di is the number of baits monopolized by the species of ant i, and Si is the number of baits that the species of ant i used but did not monopolize. Baits were considered to be monopolized when more than 5 individuals (workers and/or soldiers) of the same morphospecies were using the resource without the presence of other morphospecies. This measure (more than 5 individuals) takes into consideration that in temperate climates, ants are less abundant, and recruitment is considered weaker than in tropical environments where the index has been widely more (Santini et al., 2007). Therefore, dominant morphospecies are those that find and monopolize a larger proportion of the food resources in a given environment. The value of the index ranges from 0 (completely submissive species) to 1 (totally dominant species) and is similar to the “monopolization index” used in other studies (Dáttilo et al., 2014; Fellers, 1987; Parr & Gibb, 2012; Santini, et al., 2007). In this study, ant species with a DI lower than 0.5 were classified as submissive.

The arrival times of the ants were compared using survival curves using the survival package of the R Software program. The one-hour and 2-hour baits were used to calculate the ID of the ant species. To compare the ID between vegetation and bait types, t-tests were performed. To compare the DI between ant species, a one-way ANOVA test was performed after applying a square root transformation to meet the normality requirements. These analyses were performed using the program STATVIEW 5.0 (Abacus Concepts Inc., 1996).

Results

Fourteen morphospecies of ants were recorded in this study belonging to 11 ant genera (Table 1). For the analysis of the results from the nests and the baits placed near them, 2 morphospecies of Pheidole were identified. However, due to the low number of records for each morphospecies across different vegetation types, Pheidole sp. 1 and Pheidole sp. 2 were combined for the analysis. In the oak forest, 2 out of 19 records corresponded to Pheidole sp. 1, and in the grassland, 4 out of 15 records corresponded to Pheidole sp. 2. Since Pheidole species often share similar feeding habits and are functionally similar (Agosti, 2000; Andersen, 2000), these 2 morphospecies were grouped to increase statistical power, and the analysis was conducted at the genus level. In contrast, another Pheidole species, Pheidole sp. 3, was recorded in the baits used to determine the DI and was analyzed independently due to its different context of occurrence. The rest of the analyses were conducted at the species level.

In the case of the nests, we found 9 morphospecies distributed among 34 nests, 19 of these were located in the oak forest and 18 in the induced grassland (Table 2). No significant differences were observed in the types of food elements those individual ants transported to their nest when comparing different vegetation types (x² = 2.13, df = 2, p = 0.34). In the oak forest considering all ant species, there were no significant differences in the number of individuals carrying elements from different categories (x² = 0.03, df = 1, p = 0.86) (Fig. 2). Of the genera of ants observed carrying food elements to their nests, only 2 genera were observed with a single type of food resource: Prenolepis imparis (Say) only carried arthropods and Dorymyrmex insanus (Buckley) only seeds.

Table 1

List of ant species recorded in this study in “Flor del Bosque” State Park in Puebla, Mexico.

Formicidae/Subfamily/Ant species
Formicidae Latreille, 1809
Dolichoderinae Forel, 1878
	Leptomyrmecini Emery, 1913
		Dorymyrmex Mayr, 1866
			Dorymyrmex insanus (Buckley, 1866)
		Linepithema Mayr, 1866
			Linepithema dispertitum (Forel, 1885)
Dorylinae Leach, 1815
	Dorylini Ashmead, 1905
		Labidus Jurine, 1807
			Labidus coecus Latreille, 1802
Formicinae Latreille, 1809
	Camponotini Forel, 1878
		Camponotus Mayr,1861
			Camponotus rubrithorax Forel, 1899
	Lasiini Ashmead, 1905
		Nylanderia Emery, 1906
			Nylanderia austroccidua Trager, 1984
		Prenolepis Mayr, 1861
			Prenolepis imparis Say, 1836
Myrmicinae Lepeletier, 1835
	Attini Smith, 1858
		Pheidole Westwood, 1839
			Pheidole sp. 1
			Pheidole sp. 2
			Pheidole sp. 3
	Crematogastrini Forel, 1893
		Temnothorax Mayr, 1861
			Temnothorax augusti Baroni Urbani, 1978
	Pogonomyrmecini Ward, Brady, Fisher & Schultz, 2014
		Pogonomyrmex Mayr, 1868
			Pogonomyrmex barbatus (Smith, F., 1858)
	Solenopsidini Forel, 1893
		Monomorium Mayr, 1855
			Monomorium ebenium Forel, 1891
Pseudomyrmecinae, Smith, M.R., 1952
	Pseudomyrmecini Smith, 1952
		Pseudomyrmex Lund, 1831
			Pseudomyrmex pallidus Smith, F., 1855

Conversely, ants of the genus Pheidole spp. carried a greater quantity of plant elements compared to the other categories (x² = 10.4, df = 2, p = 0.006); only 2 of these individuals were observed carrying arthropods, and none were recorded carrying seeds. No significant differences were found for D. insanus (x² = 2, df = 2, p = 0.37), P. barbatus (x² = 1.6, df = 2, p = 0.450), nor P. imparis (x² = 4, df = 2, p = 0.135) among the 3 food categories. Only P. barbatus carried all 3 types of food (Fig. 2). 

Considering all ant species, there were no significant differences in the number of individuals carrying different food elements in the induced grassland (x² = 0.66, df = 1, p = 0.41, Fig. 2). Among the observed ant species that transported food types to their nests, P. barbatus was found to carry both plant elements and seeds. Additionally, Dorymyrmex insanus and Pheidole spp. were observed carrying food from all 3 categories (Fig. 2). Pheidole spp. (x² = 10.33, df = 2, p = 0.006) and P. barbatus (x² < 30.33, df = 2, p < 0.001), carried a higher number of plant elements. For D. insanus no significant differences were found in the number of elements of each type (x² < 3.77, df = 2, p = 0.15) (Fig. 2). 

Figure 2. Number of individuals of the different ant species recorded carrying any food resource to their nest. DORY: Dorymyrmex insanus, PHEI: Pheidole spp., POGO: Pogonomyrmex barbatus, PREN: Prenolepis impairs.

In the induced grassland, 5 ant species were attracted to both types of baits, and 1 species was found in the tuna baits. In the oak area, 5 ant species were observed in both types of baits, and 1 species was attracted to the honey bait (Fig. 3). The probability of presence of a forager during a given minute of the observation period depended on the vegetation type (x²= 8.73, df = 1, p = 0.03), the ant species (x²= 59.64, df = 7, p < 0.001), the interaction between the ant species and the bait (x²= 57.96, df = 8, p < 0.001), and the interaction between the vegetation type, the ant species, and the bait (x²= 26.12, df = 7, p < 0.001). The Tukey contrast test for the vegetation*species*bait interaction showed (Z > 3.83, p < 0.042) that the presence of Linepithema dispertitum (Forel) on the honey bait, in the oak forest, was more likely than the presence of Pheidole spp. on either the tuna or honey bait in the oak forest, Pheidole spp. on the honey bait in the induced grassland, and P. impairs on the honey bait in the oak forest (Fig. 3). The probability of Camponothus rubrithorax (Forel) being present on the honey bait in the induced grassland, was greater than that of D. insanus in either the tuna or honey bait in the induced grassland, Pheidole spp. on either the tuna or honey bait in the oak forest, as well as its presence in the honey bait in the induced grassland and P. imparis in the honey bait in the oak forest (Fig. 3). The probability that P. barbatus was present on the honey bait in the induced grassland, was greater than the presence of Pheidole spp. in either the tuna or honey bait and P. imparis in the honey bait in the oak forest (Fig. 3). Finally, the presence of Pheidole spp. on the tuna bait was more likely than its presence on the honey bait in the induced grassland (Fig. 3).

Figure 3. Probability of finding an ant forager on the honey or tuna baits at the nest of the different ant species registered in the Oak Forest (OF) and the Induced Grassland (IG) in the 5-minute period. CAMPO: Camponotus rubrithorax, DORY: Dorymyrmex insanus, LABI: Labidus coecus, LINE: Linepithema dispertitum, MONO: Monomorium ebenium, PHEI: Pheidole spp., POGO: Pogonomyrmex barbatus, PREN: Prenolepis imparis.

Figure 4. Survival curves depicting the arrival times of ant foragers on the baits (tuna or honey) during a one-hour observation.
No significant differences were found in the arrival times among ant species.

No significant differences were found in the arrival times between the vegetation types (x² = 0.4, df = 1, p = 0.5). Moreover, no significant differences were found in the arrival times among the different ant species (x² = 7.6, df = 7, p = 0.4) or the type of bait (x² = 0.01, df = 1, p = 0.9; Fig. 4).

In the dominance index baits 12 ant morphospecies were recorded and different ant morphospecies were recorded in each of the vegetation types (Table 3). When comparing the average dominance index (DI) by vegetation type (mean ± SE, n; oak forest DI = 0.372 ± 0.071, 33; induced grassland DI = 0.251 ± 0.052, 44), was not statistically different (t= 1.091; df = 75; p = 0.2787), which would suggest a submissive behavior in both types of vegetation. Average dominance index by ant species indicates that M. ebenium, P. imparis and Pheidole sp. 3 behave as submissive in both vegetation types (t = 0.590, df = 20, p = 0.5617; t = -0.237, df = 10, p = 0.8171; t = -0.315, df = 11, p = 0.7590, respectively) (Table 4). It was not possible to compare the DI between vegetation types for Temnothorax sp., T. augusti, C. rubrithorax, N. austroccidua, L. dispertitum, D. insanus or P. pallidus, because they were only present in one of the vegetation types (Table 4). 

The average dominance index for each ant species was low suggesting that all of them displayed submissive behaviors (F_8,68 = 1.949, p = 0.06). Nylanderia austroccidua and L. dispertitum had the highest dominance indices, although it should be noted that N. austroccidua only had one record and L. dispertitum had 3 records. Additionally, Monomorium ebenium showed a tendency to behave as dominant (Table 4).

When comparing the mean dominance index by bait type, no significant differences were found (t = 1.023, df = 110, p = 0.3088; mean ± SE, n; tuna = 0.382 ± 0.060, 55; honey = 0.295 ± 0.053, 57). This indicates that ants exhibited submissive behavior in both types of bait.

According to the average dominance index per ant species in relation to the bait type, Pheidole sp. 3, P. imparis, M. ebeninum, C. rubrithorax, and D. insanus exhibited a submissive behavior in both types of baits (t = 0.888, df = 13, p = 0.3907; t= 0.988, df = 15, p = 0.3388; t = 0.189, df = 32, p = 0.8513; t = -0.918, df = 24, p = 0.3675, t = 2.390, df = 5, p = 0.0624, respectively) (Table 5) and L. dispertitum exhibited a dominant behavior in both types of baits (t = 0.421, df = 3, p = 0.7021).

Temnothorax sp. and P. pallidus were not compared statistically due to the low incidence registered; nonetheless, they presented a submissive behavior index in both types of bait. Neither T. augusti nor N. austroccidua were included in this analysis, as they were only present on the honey bait. In none of the species was there a difference in behavior between the baits.

Table 2

Number of nests found in the oak forest and the induced grassland.

Species	Oak forest	Induced grassland
Camponotus rubrithorax	0	3
Dorymyrmex insanus	4	6
Labidus coecus	1	0
Linepithema dispertitum	2	0
Pheidole spp.	7	6
Pogonomyrmex barbatus	1	2
Prenolepis imparis	3	0
Monomorium ebenium	1	1

Discussion

This study found that certain ant species (e.g., C. rubrithorax, L. dispertitum, P. barbatus) exhibited a preference for honey within specific vegetation types when compared to the rest of the ant community. The compensation hypothesis (Davidson, 2005; Kaspari & Yanoviak, 2001) predicts that the attractiveness of a nutrient to an organism is higher the more limiting it is. It is possible that this result is related to the fact that sugar is generally less available on the ground strata than protein (Kaspari & Yanoviak, 2001; Kaspari et al., 2012).

However, when considering the overall community level, no noticeable food preferences were found between the ant communities in the oak forest and the induced grassland vegetation. This lack of preference was evident both in the items carried to their nests and the resources provided on the baits. There were no significant differences in arrival times at the baits between vegetation types. This similarity can be attributed to ants experiencing soil-level heterogeneity similarly in both oak forests and induced grasslands (see below).

The ant communities in both habitats were composed of ant species displaying submissive behavior. This characteristic could be attributed to stress factors such as low temperatures in the oak forest and disturbances in the induced grassland, which potentially reduce ant competition. In the oak forest, low temperatures may favor the presence of cold climate specialists that can forage at low temperatures without the need to display a dominant behavior. In the induced grassland, disturbances may favor the presence of generalist species that take advantage of the absence of dominant species such as the dominant Dolichoderinae (see below). This study emphasizes the significance of ant species’ response to their environment and their adaptability in dealing with disturbances.

Our findings do not provide support for vegetation-scale disparities in food preferences, first and second hypotheses, which proposed greater differences in the food items being carried to the nest in the oak forest and reduced food preferences in the induced grassland. There were no differences between vegetation types in the number of seeds, plants, or arthropods taken by the ant foragers to their nests. Similarly, there was not a preference for a specific food bait (honey or tuna) between the oak forest and induced grassland ant communities. These findings are consistent with previous studies conducted by Radnan et al. (2018) and Castillo-Guevara et al. (2019), which did not identify differences in food preferences or foraging strategies between natural and disturbed vegetation at the community level. However, our observations did reveal that certain ant species exhibited preferences for specific types of food, within particular vegetation types, in accordance with the first hypothesis and second hypothesis but at the species level, which proposed an influence of ant species on food preferences. The specific preferences of these ant species are discussed in detail below.

Table 3

Visit frequency (i.e., number of foragers) per species in each vegetation type and functional groups they belong to (Andersen, 2000). GM = Generalized Myrmicinae, CCS = cold climate specialist, SC = subordinate Camponotini, TCS = tropical climate specialist, O = opportunist.

Subfamily	Specie	Funtional group	Oak forest	Induced grassland
Myrmicinae	Monomorium ebenium	GM	753	167
Myrmicinae	Pheidole sp. 3	GM	14	236
Myrmicinae	Temnothorax sp.	CCS	5	0
Myrmicinae	Temnothorax augusti	CCS	0	1
Formicinae	Prenolepis imparis	CCS	289	11
Formicinae	Camponotus rubrithorax	SC	0	351
Formicinae	Nylanderia austroccidua	TCS	0	30
Dolichoderinae	Linepithema dispertitum	CCS	190	0
Dolichoderinae	Dorymyrmex insanus	O	0	26
Pseudomyrmecinae	Pseudomyrmex pallidus	TCS	0	10
Total			1,251	832

Table 4

Average numerical dominance index (DI) (mean ± SE, n) per ant morphospecies. (-) Unregistered species.

Specie	Oak forest	Induced grassland	ID
Monomorium ebenium	0.4 ± 0.1, 13	0.3 ± 0.1, 9	0.4 ± 0.0, 22
Pheidole sp.3	0.2 ± 0.2, 4	0.2 ± 0.1, 9	0.2 ± 0.1, 13
Temnothorax sp.	0 ± 0, 3	–	0 ± 0, 3
Temnothorax augusti	–	0, 1	0 ± 0, 1
Prenolepis imparis	0.3 ± 0.1, 10	0.5 ± 0.5, 2	0.3 ± 0.1, 12
Camponotus rubrithorax	–	0.1 ± 0.0, 15	0.1 ± 0.0, 15
Nylanderia austroccidua	–	1, 1	1, 1
Linepithema dispertitum	0.7 ± 0.1, 3	–	0.7 ± 0.1, 3
Dorymyrmex insanus	–	0.2 ± 0.1, 4	0.3 ± 0.1, 4
Pseudomyrmex pallidus	–	0 ± 0, 3	0 ± 0, 3

Pheidole spp., in the induced grassland, were observed carrying a greater quantity of plant elements to their nests in both vegetation types. However, it exhibited a preference for the tuna bait in induced grassland, possibly indicating a supplementary dietary preference. Given the extensive diversity within the Pheidole genus, which encompasses 900 species described worldwide (Wilson, 2003), it is not feasible to categorize them based on a specific food habit and our results suggest that it is omnivorous (Table 6). Nevertheless, these results should be interpreted with caution due to the grouping of species.

Camponotus rubrithorax and P. barbatus, in the induced grassland, as well as L. dispertitum, in the oak forest, were more frequently observed foraging the honey bait compared to other ant species. Camponotus (Mayr) is a genus known for its nectarivorous habits and consumption of other sweet secretions, such as honeydew which coincides with its preference for honey found in the present study (Nettimi & Iyer, 2015). In this study, L. dispertitum predominantly consumed honey, although it is a generalist forager species capable of consuming other types of food as well (Table 5), in the study site it is found exclusively in the oak forest (Cuautle et al., 2016).

Despite Pogonomyrmex (Mayr) has been recognized as a granivorous genus (Pirk & López-de Casenave, 2014), this study provided new insights into the foraging preferences of P. barbatus, revealing a notable inclination towards carbohydrates and plant elements within the induced grassland. Moreover, we even registered individuals transporting arthropods to their nests. These findings (Table 5) strongly indicate that certain species within the Pogonomyrmex genus exhibit a generalist foraging behavior. Other species such as D. insanus, which used resources more intensively in the induced grassland, showed no preference for any of the baits or items taken to the nest, which coincides with the generalist forager behavior, observed in open and disturbed habitats by Cuezzo and Guerrero (2012).

Table 5

Average numerical dominance index (DI) (mean ± SE, n) per ant species in both bait types (tuna, honey). (-) Unregistered species.

Species	Tuna	Honey
Monomorium ebenium	0.4 ± 0.1, 18	0.4 ± 0.1, 16
Pheidole sp.3	0.4 ± 0.2, 7	0.1 ± 0.1, 8
Temnothorax sp.	0, 1	0 ± 0, 2
Temnothorax augusti	–	0, 1
Prenolepis imparis	0.5 ± 0.1, 7	0.3 ± 0.1, 10
Camponotus rubrithorax	0.1 ± 0.0, 14	0.2 ± 0.0, 12
Nylanderia austroccidua	–	1, 1
Linepithema dispertitum	0.7 ± 0.2, 2	0.6 ± 0.3, 3
Dorymirmex insanus	0.6 ± 0.3, 3	0 ± 0, 4
Pseudomyrmex pallidus	0 ± 0, 3	0, 1

Table 6

The biology and ecology of the ant genera found in this study, in the Flor del Bosque State Park (Agosti et al., 2000; AntWiki, n.d.).

Genera	Microhabitat	Food habits
Camponotus	Ground nesting, decaying wood and in trees	Generalist foragers
Dorymyrmex	————–	Generalist foragers
Linepithema	————–	Generalist foragers
Labidus	Epigeous,bivouacs	Predators (Army ants)
Monomorium	————–	Generalized foragers, harvesters
Nylanderia	Nest in leaf litter, soil, or in rotten wood	Generalist foragers
Pheidole	Soil or decaying wood	Granivores or omnivores
Pogonomyrmex	Ground nesting	Generalist foragers and granivores
Prenolepis	————–	Generalist predators
Pseudomyrmex	Mostly arboreal (nesters and foragers), few epigaeic	Generalized predators, visit extrafloral nectaries
Temnothorax	Nesting in ground, and under stones, in wood, and in trees	Generalized foragers and parasites

The previous results (e.g., C. rubrithorax, Pheidole spp., D. insanus) support the reported food habits of specific ant genera. Nonetheless, it is also possible to interpret these results as ants taking nutrients from the baits that are not currently available or that are limited within their community to supplement their diet (Compensation hypothesis) (Davidson, 2005; Kaspari & Yanoviak, 2001). In this study, significant differences in ant species presence on baits were primarily associated with the presence of ant species on honey baits. This finding aligns with the compensation hypothesis (Davidson, 2005; Kaspari & Yanoviak, 2001), which posits that the utility of a resource remains constant across species and varies solely with availability. The hypothesis predicts a singular limiting resource that is locally in shortest supply. Consequently, habitats with relatively high protein availability should attract ants more inclined towards carbohydrates, and vice versa. Previous research has demonstrated that litter ant communities are limited by carbohydrates, whereas ant arboreal communities face protein limitations in tropical regions (Kaspari & Yanoviak, 2001; Kaspari et al., 2012). Although this study was conducted in a temperate environment, the results might be related to the usual scarcity of sugar in the ground strata compared to protein.

Our results do not align with the second hypothesis, which speculated faster arrival times in the induced grassland. We did not find significant differences in arrival times among vegetation types, ant species, or bait types. Ants could be experiencing similar heterogeneities at the soil level between the oak forest and induced grassland. While it is expected that the presence of more litter in the oak forest could hinder ant movement, within the induced grassland, the vegetation morphology itself (long grass species) could be interfering with ant movement. For instance, Hernández-Flores et al. (2016) observed that the foraging performance of P. barbatus was reduced due to the presence of herbaceous vegetation in plots where regeneration after grazing was permitted.

Dominance hierarchy. Disturbance typically favors the presence of generalist and opportunistic species, hence, we expected to find a lower dominance index within the induced grassland (third hypothesis). Nonetheless, we did not find enough evidence to support this hypothesis as we registered ant communities consisted of submissive species in both habitats. Additionally, the ants did not show a tendency to dominate a specific resource, neither at the bait level nor at the species level, suggesting a lack of food preferences. These results are consistent with the findings of Castillo-Guevara et al. (2019), who registered no significant differences in dominance indices between a native oak forest and an altered agricultural land. Moreover, the authors reported intermediate to low values of dominance within the ant communities of each vegetation type. Our findings could be attributed to the presence of other factors that could have overshadowed the role of competition in shaping the organization of ant communities. For example, a review by Parr and Gibb (2012), which encompassed data from 3 continents, indicated that the trade-off between discovery and dominance occurs primarily when parasitoids are present. In environments without parasitoids, species with high discovery abilities tend to also be dominant (Parr & Gibb, 2012). While our study did not specifically assess the discovery-dominance trade-off, the observed lack of differences in arrival times and similar dominance values among the ants suggest that dominance was not a prominent factor in our study sites. It is possible that factors such as low temperatures in the oak forest, or disturbance in the induced grassland, may have played a role in relaxing dominance. According to Andersen (2000), in disturbed vegetation like the induced grassland, it is expected to observe the presence of subdominant ant species that exploit the absence of dominant species from the native vegetation, such as dominant Dolichoderinae. The findings obtained during this study align with this prediction and provide support for it. Notwithstanding, the absence of ants with high dominance indices in the studied communities does not necessarily indicate a complete lack of dominance hierarchy. In each community, the ants can still be ordered based on their DI. For example, in the oak forest, L. dispertitum had the highest DI value (0.7 DI), M. ebenium had an intermediate value (0.4 DI) and Temnothorax spp. had the lowest value (0.0 DI).

Land use change did not seem to influence food preferences or foraging strategies at the community level. However, we observed an effect at species level, indicating that individual ant species exhibited specific food preferences. Carbohydrates could be the limiting resource in the oak forest and the induced grassland litter ant communities, as some ant species showed preference for honey baits. Although competition is typically considered a key factor in understanding food preferences among ants, it is noteworthy that both, the natural oak forest and induced grassland, were predominantly populated by submissive ant species. Therefore, it appears that other factors instead of competition may be playing a role in shaping food preferences within these communities. In conclusion, our study reveals that ant species may exhibit preferences for specific foods, which could be limited in their environment. The ability of ants to respond to available resources enables them to optimize their nutrient intake, as well as adapt and persist under variable conditions, including disturbances. Understanding the distinct dietary preferences and foraging strategies of ant species within functional groups will provide valuable insights into their ecological roles and potential impacts on ecosystem dynamics. Such investigations would enhance our ability to predict the responses of ants to diverse forms of disturbances in an anthropized world.

Acknowledgements

We appreciate the assistance provided by the authorities of “Flor del Bosque” State Park, coordinator Enrique Martínez Romero (M.S.) and director Mario Alberto Romero Guzmán (MVZ). We would also like to thank Florencio Luna Castellanos for his support with the fieldwork. This study was financed by Consejo Nacional de Humanidades Ciencias y Tecnologías (Conahcyt) as part of a grant awarded to Mariana Cuautle (223033).

References

Abacus Concepts Inc. (1996). Abacus Concepts, Stat View Reference. Berkeley, California. 

Agosti, D., Majer, J. D., Alonso, L. E., & Schultz, T. R. (2000). Ants: standard methods for measuring and monitoring bio-
diversity. Washington D.C.: Smithsonian Institution Press.

Andersen, A. (2000). A global ecology of rainforest ants: functional groups in relation to environmental stress and disturbance. In D. Agosti, J. D. Majer, L. E. Alonso, & T. R. Schultz (Eds.), Ants: standard methods for measuring and monitoring biodiversity (pp. 25–34). Washington D.C.: Smithsonian Institution Press.

Andersen, A.N. (1992). Regulation of “momentary” diversity by dominant species in exceptionally rich ant communities of the Australian seasonal tropics. The American Naturalist, 40, 401–420. https://www.journals.uchicago.edu/doi/abs/10.
1086/285419

AntWiki (n.d.). Nylanderia. AntWiki. Consulted 8/20/2024

Arnan, X., Cerdá, X., & Retana, J. (2012). Distinctive life traits and distribution along environmental gradients of dominant and subordinate Mediterranean ant species. Oecologia, 170, 489–500. https://doi.org/10.1007/s00442-012-2315-y

Castillo-Guevara, C., Cuautle, M., Lara, C., & Juárez-Juárez, B. (2019). Effect of agricultural land-use change on ant dominance hierarchy and food preferences in a temperate oak forest. PeerJ, 7, e6255. https://peerj.com/articles/6255/

Cerdá, X., Arnan, X., & Retana, J. (2013). Is competition a significant hallmark of ant (Hymenoptera: Formicidae) ecology? Myrmecological News, 18, 131–147. https://doi.org/10.25849/myrmecol.news_018:131

Cerdá, X., Retana, J., & Cros S. (1997). Thermal disruption of transitive hierarchies in Mediterranean ant communities. Journal of Animal Ecology, 66, 363–374. 

Costes-Quijano, R., Meza, A. R., Macías-Juárez, A., Berriel-Mastreta, C. A., Cortés-Atilano, B., Martínez-Romero, L. E. et al. (2006). Plan de manejo Parque Ecológico Recreativo General Lázaro Cárdenas “Flor del Bosque”. [Management Plan of the recreative ecological Park “General Lázaro Cárdenas “Flor del Bosque”]. Ciudad de México: Gobierno del Estado de Puebla/ Secretaría de Medio Ambiente y Recursos Naturales.

Cuautle, M., Vergara, C., & Badano, E. (2016). Comparison of ant community diversity and functional group composition associated to land use change in a seasonally dry oak forest. Neotropical Entomology, 45, 170–9. https://doi.org/
10.1155/2012/516058 

Cuezzo, F., & Guerrero, R. J. (2012). The Ant Genus Dorymyrmex Mayr (Hymenoptera: Formicidae: Dolicho-
derinae) in Colombia. Psyche, 51605, 1–24. https://doi.org/
10.1155/2012/516058

Dáttilo, W., Díaz-Castelazo, C., & Rico-Gray, V. (2014). Ant dominance hierarchy determines the nested pattern in ant-plant networks. Biological Journal of the Linnean Society, 113, 405–414. https://doi.org/10.1111/bij.12350.

Davidson, D. W. (2005). Ecological stoichiometry of ants in a New World rain forest. Oecologia, 142, 221–231. https://doi.org/10.1007/s00442-004-1722-0

Dejean, A., & Corbara, B. (2003). A review of mosaics of dominant ants in rainforests and plantations. In Y. Basset, V. Novotny, S. E. Miller, & R. L. Kitching (Eds.), Arthropods of tropical forests: spatio-temporal dynamics and resource use in the canopy (pp 341–347). Cambridge: Cambridge University Press.

Ellis, E. C., Goldewijk, K. K., Siebert S., Lightman, D., & Ramankutty, N. (2010). Anthropogenic transformation of
the biomes, 1700 to 2000. Global Ecology and Biogeo-
graphy, 19, 589–606. https://doi.org/10.1111/j.1466-8238.20
10.00540.x

Fellers, J. H. (1987). Interference and exploitation in a guild of Woodland ants. Ecology, 68, 1466–1478. https://doi.org/
10.2307/1939230

Fellers, J. H. (1989). Daily and seasonal activity in woodland
ants. Oecologia, 78, 69–76. https://doi.org/10.1007/BF0037
7199

Foley, J. A., DeFries, R., Asner, G. P., Barford C., Bonan, G., Carpenter, S. R. et al. (2005). Global consequences of land use. Science, 309, 570–574. https://doi/10.1126/science.1111772

Hoffmann, B. D., & Andersen, A. N. (2003). Responses of ants to disturbance in Australia, with particular reference to functional groups. Austral Ecologyl, 28, 444–464. –https://doi.org/10.1046/j.1442-9993.2003.01301.x

Houdria, M., Salas-López, A., Orivel. J., Bluthgen, N., & Menzel, F. (2015). Dietary and temporal niche differentiation in tropical ants – can they explain local ant coexistence? Biotropica, 47,208–217. https://doi.org/10.1111/btp.12184

Hernández-Flores, J., Osorio-Beristain. M., & Martínez-Garza. C. (2016). Ant foraging as an indicator of tropical dry forest restoration. Environmental Entomology, 45, 991–994. https://doi.org/10.1093/ee/nvw054

Kaspari, M., Donoso, D., Lucas, J. A., Zumbusch, T., & Kay, A. D. (2012). Using nutritional ecology to predict community structure: a field test in Neotropical ants. Ecosphere, 3,93. https://doi.org/10.1890/ES12-00136.1

Kaspari, M., & Yanoviak, S. P. (2001). Bait use in tropical litter and canopy ants-evidence of differences in nutrient limitation. Biotropica, 33,207–211. https://doi.org/10.1646/
0006-3606(2001)033[0207:BUITLA]2.0.CO;2

Łaska, G. (2001). The disturbance and vegetation dynamics: a review and an alternative framework. Plant Ecology, 157, 77–99. https://doi.org/10.1023/A:1013760320805

Lynch, J. F., Balinsky, E. C., & Vail, S. G. (1980). Foraging patterns in three sympatric forest ant species, Prenolepis imparis, Paratrechina melanderi and Aphaenogaster rudis (Hymenoptera: Formicidae). Ecological Entomology, 5,353–371. https://doi.org/10.1111/j.1365-2311.1980.tb01160.x

Mackay, W., & Mackay, E. (1989). Clave de los géneros de hormigas en México (Hymenoptera: Formicidae). El Paso, Texas: The University of Texas. 

Nettimi, R. P., & Iyer, P. (2015). Patch fidelity in Camponotus compressus ants foraging on honeydew secreted by treehoppers. Current Science, 109, 362–366. 

Parr, C. L. (2008). Dominant ants can control assemblage species richness in a South Africa savanna. Journal of Animal Ecology, 77, 1191–1198. https://doi.org/10.1111/j.13
65-2656.2008.01450.x 

Parr, L., & Gibb, H. (2012). The discovery-dominance trade-off
is the exception, rather than the rule Journal of Animal
Ecology, 81, 233–241. https://doi.org/10.1111/j.1365-2656.
2011.01899.x

Pirk, G. I., & López-de Casenave, J. (2014). Effect of harvester ants of the genus Pogonomyrmex on the soil seed bank around their nests in the central Monte desert, Argentina. Ecological Entomology, 39, 610–619. https://doi.org/10.1111/een.12140 

R Core Team (2022). R: a language and environment for statistical computing. Vienna, Austria: R Foundation
for Statistical Computing. https://www.R-project.org/

Radnan, G. N., Gibb, H., & Eldridge, D. J. (2018). Soil surface complexity has a larger effect on food exploitation by ants than a change from grassland to shrubland. Ecological Entomology, 43, 379–388. https://doi.org/10.1111/een.12510

Rico-Gray, V., & Oliveira, P. S. (2007). The ecology and evolution of ant-plant interactions. Chicago: University of Chicago Press. 

Sala, O. E., Chapin, F. S., Armesto, J. J., Berlow, E., Bloomfield, J., Dirzo. R. et al. (2000). Global biodiversity scenarios for the year 2100. Science, 287,1770–1774. https://doi.org/10.1126/science.287.5459.17

Santini, G., Tucci, L., Ottonetti, L., & Frizzi, L. (2007). Competition trade-offs in the organisation of a Medite-
rranean ant assemblage. Ecological Entomology, 32,
319–326. https://doi.org/10.1111/j.1365-2311.2007.00882.x 

Savolainen, R., Vepsäläinen, K., & Wuorenrinne, H. (1989). Ant assemblages in the taiga biome: testing the role of territorial wood ants. Oecologia, 81, 481–486. https://doi.org/10.1007/BF00378955

Louis C. Bender ^a, Jon C. Boren ^b, Shad Cox ^c, Erik Joaquín Torres-Romero ^{d, e,} *

^a New Mexico State University, Department of Extension Animal Sciences and Natural Resources, PO Box 30003 MSC 3AE, Las Cruces, New Mexico 88003, USA 

^b New Mexico State University, Cooperative Extension Service, PO Box 30003 MSC 3AE, Las Cruces, New Mexico 88003, USA

^c New Mexico State University, Corona Range and Livestock Research Center, PO Box 392, Corona, New Mexico 88318, USA

^d Universidad Politécnica de Puebla, Ingeniería en Biotecnología, Tercer Carril del Ejido, Serrano s/n, San Mateo Cuanalá, Juan C. Bonilla, 72640 Puebla, Mexico 

^e Tecnológico Nacional de México campus Zacapoaxtla, Subdirección de Investigación y Posgrado, División de Biología, Carretera Acuaco-Zacapoaxtla Km. 8, Col. Totoltepec, 73680 Zacapoaxtla, Puebla, Mexico

*Corresponding author: ejtr23@hotmail.com (E.J. Torres-Romero)

Received: 9 April 2024; accepted: 26 August 2024

Abstract

Deer in northern temperate environments show behavioral and physiological adaptations to conserve energy during winter, including decreased movements. Whether these behaviors persist in warmer temperate environments such as the arid Southwest has received little consideration. We compared daily movements as estimated by continuous-time movement models and minimum subdaily (4 h) straight-line movements of adult female mule deer between winter and spring-autumn seasons in south-central New Mexico. Deer moved less during winter daily (2.90 vs. 4.34 km/d) and subdaily (302 vs. 409 m). Similarly, for deer for which movement data for successive seasons were available, movements between successive seasons were less during the winter (daily = -1.05 km/d; subdaily = -91 m) than the following or preceding spring-autumn. Our results support conservation of decreased movements during winter in the less extreme winters of the arid Southwest. Because some proximate stimuli (i.e., deep snow, very cold temperatures) associated with energy conservation behaviors are lacking in the arid Southwest, our results further support low forage quality and availability being the primary drivers of this behavior.

Keywords: Energy conservation; Movements; Mule deer; New Mexico

© 2024 Universidad Nacional Autónoma de México, Instituto de Biología. Este es un artículo Open Access bajo la licencia CC BY-NC-ND

(http://creativecommons.org/licenses/by-nc-nd/4.0/).

Disminución de los movimientos de las hembras adultas de venado bura durante el invierno en el árido suroeste de América del Norte

Resumen

Los ciervos en entornos templados del norte muestran adaptaciones de comportamiento y fisiológicas para conservar energía durante el invierno, incluyendo una disminución en sus movimientos. Se ha explorado poco si estos comportamientos persisten en ambientes templados más cálidos, como el suroeste árido. Comparamos los movimientos diarios, mediante modelos de movimiento continuo y movimientos mínimos en línea recta subdiarios (4 horas) de hembras adultas de venado bura entre las estaciones de invierno y primavera-otoño en el centro-sur de Nuevo México. Los ciervos se movieron menos durante el invierno, tanto diario (2.90 vs. 4.34 km/día) como subdiario (302 vs. 409 m). Además, para ciervos con datos de movimiento en estaciones sucesivas, los movimientos en invierno fueron menores (diarios = -1.05 km/día, subdiarios = -91 m) en comparación con la primavera-otoño previa o siguiente. Nuestros resultados respaldan la disminución de los movimientos durante el invierno en los inviernos menos extremos del suroeste árido. Dado que algunos estímulos inmediatos (por ejemplo, nieve profunda, temperaturas muy frías) asociados con comportamientos de conservación de energía están ausentes en el suroeste árido, es evidente que nuestros resultados apoyan que la baja calidad y disponibilidad de forraje son los principales factores que impulsan este comportamiento.

Palabras clave: Conservación de energía; Movimientos; Venado bura; Nuevo México

Introduction

Deer in northern temperate environments of North America employ a complex energy conservation strategy during winter, incorporating multiple behavioral and physiological adaptations, including reducing activity (i.e., movements) and limiting feeding while relying on endogenous reserves (Alldredge et al., 1974; Short, 1981; Verme & Ullrey, 1984); growing a highly insulative pelage (Jacobsen, 1980); limiting vascular circulation to the extremities (Parker & Robbins, 1984); and lowering metabolic rate to slow the rate of loss of body reserves (Short, 1981; Silver et al., 1971; Verme & Ullrey, 1984). This strategy is considered primarily an adaptation to conserve energy in response to low forage availability and quality, as well as increased costs of movement associated with snow (Short, 1981; Verme & Ullrey, 1984). Conserving energy by minimizing radiant and convective heat loss was also believed to drive other behaviors such as yarding under dense conifer forest canopies (Marchinton & Hirth, 1984). While initially thought to be a response to extreme cold (Marchinton & Hirth, 1984), yarding likely relates to decreased costs of movement because of reduced snow depths as the presumed temperature-moderating influence of forest canopy (i.e., thermal cover) has been shown to have no real effect on deer condition (Cook et al., 1998; Freddy, 1984). Many of the behavioral aspects (at least) of the winter energy conservation strategy are not invariant, however, and can be affected by proximate stimuli. For example, both movements and feeding periods are reduced less if winter conditions are less severe (Bartmann & Bowden, 1984; Verme & Ullrey, 1984). 

Whether these energy conservation behaviors persist in warmer temperate ranges such as the arid southwestern USA and Mexico has received little consideration. In the arid Southwest, deer similarly experience winter seasonality in terms of temperature and precipitation differences (Krausman et al., 1990; McKinney, 2003; Marshal et al., 2008), which affects forage availability and quality (Kemp, 1983; Krausman et al., 1990; McKinney, 2003; Short, 1981). Consequently, mule deer lose most of their endogenous reserves over winter (Bender et al., 2012; Bender & Hoenes, 2017). However, snow is relatively rare and short-lived in much of the arid Southwest, and winter temperatures are higher than in northern environments (Table 1). Therefore, aside from decreased forage availability and quality, many of the potential proximate stimuli (i.e., deep snow, very cold temperatures) associated with energy conservation behaviors are lacking in the arid Southwest. Moreover, because winter is less extreme in terms of minimum temperatures and particularly snowfall, availability of forage may also be less limiting, although forage quality constraints are similarly severe (Bender, 2020; Kemp, 1983; Krausman et al., 1990; McKinney et al., 2003). Consequently, behavioral responses associated with the winter energy conservation strategy may be less pronounced or absent in the arid Southwest.

Table 1

Long-term range of monthly mean high and low temperatures (^oC) and monthly snowfall (cm) during Dec.-Feb. and Mar.-Nov. at the Corona Range and Livestock Research Center (CRLRC), Corona, New Mexico USA, and the Cusino Wildlife Research Station, Shingleton, Michigan USA. Cusino was selected as a comparison because of the volume and depth of deer nutritional, physiological, and behavioral research conducted there (Verme & Ullrey, 1984). Also presented are the range of conditions on the CRLRC study area for 2005-2008 study period.

Months	Climatic variable	CRLRC	Cusino	Study
Dec.-Feb.	Mean high temperatures	6.7-8.9	-2.8- -1.7	6.4-13.0
	Mean low temperatures	-5.6- -4.4	-11.7- -8.3	-6.3- -1.9
	Mean snow accumulations	13-23	69-109	7.3-17.1
Mar.-Nov.	Mean high temperatures	11.7-28.3	1.7-22.8	14.2-29.1
	Mean low temperatures	-1.1-13.3	-6.7-13.3	-0.5-13.8
	Mean snow accumulations	0-13	0-48	0-5.1

If energy conservation behaviors are maintained in the arid Southwest, mule deer (Odocoileus hemionus) should move less during the winter, conserving body reserves in the face of lower quality and less abundant forage even if the impacts of winter weather are less severe on forage availability and costs of movement. Thus, our goal was to contrast short-term movements between winter and spring-autumn seasons for adult female mule deer in a Chihuahuan desert-short grass prairie habitat in New Mexico, USA to determine whether deer reduce movements during winter as predicted by the winter energy conservation strategy. Specifically, we compared minimum daily and subdaily movement distances of adult female mule deer between winter and spring-autumn. 

Materials and methods

Our study was conducted on the Corona Range and Livestock Research Center (CRLRC; 34°15’36” N, 105°24’36” W), an 11,290-ha ranch owned and operated by New Mexico State University and located approximately 22.5 km east of Corona, New Mexico (Fig. 1). CRLRC has an average elevation of 1,900 m asl; mean annual precipitation is 40 cm, 87% of which occurs in the Mar.-Nov. period. Snowfall totals < 74 cm annually. Climate of the CRLRC shows distinct seasonality, although the magnitude of seasonal differences in winter is less than seen in northern temperate deer habitats (Table 1). 

Topography of the CRLRC is mostly rolling. Vegetation includes perennial grassland, with scattered sparse to dense pinyon (Pinus edulis) and one-seed juniper (Juniperus monosperma) woodlands and a few shrublands. Free water was abundant and comparably available in both winter and spring-autumn seasons because of numerous permanent water developments, ≥ 1 of which were present within or adjacent to annual home ranges of study deer (Fig. 1). Deer on the CRLRC do not migrate between distinct summer and winter ranges. 

We captured and collared ≥ 2.5-year-old female mule deer with GPS/VHF radio-collars (Advanced Telemetry Solution, Asanti, Minnesota, USA) programmed to record a position fix every 4 h, early-December 2005-2007, and April, 2006-2007, as part of a larger study of mule deer ecology including other VHF-only radio-collared individuals (Bender et al., 2011, 2013). Deer were captured using a helicopter by aerial net-gunning or darting with 1.5-1.8 mg of carfentanil citrate and 50-75 mg of xylazine hydrochloride per deer. We aged deer as yearling or adult by tooth wear and replacement (Robinette et al., 1957), determined lactation status (Bender et al., 2011), and treated deer with antibiotics, vitamin E/selenium, vitamin B, and an 8-way Clostridium bacterin to help alleviate capture stress. Following processing, immobilants were antagonized with naltrexone and tolazoline.

We defined seasons as winter = Dec-Feb and spring-autumn = Mar-Nov. These seasons corresponded with both typical seasonal and phenological patterns on the CRLRC and the Chihuahuan desert-short grass prairie habitats of the arid Southwest in general, as well as important periods in the annual cycle of female mule deer in the arid Southwest (Bender et al., 2011, 2012). We estimated daily movements using continuous-time movement modeling (see Fleming et al., 2014, 2016), using the Speed/Distance analysis in ctmmweb (https://ctmm.shinyapps.io/ctmmweb/) (Calabrese et al., 2016, 2021). We used only locations with 3D fixes and DOP < 2, as these had an accuracy of < 3 m in our study area. We also determined minimum subdaily movements, defined as the straight-line distance moved between successive 4 h locations, and calculated seasonal means for each deer. 

Figure 1. Topographic hillshade showing locations of annual home ranges of adult female mule deer (Bender et al., 2013), and locations of permanent water sources on the Corona Range and Livestock Research Center (CRLRC), east-central New
Mexico, USA.

We compared movement distances (km from continuous-time movement models [ctmms]; m between successive subdaily locations) between seasons using PROC GLIMMIX in SAS 9.4 (SAS, 1988), using individual deer as a random effect. We also compared mean daily and subdaily movements between successive seasons for individual deer for which we had both winter and subsequent spring-autumn data, or spring-autumn and subsequent winter data, available. We determined mean seasonal differences in movement distances for each successive time period for each individual, and used bootstrapping with N = 1,000 iterations to determine the probability that mean movement distances differed seasonally for this subset of data (see Efron & Tibshirani, 1993).

Additionally, because lactating females enter winter in poorer condition than do dry females and condition subsequently converges between the 2 classes over winter (including on the CRLRC; Bender & Hoenes, 2017), lactation status might affect desire or need to forage and thus movements of deer. However, because deer condition was very low on CRLRC during our study (i.e., lactating females were able to accrue only ≤ 5.7% percent body fat annually at the annual peak in late autumn; Bender et al., 2011, 2013), our GPS/VHF collared sample contained ≤ 2 lactating females each year, and thus we were unable to meaningfully include lactation status in our analyses. Consequently, we explored any potential effect of lactation status on movements by determining the percentile movement distances of lactating females relative to the frequency distribution of dry female movement distances to determine whether lactating females were closer to the mean or extremes of the range of dry females.

Lastly, because our GPS collared sample comprised a limited proportion of the overall radio-collared sample, we compared annual, spring-autumn, and winter home range sizes (Bender et al. 2013; N = 18-27 collared females annually) to see whether movements of our GPS sample were representative of radio-collared deer in general, comprised of N = 18-27 VHF-only collared adult females for each season (Bender et al., 2011, 2013). For this we compared annual and seasonal home range sizes of the 2 classes (GPS/VHF and VHF-only) using PROC GLM (SAS, 1988), specifically testing the year × class interaction. We visually located all deer (i.e., both GPS/VHF and VHF-only) via ground tracking of VHF signals a minimum of once per week with additional location emphasis on spring-autumn locations, and mapped locations using the Geographic Information System software package ArcGIS 10.0 (Environmental Systems Research Institute, Redlands, California, USA) (Bender et al., 2011, 2013). We constructed 95% minimum convex polygon (MCP) annual and seasonal home ranges after determining the minimum number of locations to adequately estimate seasonal home range size by plotting size as a function of number of locations (Bender et al. 2013; Kie et al. 1996). For this comparison we used only VHF visual locations of GPS/VHF collared deer so that both GPS/VHF and VHF-only samples were comprised of comparable data. 

Results

We collected GPS movement data for 37 seasonal ranges (Table 2), winter 2005-6 through winter 2007-8, from 6-10 GPS/VHF collared adult females annually (mean = 336 and 472 locations per deer for winter and spring-autumn ranges, respectively). For all deer, the OUF-anisotropic movement model provided the best fit (i.e., AICc < 2 vs. all other models) of deer movements. The OUF-anisotropic movement model is the most general of ctmms, and includes a home range, correlated positions, correlated velocities, and movements varying by direction (Calabrese et al., 2021; Fleming et al., 2014).

Table 2

Mean distance moved by adult female mule deer during winter (Dec.-Feb.) and spring-autumn (Mar.-Nov.) seasons as estimated by continuous-time movement modeling (daily) and subdaily straight line movements between successive 4 h locations (subdaily) on the Corona Range and Livestock Research Center, east-central New Mexico, 2005-2007. 

	Daily (km)	Subdaily (m)
Season	Distance	SE	N	Distance	SE	N
Winter	2.90	0.24	24	302	19	24
Spring-Autumn	4.34	0.43	13	409	44	13

Deer moved less during the winter (Table 2) for both daily ctmms (F_1,20 = 6.0; p = 0.024) and minimum subdaily straight-line distances (F_1,20 = 4.8; p = 0.041); in both cases the magnitude of differences varied among individual deer (t₁₅ < -2.72; p < 0.016). Similarly, for deer for which movement data for successive seasons were available, movements between successive seasons were always less during the winter season (p (winter < spring-autumn) = 1.000) than the following or preceding spring-autumn season (Table 3). 

Additionally, movement distances of lactating females (mean = 42^nd percentile; range = 34-55^th percentile) were always closer to the average than the extremes of the frequency distribution of movement distances of dry females for each period. Last, neither annual, spring-autumn, or winter home MCP range sizes (N = 7-16 per period) differed between GPS/VHF and VHF-only collared deer (F_5,39 ≤ 0.94; p ≥ 0.599). 

Discussion

Despite much less snow and warmer winter temperatures, mule deer females moved less during winter regardless of movement period (i.e., daily, subdaily), supporting the maintenance of this energy conservation behavior in mule deer in the arid Southwest (Verme & Ullrey, 1984). Because deep snow cover and very cold temperatures are typically lacking in the arid Southwest, our results further support low forage quality and availability being the primary drivers of this behavior in deer (Verme & Ullrey, 1984; see below). This latter was reflected by the poor nutritional condition of lactating adult females in the study area (i.e., <5.7% body fat annually at the annual peak in late autumn; Bender et al., 2011, 2013), and in the arid Southwest in general (Bender et al., 2007, 2011, 2012; Bender, 2020), as well as the significant losses in condition seen over winter (Bender & Hoenes, 2017; Hoenes, 2008), despite mule deer likely requiring lower quality forage than white-tailed deer (O. virginianus) (Staudenmaier et al., 2022).

While lack of deep snow cover may result in relative forage availability being less impacted during winter in the arid Southwest, senesced forages are still of very low quality, similar to northern environments (Kemp, 1983; Krausman et al., 1990; McKinney, 2003). Low forage quality in winter was reflected in the condition dynamics of deer on the CRLRC; for example, dry females lost on average 32% of body fat reserves and 38% of rump body condition score over winter (L. Bender, unpublished data), even though their condition in late autumn-early winter was already low relative to other southwestern populations (Bender et al., 2007, 2011, 2012). The nutritional condition of deer is driven primarily by forage quality (Bender, 2020; National Research Council, 2007; Tollefson et al., 2010; Verme & Ullrey, 1984; Wakeling & Bender, 2003), illustrating that deer in the arid Southwest face similar constraints in terms of limited nutrient gains from forage intake as do deer in northern environments during winter.

Table 3

Mean movement distances (x) and mean differences (D) in mean distance moved daily as predicted by continuous-time movement modeling (daily) and subdaily straight line movements between successive 4 h locations (subdaily) of adult female mule deer for which movement data for successive seasons were available between winter (Dec.-Feb.) and spring-autumn (Mar.-Nov.) on the Corona Range and Livestock Research Center, east-central New Mexico, 2005-2007. p = Probability that seasonal differences differ from 0; N = number of seasonal comparisons. 

Period	x Summer	x Winter	D	90% CI	p	N
Subdaily	400.2	309.6	-90.6 m	-123- -63	1.000	20
Daily	4.21	3.16	-1.05 km	-1.53- -0.63	1.000	20

Deer in northern environments do face additional energetic challenges associated with persistent snow cover, which can limit forage availability (by making location and acquisition of food more difficult; Hovey & Harstad, 1992), diet quality (due to reduced forage quality and availability; McKinney, 2003; Osborn & Jenks, 1998), and increase loss of endogenous reserves (because of increased costs of moving through snow; Bunnell et al., 1990; Mattfeld, 1973). However, while deer in the arid Southwest are less influenced by snow-depth related challenges, northern deer do not face the lack of free water experienced by most deer populations in arid environments because of persistent snow cover in northern environments. Lack of water can present an energetic cost to deer in the arid Southwest, as mule deer may increase movement distances to access water (Heffelfinger, 2006), and winter is much drier than spring-autumn in Chihuahuan desert and short-grass prairie habitats of the arid Southwest (e.g., 87% of precipitation occurs during spring-autumn on the CRLRC). 

Because of permanent water developments, water was comparably available seasonally on the CRLRC; water developments were accessible from all deer home ranges, so mule deer did not need to alter their movements in response to seasonal changes in availability of water. Hence, need or preference for free water likely had a negligible effect on deer movements on the CRLRC, unless presence of temporary sources (ephemeral pools, etc.) during the summer monsoon reduce deer movements during spring-autumn because of increased availability. Thus, despite water being effectively controlled in our study, deer still showed less movements during winter. This again supports decreased movements during winter being most influenced by the lack of energetic benefit from seeking and foraging on low quality senesced forage.

Cold temperatures are often thought to influence the winter conservation strategy, despite demonstrated lack of benefit of thermal cover in winter to deer (Cook et al., 1998; Freddy, 1984), the pronounced effect of solar radiation on warming deer (Cook et al., 1998; Parker & Gillingham, 1990; Parker & Robbins, 1984), and deer movements and tolerance of exposure to cold (including bedding in the open) except during the most extreme conditions when high quality forage is available (Moen, 1968; Verme & Ullrey, 1984). Mule deer possess a low thermal critical zone (ca. -20 ^oC) and show greater tolerance of cold than do white-tailed deer, including a lower metabolic rate response to decreasing temperatures (Mautz et al., 1985; Parker & Robbins, 1984). Hence, they are less affected by even extreme cold, which is seldom the case in the arid Southwest where mean low temperatures seldom approach their lower thermal tolerance (Table 1). Conversely, the upper thermal critical level of mule deer in winter pelage (5 ºC; Mautz et al., 1985) is lower than average high temperatures during winter on our study area (6.7-8.9 ºC; Table 1) and much of the arid Southwest. This would require active metabolic activity or behaviors (e.g., panting, etc.) by mule deer to cool themselves, or possibly limiting movements and bedding under shade (although the energetic benefits of the latter are questionable; Cook et al., 1998). Consequently, if temperature affected movements during the winter in our study area, deer would be more likely to decrease movements because of heat stress, not cold stress. 

Lastly, while at least one behavioral aspect of the winter energy conservation strategy is seen in mule deer in the arid Southwest, to what degree other adaptations are conserved is unknown. Mule deer in the arid Southwest do develop a highly insulative pelage in winter (Heffelfinger, 2006), but the extent that they may decrease metabolic rate or regulate vascular circulation to the extremities (or need to, in light of the more moderate temperatures) is unknown. Additionally, while lactating females showed a tendency to move less relative to dry females on the CRLRC, whether this is typical is unknown because of our small sample of lactating females. However, given that most females on CRLRC were in poorer condition than females elsewhere in the arid Southwest (Bender et al., 2007, 2011, 2012; Bender & Hoenes, 2017), if very low condition (such as results from lactation) increases movements during winter, this increase would likely have been seen in all CRLRC females regardless of lactation status. Moreover, GPS/VHF and VHF-only collared females showed similar movements (as indexed by home range sizes) seasonally and annually on CRLRC, indicating that movements of GPS/VHF collared females reflected females in general. Although small sample sizes (3-4) precluded including lactation status as an interactive term in the contrast of GPS/VHF and VHF-only collared females, for all females lactation status had no effect on annual or seasonal home range sizes (p ≥ 0.599).

Acknowledgments

Support for this project was provided by the U.S. Forest Service-Rocky Mountain Research Station and the New Mexico State University Cooperative Extension Service and Agricultural Experimental Station. All activities were in accordance with NMSU IACUC Permit No. 2005-023. E.J.T.-R. was supported by a postdoctoral fellowship from Consejo Nacional de Humanidades, Ciencias y Tecnologías (Conahcyt-Mexico).

References

Alldredge, A. W., Lipscomb, J.F., & Whicker, F.W. (1974). Forage intake rates of mule deer estimated with fallout
cesium-137. The Journal of Wildlife Management, 38, 508–516. https://doi.org/10.2307/3800882

Bartmann, R. M., & Bowden, D. C. (1984). Predicting mule deer mortality from weather data in Colorado. Wildlife Society Bulletin, 12, 246–248. https://doi.org/10.1002/jwmg.956

Bender, L. C. (2020). Elk, deer, and pinyon-juniper: needs, what works, and what doesn’t. In K. Malcolm, B. Dykstra, K. Johnson, D. Lightfoot, E. Muldavin, & M. Ramsey (Compilers), Symposium proceedings on piñon-juniper habitats: status and management for wildlife – 2016. Proceedings RMRS-P-77. Fort Collins, Colorado: U.S. Forest Service, Rocky Mountain Research Station.

Bender, L. C., & Hoenes, B. D. (2017). Costs of lactation to body condition and future reproduction of free-ranging mule deer Odocoileus hemionus (Cervidae). Mammalia, 81, 329–338. https://doi.org/10.1515/mammalia-2015-0143

Bender, L. C., Lomas, L. A., & Browning, J. (2007). Condition, survival, and cause-specific mortality of mule deer in northcentral New Mexico. Journal of Wildlife Management, 71, 1118–1124. https://doi.org/10.2193/2006-226

Bender, L. C., Boren, J. C., Halbritter, H., & Cox, S. (2011). Condition, survival, and productivity of mule deer in semiarid grassland-woodland in east-central New Mexico. Human-Wildlife Interactions, 5, 276–286. https://doi.org/10.26077/60n0-ks48

Bender, L. C., Hoenes, B. D., & Rodden, C. L. (2012). Factors influencing survival of desert mule deer in the greater San Andres Mountains, New Mexico. Human-Wildlife Interactions, 6, 245–260. https://doi.org/10.26077/h5bg-1829

Bender, L. C., Boren, J. C., Halbritter, H., & Cox, S. (2013). Effects of site characteristics, pinyon-juniper control, and precipitation on habitat quality for mule deer on the Corona Range and Livestock Research Center. Human-Wildlife Interactions, 7, 47–59. https://digitalcommons.usu.edu/hwi/vol7/iss1/5

Bunnell, F. L., Hovey, F. W., McNay, R. S., & Parker, K. L. (1990). Forest cover, snow conditions, and black-tailed deer sinking depths. Canadian Journal of Zoology, 68, 2403–2408. https://doi.org/10.1139/z90-333

Calabrese, J. M., Fleming, C. H., & Gurarie, E. (2016). ctmm: an R package for analyzing animal relocation data as a continuous-time stochastic process. Methods in Ecology and Evolution, 7, 1124–1132. https://doi.org/10.1111/2041-210X.12559

Calabrese, J. M., Fleming, C. H., Noonan, M. J., & Dong, X. (2021). ctmmweb: a graphical user interface for autocorrelation-informed home range estimation. Wildlife Society Bulletin, 45, 162–169. https://doi.org/10.1002/wsb.1154

Cook, J. G., Irwin, L. L., Bryant, L. D., Riggs, R. A., & Thomas, J. W. (1998). Relations of forest cover and condition of elk: a test of the thermal cover hypothesis in summer and winter. Wildlife Monographs, 141, 1–61.

Efron, B., & Tibshirani, R. J. (1993). An introduction to the bootstrap. New York: Chapman & Hall.

Fleming, C. H., Fagan, W. F., Mueller, T., Olson, K. A., Leimgruber, P., & Calabrese, J. M. (2014). From fine-scale foraging to home ranges: a semivariance approach to identifying movement modes across spatiotemporal scales. American Naturalist, 183, E154–E167. http://doi.org/10.1086/675504 

Fleming, C. H., Fagan, W. F., Mueller, T., Olson, K. A., Leimgruber, P., & Calabrese, J. M. (2016). Estimating where and how animals travel: an optimal framework for path reconstruction from autocorrelated tracking data. Ecology, 97, 576–582. https://doi.org/10.1890/15-1607.1

Freddy, D. J. (1984). Quantifying capacity of winter ranges to support deer – evaluation of thermal cover used by deer. Denver, CO: Wildlife Research Report, Colorado Division of Wildlife, USA.

Heffelfinger, J. (2006). Deer of the Southwest: a complete guide to the natural history, biology, and management of Southwestern mule deer and white-tailed deer. Texas A&M University, College Station, Texas.

Hoenes, B. D. (2008). Identification of factors limiting desert mule deer populations in the greater San Andres Mountains of southcentral New Mexico (Thesis). New Mexico State University, Las Cruces. 

Hovey, F. W., & Harestad, A. S. (1992). Estimating effects of snow on shrub availability for black-tailed deer in southwestern British Columbia. Wildlife Society Bulletin, 20, 308–313.

Jacobsen, N. K. (1980). Differences of thermal properties of white-tailed deer pelage between seasons and body regions. Journal of Thermal Biology, 5, 151–158. https://doi.org/10.1016/0306-4565(80)90014-5

Kemp, P. R. (1983). Phenological patterns of Chihuahuan desert plants in relation to the timing of water availability. Journal of Ecology, 71, 427–436.

Kie, J. G., Baldwin, J. A., & Evans, C. J. (1996). CALHOME: a program for estimating animal home ranges. Wildlife Society Bulletin, 24, 342–344.

Krausman, P. R., Ordway, L. L., Whiting, F. M., & Brown, W. H. (1990). Nutritional composition of desert mule deer forage in the Picacho Mountains, Arizona. Desert Plants, 10, 32–34.

Marchinton, R. L., & Hirth, D. H. (1984). Behavior. In L. K. Halls (Ed.), White-tailed deer ecology and management. Pennsylvania: Stackpole Books.

Marshal, J. P., Krausman, P. R., & Bleich, V. C. (2008). Body condition of mule deer in the Sonoran Desert is related to rainfall. Southwestern Naturalist, 53, 311–318. https://doi.org/10.1894/CJ-143.1

Mattfeld, G. F. (1973). The effect of snow on the energy expenditure of walking white-tailed deer. In Transactions of the 30th Northeast Fish and Wildlife Conference, Dover, Vermont, Spring, April 1973.

Mautz, W. W., Pekins, P. J., & Warren, J. A. (1985). Cold temperature effects on metabolic rates of white-tailed, mule, and black-tailed deer in winter coat. In P. F. Fennessy, & K. R. Orew (Eds.), The biology of deer. Royal Society of New Zealand Bulletin, 22, 453–457. 

McKinney, T. (2003). Precipitation, weather, and mule deer. In J. C. Jr. deVos, M. R. Conover, & N. E. Headrick (Eds.), Mule deer conservation. Issues and management strategies. Logan, Utah: Jack H. Berryman Press. 

Moen, A. N. (1968). Energy exchange of white-tailed deer, western Minnesota. Ecology, 49, 676–682. https://doi.org/10.2307/1935531

National Research Council (US). (2007). Committee on Nutrient Requirements of Small Ruminants. Nutrient requirements of small ruminants: sheep, goats, cervids, and New World camelids. Washington D.C.: National Academy Press. 

Osborn, R. G., & Jenks, J. A. (1998). Assessing dietary quality of white-tailed deer using fecal indices: effects of supplemental feeding and area. Journal of Mammalogy, 79, 437–447. https://doi.org/10.2307/1382974

Parker, K. L., & Robbins, C. T. (1984). Thermoregulation in mule deer and elk. Canadian Journal of Zoology, 62, 1409–1422. https://doi.org/10.1139/z84-202

Parker, K. L., & Gillingham, M. P. (1990). Estimates of critical thermal environments for mule deer. Journal of Range Management, 43, 73–81.  http://doi.org/10.2307/3899126  

Robinette, W. L., Jones, D. A., Rogers, G., & Gashwiler, J. S. (1957). Notes on tooth development and wear for Rocky Mountain mule deer. Journal of Wildlife Management, 21, 134–153. https://doi.org/10.2307/3797579

SAS. (1988). SAS/STAT user’s guide. SAS Institute Incorporated, Cary, North Carolina.

Short, H. L. (1981). Nutrition and metabolism. In O. C. Wallmo, R. E. McCabe, & L. R. Jahn (Eds.), Mule and black-tailed deer of North America. Lincoln: University of Nebraska Press.

Silver, H., Holter, J. B., Colovox, N. F., & Hayes, H. H. (1971). Effect of falling temperature on heat production in fasting white-tailed deer. Journal of Wildlife Management, 35, 37–46. 

Staudenmaier, A. R., Shipley, L. A., Camp, M. J., Forbey, J. S., Hagerman, A. E., Brandt, A. E. et al. (2022). Mule deer do more with less: comparing their nutritional requirements and tolerances with white-tailed deer. Journal of Mammalogy, 103, 178–195. https://doi.org/10.1093/jmammal/gyab116

Tollefson, T. N., Shipley, L. A., Myers, W. L., Keisler, D. H., & Dasgupta, N. (2010). Influence of summer-autumn nutrition on body condition and reproduction in lactating mule deer. Journal of Wildlife Management, 74, 974–986. https://doi.org/10.2193/2008-529

Verme, L. J., & Ullrey, D. E. (1984). Physiology and nutrition. In L. K. Halls (Ed.), White-tailed deer: Ecology and management. Harrisburg, Pennsylvania: Stackpole Books.

Wakeling, B. F., & Bender, L. C. (2003). Influence of nutrition on mule deer biology and ecology. In J. C. Jr. deVos, M. R. Conover, & N. E. Headrick (Eds.), Mule deer conservation. Issues and management strategies. Logan, Utah: Jack H. Berryman Press.

Fernando Arriola-Álvarez ^a, Luis Gabriel Aguilar-Estrada ^a, *, Lucía Álvarez-Castillo ^b,
Ivette Ruiz-Boijseauneau ^a y Dení Rodríguez ^a

^a Universidad Nacional Autónoma de México, Facultad de Ciencias, Laboratorio de Ficología (Biodiversidad Marina), Circuito Exterior s./n., Coyoacán, 04510 Ciudad de México, México

^b Universidad Nacional Autónoma de México, Facultad de Ciencias, Posgrado en Ciencias del Mar y Limnología, Ciudad Universitaria, Coyoacán, 04510 Ciudad de México, México

*Autor para correspondencia: lgae@ciencias.unam.mx (L.G. Aguilar-Estrada)

Recibido: 7 julio 2022; aceptado: 11 septiembre 2024

Resumen

Las macroalgas intermareales proporcionan alimento y refugio para diferentes organismos. El objetivo de este trabajo fue analizar los bivalvos asociados a macroalgas. Se realizaron muestreos en enero, mayo, julio y noviembre de 2014, se recolectaron manualmente 72 muestras de macroalgas y sus bivalvos asociados dentro de cuadros de 400 cm² en 2 localidades de Guerrero: playas El Palmar y Las Gatas. La estructura comunitaria de los bivalvos se determinó a partir de la riqueza específica, composición, abundancia, distribución e índices comunitarios: diversidad de Shannon, equidad de Pielou y dominancia de Simpson. Cada especie de macroalga (59 spp.) se asoció con la propuesta de grupos morfofuncionales. Se analizó la cobertura de macroalgas, abundancia de bivalvos y sedimento retenido. Del total de individuos (873), se reconocieron 17 especies de bivalvos. El índice de Shannon fue de 2.15 bits/individuo. Los bivalvos se asociaron a 3 grupos morfofuncionales de macroalgas. La abundancia de bivalvos y los sedimentos retenidos disminuyeron por mes, mientras que la abundancia, cobertura y sedimentos disminuyeron al aumentar el nivel de marea. Estudios como este proporcionan información importante para el conocimiento de la diversidad costera, en este caso de una zona turística en Guerrero.

Palabras clave: Malacofauna; Pacífico tropical mexicano; Riqueza específica; Sedimentos

© 2024 Universidad Nacional Autónoma de México, Instituto de Biología. Este es un artículo Open Access bajo la licencia CC BY-NC-ND

(http://creativecommons.org/licenses/by-nc-nd/4.0/).

Community structure of bivalves (Mollusca: Bivalvia) associated with intertidal macroalgae of Guerrero, Mexico

Abstract

Intertidal macroalgae provide food and shelter for different organisms. The objective of this work was to analyze the bivalves associated with macroalgae. Sampling was carried out in January, May, July, and November 2014, 72 samples of macroalgae and their associated bivalves were manually collected within 400 cm² in 2 locations in Guerrero: El Palmar and Las Gatas beaches. The community structure of bivalves was determined from specific richness, composition, abundance, distribution, and community indices: Shannon diversity, Pielou evenness and Simpson dominance. Each macroalgal species (59 spp.) was associated with proposed morphofunctional groups. Macroalgal cover, bivalve abundance and retained sediment were analyzed. Of the total number of individuals (873), 17 bivalve species were recognized. The Shannon index was 2.15 bits/individual. Bivalves were associated with 3 morphofunctional groups of macroalgae. Bivalve abundance and retained sediment decreased by month, while abundance, cover, and sediment decreased with increasing tide level. Studies like this provide important information for understanding coastal diversity, in this case of a tourist area in Guerrero.

Keywords: Malacofauna; Tropical Mexican Pacific; Species richness; Sediments

Introducción

En el Pacífico tropical mexicano se han descrito 2 patrones principales de circulación de corrientes oceánicas: primavera (marzo-abril) y otoño (septiembre-octubre), que generan variaciones espacio-temporales como las temporadas de lluvias y secas o fenómenos climatológicos (Baumgartner y Christensen, 1985; Pérez, 2013; Vega et al., 2008; Wyrtki, 1966). La zona intermareal rocosa es un sitio de transición entre los ambientes terrestre y marino, la cual se encuentra sujeta a cambios constantes de las variables ambientales como la oscilación de la marea, la intensidad lumínica, el viento, las variaciones en la salinidad y la temperatura (Salazar-Vallejo y González, 1990; Vassallo et al.,2014); ésto genera un hábitat heterogéneo con diversos microambientes en donde varios grupos de organismos pueden desarrollarse (Flores-Garza etal., 2011, 2014). Los organismos más frecuentes en esta zona son las macroalgas (Lee, 2008), mismas que proporcionan refugio y alimento para numerosos grupos de invertebrados (García-Robledo et al.,2008; Jover-Capote y Diez, 2017; Moreno, 1995; Steneck y Watling, 1982; Yang et al., 2007). Las macroalgas son un ambiente espacialmente heterogéneo, lo que hace posible que puedan albergar distintos grupos de invertebrados a lo largo del tiempo (Benedetti-Cecchi et al., 2001; Olabarria y Chapman, 2001). Los anfípodos, poliquetos y moluscos son los grupos más importantes al interior de la comunidad de macroalgas, ya que representan 70% de la abundancia en éstas (Aguilera, 2011; Colman, 1940).

Dentro del phylum Mollusca, los bivalvos son la segunda clase más representativa (Gosling, 2015). En la zona intermareal pueden vivir adheridos a diversos sustratos como rocas, arena o macroalgas (García-Cubas y Reguero, 2007). Algunos organismos que conforman la clase Bivalvia son sésiles y tienen diferentes estrategias en cuanto a sus tipos de alimentación: suspensívoros o detritívoros (Coan y Valentich-Scott, 2006). Además, generan redes de mucus para atrapar las partículas que flotan en la columna de agua (Jorgensen, 1996; Ward et al., 1998), por lo que desempeñan papeles ecológicos importantes en los cuerpos de agua; por sus hábitos de vida, son un grupo de especial interés en los estudios ecológicos (Lozada, 2015; Vega et al., 2008).

En general, los trabajos sobre moluscos en México son numerosos. En el de Sánchez (2014) se mencionó que existen alrededor de 47 contribuciones tomando en cuenta las costas del Atlántico y del Pacífico. En Guerrero se cuenta con alrededor de 30 estudios malacológicos (Gama, 2019), los cuales en su mayoría se han orientado a conocer la riqueza y composición de especies (Flores, 2004; Flores-Rodríguez et al., 2012; Lesser, 1984; López-Rojas et al., 2017). Gran parte (70%) de los trabajos en este estado se han realizado en el área de Acapulco y se han analizado distintos aspectos de la comunidad de moluscos, incluyendo a los bivalvos (Barba-Marino et al., 2010; Castro-Mondragón et al.,2016; Flores-Garza et al., 2010, 2011, 2012, 2014; Flores-Rodríguez et al.,2003; Galeana-Rebolledo et al., 2012, 2018; Garcés, 2011; Kuk-Dzul et al.,2019; Torreblanca, 2010; Torreblanca-Ramírez et al., 2012; Valdés-González etal.,2004; Villegas-Maldonado et al., 2007; Villalpando, 1986).

En la parte norte de Guerrero, que incluye Ixtapa-Zihuatanejo, se han realizado estudios ecológicos y sobre ciclos reproductivos de moluscos (Baqueiro, 1979; Flores-Rodríguez et al.,2007; Salcedo-Martínez et al.,1988), o sobre especies de importancia comercial como el de Cerros-Cornelio et al. (2021), quienes mencionaron 24 especies de moluscos, de las cuales 13 son de la clase Bivalvia. Por su parte, en la costa sur de Guerrero existen 9 trabajos malacológicos; Flores-Garza et al. (2007) analizaron la densidad de Plicopurpura columellaris (Lamarck, 1816) y su malacofauna asociada, reportando 34 especies de moluscos, de las cuales 7 fueron bivalvos.

Entre los estudios realizados en las costas de Guerrero, resalta el de Salcedo-Martínez et al. (1988) por ser el primero sobre la relación entre las macroalgas e invertebrados de Zihuatanejo, donde los moluscos, en especial la clase Gastropoda, fueron el componente mayoritario (38.72%). Existen algunos trabajos sobre la asociación alga-molusco en Ixtapa-Zihuatanejo (Aguilar-Estrada et al., 2017, 2022; Cisneros, 2016; Gama-Kwick et al., 2021; Quiroz-González et al., 2020); sin embargo, dichos estudios están enfocados a otras clases de moluscos (gasterópodos y quitones); por lo que, el conocimiento sobre la relación de las macroalgas y bivalvos es escaso. La presente contribución tiene como objetivo aportar conocimiento de la estructura comunitaria de los moluscos bivalvos asociados a macroalgas, en un ciclo anual en la zona intermareal rocosa de Ixtapa-Zihuatanejo en Guerrero.

Materiales y métodos

Se realizaron 4 salidas de campo durante enero, mayo, julio y noviembre de 2014 a Ixtapa-Zihuatanejo, Guerrero, con el propósito de observar los posibles cambios de la estructura comunitaria de los bivalvos asociados a las macroalgas en la zona. Las comunidades de macroalgas se recolectaron en la zona intermareal rocosa en 2 localidades: playa El Palmar en Ixtapa (17°39’0.4” N, 101°36’2.79” O) y el pretil de playa Las Gatas, al interior de la bahía de Zihuatanejo (17°37’22.07” N, 101°33’4.85” O) (fig. 1A). La recolección de ejemplares se realizó con un permiso otorgado por la Secretaría de Agricultura, Ganadería, Desarrollo Rural, Pesca y Alimentación (SAGARPA) (Registro Nacional de Pesca y Acuacultura -DF00000208).

Playa El Palmar es un sitio expuesto, se encuentra frente al complejo turístico Ixtapa y la playa está formada por litoral rocoso y arenoso, tiene una longitud de 2.7 km. Este lugar consta de un relieve heterogéneo, compuesto por riscos y morros de diferente tamaño y forma irregular, la zona intermareal rocosa tiene una amplitud aproximada de 2 m (Aguilar, 2017). En esta playa se han descrito patrones con una circulación del agua dominante hacia el norte a lo largo de la zona, por lo que el oleaje es intenso en la parte norte (Trasviña y Andrade, 2002) (fig. 1B).

Playa Las Gatas es un sitio protegido porque se localiza al interior de la bahía de Zihuatanejo; tiene una extensión de 350 m y se compone principalmente de arenas, fragmentos de coral y rocas (Cisneros, 2016; López, 1993). Paralelo a la línea de costa se encuentra el “pretil”, que es un conglomerado de rocas apiladas, irregulares y de tamaños variables (Aguilar, 2017; Cisneros, 2016; Urbano, 2004). La zona intermareal tiene una amplitud aproximada de 1 m (López, 1993). La presencia de este conglomerado hace que el oleaje sea menor al interior de la playa (Aguilar, 2017) (fig. 1C).

Figura 1. A, Ubicación de la zona de estudio en Ixtapa-Zihuatanejo; B, playa El Palmar; C, playa Las Gatas.

En cada localidad, 3 muestras de macroalgas se recolectaron manualmente de forma aleatoria con una espátula (Bakus, 2007), dentro de cuadros de 20 × 20 cm (0.04 m²) por cada nivel de marea: intermareal bajo (3), medio (3) y alto (3); en total, se obtuvieron 72 muestras provenientes de 36 cuadros en playa Las Gatas y de 36 cuadros en playa El Palmar. Las muestras de las comunidades de macroalgas y moluscos asociados se preservaron en una mezcla de formaldehído al 4% con agua de mar, neutralizada con borato de sodio y glicerina; después se trasladaron al Laboratorio de Ficología (Biodiversidad Marina) de la Facultad de Ciencias de la Universidad Nacional Autónoma de México (UNAM).

En el laboratorio, se calculó la cobertura de cada especie de macroalga en cm² al colocar cada muestra sobre un área delimitada de 20 × 20 cm. Para su identificación taxonómica se tomaron en cuenta las características morfológicas externas, como tipo de talo y tipo de ramificación, e internas a partir de cortes anatómicos transversales de ejes, frondas y, de estar presentes, de estructuras reproductivas; los especímenes y cortes fueron observados bajo microscopios estereoscópico y óptico (Zeiss). La identificación taxonómica de los ejemplares de macroalgas se realizó utilizando literatura especializada para macroalgas del océano Pacífico (Abbott, 1999; Abbott y Hollenberg, 1976; Dawson, 1949, 1953, 1954, 1960, 1961, 1963; Dawson y Beaudette, 1959; Rodríguez et al.,2008; Taylor, 1945). La actualización de la nomenclatura se hizo a partir de la base de datos de Algaebase (Guiry y Guiry, 2024) y, con base en ella, se elaboró una lista sistemática de las especies de macroalgas. De cada comunidad de macroalgas se separó el sedimento retenido por las algas después de la medición de la cobertura para cada especie y se midió su peso húmedo con una balanza digital, modelo OBI.

De cada muestra, se extrajeron de forma manual todos los moluscos de la clase Bivalvia. Los ejemplares con concha y parte blanda (vivos) fueron identificados al nivel taxonómico más bajo posible, género o especie, dependiendo del estado de conservación de cada ejemplar, a partir de la observación de las características morfológicas de la concha, con apoyo de un microscopio estereoscópico. La identificación taxonómica de las especies se hizo con literatura malacológica especializada para la zona del océano Pacífico oriental (Coan et al., 2000; Keen, 1971). Se elaboró una lista sistemática con base en la propuesta de Bouchet et al. (2010) para los niveles suprafamiliares y la actualización de nomenclatura se realizó a partir de la base de datos de World Register of Marine Species (WORMS) para género y especie (Horton et al., 2024).

Los ejemplares de moluscos y macroalgas recolectados fueron depositados en la colección “Invertebrados asociados a macroalgas”, en proceso de registro, con número de inventario para bivalvos (INV-1531 a INV-1638) del Laboratorio de Ficología (Biodiversidad Marina) de la Facultad de Ciencias, UNAM y las macroalgas se depositaron en la colección del Herbario de la Facultad de Ciencias (FCME) con número de catálogo (PTM-10534 a PTM-10558; PTM-10577 a PTM-10604; PTM-10613 a PTM-10630; PTM-10640 a PTM-10648).

Se elaboró una curva de acumulación de las especies de bivalvos asociados a macroalgas recolectadas en las localidades de estudio, con la finalidad de conocer la cantidad de especies que faltaría encontrar y recolectar en dichas localidades. Con los datos de riqueza de especies (S) y abundancia de bivalvos (N) de ejemplares vivos (concha y parte blanda), se estimaron los índices de Shannon (H´) y de diversidad máxima (H’ max) para cada fecha de muestreo, ya que éstos corresponden a los miembros de la comunidad en el momento de muestreo (Aguilar-Estrada et al., 2014); H´ fue expresado en bits/individuo (Magurran, 2004). Se calculó el índice de equidad de Pielou (J´) y el índice de dominancia de Simpson (D) (Moreno, 2001). Estos índices permiten hacer comparaciones cuantitativas y cualitativas entre estudios o zonas ya que se han utilizado como referente mínimo para describir la estructura comunitaria de un lugar (Aguilar-Estrada et al., 2014).

Para evaluar la normalidad de los datos de abundancia de bivalvos y de los índices de diversidad se realizaron pruebas de Shapiro-Wilk (W) (Siegel, 1990). Posteriormente, se realizaron pruebas de Levene (F) para comprobar la homogeneidad de varianzas para los datos no normales y prueba de Bartlett para los datos con distribución normal (Bartlett, 1937; Levene, 1960). Los índices de diversidad de Shannon fueron analizados usando una prueba de “t de student” para evaluar si existían diferencias estadísticamente significativas entre los meses de muestreo. Estos análisis se realizaron utilizando la paquetería Vegan Versión 2.5-6 (Oksanen et al., 2019) en el programa R Studio 2023.12.0+369 (R Core Team, 2023).

Los datos de abundancia de bivalvos no fueron normales de acuerdo con las pruebas de normalidad, por ello, se realizaron análisis de estadística no paramétrica en el software PRIMER v6 + add on Permanova (PRIMER-E Ltd., Plymouth, UK) (Anderson et al., 2008; Clarke y Gorley, 2006). Los datos fueron transformados con ⁴√ y con ellos se calculó una matriz de similitud a partir del índice de Bray-Curtis. A partir de dicha matriz, se realizó un análisis de escalamiento multidimensional no métrico (nMDS) con la finalidad de observar la distribución de las abundancias de bivalvos en el área de estudio. Se incluyó la abundancia de las especies dominantes en forma de vectores azules superpuestos en el gráfico nMDS para facilitar la interpretación de las abundancias de especies dominantes en las distintas localidades y niveles de marea, donde el círculo azul representa la variación de la abundancia. Las especies dominantes que se seleccionaron en este análisis son las que estuvieron presentes en todas las muestras con abundancias mayores a la media.

Además, se realizaron análisis de varianza permutacionales (Permanova) para determinar si existían diferencias significativas en la abundancia de bivalvos (variable dependiente) con respecto a los factores utilizados: localidad con 2 niveles (El Palmar y Las Gatas), nivel de marea con 3 niveles (alto, medio y bajo) y mes de muestreo con 4 niveles (enero, mayo, julio y noviembre). En los análisis de Permanova se utilizaron 999 permutaciones de los residuos bajo un modelo reducido. Posteriormente, para los factores donde se obtuvieron diferencias significativas, se realizaron comparaciones de pares para identificar los niveles que eran diferentes estadísticamente (Anderson et al., 2008).

Utilizando las variables numéricas, se realizaron regresiones lineales múltiples con el software SPSS Statistics v20, para evaluar el efecto de la cobertura de macroalgas y los sedimentos retenidos (variables independientes) en la abundancia de bivalvos (variable dependiente) en las playas El Palmar y Las Gatas, en los diferentes niveles de marea y en los meses de muestreo.

Por último, se determinó el grupo morfofuncional (GMF): filamentosas (Fil), foliosas (Fol), foliosas corticadas (Foc), filamentosas corticadas (Fic), coriáceas (Cor), calcáreas articuladas (Cal) para cada especie de macroalgas, con base en la propuesta de Steneck y Dethier (1994) y se asociaron con las especies de bivalvos recolectadas.

Resultados

Riqueza y composición de moluscos. En 60 de las 72 muestras recolectadas se encontraron bivalvos, 35 muestras pertenecientes a playa El Palmar y 25 a playa Las Gatas. Se obtuvieron un total de 873 individuos de la clase Bivalvia. Se identificaron 17 especies (fig. 2) agrupadas en 2 subclases, 8 órdenes, 10 familias y 15 géneros. Del total de especies recolectadas, 2 de ellas Parvilucina approximata (Dall, 1901) y Pinna rugosa G. B. Sowerby I, 1835, encontradas en playa El Palmar y playa Las Gatas, respectivamente, no se incluyeron en los análisis de estructura comunitaria, dado que solo se recolectó la concha. Se obtuvieron 13 especies en playa El Palmar y 9 en playa Las Gatas. Las familias Carditidae y Mytilidae fueron las mejor representadas con 3 (18%) y 5 (30%) especies, respectivamente.

En ambas localidades, las especies más abundantes fueron: Brachidontes adamsianus (Dunker, 1857), Leiosolenus aristatus (Dillwyn, 1817) e Isognomon janus Carpenter, 1857, con 376, 205 y 124 individuos, respectivamente, lo cual representó 91% del total de individuos. Brachidontes adamsianus obtuvo la mayor abundancia en El Palmar a lo largo del ciclo anual, con el valor más elevado (161 individuos) en enero. En Las Gatas, Leiosolenus aristatus fue la especie con la mayor abundanciacon 140 individuos en enero; sin embargo, la especie más abundante fue B. adamsianus a lo largo de los otros meses de muestreo. Mientras que las especies menos abundantes fueron: Sphenia fragilis (H. Adams & A. Adams, 1854) con 2 individuos en enero y mayo en El Palmar y 1 individuo en noviembre en Las Gatas, Crassinella ecuadoriana Olsson, 1961 con 2 individuos en mayo y julio en El Palmar, Linucula declivis (Hinds, 1843) con 2 individuos en julio y noviembre en Las Gatas, Mya sp. y Mytilus edulis Linnaeus, 1758 con 1 individuo en enero en El Palmar (tabla 1).

Tabla 1

Abundancia de las especies de bivalvos recolectadas en las localidades de estudio por mes de muestreo. El arreglo sistemático sigue la propuesta de Horton et al. (2024).

Especie	Playa El Palmar	Playa Las Gatas
	ene	may	jul	nov	ene	may	jul	nov	Total
Subclase Protobranchia
Orden Nuculida
Familia Nuculidae
Linucula declivis	–	–	–	–	–	–	1	1	2
Subclase Autobranchia
Orden Mytilida
Familia Mytilidae
Brachidontes adamsianus	161	43	18	14	41	9	34	56	376
Brachidontes semilaevis	4	–	–	1	–	–	–	–	5
Leiosolenus aristatus	2	30	1	6	140	7	5	14	205
Modiolus capax	1	2	1	1	–	–	–	–	5
Mytilus edulis	1	–	–	–	–	–	–	–	1
Orden Arcida
Familia Arcidae
Acar rostae	–	–	–	–	–	–	1	3	4
Orden Ostreida
Familia Pteriidae
Isognomon janus	51	32	8	5	2	–	10	16	124
Familia Pinnidae
Pinna rugosa	–	–	–	–	–	–	3	–	3
Orden Lucinida
Familia Lucinidae
Parvilucina approximata	1	–	–	–	–	–	–	–	1
Orden Carditida
Familia Carditidae
Carditamera affinis	5	2	–	–	1	–	–	1	9
Carditamera radiata	–	3	1	–	3	–	3	5	15
Cardites grayi	–	1	3	–	–	–	–	–	4
Familia Crassatellidae
Crassinella ecuadoriana	–	1	1	–	–	–	–	–	2
Orden Venerida
Familia Chamidae
Chama coralloides	3	2	1	–	–	–	–	7	13
Orden Myida
Familia Myidae
Mya sp.	1	–	–	–	–	–	–	–	1
Sphenia fragilis	1	1	–	–	–	–	–	1	3
Total	231	117	34	27	187	16	57	104	773

Figura 2. Especies de bivalvos asociadas a macroalgas recolectadas en Ixtapa-Zihuatanejo, Guerrero, en vista dorsal y ventral. a, Linucula declivis;b, Brachidontes adamsianus;c, Leiosolenus aristatus; d, Mytilus edulis; e, Acar rostae; f, Isognomon janus;g, Pinna rugosa; h, Parvilucina approximata; i, Carditamera radiata; j, Cardites grayi; k, Crassinella ecuadoriana; l, Chama coralloides;m, Sphenia fragilis; n, Modiolus capax;o, Brachidontes semilaevis;p, Carditamera affinis;q, Mya sp.

La riqueza de especies con relación al ciclo anual fue de 13 especies en playa El Palmar, variando entre 5 (noviembre) y 10 especies (enero) a lo largo de los meses, mientras que en playa Las Gatas se obtuvieron 9 especies, con una variación de 2 (mayo) y 9 (noviembre) especies (tabla 1). La curva de acumulación de especies mostró un comportamiento asintótico en El Palmar, mientras que Las Gatas no mostró una tendencia a ser asintótica (fig. 3).

Figura 3. Curva de acumulación de especies de bivalvos registrados en las localidades de estudio. Línea azul: playa El Palmar, línea verde: playa Las Gatas.

En playa El Palmar el nivel intermareal bajo obtuvo el mayor número de especies (13) y el menor número de especies se presentó tanto en el nivel intermareal medio como en el alto (7); mientras que en playa Las Gatas, se presentó el mayor número de especies en el nivel intermareal medio (10) y el nivel alto tuvo el menor número de especies (5). Las especies que se presentaron en ambas localidades en los 3 niveles de la zona intermareal fueron Brachidontes adamsianus, Leiosolenus aristatus e Isognomon janus.

En El Palmar se registraron valores del índice de diversidad de Shannon que fluctuaron entre 1.30 y 2.15 bits/individuo, mientras que en Las Gatas este índice varió entre 0.98 y 2.09 bits/individuo (tabla 2). En general, el índice de equidad de Pielou presentó valores mayores a 0.60 en playa El Palmar, excepto en enero, mientras que en playa Las Gatas se observó que los valores fueron bajos en enero y mayo, y mayores a 0.60 en julio y noviembre. Las pruebas de “t de student”para el índice de diversidad de Shannon entre los meses de muestreo, indicaron diferencias significativas entre todos ellos (p ≤ 0.05) en El Palmar, y Las Gatas (tabla 3).

Tabla 2

Índices comunitarios calculados para las especies recolectadas en playa El Palmar y playa Las Gatas. N = número de individuos, H’= Índice de Shannon, J’= equidad de Pielou y λ = dominancia de Simpson. Los valores máximos y mínimos por localidad para el índice de Shannon (H´) se marcan con negritas.

	Playa El Palmar	Playa Las Gatas
Mes	N	H’	J’	λ	N	H’	J’	λ
Enero	230	1.30	0.48	0.54	187	0.99	0.41	0.60
Mayo	117	2.15	0.81	0.27	16	0.98	0.41	0.50
Julio	34	2.03	0.77	0.34	54	1.63	0.68	0.44
noviembre	27	1.77	0.67	0.35	104	2.09	0.87	0.33

En el análisis de escalamiento multidimensional no métrico (nMDS), a partir de los valores de abundancia de las especies de bivalvos, se observó que las muestras se separan en 3 grupos: un primer grupo dominado por estaciones pertenecientes a playa El Palmar (P), un segundo grupo principalmente con estaciones de playa Las Gatas (G) y uno tercero con estaciones de ambas localidades que probablemente se separa por la dominancia de la especie Brachidontes adamsianus (fig. 4). El vector indica la dirección a través del plano de ordenación en la cual aumentan los valores de abundancia de las especies dominantes, la longitud de la línea indica la cantidad de variación total de cada especie, entonces, si toda la variación se explicara, la línea azul alcanzaría el círculo azul.

Tabla 3

Resultado de la prueba de “t de student”para evaluar diferencias significativas del índice de Shannon entre los meses de muestreo por localidad.

	Playa El Palmar	Playa Las Gatas
Meses	H´	t	p	gl	H´	t	p	gl
enero vs. mayo	1.30/2.15	-0.009	324	1	0.99/0.98	0.009	169	1
mayo vs. julio	2.15/2.03	0.0004	142		0.98/1.63	-19.84	4
julio vs. noviembre	2.03/1.17	0.023	60		1.63/2.09	-3.18	94

El análisis de Permanova reveló que la composición de bivalvos difiere significativamente entre las diferentes localidades muestreadas (pseudo-F = 5.2373, p = 0.001). No hubo diferencias estadísticamente significativas para los meses de muestreo (pseudo-F=1.431, p = 0.147) y el nivel de marea (pseudo-F = 1.19, p = 0.309) (tabla 4).

Tabla 4

Resultado de los análisis de Permanova para los factores de localidad, mes de muestreo y nivel de marea. GL: Grados de libertad, ms: media suma de cuadrados, perm: permutaciones, * = valores significativos.

Factor	GL	ms	Pseudo-F	p (perm)
Localidad	1	7,844.9	5.2373	0.001*
Res	57	1,497.9
Total	58
Mes	3	2,249.9	1.431	0.147
Res	55	1,572.3
Total	58
Nivel de marea	2	1,915	1.1996	0.309
Res	56	1,596.3
Total	58

Riqueza y composición de macroalgas. Se encontró un total de 3 phyla, 3 clases, 5 subclases, 14 órdenes, 20 familias, 33 géneros y 59 especies de macroalgas; de las cuales, a cada localidad le pertenecen 37 especies. De las 59 especies totales, 11 fueron Chlorophyta (18%), 42 Rhodophyta (70%) y 7 Heterokontophyta-Phaeophyceae (11%). Entre ambas localidades se compartieron 16 especies. Las familias con mayor número de especies fueron Rhodomelaceae con 12 especies, Corallinaceae con 10 y Ceramiaceae y Dictyotaceae con 5 especies cada una (tabla 5).

Asociación entre grupos morfofuncionales y moluscos. Las macroalgas se clasificaron en 6 grupos morfofuncionales: algas filamentosas, foliosas, foliosas corticadas, filamentosas corticadas, coriáceas y calcáreas articuladas (tabla 5). Los grupos dominantes en las 2 localidades fueron las algas filamentosas, filamentosas corticadas y las calcáreas articuladas. En playa El Palmar se presentó una mayor riqueza de especies (14) de bivalvos al interior de las comunidades algales que contenían el grupo de algas calcáreas articuladas. En playa Las Gatas, los grupos morfofuncionales que presentaron mayor cantidad de especies fueron las algas filamentosas, filamentosas corticadas y las calcáreas articuladas con 10 especies cada uno (tabla 6).

Cobertura de macroalgas, abundancia de moluscos y sedimentos retenidos. La regresión lineal por localidad indicó que para playa El Palmar hay una correlación baja con 24% entre los componentes: cobertura de macroalgas, abundancia de bivalvos y peso de sedimento húmedo; sin embargo, se observó una relación moderada entre los sedimentos retenidos por las macroalgas y la abundancia de bivalvos (R = 0.55 y R² = 0.24). Para playa Las Gatas la regresión lineal múltiple indicó una correlación baja con 10% entre los factores y una relación moderada entre la cobertura de macroalgas y la abundancia de bivalvos (R = 0.42 y R² = 0.10). La cobertura de macroalgas, abundancia de bivalvos y los sedimentos retenidos varió según la localidad (fig. 5A-C), mes de muestreo (fig. 5D-F) y nivel de marea (fig. 5G-I). En El Palmar la regresión lineal mostró una relación directamente proporcional entre la abundancia de bivalvos y la cobertura de macroalgas. Los valores más altos de cobertura de macroalgas (0.89 m²) se relacionaron con los valores más altos de abundancia de bivalvos.

Con respecto al mes de muestreo, las regresiones lineales múltiples indicaron una correlación alta con 98% (R = 0.99 y R² = 0.98), entre la cantidad de sedimento retenido y la abundancia de bivalvos. La regresión lineal múltiple sugirió que la abundancia de bivalvos disminuyó a lo largo del ciclo anual, donde enero presentó la mayor abundancia (de 187-230 individuos) y noviembre la menor (27 individuos). La relación de la abundancia de bivalvos con el peso húmedo de sedimento tuvo un comportamiento parecido, disminuyendo conforme al ciclo anual: en las muestras de enero se reportó mayor sedimento retenido (0.56 kg), mientras que en julio se encontraron menos sedimentos retenidos (0.17 kg).

Figura 4. Análisis de escalamiento multidimensional no métrico (nMDS) de las especies de bivalvos recolectadas en las localidades de estudio. Las líneas azules (vectores) indican las especies abundantes y frecuentes (dominantes). El círculo azul representa la variación de las abundancias.

Tabla 5

Especies de macroalgas asociadas a bivalvos, para cada localidad y grupo morfofuncional. Se marca con * la especie de macroalga presente en la localidad. Localidad: playa El Palmar (P), playa Las Gatas (G). Grupos morfofuncionales de macroalgas: filamentosas (Fil), foliosas (Fol), foliosas corticadas (Foc), filamentosas corticadas (Fic), coriáceas (Cor), calcáreas articuladas (Cal). Especies de bivalvos: 1) Linucula declivis, 2) Brachidontes adamsianus, 3) Brachidontes semilaevis, 4) Leiosolenus aristatus, 5) Modiolus capax, 6) Mytilus edulis, 7) Acar rostae, 8) Isognomon janus, 9) Carditamera affinis, 10) Carditamera radiata, 11) Cardites grayi, 12) Crassinella ecuadoriana, 13) Chama coralloides, 14) Mya sp., 15) Sphenia fragilis,16) Pinna rugosa, 17) Parvilucina approximata.

Especies	Grupo morfofuncional	Especies de bivalvos
		P	G
Chlorophyta
Ulvales
Ulvaceae
Ulva californica Wille, 1899	Fol	*	–
Ulva intestinalis Linnaeus, 1753	Fol	–	2, 4, 13
Ulva linza Linnaeus, 1753	Fol	–	2, 4, 13
Bryopsidales
Bryopsidaceae
Bryopsis pennata var. minor J. Agardh, 1887	Fil	–	2, 4, 8
Caulerpaceae
Caulerpa chemnitzia (Esper) J.V. Lamouroux, 1809	Fil	2, 4, 8, 10, 17	1, 2, 4, 7, 8, 9, 10, 15
Caulerpa sertularioides (S.G. Gmelin) M. Howe, 1905	Fil	2, 4, 8, 10, 17	1, 2, 4, 7-10, 15
Halimedaceae
Halimeda discoidea Decaisne, 1842	Cal	2, 10, 17	–
Cladophorales
Cladophoraceae
Chaetomorpha antennina (Bory) Kützing, 1847	Fil	2, 3, 7, 8-13	4, 7, 8, 11, 13
Cladophora sp.	Fil	2, 4, 5, 8, 9	–
Cladophora graminea Collins, 1909	Fil	–	16
Cladophora microcladioides Collins, 1909	Fil	2, 4, 5, 8, 9	–
Rhodophyta
Gigartinales
Cystocloniaceae
Hypnea pannosa J. Agardh, 1847	Fic	2-6, 8-15	2, 4, 8-10, 16
Hypnea spinella (C. Agardh) Kützing, 1847	Fic	2-6, 8-15, 17	2, 4, 8-10, 16
Hypnea johnstonii Setchell y N. L. Gardner, 1924	Fic	–	1, 2, 4, 7-10, 16
Phyllophoraceae
Ahnfeltiopsis gigartinoides (J. Agardh) P. C. Silva y DeCew, 1992	Fic	8	–
Gymnogongrus johnstonii (Setchell y N. L. Gardner) E.Y. Dawson, 1961	Fic	2, 8	–
Ceramiales
Ceramiaceae
Centroceras clavulatum (C. Agardh) Montagne, 1846	Fil	2-5, 8	2, 4, 10, 15

Tabla 5. Continúa
Especies	Grupo morfofuncional	Especies de bivalvos
		P	G
Ceramium sp.	Fil	–	1, 2, 4, 7, 8, 10, 13, 16
Ceramium camouii E. Y. Dawson, 1944	Fil	–	1, 2, 4, 7, 8, 10, 13, 16
Ceramium zacae Setchell y N. L. Gardner, 1937	Fil	–	1, 2, 4, 7, 8, 10, 13, 16
Gayliella flaccida (HarveyyKützing) T.O. McIvory y L.J. Cho, 2008	Fic	–	2, 4, 7, 8, 11
Delesseriaceae
Taenioma perpusillum (J. Agardh) J. Agardh, 1863	Fil	–	4, 16
Rhodomelaceae
Chondria sp.	Fic	*	*
Herposiphonia secunda (C. Agardh) Ambronn, 1880	Fic	8	2, 4, 8, 10
Herposiphonia tenella (C. Agardh) Ambronn, 1880	Fic	7	–
Melanothamnus simplex (Hollenberg) Díaz-Tapia y Maggs, 2017	Fic	2, 8	–
Melanothamnus sphaerocarpus (Borgesen) Díaz.Tapioa y Maggs, 2017	Fic	–	2, 4, 8
Polysiphonia mollis J. D. Hooker y Harvey, 1847	Fic	*	2, 4, 8, 15
Polysiphonia nathanielii Hollenberg, 1958	Fic	–	2, 4, 8, 15
Polysiphonia subtilissima Montagne, 1840	Fic	–	2, 4, 8, 15
Eutrichosiphonia confusa (Hollenberg) Savoie y G.W. Saunders, 2019	Fic	–	2, 4, 8
Laurencia sp.	Fic	2, 8,	2, 7, 8, 13
Laurencia hancockii E.Y. Dawson, 1944	Fic	2, 8	–
Laurencia subcorymbosa E.Y. Dawson, 1963	Fic	2, 8	–
Rhodymeniales
Rhodymeniaceae
Tayloriella dictyurus (J. Agardh) Kylin, 1956	Fil	2, 8	–
Lomentariaceae
Ceratodictyon tenue (Setchell y N. L. Gardner) J.N. Norris, 2014	Fil	–	2, 4, 8
Corallinales
Corallinaceae
Amphiroa beauvoisii J.V. Lamouroux, 1816	Cal	2-15, 17	1, 2, 4, 7-14, 16
Amphiroa misakiensis Yendo, 1902	Cal	2-15	1, 2, 4, 7-14, 16
Amphiroa rigida J.V. Lamouroux, 1816	Cal	2-15	1, 2, 4, 7-14, 16
Amphiroa subcylindrica E.Y. Dawson, 1953	Cal	2-15	–
Jania capillacea Harvey, 1853	Cal	–	2, 4, 7-10, 15, 16
Jania subpinnata E.Y. Dawson, 1953	Cal	2, 4, 5, 8, 9, 12, 14	2, 4, 7-10, 15, 16
Jania tenella (Kützing) Grunow, 1874	Cal	2, 4, 5, 8, 9, 12, 14	2, 4, 7-10, 15, 16
Jania tenella var. tenella	Cal	2, 4, 5, 8, 9, 12, 14	–
Gelidiales
Gelidiaceae
Gelidiella acerosa (Forsskal) Feldmann y Hamel, 1934	Fic	–	2, 7, 8
Gelidium mcnabbianum (E.Y. Dawson) B. Santelices, 1998	Fic	–	2, 4, 7-11
Gelidium pusillum (Stackhouse) Le Jolis, 1863	Fic	2, 4, 5, 8, 12, 14	2, 4, 7-11
Pterocladiaceae
Pterocladiella caloglossoides (M. Howe) Santelices, 1998	Fic	–	2, 4, 7, 8, 15
Gracilariales
Gracilariaceae
Gracilaria sp.	Fic	2, 8, 12, 14	–
Halymeniales
Halymeniaceae
Grateloupia huertana Mateo-Cid, Mendoza-González y Gavio, 2005	Fic	–	4, 8
Grateloupia versicolor (J. Agardh) J. Agardh, 1847	Fic	2, 4, 8	–
Heterokontophyta-Phaeophyceae
Ectocarpales
Scytosiphonaceae
Chnoospora minima Papenfuss, 1956	Fic	*	–
Fucales
Sargassaceae
Sargassum liebmannii Agardh, 1847	Cor	2, 4, 5, 7, 8, 10-13, 15, 17	–
Dictyotales
Dictyotaceae
Dictyota sp.	Foc	–	*
Dictyota dichotoma (Hudson) J.V. Lamouroux, 1809	Foc	2, 4, 5, 8	–
Lobophora variegata (J.V.Lamouroux)Womersley ex E.C.Oliveira, 1977	Fic	8	–
Padina mexicana var. erecta Avila-Ortiz, 2003	Foc	2, 4, 5, 8, 9	–
Padina ramonribae Avila-Ortiz, Pedroche y Díaz-Martínez, 2016	Foc	2, 4, 5, 8, 9	–

La regresión lineal por nivel de marea indicó una correlación significativa con 100% (R = 1 y R² = 1), es decir, la cobertura de algas y el sedimento retenido están relacionados con la abundancia de bivalvos en los diferentes niveles de marea. El nivel del intermareal bajo presentó los valores más elevados de cobertura de macroalgas con 0.61 m² (fig. 5G), mayor abundancia de moluscos con 275 individuos (fig. 5H) y mayor cantidad de sedimentos con 0.69 kg (fig. 5I). En el intermareal alto se registraron los valores más bajos de abundancia de bivalvos con 230 individuos (fig. 5H), cobertura de algas con 0.52 m² (fig. 5G) y peso húmedo de sedimentos con 0.30 kg (fig. 5I), es decir, conforme se asciende en el intermareal, estas variables disminuyen.

Tabla 6

Asociación de las especies de bivalvos con los grupos morfofuncionales propuestos por Steneck y Dethier (1994). Grupos morfofuncionales de macroalgas: filamentosas (Fil), foliosas (Fol), foliosas corticadas (Foc), filamentosas corticadas (Fic), coriáceas (Cor), calcáreas articuladas (Cal). × = Presencia de la especie.

Familia

Especie

Playa El Palmar

Playa Las Gatas

Fil

Fol

Foc

Fic

Cor

Cal

Fil

Fol

Foc

Fic

Cor

Cal

Nuculidae

L. declivis

–

–

–

–

–

–

×

–

–

×

–

×

Mytilidae

B. adamsianus

×

–

×

×

×

×

×

×

–

×

–

×

B. semilaevis

×

–

–

×

–

×

–

–

–

–

–

–

L. aristatus

×

–

×

×

×

×

×

×

–

×

–

×

M. capax

×

–

×

×

×

×

–

–

–

–

–

–

M. edulis

–

–

–

×

–

×

–

–

–

–

–

–

Arcidae

A. rostae

×

–

–

–

×

×

×

–

–

×

–

×

Pteriidae

I. janus

×

–

×

×

×

×

×

–

–

×

–

×

Pinnidae

P. rugosa

–

–

–

–

–

–

×

–

–

×

–

×

Lucinidae

P. approximata

×

–

–

×

×

×

–

–

–

–

–

–

Carditidae

C. affinis

×

–

×

×

–

×

×

–

–

×

–

×

C. radiata

×

–

–

×

×

×

×

–

–

×

–

×

C. grayi

×

–

–

×

×

×

×

–

–

×

–

×

Crassatellidae

C. ecuadoriana

×

–

–

×

×

×

–

–

–

–

–

–

Chamidae

C. coralloides

×

–

–

×

×

×

×

×

–

×

–

×

Myidae

Mya sp.

–

–

–

×

–

×

–

–

–

–

–

–

S. fragilis

–

–

–

×

×

×

×

–

–

×

–

×

Total

12

0

5

14

11

15

11

3

0

11

0

11

Discusión

Riqueza y composición de moluscos. De las 17 especies encontradas en el presente estudio, 16 ya habían sido observadas para el Pacífico tropical mexicano, específicamente para Ixtapa-Zihuatanejo en sustratos rocosos y arenosos (Flores-Rodríguez et al., 2007, 2012; Lesser, 1984; López-Rojas et al.,2017; Lozada, 2010; Sánchez, 2014), excepto Mytilus edulis que es un nuevo registro para Guerrero. Esta especie fue encontrada en playa El Palmar, asociada a los géneros de macroalgas Amphiroa e Hypnea. Dicha especie había sido registrada para el océano Pacífico desde las costas del Ártico, Canadá, EUA, hasta Cabo San Lucas, México, al interior del golfo de California en San Luquitas y Santa Rosalía, y para el Pacífico tropical mexicano en isla Socorro, por lo que su distribución se amplía hacia el sur, a Guerrero en ambientes litorales asociados a las comunidades de macroalgas (Cadena-Cárdenas et al., 2009; Fitch, 1953; Shaw et al., 1988). La ampliación de la distribución de esta especie puede deberse al movimiento de la corriente de California como ha sido reportado para bivalvos (Schulien et al., 2020), gasterópodos y escafópodos que ampliaron su distribución desde la provincia Californiana hacia la Panámica donde se señala la confluencia de especies de moluscos entre ambas provincias (Landa-Jaime y Arciniega-Flores 1998; Ríos-Jara et al., 2003; Gama-Kwick et al., 2021). Otra de las razones por la cual M. edulis amplió su distribución al sur podría ser por el fenómeno de El Niño que proporciona condiciones ambientales para que las especies que habitan en sitios templados y subtropicales se desplacen a sitios tropicales (Díaz y Ortlieb, 1993; Paredes et al., 1998; Velez y Zeballos, 1985).

Figura 5. Variación de cobertura de macroalgas, abundancia de bivalvos y sedimento retenido por localidades (A, B y C), mes de muestreo (D, E y F) y nivel de la zona intermareal (G, H e I). Verde = macroalgas, amarillo = bivalvos, rojo = sedimento retenido.

La familia mejor representada fue Mytilidae con 5 especies: Brachidontes adamsianus, B. semilaevis, Leiosolenus aristatus, Modiolus capax y Mytilus edulis (tabla 1), dichas especies coinciden con lo registrado por Flores-Garza et al.(2014), Galeana-Rebolledo et al.(2012) y Garcés (2011) en Acapulco, así como por López-Rojas et al. (2017) en diferentes localidades en Guerrero. Mytilidae también ha sido la familia mejor representada en otras localidades de Zihuatanejo (muelle municipal), donde se han observado 8 especies (Guzmán, 2022), así como en otras regiones tropicales como Brasil, donde esta familia también fue la mejor representada con 5 especies (Santos et al., 2020).

Una mayor riqueza de especies de la familia Mytilidae en Ixtapa-Zihuatanejo podría explicarse debido a la frecuencia y rápido asentamiento de larvas de bivalvos, ya que son organismos que están en constante reproducción a lo largo del año (Seed, 1969a). Suchanek (1978) mencionó que los miembros de la familia Mytilidae tienden a colonizar rápidamente los espacios disponibles en la zona intermareal rocosa, familia que se ha adaptado a distintos hábitats (Keen, 1971), gracias a su rápido aumento de talla y posterior asentamiento de larvas (Ceccherelli y Rossi, 1984). También se ha mencionado que los organismos de dicha familia pueden producir fibras de biso, por lo que son capaces de anclarse y permanecer en sustratos en los que otras familias de bivalvos no pueden, como rocas, arenas u otros organismos (Keen, 1971; Stella et al.,2010).

En general, se encontró una mayor riqueza de especies en El Palmar, valores que pueden explicarse debido a la heterogeneidad ambiental de una playa expuesta con oleaje intenso como dicho sitio (Morales et al., 2008). Se ha mencionado que la riqueza de especies de bivalvos es mayor en zonas expuestas en donde las olas impactan directamente, por su parte, las zonas protegidas o de menor oleaje presentan menor riqueza específica (Flores-Garza et al., 2014; Flores-Rodríguez et al., 2012; Valdés-González et al., 2004).

En El Palmar, se encontraron valores elevados para la riqueza de especies de bivalvos (17 spp.) respecto de lo encontrado en Las Gatas, lo anterior también se ha observado para otras clases de moluscos como los gasterópodos y poliplacóforos asociados a comunidades de macroalgas en playa El Palmar (Aguilar, 2017). También en esta playa se encontró menor cantidad de sedimentos retenidos por las macroalgas (fig. 5C). El Palmar está conformada por rocas de diferentes tamaños, que pueden moverse conforme se dan los cambios en el oleaje y la marea, lo que le confiere un alto dinamismo que no permite la acumulación de sedimentos (Gibbons, 1988), tal como se observó en el presente estudio. Dichas características ayudan a que se conformen una amplia gama de ambientes que pueden ser colonizados, donde se esperaría una mayor riqueza de especies como se observó para los bivalvos de El Palmar (Benedetti-Cecchi, 2001).

Como se ha mencionado anteriormente, playa Las Gatas es un sitio protegido, por esta razón es un lugar de baja energía, con características poco favorables para el establecimiento de especies de bivalvos de manera muy similar a los estudios realizados en otras zonas de Guerrero (Flores-Garza et al., 2014; Flores-Rodríguez et al., 2012; Valdés-González et al., 2004). Las fluctuaciones observadas entre ambas playas de la riqueza de especies, también podrían deberse a otro tipo de características presentes en algunas familias de bivalvos que no tienen hábitos epifaunales como los que se adhieren a las macroalgas y que presentan hábitos de vida semiinfaunales o infaunales (Garcés, 2011), por lo que estas especies podrían encontrarse en otros sitios, posiblemente en zonas más profundas en la zona submareal, o de igual forma, las especies de bivalvos se pueden ver afectadas por la contaminación al incorporar en sus tejidos bacterias patógenas (Gosling, 2015), que es común de una zona turística como Ixtapa-Zihuatanejo (IMTA, 2010; UNAM, 2013). Lo anterior podría explicar las fluctuaciones en la riqueza de especies de bivalvos observadas en Zihuatanejo.

Los valores de riqueza de especies no variaron considerablemente en las localidades estudiadas a lo largo del año. La riqueza de especies de bivalvos en un ciclo anual puede ser constante, ya que la mayoría de las especies (70%) de bivalvos del presente trabajo son generalmente inmóviles o sedentarias. Lo anterior hace que sus poblaciones se mantengan con pocos cambios con respecto a su riqueza de especies debido a sus patrones de reproducción anuales o bianuales (Baqueiro y Masso, 1988; Flores-Rodríguez et al., 2012; Seed, 1969a).

La especie más abundante en el presente trabajo fue B. adamsianus (Mytilidae). Esta especie ha sido registrada como la más abundante en otras localidades de Zihuatanejo (Guzmán, 2022). En particular, las familias Chamidae y Mytilidae se han registrado como las más abundantes para Guerrero; por su parte, la especie Isognomon janus Carpenter, 1857 (Isognomonidae) también ha presentado un gran número de organismos (Guzmán, 2022; López-Rojas et al. 2017). Las bajas abundancias de los bivalvos encontrados en este trabajo pueden estar dadas por el tipo de hábitos de vida de dichos moluscos. Por ejemplo, el género Sphenia tienen un hábito endolítico (Esqueda-González et al., 2014; Garcés, 2011; Guzmán, 2022), por lo que sería poco probable encontrarlos dentro de las comunidades de macroalgas, ya que es común encontrarlos incrustados en agujeros preexistentes de rocas, fragmentos de madera, incluso otros materiales como las conchas de moluscos, así como entre colonias de briozoos (Coan, 1999). El Permanova mostró que existen diferencias significativas entre localidades, esto puede deberse a que playa Las Gatas es considerada como un sitio protegido y playa El Palmar como un sitio expuesto, esto representó diferencias en abundancia y riqueza de especies, que pueden deberse a los requerimientos medioambientales propios de cada una de las especies de bivalvos como: sustrato, alimento, salinidad o temperatura (Borges et al., 2014; Galeana-Rebolledo et al., 2012; Seed 1969a).

En playa El Palmar, la curva de acumulación de especies fue asintótica, lo que sugiere que se encontró a la mayoría de los bivalvos asociados a macroalgas de esta localidad. Mientras que, en playa Las Gatas, la curva de acumulación de especies no fue asintótica. La cantidad de especies recolectadas en un sitio está relacionada con el esfuerzo de muestreo (Moreno, 2001). Por lo que el menor número de muestras en donde se encontraron bivalvos en playa Las Gatas (25), en comparación con las de playa El Palmar (35), podría explicarse por este hecho. Se ha mencionado que una curva de acumulación de especies en muy raras ocasiones llega a ser asintótica, ya que siempre habrá especies que no se recolecten o durante los muestreos pueden encontrarse especies raras, lo que puede estar determinado por el sitio de muestreo, temporada del año, tipo de sustrato, entre otras variables (Jiménez-Valverde y Hortal, 2003).

Además, otra de las razones para que la curva de acumulación de especies no haya sido asintótica para Las Gatas, podría deberse a que, en ambas localidades, las muestras de bivalvos fueron obtenidas de comunidades de macroalgas, que es un sustrato muy específico y es distinto a las rocas y arena que componen la zona intermareal que corresponden a los sustratos que se han estudiado en la mayoría de los distintos trabajos malacológicos de Guerrero. Por lo anterior, el inventario podría estar incompleto; sin embargo, es una buena aproximación de la biodiversidad de las especies de bivalvos asociadas a comunidades de macroalgas de la zona norte de Guerrero en Ixtapa-Zihuatanejo. Para tener una curva de acumulación de especies asintótica se recomienda muestrear todos los sustratos posibles y de esta forma, obtener casi todas las especies de bivalvos de la zona.

Se encontraron cambios con respecto a la abundancia de bivalvos asociados a macroalgas a lo largo del año, fluctuaciones que pudieran estar relacionadas con sus patrones reproductivos, lo que podría explicar la disminución en la abundancia de bivalvos en los meses cálidos (mayo y julio) y un aumento en la cantidad de individuos en los meses fríos (enero y noviembre). Los cambios en la abundancia de bivalvos se han estudiado en trabajos sobre reproducción de diferentes especies de moluscos (Baqueiro y Aldana, 2000, 2003). Algunas de las especies de bivalvos que han sido utilizadas en estudios reproductivos o sobre ciclos gonádicos son: Chione undatella (G. B. Sowerby I, 1835), Megapitaria aurantiaca (G. B. Sowerby I, 1831) (Veneridae), Brachidontes rodriguezii (d’Orbigny, 1842), Mytilus edulis y Mytilus chilensis Hupé, 1854 (Mytilidae), así como Anadara tuberculosa (G. B. Sowerby I, 1833) (Arcidae), tanto en México como otros países de Sudamérica. Dichos estudios destacan que los bivalvos se reproducen continuamente a lo largo del año y se ha mencionado que el momento del desove comprende desde la primavera hasta el otoño con temperaturas por encima de los 25 °C; se ha observado que las larvas se instalan en octubre y noviembre, con una metamorfosis de 15 a 25 días, dichas larvas permanecen adheridas a las algas de noviembre a mayo, en donde llevan a cabo su metamorfosis (Aguillón, 2011; Baqueiro y Masso, 1988; García-Domínguez et al., 2008; Hernández-Moreno et al., 2020; Oyarzún et al., 2011; Seed, 1969a; Torroglosa, 2015). Solo se encontró a Mytilus edulis como especie de importancia comercial para el presente trabajo. En especies que no se les ha considerado importantes para la industria pesquera, es muy poco lo que se conoce acerca de sus hábitos reproductivos (Aguilar, 2017).

Las 3 especies más abundantes del presente trabajo son epifaunales de hábitos filtradores por suspensión (García-Cubas, 1981), Isognomon janus (Isognomonidae) habita en la zona intermareal hasta profundidades de 20 m, Leiosolenus aristatus (Mytilidae)puede encontrarse desde el litoral hasta los 300 m sobre rocas u otros bivalvos (Coan y Valentich-Scott, 2012) y Brachidontes adamsianus (Mytilidae)habita sobre grietas de rocas grandes en zonas expuestas (Landa-Jaime et al., 2013). Las especies que no pertenecen a las familias Isognomonidae y Mytilidae, en general, son organismos infaunales de hábitos filtradores por suspensión, a excepción de la familia Myidae que son infaunales perforadores de hábitos filtradores por suspensión (García-Cubas, 1981; Coan, 1999).

Los bivalvos están presentes en diferentes lugares al interior de la zona intermareal rocosa y hay especies que tienen preferencia por algún sitio a lo largo de dicha zona (Román-Contreras et al.,1991; Sibaja-Cordero y Vargas-Zamora, 2006; Suchanek, 1978). En el presente estudio se observó un incremento en la abundancia y riqueza en los niveles intermareales medio y bajo. Este mismo patrón se ha observado en otras especies de moluscos como gasterópodos y poliplacóforos de Zihuatanejo (Aguilar, 2017; Gama-Kwick et al., 2021).

El índice de diversidad de Shannon obtenido en este trabajo fue bajo, menor a 2.16 bits/individuo para ambos sitios de muestreo, en comparación con otros trabajos sobre bivalvos asociados a macroalgas para Zihuatanejo, como el de Guzmán (2022), donde se encontraron valores de 3.41 bits/individuo en el muelle municipal. Flores-Garza et al.(2014), en su estudio en Acapulco, presentaron valores de 3.65 bits/individuo y Galeana-Rebolledo et al.(2012) reportaron valores de 3.64 bits/individuo. Las diferencias con respecto al índice de diversidad de Shannon pueden deberse a la baja riqueza de especies encontrada en Ixtapa-Zihuatanejo, respecto de otras localidades del Pacífico tropical mexicano, donde se han registrado valores superiores a 17 especies; Reguero y García-Cubas (1989) encontraron 53 especies de bivalvos para Nayarit, y Esqueda-González et al. (2014) reportaron 89 especies de bivalvos para Sinaloa. Lo anterior también puede relacionarse con los altos valores de abundancia de ciertas especies, como es el caso de Brachidontes adamsianus. Ésto también se ve reflejado en los valores calculados del índice de dominancia de Simpson, el cual dio como resultado valores moderados, ya que la abundancia de las especies de bivalvos no es homogénea a lo largo de la zona intermareal rocosa en ambas localidades del presente estudio.

Garcés (2011), describió valores similares del índice de diversidad de Shannon para especies de bivalvos en sustratos rocosos en Acapulco (2.41 bits/individuo) y mencionó que éstos varían según el tipo de sustrato, ya que la riqueza de especies fue mayor en sustratos arenosos que en sustratos rocosos, lo que ocasiona que en estos sustratos, los valores disminuyan.

En el nMDS se observan 3 agrupaciones: una con muestras de la localidad de playa El Palmar, otra con estaciones de playa Las Gatas y la tercera con estaciones de ambas localidades, donde solamente hubo presencia de la especie Brachidontes adamsianus (tabla 1). Dicha especie es una de las más comunes en varias localidades de Jalisco, Oaxaca y Guerrero, por lo que sus valores de abundancia en Ixtapa-Zihuatanejo concuerdan con lo encontrado en otras contribuciones (Castro-Mondragón et al., 2016; Galeana-Rebolledo et al., 2012; Garcés, 2011; Holguín-Quiñones y González-Pedraza 1989; Landa-Jaime et al., 2013; López-Rojas et al., 2017; Torreblanca-Ramírez et al., 2012).

Asociación entre grupos morfofuncionales y moluscos. Los bivalvos recolectados tuvieron una mayor presencia en el grupo morfofuncional de algas filamentosas en ambas localidades, este mismo resultado se observó en muelle municipal en Zihuatanejo por Guzmán (2022) para los bivalvos asociados a macroalgas. Lo anterior puede verse favorecido porque este grupo de algas proporciona refugio ante el oleaje y la desecación, de igual forma tienen un papel fundamental en el asentamiento de las larvas de bivalvos, ya que son utilizadas como un sitio para evitar la competencia entre estadios juveniles y adultos (Dobretsov, 1999; Seed, 1969a).

Las larvas de la familia Mytilidae son atraídas hacia los filamentos de algas rojas de los géneros Ceramium y Polysiphonia, en donde llevan a cabo su metamorfosis (Seed, 1969a), ambos géneros fueron encontrados en el presente trabajo. Las especies de bivalvos de la familia Mytilidae también se encontraron junto con otros géneros de macroalgas como: Bryopsis, Caulerpa, Chaetomorpha y Cladophora (Chlorophyta), dichos géneros también se han relacionado a comunidades de bivalvos en el Caribe colombiano (Quirós-Rodríguez y Campos, 2013). Los géneros Gayliella, Herposiphonia, Taenioma y Tayloriella (Rhodophyta) también se han observado junto a especies de bivalvos, en particular el género Tayloriella se ha encontrado vinculado a 9 diferentes especies de bivalvos en el muelle municipal en Zihuatanejo (Guzmán, 2022). Aun cuando es poco lo que se conoce sobre las asociaciones de bivalvos con macroalgas, es posible sugerir que las larvas de las otras familias de bivalvos encontradas en el presente trabajo (Arcidae, Carditidae, Crassatellidae, Chamidae, Myidae y Pteriidae) puedan tener un comportamiento similar al de Mytilidae y, por esta razón, se podría reconocer su presencia en las comunidades de macroalgas.

La asociación de las especies de bivalvos (B. adamsianus, B. semilaevis, L. aristatus, M. capax, M. edulis, A. rostae, I. janus, C. affinis, C. radiata, C. grayi, C. ecuadoriana, C. coralloides, Mya sp. y S. fragilis) con el grupo de algas filamentosas corticadas fue explicado para especies del género Gigartina porSeed (1969a), quien afirmó que este tipo de algas proporciona una mayor protección y partículas de sedimento retenido (alimento), comparado con las algas filamentosas. En el presente trabajo se encontraron especies de los géneros Ceratodyction, Gelidiella, Gelidium, Gracilaria e Hypnea, que podrían tener una relación similar a la dada por las especies de Gigartina, ya que comparten el mismo grupo morfofuncional propuesto por Steneck y Dethier (1994).

Todas las especies de bivalvos se asociaron al grupo de algas calcáreas articuladas. Las especies del género Corallina proporcionan condiciones ambientales adecuadas para los bivalvos como: refugio, captación de partículas y cantidad de CaCO₃ (Seed, 1969a). De esta forma, los bivalvos son capaces de alcanzar tallas mayores en menor tiempo (Seed, 1969a, b). Probablemente, la cantidad de CaCO₃ obtenida del ambiente por las macroalgas, queda disponible para que los bivalvos la aprovechen, una vez que las macroalgas mueren. En el presente trabajo, se observó a las especies de algas calcáreas articuladas de los géneros Amphiroa, Halimeda y Jania asociadas a las especies de bivalvos son muy similares a las del género Corallina y comparten el mismo grupo morfofuncional, por lo que podrían funcionar de la misma forma al ser algas que pueden fijar CaCO₃ y ser una fuente de aprovechamiento para otros organismos; estas algas calcáreas son la fuente principal de carbonatos marinos (Feely et al., 2004).

Cobertura de macroalgas, abundancia de moluscos y sedimento retenido. La cobertura de macroalgas, abundancia de bivalvos y sedimento retenido disminuyeron conforme se aumentó el nivel del intermareal, los valores más elevados para dichas variables se observaron en el nivel del intermareal bajo (fig. 5G-I). Los resultados de la regresión lineal múltiple sustentan lo anterior, ya que se encontró una relación entre la cobertura de macroalgas y abundancia de bivalvos. Lo anterior podría explicarse por las características del terreno (emersión), debido a que, en un ambiente cambiante como la zona intermareal rocosa, las macroalgas quedan expuestas a una mayor radiación solar en el nivel alto. Por lo tanto, al haber menos céspedes algales (Huovinen et al., 2006), la retención de sedimento es baja, estos 2 factores influyen en las necesidades de cada especie de bivalvo, lo que determina la supervivencia de estos moluscos en esta zona (Airoldi, 2003; Rosenberg, 1977).

Renaud et al. (1997) encontraron una relación entre la abundancia de macroalgas con respecto al sedimento, donde la cobertura de macroalgas era sistemáticamente mayor en las zonas con baja cantidad de sedimentos y un aumento notable en la abundancia de macroalgas se observó después de la remoción de sedimentos. Esto se relaciona con lo encontrado en el presente trabajo, donde se observó una menor cobertura de macroalgas asociada a valores elevados de sedimento. Con respecto de los bivalvos, Forster y Zettler (2004) observaron que la biomasa de Mya arenaria Linnaeus, 1758 se redujo con la presencia de sedimentos finos, lo que podría sugerir que en playa Las Gatas existe una menor abundancia de bivalvos (Airoldi, 2003).

El sedimento en las comunidades de macroalgas queda retenido ya que se acumula en los espacios que hay entre los talos; las cantidades de sedimento acumuladas pueden estar determinadas por la complejidad estructural del alga y las condiciones del medio (García, 2009). La variación en la tasa de sedimentación podría ser un factor de alteración y estrés sobre las comunidades de macroalgas, lo que ocasiona que la tasa de crecimiento algal se vea reducida por la falta de incidencia de luz. Por ello, la fauna que ahí se establece, como las especies de bivalvos, también es menor, ya que una gran cantidad de sedimentos retenidos provoca que los organismos puedan sofocarse y el reclutamiento larval se vea reducido (Airoldi, 2003; Rosenberg, 1977). La clase Bivalvia tiene una preferencia por sustratos arenosos, rocosos, lodosos y fangosos, sin embargo, cuando los sitios no tienen corrientes fuertes y la acumulación de sedimento es alta, la abundancia de bivalvos baja drásticamente (Gosling, 2015).

Las larvas de bivalvos se instalan sobre las algas, donde permanecen adheridas a ellas y llevan a cabo su metamorfosis entre 15 y 25 días; en este lapso, los céspedes algales les proporcionan alimento por el sedimento que retienen y protección contra depredadores y factores abióticos. Posteriormente, los bivalvos juveniles migran fuera de las algas ya que han alcanzado una talla óptima o las algas ya no les brindan suficientes recursos, y llegan a otro sitio que será el definitivo para llevar a cabo el resto de su ciclo de vida (Seed, 1969a, 1969b; Suchanek, 1978).

Olafsson (1986) encontró una relación significativa entre la abundancia de bivalvos y el sedimento, estos resultados son consistentes con lo encontrado en el presente estudio en Ixtapa-Zihuatanejo. Sin embargo, Vázquez (2009) en playa Las Gatas no encontró relación entre la riqueza de especies y el sedimento retenido. Por su parte, en playa El Palmar, que es una zona expuesta, una posible explicación para la relación entre abundancia de bivalvos y sedimento es debido al tipo de hábitos de vida de la mayoría de las especies encontradas en el presente trabajo, ya que la mayoría de ellas son epifaunales o infaunales con hábitos filtradores suspensívoros. La erosión del sustrato rocoso por la acción de las olas y por el movimiento de los sedimentos puede generar heterogeneidad ambiental (Airoldi, 2003). La mezcla constante del agua provee mayor humedad, así como nutrientes al suspender los sedimentos, lo que ocasiona una mayor supervivencia de los organismos, lo que explicaría la relación entre abundancia y sedimento (Gama-Kwick et al., 2021).

Las condiciones de playa El Palmar con un oleaje intenso (Morales et al., 2008), no permite la retención de sedimentos, su baja cantidad podría influir en la composición de las macroalgas. Por lo tanto, la composición de los invertebrados al interior de éstos puede ser muy particular en esta localidad (Gama-Kwick et al., 2021), lo que puede incrementar su diversidad (Chemello y Milazzo, 2002; Prathep et al., 2003). En contraste, playa Las Gatas es un sitio protegido con un oleaje de baja energía y alta cantidad de sedimentos, por lo que su comunidad de macroalgas y bivalvos puede ser menos diversa, abundante y compleja. Lo anterior es similar a lo encontrado en trabajos de la zona para otras clases de moluscos como gasterópodos y poliplacóforos (Aguilar, 2017).

Las comunidades de macroalgas pueden funcionar como sustrato de captura de las larvas desde la columna de agua (Aguilar, 2017). Dichos sitios constituyen zonas de crianza en parte o en todo el ciclo de vida de los bivalvos, ya que proporcionan protección y alimento (Seed, 1969a, b), que está determinado por los sedimentos que las macroalgas retienen.

El presente trabajo aporta resultados que contribuyen al conocimiento de la biodiversidad marina, en especial de los bivalvos asociados a la ficoflora del Pacífico tropical mexicano. Los esfuerzos posteriores deberían centrarse en explorar la biodiversidad de los invertebrados sobre sustratos diferentes en el litoral rocoso, posiblemente utilizando metodologías similares a las del presente estudio con la finalidad de generar comparaciones. En los estudios futuros se debe profundizar en el conocimiento de la distribución geográfica de Mytilus edulis y así corroborar los resultados de este trabajo. Se debe promover el estudio de las interacciones ecológicas, ya que muchas especies de moluscos son impactadas por diferentes factores tanto bióticos, como la depredación o la epibiosis (Aguilar-Estrada et al., 2022; García-Ibáñez et al., 2014; Quiroz-González et al., 2020), así como por distintas condiciones abióticas como circulación del agua, temperatura, pH, régimen de mareas y sedimento (López et al., 2017, 2023). Por ello, es fundamental realizar más investigaciones enfocadas hacia el estudio de la variación espacio-temporal de los organismos en periodos distintos a un ciclo anual en diferentes localidades de México, con el objetivo de sentar las bases pertinentes para su conservación y posterior manejo ante el incremento del desarrollo de infraestructuras turísticas/urbanas en zonas como Guerrero, que pueden tener un efecto desfavorable sobre las comunidades intermareales a largo plazo (Zamorano y Leyte-Morales, 2009).

Agradecimientos

Al proyecto DGAPA-PAPIIT, UNAM (IN220714), al Registro Nacional de Pesca y Acuacultura por el permiso para la recolección de material biológico (DF00000208). A Norma López por el préstamo de las instalaciones de la UMDI-Zihuatanejo, a Carlos Candelaria por su apoyo técnico en campo y a Isabel Bieler, por su apoyo en la toma de las fotografías de los ejemplares para
este estudio.

Referencias

Abbott, I. A. (1999). Marine red algae of the Hawaiian Island. Honolulu, Hawaii: Bishop Museum Press.

Abbott, I. A. y Hollenberg, G. J. (1976). Marine algae of California. Stanford, California: Sanford University Press.

Aguilar, L. G. (2017). Estructura comunitaria de los moluscos (gasterópodos pateliformes y poliplacóforos) asociados a los ensambles macroalgales en el intermareal rocoso de Ixtapa-Zihuatanejo, Guerrero, México (Tesis de maestría). Facultad de Ciencias, Universidad Nacional Autónoma de México. Ciudad de México.

Aguilar-Estrada, L. G., Ortigosa, D., Urbano, B. y Reguero, M. (2014). Análisis histórico de los gasterópodos de la laguna arrecifal de Isla Verde, Veracruz, México. Revista Mexicana de Biodiversidad, 85, 502–512. https://doi.org/10.7550/rmb.
33802

Aguilar-Estrada, L. G., Ruiz-Boijseauneau, I. y Rodríguez, D. (2017). Estadios juveniles de las especies de gasterópodos pateliformes y de poliplacóforos (Mollusca) asociados a macroalgas intermareales de Guerrero, México. Revista Mexicana de Biodiversidad, 88, 280–299. https://doi.org/10.1016/j.rmb.2017.03.021

Aguilar-Estrada, L. G., Quiroz-González, N., Ruiz-Boijseauneau, I., Álvarez-Castillo, L. y Rodríguez, D. (2022). Algal epibiont species on Chiton articulatus (Mollusca: Polyplacophora) from a rocky intertidal coast from the Mexican Tropical Pacific. Revista Mexicana de Biodiversidad, 93, e934163. https://doi.org/10.22201/ib.20078706e.2022.93.4163

Aguilera, M. A. (2011). The functional roles of herbivores in the rocky intertidal systems in Chile: A review of food preferences and consumptive effects. Revista Chilena de Historia Natural, 84, 241–261.

Aguillón, A. (2011). Variación espacio-temporal del recluta-
miento en Mollusca y Echinodermata en bahía de La Paz, Baja California Sur, México (Tesis de maestría). Ciudad de México: Centro Interdisciplinario de Ciencias Marinas, Instituto Politécnico Nacional.

Airoldi, L. (2003). The effects of sedimentation on rocky coast assemblages. Oceanography and Marine Biology an Annual Review, 41, 167–171.

Anderson, M., Gorley, R. y Clarke, K. (2008). Permanova+ for PRIMER: guide to software and statistical methods, PRIMER-E, Plymouth.

Bakus, G. J. (2007). Quantitative analysis of marine biological communities: field biology and environmental. Hoboken, New Jersey: John Wiley and Sons Inc.

Baqueiro, E. (1979). Sobre la distribución de Megapitaria aurantiaca (Sowerby), M. squalida (Sowerby) y Dosinia ponderosa (Gray) en relación a la granulometría del sedimento (Bivalvia: Veneridae). Anales del Centro de Ciencias del Mar y Limnología, Universidad Nacional Autónoma de México, 6, 25–32.

Baqueiro, E. y Aldana, D. (2000). A review of reproductive patterns of bivalve mollusks from Mexico. Bulletin of Marine Science, 66, 13–27.

Baqueiro, E. y Aldana, D. (2003). Patrones en la biología poblacional de moluscos de importancia comercial en México. Revista de Biología Tropical, 51, 97–107.

Baqueiro, E. y Masso, J. A. (1988). Variaciones poblacionales y reproducción de dos poblaciones de Chione undatella (Sowerby, 1835) bajo diferentes regímenes de pesca en la bahía de La Paz, BCS México. Ciencia Pesquera, 6, 51–67.

Barba-Marino, F., Flores-Rodríguez, P., Flores-Garza, R., García-Ibáñez, S. y Arana-Salvador, D. G. (2010). Biodiversidad y zonificación de la comunidad de moluscos, que habita el sustrato rocoso en dos sitios con distinta acción del oleaje, en la Isla “La Roqueta”, Acapulco, Guerrero, México. En L. J. Rangel, J. Gamboa, S. L. Arriaga y W. M. Contreras (Eds.), Perspectiva en malacología mexicana (pp. 44–56). Villahermosa: Universidad Juárez Autónoma de Tabasco.

Bartlett, M.S. (1937). Properties of sufficiency and statistical test. Proceedings of the Royal Society A, 160, 268–282. https://doi.org/10.1098/rspa.1937.0109

Baumgartner, T. R. y Christensen, N. (1985). Coupling of the Gulf of California to large-scale interannual climatic variability. Journal of Marine Research, 43, 825–848.

Benedetti-Cecchi, L. (2001). Variability in abundance of algae and invertebrates at different spatial scales on rocky sea shores. Marine Ecology Progress Series, 215, 79–92.

Benedetti-Cecchi, L., Rindi, F., Bertocci, I., Bulleri, F. y Cinelli, F. (2001). Spatial variation in development of epibenthic assemblages in a coastal lagoon. Estuarine, Coastal and
Shelf Science, 52, 659–668. https://doi.org/10.1006/ecss.20
01.0775

Borges, L., Merckelbach, L. M., Sampaio, Í. y Cragg, S. M. (2014). Diversity, environmental requirements, and biogeography of bivalve wood-borers (Teredinidae) in European coastal waters. Frontiers in Zoology, 11, 1–13. https://doi.org/10.1186/1742-9994-11-13

Bouchet, P., Rocroi, J. P., Bieler, R., Carter, J. G. y Coan, E. V. (2010). Nomenclator of bivalve families with a classification of bivalve families. Malacologia, 52, 1–185. https://doi.org/
10.4002/040.052.0201

Cadena-Cárdenas, L., Méndez-Rodríguez, L., Zenteno-Savín, T., García-Hernández, J. y Acosta-Vargas, B. (2009). Heavy metal levels in marine mollusks from areas with, or without, mining activities along the Gulf of California, Mexico. Archives of Environmental Contamination and Toxicology, 57, 96–102. https://doi.org/10.1007/s00244-008-9236-0

Castro-Mondragón, H., Flores-Garza, R., Valdez-González, A., Flores-Rodríguez, P., García-Ibáñez, S. y Rosas-Acevedo, J. L. (2016). Diversidad, especies de mayor importancia y composición de tallas de los moluscos en la pesca ribereña en Acapulco, Guerrero, México. Acta Universitaria, 26, 24–34. https://doi.org/10.15174/au.2016.1025.

Ceccherelli, V. U. y Rossi, R. (1984). Settlement, growth and production of the mussel Mytilus galloprovincialis. Marine Ecology Progress Series, 16, 173–184.

Cerros-Cornelio, J. C., Flores-Garza, R., Landa-Jaime, V., García-Ibáñez, S., Rosas-Guerrero, V., Flores-Rodríguez, P. et al. (2021). Composición de especies e ingreso económico por la pesca ribereña de moluscos en la Costa Grande de Guerrero México. Revista Bio Ciencias, 8, e1054. https://doi.org/10.15741/revbio.08.e1054

Chemello, R. y Milazzo, M. (2002). Effect of algal architecture on associated fauna: some evidence from phytal molluscs. Marine Biology, 140, 981–990. https://doi.org/10.1007/s00
227–002–0777–x

Cisneros, F. (2016). Estructura comunitaria de gasterópodos asociados a macroalgas en el litoral rocoso de Ixtapa-Zihuatanejo, Guerrero, México (Tesis). Facultad de Ciencias,
Universidad Nacional Autónoma de México. Ciudad de México.

Clarke K. R. y Gorley, R. N. (2006) PRIMER v6: User Manual/Tutorial. PRIMER-E Ltd, Plymouth, Reino Unido.

Coan, E. V. (1999). The eastern Pacific species of Sphenia (Bivalvia: Myidae). Nautilus-Sanibel, 113, 103–120. https://doi.org/10.5962/bhl.part.2019

Coan, E. V., Scott, P. V. y Bernard, F. R. (2000). Bivalve seashells of western North America. Santa Bárbara, California: Santa Barbara Museum of Natural History Monographs.

Coan, E. V. y Valentich-Scott, P. (2006). Marine Bivalvia. En C. F. Sturm, T. A. Pearce y A. Valdés (Eds.), The mollusks: a guide to their study, collection, and preservation (pp. 339–349). Pittsburgh, Pennsylvania: American Malacological Society.

Coan E.V. y Valentich-Scott P. (2012). Bivalve seashells of tropical West America. Marine bivalve mollusks from Baja California to Peru. Santa Barbara: Santa Barbara Museum of Natural History Monographs.

Colman, J. (1940). On the faunas inhabiting intertidal seaweeds. Journal of the Marine Biological Association of the United Kingdom, 24, 129–183. https://doi.org/10.1017/S0025315400054503

Dawson, E. Y. (1949). Resultados preliminares de un reco-
nocimiento de las algas marinas de la costa del pacífico de México. Revista de la Sociedad Mexicana de Historia Natural, 9, 215–255.

Dawson, E. Y. (1953). Marine red algae of Pacific Mexico. Part 1. Bangiales to Corallinaceae subf. Corallinoidae. Allan Hancock Pacific Expedition, 17, 1–239.

Dawson, E. Y. (1954). Marine red algae of Pacific Mexico. Part 2. Cryptonemiales (cont.). Allan Hancock Pacific Expedition, 17, 241–397.

Dawson, E. Y. (1960). Marine red algae of Pacific Mexico. Part 3. Cryptonemiales, Corallinaceae subf. Melobesioidae. Pacific Naturalist, 2, 3–125.

Dawson, E. Y. (1961). Marine red algae of Pacific Mexico. Part 4. Gigartinales. Pacific Naturalist, 2, 191–343.

Dawson, E. Y. (1963). Marine red algae of Pacific Mexico part 8. Ceramiales: Dasyaceae, Rhodomelaceae. Nova Hedwigia, 6, 401–481.

Dawson, E. Y. y Beaudette, P. T. (1959). Field notes from the 1959 Eastern Pacific Cruise of the Stella Polaris. Pacific Naturalist, 1, 1–24.

Díaz, A. y Ortlieb, L. (1993). El fenómeno “El Niño” y los moluscos de la costa peruana. Bulletin de l’Institut Français d’Études Andines, 22, 159–177.

Dobretsov, S. V. (1999). Effects of macroalgae and biofilm on settlement of blue mussel (Mytilus edulis l.) larvae, Bio-
fouling, 14, 153–165. https://doi.org/10.1080/089270199093
78406

Esqueda-González, M., Ríos-Jara, E., Galván-Villa, C. M. y Rodríguez-Zaragoza, F. A. (2014). Species composition, richness, and distribution of marine bivalve molluscs in Bahía de Mazatlán, México. Zookeys, 399, 43–69. https://doi.org/10.3897/zookeys.399.6256

Feely, R. A., Sabine, C. L., Lee, K., Berelson, W., Kleypas, J., Fabry, V. J. et al. (2004). Impact of anthropogenic CO₂ on the CaCO₃ system in the oceans. Science, 305, 362–366. https://doi.org/10.1126/science.1097329

Fitch, J. E. (1953). Common marine bivalves of California. Sacramento, California: California Department of Fish and Game.

Flores, P. (2004). Estructura de la comunidad de moluscos del mesolitoral superior en las playas de facie rocosa del estado de Guerrero, México (Tesis doctoral). Facultad de Ciencias Biológicas, Universidad Autónoma de Nuevo León. San Nicolás de los Garza.

Flores-Garza, R., Flores-Rodríguez, P., García-Ibáñez, S. y Valdés-González, A. (2007). Demografía del caracol Plicopurpura pansa (Neotaenioglossa: Muricidae) y constitución de la comunidad malacológica asociada en Guerrero, México. Revista de Biología Tropical, 55, 867–878.

Flores-Garza, R., García-Ibáñez, S., Flores-Rodríguez, P., Torreblanca-Ramírez, C., Galeana-Rebolledo, L., Valdés-González, A. et al. (2012). Commercially important marine mollusks for human consumption in Acapulco, México. Natural Resources, 3, 11–17. http://dx.doi.org/10.4236/nr.20
12.31003

Flores-Garza, R., López-Rojas, V., Flores-Rodríguez, P. y Ramírez, C. T. (2014). Diversity, distribution, and composition of the Bivalvia class on the rocky intertidal zone of marine priority region 32, Mexico. Open Journal of Ecology, 4, 961–973 https://doi.org/10.4236/oje.2014.415080

Flores-Garza, R., Torreblanca-Ramírez, C., Flores-Rodríguez, P., García-Ibáñez, S. y Galeana-Rebolledo, L. (2010). Riqueza y análisis de la comunidad malacológica en el mesolitoral rocoso de la playa Tlacopanocha, Acapulco, Guerrero. En L. J. Rangel, J. Gamboa, S. L. Arriaga y W. M. Contreras (Eds.), Perspectiva en malacología mexicana (pp. 125–138). Villahermosa: Universidad Juárez Autónoma de Tabasco.

Flores-Garza, R., Torreblanca-Ramírez, C., Flores-Rodríguez, P., García-Ibáñez, S., Galeana-Rebolledo, L., Valdés-González, A. et al. (2011). Mollusc community from a rocky intertidal zone in Acapulco, Mexico. Biodiversity, 12, 144–153. https://doi.org/10.1080/14888386.2011.625520

Flores-Rodríguez, P., Flores-Garza, R., García-Ibáñez, S. y Valdés-González, A. (2003). Riqueza y diversidad de la malacofauna del mesolitoral rocoso de la Isla la Roqueta, Acapulco, Guerrero, México. Ciencia, Revista de Investigación Científica, 11, 5–14.

Flores-Rodríguez, P., Flores-Garza, R., García-Ibáñez, S. y Valdés-González, A. (2007). Variación en la diversidad malacológica del mesolitoral rocoso en Playa Troncones, La Unión, Guerrero, México. Revista Mexicana de Biodiversidad, 78, 33–40. https://doi.org/10.22201/ib.20078706e.2007.002.298

Flores-Rodríguez, P., Flores-Garza, R., García-Ibáñez, S., Valdés-González, A., Violante-González, J., Santiago, C. E. et al. (2012). Mollusk species richness on the rocky shores of the state of Guerrero, Mexico, as affected by rains and their geographical distribution. Natural Resources, 3, 248–260. http://dx.doi.org/10.4236/nr.2012.34032

Forster, S. y Zettler, M. L. (2004). The capacity of the filter-feeding bivalve Mya arenaria L. to affect water transport in sandy beds. Marine Biology, 144, 1183–1189. https://doi.org/10.1007/s00227-003-1278-2

Galeana-Rebolledo, L., Flores-Garza, R., Torreblanca-Ramírez, C., García-Ibáñez, S., Flores-Rodríguez, P. y López-Rojas, V. I. (2012). Biocenosis de Bivalvia y Polyplacophora del intermareal rocoso en playa Tlacopanocha, Acapulco, Guerrero, México. Latin American Journal of Aquatic Research, 40, 943–954. http://dx.doi.org/10.3856/vol40-
issue4-fulltext-11

Galeana-Rebolledo, L., Flores-Garza, R., Violante-González, J., Flores-Rodríguez, P., García-Ibáñez, S., Landa-Jaime, V. et al. (2018). Socioeconomic aspects for coastal mollusk commercial fishing in Costa Chica, Guerrero, México. Natural Resources, 9, 229–241. https://doi.org/10.4236/nr.
2018.96015

García-Domínguez, F. A., Haro-Hernández, D., García-Cuellar, Á., Villalejo-Fuerte, M. y Rodríguez-Astudillo, S. (2008). Ciclo reproductivo de Anadara tuberculosa (Sowerby, 1833) (Arcidae) en Bahía Magdalena, México. Revista de Biología Marina y Oceanografía, 43, 143–152. http://dx.doi.org/10.4067/S0718-19572008000100015.

Gama, N. (2019). Aspectos ecológicos de la comunidad de gasterópodos (Mollusca: Gastropoda) asociados a ensambles macroalgales en el intermareal rocoso de Ixtapa Zihuatanejo, Guerrero, México (Tesis). Facultad de Ciencias, Universidad Nacional Autónoma de México. Ciudad de México.

Gama-Kwick, N., Aguilar-Estrada, L. G., Quiroz-González, N. y Ruiz-Boijseuneau, I. (2021). Nuevos registros de gasterópodos (Mollusca) asociados a macroalgas inter-
mareales de Guerrero, México. Revista Mexicana de Biodiversidad, 92, 1–18. http://dx.doi.org/10.22201/ib.200
78706e.2021.92.3441

Garcés, J. L. (2011). Micromoluscos bivalvos de la Bahía de Acapulco Guerrero, México: composición específica y diversidad (Tesis). Facultad de Ciencias, Universidad Nacional Autónoma de México. Ciudad de México.

García, M. (2009). Estructura comunitaria de la fauna asociada a algas submareales en tres sitios contrastantes en Zihuatanejo, Guerrero (Tesis de maestría). Facultad de Ciencias, Universidad Nacional Autónoma de México. Ciudad de México.

García-Cubas, A. (1981). Moluscos de un sistema lagunar tropical en el sur del Golfo de México (Laguna de Términos, Campeche). Publicaciones especiales-Instituto de Ciencias del Mar y Limnología, 5, 1–379.

García-Cubas, A. y Reguero, M. (2007). Catálogo ilustrado de moluscos bivalvos del Golfo de México y Mar Caribe. Ciudad de México: Instituto de Ciencias del Mar y Limnología, Universidad Nacional Autónoma de México.

García-Ibáñez, S., Flores-Rodríguez, P., Navarro, J. T. N., Garza, R. F. y Moreno, I. E. B. (2014). Respuesta del carnívoro Plicopurpura pansa (Mollusca: Gastropoda) y el herbívoro Chiton articulatus (Mollusca: Polyplacophora) a factores ambientales en Acapulco, México. CienciaUAT, 8, 11–21.

García-Robledo, E., Corzo, A., van Bergeijk, S. y Yúfera, M. (2008). Impacto de las acumulaciones de macroalgas en la comunidad biológica intermareal. Revista de la Sociedad Gaditana de Historia Natural, 8, 115–137.

Gibbons, M. J. (1988). The impact of sediment accumulation, relative habitat complexity and elevation on rocky shore meiofauna. Journal of Experimental Marine Biology Ecology, 122, 225–241. https://doi.org/10.1016/0022-0981(88)90125-6

Gosling, E. (2015). Marine bivalve molluscs. Chichester, West Sussex: John Wiley y Sons.

Guiry, M. D. y Guiry, G. M. (2024). AlgaeBase. National University of Ireland, Galway. Recuperado el 12 enero, 2024 de: http://www.algaebase.org

Guzmán, A. R. (2022). Estructura comunitaria de bivalvos asociados a macroalgas en el litoral rocoso de Muelle Municipal, Zihuatanejo, Guerrero, México (Tesis). Facultad de Ciencias. Universidad Nacional Autónoma de México. Ciudad de México

Hernández-Moreno, E. P., Romo-Piñera, A. K., Fernández-Rivera Melo, F. J., Aguilar-Cruz, C. A., Reyes-Bonilla, H. y López-Vivas, J. M. (2020). Reproductive Study of Megapitaria aurantiaca (Sowerby, 1831) (Bivalvia: Veneridae) in Puerto Libertad, Sonora, Mexico. Journal of Shellfish Research, 39, 441–447.

Holguín-Quiñones, O. E. y González-Pedraza, A. C. (1989). Moluscos de la franja costera del estado de Oaxaca, México. México D.F.: Dirección de Bibliotecas y Publicaciones, Instituto Politécnico Nacional.

Horton, T., Kroh, A., Ahyong, S., Bailly, N., Boyko, C. B. y Brandão, S. N. (2024. World Register of Marine Species. Recuperado el 14 de marzo, 2024: http://www.marinespecies.org

Huovinen, P., Gómez, I. y Lovengreen, C. (2006). A five-year study of solar ultraviolet radiation in Southern Chile (39° S): Potential impact on physiology of coastal marine algae? Photochemistry and Photobiology, 82, 515–522. https://doi.org/10.1562/2005-07-05-RA-601

IMTA (Instituto Mexicano de Tecnología del Agua). (2010). Estudio de clasificación de la Bahía de Ixtapa-Zihuata-
nejo. No. FON-CNA-2004-02-016. Informe final. Comisión Nacional del Agua.

Jiménez-Valverde, A. y Hortal, J. (2003). Las curvas de acumulación de especies y la necesidad de evaluar la calidad de los inventarios biológicos. Revista Ibérica de Aracnología, 8, 151–161.

Jorgensen, C. B. (1996). Bivalve filter feeding revisited. Marine Ecology Progress Series, 142, 287–302.

Jover-Capote, A. y Diez, Y. L. (2017). Abundancia de moluscos en mantos de macroalgas del mesolitoral rocoso en la costa suroriental de Cuba. Amici Molluscarum, 25, 27–43.

Keen, A. M. (1971). Sea shells of Tropical West America Marine mollusks from Baja California to Peru. Standford, California: Stanford University Press.

Kuk-Dzul, J. G., Padilla, J. G., Torreblanca, C., Flores, R., Flores, P. y Muñiz, X. I. (2019). Structure of molluscan communities in shallow subtidal rocky bottoms of Acapulco, Mexico. Turkish Journal of Zoology, 43, 465–479.

Landa-Jaime, V. y Arciniega-Flores, J. (1998). Macromoluscos bentónicos de fondos blandos de la plataforma continental de Jalisco y Colima, México. Ciencias Marinas, 24, 155–l67.

Landa-Jaime, V., Michel-Morfín, E., Arciniega-Flores, J., Castillo-Vargasmachuca, S. y Saucedo-Lozano, M. (2013). Moluscos asociados al arrecife coralino de Tenacatita, Jalisco, en el Pacífico central mexicano. Revista Mexicana de Biodiversidad, 84, 1121–1136. https://doi.org/10.7550/rmb.32994

Lee, R. E. (2008). Basic characteristics of the algae. En R. E. Lee (Ed.), Phycology (pp. 3-29). Nueva York: Cambridge University Press.

Lesser, H. (1984). Prospección sistemática y ecológica de los moluscos bentónicos de la plataforma continental del estado de Guerrero, México (Tesis). Facultad de Ciencias, Universidad Nacional Autónoma de México. Ciudad de México.

Levene, H. (1960) Robust tests for equality of variances. En I. Olkin (Ed.), Contributions to probability and statistics (278–292). Palo Alto, California: Standford University Press.

López, N. A., (1993). Caracterización de la ficoflora sublitoral de Acapulco y Zihuatanejo, Gro. (Tesis). Facultad de Cien-
cias, Universidad Nacional Autónoma de México. Ciudad de México.

López, N., Candelaria, C. y Ramírez-García, P. (2023). Assessment of macroalgae coverage in a scarcely studied deep rocky reef in the tropical eastern Mexican Pacific. Latin American Journal of Aquatic Research, 51, 23–33. https://doi.org/10.3856/vol51-issue1-fulltext-2920

López, N., Candelaria, C., Ramírez-García, P. y Rodríguez, D. (2017). The structure of tropical turf-forming algae assemblages. Zihuatanejo Bay, México. Latin American Journal of Aquatic Research, 45, 329–340 http://dx.doi.org/10.3856/vol45-issue2-fulltext-9

López-Rojas, V. I., Flores-Garza, R., Flores-Rodríguez, P., Torreblanca-Ramírez, C. y García-Ibáñez, S. (2017). La clase Bivalvia en sitios rocosos de las Regiones Marinas Prioritarias en Guerrero, México: riqueza de especies, abundancia y distribución. Hidrobiológica, 27, 69–86.

Lozada, O. (2010). Actualización sistemática de los bivalvos de la Colección Malacológica Dr. Antonio García-Cubas del Instituto de Ciencias del Mar y Limnología (Tesis). Facultad de Ciencias, Universidad Nacional Autónoma de México. Ciudad de México.

Lozada, O. (2015). Presencia de metales pesados en Isogno-
mon alatus (Gmelin, 1791) de la laguna Tampamachoco, Veracruz (Tesis de maestría). Facultad de Ciencias Bio-
lógicas y Agropecuarias, Universidad Veracruzana. Tuxpan, Veracruz.

Magurran, A. E. (2004). Measuring biological diversity. Malden, Massachusetts: Blackwell Publishing.

Morales, R., Vélez, H., Mejía, A. Ramírez, I., Izurierta, J. y Saldaña, P. (2008). Hidrodinámica de la Bahía de Zihuata-
nejo. XXIII Congreso Latinoamericano de Hidráulica, 2 al 6 de septiembre de 2008, Cartagena de Indias, Colombia.

Moreno, C. A. (1995). Macroalgae as a refuge from predation for recruits of the mussel Choromytilus chorus (Molina, 1782) in southern Chile. Journal of Experimental Marine Biology and Ecology, 191, 181–193.

Moreno, C. E. (2001). Métodos para medir la biodiversidad. Zaragoza, España: M&T–Manuales y Tesis SEA.

Oksanen, J., Guillaume-Blanchet, F., Friendly, M., Kindt, R. Legendre, P., McGlinn, D. et al. (2019). Vegan: community ecology package. R package versión 2.5-6. https://CRAN.R-project.org/package=vegan

Olabarria, C. y Chapman, M. G. (2001). Comparison of patterns of spatial variation of microgastropods between two contrasting intertidal habitats. Marine Ecology Progress Series, 220, 201–211. https://doi.org/doi:10.3354/meps220201

Olafsson, E. B. (1986). Density dependence in suspension-feeding and deposit-feeding populations of the bivalve Macoma balthica: a field experiment. Journal of Animal Ecology, 55, 517–526. https://doi.org/10.2307/4735

Oyarzún, P. A., Toro, J. E., Jaramillo, R., Guiñez, R., Briones, C. y Astorga, M. (2011). Ciclo gonadal del chorito Mytilus chilensis (Bivalvia: Mytilidae) en dos localidades del sur de Chile. Latin American Journal of Aquatic Research, 39, 512–525.

Paredes, C., Tarazona, J., Canahuire, E., Romero, L., Cornejo, O. y Cardoso, F. (1998). Presencia de moluscos tropicales de la provincia panameña en la costa central del Perú y su relación con los eventos “El Niño”. Revista Peruana de Biología, 5, 123–128. https://doi.org/10.15381/rpb.v5i2.8330

Pérez, L. M. (2013). Descripción espacio-temporal de la temperatura superficial del mar en el Pacífico Sur Mexicano de 1996-2009 (Tesis). Universidad del Mar. Puerto Ángel, Oaxaca.

Prathep, A., Marrs, R. H. y Norton, T. A. (2003). Spatial and temporal variations in sediment accumulation in an algal turf and their impact on associated fauna. Marine Biology, 142, 381–390. https://doi.org/10.1007/s00227–002–0940–4

Quirós-Rodríguez, J. y Campos, N. H. (2013). Moluscos asociados a ensamblajes macroalgales en el litoral rocoso de Córdoba, Caribe colombiano. Boletín de Investigaciones Marinas y Costeras, 42, 101–120.

Quiroz-González, N., Aguilar-Estrada, L. G., Ruiz-Boijseauneau, I. y Rodríguez, D. (2020). Biodiversidad de algas epizoicas en el Pacífico tropical mexicano. Acta Botanica Mexicana, 127, 1-22. https://doi.org/10.21829/abm127.2020.1645

R Core Team (2020). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/

Reguero, M., y García-Cubas, A. (1989). Moluscos de la plataforma continental de Nayarit: sistemática y ecología (cuatro campañas oceanográficas). Anales del Instituto de Ciencias del Mar y Limnología, 16, 33–58.

Renaud, P. E., Riggs, S. R., Ambrose Jr, W. G., Schmid, K. y Snyder, S. W. (1997). Biological-geological interactions: storm effects on macroalgal communities mediated by sediment characteristics and distribution. Continental Shelf Research, 17, 37–56. https://doi.org/10.1016/0278-43
43(96)00019-2

Ríos-Jara, E., López-Uriarte, E., Pérez-Peña, M. y Juárez-Carrillo, E. (2003). Nuevos registros de escafópodos para las costas de Jalisco y Colima, México. Hidrobiológica, 13, 167–170.

Rodríguez, D., López, N. y González-González, J. (2008). Gelidiales (Rhodophyta) en las costas del Pacífico mexicano con énfasis en las especies tropicales. En A. Sentíes y K. M. Dreckmann (Eds.), Monografías ficológicas, Vol. 3 (pp. 27–74). Ciudad de México: Universidad Autónoma Metropolitana.

Román-Contreras, R., Cruz-Abrego, F. M. e Ibáñez-Aguirre, A. L. (1991). Observaciones ecológicas de los moluscos de la zona intermareal rocosa de la Bahía de Chamela, Jalisco, México. Anales del Instituto de Biología, Universidad Nacional Autónoma de México, Serie Zoología, 62, 17–31.

Rosenberg, R. (1977). Effects of dredging operations on estuarine benthic macrofauna. Marine Ecology Progress Series, 62, 185–202.

Salazar-Vallejo, S. I. y González, E. (1990). Ecología costera en la región de La Mancha, Veracruz. La Ciencia y el Hombre, 6,101–120.

Salcedo-Martínez, S., Green, G., Gamboa-Contreras, A. y Gómez,
P. (1988). Inventario de macroalgas y macroinvertebrados bénticos, presentes en áreas rocosas de la región de Zihuatanejo, Guerrero, México. Anales del Instituto de Ciencias del Mar y Limnología, Universidad Nacional Autónoma de México, 15, 73–95.

Sánchez, M. (2014). Base de datos de los moluscos presentes en la colección biológica de la Secretaría de Marina-Armada de México (Tesis). Facultad de Ciencias, Universidad Nacional Autónoma de México. Ciudad de México.

Santos, L., Souza, J., Lima, S. y Guimarães, C. (2020). Diversity of bivalve mollusks associated with macroalgae on the continental shelf in the states of Alagoas, Sergipe and Bahia, Northeastern Brazil. Zoological Studies, 59, e58. https://doi.org/10.6620/ZS.2020.59-58

Schulien, J. A., Adams, J. y Felis, J. (2020). Pacific continental shelf environmental assessment (PaCSEA): characterization of seasonal water masses within the Northern California current system using airborne remote sensing off Northern California, Oregon, and Washington, 2011–2012. Camarillo, California: US Department of the Interior, Bureau of Ocean Energy Management/ Pacific OCS Region/ OCS Study BOEM.

Seed, R. (1969a). The ecology of Mytilus edulis L. (Lamellibranchiata) on exposed rocky shores. Oecologia, 3, 277–316. https://doi.org/10.1007/BF00390380

Seed, R. (1969b). The ecology of Mytilus edulis L. (Lamelli-
branchiata) on exposed rocky shores: II. Growth and mortality. Oecologia, 3, 317–350. https://www.jstor.org/stable
/4214550

Shaw, W. N., Hassler, T. J. y Moran, D. P. (1988). California sea mussel and bay mussel: species profiles. Life histories and environmental requirements of coastal fishes and invertebrates (Pacific Southwest). Washington D.C.: Fish and Wildlife Service.

Sibaja-Cordero, J. A. y Vargas-Zamora, J. A. (2006). Zonación vertical de epifauna y algas en litorales rocosos del Golfo de Nicoya, Costa Rica. Revista de Biología Tropical, 54, 49–67.

Siegel, S. (1990). Estadística no paramétrica: aplicada a las ciencias de la conducta. Ciudad de México: Trillas.

Stella, C., Vijayalakshmi, S. y Siva, J. (2010). Two new records of bivalve species of Mytilidae family from Palk Bay and Mandapam area-south-east coast of India. Global Journal of Environmental Research, 4, 40–42.

Steneck, R. S. y Dethier, M. N. (1994). A functional group approach to the structure of algal-dominated communities. Oikos, 69, 476–498. https://doi.org/10.2307/3545860

Steneck, R. S. y Watling, L. (1982). Feeding capabilities and limitation of herbivorous molluscs: a functional group approach. Marine Biology, 68, 299–319. https://doi.org/10.
1007/BF00409596

Suchanek, T. H. (1978). The ecology of Mytilus edulis L. in exposed rocky intertidal communities. Journal of Experimental Marine Biology and Ecology, 31, 105–120. https://doi.org/10.1016/0022-0981(78)90139-9

Taylor, W. R. (1945). Pacific marine algae of the Allan Hancock Expeditions to the Galapagos Islands. Allan Hancock Pacific Expeditions, 12, 1–528.

Torreblanca, C. (2010). Análisis de la diversidad y estructura de la comunidad de moluscos del mesolitoral rocoso de Acapulco, Guerrero (Tesis).Unidad Académica de Ecología Marina, Universidad Autónoma de Guerrero. Acapulco de Juárez.

Torreblanca-Ramírez, C., Flores-Garza, R., Flores-Rodríguez, P., García-Ibáñez, S. y Galeana-Rebolledo, L. (2012). Riqueza, composición y diversidad de la comunidad de moluscos asociada al sustrato rocoso intermareal de playa Parque de la Reina, Acapulco, México. Revista de Biología Marina y Oceanografía, 47, 283–294. http://dx.doi.org/10.4067/S0718-19572012000200010

Torroglosa, M. E. (2015). Biología reproductiva y crecimiento de Brachidontes rodriguezii (d´Orbign 1846) en sustratos duros artificiales en playas arenosas de la provincia de Buenos Aires (Tesis doctoral). Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires. Argentina.

Trasviña, A. y Andrade, C. A. (2002). La circulación costera del Pacifico Tropical Oriental, con énfasis en la Alberca Cálida Mexicana (ACM). En Corcas (Eds.), Circulación oceánica y climatología tropical en México y Colombia (pp. 9–37). Bogotá: Consejo Nacional de Acreditación-Ministerio de Educación.

UNAM (Universidad Nacional Autónoma de México) (2013). Laguna Las Salinas Zihuatanejo. Proyecto ejecutivo UNAM 1989-1994 basado en plano base FONATUR-SCT. Universidad Nacional Autónoma de México, México. Recuperado el 14 enero, 2024 de https://mmacreactive.wordpress.com/

Urbano, B. (2004). Estructura comunitaria de gasterópodos de Zihuatanejo, Guerrero (Tesis). Facultad de Ciencias, Universidad Nacional Autónoma de México. Ciudad de México.

Valdés-González, A., Flores-Rodríguez, P., Flores-Garza, R. y García-Ibáñez, S. (2004). Molluscan communities of the rocky intertidal zone at two sites with different wave action on Isla La Roqueta, Acapulco, Guerrero, Mexico. Journal of Shellfish Research, 23, 875–880.

Vassallo, A., Dávila, Y., Luviano, N., Deneb-Amozurrutia, S., Vital, X. G., Conejeros, C. A. et al. (2014). Inventario de invertebrados de la zona rocosa intermareal de Montepío, Veracruz, México. Revista Mexicana de Biodiversidad, 85, 349–362. https://doi.org/10.7550/rmb.42628

Vázquez, P. (2009). Efecto del sedimento retenido en la estructura de los ensambles algales submareales (Tesis). Facultad de Ciencias, Universidad Nacional Autónoma de México. Ciudad de México.

Vega, C., Olabarria, C. y Carballo, J. L. (2008). Variación espacio-temporal de moluscos y macroalgas en sustratos rocosos intermareales en la bahía de Mazatlán. Ciencia y Mar, 34, 3–16.

Velez, J. y Zeballos, J. (1985). Ampliación de la distribución de algunos peces e invertebrados durante el fenómeno “El Niño” 1982–1983. En W. Arntz, A. Landa y J. Tarazona (Eds.), “El Niño” su impacto en la fauna marina (pp. 173– 180). Callao, Perú: Boletín del Instituto del Mar del Perú.

Villalpando, E. (1986). Diversidad y zonación de moluscos de facie rocosa Isla Roqueta, Acapulco, Gro. (Tesis). Facultad de Ciencias, Universidad Nacional Autónoma de México. Ciudad de México.

Villegas-Maldonado, S., Neri-García, E., Flores-Garza, R., García-Ibáñez, S., Flores Rodríguez, P. y Arana-Salvador, D. G. (2007). Datos preliminares de la diversidad de moluscos para el consumo humano que se expenden en Acapulco, Guerrero. En E. Ríos-Jara, M. C. Esqueda-González y C. M Galván-Villa (Eds.), Estudios sobre la malacología y conquiliología en México (pp. 57-59). Guadalajara: Universidad de Guadalajara.

Ward, J. E., Sanford, L. P., Newell, R. I. E. y MacDonald, B. A. (1998). A new explanation of particle capture in suspension-feeding bivalve molluscs. Limnology and Oceanography, 43, 741–752.

Wyrtki, K. 1966. Oceanography of the eastern equatorial Pacific Ocean. Oceanography and marine biology: an annual review, 4, 33–68.

Yang, J. L., Satuito, C. G., Bao, W. Y. y Kitamura, H. (2007). Larval settlement and metamorphosis of the mussel Mytilus galloprovincialis on different macroalgae. Marine Biology, 152, 1121–1132. https://doi.org/10.1007/s00227-007-0759-0

Zamorano, P. y Leyte-Morales, G. E. (2009). Equinodermos asociados a formaciones arrecifales en Zihuatanejo y Acapulco, Guerrero, México. Boletín de Investigaciones Marinas y Costeras-INVEMAR, 38, 7–28.