The impact of environmental factors on marine turtle stranding rates

Globally, tropical and subtropical regions have experienced an increased frequency and intensity in extreme weather events, ranging from severe drought to protracted rain depressions and cyclones, these coincided with an increased number of marine turtles subsequently reported stranded. This study investigated the relationship between environmental variables and marine turtle stranding. The environmental variables examined in this study, in descending order of importance, were freshwater discharge, monthly mean maximum and minimum air temperatures, monthly average daily diurnal air temperature difference and rainfall for the latitudinal hotspots (-27°, -25°, -23°, -19°) along the Queensland coast as well as for major embayments within these blocks. This study found that marine turtle strandings can be linked to these environmental variables at different lag times (3–12 months), and that cumulative (months added together for maximum lag) and non-cumulative (single month only) effects cause different responses. Different latitudes also showed different responses of marine turtle strandings, both in response direction and timing.Cumulative effects of freshwater discharge in all latitudes resulted in increased strandings 10–12 months later. For latitudes -27°, -25° and -23° non-cumulative effects for discharge resulted in increased strandings 7–12 months later. Latitude -19° had different results for the non-cumulative bay with strandings reported earlier (3–6 months). Monthly mean maximum and minimum air temperatures, monthly average daily diurnal air temperature difference and rainfall had varying results for each examined latitude. This study will allow first responders and resource managers to be better equipped to deal with increased marine turtle stranding rates following extreme weather events.


Introduction
In recent years, tropical and subtropical regions, such as Queensland, have experienced many extreme weather events [1][2][3][4][5][6][7], including droughts, cyclones and protracted rain depressions. In Australia, during summer there is a heightened risk of extreme weather and warmer temperatures, the summer of 2010/2011 in Queensland is of particular note. During this time, cyclones and protracted rain depressions caused wide-spread flooding which in turn led to increased periods of turbid water and increased nutrient and sediment loads from freshwater run-off being dumped into all four major coastal waterways (Brisbane, Fitzroy, Burnett and Burdekin Rivers) [8]. The cyclones and floods stressed coral reefs and seagrass beds causing large-scale die-off of ecologically important seagrass species and decreased water quality intermittently along the entire length of the Queensland coastline south from Cairns [8][9][10][11]. It was postulated that within a year (short-term) of these types of catastrophes, marine megafauna show an increase in the number of stranding, mortalities and exacerbated poor health conditions [12]. In a similar ilk, it has been shown that environmental variables affect seabird wrecks numbers and locations [13].
The ongoing poor weather conditions recently experienced are unprecedented in the 35 year history of the Great Barrier Reef Marine Park [11] and Queensland in general [14] since European settlement. The magnitude and scale of the bad weather conditions experienced during early 2011 on the Great Barrier Reef have not been seen since recording began in 1918 [11].
Norman et al. [15] stated that the ability to understand and investigate marine mammal unusual mortality events and other unexpected strandings that involve substantial die-offs of the marine mammal population are important events which serve as indicators of ocean health. This may give better insight into larger environmental issues, which can have implications for human health and animal welfare. This One Health paradigm can also be applied to marine turtle strandings as marine turtles have been proposed as sentinels of environmental health [16][17][18] and, as such, an increase in the numbers of animals which strand can indicate that the environments in which they live have changed [18].
It has been suggested that marine turtle stranding numbers follow seasonal trends influenced by weather events as well as land-based and at-sea seasonal activities. There have been links made between extreme weather and increased strandings within 12 month periods as outlined by [12,[19][20][21].
Meager and Limpus [21] stated that the most plausible explanation for the high rate of strandings and mortalities of near shore green turtles during 2011 were extreme weather events that occurred in late 2010 and early 2011, which impacted on seagrass foraging areas. They linked this because most of the examined mortalities were attributed to protracted ill health/poor body condition in green turtles and dugongs; which both primarily forage on seagrass. There was evidence that seagrass pastures, mangrove forests, algal beds and coral reefs in Queensland were impacted by a combination of elevated rainfall, flooding and three cyclones (Category 5 Yasi, Category 2 Anthony and Category 1 Tasha) with a protracted low pressure system during the summer of 2010/2011. This study examined marine turtle stranding rates in relation to environmental variables (including rainfall, freshwater discharge rates and air temperature). Different latitudinal blocks, species and age classes were investigated to determine if there were different responses. We summarized and analyzed the available data to provide first responders and management agencies with information to better assist them when responding to stranding events. The databases used for this study are the most comprehensive databases available for Queensland marine turtle records and was established over 30 years ago.

Stranding data
StrandNet is the Queensland Government's Department of Environment and Heritage Protection (EHP) statewide database which records dead, sick and injured threatened marine animals for the entire coast of Queensland and adjacent Commonwealth waters. Records are received from members of the public, and employees of EHP, Queensland Parks and Wildlife (QPWS), Queensland Department of Agriculture, and Fisheries (DAF) and the Great Barrier Reef Marine Park Authority (GBRMPA). Information is collated and stored in this central database. Once reports are entered by on-ground staff the information available is verified by regional and state coordinators for standardization. Additional data is often obtained after the stranding event from veterinarians, pathologists and other biologists who complete more detailed post-mortem investigations.
As a proxy of age class, standard measurements such as curved carapace length (CCL) and tail to carapace length (TCL) were collected at the time of initial stranding [22]. This data was used to assign turtles into 3 age classes: small immatures, large immature and adult sized.
Sex was determined by gonad examination by trained personnel either onsite or using photographs or measurements [23,24].
The study area encompassed latitude -10.78˚to -28.16˚and longitude 142.15˚to 155 (  Fig 1). This part of the east coast of Queensland was selected as it has a long-term and comprehensive dataset; with data collection biased to regions of survey and higher populations. This limitation is openly acknowledged by Meager and Limpus (2012) but considered valid as a representative of a minimum recovery rate and indicative of trends occurring. As the exact location where a stranding was reported was not necessarily where the impact/incident occurred, strandings were grouped into latitudinal blocks of 1˚to account for this potential error. The areas of focus for this study were the hotspots recognized by Flint et al. [27] as -27˚, -25˚, -23˚and -19˚ (Figs 1 and 2). In addition, major embayments, irrespective of latitudinal blocks were assessed (Fig 2).
Time. The date a turtle was reported stranded was used as a proxy of time of death, grouped to a monthly scale.
Cause of stranding. The term 'stranding' is used here to incorporate all reported sick, injured, incapacitated or dead marine turtles that were either found ashore or, in rare cases, were encountered at sea. It included turtles which were entangled in fishing nets, synthetic debris or rescued from a situation where they would have died had it not been for human intervention [28].
Within StrandNet, the primary cause of death/stranding was identified based on gross examination, photograph and/or necropsy by trained personnel [12,29]. The single cause of stranding identified in StrandNet was based on the summation of information available.

Environmental data
Rainfall, freshwater discharge and air temperature were examined as environmental variables. These were selected as they provided the most comprehensive, readily available and up to date dataset of environmental conditions available. Turbidity, water temperature, pH and salinity were not used due to paucity of current available data along the Queensland coastline.
Freshwater discharge is the amount of freshwater running through a river's gauging (recording) station, measured in cumecs (cubic meter per second, m 3 .s -1

Data analysis
Data from StrandNet was grouped into 1˚latitudinal blocks from the -28˚(Queensland-New South Wales border) north to -16˚(Cape Tribulation) for each month between January 1996 and December 2013 [12]. Only natural and unknown causes of death were used for this analysis, as anthropogenic causes can be seasonal due to increased activity (eg. Fishing and boating) [12]. The "unknown cause" used as the operating practice for StrandNet was applied when there was no obvious cause of trauma or subsequent analysis done [12]. Strandings were also isolated from bays recognized from the Queensland Spatial Catalogue (http://qldspatial.information.qld.gov.au/catalogue/) under the Department of Natural Resources and Mines CC BY license. Bays were selected as representing an encapsulating body of water including some estuarine and tidal habitats, within which a population may usually reside irrespective of arbitrary coordinates. Standings were mapped using ArcGIS and then overlapped with the Bay layer. For embayment assessment, strandings were only used if they occurred within the defined bay area.

Model formation
When constructing the model, environmental discharge, air temperature and rain variables were lagged up to 12 months, with a cumulative effect. Time lag one included the environmental factor from time 0 and time -1, time lag two included the environmental factor from time 0, -1 and -2; and so on. A non-cumulative lag effect was also used for this analysis and compared against the cumulative effect.
A 12-month maximum lag time was used as there has historically been links made between marine turtle and dugong deaths occurring within this time frame of extreme weather events [12,20,30]. As seagrass loss after extreme weather events has been noted to be immeditate it is not through to delay the response observed in marine turtle stranding rates [31].
All species of marine turtle to occur within the study area were analyzed individually and collectively as a total count of strandings.
Age classes used for analysis were large immature, adult sized, small immatures, combined small immatures and large immatures, combined large immature and adult sized as well as all age classes together. Models were analyzed where sample sized allowed.
The latitudes with the most strandings (both embayments and whole blocks) were chosen to run the models. These latitudes were -27˚, -25˚, -23˚, -19˚.

Model computation
The models were run as general linear models using R [32] with the bbmle package used to calculate additional information criterion including weights and qAIC values [33,34]. The models were run a priori approach due to the complexity and number of possible models [35][36][37].
Steps followed were similar to those outlined in Bolker et al [35]. Briefly these were specifying the effects, choosing an error distribution, graphically checking variance, fit GLM model to both full model and with each factor.
The strandings data had an excess number of zeros and data was also over-dispersed, so a quasi-Poisson error distribution was used [38].

Model hypothesis
The hypotheses tested are outlined below: 1. Small minimum air temperature will cause increases in marine turtle stranding rates.
2. Maximum air temperature will not affect marine turtle stranding rates.
3. Increased rainfall will cause increased marine turtle strandings rates 7-9 months later. 4. Increased freshwater discharge will cause increased marine turtle stranding rates 7-9 months later.
5. All environmental factors combined will affect marine turtle stranding rates 7-9 months later.

Model testing
To begin with, models were run with all variables combined. These models proved non-significant (p > 0.1). After this, each environmental factor was run separately to determine the individual effect. This was done for each age class and species for each latitude chosen. A no effect model was also run for each variable, ageclass and species. In order to compare models, QAIC weights were calculated using the relative likelihood of the model. This was done following the steps outlined in Bolker [34], briefly the regular model was fit, then the over dispersion parameter was manually extracted to calculate a qAIC value. qAIC is the quasi Akaike Information Criterion (AIC). qAIC weights allow for the selection of a "best approximating model" [36]. This was then used in conjunction with the significance of the variables to determine which model most accurately explained the variance.
Strandings numbers of less than 10 over the 18-year period were excluded due to small sample size as were age class and species with less than 2 turtles per month for the 18-year period.

Numbers of animals reported stranded
The number of turtles reported stranded over the 18 years is depicted in Table 1.
Upon initial investigation green turtles were the only species which could be analyzed separately due to sample size. For the remaining sections of this study green turtles and the total number of strandings were analyzed and reported.

Green turtles
Rainfall. Table 2 summarizes the relationship between rainfall and green turtle strandings rates. In brief it shows that within the -19˚and -27˚blocks strandings decreased as rainfall increased, while the -23 and -25˚blocks showed split responses; the majority of age classes showed significant responses within the first 3 months; obvious differences between cumulative and non-cumulative effects of rainfall; different responses time noted with both embayments and whole blocks.
QAIC's for all groups assessed were different and no patterns were observed ( Table 2). In most cases, the QAICs corresponded with significant responses, with an exception for the age classes which did not produce a significant relationship.
Within each examined latitudinal block, there were no observed pattern as to which age class was the first to show significant responses (Tables 3 and 4).
These patterns did not change when comparing embayment's with whole blocks but the lag time may be extended when examining strandings within the embayment compared to whole block strandings (Tables 3 and 4

).
All examined latitudinal blocks for non-cumulative lagged effects responded similarly to cumulative effects, with non-cumulative showing responses first (Tables 3 and 4).
QAIC's for all groups assessed were different and no patterns were observed (Tables 3 and  4). In most cases, the QAICs corresponded with significant model responses, with an exception for the age classes which did not produce a significant relationship.
Peak discharge. Table 5 summarizes the relationships between green turtle stranding numbers and peak discharge. In brief, large immatures and large turtles in the -27˚block showed no significant response; in the -19˚, -23˚and -27˚degree blocks, as peak discharge Table 3. Model results for green turtles and cumulative mean discharge. " denotes increased strandings rates with increased discharge. # denotes decreased stranding rates with increased discharge. Age class abbreviations: ALL = all turtles, SI = small immature, LI = large immatures, A = adult-sized, ALL IMM = all immature sized animals (small + large), Large = all large turtles (large immatures + adult-sized). The time frame is reported in monthly ranges where responses were noted. The values reported in qAIC are the months with the most significant qAIC value.

Latitude
Age increased so did green turtle stranding rates, as a comparison, in the -25˚block showed a split response with strandings decreasing with increased discharge over the first 5 months, which then switched to increased strandings with increasing peak discharge. Table 5 also displays that within each examined latitudinal block, most age classes showed a significant stranding response to peak discharge, however, there was no observed pattern as to which age class was the first to show significant responses; all examined latitudinal blocks for non-cumulative lagged effects of peak discharge, responded similarly to cumulative effects, with non-cumulative showing responses first.
These patterns did not change when comparing embayments with whole blocks but the lag time may be extended when examining strandings within the embayment compared to whole block green turtle strandings (Table 5).
QAIC's for all groups assessed were different and no patterns were observed (Table 5). In most cases, the QAICs corresponded with significant responses, with an exception for the age classes which did not produce a significant relationship.
Monthly mean maximum air temperature. Table 6 summarizes the relationship between monthly mean maximum air temperature and green turtle stranding rates. In brief it shows that in most cases, as monthly mean maximum air temperatures increased the green turtle stranding rate decreased; there was a significant response noted within the first 4 Table 4. Model results for green turtles and mean discharge. " denotes increased strandings rates with increased discharge. # denotes decreased stranding rates with increased discharge. Age class abbreviations: ALL = all turtles, SI = small immature, LI = large immatures, A = adult-sized, ALL IMM = all immature sized animals (small + large), Large = all large turtles (large immatures + adult-sized). The time frame is reported in monthly ranges where responses were noted. The values reported in qAIC are the months with the most significant qAIC value.

Latitude
Age months; there were very obvious differences between cumulative and non-cumulative effects of monthly mean maximum air temperature, with non-cumulative effects more likely to produce split responses; and that there were similar stranding response times noted with both embayments and the whole blocks for monthly mean maximum air temperature. QAIC's for all groups assessed were different and no patterns were observed (Table 6). In most cases, the QAICs corresponded with significant responses, with an exception for the age classes which did not produce a significant relationship.
Monthly mean minimum air temperature. Table 7 summarizes the relationship between monthly mean minimum air temperature and green turtle stranding rates. In brief, it shows in most cases, as monthly mean minimum air temperatures increased the stranding rate decreased; there were very obvious differences between cumulative and non-cumulative effects, with non-cumulative effects resulted in split responses; in most cases there was a significant green turtle strandings response noted within the first 3 months of the mean minimum air temperature recorded; there were similar responses time noted with both embayments and whole blocks.
QAIC's for all groups assessed were different and no patterns were observed (Table 7). In most cases, the QAICs corresponded with significant responses, with an exception for the age classes which did not produce a significant relationship. Table 5. Model results for green turtles and peak discharge. " denotes increased strandings rates with increased discharge. # denotes decreased stranding rates with increased discharge. Age class abbreviations: ALL = all turtles, SI = small immature, LI = large immatures, A = adult-sized, ALL IMM = all immature sized animals (small + large), Large = all large turtles (large immatures + adult-sized). The time frame is reported in monthly ranges where responses were noted. The values reported in qAIC are the months with the most significant qAIC value.

Latitude
Age Environmental factors on marine turtle stranding rates Monthly average daily diurnal air temperature difference. Table 8 summarizes the relationship between monthly average daily diurnal air temperature difference and green turtle stranding rates. In brief it shows very obvious differences between cumulative and non-cumulative effects of monthly average daily diurnal air temperature differences; Non-cumulative effects resulted in split responses whereas in most cases the cumulative effects resulted in decreased stranding rate with increased mean minimum air temperature; in most cases there was a significant response noted within the first 3 months of monthly average daily diurnal air temperature difference being recorded; similar response times were noted with both embayments and whole blocks.
The exception was the -19˚block, where significant response times were varied for small immatures, immature and all green turtles and adults and large turtles within the -19˚block did not display a significant response (Table 8).
QAIC's for all groups assessed were different and no patterns were observed. In most cases, the QAICs corresponded with significant responses, with an exception for the age classes which did not produce significant relationship (Table 8).

All marine turtle strandings
Rainfall. Table 9 summarizes the relationship between rainfall and marine turtle stranding rates. In brief it shows that when comparing rainfall across all blocks, there were different Table 6. Model results for green turtles and monthly mean maximum air temperature. " denotes increased strandings rates with increased mean maximum air temperature. # denotes decreased stranding rates with increased mean maximum air temperature. Age class abbreviations: ALL = all turtles, SI = small immature, LI = large immatures, A = adult-sized, ALL IMM = all immature sized animals (small + large), Large = all large turtles (large immatures + adult-sized). The time frame is reported in monthly ranges where responses were noted. The values reported in qAIC are the months with the most significant qAIC value. patterns noted for each block; cumulative effects within the -27˚block stranding rates decreased as rainfall increased; non-cumulative effects within the -27˚block showed mixed results; within each examined latitudinal block, there were similar stranding response times noted for embayments and the whole blocks; there were very obvious differences between cumulative and non-cumulative effects of rainfall on all turtle stranding rates ( Table 9). QAIC's for all groups assessed were different and no patterns were observed ( Table 9). In most cases, the QAICs corresponded with significant responses, with an exception for the age classes which did not produce a significant relationship.

Latitude Age class
Cumulative mean and mean discharge. Tables 10 and 11 summarizes the relationship between cumulative mean, mean discharge and all marine turtle stranding rates. In brief, it shows in most cases, as cumulative mean discharge increased, the stranding rate for all turtles also increased; each examined latitudinal block, there were similar stranding response times noted for embayments and the whole blocks in respect to discharge; within each examined latitudinal block, there were also similar response times for cumulative effect vs non-cumulative effect of discharge.
Monthly mean maximum air temperature. Table 13 summarizes the relationship between monthly mean maximum air temperature and marine turtle stranding rates. In brief, it shows most cases split responses were observed; when a split response was not noted stranding rates decreased and monthly mean maximum air temperature increased; there were similar response times noted for embayments and the whole blocks for monthly mean maximum air temperature; there were also similar response times for cumulative effect vs non-cumulative effect of monthly mean maximum air temperature.
QAIC's for all groups assessed were different and no patterns were observed ( Table 13). In most cases, the QAICs corresponded with significant responses, with an exception for the age classes which did not produce a significant relationship.

Latitude Age class Cum-Whole Cum-Bay Non Cum-Whole Non Cum-Bay
Monthly average daily diurnal air temperature difference. Table 15 summarizes the relationship between monthly average daily diurnal air temperature difference and marine turtle stranding rates. In brief it shows that in most cases, a split response was observed for monthly average daily diurnal air temperature difference; adults and large turtles from the whole -19˚block did not show significant stranding responses; within each examined latitudinal block, there were similar response times noted for embayments and the whole blocks.
QAIC's for all groups assessed were different and no patterns were observed (Table 15). In most cases, the QAICs corresponded with significant responses, with an exception for the age classes which did not produce a significant relationship. Table 10. Model results for all turtles and cumulative mean discharge. " denotes increased strandings rates with increased discharge. # denotes decreased stranding rates with increased discharge. Age class abbreviations: ALL = all turtles, SI = small immature, LI = large immatures, A = adult-sized, ALL IMM = all immature sized animals (small + large), Large = all large turtles (large immatures + adult-sized). The time frame is reported in monthly ranges where responses were noted. The values reported in qAIC are the months with the most significant qAIC value.

Discussion
This is the first study of its kind to elucidate the effects of individual environmental variables on the stranding rates of coastal marine turtle populations and provides a baseline for future predictive models that can be used as real-time management tools. We found that strandings occurred after a lag phase, with water discharge having the greatest effect on stranding numbers. This study found that the cumulative effects of freshwater discharge in all latitudes resulted in increased strandings 7-12 months later (Tables 3-5, Tables 10-12). The cumulative effects of mean maximum and minimum air temperature resulted in decreased stranding rates immediately through to a lag of 9 months (Tables 6 and 7, Tables 13 and 14). Monthly average daily diurnal air temperature difference resulted in increased strandings immediately through to a lag of 8 months (Tables 8-15). There was an overall decrease in stranding rates 2-8 months after high rainfall events, although the relationship was less clear (Tables 2-9). When comparing cumulative effects against non-cumulative effects, non-cumulative effects were more likely to produce split or dual responses (Tables 2-15). This could be due to cumulative effects having a more lasting, stronger effect. The cumulative effect of multiple months of increased discharge and rainfall potentially does not allow time for the seagrasses to recover and hence have a stronger effect on marine turtles through their diets. Table 11. Model results for all turtles and mean discharge. " denotes increased strandings rates with increased discharge. # denotes decreased stranding rates with increased discharge. Age class abbreviations: ALL = all turtles, SI = small immature, LI = large immatures, A = adult-sized, ALL IMM = all immature sized animals (small + large), Large = all large turtles (large immatures + adult-sized). The time frame is reported in monthly ranges where responses were noted. The values reported in qAIC are the months with the most significant qAIC value.

Latitude
Age Environmental factors on marine turtle stranding rates When analyzing latitude along the Queensland coastline, there was no evidence that stranding rates were different in different latitude, although there were some noticeable differences (Tables 2-15). When comparing the effect of latitude for discharge, the -25˚block produced split responses whereas the other blocks produced single responses ( Tables 3-5, Tables 10-12). Although stranding rates at different latitudes responded differently the overall pattern of lagged stranding was similar, suggesting the increase in marine turtle stranding was not just a local issue rather, at least, a state-wide issue that occurred and warranted consideration at a state or larger regional level. Given the migratory pattern of marine turtles and their ability to move within the broader range of their individual home sites [39], mitigation needs to consider widespread impacts and not just local habitats of known marine turtle populations.
When analyzing age classes across the variables, there were no observed patterns in relation to which group responded first for each variable (Tables 2-15). This was not expected, it was expected that small immatures would be more susceptible to changes in dietary availability and would show responses before other age classes.
Embayments, when compared to the whole latitudinal block, did not influence the pattern of strandings but did decrease the lag phase for each examined environmental variable (Tables  2-15). This could indicate that embayments are areas of concentrated discharge which is not Table 12. Model results for all turtles peak discharge. " denotes increased strandings rates with increased discharge. # denotes decreased stranding rates with increased discharge. Age class abbreviations: ALL = all turtles, SI = small immature, LI = large immatures, A = adult-sized, ALL IMM = all immature sized animals (small + large), Large = all large turtles (large immatures + adult-sized). The time frame is reported in monthly ranges where responses were noted. The values reported in qAIC are the months with the most significant qAIC value.
As with any exploratory modeling, we identified several limitations that may influence the accuracy of any developed model including distributed sample equality, equal adequate sample sizes for each species, availability of environmental data such as seagrass abundance, habitat type and the distance offshore that an event was recorded. One of the limitations of these models is that the stranding sample size was different for each examined latitudinal block. Larger sample size may make the relationships more noticeable than smaller sample sizes, but as this used one of the longest running and largest datasets available, this may be difficult to correct. The -27˚block recorded the most number of strandings over the study period ( Table 1). This latitudinal block encompasses Moreton Bay which is known to support large fields of seagrass and other aquatic vegetation. The -23˚block recorded the least number of strandings over the study period for a recognized hotspot ( Table 1). Table 13. Model results for all turtles and monthly mean maximum temperature. " denotes increased strandings rates with increased mean maximum temperature. # denotes decreased stranding rates with increased mean maximum temperature. Age class abbreviations: ALL = all turtles, SI = small immature, LI = large immatures, A = adult-sized, ALL IMM = all immature sized animals (small + large), Large = all large turtles (large immatures + adult-sized). The time frame is reported in monthly ranges where responses were noted. The values reported in qAIC are the months with the most significant qAIC value.

Latitude
Age Environmental factors on marine turtle stranding rates The model may have been strengthened by the use of food availability/viability as a factor. However, due to the paucity of data, it was decided to use weather as a proxy to this as weather data is available in immediate time. There is evidence that discharge and rainfall are adequate proxies for seagrass abundance as large-scale seagrass die-off have been closely associated in time and intensity to flooding [40][41][42]. This study may also have been strengthened by determining if different species showed different responses times and directions. This was not possible due to the small sample sizes of the other species occurring within the study location.
Within coastal waters green turtles are almost exclusively herbivorous, feeding principally on seagrass and a wide range of algae and mangrove fruits [43,44]. Occasionally, green turtles feed on macroplankton, including jellyfish, bluebottles, small crustaceans and dead fish [43,44]. Brand-Gardner et al. [45] found that within Moreton Bay small immature green turtles forage selectively on plants with higher nitrogen levels and lower levels of fiber (such as Gracilaria sp.). Due to this strong dependency on aquatic vegetation, it has meant that green turtles living within inshore coast habitats where aquatic vegetation is a large component of their diet have suffered during and post the extreme weather events, such as the flooding in Queensland in 2010-11. Table 14. Model results for all turtles and monthly mean minimum temperature. " denotes increased strandings rates with increased monthly mean minimum temperature. # denotes decreased stranding rates with increased monthly mean minimum temperature. Age class abbreviations: ALL = all turtles, SI = small immature, LI = large immatures, A = adult-sized, ALL IMM = all immature sized animals (small + large), Large = all large turtles (large immatures + adult-sized). The time frame is reported in monthly ranges where responses were noted. The values reported in qAIC are the months with the most significant qAIC value.