ORIGINAL RESEARCH article

Front. Ecol. Evol., 18 April 2022

Sec. Behavioral and Evolutionary Ecology

Volume 10 - 2022 | https://doi.org/10.3389/fevo.2022.839423

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

  • 1. Department of Environment, Parks and Water Security, Northern Territory Government, Darwin, NT, Australia

  • 2. Complex Systems Consulting, Brisbane, QLD, Australia

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 km2 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 km2 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 (; ). Projections of continued reductions in arctic ice area indicate that the polar bear will experience severe population declines by 2100 (; ). 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 (). 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 (). 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 (). Similar results show that increased temperature led to shorter incubation periods for the same species in Mexico (). Moreover, the sex ratio is determined, in crocodilian species, by incubation temperature, and climate warming is expected to interfere with this (; ; ). 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 ().

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 (; ; ). 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 (; ).

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 (). 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

TABLE 1

VariableDescriptionSource
1ElevationMeters above the average sea levelDEM 3 s (Geoscience Australia, 2010)
2BIO06Minimum temperature of the coldest month in a yearBIOCLIM (WorldClim, 2021)
3BIO15Coefficient of variation in monthly precipitation expressed as a percentageBIOCLIM (WorldClim, 2021)
4Early wet season rainfallMean total rainfall (mm) in October–DecemberBIOCLIM (WorldClim, 2021)
5FloodplainDistance (km) to the closest land subject to inundationGEODATA TOPO 250 series 3 (Geoscience Australia, 2007)
6Perennial lakesDistance (km) to the closest perennial lakesGEODATA TOPO 250 series 3 (Geoscience Australia, 2007)
7Perennial streamsDistance (km) to the closest perennial watercourseGEODATA TOPO 250 series 3 (Geoscience Australia, 2007)
8Freshwater swampCategorical value of freshwater swamp vs. other landformGEODATA TOPO 250 series 3 (Geoscience Australia, 2007)

Biological relevance

1Nesting habitat are found at low elevation. *,**
2Distribution is limited by the minimum temperature.*
3Distribution is limited to the monsoonal climate and breeding occurs in the wet season.*
4Higher recruitment occurs after the west season with more rainfall in October–December.***
5Freshwater floodplains are major nesting habitat. *
6Perennial lakes are major residential habitat for females, connecting to breeding sites.*
7Perennial watercourses are major residential habitat for females, connecting to breeding sites.*
8Patchy freshwater swamps are residential and nesting habitat.*

Environmental variables and their attributes used for the Maxent models.

*, **, ***.

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 (). The species is physiologically adapted to both the freshwater and saline environment (; ; ) although they require freshwater for breeding (, ). 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 (; ). 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 (, ; ). Eggs are incubated by the heat generated by the decomposition of the vegetation material and hatch after approximately 75–95 days (; ). 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 (). 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, ). 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 (, ).

Nesting Habitat Modeling

We used the software package Maxent version 3.4.4 () 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 () and has frequently outperformed other distribution modeling techniques (; ).

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 (). 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 (). 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 (). 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 (; ); and five variables from a Northern Territory 1:250,000 topographic map (), 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 , 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 () 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 (), we also selected a threshold (Maximum training sensitivity plus specificity; ) 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 under two climate scenarios based on low Shared Socio-economic Pathways (SSP 126) and high (SSP 585) future emissions (; ).

Sea Level-Rise in Kakadu

Detailed tidally driven, hydrodynamic models of SWI driven by SLR have previously been developed for the Kakadu region (), 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 , 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 (). We used the SWI simulation with 1.1 m SLR from 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 (). 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 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 km2 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 km2 under the low emissions scenario (SSP 126) and 42,740.78 km2 under the high emissions scenario (SSP 585) across the study area.

Within the Kakadu Region, a total area of 4232.08 km2 was assigned as suitable nesting habitat by our current model (Figure 3A). However, the hydrodynamic models by 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 km2 remains suitable for nesting (Figure 3B). After 1.1 m SLR (that is, sea levels forecast for 2100), only 1775.70 km2 remains suitable (Figure 3C). This represents a 58.04% reduction in predicted suitable nesting habitat over 87 years, between 2013 and 2100).

FIGURE 3

, and (C) 1.1 m SLR and excluding areas with > 0.25 m SWI as simulated for 2100 by .

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 (, ; ). The perennial lakes appear to be important as females guarding their nests during the incubation period, typically 2–3 months, require access to freshwater (; ). 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 ().

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 km2 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 (, ). 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 (; ).

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. 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 () 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 (). 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 () would facilitate more accurate forecasts of potential future suitable habitat.

Another source of uncertainty is the degree of future SLR. adopted a value of 1.1 m in 2001 based on a range of published projections (; ). 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 (). 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 (; ), 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 (). 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 () and the global mean SLR is 4.8 mm annually (), even infrequent SWI can lead to sudden changes in salinity (; ). Freshwater plant taxa have a unique tolerance to salinity and water depth (), and most freshwater plant species are unable to tolerate salinity much in excess of 1 ppt (). 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 (), and Oryza species are unable to tolerate salinity above 1 ppt (). In contrast, mangrove species can tolerate salinity to 50 ppt (), and are not obligate halophiles, that is, they are also able to live in freshwater conditions ().

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 (; ; ), 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.

Statements

Data availability statement

The original contributions presented in the study are publicly available. This data can be found here: .

Ethics statement

Ethical review and approval was not required for the animal study because this study did not involve live animals and only used the location data of crocodile nests that were reported to the Northern Territory Government, Australia.

Author contributions

YF, PM, and BC contributed to conception and design of the study, and wrote the manuscript. YF organized the data access. YF and PM performed the spatial analysis. All authors contributed to manuscript revision, read, and approved the submitted version.

Funding

Department of Environment, Parks and Water Security provided financial support.

Acknowledgments

Crocodile Farmers Association of the Northern Territory provided the location data of the crocodile nests used in the analyses. These location data were collected in accordance to the Code of practice for the humane treatment of wild and farmed Australian crocodiles () as part of the sustainable harvest program authorized by the Northern Territory Government (). Alaric Fisher, Glenn Edwards, Tim Clancy, and Tony Griffiths supported the study.

Conflict of interest

BC was affiliated with Complex Systems Consulting. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

Summary

Keywords

climate change, sea level rise (SLR), crocodile, habitat suitability analysis, Maxent, saltwater incursion, Kakadu National Park, Australia

Citation

Fukuda Y, McDonald PJ and Crase B (2022) Lost to the Sea: Predicted Climate Change Threats to Saltwater Crocodile Nesting Habitat. Front. Ecol. Evol. 10:839423. doi: 10.3389/fevo.2022.839423

Received

20 December 2021

Accepted

10 February 2022

Published

18 April 2022

Volume

10 - 2022

Edited by

Jeanine M. Refsnider, The University of Toledo, United States

Reviewed by

Lin Zhang, Hubei University of Chinese Medicine, China; Paula A. White, University of California, Los Angeles, United States

Updates

Copyright

*Correspondence: Yusuke Fukuda,

This article was submitted to Behavioral and Evolutionary Ecology, a section of the journal Frontiers in Ecology and Evolution

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics