Community based actions save Yellow-spotted river turtle (Podocnemis unifilis) eggs and hatchlings flooded by rapid river level rises

The conservation and recovery of increasingly threatened tropical freshwater turtle populations depends on effective management plans and actions. Here we show that community-based actions saved Yellow-spotted river turtle (Podocnemis unifilis) eggs submerged by unseasonal flooding and ensured the release of hatchlings. We recovered 926 eggs and 65 premature hatchlings from 74 submerged nests at 16 flooded nesting areas along 75 km of waterways. The rescued eggs were transferred to a rearing center and incubated. Hatchlings emerged from eggs that had remained underwater for up to two days. Hatchlings were maintained in 250–500 L nursery tanks until yolk sac scars had closed. Healthy hatchlings were then immediately released around the original nesting areas. We released 599 hatchlings (60.4%) from 991 submerged eggs and hatchlings. Egg survival (61.7% (571/926)) was substantially less than hatchling survival (94.2% (599/636)) but within the expected range of values reported for this species. These findings suggest that Yellow-spotted river turtle eggs and embryos are resistant to short-term submersion, which could help explain the widespread distribution of this species across highly seasonal Amazonian rivers. Management plans should take the possible survival of submerged eggs into consideration as part of species conservation and recovery actions.


INTRODUCTION
Like many tropical species Amazonian freshwater turtles are threatened by deforestation (Fagundes et al., 2018), climate change (Eisemberg et al., 2016), and unsustainable exploitation (Rachmansah, Norris & Gibbs, 2020;Smith, 1979). Precautionary estimates suggest that populations of the once abundant Yellow-spotted river turtle (Podocnemis our results to discuss how management plans could include rescue of submerged eggs and premature hatchlings as part of species conservation and recovery actions.

Ethics statement
All methods were carried out in accordance with relevant guidelines and regulations. The actions we report were developed and applied following international (Beaupre et al., 2004) and national (Balestra, 2016) guidelines. Fieldwork and data collection were conducted under research permit numbers SISBIO 63668-1 and 63668-2 to DN, issued by the Brazilian Ministério do Meio Ambiente ("MMA"). Fieldwork with turtles was approved by the Animal Use Ethics Committee from the Federal University of Amapá (UNIFAP) (CEUA-UNIFAP approval 009/2017).

Study area
The study was conducted in the Araguari river basin, located in the Brazilian State of Amapá (Fig. 1). Nesting of Podocnemis unifilis has been monitored since 2011 in the study area and here we present a brief summary of the previously described study area and P. unifilis nesting areas (Michalski et al., 2020;Norris & Michalski, 2013;Norris, Michalski & Gibbs, 2018a, 2018bQuintana et al., 2019). Climate in the area is characterized as equatorial monsoon (Kottek et al., 2006) with an annual rainfall greater than 2,000 mm. Rainfall is strongly seasonal, with the dry season extending between September and November (total monthly rainfall <150 mm), and the wet season from February to April (total monthly rainfall >300 mm (De Oliveira, Norris & Michalski, 2015;Paredes et al., 2017)).

Data collection
As part of a collaborative monitoring scheme with the local community, we monitored 29 P. unifilis nesting areas from September to December 2018. This is a representative subset of the available nesting areas described previously (Norris, Michalski & Gibbs, 2018a;Quintana et al., 2019) that were chosen in 2018 to encompass broad gradients (environmental and anthropogenic) and include areas with the majority of nests (including nesting areas with >90% nests recorded in previous years). Nesting areas were visited at least once a month to search for turtle nests and to monitor conditions of nests that were individually marked with a small stick. Nests were located by following turtle tracks on the sandy/gravel substrates and systematic substrate searches with blunt wooden sticks and fingertips (Escalona & Fa, 1998;Norris, Michalski & Gibbs, 2018a;Quintana et al., 2019). Searches were conducted by teams of two to three observers along lengthwise transects throughout each area at a standardized speed (mean 0.8, range 0.2-1.3 km per hour); time spent searching nesting areas ranged from 10 to 97 min depending on the size of the area (Norris, Michalski & Gibbs, 2018a).

Recovery of submerged eggs
In 2018, the work with local communities could be classified as collaborative monitoring with external data interpretation (Danielsen et al., 2009). In addition to monitoring nesting areas, members of the local communities also participated in management actions including installation and monitoring of nest predator exclusion devices (Norris, Michalski & Gibbs, 2018b). Previously, strategically targeted community-based actions at only two nesting areas enabled the successful protection of 75% of nests along a 33 km stretch of river (Norris, Michalski & Gibbs, 2018b). In 2018 actions were expanded to include two monitoring teams working simultaneously (one team on each of the Araguari and Falsino rivers).
In 2018 there was an unprecedented and rapid rise in river levels ( Fig. 2). On the 29 November 2018 the field team from the Falsino river reported rising river levels due to unseasonal early rains and that on 28 November some nests had been moved to higher points on at least one nesting area. Based on our knowledge from previous years in the region there was no expectation that the river water level would continue to rise. However, after assessing the river level at the field site on 30 November it became clear that the rapid water level rise was flooding nests along both Araguari and Falsino rivers. Members of the local community living in the region for more than 30 years remember early rain, but do not recall river levels rising so quickly. Hence, as this was an unprecedented event, our team took action to try and save at least some of the nests.
To reduce the loss of eggs and hatchlings our plan of action followed these steps (Fig. 3): 1. Move nests to higher locations that would not be flooded keeping them in the same nesting area where eggs were laid.
3. Transport and dry sand from the nesting beaches (needed for incubating eggs and strengthening of hatchlings).

Set up incubators for eggs.
5. Transfer hatchlings to a "nursery" for strengthening.
6. Release hatchlings in the same region where they were rescued.
By 1 December 2018 it was clear that the river levels would not decline and that we would need to take additional steps (i.e., moving nests to higher points within the same nest site was no longer a viable option). We established an improvised rearing facility located on the property of a local landowner. We created incubators for eggs and temporary nursery tanks for premature hatchlings. Eggs were placed in incubators (50-80 L polystyrene foam boxes) filled with sand from the submerged beaches and heated by lamps with the aid of a generator that was turned on at least four hours a day (Fig. 3). Hatchlings were transferred to water tanks (250 or 500 liter) that served as nursery areas for strengthening and to enable complete absorption of the yolk sac ( Fig. 3). In the nursery tanks fresh vegetation (cassava (Manihot esculenta) leaves) was also provided for both food and shelter. Hatchlings were periodically (every 2-3 days) weighed and measured (straight-line carapace length and width) to monitor their health and development.
Hatchlings were released when umbilical scars were closed. Prior to release hatchlings were examined for any external sign of disease (i.e., fungus, botfly larvae) and weighed. Although all hatchlings were examined, due to logistical difficulties it was not always possible to weigh hatchlings prior to release. We were able to weigh 409 (68%) of the 599 released hatchlings (145 and 264 from Araguari and Falsino rivers, respectively). Healthy hatchlings were then released close to the original nesting area. By the time of hatchling release (up to 1 month after rescue) all nesting areas had been flooded and were totally submerged. We therefore selected release areas that were upstream and close (within 200 m) to the original nesting areas with relatively calm water and access to shelter and food along the river banks (i.e., small perennial streams that join the main rivers).

Data analysis
All statistical analysis and graphics production were undertaken within the R language and environment for statistical computing (R Development Core Team, 2020). One-way Analysis of Variance (ANOVA) was used to test for differences in body mass and condition between hatchlings from Araguari and Falsino rivers. We calculated two body condition indexes for released hatchlings to evaluate the relative condition on the day they were released. The indexes were obtained from the residuals of Generalized Additive Models (GAMs, Wood, 2011Wood, , 2017 of the response of ln transformed body mass predicted separately by ln transformed straight-line carapace length and ln transformed carapace area. The carapace area (cm 2 ) was calculated as an ellipse by standard formula from the radius of the major (ra) and minor (rb) axis (expressed in cm) as follows: We used GAMs to avoid known problems of residual body condition indexes obtained from linear (ordinary least square) regression, namely the assumption that mass and body size relationships are linear (Green, 2001). GAMs include non-parametric smooth terms and are able to model non-linear relationships (Wood, 2011) such as those

RESULTS
The river water level rose 2.2 m (from 343 cm to 561 cm) over 16 days between 24 November and 10 December 2018 (Fig. 2). On average P. unifilis nests were laid at a height equivalent to the river level of 472 cm (Fig. 2). The highest nests were laid at 526 cm, which means that 100% of nests would have been entirely submerged by 8 December (Fig. 2). Eggs and premature hatchlings were recovered from a total of 74 submerged nests in 16 nesting areas along 74.6 km of rivers (Table 1). On average nests were laid at similar heights above the water along both rivers (Table 1), but submerged nests were detected and recovered earlier along the Araguari river (between 30 November and 2 December) compared with the Falsino river (between 1 and 5 December), which had the highest nesting areas with nests laid 69 cm above the maximum height on the Araguari river (Table 1). Nests were submerged for 0 to two days (Table 1). Premature hatchlings were only recovered alive from waterlogged nests that had water infiltrating from below (i.e., 0 days submerged). Eggs were recovered from nests that were waterlogged (0 days submerged) to two days submerged. At the time of recovery the depth of water ranged between 0.05 to 1.30 m above the nests. The overall survival rate of eggs and pre-mature hatchlings was 60.4%, with 599 hatchlings released from 991 eggs and pre-mature hatchlings ( Table 2). A total of 582 hatchlings emerged from 926 eggs that had been submerged for up to 2 days. Hatchling emerged more than a month after submersion, with the last hatchlings emerging on 10 January 2019 (36 to 40 days after recovery). The vast majority (94.2%) of emerged and pre-mature hatchlings survived (Table 2) and all surviving hatchlings were released by 31 January 2019.
The body mass of released hatchlings varied two fold (from 9.8 to 23.2 g, Table 2). Although mean body mass was significantly greater for hatchlings released on the Falsino river (Table 2, one-way ANOVA, F 1,407 = 8.00, P = 0.005), the increase of 0.6 g compared with body mass of Araguari hatchlings represented only 3% of the mean body mass (18.0 g) from all released hatchlings (Table 2). The relationships between ln transformed body mass and ln transformed body size (carapace length and area) were not linear (Fig. 4, GAM effective degrees of freedom (EDF) = 5.5 and 7.1 for straight-line carapace length  and carapace area respectively). The body condition index values ranged from −0.15 to 0.08 and −0.36 to 0.26 from the straight-line carapace length and carapace area models respectively (Fig. 4). Hatchling body condition increased slightly in Falsino hatchlings (Fig. 4). Body condition based on ln straight-line carapace length differed weakly between rivers (one-way ANOVA, F 1,407 = 5.51, P = 0.019) and a similar pattern was found with the residual body condition index derived from ln carapace area (one-way ANOVA, F 1,407 = 11.47, P = 0.0008). Although statistically significant, there were only small differences in mean condition index values between rivers (Araguari (−0.008, −0.011) and Falsino (0.004, 0.006), mean values from straight-line carapace length and carapace area models respectively) and the distributions of body condition values were also similar (Fig. 4). Additionally, the interquartile range of body condition index values for hatchlings from both rivers overlapped zero (Fig. 4).

DISCUSSION
To our knowledge we demonstrate for the first time that P. unifilis embryos are resistant to short-term submersion. We also show that healthy hatchlings survived and were successfully released following egg submersion. Early stages of freshwater turtles are known to suffer elevated mortality (Iverson, 1991). The mortality of P. unifilis eggs and hatchlings is usually elevated and varies widely, typically with 60-100% mortality in natural settings (Caputo, Canestrelli & Boitani, 2005;Escalona & Fa, 1998;Ferreira Júnior & Castro, 2010). Additionally, previous studies show that human harvest is also a major cause of nest loss in our area (Norris, Michalski & Gibbs, 2018b;Quintana et al., 2019). As such the release of hatchlings from over half of the submerged eggs (that would otherwise have suffered 100% loss) is within the expected survival of natural nests. The average body mass (18.0 g) of released hatchlings was similar to values reported elsewhere for this species. Previous studies report mean hatchling weights of 17.8 to 20.0 g in Colombia (Páez & Bock, 2004), 17.3 g in Ecuador (Caputo, Canestrelli & Boitani, 2005) and 20.8 g from hatchlings (n = 99) along the Falsino River in 2011 (Arraes, 2012). Although we obtained statistical differences in body mass between hatchlings from Araguari and Falsino rivers, such a small difference (0.6 g) is unlikely to be biologically relevant. We were unable to monitor the long-term survival of the released hatchlings; however, there is no reason to expect that they would experience any negative impacts of short-term rearing. We believe that we effectively reduced known issues of captive rearing (Burke, 2015) by maintaining hatchlings at low densities and for a short period of time.
Different eggshell types may provide different levels of protection to embryos (Deeming, 2018). The rigid cased P. unifilis eggshells may help to explain the resistance of embryos to short-term submersion (Packard, Packard & Boardman, 1982;Winkler, 2006). Although P. unifilis eggs are rigid, the shells remain somewhat permeable, for example allowing passage of agro-toxins (Hirano et al., 2019). Yet P. unifilis eggshells have a low pore density (Winkler, 2006;Winkler & Sánchez-Villagra, 2006), which is likely to limit exchange of water with the environment (Packard, Packard & Boardman, 1982).
Although rigid eggshells can resist hydric changes, embryos are susceptible to changes in temperature associated with submersion. As P. unifilis only occurs within the tropics Rachmansah, Norris & Gibbs, 2020), temperature is likely to be a key determinant for the survival of all stages and it is unlikely that embryos have any adaptation to resist low temperatures such as those experienced under flowing rivers. The temperature range for successful incubation of P. unifilis nests is typically between 27-32 C (De Souza & Vogt, 1994;Ferreira Júnior & Castro, 2006;Páez & Bock, 1998;Páez & Bock, 2004). Although temperature dependent sex determination has been well studied in P. unifilis (De Souza & Vogt, 1994), other aspects such as embryo survival in relation to temperature remain poorly investigated. For example, in the pivotal study for the species, De Souza & Vogt (1994) state "there was high mortality among the eggs switched from 34 C to 28 C", providing no information on when the change was made or quantity (n or %) that died.
Intervening to prevent losses is often a necessary component of wildlife conservation projects. Yet, there are ethical and moral concerns associated with these actions. Turtles can be considered a special case in several ways. All turtles are oviparous and P. unifilis do not provide any post-oviposition parental care, therefore, no surrogates for parental care and training are needed to rear eggs and hatchlings (Burke, 2015). Indeed rearing P. unifilis hatchlings has become an integral part of management approaches across the species range (Harju, Sirén & Salo, 2018;Páez et al., 2015;Sinovas et al., 2017). We took reasonable precautions to address ethical concerns, including providing appropriate temperature, water and food (reflected in the increase in body weight), avoiding risks of disease spread from exotic species, and ensuring that premature hatchlings were retained for the minimal time necessary before release to the wild. We suggest that the improvised response to unprecedented river level rise was appropriate and provided a positive contribution to the conservation of P. unifilis populations along the rivers. Continued long-term monitoring is required to establish if such interventions enable freshwater turtle populations to persist and even recover across Amazonian rivers increasingly impacted by anthropogenic harvests  and flow-rate alterations (Castello et al., 2013).
Considering the drastic changes to Amazonian river flow regimes it seems reasonable for management actions to include contingency plans for recovery and release of submerged P. unifilis eggs and hatchlings. A total of US$5007 was spent (exchange rate of 3.78 Brazilian R$ to 1 US$, food, gasoline, daily stipends, rearing equipment) during the two month improvised rearing process (28 November 2018-31 January 2019, Fig. 3), representing an overall financial investment of US$8.36 for each P. unifilis hatchling released. In comparison, the community based monitoring implemented with two teams required on average US$2302 per month (food, gasoline, daily stipends) that is, US$4604 over a two month period. Therefore improvised rearing does not necessarily represent an excessive financial investment beyond that expected as part of the actions for conservation and monitoring of the species.

CONCLUSION
The resistance of P. unifilis eggs and embryos to short-term submersion could help explain their widespread distribution across tropical South America. The continued decline of the species across its range shows that such adaptations are insufficient to respond to unprecedented 21st century changes across riverine habitats. Our successful rescue of eggs and premature hatchlings from flooded nests and subsequent release of healthy hatchlings demonstrates the viability of recovery actions. Such actions are likely to be increasingly necessary across rivers where not only are flow rates affected by damming but also anthropogenic harvests can drive unsustainable harvests of turtle nests. We suggest that such recovery actions should be included in action plans for species conservation and recovery. Additionally, hydroelectric developments should evaluate localized changes in flow rates to avoid seasonal flooding of nesting areas and pro-actively avoid the need for recovery actions. James P. Gibbs conceived and designed the experiments, authored or reviewed drafts of the paper, and approved the final draft.

Animal Ethics
The following information was supplied relating to ethical approvals (i.e., approving body and any reference numbers): The Animal Use Ethics Committee from the Federal University of Amapá (UNIFAP) approved the study (CEUA-UNIFAP approval 009/2017).

Data Availability
The following information was supplied regarding data availability: The raw measurements used in analysis and figures are available in the Supplemental Files.