Data collection
In the present study, turtle captures from 1981 to 2022 occurred mostly in the lagoon and near the channels of Aldabra (Fig. 1), but tended to be concentrated on the western part of the atoll, near the research station. Sampling sessions were opportunistic and not standardized spatially or temporally. Turtles were either hand-captured by walking through tidal pools and partially exposed seagrass and macroalgae mats at low tide, or by jumping onto or next to individuals from a boat during high tide (i.e., ‘rodeo’ technique) (Eckert et al. 1999). For logistical reasons, captures by boat occurred in sandy areas of the lagoon near the research station, with several focal capture sites within this area (Fig. 1).
Sampling locations were determined by sea conditions and tide height. The first turtle sighted at the location was pursued until caught, or a new turtle followed if the first escaped capture. Once captured, turtles were tagged on the trailing edges of both front flippers (or existing tags recorded), following standard procedures (Balazs 1999). Monel (National Band & Tag Co. styles 681, 49) and plastic yellow Roto-tags (Dalton Supplies Ltd.) were used from the 1980s through mid-1990s, followed by Inconel tags (National Band & Tag Co. style 681) from the mid-1990s to present, with larger Titanium tags (Aust. ‘Turtle’/ Titanium Stockbrands) applied to larger turtles (> 15 kg before 2020, > 10 kg since 2020). Carapace measurements were recorded including minimum (CCLmin) and notch-to-tip (CCLn-t) (as defined by Bolten 1999), and both straight carapace (SCL) and curved carapace lengths (CCL) were taken.
In cases where CCLn-t measurements were not taken and CCLmin measurements existed for an individual, a CCLmin conversion factor was used to obtain the CCLn-t (cm), to maximize capture measurements (Table S1). All length measurement comparisons are reported in CCLn-t. Since measurements were taken by several observers across the decades, measurement error was calculated by comparing the consistency in repeat measurements on the same individuals that were collected within two weeks of each other (Braun-McNeill et al. 2008); mean measurement absolute error was calculated by taking the difference between measurements and calculating absolute mean. It was estimated as 0.52 ± 1.4 cm SD (-4.7–0.8 cm; n = 15).
Data analysis
Analyses were performed in RStudio Version 1.4.1106 and R 4.0.4 (R Core Team 2021). Generalized Additive Models (GAMs) were implemented with package ‘mgcv’ (Wood 2011). Thin plate regression splines and smoothness parameters were estimated by restricted maximum likelihood, REML (Wood 2011; Pedersen et al. 2019). To assess the K (the number of knots used in the fitted splines), the functions Gam.check and k.check were used in the mgcv package (Pedersen et al. 2019). GAMs and confidence intervals were plotted with LOESS smooths using R package ‘ggplot2’ (Wickham 2016). All data were inspected for normality with a Shapiro-Wilk test and non-parametric analyses were used when normality was not met.
Size distributions were estimated for first time captures only. Study day for each capture was defined as the number of days since the beginning of the growth dataset. Measurements resulting in negative growth, due to either measurement error or carapace damage, were included to avoid bias (Bjorndal et al. 2016, 2017).
To mute seasonal peaks or troughs related to condition and minimize growth measurement errors, a recapture dataset was created using a minimum interval of 11 months between captures (Bjorndal et al. 2017). If the period between recaptures was less than 11 months, the next chronological recapture was used. When available, multiple recapture intervals of the same individual were used. Annual growth rate was calculated as the change in CCL during the growth interval, divided by the number of years in the growth interval (Casale et al. 2009; Colman et al. 2015; Bjorndal et al. 2017; Bellini et al. 2019).
To assess if growth rate has changed over time, the 30–40 cm and 40–50 cm CCL size classes, for which sufficient samples sizes were available for a meaningful temporal comparison, were analyzed with year. Growth rate was investigated as a function of year (mean year was binned into years) and turtle ID was used as a random, independent variable to account for repeated measurements of the same individual and for individual-specific heterogeneity: gam(GR ~ s(year), random = list(id ~ 1), method = “REML”, family = gaussian) turning it into a generalized additive mixed model (GAMM). A Kruskal Wallis test with pairwise comparisons was used if the GAMM was significant.
We explored the relationship between age and size in three steps. First, the time spent by turtles at Aldabra was estimated. This was done by investigating growth rate through a GAMM as a function of the following continuous covariates: mean size, mean sampling year, and recapture interval (RI; in years). ID was used as a random, independent variable, as done when assessing growth rate over time (above). The equation was: gam(GR ~ s(CCL) + s(year) + s(RI), random = list(id ~ 1), method = “REML”, family = Gaussian; see also Casale et al. 2009; Colman et al. 2015; Bellini et al. 2019). GAMMs were run separately for the two species. Size specific predictions of the growth rate function were extracted (with 95% confidence intervals) for each capture-mark-recapture (CMR) record, ordered by ascending CCL values. (i.e., the first and last records are the smaller and larger turtles, respectively).
Since it is unknown how old turtles were before recruiting to Aldabra as a foraging ground, we estimated the age of turtles, since recruitment, with year 0 associated with the smallest median size for each species in the data. The integration equation, y(CCLi) = y(CCLi−1) + (CCLi – CCLi−1)/ r(CCLi), was used, where: y(CCLi) represents the age at the initial mean size, y(CCLi−1) is the age of the previous recorded mean size, and r(CCLi) is the individual annual growth rate (Colman et al. 2015; Bellini et al. 2019). The integration was also performed to obtain confidence intervals from the GAMMs. The smallest median size for each species was used as the starting CCLi−1. The last predicted value represented the time spent at Aldabra.
Second, the age of turtles at Aldabra was estimated by adding the estimated time spent at Aldabra (calculated above) to the age of recruitment, which was unknown. To estimate the age of turtles when they first recruit to Aldabra, we compiled figures from the literature. For green turtles, recruitment to neritic habitat has been considered 2–7 years for 30–40 cm CCL (Brazil; Lenz et al. 2017) or 6–9 years for 35–47 cm CCL (south Great Barrier Reef; Limpus and Chaloupka 1997). Using skeletochronology and straight carapace length (SCL), it was estimated as 5–6 years in Hawaii (Zug et al. 2002), 3–6 years in Florida (Zug and Glor 1998) and 1–7 years on the US east coast (Goshe et al. 2010). For hawksbill turtles, recruitment was estimated at 4.5–5.5 years for 33.5 cm SCL in Ascension (Putman et al. 2014; Weber et al. 2017), 1–3 years for 20–25 cm SCL in the Virgin Islands (Boulon Jr 1994), or 2–4 years with SCL by skeletochronology in Hawaii (Snover et al. 2013). Ranges of 2–8 years and 2–5 years were used as the age of recruitment (to be associated with the smallest captured sizes per species) to Aldabra for green turtles and hawksbill turtles, respectively, in this study.
To attempt an age at sexual maturity, we estimated the time needed to grow from the largest observed size class in our dataset to the mean adult nesting size. Published adult growth rates in peer-reviewed literature from Aldabra (for green turtles) and rookeries in the Seychelles (for hawksbill turtles) were used in the model to estimate age at sexual maturity. For green turtles, we used the mean nesting size (108.9 cm CCLn-t) on Aldabra from 1995–2016 (n = 4635; SIF unpubl. data, same dataset from Mortimer et al. 2022, and a female adult growth rate of 0.14 cm year-1 from Aldabra nesters (Mortimer et al. 2022) in the model. For hawksbill turtles, we used the mean nesting size (86.4 cm CCLn-t) from Cousine Island, Seychelles (Gane et al. 2020), and an adult growth rate of 0.17 cm year-1 for females > 80 cm CCL (Coral Sea, Queensland, Australia; Bell and Pike 2012).
A size-specific growth rate function for Aldabra with the addition of the adult growth rate was derived through GAMMs with the model: gam(GR ~ s(CCL), random = list(id ~ 1), method = “REML”, family = gaussian). The GAMM-predicted growth rates were then run through an integration equation (Colman et al. 2015, Bellini et al. 2019) from the largest size observed per species to the mean adult nesting size. Due to a gap in body length measurements between the largest immature captures in the Aldabra dataset and the mean nesting size, the equation was modified; growth rates from the GAMM were predicted at intervals of 1-cm CCL instead of by the mean CCL. The integration equation was applied to a size range of 67–109 cm CCL for green turtles and 73–87 cm CCL for hawksbill turtles. The adult growth rate (as above) was conservatively assumed to be the minimum growth and any estimated growth rates with values lower than the adult growth rate (only applicable for hawksbill turtles in the 84–86 cm CCL range) were substituted with the adult growth rate. An estimate of age of sexual maturity was calculated by adding this time (from the largest size in the dataset to the mean nesting size) to the estimated age at the largest size in the dataset.