Acessibilidade / Reportar erro

Seasonal factors affecting sea turtle nesting in the Southeastern Caribbean Sea (Gulf of Paria, Venezuela)

Abstract

The nesting characteristics (number of nests and eggs, time of year, nesting initiation, and nesting length) of leatherback (Dermochelys coriacea) and hawksbill (Eretmochelys imbricata) sea turtles of the southern Caribbean Sea (specifically in the Gulf of Paria in Venezuela), were examined in association with weekly precipitation averages and number of rainy days per week during the period between 2009 and 2018. We hypothesized about the influence of rainfall intensity and patterns as the main abiotic factor for sea turtle nesting. On average, leatherbacks preferred nesting during the drier season of each year (March, April, and May), while hawksbills nested during the rainy season (June to September). For both species, we found few significant correlations between the number of nests or clutch size (number of eggs per nest) and weekly averages of seasonal precipitation rates in the region. Average hawksbill clutch sizes were not correlated with average precipitation rates but were positively correlated with the number of rainy days per week (r=0.66, P≤0.05). Average hawksbill clutch sizes decreased each year on average (-3.3 eggs/year, r=-0.88, P≤0.001), which coincided with a negative long-term trend in the number of rainy days (-0.11 rainy days/week, r=-0.69, P≤0.05). During the study period, nesting activities for both leatherback and hawksbills started progressively later (0.9 and 0.6 weeks/year, respectively p≤0.05) and were shorter (-0.9 and -0.8 weeks /year, P≤0.1 and P≤0.05, respectively).

Descriptors:
Rainfall; Sea turtles; Remote sensing; Seasonal nesting; Fecundity

INTRODUCTION

Sea turtles need specific environmental conditions for nesting in order to increase the chances of successful incubation of the eggs they bury in sandy beaches (Mortimer 1990MORTIMER, J. A. 1990. The influence of beach sand characteristics on the nesting behavior and clutch survival of green turtles (Chelonia mydas). Copeia, 1990(3), 802-817., Wood & Bjorndal 2000WOOD, D. W. & BJORNDAL, K. A. 2000. Relation of temperature, moisture, salinity, and slope to nest site selection in loggerhead sea turtles. Copeia, 2000(1), 119-128, Pike 2013PIKE, D. A. 2013. Climate influences the global distribution of sea turtle nesting. Global Ecology and Biogeography, 22(5), 555-566.). Turtles seek specific beach sand characteristics, but nesting is also affected by the seasonal conditions such as temperature and rainfall (Maulany et al. 2012MAULANY, R. I., BOTH, D. T. & BAXTER, G. S. 2012. The effect of incubation temperature on hatchling quality in the olive ridley turtle, Lepidochelys olivacea, from Alas Purwo National Park, East Java, Indonesia: Implications for hatchery management. Marine Biology, 159, 2651-2661, DOI: https://doi.org/10.1007/s00227-012-2022-6
https://doi.org/10.1007/s00227-012-2022-...
, Rafferty et al. 2017RAFFERTY, A., JOHNSTONE, C. P., GARNER, J. A. & REINA, R. D. 2017. A 20-year investigation of declining leatherback hatching success: implications of climate variation. Royal Society Open Science, 4(10), 170196, DOI: http://doi.org/10.1098/rsos.170196
http://doi.org/10.1098/rsos.170196...
). Understanding this dependence on specific environmental conditions is important for planning management strategies under future climate scenarios (Patricio et al. 2021PATRICIO, A., HAWKES, L. A., MONSINJON, J. R., GODLEY, B. J. & FUENTES, M. P. 2021. Climate change and marine turtles; recent advances and future directions. Endangered Species Research, 44, 362-395, DOI: https://doi.org/10.3354/esr01110
https://doi.org/10.3354/esr01110...
).

Sand moisture and precipitation are likely to be important factors in determining when loggerhead (Caretta caretta) and hawkbill (Eretmochelys imbricata) sea turtles nest (McGehee 1990MCGEHEE, M. A. 1990. Effects of moisture on eggs and hatchlings of loggerhead sea turtles (Caretta caretta). Herpetologica, 46, 251-258., Hitchins et al. 2003HITCHINS, P., BOURQUIN, O., HITCHINS, S. & PIPER, S. 2003. Factors influencing emergences and nesting sites of hawksbill turtles (Eretmochelys imbricata) on Cousine Island, Seychelles, 1995-1999. Phelsuma, 11, 59-69.). Hatchling success for leatherback turtle (Dermochelys coriacea) clutches is negatively affected by intense rainfalls (Houghton et al. 2007HOUGHTON, J., MYERS, A. E., LLOYD, C., KING, R. S., ISAACS, C. & HAYS, G. C. 2007. Protracted rainfall decreases temperature within leatherback turtle (Dermochelys coriacea) clutches in Grenada, West Indies: ecological implications for a species displaying temperature dependent sex determination. Journal of Experimental Marine Biology and Ecology, 345(1), 71-77.), while leatherback hatchling success is also affected by rainfall on a global scale (Santidrian et al. 2015SANTIDRIAN, P., SABA, V., LOMBARD, C. D., VALIULIS, J. M., ROBINSON, N. J., PALADINO, F. V., SPOTILA, J. R., FERNANDEZ, T., RIVAS, M. L., TUCEWK, J., NEL, R. & ORO, D. 2015. Global analysis of the effect of local climate on the hatchling output of leatherback turtles. Science Reports, 5, 16789, DOI: https://doi.org/10.1038/srep16789
https://doi.org/10.1038/srep16789...
). Saragoça et al. (2020)SARAGOÇA, R. B., RESTREPO, J. A. & VALVERDE, R. A. 2020. Effects of El Niño Southern Oscillation and local ocean temperature on the reproductive output of green turtles (Chelonia mydas) nesting at Tortuguero, Costa Rica. Marine Biology, 167(9), 128, DOI: https://doi.org/10.1007/s00227-020-03749-z
https://doi.org/10.1007/s00227-020-03749...
found a positive correlation between precipitation during February (dry season) between 1987 and 2015 and increased green turtle (Chelonia mydas) clutch numbers in Tortuguero, Costa Rica. This positive relationship was accentuated by higher/lower precipitation and larger/smaller clutch numbers during El Niño/La Niña events.

Rain directly affects sand temperature, refreshing and cooling the substrate (McGehee 1990MCGEHEE, M. A. 1990. Effects of moisture on eggs and hatchlings of loggerhead sea turtles (Caretta caretta). Herpetologica, 46, 251-258.). A consequence of temperature change can be alterations in the hatchling sex ratio of loggerhead turtles (Caretta caretta) (Lolavar & Wyneken 2015LOLAVAR, A. & WYNEKEN, J. 2015. Effect of rainfall on loggerhead turtle nest temperatures, sand temperatures and hatchling sex. Endangered Species Research, 28(3), 235-247., Lolavar & Wyneken 2017LOLAVAR, A. & WYNEKEN, J. 2017. Experimental assessment of the effects of moisture on loggerhead sea turtle hatchling sex ratios. Zoology, 123, 64-70, DOI: http://dx.doi.org/10.1016/j.zool.2017.06.007
http://dx.doi.org/10.1016/j.zool.2017.06...
). Heavy rainfalls can lead to lower incubation temperatures, as observed in a South Pacific green turtle rookery which produced 54% male hatchlings between season 2007 until season 2019 (Laloë et al. 2020LALOË, J. O., MONSINJON, J., GASPAR, C., TOURON, M., GENET, Q., STUBBS, J., GIRONDOT, M. & HAYS, G. C. 2020. Production of male hatchlings at a remote South Pacific green sea turtle rookery: conservation implications in a female dominated world. Marine Biology, 167, 70, DOI: https://doi.org/10.1007/s00227-020-03686-x
https://doi.org/10.1007/s00227-020-03686...
). A similar pattern was observed for hawksbills in the northern Australia Great Barrier Reef (Staines et al. 2020STAINES, M. N, BOOTH, D. T., MADDEN, C. A. & HAYS, G. C. 2020. Impact of heavy rainfall events and shading on the temperature of sea turtle nests. Marine Biology, 167, 190, DOI: https://doi.org/10.1007/s00227-020-03800-z
https://doi.org/10.1007/s00227-020-03800...
).

In this work, we hypothesized about the possible influence of rainfall intensity and patterns on the number of nests and clutch size for hawksbill and leatherback sea turtles along the southern shores of the Paria Peninsula (Venezuela), in the southeastern Caribbean Sea. The prediction was that certain pluviometry variables correlate with nest numbers and clutch size for each species. This study complements an earlier study on the impact of human activities on nesting variability around the periphery of the Gulf of Paria, in which only nest numbers were considered (Balladares & Quintero-Torres 2019BALLADARES, C. & QUINTERO-TORRES, A. 2019. Is a small sea turtles rookery doomed to local extinction? Decreasing nesting trends at the Paria Gulf, Venezuela. Marine Ecology, 40(5), e12562, DOI: https://doi.org/10.1111/maec.12562
https://doi.org/10.1111/maec.12562...
). Finally, here we consider only rainfall, not temperature, as the primary abiotic factor affecting nesting.

METHODS

Study Area

The Gulf of Paria is an estuarine region in the southeastern Caribbean Sea, with the Orinoco River delta to the south. It is located south of the Paria peninsula, between Venezuela and Trinidad (Figure 1; Rincon et al. 2008RINCÓN, F., ASTOR, Y., MULLER-KARGER, F., VARELA, R. & ODRIOZOLA, A. 2008. Características oceanográficas del flujo en Boca de Dragón, Venezuela. Memorias de la Fundación La Salle de Ciencias Naturales, 168, 7-24.). Along the southern border of the Paria peninsula, there are several small beaches that are regularly used by sea turtles for nesting (Figure 1; Balladares & Dubois 2010). These beaches are: 1) Macurito, 2) Los Garzos, 3) Silvano, 4) Obispo, and 5) Cerezo. The beaches are narrow, as the coast features a steep terrain. Typical beaches are approximately less than 300 m long and 50 m wide. The fives beaches sampled are sandy with low berm, and covered with supra-littoral vegetation of manchineel (Hippomane mancinella L.) and portia trees (Thespesia populnea L.). Beaches 2, 3, 4, and 5 are inside the Peninsula de Paria National Park. Macurito beach is outside the park and near the fishing town of Macuro.

Figure 1
Location of the main nesting beaches along the southeastern Caribbean Sea in the Gulf of Paria, Venezuela. The beaches monitored in this study were (1) Macurito, (2) Los Garzos, (3) Silvano, (4) Obispo, and (5) Cerezo.

Field work

Surveys began in 2003, from mid-March though the end of October, and were conducted each year through 2021. Each site was visited four times per week during the day, and during the night (between 8 pm and 4 am) at Los Garzos beach. Beach patrols were performed with at least two local trained patrollers assigned to the main nesting grounds of playa Los Garzos and Macurito. All females emerging from the sea to nest were identified by species, measured, and photographed. Turtles were also tagged when possible. Special care was taken in counting eggs in new nests each night. Starting in August 2006, once a proper and secure hatchery was built in the town of Macuro, we moved all nests from their original beaches to the hatchery to avoid predation and poaching. By 2007, we had expanded the search to all five beaches (Balladares & Dubois 2014BALLADARES, C. & DUBOIS, E. 2014. Saqueo y depredación de nidadas de tortugas marinas, durante las temporadas 2003 a 2012, en seis playas del Golfo de Paria, Venezuela. Cuadernos de Investigación UNED, 6(2), 239-243.). During 2012, 2017, and 2018, we were unable to translocate the eggs due to maintenance and expansion of the hatchery. During those years, we intensified patrolling of the beaches. Natural and translocated nests were counted for clutch size and emergence success. Hatchlings were then set free on the beaches where the eggs had been laid if possible, or on Macuro, where the hatchery was located.

Precipitation Datasets and Processing

Precipitation records from regional meteorological stations are irregular and with an unknown uncertainty because of continuously changing observation methods and poor stations maintenance. We therefore used rainfall estimates from the Multi-Satellite Precipitation Analysis (TMPA V7), which merges rainfall estimates from several satellites with gauge data (Huffman et al. 2007HUFFMAN, G. J., ADLER, R. F., CURTIS, S., BOLVIN, D. T. & NELKIN, E. J. 2007. Global rainfall analyses at monthly and 3-h time scales. In: LEVIZZANI, V., BAUER, P. & TURK, J. (ed.). Measuring precipitation from space. Dordrecht: Springer, pp. 291-305.). We specifically downloaded the TMPA-V7 3B42 product from NASA (https://pmm.nasa.gov/data-access/downloads/trmm). This provided daily-accumulated precipitation with a spatial resolution of 0.25 degrees (~28 km). We extracted the precipitation data over the Gulf of Paria, off the distal part of the Paria peninsula (10.625°N - 61.875°W). We calculated weekly precipitation averages and counted the number of rainy days during the week to compare with the weekly sea turtle nest data.

The TMPA-V7 3B42 dataset provides three products: ‘Precipitation’ (Surface Precipitation Estimate), ‘IR-precipitation’ (Infrared Precipitation Estimate), and ‘HQ-precipitation’ (Microwave Precipitation Estimate). In order to identify which product was closest to in-situ precipitation data, we compared weekly averages from each satellite products with the weekly average precipitation from the Piarco Meteorological Station in Trinidad for the period 1998-2015 (10.61°N and 61.35°W) located ~60 km east from the center of the grid cell of the satellite dataset (Piarco meteorological data moved to http://www.metservice.gov.tt/climate_data/Piarco_view.php). The highest correlation was found with ‘Precipitation’ (r = 0.58, R2=0.34, slope=0.43, P<0.001, n=902; ‘IR-precipitation’ and ‘HQ-precipitation’ showed R2 = 0.30 and 0.21, respectively). The weekly averages of TMPA ‘Precipitation’ shows the same seasonality as the in-situ Piarco data (Fig. 2a). The TMPA ‘Precipitation’ averages, however, underestimated the magnitude of in-situ precipitation (weekly satellite precipitation is approximately 63% of the weekly in-situ precipitation at Piarco). Long-term averages and standard deviations for the weekly time series were 126.7 ± 133.6 mm/day for TMPA ‘Precipitation’. The similarity in their seasonality and the high standard deviations of both datasets are the reason why only 10 of 52 climatological weekly averages showed significant differences between the in-situ and satellite precipitation (paired student’s t-test P<0.05, Fig. 2a). The precipitation product TMPA-V7 has been found to systematically underestimate or overestimate rainfall in different climatic regions (Yong et al. 2014YONG, B., CHEN, B., GOURLEY, J. J., REN, L., HONG, Y., CHEN, X., WANG, W., CHEN, S. & GONG, L. 2014. Intercomparison of the Version-6 and Version-7 TMPA precipitation products over high and low latitudes basins with independent gauge networks: Is the newer version better in both real-time and post-real-time analysis for water resources and hydrologic extremes? Journal of Hydrology, 508, 77-87., Qiao et al. 2014QIAO, L., HONG, Y., CHEN, S., ZOU, C. B., GOURLEY, J. J. & YONG, B. 2014. Performance assessment of the successive Version 6 and Version 7 TMPA products over the climate-transitional zone in the southern Great Plains, USA. Journal of Hydrology, 513, 446-456., Guo et al. 2015GUO, H., CHEN, S., BAO, A., HU, J., GEBREGIORGIS, A., XUE, X. & ZHANG, X. 2015. Inter-comparison of high-resolution satellite precipitation products over Central Asia. Remote Sensing, 7(6), 7181-7211.). At other locations close to our study area, such as off French Guiana and the north region of Brazil, the satellite precipitation products also tend to underestimate precipitation. However, the TMPA-V7 product performed better than other satellite precipitation products (Ringard et al. 2015RINGARD, J., BECKER, M., SEYLER, F. & LINGUET, L. 2015. Temporal and spatial assessment of four satellite rainfall estimates over French Guiana and North Brazil. Remote Sensing, 7(12), 16441-16459.). Given these results, we chose the ‘Precipitation’ product for our analyses.

Figure 2
Long term weekly averages of (a) satellite average precipitation and number of rainy days per week over the study area (10.625°N - 61.875°W, 1998-2018), and in-situ precipitation from a meteorological station located 60 km east of the study area (Piarco, Trinidad & Tobago, 1998-2015). Comparison between the weekly time series of satellite and in-situ precipitation (paired Student’s t-test P<0.05) showed only 10 weeks with significant differences (white circles). Weekly averages of the number of nests and the average clutch size per nest (2009-2018) for (b) leatherback, Dermochelys coriacea and (c) hawksbill, Eretmochelys imbricata. For each species, the total time series clutch size average per nest ± 1 SD are shown for comparison with the weekly averages. The smallest clutch size of leather back during the middle of March comes from only one nest during the entire time-series.

Statistical Analyses

For the analysis we used the number of nests per beach (excluding poached or lost) and clutch size for leatherbacks and hawksbills in new nests, counted daily. We analyzed sea turtle data from 2009 to 2018 due to the more systematic sampling of all five beaches during this period. Weekly averages and climatology (2009-2018) were calculated for the number of nests and clutch sizes. Precipitation was averaged weekly for the same decade, and the number of days with rainfall for each week was recorded. Annual seasonal summaries for all variables during the nesting seasons of leatherbacks and hawksbills are presented in Tables 1 and 2.

Table 1
Summary of variables obtained during the leatherback (Dermochelys coriacea) nesting season from mid-March to mid-July. The earliest and latest nesting weeks for the study period defined the “nesting season” (weeks 14th to 43th). Values shown are for each year and for the entire time series (2009-2018). Weekly rainy days averages and precipitation averages (mm d-1) were calculated for the entire season, and for the weeks with and without. (*) Where there were more than one week with the highest seasonal number of nests, the values shown are for the first week.
Table 2
Summary of variables obtained during the hawksbill (Eretmochelys imbricata) nesting season from April to mid-October. The earliest and latest nesting weeks for the whole study period was consider the “nesting season” (weeks 14th to 43th). Values shown for each year and for the whole time series (2009-2018). Weekly rainy days averages and precipitation averages (mm d-1) were calculated for the whole season, and for the weeks with nests and the weeks with no-nests. (*) In case there were more than one week with the highest seasonal number of nests, the values shown are for the first week.

The first and the last nesting weeks of the entire study period were used to define the nesting season for each species. For leatherbacks, this was on average from mid-March to mid-July (weeks 11 to 28), while for hawksbills, this was from April to mid-October (weeks 14 to 43).

Possible associations between each species’ nesting parameters and the precipitation variables were explored using linear partial correlation models in MATLAB (2014)MATHWORKS (US). 2014. MATLAB and Statistics Toolbox Release. Natick: The MathWorks, Inc.. We used partial correlation to check the association between the dependent variable with each predictor after adjusting for the other variables, to avoid misleading results should there be a confounding variable related to both the dependent and independent variables. We analyzed the possible relationships between the number of nests and clutch size (which reflects fecundity) and between these variables and weekly average rainfall and number of rainy days during the week. We used regression and correlation (Zar 2010ZAR, J. H. 2010. Biostatistical analysis. New Jersey: Pearson Prentice Hall., Legendre & Legendre 2012LEGENDRE, P. & LEGENDRE, L. 2012. Numerical ecology. 3rd ed. Amsterdam: Elsevier.) to analyze these relationships, producing estimates for the total time series and for each year studied.

RESULTS

General seasonal pattern for each species and satellite data

The satellite precipitation records (Figure 2a) show a short dry season between February and early April (Weeks 4-15) for the study area. There is a transition period between mid-April to May with increasing precipitation, and a lengthy rainy season between mid-June and November (Weeks 24-49) with some variability. The number of rainy days during the week was proportional to precipitation with similar seasonality: up to 1 rain day per week during the dry season, and 2.5-3 rain days per week during the peak of the rainy season.

Leatherbacks had an earlier and shorter nesting season compared to hawksbills. The weekly nesting climatology (Figure 2b) shows the start of the nesting season for the leatherbacks as early as Week 11 and up to Week 28 (March to mid-July). Leatherback nesting peaks between mid-April and mid-May (Weeks 15-20; ~1.5-2 nests per week). The leatherback nesting season peaks during the transition from dry to rainy season; when precipitation and number of rainy days peaks, leatherback nesting season is ending. The long-term (2009-2018) weekly average number of nests was 2.5 ± 2.1SD per week (Table 1). The clutch size of leatherbacks did not show seasonality, though an apparent higher variability in weekly clutch size averages occurred at the beginning and end of the nest season, when there were fewer nests (Table 1). In general, the weekly averages of leatherbacks clutch sizes per nest are within 110.8 ± 17.7 eggs per nest (Figure 2b and Table 1).

Hawksbill nesting season length was over twice as long as that of leatherbacks (Table 2), lasting from 19 weeks (2018) to 29 weeks (2013), with an average nesting season of 22.8 weeks. On average, ~83% of the weeks during the hawksbill nesting season saw new nests. The start and end of hawksbill nesting season was less variable, with ranges of 8 and 6 weeks, respectively (Table 2). Nesting season began between Week 14 (2010) and Week 22 (2016), with an average starting week of 17.6. The season ended between Week 37 (2011) and Week 43 (2013), with an average ending week of 39.4. On average, new nests were identified during 83% of the weeks during the nesting season. The beginning of hawksbill nesting season coincided with increasing precipitation and number of rainy days during the transition period, and the nesting peak coincided with the first two peaks of the rainy season (Fig. 2c). The weekly climatology shows that hawksbills started nesting as early as Week 14 week and up to Week 43 (mid-March to mid-October, Fig. 2c), with a nesting peak of approximately eight nests per week occurring between June-August (Weeks 22-34). The long-term (2009-2018) weekly average number of nests during hawksbill nesting season was 5.8 ± 4.3 per week (Table 2). Similar to leatherbacks, hawksbill clutch sizes did not show seasonality and showed higher apparent variability at the beginning and end of the nest season, when there were fewer nests (Table 2). However, those weekly averages are within 1.5 standard deviations of the long-term average (154.0 ± 19.6 eggs per nest, Table 2).

The weekly time series of satellite-derived precipitation observations (Figure 3a) showed values oscillating between 0 to 640 mm d-1, with high precipitation variability during the rainy season. No long-term trends for weekly precipitation nor weekly number of rainy days were observed.

Figure 3
Time-series (2009-2018) of (a) weekly averages of satellite precipitation and weekly number of days with rain, and five-week moving-averages for precipitation (red) and rainy days (blue). Number of nests per week and weekly averaged clutch size for (b) leatherback and (c) hawksbill. The clutch size of hawksbill showed a significant negative long-term trend of -3.3 eggs per year (dashed line, R2 0.21, P<0.01, n=191); the other parameters did not show any long-term trend.

Figure 3b and 3c illustrate the annual number of nests and clutch size for 2009-2018. Leatherbacks laid over twice the number of nests during 2009 than during the rest of the series. We did not detect nests during 2017, but there was no unusual precipitation during that period (Figure 3a). The weekly time series of hawksbill clutch size showed a significant negative long-term trend of -3.3 eggs per year (Figure 3c, r= -0.45, R2=0.21, P<0.01, slope= -3.29, n= 191), while none of the other weekly time series parameters showed any significant trend.

Annual correlations of nesting vs. rainfall

For each study year, the partial correlations did not show any significant associations between the number of nests or clutch sizes of leatherbacks with weekly rainfall or number of rainy days (Table 3). During 2009, leatherbacks had an exceptionally long nesting season (12 weeks compared to an average of 5.5 weeks during the rest of the study period); it follows that 2009 was the year with the largest number of nests. However, there was no apparent relationship between number of nests and rainfall in any year. Clutch size and precipitation were only correlated during 2010 (negatively) and 2014 (positively), albeit with opposite signs. Clutch size was significantly correlated to number of rainy days per week during 2012 (negatively) and 2013 (positively). Again the signs of the correlations were opposite in those years.

Table 3
Partial correlations in a probability distribution of four variables: mean daily precipitation in a week, number of rainy days in a week, number of nests per week, and average clutch size per week, for Leatherback (Dermochelys coriacea). The values shown in the cells are the partial correlation coefficients. The asterisks indicate the level of significance of the two-tailed test; * P < 0.1, ** P < 0.05, *** P < 0.01. The values used for the calculations were those of the weeks with nests; the number of these weeks is n. Years 2015 to 2018 did not have enough weeks with nests to calculate partial correlations.

Hawkbill nest count per week and clutch size were positively correlated with precipitation during 2009 and with rainy days during 2010. Clutch size was positively correlated with precipitation only in 2014 and negatively correlated with rainy days (Table 4).

Table 4
Partial correlations in a probability distribution of four variables: mean daily precipitation in a week, number of rainy days in a week, number of nests per week, and average clutch size per week, in the Hawksbill (Eretmochelys imbricata). The values shown in the cells are the partial correlation coefficients. The asterisks indicate the level of significance of the two-tailed test; *P < 0.1, ** P < 0.05, *** P < 0.01. The values used for the calculations were those of the weeks with nests; the number of these weeks is n.

Tables 3 and 4 also show the partial correlations between number of nests and clutch size for each species. There were no significant correlations between those two variables in any season nor for the entire study period. This means that the hypothesis of a relationship between these fecundity parameters and precipitation is not supported.

Seasonal correlations

We applied multiple correlation analysis on the seasonal average precipitation (and number of rainy days) averaged for each species nesting season and for the weeks with nests and the weeks of no nesting (Tables 3 and 4), but found no statistical differences between those averages using the Welch’s t-test, implying that rainy conditions during the rainy seasons were similar for weeks with and without nesting (Table 5a).

Table 5
(a) Partial correlations in a probability distribution of five variables: mean daily precipitation in a week, number of rainy days in a week, length of nesting season, number of nests per week, and average clutch size per week, in the Leatherback (Dermochelys coriacea) and the Hawksbill (Eretmochelys imbricata). The replicates are the mean values of each year. The values shown in the cells are the partial correlation coefficients. The asterisks indicate the level of significance of the two-tailed test; *P < 0.1, ** P < 0.05, *** P < 0.01. The values used for the calculations were the means of those of the weeks with nests; the number of these weeks is n. (b) Also shown are the simple correlations between some variables and year.

The analysis of seasonal parameters showed statistically significant long-term trends for these and other variables (Table 5b). The hawksbill clutch size showed a negative long-term trend (Figure 3c; regression parameters b= -3.35 eggs per year, r= -0.88, P ≤0.01, n=10) when using the complete weekly time series; the regression parameters were b= -3.29, r= -0.45, P ≤0.01, n= 191. Long-term trends calculated with seasonal parameters (i.e., averages for each year’s nesting season), which have fewer data points, rendered practically the same results as the trend calculated with the longer time series.

For both species, the start and duration of nesting season shifted later and shorter, respectively (Tables 1 and 2). For 2009-2018, long-term trends for leatherbacks showed a later start of the nesting season (First week with nests: correlation slope b= 0.90 weeks/year, r= 0.69, P ≤0.05) and a shorter nesting period (Season length: b= -0.85 weeks/year, r= -0.64, p ≤0.1). Similarly, for hawksbills, the nesting seasons shifted to later starts (b= 0.59 weeks/year, r= 0.73, P ≤0.05), and shorter durations (b= -0.76 weeks/year, r= 0.66, P ≤0.05). During the leatherback nesting period, which coincides with the dry season, there were no significant trends for any of the precipitation parameters. However, during the hawksbill nesting season, there was a long-term trend of fewer rainy days per week (b= -0.11 rainy days per week/year, r= -067 P ≤0.05).

DISCUSSION

General seasonal pattern for both species

There is a clear separation of the nesting season of leatherback and hawksbill sea turtles in the Gulf of Paria. Leatherbacks prefer to nest late in the dry season, allowing hatchlings to emerge during the rainy season. This is consistent with reports by Santidrian et al. (2015)SANTIDRIAN, P., SABA, V., LOMBARD, C. D., VALIULIS, J. M., ROBINSON, N. J., PALADINO, F. V., SPOTILA, J. R., FERNANDEZ, T., RIVAS, M. L., TUCEWK, J., NEL, R. & ORO, D. 2015. Global analysis of the effect of local climate on the hatchling output of leatherback turtles. Science Reports, 5, 16789, DOI: https://doi.org/10.1038/srep16789
https://doi.org/10.1038/srep16789...
regarding the greater number of leatherback newborns during increased precipitation conditions in areas with dry climates in four large rookeries worldwide. Conversely, hawksbills nest and emerge during the rainy period. The separation of nesting seasons between the two species could allow for the maximization of nesting on small rookeries such as ours.

Rainfall is deemed a positive factor for sea turtle reproduction (McGehee 1990MCGEHEE, M. A. 1990. Effects of moisture on eggs and hatchlings of loggerhead sea turtles (Caretta caretta). Herpetologica, 46, 251-258., Bezy et al. 2016BEZY, V., GIRONDOT, M. & VALVERDE, R. 2016. Estimation of the net nesting effort of olive Ridley Arribada sea turtles based on nest densities at Ostional Beach, Costa Rica. Journal of Herpetology, 50(3), 409-415.). However, for the period examined, precipitation showed no trends, and hawksbill fecundity decreased during the study period, possibly due to anthropogenic impacts (Balladares & Quintero-Torres 2019BALLADARES, C. & QUINTERO-TORRES, A. 2019. Is a small sea turtles rookery doomed to local extinction? Decreasing nesting trends at the Paria Gulf, Venezuela. Marine Ecology, 40(5), e12562, DOI: https://doi.org/10.1111/maec.12562
https://doi.org/10.1111/maec.12562...
). Our results suggest that nest numbers depend more on human-related impacts and on biotic factors such as feeding and migration than on precipitation patterns.

Seasonal Correlations and Final remarks about Abiotic conditions vs. Fecundity of sea turtles in the region

Our analysis focused on the fecundity metrics nest numbers and clutch size. We did not consider hatching nor emergence success. Regarding this parameter, Olive Ridley turtles Lepidochelys olivacea in Costa Rica (Dornfeld et al. 2015DORNFELD, T., ROBINSON, N. J., SANTIDRIÁN, P. S. & PALADINO, F. V. 2015. Ecology of solitary nesting olive ridley sea turtles at Playa Grande, Costa Rica. Marine Biology, 162(1), 123-139.; Bezy et al. 2016BEZY, V., GIRONDOT, M. & VALVERDE, R. 2016. Estimation of the net nesting effort of olive Ridley Arribada sea turtles based on nest densities at Ostional Beach, Costa Rica. Journal of Herpetology, 50(3), 409-415.), had higher hatchling success during the rainy season than during the dry season. Rafferty et al. (2017)RAFFERTY, A., JOHNSTONE, C. P., GARNER, J. A. & REINA, R. D. 2017. A 20-year investigation of declining leatherback hatching success: implications of climate variation. Royal Society Open Science, 4(10), 170196, DOI: http://doi.org/10.1098/rsos.170196
http://doi.org/10.1098/rsos.170196...
reported an increase in leatherback hatchling mortality with decreasing precipitation at Sandy Point National Wildlife Refuge (SPNWR). However, no such trend in local temperature or rainfall was observed between 1990-2010 at SPNWR, according to the same authors.

We did not find a statistically significant relationship between nest count and clutch size for either leatherbacks or hawksbills. If the number of nests is indicative of nesting females (Robinson et al. 2014ROBINSON, N. J., VALENTINE, S. E., SANTIDRIÁN, P. S., SABA, V. S., SPOTILA, J. R. & PALADINO, F. V. 2014. Multidecadal trends in the nesting phenology of Pacific and Atlantic leatherback turtles are associated with population demography. Endangered Species Research, 24(3), 197-206.), the absence of relation between nest count and size suggests that there is no density-dependence. Said differently, the density of mature females did not affect their fecundity in our study area. Halley et al. (2018)HALLEY, J. M., VAN HOUTAN, K. S. & MANTUA, N. 2018. How survival curve affects populations’ vulnerability to climate change. PloS One, 13(9), e0203124. suggested that leatherbacks follow a Type III life history strategy of high fecundity, low juvenile survival, and increasing survival of later age classes, and the author concluded that organisms with this strategy are more strongly influenced by environmental variability and climatic changes. Halley et al. (2018)HALLEY, J. M., VAN HOUTAN, K. S. & MANTUA, N. 2018. How survival curve affects populations’ vulnerability to climate change. PloS One, 13(9), e0203124. assumed a density-independent model for their study, suggesting that the role of density-dependence is still not well understood. In our case, while it is possible that density effects are present due to intraspecific competition, these are not expressed as a decrease in clutch size per female when density increases.

In the Gulf of Paria, the longest leatherback nesting season seen in 2009 was not related to precipitation variables but instead to a larger number of nesting females. Regarding hawksbills, during 2014 and 2010, significant correlations were found between clutch size and the number of rainy days per week. However, that relationship was not evident when analyzing the seasonal averages of all years together. Both parameters showed long-term trends of smaller clutch sizes and fewer rainy days per week. Although possible relationships between those two parameters are unclear, Balladares et al. (2020)BALLADARES, C., GONZÁLEZ, M. F. & RODRÍGUEZ, D. 2020. A matrix population model for the hawksbill sea turtle (Eretmochelys imbricata) in the Gulf of Paria, Venezuela. LatinAmerican Journal of Aquatic Research, 48(5), 739-748, DOI: https://doi.org/10.3856/vol48-issue5-fulltext-2476
https://doi.org/10.3856/vol48-issue5-ful...
found that larger hawksbill females in the region have more than twice the fecundity of younger reproductive females.

This work aimed to understand how precipitation, which in other regions regulates sea turtle nesting (Lolavar & Wyneken 2015LOLAVAR, A. & WYNEKEN, J. 2015. Effect of rainfall on loggerhead turtle nest temperatures, sand temperatures and hatchling sex. Endangered Species Research, 28(3), 235-247., Bezy et al. 2016BEZY, V., GIRONDOT, M. & VALVERDE, R. 2016. Estimation of the net nesting effort of olive Ridley Arribada sea turtles based on nest densities at Ostional Beach, Costa Rica. Journal of Herpetology, 50(3), 409-415., Saragoça et al. 2020SARAGOÇA, R. B., RESTREPO, J. A. & VALVERDE, R. A. 2020. Effects of El Niño Southern Oscillation and local ocean temperature on the reproductive output of green turtles (Chelonia mydas) nesting at Tortuguero, Costa Rica. Marine Biology, 167(9), 128, DOI: https://doi.org/10.1007/s00227-020-03749-z
https://doi.org/10.1007/s00227-020-03749...
), affected the nesting of leatherback and hawksbill sea turtles in the southeastern Caribbean Sea. Aside from clear seasonality in both species, we could not demonstrate a clear relationship between variation or trends in precipitation and reproductive output during the study period (2009-2018). However, we did find that the nesting seasons for both leatherbacks and hawksbills are tending to start later and to be shorter. The drivers of these trends are worth exploring and could be related to other environmental variables (retaste temperatures or humidity inside the micro nesting environment), or to ecological factors such migrations, interspecific competition, feeding habits, or human activities.

It is important to sustain monitoring of the nesting characteristics of sea turtles to understand whether these patterns reflect short- or long-term changes associated with climate change (McGehee 1990MCGEHEE, M. A. 1990. Effects of moisture on eggs and hatchlings of loggerhead sea turtles (Caretta caretta). Herpetologica, 46, 251-258., Hawkes et al 2009HAWKES, L., BRODERICK, A., GODFREY, M. & GODLEY, B. 2009. Climate change and marine turtles. Endangered Species Research, 7, 137-154.). However, though important to understanding climate impacts, long-term monitoring projects are challenging to maintain given limited resources, especially in less developed countries that provide critical habitats for endangered species such as sea turtles.

CONCLUSIONS

In the region of the Gulf of Paria in Venezuela, there is a clear time separation between the nesting seasons of leatherback and hawksbill sea turtles, with both seasons occurring during rainy months. During the study period (2009-2018) the nest count per year decreased in both species. Although in some years there seems to be some influence of precipitation on clutch size (leatherback) and on nests count and clutch size (hawksbill), the overall data indicate that precipitation does not influence these fecundity parameters in this region of the Southern Caribbean Sea.

ACKNOWLEDGMENTS

The first author (CB) thanks the people of Macuro, especially field assistants Modesto Garcia and Jose Miguel Reinosa, the Venezuelan Environment Ministry, and the financial support of Timberland Venezuela, C.A. during 2016 -2019. CB also acknowledges the effort of coauthors DRR, DR, FMK, and HBG. CB is also grateful for the support provided to this program from: SEE Turtles (Brad Nahill), Global Conservation (Mr. Jeff Morgan), and the Marine Conservation Action Fund from the New England Aquarium (Elizabeth Stephenson, Bess Edwards, and Emily Duwan) in the United States, and to Provita, A.C. (Bibiana Sucre, Juan Carlos Amilibia, Administration team, and Dr. Jon Paul Rodriguez) in Venezuela. Finally, CB acknowledges Dr. Bryan Wallace and Dra. Michelle Maria Early Capistrán, who definitively reviewed this manuscript and the English grammar.

REFERENCES

  • BALLADARES, C. & DUBOIS, E. 2014. Saqueo y depredación de nidadas de tortugas marinas, durante las temporadas 2003 a 2012, en seis playas del Golfo de Paria, Venezuela. Cuadernos de Investigación UNED, 6(2), 239-243.
  • BALLADARES, C., GONZÁLEZ, M. F. & RODRÍGUEZ, D. 2020. A matrix population model for the hawksbill sea turtle (Eretmochelys imbricata) in the Gulf of Paria, Venezuela. LatinAmerican Journal of Aquatic Research, 48(5), 739-748, DOI: https://doi.org/10.3856/vol48-issue5-fulltext-2476
    » https://doi.org/10.3856/vol48-issue5-fulltext-2476
  • BALLADARES, C. & QUINTERO-TORRES, A. 2019. Is a small sea turtles rookery doomed to local extinction? Decreasing nesting trends at the Paria Gulf, Venezuela. Marine Ecology, 40(5), e12562, DOI: https://doi.org/10.1111/maec.12562
    » https://doi.org/10.1111/maec.12562
  • BEZY, V., GIRONDOT, M. & VALVERDE, R. 2016. Estimation of the net nesting effort of olive Ridley Arribada sea turtles based on nest densities at Ostional Beach, Costa Rica. Journal of Herpetology, 50(3), 409-415.
  • DORNFELD, T., ROBINSON, N. J., SANTIDRIÁN, P. S. & PALADINO, F. V. 2015. Ecology of solitary nesting olive ridley sea turtles at Playa Grande, Costa Rica. Marine Biology, 162(1), 123-139.
  • GUO, H., CHEN, S., BAO, A., HU, J., GEBREGIORGIS, A., XUE, X. & ZHANG, X. 2015. Inter-comparison of high-resolution satellite precipitation products over Central Asia. Remote Sensing, 7(6), 7181-7211.
  • HALLEY, J. M., VAN HOUTAN, K. S. & MANTUA, N. 2018. How survival curve affects populations’ vulnerability to climate change. PloS One, 13(9), e0203124.
  • HAWKES, L., BRODERICK, A., GODFREY, M. & GODLEY, B. 2009. Climate change and marine turtles. Endangered Species Research, 7, 137-154.
  • HITCHINS, P., BOURQUIN, O., HITCHINS, S. & PIPER, S. 2003. Factors influencing emergences and nesting sites of hawksbill turtles (Eretmochelys imbricata) on Cousine Island, Seychelles, 1995-1999. Phelsuma, 11, 59-69.
  • HOUGHTON, J., MYERS, A. E., LLOYD, C., KING, R. S., ISAACS, C. & HAYS, G. C. 2007. Protracted rainfall decreases temperature within leatherback turtle (Dermochelys coriacea) clutches in Grenada, West Indies: ecological implications for a species displaying temperature dependent sex determination. Journal of Experimental Marine Biology and Ecology, 345(1), 71-77.
  • HUFFMAN, G. J., ADLER, R. F., CURTIS, S., BOLVIN, D. T. & NELKIN, E. J. 2007. Global rainfall analyses at monthly and 3-h time scales. In: LEVIZZANI, V., BAUER, P. & TURK, J. (ed.). Measuring precipitation from space Dordrecht: Springer, pp. 291-305.
  • LALOË, J. O., MONSINJON, J., GASPAR, C., TOURON, M., GENET, Q., STUBBS, J., GIRONDOT, M. & HAYS, G. C. 2020. Production of male hatchlings at a remote South Pacific green sea turtle rookery: conservation implications in a female dominated world. Marine Biology, 167, 70, DOI: https://doi.org/10.1007/s00227-020-03686-x
    » https://doi.org/10.1007/s00227-020-03686-x
  • LEGENDRE, P. & LEGENDRE, L. 2012. Numerical ecology 3rd ed. Amsterdam: Elsevier.
  • LOLAVAR, A. & WYNEKEN, J. 2015. Effect of rainfall on loggerhead turtle nest temperatures, sand temperatures and hatchling sex. Endangered Species Research, 28(3), 235-247.
  • LOLAVAR, A. & WYNEKEN, J. 2017. Experimental assessment of the effects of moisture on loggerhead sea turtle hatchling sex ratios. Zoology, 123, 64-70, DOI: http://dx.doi.org/10.1016/j.zool.2017.06.007
    » http://dx.doi.org/10.1016/j.zool.2017.06.007
  • MATHWORKS (US). 2014. MATLAB and Statistics Toolbox Release. Natick: The MathWorks, Inc.
  • MAULANY, R. I., BOTH, D. T. & BAXTER, G. S. 2012. The effect of incubation temperature on hatchling quality in the olive ridley turtle, Lepidochelys olivacea, from Alas Purwo National Park, East Java, Indonesia: Implications for hatchery management. Marine Biology, 159, 2651-2661, DOI: https://doi.org/10.1007/s00227-012-2022-6
    » https://doi.org/10.1007/s00227-012-2022-6
  • MCGEHEE, M. A. 1990. Effects of moisture on eggs and hatchlings of loggerhead sea turtles (Caretta caretta). Herpetologica, 46, 251-258.
  • MORTIMER, J. A. 1990. The influence of beach sand characteristics on the nesting behavior and clutch survival of green turtles (Chelonia mydas). Copeia, 1990(3), 802-817.
  • PATRICIO, A., HAWKES, L. A., MONSINJON, J. R., GODLEY, B. J. & FUENTES, M. P. 2021. Climate change and marine turtles; recent advances and future directions. Endangered Species Research, 44, 362-395, DOI: https://doi.org/10.3354/esr01110
    » https://doi.org/10.3354/esr01110
  • PIKE, D. A. 2013. Climate influences the global distribution of sea turtle nesting. Global Ecology and Biogeography, 22(5), 555-566.
  • QIAO, L., HONG, Y., CHEN, S., ZOU, C. B., GOURLEY, J. J. & YONG, B. 2014. Performance assessment of the successive Version 6 and Version 7 TMPA products over the climate-transitional zone in the southern Great Plains, USA. Journal of Hydrology, 513, 446-456.
  • RAFFERTY, A., JOHNSTONE, C. P., GARNER, J. A. & REINA, R. D. 2017. A 20-year investigation of declining leatherback hatching success: implications of climate variation. Royal Society Open Science, 4(10), 170196, DOI: http://doi.org/10.1098/rsos.170196
    » http://doi.org/10.1098/rsos.170196
  • RINCÓN, F., ASTOR, Y., MULLER-KARGER, F., VARELA, R. & ODRIOZOLA, A. 2008. Características oceanográficas del flujo en Boca de Dragón, Venezuela. Memorias de la Fundación La Salle de Ciencias Naturales, 168, 7-24.
  • RINGARD, J., BECKER, M., SEYLER, F. & LINGUET, L. 2015. Temporal and spatial assessment of four satellite rainfall estimates over French Guiana and North Brazil. Remote Sensing, 7(12), 16441-16459.
  • ROBINSON, N. J., VALENTINE, S. E., SANTIDRIÁN, P. S., SABA, V. S., SPOTILA, J. R. & PALADINO, F. V. 2014. Multidecadal trends in the nesting phenology of Pacific and Atlantic leatherback turtles are associated with population demography. Endangered Species Research, 24(3), 197-206.
  • SANTIDRIAN, P., SABA, V., LOMBARD, C. D., VALIULIS, J. M., ROBINSON, N. J., PALADINO, F. V., SPOTILA, J. R., FERNANDEZ, T., RIVAS, M. L., TUCEWK, J., NEL, R. & ORO, D. 2015. Global analysis of the effect of local climate on the hatchling output of leatherback turtles. Science Reports, 5, 16789, DOI: https://doi.org/10.1038/srep16789
    » https://doi.org/10.1038/srep16789
  • SARAGOÇA, R. B., RESTREPO, J. A. & VALVERDE, R. A. 2020. Effects of El Niño Southern Oscillation and local ocean temperature on the reproductive output of green turtles (Chelonia mydas) nesting at Tortuguero, Costa Rica. Marine Biology, 167(9), 128, DOI: https://doi.org/10.1007/s00227-020-03749-z
    » https://doi.org/10.1007/s00227-020-03749-z
  • STAINES, M. N, BOOTH, D. T., MADDEN, C. A. & HAYS, G. C. 2020. Impact of heavy rainfall events and shading on the temperature of sea turtle nests. Marine Biology, 167, 190, DOI: https://doi.org/10.1007/s00227-020-03800-z
    » https://doi.org/10.1007/s00227-020-03800-z
  • WOOD, D. W. & BJORNDAL, K. A. 2000. Relation of temperature, moisture, salinity, and slope to nest site selection in loggerhead sea turtles. Copeia, 2000(1), 119-128
  • YONG, B., CHEN, B., GOURLEY, J. J., REN, L., HONG, Y., CHEN, X., WANG, W., CHEN, S. & GONG, L. 2014. Intercomparison of the Version-6 and Version-7 TMPA precipitation products over high and low latitudes basins with independent gauge networks: Is the newer version better in both real-time and post-real-time analysis for water resources and hydrologic extremes? Journal of Hydrology, 508, 77-87.
  • ZAR, J. H. 2010. Biostatistical analysis New Jersey: Pearson Prentice Hall.

Edited by

Editor: Rubens M. Lopes

Publication Dates

  • Publication in this collection
    03 Oct 2022
  • Date of issue
    2022

History

  • Received
    12 May 2022
  • Accepted
    10 Aug 2022
Instituto Oceanográfico da Universidade de São Paulo Praça do Oceanográfico 191, CEP: 05508-120, São Paulo, SP - Brasil, Tel.: (11) 3091-6501 - São Paulo - SP - Brazil
E-mail: diretoria.io@usp.br