Models of primary sex ratios at a major flatback turtle rookery show an anomalous masculinising trend
© Stubbs et al.; licensee BioMed Central Ltd. 2014
Received: 1 July 2014
Accepted: 29 July 2014
Published: 22 October 2014
Quantifying primary sex ratios is essential for assessing how global warming will influence the population dynamics of species with temperature-dependent sex determination (TSD). Process-explicit (mechanistic) models can accurately estimate primary sex ratios but require the resolution of the key physiological parameters that influence sex determination and validation of the model by testing predictions against empirical data.
To address these goals, we conducted incubation experiments on flatback sea turtle (Natator depressus) embryos from a large winter-nesting rookery at Cape Domett in the East Kimberley region of Western Australia. A TSD model fitted to laboratory and field nest data indicated that the pivotal temperature producing equal sex ratios was 29.4°C, with males produced below 27.7°C and females produced above 31.1°C. Back-switch experiments revealed that the thermosensitive period (TSP), when gonads differentiate into testes or ovaries, occurs between 43% and 66% of development to hatching. Integrating this new information with sand temperatures reconstructed from 23 years of historical climate data shows that male-biased sex ratios are likely if the TSP falls during the Austral winter. Annual variation in the simulated sand temperatures increased from 1990 to 2013, with cooler winters producing conditions that favoured male hatchlings for longer periods. The same model projected to 2030 and 2070 suggests that female-biased primary sex ratios will become more prevalent over time.
Our results show that accurate modelling of primary sex ratios depends on quantifying the thermal biology of embryos and on parameterising mechanistic models of sand temperatures with site-specific climate data.
KeywordsTemperature-dependent sex determination Thermosensitive period Flatback turtle Natator depressus Sex ratio Sand temperature Climate change Cape Domett Western Australia
In most animals, sex is determined at conception by inherited sex chromosomes  but in many reptiles, the sex of an embryo is determined by the temperatures experienced during incubation (temperature-dependent sex determination or TSD) . Primary sex ratios are of particular interest in the context of climate change, as a warming climate will increase nest temperatures and could create or exacerbate existing sex ratio biases in hatchling cohorts. Persistent biases in the primary sex ratio can potentially lead to demographic collapse or localised extinction, as has recently been shown in an island population of tuatara .
Marine turtles are a major lineage of reptiles, and all extant species have TSD. In each species, females are produced at temperatures above a threshold and males below . There is an increasing focus on modelling primary sex ratios in these taxa . In part, this is due to the challenges of fieldwork at remote rookeries and also to difficulties in sexing hatchlings using non-invasive methods . Moreover, primary sex ratios vary within and across populations , and measuring sex ratios at one site, or over a single nesting season, will not accurately reflect the primary sex ratio of a population over time. Methods for accurately predicting hatchling sex ratios in the complex thermal environment of natural nests  are needed to reduce the need for broad-scale collection of hatchings from natural nests and for predicting sex ratios under future climates.
Earlier attempts to predict hatchling sex ratios in reptiles with TSD were based on correlations between sex ratios and air or nest temperatures (e.g. –). Correlative models have a key constraint in that they should not be extrapolated to predict sex ratios when temperatures exceed the range used to generate the model , and this constraint can limit their utility for predicting the impacts of global warming. More recently, mechanistic (or ‘process-explicit’) models have been used to estimate nest temperatures for species with TSD ,. A further innovation has been to link predicted nest temperatures to a developmental model of an embryo, where key traits such as hatchling sex, development time and heat-induced mortality can be estimated after quantifying the physiological responses of an embryo to temperature . An advantage of the mechanistic approach is that behavioural responses to counter the effects of global warming, or processes such as metabolic heating, can be simulated by adjusting model inputs or constructs.
Both correlative and mechanistic models require knowledge of the pivotal temperature (Tpiv—the temperature that produces a 1:1 sex ratio) and the transitional range of temperatures (TRT) where both sexes can be produced for the population under study. Mechanistic models further require estimation of the rate of embryonic development across the broad range of current and future environmental temperatures, as well as identification of the period of development where indifferent gonads are sensitive to temperature, known as the thermosensitive period (TSP). In the reptile species so far examined, the TSP broadly falls within the middle half to middle third of incubation, and in most species, it has not been resolved any further ,. Delineating the TSP more precisely allows the identification of the particular portion of the thermal profile of a nest that influences sex determination. In general, a subset of nest temperatures will be required for sex ratio estimation, and selecting the correct subset enhances the accuracy with which nest sex ratios can be predicted.
The flatback turtle (Natator depressus, Garman) is endemic to Australian waters and is listed as ‘data deficient’ by the International Union for the Conservation of Nature ,. While pivotal temperatures have been estimated for both Western and Eastern Australian populations of N. depressus–, other TSD parameters are poorly resolved. Cape Domett in the East Kimberley region of Western Australia hosts one of Australia’s largest N. depressus rookeries where females have an extended nesting season from April to November that peaks in the Austral winter . Most other populations of N. depressus show a nesting peak in summer; hence, a winter nesting population provides an important contrast when considering the relative vulnerability of flatback populations to climate change.
2.1 Incubation times and embryonic development rate
Of 296 eggs incubated in the laboratory under a range of thermal regimes (see ‘Methods’ section), 156 hatched and almost all eggs that did not hatch showed no signs of development (confirmed by dissection). All 156 hatchlings plus 138 late-stage embryos sampled from field nests were sexed using histological techniques, and the classification of sex was highly repeatable (binomial test; p < 0.0001, probability of the same classification = 0.97).
2.2 Thermosensitive period (TSP)
Design of the back-switch experiments used to delineate the thermosensitive period
Starting day of switcha
Length of switcha(d)
Proportion developed at start
Proportion developed at end
28°C–28°C–28°C; M control
32.5°C–32.5°C–32.5°C; F control
2.3 Pivotal temperature (Tpiv) and transitional range of temperatures (TRT)
Variation in the Tpiv and TRT based on data from different subsets of field nests
Nests used to fit model
Lower bound of TRT (°C)
Upper bound of TRT (°C)
Mean residual ± SE (°C)
0.067 ± 0.02
0.102 ± 0.03
0.091 ± 0.03
0.094 ± 0.03
0.081 ± 0.02
0.108 ± 0.04
0.108 ± 0.03
0.082 ± 0.03
0.080 ± 0.03
0.084 ± 0.03
2.4 Validation of the sand temperature model
2.5 Validation of the physiological model
CTEs and sex ratios predicted for field nests based on measured and estimated sand temperatures
Sex of embryos sampled from nest
Range of Constant Temperature Equivalents (CTEs, °C) based on:
Predicted sex ratio based on:
Measured nest temps
Modelled sand temps
Weather station sand CTEs
AWAP sand CTEs
2.6 Historical sand temperatures and sex ratios
2.7 Sand temperatures and sex ratios under climate change
Sand properties required by NicheMapR’s microclimate model
Sand Depth (cm)
Bulk density (g/cm3)
Saturated water content (cm3/cm3)
Mineral thermal conductivity (W/mK)
Mineral heat capacity (J/KgK)
Mineral density (g/cm3)
% sand moisture
Cape Domett hosts the major winter nesting aggregation of N. depressus in Northern Australia . Long-term data on nesting phenology does not exist for this rookery, but in the single year that this was assessed (April 2006–April 2007), peak nesting occurred in August and September . Our sand temperature models show that from 1990 to 2013, the rookery experienced progressively cooler incubation temperatures in August and September (Figure 6), which are likely to have led to the production of predominantly male hatchings. This is contrary to a warming (and thereby feminising) trend that has been either documented for sea turtle beaches in recent years or else has been estimated from correlative models ,,,. Our anomalous result is explained by Cape Domett falling within a small region of north-western Australia where air temperatures have declined relative to historical benchmarks, due to large increases in rainfall and increase in atmospheric aerosols . This trend is captured in the AWAP climate database that forced our model of historical sand temperatures (Figure 1). A change in peak nesting to spring months could compensate for the recent cooling of sand temperatures in winter. Once long-term data on nesting seasonality becomes available for this rookery, our mechanistic model could be used to assess whether phenological shifts would neutralise a male sex ratio bias (e.g. ).
Sand temperatures at 50 cm depth reconstructed over a 23-year period (1990–2013; Figure 6) showed similar seasonal patterns to the two recent years that we depict in greater detail (2011–2012; Figure 7). In most years, sand temperatures showed annual variation of approximately 10°C (Figure 6). As female N. depressus produce an average of 2.8 clutches per season and have a renesting interval of 12–16 days ,, most females should produce offspring of both sexes in each year. Based on the new physiological data reported here, the TSP begins approximately 22 days after oviposition at the Tpiv of 29.4°C, suggesting that the peak nesting period identified in 2006–2007  would have produced predominantly male or mixed sex hatchlings (Figure 6). As our field work occurred during September and October 2012, it was likely to have occurred towards the end of the peak nesting period. Sand temperatures are several degrees centigrade higher in these months than they are in August (Figure 7), and most embryos sampled from nests laid at this time were female (Table 3). Our sand temperatures and sex ratio predictions in 2012 showed good agreement with our empirical data (Table 3).
A trend evident in Figure 6 was an increase in the annual range of sand temperatures, driven primarily by the decreased winter temperatures. This pattern led to a greater proportion each nesting season where sand temperatures favour the production of male hatchlings. If these trends continue without a consequent shift in the timing of peak nesting, then primary sex ratios are likely to become more male-biased at this rookery, contrary to our original hypothesis. The increasing sand temperature range over the 23-year period also coincided with warmer summer temperatures in later years (e.g. 2005–2012; Figure 6D,E). These recent summer sand temperatures exceed the upper threshold of ~35°C for survival of marine turtle embryos , while temperatures predicted under future climates peak at 38.3°C in December (Figure 8D). Hence, year-round nesting would come under strong negative selection due to heat-induced mortality of embryos. Cooler winter and warmer summer temperatures could result in spring and autumn months becoming the most suitable for producing viable hatchlings of both sexes. Potentially, under future climates, this population could shift from predominantly winter nesting to bimodal nesting with peaks in both spring and autumn.
3.1 Validation of sand temperature predictions
The validity of the predicted sand temperatures is strengthened by comparison with actual measurements of nest temperatures at equivalent depths (Additional file 1: Figure S1 and Figure 5). NicheMapR’s microclimate model produced accurate estimates of nest temperatures (without any data-based calibration) across a wide range of depths, except towards the end of incubation when metabolic heating is likely to have increased nest temperatures above that of the surrounding sand. For example, in green (Chelonia mydas) and loggerhead turtle (Caretta caretta) nests, a temperature differential between 0.07°C and 2.86°C has been attributed to metabolic heating ,. A divergence in measured and predicted temperatures towards the end of incubation would not affect our sex ratio estimate, as metabolic heating should be most pronounced after the TSP. For example, in a direct comparison of sand and nest temperatures at identical depths in a study of N. depressus in the Pilbara, Box (2010) concluded that only two of seven nests showed any detectable sign of metabolic heating . Nevertheless, it would be advisable to include a metabolic heating process (e.g. ) into a mechanistic model if traits such as hatchling mortality or incubation duration were of interest.
We did not examine the influence of fine-scale variation in nest orientation, shade or sand properties on nest temperatures in this study, as the Cape Domett rookery is fairly homogenous in these respects. The location of the nest within the beach profile would likely produce the most sand temperature variation due to variable sand moisture and potential shading from vegetation. These parameters can all be modified and mapped spatially using the NicheMapR system (e.g. ). Rapid changes in beach topography can make fine-scale modelling of nest shade and exposure challenging, but high-resolution, high frequency satellite and LIDAR imaging could ultimately allow the dynamic rendering of beach temperatures and sex ratios.
3.2 Delineation of the thermosensitive period (TSP)
The precise delineation of the TSP is important for determining the sequence of incubation temperatures that contribute to sex determination . The TSP is the least researched of all TSD parameters and most reptile studies delineate the TSP no more precisely than the middle half to the middle third of development ,. We showed that the TSP of N. depressus spans a period of about 23% of embryonic development (43%–66%), supporting our hypothesis that the TSP would be less than 33%–50% of the total development period. Hewavisenthi and Parmenter  provide the only other estimate of the TSP in N. depressus, concluding that it occurs between days 32 and 40 of incubation at 26°C. Based on incubation temperatures reported in the Hewavisenthi and Parmenter study, and the development rate model developed here (Figure 2), the proportion of development that would have occurred on days 32 and 40 would have been ~44% and ~54%, respectively. Both studies therefore infer that the TSP starts at about the same developmental stage. However, this study infers that the TSP is longer; encompassing approximately 23% rather than 10% of embryonic development shown in the Hewavisenthi and Parmenter study.
All hatchlings (except one) produced from female–male–female (FMF) treatments were female, indicating that either a masculinising switch does not fall within the TSP we identified, or else that our switch periods were too short or not sufficiently cool to trigger the development of males. The switch periods were designed to encompass the same proportions of development for switches in both directions. Instead, development of N. depressus embryos from Cape Domett occurred at slightly faster rates than for embryos collected from the Pilbara . This caused a greater proportion of development to be completed by the start of the FMF switches than for the corresponding male–female–male (MFM) switches.
Switches to 32.5°C over as little as only 9% of development produced a feminising effect (Figure 3; treatment 5). Experiments on the pond turtle (Emys orbicularis) show that bipotential gonads develop into testes in the absence of a feminising temperature trigger during the TSP . In contrast, when this trigger occurs, the enzyme aromatase acts on gonads to stimulate the production of oestrogen and development of ovaries ,,. In N. depressus, we hypothesise that the TSP is similar for triggering testicular or ovarian development, but that the proportion of development at feminising temperatures may be 9% (or less) to trigger ovarian development. Hence, if gonads are to develop into testes, the majority of the TSP needs to occur at a masculinising temperature. This is a possible reason for the lack of masculinisation seen in the FMF back-switch experiments (treatments 7–11), as experimental switch periods may have been too short for the majority of the TSP to occur at the masculinising temperature. Further experimentation is required to test this hypothesis, and its resolution will have important implications for estimating hatchling sex ratios in natural nests that fluctuate between male- and female-producing temperatures.
3.3 TSD parameters and embryonic development rate
Identifying the Tpiv, TRT and embryonic development rate were not major foci of this study, but the ability of our physiological model (Figure 1) to provide good estimates of hatchling sex ratios at Cape Domett (Table 3) suggests that these parameters and the boundaries of the TSP were reasonably accurate. The Tpiv estimated in this study of 29.4°C for the Cape Domett population is similar to the 29.3°C estimated by Limpus (1995, cited in ) and the 29.5°C estimated for Eastern Australian populations . It is, however, lower than the 30.1°C estimated for a population of N. depressus from the Pilbara region of Western Australia . Differences in the Tpiv estimates for the two Western Australian populations are concordant with the genetic differentiation reported between summer nesting populations from the Pilbara and winter nesting populations in Northern Australia .
We used sex ratio data from field nests to improve our estimates of the Tpiv and TRT, but field data is not an ideal substitute for carefully controlled laboratory experiments. In addition, we have not assessed sex ratios close to the estimated Tpiv (29.4°C). Consequently, constant temperature incubation experiments at temperatures at and near the Tpiv are needed to further refine the TSD parameters for the Cape Domett population. Incubation experiments at higher temperatures would also be highly relevant when considering impacts of global warming, as development rate increases with temperature only up to a certain point, after which it rapidly decreases due to enzyme deactivation . High nest temperatures are a documented cause of mortality in marine turtle nests , and well-resolved thermal response curves will be essential for predicting embryonic survival under climate change.
3.4 Implications for population viability and management
Marine turtles have adapted to past climate change and may also adapt to accelerated rates of climate change associated with anthropogenic global warming . Adaptations to counter increasing temperatures include adjusting the timing of nesting, altering nest characters such as shade or depth and altering nesting latitude ,. As peak nesting at Cape Domett coincides with the coolest sand temperatures, it constrains opportunities for females to temporally adjust peak nesting times to counteract increases in incubation temperatures. Another behavioural characteristic—nest site choice with respect to shade cover—has shown some degree of plasticity in other reptiles, where females select sites of differing shade availability in different climates ,. The effectiveness of this response may be limited at Cape Domett, as shaded nest sites are scarce. Alternatively, females could dig deeper nests, where sand temperatures are usually cooler and more stable , but nest depth is limited primarily by limb length in turtles, meaning that nest depths are relatively fixed relative to those of other reptiles . Similarly, cooler nest locations low on the beach risk flooding and suffocation of embryos.
The heritable variation in TSD traits presents an alternative avenue for adaptation to increasing temperatures. Microevolution of the Tpiv and/or TRT in response to climate change is possible ,, and, unlike behavioural adaptations which alter temperatures experienced by embryos, changes in the Tpiv or TRT would alter the sex ratio produced at a given temperature. In general, data on the heritability of physiological traits in reptiles is scarce and as such provides an important avenue for future research . Mechanistic models such as those developed here readily allow the effectiveness of microevolutionary change to be examined.
The overwhelming consensus that anthropogenic climate change will increase global temperatures means that female biases in primary sex ratios could become more prevalent in marine turtles ,. Consequently, population viability may decline and in extreme cases, localised extinctions could occur . Projected global warming in Queensland suggests that within 50 years, ratios of one male to four females are possible for many marine turtle rookeries in the Great Barrier Reef , and such biases may not be sustainable when combined with other anthropogenic stressors such as habitat modification and loss of rookeries ,. However, polyandrous mating systems may allow marine turtles to persist under female-biased operational sex ratios (e.g. ). Understanding responses to operational sex ratio biases and accurately predicting primary sex ratios will both be important when assessing population viability under a warming climate.
Our reconstruction of cooling in sand temperatures at Cape Domett in recent decades (Figure 6) contrasts with our projections of warmer sand temperatures under future climates (Figure 8). While our projections of future sand temperatures could be further refined by applying downscaled global climate models that anticipate changes in rainfall and cloud cover for this region , it is clear that warmer sand temperatures would lead to more female hatchlings being produced at Cape Domett relative to the primary sex ratios produced historically. Furthermore, the periods of year when sand temperatures remain below an upper limit for embryonic viability are reduced appreciably. If the trend of an increasing range in annual sand temperature continues (as seen in the 23-year reconstruction in Figure 6), then the winter sand temperatures that currently coincide with peak nesting could decrease further, potentially favouring male hatchlings. The fact that Cape Domett falls within a small region of Australia where annual minimum air temperatures have decreased  means that monitoring and conservation of this rookery should be a key priority, as rookeries that produce predominantly males can balance a feminising trend elsewhere in the nesting distribution.
5.1 Ethics Statement
All procedures were reviewed and approved by the University of Western Australia animal ethics committee (UWA animal ethics permit RA/3/100/1145) and the Department of Parks and Wildlife (scientific licence SF008844).
5.2 Study site and egg collection
The N. depressus rookery at Cape Domett occurs on a gently sloping, 1.9 km-long, north-west facing beach in the East Kimberley region of Western Australia, Australia (14.798° S, 128.415° E). Freshly oviposited eggs were collected for transport to a laboratory in Perth, Western Australia on September 20 2012. Approximately, whole clutches (52 eggs) were collected from each of six nesting females (total 296 eggs) and stored in damp vermiculite within a portable refrigerator (model: Engel MT45FP) set to 8°C to suspend embryonic development thereby facilitating egg viability during transport (Harry & Limpus 1989). Eggs were transported by helicopter to Kununurra where they awaited further transport via a commercial airline to Perth. Inadvertently, eggs were held over in Kununurra due to a lack of cargo space, and consequently eggs reached the University of Western Australia about 100 h after oviposition. This time period fell outside the 72-h window of viability defined by Harry and Limpus .
In situ incubation temperatures were recorded in 11 nests on the nights of September 16 and 17 in 2012. Four Thermochron® iButtons (Maxim Integrated Products; DS1921H; accuracy ±1°C, precision 0.125°C) were placed at intervals into each nest chamber during oviposition, with the goal being to record temperatures at different depths within the nest chamber. The nests were situated at various locations within the ‘pink sand’ sections of the beach  and were marked with a GPS device (Garmin eTrex Vista HCx) and photographed for ease of relocation. A weather station (WeatherHawk, Signature Series 232) was erected midway along the beach on a small dune above the high tide mark. Weather data (air temperature, humidity, wind speed, solar radiation) were recorded at 30-min intervals.
The rookery was revisited 40 days later (25–26 October 2012) and all nests were relocated and excavated. Between 12 and 16 eggs that were closest to each of the four iButtons were removed from each nest (total of 138 from all nests). The depths of iButtons below the sand surface were measured to the nearest 10 mm before removal. All remaining eggs were reburied. The embryos inside the sampled eggs were euthanised by chilling for approximately 24 h. Chilling was achieved on the beach using ice, and eggs were relocated to a portable refrigerator set to 2–4°C within 2–3 h of collection. Chilling was maintained during transport by boat and road back to Kununurra. There, dead embryos were removed from their shells and fixed in 10% buffered neutral formalin for road transport to Perth. Ten of the eleven nests (nest IDs 1, 2, 3, 4, 5, 6, 10, 12, 13 and 15) contained viable embryos; eggs from the remaining nest showed no signs of development.
5.3 Incubation experiments
The six live clutches returned to Perth were either used in ‘back-switch’ experiments designed to delineate the TSP, or else were incubated at a constant temperature close to the estimated pivotal temperature. Ten eggs (2 each from 5 clutches) that were not used in back-switch experiments were incubated at a constant temperature of 30.1°C (Model i180), which was the estimated pivotal temperature for a population of N. depressus from the Pilbara region of Western Australia .
For the back-switch experiments, one egg from each clutch was buried in damp sand within one of 48 partially sealed containers. Containers were placed in one of four Steridium incubators (models: i180 and i500) set to either a male-producing temperature (28°C; 24 containers) or a female-producing temperature (32.5°C; 24 containers). During development, each container of 5–6 eggs was switched to the opposite temperature treatment for either 10% or 15% of embryonic development (96–240 h, depending on the time and temperature), after which time they were switched back to their original temperature (Table 1). Each back-switch treatment was replicated within a pair of incubators (one pair of 180 L incubators and one pair of 500 L incubators), and iButtons were placed next to an egg when each container underwent a switch to record the temperatures experienced by the embryos.
We used a development rate function developed for N. depressus by Box (2010) to calculate the day that embryos should have completed 35% of development to hatching when incubated at either 28°C or 32.5°C. Switches then began on this day (or when embryos were further developed—refer Table 1) and concluded when a further 10% or 15% of development to hatching should have been completed at the switch temperature (28°C or 32.5°C). Procedural controls were also included for the earliest switches, in which case embryos incubated at 28°C were switched to a different incubator at 28°C for the appropriate period, and then back again to the original incubator. A procedural control was also applied to eggs incubated at 32.5°C.
Eggs in both the back-switch and Tpiv experiments were weighed at the start of incubation and thereafter every seven days to assess embryonic viability. Unviable eggs (those that lost weight and became discoloured) were removed from incubation boxes. On each occasion, containers were repositioned randomly in the incubators to minimise any impact of subtle temperature gradients within incubators. Containers were checked daily after 35 days of incubation (one week before the earliest expected hatch date) and incubation time was recorded as the time from the start of incubation to when pipping commenced. The four days between oviposition and the start of incubation were not included in the incubation time, as eggs were cooled during this period and embryonic development was suspended . Once pipping commenced, hatchlings were weighed and euthanised by an intracoelomic injection of 0.4 mL sodium pentabarbitone at a concentration of 160 mg/kg. Hatchlings were labelled and stored in 10% buffered neutral formalin for later dissection.
5.4 Identification of hatchling sex
The gonads of marine turtle hatchlings are small (<500 μm) and attached to the kidney . Consequently, the left kidney of each preserved individual was removed. Kidneys were then prepared as paraffin-embedded sections and stained with haematoxylin and eosin for light microscopy . The sex of each specimen was identified based on the differentiation of gonadal medulla and cortex, where males have seminiferous tubules in the medulla and a regressed cortex, and females have a disorganised medulla and a thick, well-developed cortex ,. Each specimen was classified as male, female or unknown on three separate occasions, without reference to previous assessments, and a repeatability analysis was conducted to determine the reliability of the classification. Any specimens where a gonad was not visible or was unable to be classified were resectioned and reexamined until each individual could be classified as male or female.
5.5 Calculation of a non-linear function for embryonic development rate
The parameters fitted by DEVARA describe the maximum development rate (b1), its corresponding temperature (b3) and the temperature at which development rate approaches zero (b2). Parameters b4 and b5 control the asymmetry and steepness of the curve and were fixed at 6 and 0.4, respectively, as is recommended when development rates at extreme incubation temperatures are unknown .
5.6 Delineation of the thermosensitive period (TSP)
Using the parameters estimated for the non-linear development rate function defined above, hourly temperature records from the back-switch experiments were converted into developmental increments. These increments were integrated to determine the cumulative proportion of development that was completed on each day. The length of each switch window was thereafter described as a developmental proportion (e.g. 12% of development). If any back-switch regime produced a sex ratio different to that expected from the dominant incubation temperature, then the TSP was assumed to fall within the portion of development where the temperature was switched.
5.7 Determination of the pivotal temperature (Tpiv) and transitional range of temperatures (TRT)
The relationship between incubation temperature and sex ratio was calculated using TSD software version 4.0.3 [http://max2.ese.u-psud.fr/epc/conservation/TSD/index.html] developed by . This program is widely used to compare the fit of up to five threshold models using maximum likelihood (Richards/a-logistic, Weibull, Weibull*, Hill and Hill*) where * indicates the a-logistic or ‘asymmetrical’ version of the model. Data from constant temperature incubation were intended to be used for this analysis, but the program requires that three or more temperatures produce mixed sexes in order for the fit of different models to be compared . Of the five constant incubation temperature regimes applied in this study (Tpiv experiment, male back-switch controls and female back-switch controls), only one produced mixed sexes. Hence, sex ratio and temperature data from field nests were added to the laboratory data to produce a more rigorous estimate of a sex ratio function.
As the TSD software is designed for constant temperature data, temperature records from field nests were converted to constant temperature equivalents (CTEs; ). Effectively, the CTE is the temperature above which half of development occurs (i.e. a developmental median) and thereby the relationship between development rate and temperature must first be established. Furthermore, when incubation temperatures are variable, the CTE during the TSP (and not the entire incubation period) is most relevant for fitting the TSD function . Hence, we used our development rate function (Figure 2) to convert the nest temperatures that fell within the TSP into hourly developmental increments. These increments were ranked and integrated to determine the temperature above which half of development occurs (the CTE).
where sr(t) is the sex ratio at a given temperature t, P is the pivotal temperature, S describes the steepness of the transition from male to female producing temperatures and K is a parameter that describes the asymmetrical shape of the function. The fitted values of P, S and K and the corresponding TRT were noted and used as reference points for a further ten models, which were fitted from random subsets of five field nests (Table 2). By using only five of the ten field data points available, we were able to meet our goal of retaining an independent data set to test our mechanistic framework (Figure 1). For each of the ten subsets, there was more than one equation that produced an equally good fit (∆AIC <2); however, the Hill equation consistently produced the lowest AIC value across all subsets. Hence, we used the Hill equation to make a standardised comparison of the TSD equation generated by each of the ten datasets, and used the mean residuals to select a consensus model that could be used to estimate the Tpiv and TRT of the Cape Domett population.
5.8 Sand temperature reconstruction and projection
NicheMapR is an R version of the mechanistic modelling program NicheMapperTM;  and was used to predict hourly sand temperatures by simultaneously solving heat and mass balance equations based on physical properties of beach sand and on regional climate data (e.g. ,). The parameters of the microclimate model were not ‘tuned’ to the observed sand temperatures. Rather, the sand properties for the Cape Domett rookery were estimated from relevant literature (e.g. ,) or were measured empirically (Table 5). Percent sand moisture was measured by determining the wet and dry weights of sand samples collected from a range of depths from each field nest during excavation in October. Solar reflectance of sand sampled from two nests in the wavelength range 300–2,100 nm was measured using two spectrometers (Ocean Optics USB2000 for the UV-visible range and NIRQuest for the NIR range) and two light sources (Ocean Optics PX-2 pulsed xenon light for the UV-visible range and HL-2000 tungsten halogen light for the visible-NIR range), all connected with a quadrifurcated fibre optic. The probe on the end of the fibre optic was held within an ocean optics RPH-1 probe holder at a constant angle (45°) and distance from the surface, and each measurement was expressed relative to a Spectralon 99% white reflectance standard (Labsphere, Inc., North Sutton, NH, USA), and weighted by solar irradiance. There was no significant difference in the solar reflectance of the two samples (unpaired t-test; p = 0.959).
Daily maximum and minimum temperatures, relative humidity and wind speeds were collated from data generated by the weather station we installed at the rookery between September 18 and October 24, 2012. Rainfall data during this period were obtained from the nearest Bureau of Meteorology weather station at Wyndham (http://www.bom.gov.au). These weather values were used as inputs into NicheMapR, and sand temperatures were predicted for user-specified depths (range 36–80 cm) using the parameters in Table 2 and assuming 0% shade. Predicted temperatures at a particular depth were then compared to actual nest temperatures measured at the same depth, using the r2 statistic.
Historical (1990–2013) climate data for Cape Domett (daily maximum and minimum temperatures, relative humidity, rainfall and solar radiation) were obtained from the Australian Water Availability Project (AWAP). Gridded long-term average wind speed data were obtained from Australian National University Climate software package (ANUCLIM; ). Both the AWAP and ANUCLIM databases are derived from interpolated data records from weather stations across Australia. These historical climate data were interpolated at a point approximately 25 km south of the Cape Domett rookery (15.003° S, 128.383° E).
The AWAP climate data, combined with the shade and sand parameters defined earlier, were used as inputs into the microclimate model within NicheMapR to estimate sand temperatures at 50 cm depth for the 23 years from 1990 to 2013. To investigate the influence of climate change on sand temperatures, the 2007 air temperatures from the AWAP database were adjusted in accordance with the Commonwealth Scientific and Industrial Research Organisation’s (CSIRO) projections of future air temperatures for Australia . Under a low emissions scenario, AWAP air temperatures were increased by 0.6°C or 1.8°C for the years 2030 and 2070, respectively. Under a high emissions scenario, air temperatures were increased by 1.5°C or 3.4°C for the same years.
5.9 Model validation
The DEVARA function was used to convert the hourly sand temperatures predicted by NicheMapR (or the actual temperatures measured in nests) into developmental increments for the five nests not used to fit the TSD function. Development was assumed to have started on the date of oviposition; hence, the first developmental increment we calculated was for the hour immediately following oviposition. All the hourly development increments calculated for a particular temperature record were integrated to determine the proportion of development completed on each day and thereby identify the dates that formed the boundaries of the TSP. We calculated a CTE for all temperatures that fell within the TSP (as described earlier) and CTEs were inputted into the TSD function (Figure 4B) to estimate the sex ratio of the embryos sampled at the various depths in each nest. Sex ratios were classified as female (<5% males), mixed (5–95% males) or male (>95% males). Predicted sex ratios were then compared to the actual sex ratios of the sampled embryos, which allowed us to assess how effectively our physiological model estimated sex ratios from hourly temperature records (Figure 1).
temperature-dependent sex determination
transitional range of temperatures
constant temperature equivalent
Australian water availability project
Australian National University Climate software package
We are grateful to the Traditional Owners of Miriuwung Gajerrong land, in particular Lawrence Carter, for allowing access to the rookery and for sharing their knowledge. Thank you also to Leia Howes for assistance with egg collection and marking nests, to Jamie Tedeschi and Lorian Woolgar for help in the laboratory, and to Tom Stewart for preparing and sectioning histological samples.
The research was funded by the School of Animal Biology at the University of Western Australia (UWA) and the Western Australian Marine Science Institution. Support and accommodation for our research was provided by Daryl Moncrieff, Allan Thomson, Ben Corey and boat skippers from the Department of Parks and Wildlife (DPaW) in Western Australia.
- Pieau C, Dorizzi M, Richard-Mercier N: Temperature-dependent sex determination and gonadal differentiation in reptiles. Cell Mol Life Sci 1999, 55: 887–900. 10.1007/s000180050342View ArticlePubMedGoogle Scholar
- Valenzuela N, Lance VA: Temperature-Dependent Sex Determination. Smithsonian Institution Press, Washington D C; 2004.Google Scholar
- Grayson KL, Mitchell NJ, Monks JM, Keall SN, Wilson J, Nelson NJ: Sex ratio bias and extinction risk in an isolated population of tuatara (Sphenodon punctatus) . PLoS One 2014, 9: 1–10.Google Scholar
- Ewert MA, Etchberger CR, Nelson CE: Turtle Sex-Determining Modes and TSD Patterns, And Some TSD Pattern Correlates. In Temperature-Dependent Sex Determination In Vertebrates. Edited by: Valenzuela N, Lance VA. Smithsonian Books, Washington; 2004:21–32.Google Scholar
- Fuller W, Godley B, Hodgson D, Reece S, Witt M, Broderick A: Importance of spatio-temporal data for predicting the effects of climate change on marine turtle sex ratios. Mar Ecol Prog Ser 2013, 488: 267–274. 10.3354/meps10419View ArticleGoogle Scholar
- Delmas V: A mechanistic model of temperature-dependent sex determination in a chelonian: the European pond turtle. Funct Ecol 2007, 22: 84.Google Scholar
- Standora EA, Spotila JR: Temperature dependent sex determination in sea turtles. Copeia 1985, 1985: 711–722. 10.2307/1444765View ArticleGoogle Scholar
- Mitchell NJ, Kearney MR, Nelson NJ, Porter WP: Predicting the fate of a living fossil: how will global warming affect sex determination and hatching phenology in tuatara? P Roy Soc Lond B Bio 2008, 275: 2185–2193. 10.1098/rspb.2008.0438View ArticleGoogle Scholar
- Hays GC, Broderick AC, Glen F, Godley BJ: Climate change and sea turtles: a 150-year reconstruction of incubation temperatures at a major marine turtle rookery. Glob Change Biol 2003, 9: 642–646. 10.1046/j.1365-2486.2003.00606.xView ArticleGoogle Scholar
- Hawkes LA, Broderick AC, Godfrey MH, Godley BJ: Investigating the potential impacts of climate change on a marine turtle population. Glob Change Biol 2007, 13: 923–932. 10.1111/j.1365-2486.2007.01320.xView ArticleGoogle Scholar
- Fuentes M, Porter W: Using a microclimate model to evaluate impacts of climate change on sea turtles. Ecol Model 2013, 251: 150–157. 10.1016/j.ecolmodel.2012.12.020View ArticleGoogle Scholar
- Maulany RI, Booth DT, Baxter GS: Emergence success and sex ratio of natural and relocated nests of Olive Ridley turtles from Alas Purwo National park, East Java, Indonesia. Copeia 2012, 2012: 738–747. 10.1643/CH-12-088View ArticleGoogle Scholar
- Janzen FJ: Climate change and temperature-dependent sex determination in reptiles. Proc Natl Acad Sci U S A 1994, 91: 7487–7490. 10.1073/pnas.91.16.7487View ArticlePubMed CentralPubMedGoogle Scholar
- Booth D, Freeman C: Sand and nest temperatures and an estimate of hatchling sex ratio from the Heron Island green turtle (Chelonia mydas) rookery, Southern Great Barrier Reef. Coral Reefs 2006, 25: 629–633. 10.1007/s00338-006-0135-4View ArticleGoogle Scholar
- Georges A, Beggs K, Young JE, Doody JS: Modelling development of reptile embryos under fluctuating temperature regimes. Physiol Biochem Zool 2005, 78: 18–30. 10.1086/425200View ArticlePubMedGoogle Scholar
- Pendoley K, Bell C, McCracken R, Ball K, Sherborne J, Oates J, Becker P, Vitenbergs A, Whittock P: Reproductive biology of the flatback turtle Natator depressus in Western Australia. End Spec Res 2014, 23: 115–123. 10.3354/esr00569View ArticleGoogle Scholar
- Limpus C: The flatback turtle, Chelonia depressa , in Queensland: reproductive periodicity, philopatry and recruitment. Aust Wildlife Res 1984, 11: 579–587. 10.1071/WR9840579View ArticleGoogle Scholar
- Hewavisenthi S, Parmenter C: Hydric environment and sex determination in the flatback turtle ( Natator depressus Garman) (Chelonia: Cheloniidae). Aust J Zool 2000, 48: 653–659. 10.1071/ZO00049View ArticleGoogle Scholar
- Hewavisenthi S, Parmenter C: Thermosensitive period for sexual differentiation of the gonads of the flatback turtle ( Natator depressus Garman). Aust J Zool 2002, 50: 521–527. 10.1071/ZO02014View ArticleGoogle Scholar
- Box T: A Mechanistic Approach To Predicting Hatchling Sex Ratios Of The Flatback Turtle ( Natator Depressus ). In Honours thesis. The University of Western Australia, Animal Biology, Western Australia; 2010.Google Scholar
- Whiting A, Thomson A, Chaloupka M, Limpus C: Seasonality, abundance and breeding biology of one of the largest populations of nesting flatback turtles, Natator depressus : Cape Domett, Western Australia. Aust J Zool 2009, 56: 297–303. 10.1071/ZO08038View ArticleGoogle Scholar
- Kearney M, Shamakhy A, Tingley R, Karoly D, Hoffmann A, Briggs P, Porter W: Microclimate modelling at macro scales: a test of a general microclimate model integrated with gridded continental-scale soil and weather data. Meth Ecol Evol 2014, 5: 273–286. 10.1111/2041-210X.12148View ArticleGoogle Scholar
- Jones DA, Wang W, Fawcett R: High-quality spatial climate data-sets for Australia. Aust Meteorol Oceanogr J 2009, 58: 233.Google Scholar
- Wibbels T, Hillis-Star Z, Phillips B: Female-biased sex ratios of hatchling hawksbill sea turtles from a Caribbean nesting beach. J Herpetol 1999, 33: 142–144. 10.2307/1565556View ArticleGoogle Scholar
- Hanson J, Wibbels T, Martin R: Predicted female bias in sex ratios of hatchling loggerhead sea turtles from a Florida nesting beach. Can J Zool 1998, 76: 1850–1861. 10.1139/z98-118View ArticleGoogle Scholar
- Rotstayn LD, Jeffrey SJ, Collier MA, Dravitzki SM, Hirst AC, Syktus JI, Wong KK: Aerosol- and greenhouse gas-induced changes in summer rainfall and circulation in the Australasian region: a study using single-forcing climate simulations. Atmos Chem Phys 2012, 12: 6377–6404. 10.5194/acp-12-6377-2012View ArticleGoogle Scholar
- Ackerman RA: The Nest Environment And The Embryonic Development Of Sea Turtles. In The Biology of Sea Turtles. Edited by: Lutz PL, Musick JA, Wyneken J. CRC Press, Florida; 1997:83–106.Google Scholar
- Broderick AC, Godley BJ, Hays GC: Metabolic heating and the prediction of sex ratios for green turtles ( Chelonia mydas ). Physiol Biochem Zool 2001, 74: 161–170. 10.1086/319661View ArticlePubMedGoogle Scholar
- Zbinden J, Margaritoulis D, Arlettaz R: Metabolic heating in Mediterranean loggerhead sea turtle clutches. J Exp Mar Biol Ecol 2006, 334: 151–157. 10.1016/j.jembe.2006.01.021View ArticleGoogle Scholar
- Sandoval S, Gómez-Muñoz V, Gutiérrez J, Porta-Gándara MÁ: Metabolic heat estimation of the sea turtle Lepidochelys olivacea embryos. J Therm Biol 2011, 36: 138–141. 10.1016/j.jtherbio.2010.12.007View ArticleGoogle Scholar
- Woolgar L, Trocini S, Mitchell N: Key parameters describing temperature-dependent sex determination in the southernmost population of loggerhead sea turtles. J Exp Mar Biol Ecol 2013, 449: 77–84. 10.1016/j.jembe.2013.09.001View ArticleGoogle Scholar
- Mrosovsky N, Pieau C: Transitional range of temperature, pivotal temperatures and thermosensitive stages for sex determination in reptiles. Amphibia-Reptilia 1991, 12: 169–179. 10.1163/156853891X00149View ArticleGoogle Scholar
- Pieau C, Dorizzi M: Oestrogens and temperature-dependent sex determination in reptiles: all is in the gonads. J Endocrin 2004, 181: 367–377. 10.1677/joe.0.1810367View ArticleGoogle Scholar
- Pieau C, Dorizzi M, Richard-Mercier N, Desvages G: Sexual differentiation of gonads as a function of temperature in the turtle Emys orbicularis : endocrine function, intersexuality and growth. J Exp Zool 1998, 281: 400–408. 10.1002/(SICI)1097-010X(19980801)281:5<400::AID-JEZ5>3.0.CO;2-SView ArticlePubMedGoogle Scholar
- Blamires S, Guinea M: The Influence Of Temperature On Egg Mortality, Emergence Success And Hatchling Sex Ratio For Flatback Sea Turtles ( Natator Depressus) At Fog Bay, Northern Territory, Australia. In Sea Turtles of the Indo-Pacific: Research Management and Conservation. Edited by: Pilcher N, Ismail G. ASEAN Academic Press, London; 2000.Google Scholar
- Pittart S: Genetic Population Structure Of The Flatback Turtle ( Natator depressus ): A Nuclear And Mitochondrial DNA Analysis. In Honours thesis. Institute for Applied Ecology, The University of Canberra; 2010.Google Scholar
- Matsuzawa Y, Sato K, Sakamoto W, Bjorndal K: Seasonal fluctuations in sand temperature: effects on the incubation period and mortality of loggerhead sea turtle ( Caretta caretta ) pre-emergent hatchlings in Minabe. Japan Mar Biol 2002, 140: 639–646. 10.1007/s00227-001-0724-2View ArticleGoogle Scholar
- Hamann M, Limpus C, Read MA: Vulnerability Of Marine Reptiles To Climate Change In The Great Barrier Reef, Climate Change And The Great Barrier Reef A Vulnerability Assessment. In Climate Change and the Great Barrier Reef a vulnerability assessment. Edited by: Johnson JE, Marshall PA. Great Barrier Reef Marine Park Authority and Australian Greenhouse Office, Australia; 2007.Google Scholar
- Hawkes LA, Broderick AC, Godfrey MH, Godley BJ: Climate Change And Marine Turtles. End Spec Res 2009, 7: 137–154. 10.3354/esr00198View ArticleGoogle Scholar
- Refsnider JM, Janzen FJ: Behavioural plasticity may compensate for climate change in a long-lived reptile with temperature-dependent sex determination. Biol Conserv 2012, 152: 90–95. 10.1016/j.biocon.2012.03.019View ArticleGoogle Scholar
- Doody JS, Guarino F, Georges A, Corey B, Murray G, Ewert MW: Nest site choice compensates for climate effects on sex ratios in a lizard with environmental sex determination. Evol Ecol Res 2006, 20: 307–330. 10.1007/s10682-006-0003-2View ArticleGoogle Scholar
- Hulin V, Delmas V, Girondot M, Godfrey M, Guillon J-M: Temperature-dependent sex determination and global change: are some species at greater risk? Oecologia 2009, 160: 493–506. 10.1007/s00442-009-1313-1View ArticlePubMedGoogle Scholar
- Mitchell NJ, Janzen FJ: Temperature-dependent sex determination and contemporary climate change. Sex Devel 2010, 4: 129–140. 10.1159/000282494View ArticleGoogle Scholar
- Broderick AC, Godley BJ, Reece S, Downie JR: Incubation periods and sex ratios of green turtles: highly female biased hatchling production in the eastern Mediterranean. Mar Ecol Prog Ser 2000, 202: 273–281. 10.3354/meps202273View ArticleGoogle Scholar
- Fuentes M, Hamann M, Limpus CJ: Past, current and future thermal profiles of green turtle nesting grounds: implications from climate change. J Exp Mar Biol Ecol 2010, 383: 56–64. 10.1016/j.jembe.2009.11.003View ArticleGoogle Scholar
- Jensen M, Abrue-Grobois F, Frydenberg J, Loeschcke V: Microsatellites provide insight into contrasting mating patterns in arribada vs. non‐arribada olive ridley sea turtle rookeries. Mol Ecol 2006, 15: 2567–2575. 10.1111/j.1365-294X.2006.02951.xView ArticlePubMedGoogle Scholar
- Harry JL, Limpus CJ: Low-Temperature Protection of Marine Turtle Eggs During Long-Distance Relocation. Wildl Res 1989, 16: 317–320. 10.1071/WR9890317View ArticleGoogle Scholar
- Miller J, Limpus C: Ontogeny Of Marine Turtle Gonads. In The Biology of Sea Turtles. Volume 2 edition. Edited by: Lutz P, Musick J, Wyneken J. CRC Press, Boca Raton, FL; 2003:199–224.Google Scholar
- Ceriani S, Wyneken J: Comparative morphology and sex identification of the reproductive system in formalin-preserved sea turtle specimens. Zoology 2008, 11: 179–187. 10.1016/j.zool.2007.07.007View ArticleGoogle Scholar
- Miller J, Limpus C: Incubation Period and Sexual Differentiation in the Green Turtle Chelonia mydas L. In Proceedings Of The Melbourne Herpetological Symposium. Edited by: Banks C, Martin A. The Royal Melbourne Zoological Gardens, Melbourne; 1980.Google Scholar
- Dallwitz MJ, Higgins JP: Devar: A Computer Program For Estimating Development Rate As A Function Of Temperature. In 1993rd edition. ; 1992., [http://delta-intkey.com/devar/]
- Girondot M: Statistical description of temperature-dependent sex determination using maximum likelihood. Evol Ecol Res 1999, 1: 479–486.Google Scholar
- Georges A: The Influence Of Fluctuating Temperatures On Hatchling Sex Ratios - A Model And Proposed Test Using Caretta Caretta . In Proceedings of the Australian Marine Turtle Conservation Workshop; Gold Coast. Edited by James R. Australian Nature Conservation Agency; 1994:156–162.
- Porter WP, Mitchell JW: System and method for calculating the spatial-temporal and other effects of climate and other environmental conditions on animals. US Patent 7 2006, 155: 377.Google Scholar
- Marshall T, Holmes J, Rose C: Soil Physics. Cambridge University Press, New York; 1996.View ArticleGoogle Scholar
- Smits K, Sakaki T, Linsuwat A, Illangasekare TH: Thermal conductivity of sands under varying moisture and porosity in drainage-wetting cycles. Vadose Zone J 2010, 9: 172–180. 10.2136/vzj2009.0095View ArticleGoogle Scholar
- Houlder DJ, Hutchinson MF, Nix HA, McMahon JP: ANUCLIM User Guide, Version 5.0. In Book ANUCLIM User Guide, Version 5.0. Centre for Resource and Environmental Studies, Australian National University, Canberra; 1999.Google Scholar
- Climate Change in Australia, Technical Report In Climate Change in Australia, Technical Report. CSIRO and the Australian Bureau of Meteorology, Canberra; 2007.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.