Lost to the Sea: Predicted Climate Change Threats to Saltwater Crocodile Nesting Habitat

Abstract

Climate change is predicted to have devastating impacts on apex predators such as eliminating their required habitats. Crocodilians are no exception as most species require freshwater for nesting, and such freshwater habitats are particularly vulnerable to saltwater inundation (SWI) caused by the sea level rise (SLR) from global warming. Here, we examined the impacts of climate change on saltwater crocodiles Crocodylus porosus in terms of the potential loss of nesting habitat in the Northern Territory, Australia; an area that contains the world’s most extensive nesting habitat for the species. Our spatial model, derived from 730 nest locations and selected environmental features, estimated a total of 32,306.91 km² of current suitable habitat across the study region. The most important variable was distance to perennial lakes (71.0% contribution, 87.5% permutation importance), which is negatively correlated with nesting habitat suitability. We found that projected changes in temperature and rainfall by 2100 could impact the area of suitable nesting habitat negatively or positively (0.33% decrease under low future emission climate scenario, and 32.30% increase under high emission scenario). Nevertheless, this can be canceled by the strong negative impact of SLR and concomitant SWI on nesting areas. A portion (16.40%) of the modeled suitable habitat for a subsection of our study area, the Kakadu Region, were already subject to > 0.25 m SWI in 2013. The suitable area for nesting in this region is predicted to be further reduced to 1775.70 km² with 1.1 m SLR predicted for 2100, representing 49.81% loss between 2013 and 2100. Although the estimates of habitat loss do not account for the potential creation of new habitat, nor for the uncertainty in the degree of future SLR, our results suggest that SLR driven by continuing global warming can be the major threat to mound-nest-building crocodilians including C. porosus, rather than direct impacts from changes in temperature and rainfall. The degree of impact on saltwater crocodiles will be determined by the interplay between the loss of nesting habitat, which would appear inevitable under current global warming, and the ability to expand into new areas created by the expansion of the tropics.

Introduction

Apex predators include some of the world’s most imperiled species and climate change is predicted to have devastating impacts on some members of this important trophic group. Perhaps the best publicized example of climate impacts on an apex predator is the polar bear (Ursus maritimus), which undergoes population decline in years with low levels of arctic ice coverage due to reduced reproduction and adult female survivorship (Hunter et al., 2010; Laidre et al., 2020). Projections of continued reductions in arctic ice area indicate that the polar bear will experience severe population declines by 2100 (Hunter et al., 2010; Molnár et al., 2011). Climate change is also predicted to impact tropical apex predators. For example, in the mangrove Sundarbans of southern Bangladesh, sea-level rise is predicted to eliminate all suitable habitat for the Bengal tiger (Panthera tigris tigris) by 2070 (Mukul et al., 2019). Climate change impacts are also likely to extend to reptile apex predators, including the world’s largest lizard, the Komodo dragon (Varanus komodoensis). The Komodo dragon occurs on five islands in Indonesia and climate change under moderate emissions scenarios is predicted to reduce suitable habitat by ∼90% by 2050 (Jones et al., 2020). From these examples, it is apparent that climate change will impact apex predators via several pathways and across multiple biomes.

Crocodilians, with 25 extant species currently recognized, function as a crucial apex predator in semiaquatic ecosystems, and are also anticipated to be negatively impacted by climate change. For example, 37-year observations in Florida showed that Crocodylus acutus hatching shifted to earlier dates by 1.5 days, every 2 years with increased sea surface temperature (Cherkiss et al., 2020). Similar results show that increased temperature led to shorter incubation periods for the same species in Mexico (Charruau et al., 2017). Moreover, the sex ratio is determined, in crocodilian species, by incubation temperature, and climate warming is expected to interfere with this (Maciejewski, 2006; González et al., 2019; Bock et al., 2020). For Alligator mississippiensis in Florida, it is estimated that a temperature rise by 1.1–1.4°C in 2040–2050 may skew the sex ratio initially to 95.6% males and then to 97.8% females with a temperature rise by 1.6–3.2°C in 2090–2100 (Bock et al., 2020).

Global warming also can be a threat to crocodilians by destroying their habitat, in particular freshwater swamps or floodplains, through saltwater inundation (SWI) as a result of the sea level rise (SLR). Most species require freshwater habitat for breeding and nesting and such areas typically lie at a low elevation along coasts or rivers and are, therefore, vulnerable to imminent SWI (Pezeshki et al., 1990; Mulrennan and Woodroffe, 1998; Pettit et al., 2018). Despite the adaptation to the saline environment, as implied by their common name, saltwater crocodiles, C. porosus is one such species and requires constant or regular access to freshwater for breeding.

Here, we examine the impacts of climate change on C. porosus, the largest extant crocodilian species, and quantify the potential loss of nesting habitat through SWI and SLR in the Northern Territory (NT), Australia. This region is of global significance, as it supports the most extensive freshwater wetlands and floodplains and contains the largest population of this apex predator in the world (Webb et al., 2010; Fukuda et al., 2021).

Materials and Methods

Study Area

The study area is the northern coastal region of the NT, Australia, within the latitude range −11.0 and −17.0°, called the Top End (Figure 1). The Top End includes the Kakadu Region, which largely consists of the Kakadu National Park. Four major tidal rivers (East, South, West Alligator Rivers and Wildman River) feed into extensive freshwater floodplains contained within the Kakadu National Park. The climate is tropical and monsoonal with distinct dry (May–October) and wet (November–April) seasons. In the coastal areas of the study region, the daily rainfall can exceed 200 mm, and averages approximately 25 mm daily at the peak wet season (Bureau of Meteorology, 2020). The annual rainfall typically ranges between 1,500 and 2,000 mm. The mean maximum and minimum monthly ambient temperature is approximately 17 and 34°C, respectively.

FIGURE 1

Study Species

C. porosus is the largest extant crocodilian species, with the largest individuals exceeding a total length of 6 m and weighing over 1,000 kg (Britton et al., 2012). The species is physiologically adapted to both the freshwater and saline environment (Grigg et al., 1980; Grigg, 1981; Taplin and Grigg, 1981) although they require freshwater for breeding (Webb et al., 1977, 1983). They are found in many different waterbodies including beaches, estuaries, lakes, rivers, and swamps. Some individuals have been reported in the sea far from the shore (Brackhane et al., 2018; Spennemann, 2021). Females build a mound-like nest from vegetation such as tall grasses in the freshwater floodplains or swamps and lay typically 40–55 eggs inside the nest (Webb et al., 1977, 1983; Fukuda and Cuff, 2013). Eggs are incubated by the heat generated by the decomposition of the vegetation material and hatch after approximately 75–95 days (Webb et al., 1983; Richardson et al., 2002). Their breeding is annual and highly seasonal, and is restricted to the wet season.

The species was heavily hunted for commercial use between the 1940s and 1960s, but since legislative protection in 1971, has substantially recovered (Fukuda et al., 2011). They are not considered threatened at any level in the NT and are categorized as Least Concern under the Territory Parks and Wildlife Conservation Act 1976. The Australian population is listed in Appendix II under the Convention on the International Trade in Endangered Species (CITES, IUCN, 2012). As part of the sustainable harvest program implemented in the NT, up to 70,000 eggs are collected annually across the study area for the commercial ranching program (Saalfeld et al., 2015, 2016).

Nesting Habitat Modeling

We used the software package Maxent version 3.4.4 (Phillips et al., 2021) to estimate the current and future areas of suitability for saltwater crocodile nesting in the Top End. Maxent is a presence-only model that minimizes the relative entropy of estimated probability densities between presence points and the background landscape (Elith et al., 2011) and has frequently outperformed other distribution modeling techniques (Hernandez et al., 2006; Wisz et al., 2008).

For the model of suitable nesting habitat at the present time, we collated 730 individual nest locations from across the study area. The crocodile nests were located and harvested by multiple commercial operators during the wet season in 2019 (November–December) and 2020 (January–May) as part of the sustainable harvest program (Saalfeld et al., 2015). We started with an in initial list of 15 environmental variables anticipated to limit saltwater crocodile nesting habitat based on the literature (Table 1). We converted all vector topographic variables to rasters, matching the grain size of a 3 s (∼90 m) digital elevation model (DEM), and masked existing raster layers to the same grain size as the DEM, in ArcGIS version 10.6.1 (esri, 2021). We screened variables for collinearity using the correlation matrix in ArcMap, and excluded one of variables in a pair with a correlation coefficient of ≥ 0.7 (Merow et al., 2013). This left us with the following eight variables: (1) “elevation,” 3 s (∼90 m) digital elevation model; three climate variables from WorldClim downloaded at 30 s resolution including (2) “BIO06,” minimum temperature of coldest month; (3) “BIO15,” precipitation seasonality (coefficient of variation); and (4) “early wet season rainfall,” October–December rainfall (Booth et al., 2014; WorldClim, 2020); and five variables from a Northern Territory 1:250,000 topographic map (Geoscience Australia, 2006), including (5) “floodplain,” distance to land subject to inundation; (6) “perennial lakes,” distance to perennial lakes; (7) “perennial streams,” distance to perennial watercourse; and (8) “freshwater swamp,” categorical value of freshwater swamp (1) vs. other landform type (0) (Table 1). We made this variable categorical because, according to the definitions by Geoscience Australia (2006), the freshwater swamps are small and sparsely distributed in some catchments across the study area. We fitted a nesting habitat suitability model in Maxent with only the linear relationships feature option selected (Merow et al., 2013) and retained 20% of records for validation testing. In addition to the standard cloglog suitability raster output, which we interpreted as relative nesting habitat suitability (Phillips and Dudík, 2008), we also selected a threshold (Maximum training sensitivity plus specificity; Bean et al., 2012) to assign areas of the study area as suitable vs. unsuitable based on the fitted model. To predict suitable nesting habitat in the future, we applied a projection in Maxent using the WorldClim future weather climate projections for the three climate variables. We selected the period 2081–2100 using BCC-CSM2-MR, a medium-resolution global climate model developed by Wu et al. (2019) under two climate scenarios based on low Shared Socio-economic Pathways (SSP 126) and high (SSP 585) future emissions (Hausfather, 2018; WorldClim, 2020).

Sea Level-Rise in Kakadu

Detailed tidally driven, hydrodynamic models of SWI driven by SLR have previously been developed for the Kakadu region (Bayliss et al., 2016), enabling us to make detailed predictions for changes in habitat suitability in this important region. Our thresholded Maxent model produced a raster with binary nesting suitability (0 for unsuitable and 1 for suitable). We overlaid this raster of the thresholded model with the raster datasets of the coastal and river freshwater floodplains simulated by Bayliss et al. (2016), using ArcGIS version 10.6.1. The sea level around the Australian coastlines is expected to rise in a range of 0.75–1.90 m with the mid-range value of 1.1 by 2100 (Short and Woodroffe, 2009). We used the SWI simulation with 1.1 m SLR from Bayliss et al. (2016) to estimate how much of the suitable nesting habitat in the Kakadu Region would be affected by SWI in 2100. Although C. porosus prefers the freshwater environment for nesting, some areas with saline vegetation such as Halosarcia, Tecticornia, and Suaeda are sometimes used because they largely become freshwater in the breeding wet season due to the large input of flushing rainwater during monsoonal events (Fukuda and Cuff, 2013). Thus, we considered that floodplains with less than 0.25 m SWI would remain as habitat suitable for nesting and those areas with SWI of more than 0.25 m are unsuitable. We used the raster predictions from Bayliss et al. (2016) for > 0.25 m SWI to identify areas that will be lost from our thresholded model of current suitable nesting habitat by the year 2100.

Results

Our model of current saltwater crocodile nesting habitat performed well (test AUC of 0.958 ± SD 0.003) and showed suitable areas around the coasts and floodplains of the Top End of the NT (Figure 1). The most important variable was distance to perennial lakes (71.0% contribution, 87.5% permutation importance), with a negative logistic relationship with nesting habitat suitability (Figure 2A). Early wet season rainfall was the next most important variable (12.9% contribution; 0.9% permutation importance), with a positive logistic relationship between Early wet season rainfall and nesting habitat suitability (Figure 2B). The freshwater swamp variable (8.3% contribution; 1.2% permutation importance) showed much higher suitability than other landform types (Figure 2C). Seasonality in rainfall (BIO15) showed 6.6% contribution; 3.0% permutation importance, with a negative logistic relationship to nesting habitat suitability (Figure 2D). All other variables contributed < 5% to the model.

FIGURE 2

For the whole of the Top End, our thresholded Maxent model predicted a total of 32,306.91 km² of suitable nesting habitat across the study area for the current time. For the period 2081–2100, using the predicted future temperature and rainfall patterns derived from the BCC-CSM2-MR climate model but retaining the other variables at their present values, our Maxent model predicted an area of suitable habitat of 32,199.32 km² under the low emissions scenario (SSP 126) and 42,740.78 km² under the high emissions scenario (SSP 585) across the study area.

Within the Kakadu Region, a total area of 4232.08 km² was assigned as suitable nesting habitat by our current model (Figure 3A). However, the hydrodynamic models by Bayliss et al. (2016) show that as at 2013 sea level, approximately 16.40% of the suitable habitat was already affected by < 0.25 m SWI, and 3538.10 km² remains suitable for nesting (Figure 3B). After 1.1 m SLR (that is, sea levels forecast for 2100), only 1775.70 km² remains suitable (Figure 3C). This represents a 58.04% reduction in predicted suitable nesting habitat over 87 years, between 2013 and 2100).

FIGURE 3

Discussion

Our model of suitable nesting habitat for the saltwater crocodile revealed a dominant importance of abundant, perennial lakes in Australia’s Top End. The estimate of nesting suitability rapidly dropping beyond a 50 m range from perennial lakes (Figure 2A) is consistent with field observations that a vast majority of nests are made around the edges of waterbodies during the breeding season (Webb et al., 1977, 1983; Fukuda and Cuff, 2013). The perennial lakes appear to be important as females guarding their nests during the incubation period, typically 2–3 months, require access to freshwater (Webb et al., 1977; Magnusson, 1980). The model also identified the importance of early rains in the wet season (Figure 2B), supporting previous findings that rains in early wet season (Oct–Dec) trigger the highly seasonal breeding of C. porosus, and higher precipitation during this period results in higher abundance of hatchlings in the following dry season (Fukuda and Saalfeld, 2014).

We found that projected changes in temperature and rainfall under climate change could impact the area of suitable nesting habitat negatively or positively (0.33% decrease under a low emission scenario SSP 126 and 32.30% increase under a high emission scenario SSP 585). It is worth mentioning that, apart from the adverse impact of SLR, the projected nesting habitat would otherwise be increased to some extent under the higher future emission scenario, because of some positive effects on the climate variables such as increased early rains and decreased seasonality. Nevertheless, these effects will be canceled by the much higher, negative impact by SLR.

Most importantly, our analysis for the Kakadu Region showed that almost 50% of the suitable nesting habitat in 2013 would be lost to the 1.1 m SLR by 2100 (Figure 3). While fine-scale SLR data and forecasts are not currently available outside of the Kakadu region, if a similar proportion of the habitat is affected across the larger study area, 16,522 km² would be lost to SLR across the Top End. This is of significant concern as the majority of the Australian population of C. porosus resides in this area (Fukuda et al., 2007, 2021). This suggests that SLR by climate change may represent the major threat to the species. Other crocodilians that build mound nests in freshwater wetlands or floodplains, such as A. mississippiensis and C. moreletti, could also be at similar risk, as their nests are prone to flooding (Kushlan and Jacobsen, 1990; Platt et al., 2021).

One important caveat to the predicted loss of saltwater crocodile nesting habitat is that the creation of new freshwater habitats in the Kakadu region has not been accounted for. Bayliss et al. (2016) were unable to predict new areas of freshwater habitat created through SLR because the Lidar-derived fine-scale elevation data did not extend beyond the current floodplain. Undoubtedly, SLR will result in the formation of some areas of new habitat suitable for saltwater crocodile nesting which may somewhat offset the areas lost. However, the sandstone escarpments that occur further inland (Northern Territory and Geological Survey, 2006) form a physical barrier to the formation of extensive freshwater floodplains. In addition, while small billabongs or lakes may form on the escarpment, the steep and rocky terrain curtails the dispersal of saltwater crocodiles, and to date, they have never been observed in these areas (Letnic and Connors, 2006). For areas outside of the sandstone escarpments, the potential creation of new floodplain habitat is an important avenue for further research. Additional data such as LiDAR (Dong and Chen, 2017) would facilitate more accurate forecasts of potential future suitable habitat.

Another source of uncertainty is the degree of future SLR. Bayliss et al. (2016) adopted a value of 1.1 m in 2001 based on a range of published projections (Solomon et al., 2007; Jevrejeva et al., 2010). The IPCC recently released their 6th Assessment Report which includes predictions of 0.63–1.01 m by 2100 under the moderate emission scenario and 0.66–1.33 m by 2150 under the very high emissions scenario (Masson-Delmotte et al., 2021). Therefore, the 1.1 m SLR used in our study represents likely SLR between 2100 and 2150. However, several studies document the acceleration of SLR in recent decades (Nerem et al., 2018; Dangendorf et al., 2019), and in addition, the IPCC stated that due to the substantial uncertainty in global ice sheet dynamics, SLR approaching 2 m by 2100 and 5 m by 2150 cannot be ruled out (Masson-Delmotte et al., 2021). The substantial uncertainty in the potential creation of new habitat and the degree of SLR means that estimates of the loss of saltwater crocodile nesting habitat should be recalculated as new data become available.

It should be noted that the change from a freshwater to more saline habitat is not anticipated to be monotonic, and gradual replacement of freshwater plant species with those that are more saline tolerant is the likely scenario. While the floodplain vegetation is determined by fine-scale variation in topography (Finlayson et al., 2006) and the global mean SLR is 4.8 mm annually (Dangendorf et al., 2019), even infrequent SWI can lead to sudden changes in salinity (Finlayson et al., 2006; Bayliss et al., 2016). Freshwater plant taxa have a unique tolerance to salinity and water depth (Cowie, 2003), and most freshwater plant species are unable to tolerate salinity much in excess of 1 ppt (Pettit et al., 2018). Previous work has shown an association between crocodile nesting and Oryza dominated tall tussock and Melaleuca open grassland with an understory comprised of species also present in the Oryza dominated grasslands (Fukuda and Cuff, 2013), and Oryza species are unable to tolerate salinity above 1 ppt (Pettit et al., 2018). In contrast, mangrove species can tolerate salinity to 50 ppt (Ball, 1998), and are not obligate halophiles, that is, they are also able to live in freshwater conditions (Woodroffe, 1988).

Although fossil records and molecular analysis indicate that extant crocodilians were capable of tracking changes in their distribution in response to drastic climate changes in the last 100 million years (Roos et al., 2007; Brochu et al., 2009; Ryberg and Lawing, 2018), whether they can continue to track changes in the spatial distribution and quality of suitable nesting habitat is unknown. Crocodylus porosus is one of the oldest extant species in the world, little changed since dinosaurs roamed the planet, yet their survival during the Anthropocene may now depend on how quickly humanity can transition away from the burning of fossil fuels and thus stabilize global warming.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

• 1

  BallM. (1998). Mangrove species richness in relation to salinity and waterlogging: a case study along the Adelaide River floodplain, northern Australia.Glob Ecol. Biogeogr. Lett.773–82. 10.1111/j.1466-8238.1998.00282.x

  • CrossRef
  • Google Scholar

• 2

  BaylissP.SaundersK.DutraL. X. C.MeloL. F. C.HiltonJ.PrakashM.et al (2016). Assessing sea level-rise risks to coastal floodplains in the Kakadu Region, northern Australia, using a tidally driven hydrodynamic model.Mar. Freshwater Res.691064–1078. 10.1071/MF16049

  • CrossRef
  • Google Scholar

• 3

  BeanW. T.StaffordR.BrasharesJ. S. (2012). The effects of small sample size and sample bias on threshold selection and accuracy assessment of species distribution models.Ecography35250–258. 10.1111/j.1600-0587.2011.06545.x

  • CrossRef
  • Google Scholar

• 4

  BockS. L.LowersR. H.RainwaterT. R.StolenE.DrakeJ. M.WilkinsonP. M.et al (2020). Spatial and temporal variation in nest temperatures forecasts sex ratio skews in a crocodilian with environmental sex determination.Proc.R. Soc. B.Biol. Sci.287:20200210. 10.1098/rspb.2020.0210

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5

  BoothT. H.NixH. A.BusbyJ. R.HutchinsonM. F. (2014). bioclim: the first species distribution modelling package, its early applications and relevance to most current MaxEnt studies.Divers. Distrib.201–9. 10.1111/ddi.12144

  • CrossRef
  • Google Scholar

• 6

  BrackhaneS.WebbG.XavierF. M. E.GusmaoM.PechacekP. (2018). When Conservation Becomes Dangerous: Human-Crocodile Conflict in Timor-Leste.J. Wildl. Manag.821332–1344. 10.1002/jwmg.21497

  • CrossRef
  • Google Scholar

• 7

  BrittonA. R. C.WhitakerR.WhitakerN. (2012). Here be a dragon: exceptional size in a saltwater crocodile (Crocodylus porosus) from the Philippines.Herpetol. Rev.43541–546.

  • Google Scholar

• 8

  BrochuC. A.WagnerJ. R.JouveS.SumrallC. D.DensmoreL. D. (2009). A Correction Corrected: Consensus Over the Meaning of Crocodylia and Why It Matters.Syst. Biol.58537–543. 10.1093/sysbio/syp053

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9

  Bureau of Meteorology. (2020). Climate Data Online. Available online at: http://www.bom.gov.au/climate/data/index.shtml(Accessed on February 20, 2020)

  • Google Scholar

• 10

  CharruauP.CantónD.Méndez-de-la-CruzF. (2017). Additional details on temperature-dependent sex determination in Crocodylus acutus.Salamandra53304–308.

  • Google Scholar

• 11

  CherkissM. S.WatlingJ. I.BrandtL. A.MazzottiF. J.LindsayJ.BeauchampJ. S.et al (2020). Shifts in hatching date of American crocodile (Crocodylus acutus) in southern Florida.J. Thermal Biol.88:102521. 10.1016/j.jtherbio.2020.102521

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12

  CowieI. (2003). Freshwater aquatic plants of Darwin Harbour catchments. in Proceedings of Darwin Harbour Region: Current Knowledge and Future Needs. Public Presentations. (Darwin, AUS: Working Group for the Darwin Harbour Advisory Committee), 160–177.

  • Google Scholar

• 13

  DangendorfS.HayC.CalafatF. M.MarcosM.PiecuchC. G.BerkK.et al (2019). Persistent acceleration in global sea-level rise since the 1960s.Nat. Clim. Chang.9705–710. 10.1038/s41558-019-0531-8

  • CrossRef
  • Google Scholar

• 14

  Department of the Environment and Energy. (2013). Code of Practice for the Humane Treatment of Wild and Farmed Australian Crocodiles. Available online at: http://www.environment.gov.au/resource/code-practice-humane-treatment-wild-and-farmed-australian-crocodiles(Accessed on February 28, 2014)

  • Google Scholar

• 15

  DongP.ChenQ. (2017). LiDAR Remote Sensing and Applications. (Boca Raton, FL: CRC Press).

  • Google Scholar

• 16

  ElithJ.PhillipsS. J.HastieT.DudíkM.CheeY. E.YatesC. J. (2011). A statistical explanation of MaxEnt for ecologists.Divers. Distrib.1743–57. 10.1111/j.1472-4642.2010.00725.x

  • CrossRef
  • Google Scholar

• 17

  esri. (2021). ArcGIS. Available online at: https://www.esri.com/en-us/home(accessed February 17, 2022).

  • Google Scholar

• 18

  FinlaysonC. M.LowryJ.BellioM. G.NouS.PidgeonR.WaldenD.et al (2006). Biodiversity of the wetlands of the Kakadu Region, northern Australia.Aquat. Sci.68374–399. 10.1007/s00027-006-0852-3

  • CrossRef
  • Google Scholar

• 19

  FukudaY.CuffN. (2013). Vegetation communities as nesting habitat for the saltwater crocodiles in the Northern Territory of Australia.Herpetol. Conservat. Biol.8641–651.

  • Google Scholar

• 20

  FukudaY.McDonaldP. (2022). Maxent_results_crocodile_nesting_suitability [Dataset]. Figshare. Available online at: https://doi.org/10.6084/m9.figshare.17292191.v1(accessed February 18, 2022).

  • Google Scholar

• 21

  FukudaY.SaalfeldK. (2014). Abundance of Saltwater Crocodile Hatchlings is Related to Rainfall in the Preceding Wet Season in Northern Australia.herp.70439–448. 10.1655/HERPETOLOGICA-D-13-00096R3

  • CrossRef
  • Google Scholar

• 22

  FukudaY.WebbG.EdwardsG.SaalfeldK.WhiteheadP. (2021). Harvesting predators: simulation of population recovery and controlled harvest of saltwater crocodiles Crocodylus porosus.Wildl. Res.48252–263. 10.1071/WR20033

  • CrossRef
  • Google Scholar

• 23

  FukudaY.WebbG.ManolisC.DelaneyR.LetnicM.LindnerG.et al (2011). Recovery of saltwater crocodiles following unregulated hunting in tidal rivers of the Northern Territory.Aus. J. Wildl. Manag.751253–1266. 10.1002/jwmg.191

  • CrossRef
  • Google Scholar

• 24

  FukudaY.WhiteheadP.BoggsG. (2007). Broad-scale environmental influences on the abundance of saltwater crocodiles (Crocodylus porosus) in Australia.Wildl. Res.34167–176. 10.1071/WR06110

  • CrossRef
  • Google Scholar

• 25

  Geoscience Australia. (2006). GEODATA TOPO 250K Series 3. (Canberra, AUS: Australian Government)

  • Google Scholar

• 26

  GonzálezE. J.Martínez-LópezM.Morales-GarduzaM. A.García-MoralesR.CharruauP.Gallardo-CruzJ. A. (2019). The sex-determination pattern in crocodilians: A systematic review of three decades of research.J. Animal Ecol.881417–1427. 10.1111/1365-2656.13037

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27

  GriggG. C. (1981). Plasma homeostasis and cloacal urine composition in Crocodylus porosus caught along a salinity gradient.J. Comp. Physiol. B.144261–270. 10.1007/BF00802765

  • CrossRef
  • Google Scholar

• 28

  GriggG. C.TaplinL. E.HarlowP.WrightJ. (1980). Survival and growth of hatchling Crocodylus porosus in saltwater without access to fresh drinking water.Oecologia47264–266. 10.1007/BF00346830

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29

  HausfatherZ. (2018). Explainer: How ‘Shared Socioeconomic Pathways’ explore future climate change.Carbon Brief. Available online atat: https://www.carbonbrief.org/explainer-how-shared-socioeconomic-pathways-explore-future-climate-change(Accessed on December 13, 2021)

  • Google Scholar

• 30

  HernandezP. A.GrahamC. H.MasterL. L.AlbertD. L. (2006). The effect of sample size and species characteristics on performance of different species distribution modeling methods.Ecography29773–785. 10.1111/j.0906-7590.2006.04700.x

  • CrossRef
  • Google Scholar

• 31

  HunterC. M.CaswellH.RungeM. C.RegehrE. V.AmstrupS. C.StirlingI. (2010). Climate change threatens polar bear populations: a stochastic demographic analysis.Ecology912883–2897. 10.1890/09-1641.1

  • CrossRef
  • Google Scholar

• 32

  IUCN (2012). Guidelines for application of IUCN Red List criteria at regional and national levels: version 4.0. (Gland, Switzerland: IUCN Species Survival Commission). Available at: http://www.iucnredlist.org/documents/reg_guidelines_en.pdf(Accessed on January 22, 2014)

  • Google Scholar

• 33

  JevrejevaS.MooreJ. C.GrinstedA. (2010). How will sea level respond to changes in natural and 394 anthropogenic forcings by 2100?Geophys. Res. Lett.37:L07703. 10.1029/2010GL042947

  • CrossRef
  • Google Scholar

• 34

  JonesA. R.JessopT. S.AriefiandyA.BrookB. W.BrownS. C.CiofiC.et al (2020). Identifying island safe havens to prevent the extinction of the World’s largest lizard from global warming.Ecol. Evol.1010492–10507. 10.1002/ece3.6705

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  KushlanJ. A.JacobsenT. (1990). Environmental Variability and the Reproductive Success of Everglades Alligators.J. Herpetol.24176–184. 10.2307/1564225

  • CrossRef
  • Google Scholar

• 36

  LaidreK. L.AtkinsonS.RegehrE. V.SternH. L.BornE. W.WiigØet al (2020). Interrelated ecological impacts of climate change on an apex predator.Ecol. Appl.30:e02071. 10.1002/eap.2071

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37

  LetnicM.ConnorsG. (2006). Changes in the distribution and abundance of saltwater crocodiles (Crocodylus porosus) in the upstream, freshwater reaches of rivers in the Northern Territory.Aus. Wildl. Res.33529–538.

  • Google Scholar

• 38

  MaciejewskiK. (2006). Temperature-dependent sex determination in the Nile crocodile, Crocodylus niloticus, in the Okavango River, Botswana, and the effect of global climate change. Available online at: https://scholar.sun.ac.za:443/handle/10019.1/50648(Accessed on October. 29, 2021).

  • Google Scholar

• 39

  MagnussonW. E. (1980). Habitat Required for Nesting by Crocodylus porosus (Reptilia?: Crocodilidae).Wildl. Res.7149–156. 10.1071/wr9800149

  • CrossRef
  • Google Scholar

• 40

  Masson-DelmotteV.ZhaiP.PiraniA.ConnorsS. L.PéanC.BergerS.et al (2021). Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change.(Cambridge: Cambridge University Press).

  • Google Scholar

• 41

  MerowC.SmithM. J.SilanderJ. A.Jr. (2013). A practical guide to MaxEnt for modeling species’ distributions: what it does, and why inputs and settings matter.Ecography361058–1069. 10.1111/j.1600-0587.2013.07872.x

  • CrossRef
  • Google Scholar

• 42

  MolnárP. K.DerocherA. E.KlanjscekT.LewisM. A. (2011). Predicting climate change impacts on polar bear litter size.Nat. Commun.2:186. 10.1038/ncomms1183

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43

  MukulS. A.AlamgirM.SohelM. S. I.PertP. L.HerbohnJ.TurtonS. M.et al (2019). Combined effects of climate change and sea-level rise project dramatic habitat loss of the globally endangered Bengal tiger in the Bangladesh Sundarbans.Sci. Total Environ.663830–840. 10.1016/j.scitotenv.2019.01.383

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44

  MulrennanM. E.WoodroffeC. D. (1998). Saltwater intrusion into the coastal plains of the Lower Mary River. Northern Territory, Australia.J. Environ. Manag.54169–188. 10.1006/jema.1998.0229

  • CrossRef
  • Google Scholar

• 45

  NeremR. S.BeckleyB. D.FasulloJ. T.HamlingtonB. D.MastersD.MitchumG. T. (2018). Climate-change–driven accelerated sea-level rise detected in the altimeter era.PNAS1152022–2025. 10.1073/pnas.1717312115

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46

  Northern Territory, and Geological Survey. (2006). Geological Map of the Northern Territory. Available online at: https://geoscience.nt.gov.au/gemis/ntgsjspui/handle/1/82057(Accessed on November 2, 2021).

  • Google Scholar

• 47

  PettitN. E.BaylissP.BartoloR. (2018). Dynamics of plant communities and the impact of saltwater intrusion on the floodplains of Kakadu National Park.Mar. Freshwater Res.691124–1133. 10.1071/MF16148

  • CrossRef
  • Google Scholar

• 48

  PezeshkiS. R.DelauneR. D.PatrickW. H. (1990). Flooding and saltwater intrusion: Potential effects on survival and productivity of wetland forests along the U.S.Gulf Coast. Forest Ecol.Manag.3287–301. 10.1016/0378-1127(90)90199-L

  • CrossRef
  • Google Scholar

• 49

  PhillipsS. J.DudíkM. (2008). Modeling of species distributions with Maxent: new extensions and a comprehensive evaluation.Ecography31161–175. 10.1111/j.0906-7590.2008.5203.x

  • CrossRef
  • Google Scholar

• 50

  PhillipsS. J.DudikM.SchapireR. E. (2021). Maxent software for modeling species niches and distributions. (New York, NY:American Museum of Natural History).

  • Google Scholar

• 51

  PlattS. G.RainwaterT. R.McMurryS. T. (2021). Fauna associated with the nests of Crocodylus moreletii and Crocodylus moreletii × acutus in Belize.J. Nat. His.55133–149. 10.1080/00222933.2021.1895350

  • CrossRef
  • Google Scholar

• 52

  RichardsonK. C.WebbG. J. W.ManolisS. C. (2002). Crocodiles: Inside Out. A Guide to the Crocodilians and their Functional Morphology. (Sydney, Aus: Surrey Beatty & Sons).

  • Google Scholar

• 53

  RoosJ.AggarwalR. K.JankeA. (2007). Extended mitogenomic phylogenetic analyses yield new insight into crocodylian evolution and their survival of the Cretaceous–Tertiary boundary.Mol. Phylogene. Evol.45663–673. 10.1016/j.ympev.2007.06.018

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54

  RybergW.LawingM. (2018). “Genetic Consequences and Management Implications of Climate Change for the American Alligator (Alligator mississippiensis),” in American Alligators: Habitats, Behaviors, and Threats, edsHenkeS.EversoleC.123–153. (Hauppauge, NY: Nova Science Publishers).

  • Google Scholar

• 55

  SaalfeldK.FukudaY.DuldigT.FisherA. (2015). Wildlife Trade Management Plan for the Saltwater Crocodile (Crocodylus porosus) in the Northern Territory of Australia, 2016-2020. (Darwin, AUS: Northern Territory Department of Land Resource Management).

  • Google Scholar

• 56

  SaalfeldK.FukudaY.DuldigT.FisherA. (2016). Management Program for the Saltwater Crocodile in the Northern Territory of Australia, 2016-2020. (Darwin, Australia: Northern Territory Department of Environment and Natural Resources). Available at: https://nt.gov.au/__data/assets/pdf_file/0007/443581/crocodile-management-program.pdf(accessed February 17, 2022).

  • Google Scholar

• 57

  ShortA.WoodroffeC. (2009). The coast of Australia. Available online at: https://ro.uow.edu.au/scipapers/4162(accessed February 17, 2022).

  • Google Scholar

• 58

  SolomonS.IpccIntergovernmental Panel on Climate Change. (2007). Climate change 2007: the physical science basis: contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. (Cambridge?; NY: Cambridge University Press)

  • Google Scholar

• 59

  SpennemannD. H. R. (2021). Cruising the currents: Observations of extra-limital saltwater crocodiles (Crocodylus porosus Schneider, 1801) in the Pacific Region.Pacific Science74211–227. 10.2984/74.3.1

  • CrossRef
  • Google Scholar

• 60

  TaplinL. E.GriggG. C. (1981). Salt Glands in the Tongue of the Estuarine Crocodile Crocodylus porosus.Science2121045–1047. 10.1126/science.212.4498.1045

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61

  WebbG. J. W.ManolisS. C.BrienM. L. (2010). “Saltwater Crocodile Crocodylus porosus,” in Crocodiles Status Survey and Conservation Action Plan, edsManolisS. C.StevensonC. (Darwin, AUS: Crocodile Specialist Group), 99–113. Available online at: http://www.iucncsg.org/365_docs/attachments/protarea/18%20–8088e67a.pdf(accessed February 17, 2022).

  • Google Scholar

• 62

  WebbG. J. W.ManolisS. C.BuckworthR.SackG. C. (1983). An Examination of Crocodylus porosus nests in two northern Australian freshwater swamps, with an analysis of embryo mortality.Wildl. Res.10571–605. 10.1071/wr9830571

  • CrossRef
  • Google Scholar

• 63

  WebbG. J. W.MesselH.MagnussonW. E. (1977). The nesting biology of Crocodylus porosus in Arnhem Land, northern Australia.Copeia1977238–249. 10.2307/1443905

  • CrossRef
  • Google Scholar

• 64

  WiszM. S.HijmansR. J.LiJ.PetersonA. T.GrahamC. H.GuisanA.et al (2008). Effects of sample size on the performance of species distribution models.Divers. Distrib.14763–773. 10.1111/j.1472-4642.2008.00482.x

  • CrossRef
  • Google Scholar

• 65

  WoodroffeC. D. (1988). Relict Mangrove Stand on Last Interglacial Terrace. Christmas Island, Indian Ocean.J. Tropical Ecol.41–17.

  • Google Scholar

• 66

  WorldClim. (2020). Bioclimatic variables — WorldClim 1 documentation. Available online at: https://www.worldclim.org/data/bioclim.html(Accessed on October 18, 2021)

  • Google Scholar

• 67

  WuT.LuY.FangY.XinX.LiL.LiW.et al (2019). The Beijing Climate Center Climate System Model (BCC-CSM): the main progress from CMIP5 to CMIP6.Geosci. Model Devel.121573–1600. 10.5194/gmd-12-1573-2019

  • CrossRef
  • Google Scholar