Mangrove dispersal disrupted by projected changes in global seawater density
Mangrove forests thrive along tropical and subtropical shorelines and their distribution extends to warm temperate regions1. They are globally recognized for the valuable ecosystem services they provide2 but are expected to be substantially influenced by climate change-related physical processes in the future3,4. Under warming winter temperatures, poleward expansion is predicted for mangroves5,6, with potential implications for ecosystem structure and functioning, as well as human livelihoods and well-being7,8. The global distribution, abundance and species richness of mangroves is governed by a broad range of biotic and environmental factors, including temperature and precipitation9 and diverse geomorphological and hydrological gradients10. Climate and aspects related to coastal geography (for example, floodplain area) determine the availability of suitable habitat for establishment11,12. However, the potential for mangroves to track changing environmental conditions and expand their distributions ultimately depends on dispersal11,13. The importance of dispersal in controlling mangrove distributions has been demonstrated by mangrove distributional responses to historical climate variability14, past mangrove (re)colonization of oceanic islands15 and from the long-term survival of mangrove seedlings planted beyond natural range limits16. As such, quantifying changes in the factors that influence dispersal is important for understanding climate-driven distributional responses of mangroves under future climate conditions.In mangroves, dispersal is accomplished by buoyant seeds and fruits (hereafter referred to as ‘propagules’). In combination with prevailing currents, the spatial scale of this process, ranging from local retention to transoceanic dispersal over thousands of kilometres13, is determined by propagule buoyancy17, that is, the density difference between that of propagules and the surrounding water. Hence, the course of dispersal trajectories for propagules from these species depends on the interaction between spatiotemporal changes in both propagule density and that of the surrounding water, rendering this process sensitive to climate-driven changes in coastal and open-ocean water properties. The biogeographic implications of such density differences were recognized more than a century ago by Henry Brougham Guppy, who discussed18 ‘the far-reaching influence on plant-distribution and on plant-development that the relation between the specific weight of seeds and fruits and the density of sea-water must possess’.Since the time of Guppy’s early observations, climate change from human activities has driven pronounced changes in ocean temperature and salinity, with further changes predicted throughout the twenty-first century19. Ocean density is a nonlinear function of temperature, salinity and pressure20; therefore, these changes may influence dispersal patterns of mangrove propagules by altering their buoyancy and floating orientation. As Guppy noted18, ‘[for] plants whose seeds or fruits are not much lighter than seawater […] the effect of increased density of the water is to extend the flotation period’ or ‘to increase the number that floated for a given period’. Guppy also reported that the seedlings of the widespread mangrove genera Rhizophora and Bruguiera present exceptional examples of propagules with densities somewhere between seawater and freshwater18. Previous studies of the impacts of climate change on mangroves have focused on factors such as sea level rise, altered precipitation regimes and increasing temperature and storm frequency4,21,22,23 but the potential impact of climate-driven changes in seawater properties on mangroves has not yet been examined. This is somewhat surprising, as the ocean is the primary dispersal medium of this ‘sea-faring’ coastal vegetation and dispersal is a key process that governs a species’ response to climate change by changing its geographical range. This knowledge gap contrasts with recent efforts to expose links between climate change and dispersal in other ecologically important marine taxa such as zooplankton and fish species24,25,26,27.In this study, we investigate predicted changes in sea surface temperature (SST), sea surface salinity (SSS) and sea surface density (SSD) for coastal waters bordering mangrove forests (hereafter referred to as ‘coastal mangrove waters’), over the next century. Using a biogeographic classification system for coastal and shelf areas28, we examine spatiotemporal changes in these surface ocean properties, with a particular focus on the world’s two major mangrove diversity hotspots: (1) the Atlantic East Pacific (AEP) region, including all of the Americas, West and Central Africa and (2) the Indo West Pacific (IWP) region, extending from East Africa eastwards to the islands of the central Pacific1. Finally, we synthesize available data on the density of mangrove propagules for different mangrove species and explore the potential impact of climate-driven changes in SSD on propagule dispersal.To assess changes in SST and SSS throughout the global range of mangrove forests, we used present (2000–2014) and future (2090–2100) surface ocean properties from the Bio-ORACLE database29,30. SSD estimates were derived from these variables using the UNESCO EOS-80 equation of state polynomial for seawater31. Changes in SST, SSS and SSD (Fig. 1) were calculated for four representative concentration pathways (RCPs) and derived for coastal waters closest to the 583,578 polygon centroids from the 2015 Global Mangrove Watch (GMW) database32. After removing duplicates, our dataset contained 10,108 unique mangrove occurrence locations, with corresponding present conditions and predicted future changes in mean SST, SSS and SSD. Under the low-warming scenario RCP 2.6, mean SST of coastal mangrove waters is predicted to change by +0.64 (±0.11) °C and mean SSS by −0.06 (±0.25) practical salinity units (PSU). Combined, this results in an average change in mean SSD of −0.25 (±0.20) kg m−3 in coastal mangrove waters by the late twenty-first century (Supplementary Table 1). These values roughly double under RCP 4.5 (Supplementary Table 2), while under RCP 6.0, a change of +1.69 (±0.14) °C in mean SST, −0.21 (±0.42) PSU in mean SSS and −0.71 (±0.32) kg m−3 in mean SSD is predicted (Supplementary Table 3). Under RCP 8.5, our study predicts a change in SST of +2.84 (±0.21) °C (range 2.11–4.01 °C), a change in SSS of −0.30 (±0.74) PSU (−2.01–1.26 PSU) and a corresponding change in SSD of −1.17 (±0.56) kg m−3 (−2.53–0.03 kg m−3) (Supplementary Table 4).Fig. 1: Global map showing the change in sea surface variables across mangrove bioregions under RCP 8.5.a–c, Change in SST (a), SSS (b) and SSD (c). Changes in SST and SSS are based on present-day (2000–2014) and future (2090–2100) marine fields from the Bio-ORACLE database29,30, from which SSD data were derived. The vertical line (19° E) separates the two major mangrove bioregions: the AEP and IWP.Full size imageSpatial variability in predicted surface ocean property changes was examined by considering the two major mangrove bioregions (AEP and IWP) (Fig. 2) and using the Marine Ecoregions of the World (MEOW) biogeographic classification28 (Fig. 3). Both the range and changes in mean SST were comparable for the AEP and IWP mangrove bioregions, for all respective RCP scenarios (Fig. 2a and Supplementary Tables 1–4). Under RCP 8.5, mean SST in both mangrove bioregions is predicted to warm ~2.8 °C by 2100, which is roughly 4.5 times the predicted increase in mean SST under RCP 2.6 (Supplementary Tables 1 and 4). Predictions for the RCP 8.5 scenario are generally consistent with reported global ocean temperature trends33 and show that the greatest warming occurs in coastal waters near the Galapagos Islands (change in mean SST of 3.92 ± 0.06 °C). Pronounced SST increases are also predicted for Hawaii (change in mean SST of 3.36 ± 0.05 °C), the Southeast Australian Shelf (3.30 ± 0.25 °C), Northern and Southern New Zealand (3.25 ± 0.07 °C and 3.34 ± 0.02 °C, respectively), Warm Temperate Northwest Pacific (3.27 ± 0.16 °C), the Red Sea and Gulf of Aden (3.24 ± 0.08 °C), Somali/Arabian Coast (3.23 ± 0.15 °C), South China Sea (3.07 ± 0.10 °C), the Tropical East Pacific (3.09 ± 0.15 °C) and the Warm Temperate Northwest Atlantic (3.14 ± 0.13 °C) (Fig. 3b and Supplementary Tables 4).Fig. 2: Change in surface ocean properties for coastal waters bordering mangrove forests and in the two major mangrove bioregions, the AEP and IWP, for different RCPs.a–c, Variation in SST (a), SSS (b) and SSD (c) under various RCP scenarios. Grey indicates global distribution (n = 10,108), orange denotes AEP (n = 3,190) and green represents IWP (n = 6,918). Data for SST and SSS consist of present-day (2000–2014) and future (2090–2100) marine fields from the Bio-ORACLE database29,30, from which SSD data were derived. The cat-eye plots50 show the distribution of the data. Median and mean values are indicated with black and white circles, respectively, and the vertical lines represent the interquartile range.Full size imageFig. 3: Global spatial variability in SST, SSS and SSD for coastal waters bordering mangrove forests under RCP 8.5.a, Global map showing the provinces (colour code and numbers) from the MEOW database28 used to investigate spatial patterns in mangrove coastal ocean water changes by 2100. b–d, Longitudinal gradient of the change in SST (b), SSS (c) and SSD (d) under RCP 8.5 in the AEP and the IWP mangrove bioregions; circles are coloured according to the MEOW province in which respective mangrove sites are located.Full size imagePredicted SSS changes exhibit an opposite trend in the AEP and IWP bioregions, with increased salinity in the AEP and reduced salinity in the IWP under global warming (RCP 2.6–RCP 8.5; Fig. 2b); this is reflected in contrasting SSD changes in both mangrove bioregions (Fig. 2c) and associated with predicted global changes in precipitation, with extensions of the rainy season over most of the monsoon domains, except for the American monsoon34. Under RCP 8.5, the spatially averaged change in mean SSS is +0.51 (±0.57) PSU in the AEP and −0.68 (±0.44) PSU in the IWP region. The maximum decrease in mean SSS (−2.01 PSU) is predicted for the Gulf of Guinea in the AEP bioregion (Fig. 3c and Supplementary Table 4). Within the IWP, the Western Indian Ocean region shows little or no changes in SSS, which contrasts with the pronounced freshening trends predicted in the eastern part of this ocean basin and the Tropical West Pacific (Figs. 1b and 3c). Increased freshening is predicted in the Bay of Bengal (SSS change: −1.17 ± 0.43 PSU), the Sunda Shelf (SSS change: −1.21 ± 0.29 PSU) and the Western Coral Triangle province (mean SSS change: −0.80 ± 0.17 PSU) (Fig. 3c and Supplementary Table 4). Within the AEP, salinity increases exceed +0.96 PSU in the Tropical Northwestern Atlantic, +0.80 in the Warm Temperate Northwest Atlantic and +0.68 in the West African Transition (Fig. 3c and Supplementary Table 4). The spatial heterogeneity in SSS across the global range of mangrove forests corresponds with observed changes in SSS35. Trends in SSD (Fig. 3d) strongly track changes in SSS (Fig. 3c) rather than SST. All RCP scenarios predict an overall decrease in SSD for both mangrove bioregions; however, the predicted decrease in SSD in the IWP region was a factor of 2 (RCP 6.0) and 2.5 (RCP 2.6, RCP 4.5 and RCP 8.5) stronger than in the AEP (Figs. 2 and 3d and Supplementary Tables 1–4).Propagule density values from our literature survey range from 1,080 kg m−3 for different mangrove species (Fig. 4 and Supplementary Table 5). The low densities reported for Heritiera littoralis propagules provide a strong contrast with the near-seawater propagule densities reported for Avicennia and members of the Rhizophoraceae (Bruguiera, Rhizophora and Ceriops). Floating characteristics of the latter may be particularly sensitive to changes in SSD. To illustrate the potential influence of changing ocean conditions on mangrove propagule dispersal, we considered threshold water density values (1,020 and 1,022 kg m−3) that are within the range where elongated propagules of important mangrove genera tend to change floating orientation (Fig. 4a). More specifically, we determined the ocean surface area with an SSD below or equal to these thresholds under different climate change scenarios (Fig. 5). Under RCP 8.5, the ocean surface covered by mangrove coastal waters (coastal waters bordering present mangrove forests) with a density ≤1,020 kg m−3 increases ~27% by 2100, notably more so in the IWP (~37%) than in the AEP (~6%) (Supplementary Table 6). A threshold of 1,022 kg m−3 results in increases of roughly +11% (global), +12% (IWP) and +8% (AEP) (Supplementary Table 7). Similar spatial patterns are observed for open-ocean waters within the global latitudinal range of mangroves (Fig. 5 and Supplementary Figs. 1 and 2).Fig. 4: Potential effect of future declines in SSD on mangrove propagule dispersal.a, Range of reported propagule density values for wide-ranging mangrove species and present and future range of SSD for coastal waters along the range of those mangrove species. Mangrove propagule data are extracted from the literature (Supplementary Table 5). H. lit, Heritiera littoralis; X. gra, Xylocarpus granatum; A. ger, Avicennia germinans; A. mar, Avicennia marina; B. gym, Bruguiera gymnorrhiza; C. tag, Ceriops tagal; R. man, Rhizophora mangle; R. muc, Rhizophora mucronata. Bottom part adapted from ref. 51. b, Conceptual figure of the potential effects of ocean warming and freshening on mangrove propagule dispersal. Ocean warming and freshening drive changes in SSD and may reduce the timeframe for opportunistic colonization. For a propagule with a specific density and floating profile under present surface ocean conditions, reduced SSD of coastal and open-ocean waters may reduce floatation time (shaded area) and hence, reduce the proportion of long-distance dispersers. For simplicity, the density of propagules is assumed to increase linearly over time, although the actual increase may be nonlinear.Full size imageFig. 5: Future changes in SSD.a–d, Spatial extent of coastal and open-ocean surface waters with a density ≤1,020 kg m−3 (a,b) and 1,022 kg m−3 (c,d), for present (2000–2014) (a,c) and future (2090–2100; RCP 8.5) (b,d) scenarios. Data are shown for surface ocean waters within the global latitudinal range of mangrove forests (between 32° N and 38° S). The two density thresholds considered are within the range of densities at which mangrove propagule buoyancy and floating orientation of several mangrove genera change, as reported in available literature. Black dots along the coast represent the global mangrove extent from the 2015 GMW dataset32. Magenta-coloured circles represent SSD values More