Substrate influences human removal of freshwater turtle nests in the eastern Brazilian Amazon

Substrate type determines nesting success and fitness in turtles and is a critical consideration for nesting area protection and restoration. Here, we evaluated the effect of substrate on nest removal by humans in the eastern Brazilian Amazon. We analyzed substrate composition and fate of 216 P. unifilis nests along 88 km of rivers. River segment and substrate type were the most important predictors of nest removal by humans. We found up to 36% lower removal of nests in fine sand and experimental results support the hypothesis that wind more often obscures tracks of nesting females in substrates with more (>66%) fine sand. Our findings are useful for informing the restoration of artificial nesting areas across the Amazon, as nesting area substrates should be selected not only to maintain hatchling sex ratios, but also to minimize nest removal by humans.


Results
A total of 216 nests (95 in the Falsino and 121 in the Araguari river) were found in 43 nesting areas (Table 1), yielding an average of 5.02 nests per nesting area (average ± SD = 5.02 ± 3.90, range = 1-17 nests) in both river segments. More than half (n = 127, 58.80%) of all nests were removed by humans, with the highest proportion of nest removal (75.21%) in the Araguari segment ( Table 1).
The grain size distribution of nest substrates varied greatly, from nests laid in 76% fine sand to 93% gravel ( Table 2, Figs. 1 and S1). However, the majority of nests were laid in substrates with higher proportions (> 50%) of sand ( Table 2, Table S1). Although there was variation in substrate grain sizes between nests, there was no evidence to suggest that grain sizes differed between the Araguari and Falsino river segments (Fig. S1).
Of the eight environmental, anthropogenic and spatial variables considered, nest removal by humans was affected (P < 0.05) by river segment and substrate type (Table 3). Turtle nests in the Falsino segment were removed less by humans when compared to nests located in the Araguari segment (Table 3). Marginal predictions from the GAM (where effects of other variables were controlled using mean values of continuous variables and factor reference levels) showed an average threefold difference (28.4% to 88.9%) in removal across the different substrates (range of PC1 values). Indeed, nests were removed less by humans when located on substrates with www.nature.com/scientificreports www.nature.com/scientificreports/ higher proportions of fine sand, the substrate type generally represented by our first PCA axis (Table 3, Fig. 2). This pattern of reduced removal with increasing proportions of fine sand was consistent across both river sections, with predicted removal declining in parallel (Fig. 2).
In laboratory experiments, the conspicuousness (and therefore detectability) of tracks declined abruptly at a threshold of 66-76% fine sand (Fig. 3). The proportion of visible tracks declined with increasing wind speed such that at speeds of 5.6 m/s (equivalent to "moderate breeze" on the Beaufort scale) no tracks were visible in substrates with 66-76% fine sand (Fig. 3). Yet, the majority of tracks were still visible in substrates with lower proportions of fine sand, with 100% and 75% of tracks visible in substrates with 0-10% and 36-46% of fine sand, respectively (Fig. 3).

Discussion
In this study across a large river extension of eastern Brazilian amazon we (1) quantified the grain size distribution of substrate types used by Podocnemis unifilis to lay nests and (2) showed that substrate grain size affects human removal of turtle nests. We first turn to explore the flexibility showed by P. unifilis in nesting areas in our study region and then demonstrate that nest-site substrate type affects nest removal rates by humans. Finally, we provide evidence to guide nest-site restoration by identifying which substrate type reduces nest removal by humans and will likely maximize turtle reproductive success along rivers accessible to humans. Overall, Podocnemis unifilis laid nests in substrates with a wide variety of grain sizes, ranging from very fine sand to gravel. This diversity of grain sizes was expected, as compared with other freshwater turtles P. unifilis is considered a generalist species in terms of nest-site selection 17,34,40 . Thus, our results demonstrating that nests can be found in substrates ranging from over 76% of fine sand to over 93% of gravel across 88 km of waterways corroborates the generalist behavior for nest-site selection of this species.
We found that river segment (Araguari versus Falsino) and substrate grain size (represented by the first axis of the Principal Component Analysis) were the strongest predictors of nest removal by humans across all turtle nest-sites. Nests laid along the Araguari river were more likely to be removed by humans than nests along the Falsino river. This result is not surprising as the Araguari river has a larger number of houses and increased boat traffic compared with the Falsino 30,41,42 . Our findings highlight that even around protected areas with relatively low human disturbance, it is possible to have markedly different anthropogenic impacts and these differences should be considered in management plans and conservation initiatives in and around sustainable-use reserves.
Our study is the first to quantify how substrate type affects nest removal rates by humans. A study of marine turtle (Caretta caretta and Chelonia mydas) nesting in Syria identified that the visibility of tracks from these large (~160 kg) turtles depended on the texture of the beach substrate, with tracks hardly visible in dry coarse gravel 39 . The same study found that the longevity of a turtle track also depended on the substrate, with even a light breeze   www.nature.com/scientificreports www.nature.com/scientificreports/ rapidly obliterating tracks on loose sand 39 . Our finding that turtle nests were less likely to be removed by humans in substrates with higher proportions of fine sand, also suggests that different substrates may produce differences in human detections of turtle tracks and nest marks. Our experimental manipulation of track longevity in different substrates exposed to different wind speeds corroborates this hypothesized link between human removal of turtle nests and track detectability. The reduced levels of nest removal that we found in nesting sites with more (>66%) fine sand could therefore be linked with the rapid obliteration of turtle tracks and nesting marks in the finer sand by even moderate breezes (i.e. of sufficient speed to raise dust or loose paper). Independent of the reason for reduced removal of fine-sand nests our results present compelling evidence that substrate grain size does affect turtle nest removal by humans and hence, must be considered an important variable in the management of turtle nesting areas with human activity.
Finally, due to the increasing number of new dams in the Amazon 19,22 and the likely submersion of many freshwater turtle nesting areas 25 baseline information for the construction of effective artificial nesting areas is required. Our results provide compelling evidence that substrate grain size does affect human removal of turtle nests and should be considered in the design of restored nesting areas. We suggest that using finer sand compared with coarse sand or gravel in restored nesting areas can increase nesting success along 164,971 km of rivers accessible to humans throughout the species range 28 . For example, as wind can obliterate tracks in loose sand 39 , combining the use of a thin layer of fine sand covering turtle nesting areas, could potentially reduce human removal. Also, as substrate texture of freshwater turtle nest-sites was related with nest temperature and sex ratio determination 38 , maintaining heterogeneous substrate types (i.e., mix of fine and coarse sand) in nesting areas is important. Substrate type has also been shown to be an important factor for hatchling success in some turtle species 35,43 . We were unable to assess hatchling success but the results from previous studies enable us to generate informed expectations for likely outcomes. Differences in grain size did not influence the hatchling success of P. unifilis in previous studies 35 , but showed to be important for the larger congener P. expansa, with nests placed on fines sediments having improved hatching success 35 . Additionally, in loggerhead turtles (Caretta caretta), hatchlings from fine-grain sand nests had higher crawling and swimming performance, and fine-grain sand was recommended for use in hatchery incubation beds to increase hatchling survival and performance 43 . It therefore seems likely that the use of fine-grain sand substrates would not negatively affect P. unifilis hatchling success. It will be important for future studies to carefully monitor egg and hatchling survival as part of any conservation and restoration actions involving nesting area substrates.
Our results demonstrate that although P. unifilis uses a variety of substrate types, conservation initiatives must take nest substrate into consideration not only for gender determination, but also to reduce human removal of nests. Proportions of 66% or more of fine sand in nesting areas can directly reduce human detection of turtle www.nature.com/scientificreports www.nature.com/scientificreports/ nests. Therefore, our results can help to inform the successful restoration of P. unifilis nesting areas along the 164,971 km of rivers accessible to humans across the species range 28 .

Methods
Study area. The study was conducted in the Araguari river basin, located in the Brazilian State of Amapá ( Fig. 4). Climate in the area is characterized as equatorial monsoon 44 with an annual rainfall greater than 2000 mm. The dry season extends between September and November (total monthly rainfall <150 mm), and the wet season from February to April (total monthly rainfall >300 mm) 45 .
This study was conducted along a total of 88 km of contiguous segments of Falsino and Araguari rivers, 42 km and 46 km, respectively (Fig. 4). These river segments were located between two sustainable-use protected areas, the National Forest of Amapá (hereafter FLONA), and the Amapá State Forest (hereafter FLOTA). Both segments can be classified as "clear-water", with a low density of suspended particles 46 , but they experience markedly different human activity 41 . The Araguari river is characterized by higher human density (we counted 13 houses in 2017) and more intense boat and fishing activities (including commercial fishing) 41,42 . In contrast there is no commercial fishing and only 4 houses along the Falsino river.
Turtle nest surveys and sampling. Monthly boat surveys were used to locate nesting areas. Nest searching was conducted between October and December 2017 during the nesting season of Podocnemis unifilis in the region 23,25 . While navigating along the rivers in a motorized boat at a constant speed (ca. 10 km/h) we performed an extensive search for nesting areas that involved identifying all suitable nesting areas through visually searching river banks, circling islands, stopping to search among boulders and rapids. Nesting areas were identified as areas>5 m 2 of exposed substrate sufficiently raised above the river level not to be waterlogged at a depth of 15 cm, a representative depth that females dig when nesting 25,34,47 . Although we cannot guarantee that all nesting areas were located across 88 km of rivers, we assume that the comprehensive survey effort enables us to provide a robust and representative evaluation of the nest-site substrates.
In all nesting areas, a search for nests was done together with local residents with over 30 years of knowledge on nesting sites 25 . Searches were conducted by a team of three observers 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. To minimize possible observer-related biases, at least one surveyor was constantly present with every search team throughout the entire study period. Logistical limitations hampered our capacity to visit all nesting areas more than once along 88 km of rivers. In order to standardize the number of visits for each nesting area, we included only data from a single visit during the first fortnight of November because this corresponded to the period when both turtle nesting and removal of nests by humans reached their peak during the nesting season 30 .

Substrate collection and grain size determination.
At each turtle nest we collected 200 g of substrate from the nest surface (1 -3 cm), which was stored in plastic bags labeled with the identifications of the nesting area and the nest. The location of each turtle nest was obtained in the field with a handheld Global Positioning System (GPS) accurate to +2 m. We could not evaluate hatchling success from turtle nests surveyed due to high rates of nest harvest by humans in the study area. For example, in 2017, 42% of the nesting areas along 33 km of the Falsino River (the river segment with the lowest anthropogenic disturbance) were harvested by humans 23 . Thus, hatchling success in our study region could not be linked directly with nest substrate.
All substrate samples were analyzed at the Laboratório de Análises de Sedimentos do Instituto de Pesquisa Cientificas e Tecnológicas Estado do Amapá -IEPA (Macapá, Brazil). To obtain grain size frequency distribution of substrate samples we used dry sieving, which is a well-established technique 48 . In the laboratory, all organic material (i.e., vegetation) was removed manually from samples with tweezers. The samples were dried at 50 °C for 24 hours, and were shaken by a mechanical shaker for 10 minutes through a stack of 12 sieves 49  Response variable. To quantify nest removal by humans we identified if each nest was removed or not by humans 30 , retaining individual nests as the sample unit in the analysis. Human removal of eggs was identified when an open nest was found with a typical depth between 10-15 cm, but lacking eggs. Signs of human activities in the nesting areas, such as footprints, charred wood from fires, campsite or trash were used for corroboration. Nests depredated by wild animals (identified by the presence of broken eggshells and/or partially eaten eggs around the nest, with animal excavation marks were present) though damaged were considered not removed by humans 30 . Explanatory variables. We examined a total of 8 (4 environmental, 2 anthropogenic and 2 spatial) variables that were likely to influence nest removal by humans (Table 4, Tables S1-S2). We measured four variables at each nesting area: (1) nesting area type (environmental: either island or bank); (2) river segment (anthropogenic: either Falsino or Araguari), (3) nest density (environmental: expressed as nests per m 2 ); (4) distance to the nearest house (anthropogenic: km). We also recorded two variables at each individual nest: (1) distance to water (environmental: m), and (2) substrate grain sizes (environmental: PC1). We also included geographic coordinates of the nests and the identity of the nesting area to control for possible spatial autocorrelation in the models.
Variables were selected based on previous studies conducted with P. unifilis, as well as our previous knowledge about the study area and species 23,25,30,42 . For example, nesting area type (island or bank) was found to be an important environmental variable for nest removal by humans at the scale of nesting areas in our study region 30 . River segment represents different intensities of human activity (i.e., high for Araguari and low for Falsino) 41 , and it was also highlighted as important predictor for nest removal by humans in previous studies 30 . Nest density was included as based on classical optimal foraging theory 50,51 we would expect relatively increased effort and increased removal where nest density is greater and nests therefore more frequently/easily found. Distance to houses was an important variable to explain nest removal in our study area 30 , being more important than number of houses within radius of 1 or 5 km for explaining removal between nesting areas. Finally, nest distance to water and proportion of substrate grain size could both affect the detectability by humans. As all nesting sites were located along navigable rivers, and the main transport for riverine community is by boats and canoes, distance of nests from water edge as well as different substrate grain size may result in different detectability and consequent different rates of human removal.
Density of nests was obtained by dividing the total number of nests by the surface nesting site area. This was measured by mapping the surface nesting area using a GPS handheld in situ. Rocks and dense vegetation (non-nesting zones) were excluded from density calculations. Distance to the nearest house was obtained from the centroid point of each nest site to the nearest house (based on GPS fixes obtained in situ at the sampling period). Nest distance to the water's edge was measured with a measuring tape at the day of the survey. Nesting area size was not included in our analysis as human removal of nests was not influenced by this variable in our study area 30 . Substrate type and track visibility. As humans often search for turtle nests by following the tracks of nesting females, previous studies suggest that differences in human removal of turtle nests can be explained by differences in detectability of the tracks in different substrates 39 . We experimentally tested different substrates to evaluate how the visibility of tracks was affected by different wind speeds. Following the dry sieving, we established three classes of substrate to represent the gradient of fine sand in the 216 nests: 0-10%, 36-46%, and 66-76%. We then selected five nests at random in each class.
In the laboratory samples were placed in tray providing a substrate depth of 6 mm. We drew four tracks in the tray, with size and depth following those recorded from tracks of adult females (Fig. S2). Based on the quantity of substrate collected and size of tray required to maintain a depth of 6 mm it was only possible to draw four tacks. Once drawn, tracks were exposed to different wind speeds (1.5, 3.5 and 5.6 m/s, measured using a digital anemometer) generated by a multispeed electric fan fixed at a constant height and distance from the tray. These speeds represent "light air", "gentle breeze", and "moderate breeze" on the Beaufort scale and are typical of those Data analysis. All statistical analysis and graphics production were undertaken within the R language and environment for statistical computing 52 . The analysis of grain size samples followed the method of Folk and Ward 48 with the packages "rysgran" 53 and "soiltexture" 54 . Principal Component Analysis (PCA) of our seven substrate grain sizes (Table S3) was used to summarize the information of substrate grain size sampled on P. unifilis nests to avoid multicolinearity in subsequent analysis. The PCA was obtained from a scaled and centered data matrix of our seven substrate grain sizes, and axes derived from the squared correlations coefficients with the variables. Based on the Kaiser-Guttman criterion only the first three principal components (with eigenvalues >1) should be interpreted 55 (Table S3). Component loadings, summary of principal components after varimax rotation (Table S4), and inspection of correlations revealed that the three axes with eigenvalues >1 (PC1, PC2, and PC3, representing 90% of the variance) were generally represented by proportion of fine sand (negative correlation), medium sand (negative correlation), and gravel (positive correlation), respectively (Fig. S1).
We used Generalized Additive Models (GAMs) to assess the effects of environmental and anthropogenic variables on turtle nests removed by humans. This approach allows the shape of the relationship between the response and the explanatory variables to be determined from the data, rather than following a prescribed functional form 56 . Thus, we analyzed the response of the nests removed by humans (binary variable), and examined the influence of the explanatory variables (nesting area type, river segment, nest density, distance to the nearest house, distance to water, and scores of PCA first axis related to substrate grain sizes) during the peaks of turtle nesting season and nest removal. A two-step process was adopted with GAMs run (with default settings) using the R package "mgcv" 57 . Firstly we established whether variables should be include as non-parametric (non-linear) smoothed terms by checking estimated degrees of freedom (EDF) values (typically an EDF value close to or less than 1 suggests linear relationship) and standard diagnostic plots 56 . Based on results from this preliminary model, we then re-fitted the model including variables as parametric (i.e. linear) terms as appropriate.
To control for spatial autocorrelation within and between nesting areas we included the geographic coordinates of the nests as a non-parametric term and the nesting area ID as a random effect (penalized smoothed regression term) 56 .

Data availability
The raw data was supplied as a Supplementary Information File.   Table 4. Working hypotheses and variables used to explain P. unifilis nest removal by humans in the eastern Brazilian Amazon. a Strength of variable support from our information theoretic analysis. Unsupported = "−", weakly supported "+", and strongly supported = "++".