Introduction

Ghost crabs are among the most common burrowing organisms on sandy beaches from the tropics to temperate latitudes (Lucrezi and Schlacher 2014; Jones 1972; Ma et al. 2019; Wong et al. 2012). These crabs use self-constructed burrows for shelter (Christoffers 1986; Kristensen et al. 2012), mating (Linsenmair 1967; Hughes 1973), egg development (Haley 1973), molting (Christoffers 1986), sex-specific signaling (Linsenmair 1967; Hartnoll 1969) and feeding (Crane 1941; Robertson and Pfeiffer 1982; Trott 1987). In addition, the tunnel diameter and depth is sufficiently small to prevent deformation and collapse of their burrows (Shinoda et al. 2019). With the exception of five or six zoea stages, which are planktonic (Díaz and Costlow 1972; McDermott 2009; Jiang et al. 2014), most ghost crabs spend their life in their self-constructed burrow. In this regard, digging activity of ghost crabs has important bioturbation effects on sediment characteristics, biogeochemistry and substrate (Schlacher et al. 2011). Thus, the burrowing behaviour of these crabs leads to improving aeration, substrate oxidation–reduction potential, and increases the complexity of the substrate (Chan et al. 2006; Lucrezi and Schlacher 2014).

Many studies have done on burrows of Ocypode crabs in recent years which have investigated morphological characteristics, distribution and architecture of their burrows within the specified time of the year (e.g. on O. ryderi (Vannini 1980), on O. saratan (Ali-Adnan 1985; Eshky 1985), on O. gaudichaudii (Schober and Christy 1993), on O. ceratophthalama (Chakrabarti 1981; Brooke 1981; Chan et al. 2006; Seike and Nara 2008; Lim et al. 2011), on O. sinensis (Seike and Nara 2008), on O. quadrata (McDermott 2009; Branco et al. 2010; Correa et al. 2014), on O. macrocera (Haque and Choudhury 2014)).

To date, two species of ghost crab have previously been reported on Iranian sandy beaches, namely, O. rotundata (Naderloo et al. 2015( and O. sinensis)Naderi et al. 2018a). Our target species, O. rotundata is commonly found in sandy beaches along the East coast of the Arabian Peninsula (Oman) to the west coast of India (Mumbai), including the Persian Gulf (Iran, Kuwait, Saudi Arabia, Bahrain, U.A.E.) (Naderloo and Türkay 2012; Sakai and Türkay 2013).This crab species has a life span of four years and can reach sizes of up to 60.02 mm carapace width (Naderi et al. 2019), with the main breeding season from March to October, e.g., in the Qeshm Island (Naderi et al. 2018b). Despite this considerable distribution range, few works have been conducted on the biology of O. rotundata (in Chabahar Bay, Najafi 2014; in Jask harbour, Zahedi et al. 2014, in Qeshm Island, Naderi et al. 2018a and 2019).

Detailed studies on the burrowing activity of these species are scarce. Burrows constructed by O. rotundata included Y-shaped, J-shaped, single Shaft, anchor, and spiral shapes (Zahedi 2014; Naderi and Pishehvarzad 2019). According to mentioned contents and dependence of O. rotundata on sandy beaches, this study will focus on the influence of environmental conditions (organic matter content of sediment, substrate moisture percentage, sand granulometry) and temporal variability (reproduction and non-reproduction period) on the shape and morphology of burrow construction of O. rotundata on foreshore and backshore of sandy beaches in Qeshm Island. Such studies could be important from point ichnology. So that, exegesis of the burrowing behaviour can help to better commentaries of the depositional conditions which ichnology is greatly applied. Also, study of trace fossils are useful as indicators of past sea-level position (Curran and White 1991) that it would first be necessary to do detailed morphological analyses of burrow types on sandy beaches of the world. On the other hand, our results will provide relevant information on the key features of a sensitive species inhabitant of Iranian sandy beaches; this could help in the design of conservation strategies to protect and manage a key coastal habitat, where human pressure and climate change will jeopardize its ecological values over the next few decades.

Methods

Study Area

This study was conducted during the reproduction period of the ghost crab O.rotundata (April to September) and during the non-reproduction period (October to March) in Salakh sandy beach, Qeshm Island, Iran (Fig. 1) in February 2015, May, August and November 2016.

Fig. 1
figure 1

Location of study site at Qeshm Island for Ocypode rotundata

Based on topographic features, three areas were demarcated for Salakh beach which were foreshore, backshore, and dune zones. In this study, burrows were investigated on two zones: foreshore and backshore. The foreshore had a relatively flat surface with a seaward slope of 5°. The backshore was also flat, with a slope 2° landward. The dune zone was covered by vegetation, e.g., Halopayrum macronatum and Salsola vermiculata. Burrow distribution of O. rotundata was observed from the upper foreshore to the lower dune zone.

Experimental Design

The burrows of O. rotundata were distributed in 4000 m2 of sandy beach that along four transects (10 m width, positioned perpendicularly to sea, 10 m length), galleries with different diameters were chosen randomly; next, gypsum slurry (2:1 volume ratio of gypsum to water) was poured into the holes to obtain molds of the galleries` morphology. The distance of each burrow to the shoreline was noted. Crabs emerged from the burrows when the slurry poured and kept in plastic bags separately and labelled for further analysis. When the gypsum is set, the casts were carefully excavated by hand. The parameters of burrow morphology were measured, which included burrow opening diameter, length, depth, horizontal length and α and β angles (Fig. 2A, B). Out of the 106 burrows that were cast, 100 shapes were completed (56 burrows in non-reproduction period and 44 burrows in reproduction period), which were categorized into nine morphologies, as follows: single Shaft burrows (Fig. 3A), J-shaped (Fig. 3B), Y-shaped with short branch (Fig. 3C), Y-shaped (Fig. 3D, E), spiral burrow (Fig. 3F, G), complex burrows (Fig. 3H-J), J-shaped with extra branch (Fig. 3K), Y-shaped with an extra branch and anchor-shaped (Fig. 3L) and (Fig. 4). Regardless of this classification, the number of observed Y-shaped burrows with short branch, anchor-shaped, J-shaped with extra branch, Y-shaped with extra branch were very low, and thus were considered to be in the category of Y-shaped, complex burrows, that are J-shaped and Y-shaped, respectively. The carapace width and length of each crab were measured using vernier calipers, then weighed (analytical balance, ± 0.01 g, AND, EK-300i, Japan) and sexed according to abdomen, gonopod and pleopod morphology. According to size at the onset of maturity (Naderi et al. 2018b), all collected crabs were divided into two size classes of juvenile (mean: 25.52 ± 9.39 of carapace width) and adults crabs (mean: 44.55 ± 3.88 of carapace width).

Fig. 2
figure 2

A: a: length of burrow cast, b: depth of burrow cast; c: length of primary arm of burrow cast; d: length of secondary arm of burrow cast; e: horizontal length of burrow cast; B: α: ground angle of open arm, β: ground angle of blind-ending arm in Ocypode rotundata

Fig. 3
figure 3

Burrow casts of Ocypode rotundata on Salakh beach, A: single shaft; B: J-shaped; C: Y-shaped with short branch; D-E: y-shaped; F-G: spiral; H-I-J: complex; K: J-shaped with extra branch; L: anchor

Fig. 4
figure 4

Types of constructed burrow in reproduction and non-reproduction period by Ocypode rotundata

Sediment samples with three repetitions (500 g) were collected at 0–50 cm and 50–100 cm depth intervals which used for analyzing the vertical differences in sediment composition.

The organic matter content of the sediment was calculated by loss of weight on ignition at 500 °C. To do this, three 10 g subsamples were placed in porcelain capsules and incinerated at 500 °C for 3 h (Mantelatto and Fransozo 1997).

$$OM=\frac{W2-W3}{W2-W1}\times 100$$

where OM is organic matter content, W1 is Weight of tin (g), W2 is Weight of dry sediment + tin (g) and W3 is Weight of sediment after ignition + tin (g).

For calculating sediment moisture percentage, three 10 g subsamples were weighed and then dried at 105 °C in oven until a constant weight was attained. The following formula used to calculating sediment moisture percentage:

$$MC=\frac{w2-w3}{w3-w1}\times 100$$

where MC is sediment moisture percentage, W1 is Weight of tin (g), W2 is Weight of moist sediment + tin (g) and W3 is Weight of dried sediment + tin (g).

For the study of sand granulometry, sand samples (100 g) were taken from a depth of 0–50 cm and 50–100 cm in the last burrow observed both seaward and landward. Seven grades of sand were obtained by dry sieving, as follows: Coarse silt (31–62 µm), very fine sand (62–125 µm), fine sand (125–250 µm), medium sand (250–500 µm), coarse sand (500–1,000 µm), very coarse sand (1–2 mm) and granules (2–4 mm) (Buchanan 1984).

For determination of sediment pH, three 70 g subsamples were placed into the glass beaker, added 70 cc distilled water for 24 h, and then pH of sediment was obtained by pH meter (Metrohm, ‌827, Switzerland).

GRADISTAT softwar (Version 4.0) was used to determine grain size (Blott and Pye 2001).

Statistical Analyses

Pearson correlation coefficient was estimated to determine the relationship between burrow opening diameter and carapace width. t-test was also performed to compare between length of primary and secondary branch of Y-shaped burrows. The data were analyzed by using statistical software R version 3.5.0. The proportions of constructed burrows were compared by multinomial proportions analysis.

Results

Overall, 100 shapes were completed. On the other hand, 30 crabs emerged from the burrows were collected (17 females and 13 males). Relationship between carapace width and burrow diameter (Pearson’s correlation coefficient = 0.57, t = 3.48, df = 25, P = 0.001), burrow depth (Pearson’s correlation coefficient = 0.37, t = 2.13, df = 28, P = 0.04), were positively correlated (Fig. 5A, B). Relationship between carapace width and horizontal length of Y shaped burrow (Pearson’s correlation coefficient = 0.57, t = 1.7, df = 6, P = 0.13), α angle (Pearson’s correlation coefficient = 0.03, t = 0.03, df = 26, P = 0.96) were uncorrelated which were showed in Fig. 5C, D. Also, the scatterplot and the regression line of carapace width against ratio of burrow diameter/carapace was showed in Fig. 6. There was not any significant correlation between carapace width and ratio of burrow diameter/carapace (Pearson’s correlation coefficient = 0.27, t = -1.48, df = 26, P = 0.14). According to results, constructed single Shaft burrows by host crabs (28.76 ± 11.04 mm of carapace width) had a mean opening diameter of 28.08 mm and depth of 44.14 mm. Spiral burrows had a mean opening diameter of 31.05 mm and were produced by crabs with mean carapace width of 37.36 mm which had longer depth with mean of 67.75 mm than other three types of burrow. J-shaped burrow with mean opening diameter of 32.14 mm were produced by crabs with mean carapace width of 32.22 mm. Y-shaped burrow had the smallest opening diameter with a mean depth of 43.87 mm which had constructed by crabs with mean carapace width of 26.6 mm. The single Shaft entrance was more inclined than other burrows (Fig. 5E). Only collected male crab from complex burrow had 48.4 mm and 42.8 mm width and length of carapace respectively which made burrow with opening diameter of 60.8 mm and depth of 37 mm.

Fig. 5
figure 5

Regression relationship between carapace width and A burrow diameter; B burrow depth. Spot plot of carapace width and A horizontal length (just in Y shaped burrow); D α angle. E The variations of carapace width of crabs, opening diameter, depth and inclination angle among the J-shaped, Y-shaped, spiral and single tube burrows in Ocypode rotundata

Fig. 6
figure 6

Spot plot of carapace width and ratio of burrow diameter/carapace in Ocypode rotundata

Based on carapace width, percentage of J-shaped and Y-shaped burrow were more in female young crabs (Fig. 7A). Whilst, percentage of constructed single Shaft, J-shaped burrows were more frequent in male adult crabs. Also, created complex burrow was observed just in male adult crabs (Fig. 7B).

Fig. 7
figure 7

Percentage of constructed burrows by sex of A young crabs and B adult crabs in Ocypode rotundata

Y shaped burrow was had most frequency (22 burrows) in comparison to others. On the other hand, frequency of Y (19 burrows), J-shaped (16 burrows) and spiral (11 burrows) were more observed in non-reproduction period whilst most number of single Shaft burrows (14 burrows) were collected in reproduction period. Furthermore, percentage of complex burrow was significantly low (P < 0.05). Conversely, there was not observed significant difference between other burrows (Table 1). Also, results showed that percentage of Y shaped and complex burrow was significantly lower and higher in reproduction period respectively. Moreover, there was not any significant difference between other constructed burrows in reproduction and non-reproduction period (Table 2).

Table 1 Comparison of constructed burrow types percentage by Ocypode rotundata in reproduction and non-reproduction period by proportion test
Table 2 Comparison of each constructed burrows percentage Ocypode rotundata between reproduction and non-reproduction period by proportion test

Single Shaft, J-shaped, Y-shaped and spiral burrows were observed in months of February, May, August and November, whilst complex burrows occurred only in May and August (reproduction period) (Fig. 4).

Figure 8 shows Percentage of the different burrow morphologies.Based on sex, male crabs were not only more prompted to make single Shaft and spiral burrows, but also constructed more complex burrows. Y-shaped and J-shaped burrows were more frequent in female crabs.

Fig. 8
figure 8

Percentage of the different burrow morphologies of Ocypode rotundata

Figure 8 shows percentage of made burrows by sex during reproduction (5 males and 4 females) and non-reproduction period (8 males and 13 females).

Some differences were observed in frequency and morphology of burrows between foreshore and backshore. Most of the single Shaft burrows occurred in the foreshore zone during the reproduction period, with a mean opening diameter of 15.8 ± 6.0 mm, length of 265 ± 78.6 mm and α angle of 82.5°. J-shaped burrows had a mean opening diameter of 25 ± 14.5 mm, length of 537 ± 230.5 mm, depth of 404 ± 153 mm and α angle of 53.6°. For spiral burrows, opening diameter was 54.5 ± 7.91 mm, in length was 700 ± 149.9 mm, and depth reached to 660 ± 143.4 mm, with α angle being 57.3°. The last single Shaft burrow seaward distance was 480 cm away from shoreline in average 168.7 cm. This distance was 550 cm, 430 cm and 116 cm for last J-shaped, spiral and complex burrows seaward, respectively. The diameter, length and depth of single Shaft, J-shaped and Y-shaped burrows decreased with increasing distance from the shore during the reproductive period. Conversely, these parameters decreased increasing distance from the shore during non-reproductive period. This trend was similar to that observed in spiral burrows in the reproductive and non-reproductive period. Accordingly, diameter, length and depth of O. rotundata burrows increased landward. A similar result occurred for complex burrow (Table 3).

Table 3 Measurement data for Ocypode rotundata burrows within the foreshore and backshore zones

Most of the Y-shaped burrows were observed on the foreshore zone in the non-reproductive period, having a mean opening diameter of 15.8 ± 7.14 mm, length of 471 ± 277.8 mm, depth of 445 ± 241.3 mm, horizontal length of 176 ± 157.1 mm and α-angle of 61° and β-angle of 51.18°. The primary branch of all Y-shaped burrows was significantly longer than the secondary one, (t-statistics = 4.23, df = 16, P = 0.001), which joined together into a straight shaft with a mean length of 156 ± 53.3 mm. The secondary branch with a spherical end did not extend up to the surface. Primary branches faced seaward, whilst secondary faced landward. In most of the complex burrows, main entrance was divided into 2 shafts: one extended straight down and another turning right or left contrariwise to the first branch (Fig. 3J). These burrows had the largest mean diameter, with 65.8 ± 6.94 mm. Most of these burrows occurred on the foreshore zone during reproduction period, with a mean opening diameter 64.3 ± 7.06 mm. Ground angle for open branch was 50.66°.

The percentage of organic matter and moisture in relation to depth are shown in Tables 4 and 5, respectively. The results showed a decrease of moisture content from foreshore to backshore. The same result was observed for the percentage of organic matter in depths of 50–100 cm (Table 4). Maximum moisture content was obtained in depths of 50–100 cm for three studied zones (Table 5). Grain size showed variations in sediment composition, from fine to medium and coarse sand (Table 6). The mean value of pH was 8.1 where crabs were distributed.

Table 4 percentage of organic matter on Salakh beach in different depths (cm)
Table 5 moisture content on Salakh beach in different depths
Table 6 Granulometry data relating to the sand samples collected from Salakh beach

Discussion

The trace fossil O. rotundata commonly occurred from the backshore to the foreshore of Salakh sediment with single shaft, J, Y, spiral and complex shaped. Increasing size of O. rotundata accompanied by morphological changes (e.g., shape, diameter, length, and orientation) of burrows. In this study, the burrows with smaller diameter mostly occurred on upper foreshore (Table 7), which can be due to inability of their gills to tolerate long time of air exposure. This pattern of burrow distribution was similar to that reported for O. cursor (Shuchman and Warburg 1978), O. ceratophthalmus (Chakrabarti 1981; Chan et al. 2006). The mean depth of single Shaft and J-shaped burrows on foreshore zone were 309 and 400 mm respectively which were more than the depths measured in O. quadrata on the coast of Texas (Hill and Hunter 1973), O. ceratophthalma in Hong Kong ‌(Chan et al. 2006), India (Chakrabarti 1981) and Taiwan (Takahasi 1932). In our study, J-shaped burrows of juvenile crabs lacked any extra branch, as observed in India, Taiwan, America and Hong Kong (Takahasi 1932; Chakrabarti 1981; Chan et al. 2006). However, only two burrows had an extra branch, which is constructed by adult crabs. The orientation of the main branches was seaward, whilst the secondary branch was landward. This pattern was similar to crab populations in Hong Kong and India (Chakrabarti 1981; Chan et al. 2006). Conversely, none of the secondary branches of Y-shaped burrows opened to surface on Salakh beach, whilst a limited number of the secondary branches opened to surface on the Texas, India and Hong Kong coasts (Hill and Hunter 1973; Chakrabarti 1981; Chan et al. 2006). We hypothesize that O. rotundata constructed secondary branches in order to further safety and prevent predation which can be reason for low abundance of Y-shaped burrows in reproduction period. At the beginning of the reproduction period, from March to October (Naderi et al. 2018b), male crabs excavated complex burrows. Also, males created sand pyramids facing seaward in front of the burrows during daylight in high tide. Pyramids in front of the burrows of O. saratan and O. ceratophthalma were similar to those observed by Takahasi (1932) and Linsenmair (1967). Farrow (1971) stated that female crabs have Y-shaped burrows whilst males have spiral burrows; however, crabs of both sexes in O. rotundata were found to construct both Y-shaped and spiral burrows, as found in the study of Chan et al. (2006) on O. ceratophthalma. In our study, the trend of female crabs constructing mainly Y-shaped burrows whilst males made spiral burrows was different to the pattern noted by Chakrabarti (1981) in O. ceratophthalma.

Table 7 Measurement diameter burrow of Ocypode rotundata within the foreshore and backshore zones

With concerning to recruitment of O. rotundata from March to October (Naderi et al. 2018b), the main reason for high created single shaft burrow abundance is strong presence of megalopae (stage marks the passage from the sea to the beach (Smith 1873)) on foreshore zone where substrate moisture content is suitable for juvenile crabs.

During burrowing activities of O. rotundata, 4 types of sand disposal behaviour were observed: 1) Knocking sand with walking legs, 2) Throwing sand, 3) Stacking sand, which is similar to the observations of Schober and Christy (1993), Crane (1941), Vannini (1980) and Linsenmair (1967) and 4) Slamming sand with both chelae (especially major chelae) (Fig. 9) which done by male crabs. Like sand disposal behaviour of slamming, stacking sand was only made by male crabs, which is a sign for attraction of female crabs during the mating period (Schober and Christy 1993).

Fig. 9
figure 9

Types of sand disposal behavior in adult crabs of Ocypode rotundata: A Knocking, B throwing, C Slamming, D stacking

Salakh beach included coarse, medium and fine sand. The coarser granulometry between 0 and 50 cm depth on the shoreline zone could represent the accumulation of chopped oyster shells, coarse particles, etc. The results showed three grain size classes: symmetrical (foreshore and shoreline between 0–50 cm depth), coarse skewed (backshore between 0–50 cm depth) and very coarse skewed (shoreline and backshore between 50–100 cm depth). According to results, mating and reproductive behaviour of O. rotundata mainly occurred on the upper foreshore, with increasing diameter and length of burrows during the reproductive season. The reason for crabs migrating to the foreshore during reproductive period is related with the higher percentage of substrate moisture in comparison with backshore zone, which provide suitable substrate for burrowing and surpassing most of the egg development phase (Fellows 1973; Trott 1998; Negreiros-Fransozo et al. 2002). This could be an explanation for the relatively low number of ovigerous females found on the beach surface in several studies (Brooke 1981; Negreiros-Fransozo et al. 2002; Correa et al. 2014; Naderi et al. 2018b). Thus, the sand moisture could be a limiting factor in the distribution and abundance of ghost crabs, and hence an explanation of the low abundance of O. rotundata on the backshore zone throughout the year. On the other hand, ghost crabs are able to take up water from damp sand through the trichomes of setae, which are located between the third and fourth pereiopods (Wolcott 1978; Eshky 1985). In addition, none of the burrows investigated penetrate into the water table. As expected, the water content of sand decreased with distance from the sea, which was associated with an increase in burrow depth. Similar results were found on O. cursor (Hayasaka 1935; Warburg and Shuchman 1979; Williams 1983; Antia 1989; Strachan et al. 1999). Furthermore, Shinoda et al. (2019) found that crabs construct their burrows in appropriately wet zones.

The reason for curved burrows of O. rotundata depends on the heavy use of the same side for digging activity (Ali-Adnan 1985; Atkinson 1974). Conversely, Naderi et al. (2018b) obtained an equal number of left and right directions in burrows of O. rotundata. Other studies show similar results for O. ceratophthalma (Barrass 1963), O. saratan (Linsenmair 1967; Eshky 1985) and O. cursor (Strachan et al. 1999). We observed that the major claw was always the last part to access the burrows, which probably means that this part is used as a protective shield during resting time. In addition, Vannini (1980) observed this pattern for O. ryderi.

Burrows constructed (e.g., I, J, Y and spiral in shape) by O. rotundata in this study had a similar architecture to the results obtained in several studies on O. ceratophthalmus (Hayasaka 1935; Takahasi 1935; De 2005; Chan et al. 2006; Seike and Nara 2008; Lim et al. 2011), O. sinensis (Seike and Nara 2008), and O. rotundata (Zahedi et al. 2014), although some variability in shapes was evident. Lim et al. (2011) hypothesized that burrow shape was affected by sediment grain size. Also, they noted that Y-shaped burrows tend to be constructed in fine sediment, whereas J-shaped burrows are common in coarse granulometry. In our study, the casts of our burrows presented all kinds of shapes, occurring on the foreshore and backshore zone of Salakh beach, having coarse, medium and fine sediment. A similar pattern was observed on some beaches in Qeshm island (Kani, khorbas cave and Souza beaches) (Naderi unpublished data).

According to our results, O. rotundata individuals excavate single Shaft, J-shaped and Y-shaped burrows at the different growth stages, which was similar to observations reported by Lim et al. (2011). On the other hand, spiral and complex burrows only occurred in the adult crabs.

O. rotundata excavates burrows with different angles throughout the year cycle. The entrance branch of all burrow shapes was steeper during the reproduction period (Table 8), which can be a strategy of crabs for more safety and protection against predators. Also, the depth of all constructed burrows was bigger in non-reproduction period. In our study, we hypothesize that the increase in all depths of created burrows, in November and February, could be related to the reduction in weather. Lim et al. (2011) noted that steeper-sloped burrows in the opened branch of the Y-shaped burrows may provide better protection from predators since they provide a rapid retreat to deeper depths. This description is inconsistent as it fails to explain why J-shaped burrows do not present the same pattern. Results from our study showed that temporal variation had a significant impact on burrow slope. Also, Takahasi (1932) and Fellows (1966) pointed out that the tunnel in burrows becomes more horizontal with the increase in weight of the crab and substratum instability. In this study, we noted no relationship between those parameters and the angle of burrow entrance (Fig. 10).

Table 8 angles of constructed burrows of Ocypode rotundata during reproductive and non-reproductive periods on Salakh beach
Fig.10
figure 10

Spot plot of mass and slope burrow in Ocypode rotundata

According to semi-terrestrial O. rotundata, they closed their entrance of burrow as dome shaped before tidal inundations to provide a chamber of air at the top of the burrow. It seems that this activity could be helped to protect sand moisture in the burrow in order to renew their respiratory water (Fig. 11).

Fig.11
figure 11

plugged burrow by Ocypode rotundata

Generally, our results showed that temporal variation had significant impact on burrow slope, so that the entrance branch of all burrow shapes was clearly steeper during the reproduction period in contrast to the same burrows created in the non-reproduction period. Also, O. rutondata selects upper foreshore zone for mating and reproduction behavior.