Hostname: page-component-78c5997874-t5tsf Total loading time: 0 Render date: 2024-11-15T03:24:40.635Z Has data issue: false hasContentIssue false

Protecting nests of the Critically Endangered South Pacific loggerhead turtle Caretta caretta from goanna Varanus spp. predation

Published online by Cambridge University Press:  29 November 2019

Christine A. Madden Hof*
Affiliation:
World Wide Fund for Nature Australia, Level 1, 17 Burnett Lane, Brisbane, QLD 4000, Australia.
Gabriela Shuster
Affiliation:
World Wide Fund for Nature Australia, Level 1, 17 Burnett Lane, Brisbane, QLD 4000, Australia.
Nev McLachlan
Affiliation:
TurtleCare Volunteers Queensland Inc., Buderim, Australia
Bev McLachlan
Affiliation:
TurtleCare Volunteers Queensland Inc., Buderim, Australia
Saranne Giudice
Affiliation:
Burnett Mary Regional Group for Natural Resource Management, Bundaberg, Australia
Colin Limpus
Affiliation:
Queensland Department of Environment and Science, Brisbane, Australia
Tomoharu Eguchi
Affiliation:
Marine Mammal and Turtle Division, Southwest Fisheries Science Center, National Marine Fisheries Service, National Oceanic and Atmospheric Administration, La Jolla, USA
*
(Corresponding author) E-mail [email protected]

Abstract

The South Pacific subpopulation of the loggerhead turtle Caretta caretta is categorized as Critically Endangered on the IUCN Red List because of significant population declines. Five Queensland beaches support high-density nesting of this subpopulation, but egg and hatchling survival are low at some beaches because of feral and native terrestrial predators. We quantified predation of loggerhead turtle eggs by two species of goanna, Varanus panoptes and Varanus varius, at Wreck Rock beach, one of the turtle's major nesting beaches. In addition, we conducted an experiment to determine the efficacy of a nest protection device. Predation rates at Wreck Rock beach were 15.2% for treatment and 45.8% for non-treatment clutches during the 2013–2014 nesting season. A higher probability of predation (64%) was predicted for the northern beach. Although nests were only partially predated (16.4% of the total number of eggs), nest loss to predators and beach erosion (caused by a cyclone) was 91.7%. If left unmanaged, the cumulative impact of predation and other threats, including those exacerbated by climate change, can cause unsustainable loss of loggerhead turtle nests. This study provides one of the first quantitative data sets on rates of loggerhead turtle clutch predation in the South Pacific. It enhances our understanding of goanna predation impacts and identifies an efficient predator exclusion device for mitigating the effects of terrestrial predators at Wreck Rock beach, and for protecting marine turtle nests across northern Australia and globally.

Type
Article
Copyright
Copyright © 2019 Fauna & Flora International

Introduction

There is one genetic subpopulation for the loggerhead turtle Caretta caretta in the South Pacific Ocean (Bowen et al., Reference Bowen, Kamezaki, Limpus, Hughes, Meylan and Avise1994; Limpus et al., Reference Limpus, Boyle and Sunderland2006; Boyle et al., Reference Boyle, FitzSimmons, Limpus, Kelez, Velez-Zuazo and Weycott2009) and most breeding females from this subpopulation nest on beaches in eastern Australia and New Caledonia (Limpus & Limpus, Reference Limpus, Limpus, Bolten and Witherington2003; Limpus, Reference Limpus2008). Post-hatchlings emerging at these beaches undertake trans-Pacific migrations to the west coast of South America, reaching the coastal waters of Peru and Chile (Kelez et al., Reference Kelez, Mamrique and Velez-Zuazo2005; Boyle et al., Reference Boyle, FitzSimmons, Limpus, Kelez, Velez-Zuazo and Weycott2009; Donoso & Dutton, Reference Donoso and Dutton2010; Alfaro-Shigueto et al., Reference Alfaro-Shigueto, Mangel, Bernedo, Dutton, Seminoff and Godley2011). Along the east coast of Australia, sandy beaches in Queensland support the majority of nesting, with c. 80% occurring at five beaches: Wreck Island, Woongarra Coast, Tryon Island, Erskine Island and Wreck Rock (Limpus, Reference Limpus2008; Limpus & Casale, Reference Limpus and Casale2015).

Despite its protected status and several decades of conservation effort, particularly in Queensland, the South Pacific subpopulation of the loggerhead turtle continues to decline (Limpus et al., Reference Limpus, Parmenter and Chaloupka2013; Limpus & Casale, Reference Limpus and Casale2015) and was categorized as Critically Endangered on the IUCN Red List in 2015 (Limpus & Casale, Reference Limpus and Casale2015). In recent years, although the number of nesting loggerhead turtles has been increasing (Limpus et al., Reference Limpus, Parmenter and Chaloupka2013), the escalated mortality of post-hatchlings from synthetic debris ingestion in the first few months after leaving the nesting beaches (Boyle & Limpus, Reference Boyle and Limpus2008) and bycatch mortality of large post-hatchlings caught in long-line fisheries in Peru and Chile (Donoso & Dutton, Reference Donoso and Dutton2010; Alfaro-Shigueto et al., Reference Alfaro-Shigueto, Mangel, Bernedo, Dutton, Seminoff and Godley2011) have hampered population recovery.

In this context, reducing mortality from other causes at nesting beaches could support population growth. Such causes include the predation of eggs and hatchlings by introduced and native predators, the hatchling loss associated with light pollution, and beach degradation caused by anthropogenic and natural events (Limpus, Reference Limpus2008; Berry et al., Reference Berry, Booth and Limpus2013; Limpus & Casale, Reference Limpus and Casale2015).

Increasing hatchling production by reducing clutch predation is a recommended and common practice for marine turtles (Mroziak et al., Reference Mroziak, Salmon and Rusenko.2000; Engeman et al., Reference Engeman, Martin, Constantin, Noel and Woolard.2003, Reference Engeman, Martin, Smith, Woolard, Crady and Constantin2006; O'Connor et al., Reference O'Connor, Limpus, Hofmeister, Allen and Burnett2017). Native and feral predators include the red fox Vulpes vulpes, raccoon Procyon lotor and pig Sus scrofa. Strategies to reduce depredation of marine turtle eggs have been used with varying levels of success; e.g. predator control through shooting or baiting, deployment of exclusion devices (fencing, cages) and scent deterrents, and clutch relocation (Stancyk et al., Reference Stancyk, Talbert and Dean1980; Addison & Henricy, Reference Addison and Henricy1994; Addison, Reference Addison1997; Ratnaswamy et al., Reference Ratnaswamy, Warren, Kramer and Adam1997; Yerli et al., Reference Yerli, Canbolat, Brown and MacDonald1997; Blamires & Guinea, Reference Blamires and Guinea2003; Baskale & Kaska, Reference Baskale and Kaska2005; Norris et al., Reference Norris, Low, Gordon, Saunders, Lapidge and Lapidge2005; Kurz et al., Reference Kurz, Straley and DeGregorio2011; Engeman et al., Reference Engeman, Martin, Woolard, Stahl, Pelizza, Duffiney and Constantin2012; Lamarre-DeJesus & Griffin, Reference Lamarre-DeJesus and Griffin2013). However, there are few in situ studies examining the effectiveness of these strategies for loggerhead turtles (Schroeder, Reference Schroeder1981; Macdonald et al., Reference Macdonald, Brown, Yerli and Canbolat1994; Yerli et al., Reference Yerli, Canbolat, Brown and MacDonald1997; O'Connor et al., Reference O'Connor, Limpus, Hofmeister, Allen and Burnett2017).

Wreck Rock beach supports the second largest number of loggerhead turtles nesting on the eastern Australian mainland, with c. 400 nests per season (Limpus, Reference Limpus2008). Clutch loss on Wreck Rock beach results from two primary sources: storm surge erosion and flooding of nests, which varies widely between years (C.J. Limpus, 2017, pers. comm.), and depredation by goannas Varanus spp., which may now have the greatest impact on hatchling production at this beach (McLachlan et al., Reference McLachlan, McLachlan, Hof, Giudice, Shuster, Bunce, Eguchi, Whiting and Tucker2015). Historically, predation by feral foxes was a major threat, with predation levels reaching c. 90–95% of clutches laid prior to the mid 1980s (Limpus, Reference Limpus2008). A baiting programme was introduced in 1987 and led to the near-complete elimination of clutch predation by foxes. Although clutch loss to native goannas was considered negligible historically, anecdotal evidence now suggests reduced fox numbers have shifted the balance of predator–prey interactions, resulting in increased goanna predation of turtle nests at Wreck Rock beach. Although goannas have long been known to consume marine turtle eggs (Marquez, Reference Marquez1990, p. 51), their predation rates on loggerhead turtle clutches at Wreck Rock beach remained unquantified (Lei & Booth, Reference Lei and Booth2017a). In addition, there is a dearth of information on varanid ecology, behaviour and predation cues to develop effective goanna management strategies for this region.

Here we report on the outcome of experiments that successfully reduced goanna predation on loggerhead turtle clutches at Wreck Rock beach. We quantify goanna predation rates, provide an effective method to decrease predation and present alternative management interventions to reduce the threat of goanna predation on clutches, which can be applied widely at nesting beaches in northern Australia and globally.

Study area

The study site at Wreck Rock (Fig. 1) comprises 22 km of east-facing beaches from Broadwater Creek to Red Rock. The site lies at the southern end of the Great Barrier Reef Coast Marine Park (Queensland State Government) and Great Barrier Reef Marine Park (Commonwealth), situated within and adjacent to various land tenure (Fig. 1). With a coastline exposed to onshore winds and largely unprotected by the reef, the coastal dunes are a naturally dynamic ecosystem characterized by coastal scrubs, eucalypt woodlands, wet heaths and sedge lands. Delineated by camp and beach access points, Wreck Rock beach is divided into north and south beach sectors and marked every 100 m with timber pegs. Pegs on the north beach sector are numbered 0–65 and on the south beach sector 66–210. The north beach is located adjacent to Deepwater National Park. There are two public access campsites (Middle Rock and Wreck Rock) within the Park (Fig. 1).

Fig. 1 Study site at Wreck Rock beach, land tenure and location of camp sites. Markers pegs 0–65 are on the north beach and 66–210 on the south beach. C and D indicate the numbers of control and treatment (deployment of a predator exclusion device) loggerhead turtle Caretta caretta clutches, respectively, on a stretch of beach between the marker pegs shown. The size of the circles represents the total nesting activity, including failed nesting attempts, per year on the respective beach stretches.

Methods

Data collection

The Wreck Rock Turtle Research Team, under the guidance of the Queensland Government's Queensland Turtle Research Programme, has monitored the Wreck Rock beach study site for loggerhead turtle nesting activity for more than 40 years, since 1977. During the 2013–2014 nesting season (1 December 2013–12 March 2014), in addition to carrying out the standard nesting and hatchling emergence census for all turtle species (loggerhead, green Chelonia mydas and flatback Natator depressus turtles), we monitored predator activity at 57 experimental loggerhead turtle clutches (33 treatment and 24 control) daily. Following oviposition, we marked each nest with flagging tape during the night patrol, and designated them for control or treatment plots (with an exclusion device installed at the latter) the following day. Treatment and control nests were spread throughout the study area to ensure sufficient representation of both on the north (control = 14, treatment = 22) and south beaches (control = 10, treatment = 11; Fig. 1). It was not possible to place control and treatment plots randomly in space and time.

We constructed predator exclusion devices, modelled after those developed and deployed successfully on the Sunshine Coast (O'Connor et al., Reference O'Connor, Limpus, Hofmeister, Allen and Burnett2017), from interlocking aluminium mesh panels. We chose aluminium mesh because it has been used successfully in previous studies and is unlikely to disrupt the hatchlings’ magnetic imprinting (Irwin & Lohmann, Reference Irwin and Lohmann2003). Each exclusion device (cage) consisted of a 1 m2 top panel and four side flaps of 20–25 cm width (Supplementary Plate 1). The mesh size of 70 mm was sufficiently large to allow hatchlings to emerge, but small enough to prevent access by varanid and mammalian predators. We placed predator exclusion devices at a depth of c. 20 cm below the sand surface over 33 treatment nests and left 24 control nests unprotected.

We conducted predator activity surveys of the treatment and control plots during daylight hours at low tide using a pro-forma datasheet. For reasons of practicality, data were recorded within a 152.4 cm (total length of the measuring tape used) radius from each marked nest centre as a standard measure and distance of likely plot predation. We collected the following data: sand disturbance or predator tracks (species identified as goanna, fox, dog, ghost crab); attempted predation (holes had been dug but no empty egg shells were found); predation: holes dug, number of empty egg shells (if greater than half a shell) and number of undeveloped or unhatched eggs on sand surface; and the distance between the nest centre and the egg remains farthest away from it. After recording the data, we manually removed signs of predator activity using footwear, so as to minimize human scent at the nest site.

In addition, we monitored 12 nests (10 treatment and two control nests) for predator activity throughout the incubation period, using infrared wildlife motion-sensor cameras (Reconyx Hyperfire HC600, Reconyx, Holmen, USA). Cameras were secured onto timber marker posts using bungee cords and positioned to aim at the centre of each nest. The cameras recorded ambient air temperatures and still images whenever movement was detected. If a camera-monitored nest was destroyed by predators, we moved the camera to another nest. We identified species that visited the clutches from the recorded camera images. Goannas were not individually identified and we assumed a new visitation event after 60 s of no recording.

We excavated the clutches when the hatchlings emerged, or after 60 days of incubation, to determine hatching success. We calculated hatching success using data from the unpredated clutches and the following equation:

Hatching success = (no. of empty shells)/(no. of unhatched eggs + no. of undeveloped eggs + no. of empty shells). We then used the median hatching success to estimate the total hatchling production of the beach throughout the nesting season.

Statistical analysis

We calculated clutch loss as the proportion of clutches lost to all causes (predation and erosion) out of the total number of control clutches. Predation rates were computed as the proportion of clutches lost to predation out of the total number of control clutches. Nest monitoring was interrupted during 26 January–7 February 2014 because a cyclone made landfall nearby. Therefore, we calculated loss of nests in two distinct periods: before (quantified predation rate) and after (quantified predation + cyclone-caused nest loss) the cyclone event.

To determine the effectiveness of the exclusion devices on loggerhead turtle clutches, we considered five logistic regression models, using three predictor variables: days since the beginning of the experiment (Days), location of the nests measured by the peg numbers (Location), and treatment with a predator exclusion device (Treatment; Table 1). We compared models using the Akaike information criterion (AIC; Akaike, Reference Akaike1974). Using the best model from the analysis on the effectiveness of the treatment, we predicted the probability of predation at all observed nest locations along the entire beach.

Table 1 Logistic regression models used to determine the effectiveness of predator exclusion devices on loggerhead turtle Caretta caretta clutches, showing predictor variables and difference in Akaike information criterion values from best-performing model (ΔAIC).

1 Days, days since the beginning of the experiment; Location, location of the nests measured by the peg numbers along the beach; Treatment, clutch treatment with a predator exclusion device.

To determine whether the exclusion device delayed predation, we compared the mean number of days from laying to predation between control and treatment plots using a Bayesian t test. We hypothesized that the mean number of days between laying and predation would be greater for the treatment plots than for the control plots.

To estimate the total number of hatchlings on the entire beach, we used a parametric bootstrap procedure. For all observed clutches at the nesting beach, we estimated the probability of predation using the logistic regression models. The median nest size and median hatching success rate were used to compute the total number of hatchlings for the unpredated nests. We repeated this process 7,000 times to obtain the uncertainty in the estimate from the predation probability.

We conducted all statistical analyses in R 3.1.2 (R Core Team, 2015) with packages ggplot2 (Wickham, Reference Wickham2009) and gamm4 (Wood & Scheipl, Reference Wood and Scheipl2014). Bayesian computations were conducted using jags 3.3.0 (Plummer, Reference Plummer2015) through rjags (Plummer, Reference Plummer2015) in R.

Results

During the 2013–2014 nesting season (1 December 2013–12 March 2014) we recorded 78 loggerhead, 30 green and 20 flatback turtles laying one or more clutches at Wreck Rock beach (Supplementary Table 1). Of the loggerhead turtles, 16 were new recruits and 62 were recaptures from previous nesting seasons. During the nesting season, we recorded 218 loggerhead nesting activities, with 195 successful nests and 23 failed attempts. Of the 195 successful nesting events, 148 clutches were laid on the north and 47 on the south beach. We used 57 of the 195 successfully laid nests in our study, deployed exclusion devices on 33 treatment clutches and used 24 as control. During the study, 19 plots (eight treatment and 11 controls) were lost to erosion when cyclone Dylan made landfall near Wreck Rock beach on 31 January 2014.

Clutch loss

During the 2013–2014 season, 11 of 24 control clutches were depredated (mean = 45.8%, 95% CI = 26.2–66.8%), an additional 11 were lost to erosion. In total, clutch loss to predators and erosion events (exacerbated by a cyclone) was 91.7% (95% CI = 71.5–98.5%). The predation rate was greater prior to the cyclone (41.67%, 10 of 24 control plots depredated; 95% CI = 22.8–63.1%) than after the cyclone (7.7%, 1 of 13 control plots depredated; 95% CI = 0.4–37.9%). The majority of predation of control nests occurred on the northern beach (9/14 = 64%; 95% CI = 35.6–86.0%).

Exclusion devices appeared to be effective in reducing clutch predation. With exclusion devices deployed at 33 nests, only one hatchling became entrapped during emergence. Compared with the control plots (11 of 24 control nests; 45.8%; 95% CI = 26.2–66.8%), fewer treatment plots were depredated during the nesting season (5 of 33 treatment nests; 15.2%; 95% CI = 5.7–32.7%). The point estimate of predation with exclusion devices was 6% (2/33; 95% CI = 0.7–20.2%) prior to the cyclone, and 12% (3/25; 95% CI = 3.2–32.3%) after the cyclone. Although the point estimate was greater for the treatment plots than control (12.0% vs 7.7%), the binomial 95% CI overlapped (3.2–32.3% vs 0.3–37.9%). Because of the large uncertainties associated with small sample sizes, they were not statistically different.

Amongst the models used to determine the effects of the exclusion device on the predation of loggerhead turtle clutches, Model 4 had the smallest AIC value, although differences in AIC values between the top 3 models were < 3, indicating that these models were not substantially different (Table 1). In the following analyses, we use the best model, which included treatment, location, and their interactions.

The interaction term was significant at α = 0.10 (Table 2), indicating that the effectiveness of the treatment varied significantly along the beach. This can be seen in Fig. 2, where the predicted probability of predation decreases with increasing peg number (i.e. from north to south along the beach) for the control plots, but increases for the treatment clutches. However, because of the small sample sizes in the area farthest from the field station, on the south beach, confidence intervals are wide in this area. For the area with larger sample sizes, on the north beach, the predicted probability of predation was smaller for the treatment than for the control plots.

Table 2 Estimated linear model coefficients of the best-performing model (Model 4), with standard deviations and P values. The model was Predation ~ Intercept + Treatment + Location + Treatment × Location.

Fig. 2 Predicted probability of loggerhead turtle nest predation as a function of location on the beach and control vs predator exclusion device treatment. Shaded areas correspond to 95% confidence intervals.

The main effect of the exclusion device was also significant at α = 0.10 (Table 2), indicating that the devices significantly reduced the predation of loggerhead turtle clutches. This is supported by the percentage of predated clutches in control vs treatment plots (47.8 vs 15.2%, respectively).

Hatching success

The mean clutch size was 113.6 ± SE 7.7 (n = 31; median 118) and the mean hatching success was 0.779 ± SE 0.056 (n = 42; median 0.894). Using the most parsimonious logistic regression model (Model 4) from the previous analysis, we predicted the probability of predation at all observed nest locations; i.e. extrapolated to the entire beach (Fig. 3). Using the probability of predation, median hatching success, median clutch size, and a parametric bootstrap analysis, we estimated the total number of hatchlings produced at Wreck Rock beach during the nesting season to be 11,599 with a 95% CI of 10,439–12,865 (Supplementary Fig. 1). Some loggerhead turtle clutches were only partially predated. Although sample sizes were small, the mean number of predated eggs per nest was c. 16 ± SE 4.5 (16.4%, n = 18). Undeveloped eggs (unhatched eggs with no obvious embryo) also accounted for 16% of the total eggs in nests.

Fig. 3 Predicted predation probabilities with respect to beach marker pegs (location) and time since the beginning of the experiment (days). The size of the points corresponds to the probabilities.

The Bayesian equivalent of a two-sampled t test indicated there was no significant difference in the timing of predation between treatment and control nests (Supplementary Fig. 2). Further, the 95% posterior interval (Gelman et al., Reference Gelman, Carlin, Stern, Dunson, Vehtari and Rubin2014) of the contrast coefficient (control vs treatment effects) included zero (−1.973, 0.724). For both control and treatments, however, predation seemed to occur either in the early or late stages of incubation (Supplementary Fig. 2).

Predation

With individual cameras operating for 51–62 days, the total camera-trapping time was 13,800 hours over a period of 79 days. Species recorded by automatic cameras included hares Lepus sp., kites Milvus sp., cats Felis sp., foxes Vulpes sp., emus Dromaius sp., pied butcherbirds Cracticus nigrogularis and wallabies Macropus sp. The two species of goannas recorded near turtle nests were the yellow spotted or Argus monitor Varanus panoptes and the lace monitor Varanus varius (A. Amey, 2014, pers. comm.). A third goanna Varanus gouldii may be present, but the species identification could not be confirmed. As expected, goanna activity tended to occur during daylight hours (6.00–18.00).

Based on tracks, predation at both treatment and control plots was carried out by goannas, with only one incidence of fox predation recorded, at a treatment plot. On one occasion, more than six goannas were observed predating on loggerhead turtle hatchlings as they emerged from a nest. This was the first recorded observation of goanna predation on emerging hatchlings at Wreck Rock beach. On several occasions we observed predation on nests adjacent to control plots (i.e. on nests not included in this study).

We recorded 177 goanna visitation events to 12 loggerhead turtle nests monitored by automated cameras. The yellow spotted monitor was found most frequently (53 recorded visits) followed by the lace monitor (37 visits). Of the 12 monitored clutches, 8 were visited multiple times by the same species on 74 occasions (46 times by yellow spotted monitors and 28 times by lace monitors). The remaining four nests were visited once by one goanna species, and at two of these nests the second species continued to visit multiple times after the first species had visited.

Ambient air temperature during goanna visits to the clutches was 20–47 °C. We found differences between the mean temperatures at visits by the two goanna species and other species (birds and mammals; Supplementary Fig. 4) Lace monitors visited turtle nests at a lower ambient temperature (mean = 32.7 ± SE 0.44 °C, n = 46) than yellow spotted monitors (36.8 ± SE 0.55 °C, n = 53).

Discussion

As the South Pacific loggerhead turtle subpopulation continues to decline, it requires renewed attention to decrease mortality at all life history stages. Nest predation by terrestrial predators may significantly reduce recruitment and could lead to additional declines in the already depleted population (CMS, 2014). This study provides new information about goanna depredation and demonstrates that exclusion devices effectively reduce predation. We also provide a recent estimate of the predation rate on loggerhead turtle clutches at a major nesting beach in the South Pacific.

The predation rate of loggerhead turtle clutches (45.8%) at Wreck Rock beach was lower than that reported at other nesting beaches. For example, 52% of flatback turtle clutches were predated at Fog Bay (Blamires, Reference Blamires1999; Blamires & Guinea, Reference Blamires and Guinea2003), 65–70% of flatback turtle clutches at Pennefather beach (J. Doherty and Cape York Peninsula Development Association, unpubl. data in Whytlaw et al., Reference Whytlaw, Edwards and Congdon2013), 64–89% of green and loggerhead turtle clutches in Turkey (Macdonald et al., Reference Macdonald, Brown, Yerli and Canbolat1994), and 67–95% of loggerhead turtle clutches in North America (Stancyk et al., Reference Stancyk, Talbert and Dean1980; Engeman et al., Reference Engeman, Martin, Constantin, Noel and Woolard.2003). During the 2014–2015 and 2015–2016 nesting seasons at Wreck Rock beach, Lei & Booth (Reference Lei and Booth2017a) reported 57.7% and 17.4% of clutches predated by goannas, respectively. However, their study was confined to a 6-km stretch of beach, whereas we monitored the entire 22-km stretch. The difference in results may be a result of this spatial difference in sampling effort and the heterogeneity of predation across a rookery (Blamires, Reference Blamires1999). It is also possible that predation is temporally variable because of fluctuations in predator and prey density.

In recent years, as also evidenced in Lei & Booth (Reference Lei and Booth2017a), goannas have been the primary predators of turtle eggs and hatchlings at Wreck Rock beach. For the first time, we observed predation of hatchlings by multiple goannas. Foxes are now minor predators of loggerhead turtle clutches, presumably the result of the long-term fox baiting programme. Yellow spotted and lace monitors were the most prominent visitors to turtle nests at the study site. A third species, Varanus gouldii, is frequently found in close proximity to yellow spotted monitors (Shine, Reference Shine1986; A. Amey, 2014, pers. comm.) but further studies are needed to confirm its presence.

In this study, the aluminium cage predator exclusion devices were effective in reducing predation of loggerhead turtle nests at Wreck Rock beach. With only one entrapped hatchling amongst 33 cage deployments, the devices were deemed successful for reducing predation and letting hatchlings pass through. The same conclusion was also reached by Lei & Booth (Reference Lei and Booth2017b) in subsequent years. Given Wreck Rock beach is predicted to produce c. 12,000 hatchlings per season, the use of these anti-predator devices should result in increased hatchling production.

The loss of clutches observed in this study (from predation and other causes combined) was 91.7%, which is a cause for concern. Historically only quantified for beaches at Mon Repos (mean 13% loss of entire clutches from natural causes, range 8.4–83.0%; Limpus, Reference Limpus2008), clutches are regularly lost to natural erosion and flooding at Wreck Rock beach (N. McLachlan, unpubl. data). These events will continue to be exacerbated with severe cyclones, floods and storms predicted to increase in frequency with climate change (IPCC, 2007). During the study, a natural tide storm surge and wind additionally exposed turtle nests to predation and increased the total number of clutches lost. This loss was greater for the treatment nests, probably because the exclusion devices provided visual cues for the predators. In addition, hatchlings are affected by light pollution associated with coastal development. Overall, Wreck Rock beach grossly exceeds a sustainable level of annual clutch loss of c. 30% (CMS, 2014). Consequently, it is important to develop ongoing predator management (cognizant of the threat posed by extreme weather events) as this may be the only threat that can be addressed effectively in a timeframe that allows the population to recover.

The deployment of exclusion devices, although effective, may not be feasible for nesting beaches with moderate to high turtle numbers (e.g. > 20 turtles per night). Devices may impede nesting attempts (Kurz et al., Reference Kurz, Straley and DeGregorio2011), and resources are often limited (Lei & Booth, Reference Lei and Booth2017b; Lei et al., Reference Lei, Booth and Dwyer2017). A suite of alternative predator control methods (some of which are already being applied at Wreck Rock beach) should be included in future studies to develop a conservation strategy considerate of the main predator species, habitat and climatic conditions, and budgetary and logistical constraints. Additional types of predator exclusion devices that could be deployed include plastic mesh nest protectors (Lei & Booth, Reference Lei and Booth2017b) and flat chain link screening or less rigid wire mesh cages, if not disruptive to hatchlings magnetic imprinting (Addison & Henricy, Reference Addison and Henricy1994; Addison, Reference Addison1997). The application of scent deterrents such as habanero powder (see Ratnaswamy et al., Reference Ratnaswamy, Warren, Kramer and Adam1997; Lamarre-DeJesus & Griffin, Reference Lamarre-DeJesus and Griffin2013; Lei et al., Reference Lei, Booth and Dwyer2017), human scent (see Burke et al., Reference Burke, Schneider and Dolinger2005), visual disturbance (such as flags, see Burke et al., Reference Burke, Schneider and Dolinger2005), or the use of pheromones (goanna's own or other species such as cane toads) may be effective in deterring goannas. Nest relocation (including into hatcheries) should only be trialled as a last resort because it could reduce hatchling imprinting and success (see Kornaraki et al., Reference Kornaraki, Matossian, Mazaris, Matsinos and Margaritoulis2006).

Culling of goannas to protect turtle nests is not viable as they are protected in Australia. Coastal populations of the yellow spotted monitor may be the species’ last stronghold, as many inland populations have declined as a result of a toxic diet of cane toads (Ujvari & Madsen, Reference Ujvari and Madsen2009; Shine, Reference Shine2010). Because the species is not categorized as threatened on the IUCN Red List, Lei & Booth (Reference Lei and Booth2017b; Lei et al., Reference Lei, Booth and Dwyer2017) suggest the temporary removal of male yellow spotted monitors (the primary predators of turtle eggs) during the turtle nesting season. However, given the monitor's ecological role as a mesopredator, resource managers must first understand the local predator–prey interactions and the ecosystem-level effects of any predator control method selected (Prugh et al., Reference Prugh, Stoner, Epps, Bean, Ripple, Laliberte and Brashares2009; Welicky et al., Reference Welicky, Wyneken and Noonburg2012). Studies on the raccoon Procyon lotor, another mesopredator that targets marine turtle clutches, showed its removal had no effect on mainland clutch depredation (Ratnaswamy et al., Reference Ratnaswamy, Warren, Kramer and Adam1997; Barton & Roth, Reference Barton and Roth2007). Tsellarius et al. (Reference Tsellarius, Men'shikov and Tsellarius2011) suggested goannas scent-mark territorial boundaries to keep strangers out and control local population sizes, so the removal of individual goannas from a local populations could potentially lead to an increase in goanna numbers.

Goanna predation and nest selectivity is probably driven by olfactory and visual cues and influenced by spatial and temporal nest deposition (Blamires & Guinea, Reference Blamires and Guinea2003; Blamires, Reference Blamires2004; Welicky et al., Reference Welicky, Wyneken and Noonburg2012), proximity to urbanized areas (Smith & Engeman, Reference Smith and Engeman2002; Blamires & Guinea, Reference Blamires and Guinea2003; Prange et al., Reference Prange, Gehrt and Wiggers2004) and human and goanna conspecific activity (Ferreira, Reference Ferreira2012). In contrast to Welicky et al. (Reference Welicky, Wyneken and Noonburg2012), predation risk in this study was more probable on the north beach in areas heavily utilized by humans, with campsites and food waste. This raises the question of why goannas have become more abundant at Wreck Rock beach (they were considered uncommon in the 1970s; C. Limpus, 2017, pers. comm.) and what drives their behaviour. It is possible that the control of apex predators (dingoes and foxes) has altered predator–prey relationships and reduced competition for goannas. This, together with increased human activity and camp site waste (particularly on the north beach), could have increased food availability for goannas, supporting them in greater numbers (Prugh et al., Reference Prugh, Stoner, Epps, Bean, Ripple, Laliberte and Brashares2009).

Research with a focus on goanna ecology and loggerhead turtle clutch predation could improve management interventions for Wreck Rock beach. Such studies should include research into the biology and behaviour of goannas, examination of long-term predation effects and human impacts on predator behaviour. An examination of clutch predation could uncover vital information regarding the timing of predation (Ferreira, Reference Ferreira2012; Welicky et al., Reference Welicky, Wyneken and Noonburg2012), repeated visitation by predators, and complete (Limpus, Reference Limpus1971) vs partial clutch loss (Chatto & Baker, Reference Chatto and Baker2008). For example, further research could provide insights into why control clutches were predated less than anticipated, and why nests not included in this study but directly adjacent to control plots were also predated (as per the simulation in Blamires & Guinea, Reference Blamires and Guinea2003). Although the lace monitor was not previously considered a predator of turtle eggs, this species accounted for more than one-third of total predation and visited turtle nests at a lower temperature than the yellow spotted monitor (33 vs 37 °C, respectively).

Until further studies are undertaken, interim management is recommended particularly for the north beach. Camp site usage and waste management should be reviewed as a priority to reduce goanna activity in this area. This could be achieved with relatively little effort through education and enforcement. Future site-specific management is necessary to maximize hatchling success. The exclusion devices deployed in this study are effective in reducing depredation on marine turtle nesting beaches, but are not cost-effective for Wreck Rock beach because of high nesting numbers and limited resources. An alternative predator control programme should be explored, to provide the most efficient allocation of resources and management. This will be the most robust strategy to maximize hatchling production and thus contribute to future recruitment of the declining and Critically Endangered South Pacific subpopulation of the loggerhead turtle.

Acknowledgements

We thank the volunteers of TurtleCare's Wreck Rock Turtle Research Team and the Gidarjil Aboriginal Corporation for their assistance with field work; and Kirsten Wortel and the Burnett Mary Regional Group (BMRG) for assistance with mapping. We acknowledge funding from WWF-Australia and BMRG.

Author contributions

Study design: all authors; data collection: CAMH, GS, NM, BM, SG; data analysis: CAMH, GS, TE; writing: CAMH; revisions: all authors.

Conflicts of interest

None.

Ethical standards

This research abided by the Oryx guidelines on ethical standards. Research protocols for this study were approved by an authorised ethics committee (SA 2015/11/531) and authority under the Nature Conservation Act 1994.

Footnotes

Supplementary material for this article is available at https://doi.org/10.1017/S0030605318001564

References

Addison, D.S. (1997) Galvanized wire cages can protect nest depredation. Marine Turtle Newsletter, 76, 811.Google Scholar
Addison, D.S. & Henricy, S. (1994) A comparison of galvanized wire mesh cages vs flat chain-link screen in preventing Procyon lotor depredation of Caretta caretta nests. NOAA Technical Memorandum, 351, 174.Google Scholar
Akaike, H. (1974) A new look at the statistical model identification. IEEE Transactions on Automatic Control, 19, 716723.CrossRefGoogle Scholar
Alfaro-Shigueto, J., Mangel, J., Bernedo, F., Dutton, P.H., Seminoff, J.A. & Godley, B.J. (2011) Small-scale fisheries of Peru: a major sink for marine turtles in the Pacific. Journal of Applied Ecology, 48, 14321440.CrossRefGoogle Scholar
Barton, B.T. & Roth, J. (2007) Raccoon removal on sea turtle nesting beaches. The Journal of Wildlife Management, 71, 12341237.CrossRefGoogle Scholar
Baskale, E. & Kaska, Y. (2005) Sea turtle nest conservation techniques on southwestern beaches in Turkey. Israel Journal of Zoology 51, 1326.CrossRefGoogle Scholar
Berry, M., Booth, D.C. & Limpus, C.J. (2013) Artificial lighting and disrupted sea-finding behaviour in hatchling loggerhead turtles (Caretta caretta) on the Woongarra coast, south-east Queensland, Australia. Australian Journal of Zoology, 61, 137145.CrossRefGoogle Scholar
Blamires, S.J. (1999) Quantifying the impact of predation on sea turtle nests by varanids at Fog Bay. MSc thesis. Northern Territory University, Darwin, Australia.Google Scholar
Blamires, S.J. (2004) Habitat preferences of coastal goannas (Varanus panoptes): are they exploiters of sea turtle nests at Fog Bay, Australia? Copeia, 2, 370377.CrossRefGoogle Scholar
Blamires, S.J. & Guinea, M.L. (2003) Emergence success of flatback Sea turtles (Natator depressus) at Fog Bay, Northern Territory, Australia. Chelonian Conservation and Biology, 4, 548556.Google Scholar
Bowen, B.W., Kamezaki, N., Limpus, C.J., Hughes, G.R., Meylan, A.B. & Avise, J.C. (1994) Global phylogeography of the loggerhead turtle (Caretta caretta) as indicated by mitochondrial DNA haplotypes. Evolution, 48, 18201828.Google ScholarPubMed
Boyle, M.C. & Limpus, C.J. (2008) The stomach contents of post-hatchling green and loggerhead sea turtles in the southwest Pacific: an insight into habitat association. Marine Biology, 155, 233241.CrossRefGoogle Scholar
Boyle, M.C., FitzSimmons, N.N., Limpus, C.J., Kelez, S., Velez-Zuazo, X. & Weycott, M. (2009) Evidence for transoceanic migrations by loggerhead sea turtles in the southern Pacific Ocean. Proceedings of the Royal Society B: Biological Sciences, 276, 19931999.CrossRefGoogle ScholarPubMed
Burke, R.L., Schneider, C.M., Dolinger, M.T. (2005) Cues used by raccoons to find turtle nests: effects of flags, human scent, and diamond-backed terrapin sign. Journal of Herpetology, 39, 312315.CrossRefGoogle Scholar
Chatto, R. & Baker, B. (2008) The Distribution and Status of Marine Turtle Nesting in the Northern Territory. Parks and Wildlife Service of the Northern Territory, Palmerston, Australia.Google Scholar
CMS (Convention on the Conservation of Migratory Species of Wild Animals) (2014) Single Species Action Plan for the Loggerhead Turtle, Caretta caretta, in the South Pacific Ocean. https://www.cms.int/sites/default/files/document/COP11_Doc_23_2_2_Rev1_Annex_2_SSAP_Loggerhead_Turtle_E_0.pdf [accessed 17 May 2019].Google Scholar
Donoso, M. & Dutton, P. (2010) Sea turtle bycatch in the Chilean pelagic longline fishery in the south-eastern Pacific: opportunities for conservation. Biological Conservation, 143, 26722684.CrossRefGoogle Scholar
Engeman, R.M., Martin, R.E., Constantin, B., Noel, R. & Woolard., J. (2003) Monitoring predators to optimize their management for marine turtle nest protection. Biological Conservation, 113, 171178.CrossRefGoogle Scholar
Engeman, R.M., Martin, R.E., Smith, H.T., Woolard, J., Crady, C.K., Constantin, B. et al. (2006) Impact on predation of sea turtle nests when predator control was removed midway through the nesting season. Wildlife Research, 33, 187192.CrossRefGoogle Scholar
Engeman, R., Martin, R.E., Woolard, J., Stahl, M., Pelizza, C., Duffiney, A. & Constantin, B. (2012) An ideal combination for marine turtle conservation: exceptional nesting season, with low nest predation resulting from effective low-cost predator management. Oryx, 46, 229235.CrossRefGoogle Scholar
Ferreira, M.B.M.d.S. (2012) Nesting habitat preferences and nest predation of green turtles (Chelonia mydas) in the Bijagós Archipelago, Guinea Bissau. MSc thesis, Universidade de Lisboa, Lisbon, Portugal.Google Scholar
Gelman, A., Carlin, J.B., Stern, H.S., Dunson, D.B., Vehtari, A. & Rubin, D.B. (2014) Bayesian Data Analysis. CRC Press, Boca Raton, USA.CrossRefGoogle Scholar
Intergovernmental Panel on Climate Change (2007) Climate Change 2007–The Physical Science Basis: Working Group I Contribution to the Fourth Assessment Report of the IPCC. Cambridge University Press, Cambridge, UK.Google Scholar
Irwin, W.P. & Lohmann, K.J. (2003) Magnet-induced disorientation in hatchling loggerhead sea turtles. Journal of Experimental Biology, 206, 497501.CrossRefGoogle ScholarPubMed
Kelez, S., Mamrique, C. & Velez-Zuazo, X. (2005) Conservation of Sea Turtles Along the Coast of Peru. Unpublished Report to UNEP/CMS from Asociación Peruana para la Conservación de la Naturaleza y Grupo de Tortugas Marinas—Perú. https://www.cms.int/sites/default/files/publication/Tortugas_Peru_informe_final_CMS_En.pdf [accessed 3 July 2019].Google Scholar
Kornaraki, E., Matossian, D.A., Mazaris, A.D., Matsinos, Y.G. & Margaritoulis, D. (2006) Effectiveness of different conservation measures for loggerhead sea turtle (Caretta caretta) nests at Zakynthos Island, Greece. Biological Conservation, 130, 324330.CrossRefGoogle Scholar
Kurz, D.J., Straley, K.M. & DeGregorio, B.A. (2011) Out-foxing the red fox: how to best protect the nests of the Endangered loggerhead marine turtle Caretta caretta from mammalian predation? Oryx, 46, 223228.CrossRefGoogle Scholar
Lamarre-DeJesus, A.S. & Griffin, C.R. (2013) Use of habanero pepper powder to reduce depredation of loggerhead sea turtle nests. Chelonian Conservation and Biology, 12, 262267.CrossRefGoogle Scholar
Lei, J. & Booth, D.T. (2017a) Who are the important predators of sea turtle nests at Wreck Rock beach? Peer J, 5, e3515.CrossRefGoogle Scholar
Lei, J. & Booth, D.T. (2017b) How best to protect the nests of the Endangered loggerhead turtle Caretta caretta from monitor lizard predation? Chelonian Conservation and Biology, 16, 246249.CrossRefGoogle Scholar
Lei, J., Booth, D.T. & Dwyer, R.G. (2017) Spatial ecology of yellow-spotted goannas adjacent to a sea turtle nesting beach. Australian Journal of Zoology, 65, 7786.CrossRefGoogle Scholar
Limpus, C.J. (1971) The flatback turtle, Chelonia depressa Garman in Southeast Queensland, Australia. Herpetologica, 27, 431446.Google Scholar
Limpus, C.J. (2008) A Biological Review of Australian Marine Turtles. 1. Loggerhead Turtle Caretta caretta (Linneaus). Queensland Environment Protection Agency, Brisbane, Australia.Google Scholar
Limpus, C. & Casale, P. (2015) Caretta caretta (South Pacific subpopulation). In The IUCN Red List of Threatened Species 2015: e.T84156809A84156890. http://dx.doi.org/10.2305/IUCN.UK.2015-4.RLTS.T84156809A84156890.en [accessed 8 January 2017].CrossRefGoogle Scholar
Limpus, C.J. & Limpus, D.J. (2003) Loggerhead turtles in the equatorial and Southern Pacific Ocean: a species in decline. In Loggerhead Sea Turtles (eds Bolten, A.B. & Witherington, B.E.), pp. 199209. Smithsonian Institution, Washington, DC, USA.Google Scholar
Limpus, C.J., Boyle, M. & Sunderland, T. (2006) New Caledonian loggerhead turtle population assessment: 2005 pilot study. In Kinan, I. Proceedings of Second Western Pacific Sea Turtle Cooperative Research and Management Workshop. Volume II. North Pacific Loggerhead Sea Turtles, pp. 7792. Western Pacific Regional Fisheries Management Council, Honolulu, Hawaii, USA.Google Scholar
Limpus, C.J., Parmenter, C.J. & Chaloupka, M. (2013) Monitoring of Coastal Sea Turtles: Gap Analysis 1. Loggerhead Turtles, Caretta caretta, in the Port Curtis and Port Alma Region. Report Produced for the Ecosystem Research and Monitoring Program Advisory Panel as Part of Gladstone Ports Corporation's Ecosystem Research and Monitoring Program. Gladstone Ports Corporation, Gladstone, Australia.Google Scholar
Macdonald, D.W., Brown, L., Yerli, S. & Canbolat, A-F. (1994) Behavior of red foxes, Vulpes vulpes, caching eggs of loggerhead turtles, Caretta caretta. Journal of Mammalogy, 75, 985988.CrossRefGoogle Scholar
Marquez, M.R. (1990) Sea Turtles of the World. Vol. 11. An Annotated and Illustrated Catalogue of Sea Turtle Species Known to Date. Food and Agriculture Organisations of the United Nations, Rome, Italy.Google Scholar
McLachlan, N., McLachlan, B., Hof, C., Giudice, S., Shuster, G., Bunce, A. & Eguchi, T. (2015) Predator reduction strategies for protecting loggerhead turtle nests at Wreck Rock beach in Queensland. In: Proceedings of the Second Australian and Second Western Australian Marine Turtle Symposia, Perth 25-27 August 2014 (compiled by Whiting, S.D. & Tucker, A.), p. 57. Science Division, Department of Parks and Wildlife, Perth, Australia.Google Scholar
Mroziak, M.L., Salmon, M. & Rusenko., K. (2000) Do wire cages protect sea turtles from foot traffic and mammalian predators? Chelonian Conservation Biology, 3, 693698.Google Scholar
Norris, A., Low, T., Gordon, I., Saunders, G., Lapidge, S., Lapidge, K. et al. (2005) Review of the Management of Feral Animals and their Impact on Biodiversity in the Rangelands: a Resource to Aid NRM Planning. Pest Animal Control Cooperative Research Centre, Canberra, Australia.Google Scholar
O'Connor, J.M., Limpus, C.J., Hofmeister, K.M., Allen, B.L., Burnett, S.E. (2017) Anti-predator meshing may provide greater protection for sea turtle nests than predator removal. PLOS ONE, 12, e0171831.CrossRefGoogle ScholarPubMed
Plummer, M. (2015) rjags: Bayesian Graphical Models using MCMC. R package version 3-15. http://CRAN.R-project.org/package=rjags [accessed 17 May 2019].Google Scholar
Prange, S., Gehrt, S.D. & Wiggers, E.P. (2004) Influences on anthropogenic resources on raccoon (Procyon lotor) movements and spatial distribution. Journal of Mammology, 85, 483490.CrossRefGoogle Scholar
Prugh, L.R., Stoner, C.J., Epps, C.W., Bean, W.T., Ripple, W.J., Laliberte, A.S. & Brashares, J.S. (2009) The rise of the mesopredator. Bioscience, 59, 779791.CrossRefGoogle Scholar
Ratnaswamy, M.J., Warren, R.J., Kramer, M.T. & Adam, M.D. (1997) Comparisons of lethal and nonlethal techniques to reduce raccoon depredation of sea turtle nests. The Journal of Wildlife Management, 61, 368376.CrossRefGoogle Scholar
R Core Team (2015) R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.Google Scholar
Schroeder, B.A. (1981) Predation and nest success in two species of marine turtles (Caretta caretta and Chelonia mydas) at Merritt Island, Florida. Florida Scientist, 44, 35.Google Scholar
Shine, R. (1986) Food habits, habitats and reproductive biology of four sympatric species of varanid lizards in tropical Australia. Herpetologica, 42, 346360.Google Scholar
Shine, R. (2010) The ecological impact of invasive cane toads (Bufo marinus) in Australia. The Quarterly Review of Biology, 85, 253291.CrossRefGoogle ScholarPubMed
Smith, H.T. & Engeman, R.M. (2002) An extraordinary raccoon density at an urban park in Florida. Canadian Field Naturalist, 116, 636639.Google Scholar
Stancyk, S.E., Talbert, O.R. & Dean, J.M. (1980) Nesting activity of the loggerhead turtle Caretta caretta in South Carolina II. Protection of nests from raccoon predation by transplantation. Biological Conservation, 18, 289298.CrossRefGoogle Scholar
Tsellarius, A.Y., Men'shikov, Y.G. & Tsellarius, E.Y. (2011) Spacing pattern and reproduction in Varanus griseus of Western Kyzylkum. Russian Journal of Herpetology, 2, 153165.Google Scholar
Ujvari, B. & Madsen, T. (2009) Increased mortality of naive varanid lizards after the invasion of non-native cane toads (Bufo marinus). Herpetological Conservation and Biology, 4, 248251.Google Scholar
Welicky, R.L., Wyneken, J. & Noonburg, E.G. (2012) A retrospective analysis of sea turtle nest depredation patterns. The Journal of Wildlife Management, 76, 278284.CrossRefGoogle Scholar
Whytlaw, P.A., Edwards, W. & Congdon, B.C. (2013) Marine turtle nest depredation by feral pigs (Sus scrofa) on the Western Cape York Peninsula, Australia: implications for management. Wildlife Research, 40, 377384.CrossRefGoogle Scholar
Wickham, H. (2009) Ggplot2: Elegant Graphics for Data Analysis. Springer, New York, USA.CrossRefGoogle Scholar
Wood, S. & Scheipl, F. (2014) gamm4: Generalized additive mixed models using mgcv and lme4. R package version 0.2-3. http://CRAN.R-project.org/package=gamm4 [accessed 17 May 2019].Google Scholar
Yerli, S., Canbolat, A.F., Brown, L.J. & MacDonald, D.W. (1997) Mesh grids protect loggerhead turtle (Caretta caretta) nests from red fox (Vulpes vulpes) predation. Biological Conservation, 82, 109111.CrossRefGoogle Scholar
Figure 0

Fig. 1 Study site at Wreck Rock beach, land tenure and location of camp sites. Markers pegs 0–65 are on the north beach and 66–210 on the south beach. C and D indicate the numbers of control and treatment (deployment of a predator exclusion device) loggerhead turtle Caretta caretta clutches, respectively, on a stretch of beach between the marker pegs shown. The size of the circles represents the total nesting activity, including failed nesting attempts, per year on the respective beach stretches.

Figure 1

Table 1 Logistic regression models used to determine the effectiveness of predator exclusion devices on loggerhead turtle Caretta caretta clutches, showing predictor variables and difference in Akaike information criterion values from best-performing model (ΔAIC).

Figure 2

Table 2 Estimated linear model coefficients of the best-performing model (Model 4), with standard deviations and P values. The model was Predation ~ Intercept + Treatment + Location + Treatment × Location.

Figure 3

Fig. 2 Predicted probability of loggerhead turtle nest predation as a function of location on the beach and control vs predator exclusion device treatment. Shaded areas correspond to 95% confidence intervals.

Figure 4

Fig. 3 Predicted predation probabilities with respect to beach marker pegs (location) and time since the beginning of the experiment (days). The size of the points corresponds to the probabilities.

Supplementary material: PDF

Madden Hof et al. supplementary material

Madden Hof et al. supplementary material

Download Madden Hof et al. supplementary material(PDF)
PDF 283.1 KB