Abstract
Exposure to coastal flooding is increasing due to growing population and economic activity. These developments go hand-in-hand with a loss and deterioration of ecosystems. Ironically, these ecosystems can play a buffering role in reducing flood hazard. The ability of ecosystems to contribute to reducing coastal flooding has been emphasized in multiple studies. However, the role of ecosystems in hybrid coastal protection (i.e. a combination of ecosystems and levees) has been poorly quantified at a global scale. Here, we evaluate the use of coastal vegetation, mangroves, and marshes fronting levees to reduce global coastal protection costs, by accounting for wave-vegetation interaction.The research is carried out by combining earth observation data and hydrodynamic modelling. We show that incooperating vegetation in hybrid coastal protection results in more sustainable and financially attractive coastal protection strategies. If vegetated foreshore levee systems were established along populated coastlines susceptible to flooding, the required levee crest height could be considerably reduced. This would result in a reduction of 320 (range: 107-961) billion USD2005 Power Purchasing Parity (PPP) in investments, of which 67.5 (range: 22.5- 202) billion USD2005 PPP in urban areas for a 1 in 100-year flood protection level.
Similar content being viewed by others
Introduction
Globally, about 600 million people are at risk of coastal flooding, of which 320âmillion are in urban areas1. Future population growth and urbanization will expose an increasing number of people and amount of assets to coastal flooding2,3. Originally, coastal areas are solely protected by natural features. However, coastal areas where these features alone are insufficient, where occupied land is low lying or where people have encroached coastwards are often protected from floods using human-made structures. These structures are also known as âgreyâ coastal protection, such as seawalls and levees. Nowadays Nature-based Solutions (NbS) are considered potentially sustainable and cost-effective complements to engineered flood defences4,5,6. For example, foreshores vegetated with salt marshes and mangroves can substantially reduce incoming wave heights7,8,9,10, meaning that levees protecting the hinterland can be lower than âgreyâ coastal protection, resulting in reduced initial investment costs and maintenance costs. Besides flood hazard reduction during storm events, marshes and mangroves provide many other ecosystem services such as: carbon storage, habitats for fish and birds, improved water quality and accumulation of sediments11,12. These wetland areas are decreasing globally13,14. Coastlines are inherently dynamic and are either accreting, eroding or stable15,16. In the short term (coming decade), the structural loss of wetlands is likely to result in extra investment costs for coastal infrastructure, and in the long term (this century) it may lead to widespread loss of coastal lands and relocation of millions of people inhabiting coastal areas.
As flood risk is expected to increase in the future as a result of rising extreme sea levels and socioeconomic developments, there is a need to increase coastal resilience accompanied with a strong demand for coastal flood protection measures17,18. The role that coastal vegetation can play in reducing coastal flood risk has been quantified in several local, regional, national and global studies5,19,20,21,22. Only few global studies have been based on process-based wave modelling, and those that have taken this approach focused on mangroves only23. While these studies exemplify the role of coastal vegetation as ecosystem services to reduced flood risk, the potential of combining coastal vegetation with traditional flood defence measures so-called hybrid or green/grey protection)âhas not been assessed on a global scale. Such an assessment is important because hybrid protection can be very effective. Small- to medium-sized coastal vegetation belts alone cannot prevent inundation21, while levees without protective coastal vegetation in front require larger dimensions, and hence investments, for the same level of safety compared to hybrid protection. Thus conserving coastal ecosystems and accounting for coastal ecosystem presence fronting âgreyâ coastal protection could be cost-effective. Considering hybrid protection on a global scale is unique as it bridges the gap between studies on ecosystem services and flood risk reduction.
A vegetated foreshore in front of a levee reduces wave height and thereby wave run-up and wave overtopping10,24,25. Consequently, integration of vegetated foreshores in coastal protection systems, as a supplement to seawalls and levees, allows for lower crest heights24 and results in hybrid coastal protection systems that are more adaptive to new information and conditions that may emerge, such as sea-level rise (SLR)25. In this study we focus on wave-vegetation interaction and as a conservative approach we neglect the effect of coastal vegetation on storm surge levels, because surge reduction is typically only relevant for very extensive coastal vegetation belts26 and is largely dependent on the local spatial configuration27. In addition, coastal ecosystems have the ability to build up vertically due to the accumulation of sediments caused by biophysical feedback mechanisms28. Hereby, coastal ecosystems can alter the intertidal elevation profile. This study focusses on the present situation, but (future) changes in intertidal elevation (e.g. due to reduced sediment supply or SLR) will influence propagation of both waves and storm surges29.
To assess the effect of coastal vegetation on reducing levee crest height, we assumed the presence of a levee at the back of the vegetated foreshore. Global information on the presence of coastal levees is not available. In addition to available global elevation and bathymetric maps more accurate intertidal elevation data (20âm horizontal resolution, 0.52âm RMSE vertical accuracy) are obtained by creating a new global data layer based on time-ensemble average satellite images of the probability of inundation (Methods). Similarly, to acquire high-resolution data on the presence of coastal vegetation, Sentinel-2âA and Landsat-8 images were used to construct a global vegetation map using NDVI values (Methods). This map was combined with existing vegetation maps14,30,31,32 to add vegetation type, focusing on marshes and mangroves. As forcing, we used wave heights and periods from a re-analysis of ERA-Interim33 and extreme water levels from a global tide and surge model34.
Coast-normal transects were defined for global coastlines ranging between 66° N and â60° S, with alongshore distances of ~1âkm. For each transect, a bathymetric profile, vegetation cover and hydrodynamic boundary conditions were determined based on the described data sources. The boundary conditions were derived for nine flood level return periods (2, 5, 10, 25, 50, 100, 250, 500 and 1000 years). We translated offshore wave heights to nearshore conditions and calculated foreshore wave propagation (Methods). Wave damping by vegetation was obtained by comparing wave propagation over the transect with and without vegetation. For both situations, we determined the required levee crest height to prevent flooding, assuming the same landform and a levee at the back of the vegetated foreshore (Supplementary Fig. 1). The difference in crest height between the two situations was used to monetize the effect of coastal vegetation, using unit investment costs of levees that are corrected for differences in construction costs across countries (Methods). We identified populated coastal areas susceptible to flooding based on inundation using flood maps of 1âkm resolution and different population density classes (Methods). Finally, a distinction is made between urban and rural areas (Methods).
Here, we show the first global overview of the protective value of hybrid coastal protection, in which coastal vegetation is integrated in vegetated foreshore-levee systems. We assess the potential reduction in costs compare to traditional âgreyâ coastal protection. The resulting cost saving represents a reduction in coastal protection investment costs. This also applies to areas where levees already exist, as including the effects of vegetation may imply that costs for future levee heightening, for example due to projected sea-level rise, are not necessary. For areas where no levees currently exist, this method gives an approximation of the costs that could potentially be saved. In the current study we use open source earth observation (EO) data and tidal statistics to produce unprecedented high-resolution global intertidal elevation maps (20âm horizontal resolution, 0.52âm RMSE vertical accuracy) and vegetation maps (10âm resolution) (Methods) to overcome data scarcity in the intertidal zone. In combination with other open source global datasets14,30,31,32,35,36, we use these data in a numerical model of wave attenuation to obtain the first process-based global assessment of flood hazard reduction by coastal vegetation in hybrid coastal protection systems. This assessment offers valuable insights in locations where coastal vegetation is of great importance and where hybrid coastal protection systems could be applicable.
Results
Coastal vegetation presence
We find that 18.5% of the global coastline is covered by marsh or mangrove vegetation. This is mostly a fairly narrow (25â250âm) vegetated belt, which covers 6.3% of the global coastline. In addition, 3.9% and 3.3% is covered, respectively, with a vegetation belt of 250â500âm and 500â1000âm. Furthermore, a vegetation belt exceeding 1000âm width covers 5% of the global coastline (Fig. 1a). Wide vegetation belts occur near the equator, where extensive mangrove forests still exist. Of the worldâs coastline, 17% is populated (population density greater than 1 per km2) and susceptible to coastal flooding (hereafter referred to as the susceptible coastlines). Of the rural and urban populated coastlines susceptible to flooding, 25% and 18.5% are vegetated, respectively (Supplementary Fig. 3). Wider vegetation belts generally occur along rural coastlines. On average, the width of mangrove belts in populated and flood susceptible locations is 858âm in urban areas and 1233âm in rural areas. For marshes, the average width of the vegetation belt in populated and flood susceptible locations is 483âm (urban) and 613âm (rural) (Supplementary Fig. 3).
a Global distribution of coastal vegetation belt width in metres. b Reduced wave transmission solely due to wave-vegetation interaction based on storm conditions with a return period of 100 years along global coastlines susceptible to flooding. Map is created with Python 3.8.10 (https://python.org) using Cartopy (v0.18.0. Met Office UK. https://pypi.python.org/pypi/Cartopy/0.18.0), GeoPandas v0.8.1 (https://geopandas.org) and Matplotlib v3.3.497.
Wave attenuation by coastal vegetation
We find that for 11.5% of the global coastline, vegetated foreshores can reduce incoming wave heights by more than 25%, solely due to wave-vegetation interaction (Fig. 1b). This figure is calculated by comparing wave heights for a return period of 100 years between bare and vegetated foreshores. Despite the large mean width of the vegetation belt in rural areas, the difference in wave reduction between rural and urban areas is limited. This is mainly due to the non-linear relationship between wave attenuation and vegetation width37 (Supplementary Fig. 3).
Impact on required levee crest heights
For 27.6% of populated susceptible coastlines, the current presence of coastal vegetation allows for lower levee crest heights while maintaining the same protection standard. The required crest heights for a 100-year protection standard are reduced by 25â50âcm for 5.8% of the susceptible coastline, 50â75âcm for 9.5% and by more than 75âcm for 12.3%. Across all vegetated susceptible coastlines, the mean crest height reduction for a return period of 100 years is 96âcm. For return periods of 2 and 1000 years, the mean crest height reductions are 82 and 104âcm, respectively (Supplementary Fig. 5a). For 22.1% of the susceptible coastlines, the current presence of coastal vegetation allows for a reduction in the required crest height equal to or greater than the projected sea-level rise of 0.49âm (RCP4.5) by the end of the 21st century38.
Results on country level
If levees with a 100-year protection standard were constructed along all populated susceptible coastlines, the potential reduction in costs resulting from the presence of current foreshore vegetation presence is 320.2 billion USD (all monetary values in USD2005 PPP), of which 67.5 billion USD is in urban areas (Fig. 2a). The ten countries with the highest reductions in cost are shown in Fig. 2a. These ten countries all have relatively long coastlines and account for over 65% of the total potential reduction in costs, with a combined potential of 208.6 billion USD. The average levee crest height reduction per kilometre for susceptible coastlines is 0.24âmâkmâ1 globally, with the highest average value in Guinea-Bissau (0.77âmâkmâ1) (Supplementary Fig. 4). Expressing the potential saved investment costs as a percentage of the countryâs gross domestic product (GDP) highlights the importance of current coastal vegetation for Small Island Developing States (SIDS) (Fig. 2b). For many of these countries, âgreyâ coastal protection along their full coastline is not socially and financially feasible39. Hence, coastal vegetation as part of hybrid coastal protection offers opportunities at these places.
a Potential levee cost reduction by coastal vegetation per country in billion USD2005 PPP. b Potential levee cost reduction by coastal vegetation in urban areas per country as percentage of GDP2005. Values in these maps apply to a 100-year protection standard and a âMediumâ levee unit cost scenario. Non-benefiting countries are indicated with light grey colour. Map is created with Python 3.8.10 (https://python.org) using Cartopy (v0.18.0. Met Office UK. https://pypi.python.org/pypi/Cartopy/0.18.0), GeoPandas v0.8.1 (https://geopandas.org), Matplotlib v3.3.497 and GADM v2 (https://gadm.org) administrative boundaries.
The majority of the coastlines along which we find potential to implement hybrid coastal protection have a low to medium exposed population (1â100 people kmâ2) (Fig. 3). Eight out of the top ten countries with the potential to reduce coastal protection costs by applying hybrid coastal protection as Nature-based Solutions in higher populated areas (>100 people km-2) can be found in Asia (Philippines, China, Indonesia, India, Vietnam, Myanmar), Africa (Madagascar, Nigeria, Guinea) and Europe (United Kingdom).
Salt marsh (a) and mangrove (b) width of vegetation belt for four classes of population density. Wave attenuation on foreshores vegetated by salt marshes (c) and mangroves (d) in comparison to bare foreshores. The bars represent the upper 85%-ile and the 15%-ile. Levee cost reduction for three levee unit cost scenarios, three standards of protection and four classes of population density for both salt marshes (e) and mangroves (f).
On average the width of the vegetation belt decreases with increasing population density (Fig. 3), although this trend is not statistically significant given the large error margins. This can be related to the fact that vegetation extent has been decreasing for decades as a result of conversion of green belts for purposes such as agriculture or urban land use40. As a result, the potential wave height reduction by mangroves and salt marshes decreases for vegetation belts with lower widths. The hydrodynamic interaction between waves and vegetation depends on physical parameters such as the bathymetry profile, surge level, wave height, vegetation width and height. While these parameters are important for the correct calculation of wave height reduction, uncertainty in the levee cost reduction results are mainly determined by uncertainty in levee construction costs (Supplementary Fig. 10). To account for this uncertainty we applied three different standard levee unit cost scenarios (low, mid, high) in line with previous research41,42. The difference between the high and low levee cost scenarios is far greater than the difference between a low standard of protection (RP5) and a high standard of protection (RP1000) (Fig. 3).
Discussion
This global study highlights where including coastal vegetation in hybrid coastal protection schemes is applicable and provides insight in the total potential cost reduction of these Nature-based Solutions. Existing coastal vegetation can contribute substantially to cost-effective flood risk reduction strategies for countries along a substantial part of the global coastlines. The ability of coastal vegetation to reduce wave impact is vital for countries with long rural coastlines where conventional protection by hard infrastructure alone might be economically unfeasible (Fig. 3e). In addition, several small island developing states receive major benefits from the wave-reducing abilities of coastal vegetation (Fig. 2b). Hybrid coastal protection schemes combining coastal vegetation and levees or seawalls, have large potential in countries such as China, United Kingdom and Indonesia, where coastal vegetation helps to protect the highest number of people and amount of assets. The use of hard flood risk protection structures is widespread around the globe. With SLR and growing population in vulnerable areas, future economic benefits and the need for protective infrastructure are expected to be larger17,28,43,44. Results of the current study should be considered as a warning not to remove coastal vegetation and to regulate coastal developments strongly, both along rural and urban coastlines. Mangroves and salt marshes play an important role in keeping these coastlines stable and could considerably lower the costs of coastal protection infrastructure. Moreover, next to the benefits that coastal vegetation can provide to coastal flood protection, other benefits of coastal vegetation include positive impacts of tourism45, regulating fish and shellfish stocks and playing a role in the sequestration of carbon12.
Global mapping studies that combine models and different data sources inevitably accumulate uncertainties. For example, the limited spatial resolution affects the accuracy of the metrological forcing data used from ERA-I, the Global Tide and Surge Model (GTSM), and the resulting estimates of extreme waves and water levels. ERA-I waves have a bias between 0 and -0.2 metres, which results in an expected error (scatter index) between 16 and 20%33. Extreme water levels show an average relative error in the range of 11â14%34. Further, the newly developed intertidal elevation data constructed using EO data and tidal epoch data has a vertical accuracy of 0.52âm RMSE (Supplementary Fig. 7). Note that to be conservative we did assume the same landform for the simulations with and without vegetation. For this study, each data layer was validated with local data or local models, and a qualitative reliability analysis is performed (Supplementary Figs. 6, 7, 8, 9). We performed a sensitivity analysis (Methods) to identify the factors that have the largest influence on the coastal reduction costs (Supplementary Fig. 10). Despite the large uncertainty in the hydrodynamic data and topographic data used, our analysis shows that the levee costs reduction potential is influenced most by assumptions on the required coastal protection needs (people exposed) and associated levee construction costs (Fig. 3, Supplementary Fig. 10). This implies that the current study can be improved by further limiting the uncertainty in the levee costs and requirements, activities that are typically undertaken in a local assessment. Deviations in the width of the vegetation belt play a minor role in comparison to the aforementioned variables, because of the non-linear relationship between wave height reduction and vegetation width37. Global assessments have their limitations, and mismatches between model abilities and user expectations should be prevented46. The current study provides a global overview and is meant to stimulate detailed local assessments in areas where it matters most. We stress that local assessments are required for the correct design and implementation of Nature-based Solutions.
The percentage of the susceptible populated coastline benefitting differs for the various protection standards (Supplementary Fig 5), with 19.2% and 23.2% for a 2 and 1000 years protection standard respectively using a critical wave overtopping rate of 1âlâsâ1âmâ1. Here, benefitting is defined as having a reduction in required levee crest height of at least 50âcm due the presence of coastal vegetation. This spread (19.2â23.2%) is smaller in comparison to the outcomes for various critical wave overtopping rates. For rates of 0.1, 1.0 and 10âlâsâ1âmâ1 the percentage of coastlines benefitting is 25.6%, 21.8% and 11.1%, respectively, for a 1 in 100-year protection standard. These results point out that the reduction of the required crest height is largest for simple earthen levees, as the applicable critical wave overtopping rate depends mainly on the quality of the levee47. Consequently, reduced levee heights in countries with more strict levee quality are smaller, as the reduction in required levee height due to the presence of coastal vegetation is less. However, the associated cost savings of high quality and low quality levees might be comparable. A smaller levee crest height reduction for high quality levees might be balanced out, because high quality levees may actually be more expensive.
For this study the effects of vegetation on storm surges were not included. Although several studies mention the positive effects of marshes on surge reduction24,29, other studies emphasize that surge reduction depends strongly on storm duration and intensity, and that effects are only significant for mild and short storms48 and that surge levels of severe storms are comparable for bare tidal flats versus marshes with widths below 1200âm6. This suggests that surge propagation is context specific and influenced by local features, specificially nearshore geomorphological configuration. Conservatively, we assumed that to reduce surges, vegetated foreshores of multiple kilometres are required. Such large wetlands only occur on a few deltas around the world, with only 4 out of 11 providing effective protection26. Most of the vegetation belts along populated susceptible coastlines have an insufficient width (median belt width of 525âm) to have substantial effects on surge reduction. Vegetation belts with widths exceeding 2âkm only exist for 3.0% of the vegetated transects in urban areas susceptible to flooding.
The resilience of coastal vegetation under severe conditions is uncertain49. Damage reports following severe events, such as hurricane Haiyan, show destruction of planted mangroves that were in the path of the typhoon50. Uprooting and shear strength influence the ability of coastal vegetation to reduce wave heights during storm conditions, but are not captured by the current model. For this study we determined coastal vegetation belt widths, but we used spatially constant (conservative) vegetation characteristics (Methods). For salt marshes, we used a winter state as found in NW Europe, and for mangroves we used characteristics of young fringing pioneer mangroves (Methods). In reality salt marshes and mangroves are complex habitats with considerable spatial variation due to creeks, differing vegetation cover and sediment characteristics. Performing the current analysis with spatially varying vegetation characteristics (such as mangrove canopy height51 and root density), based on new EO data techniques, would be an interesting topic for future study. In the current study, we tried to make a global assessment while approaching reality as much as possible by focusing on the main factors of relevance and reducing uncertainty by developing improved global data layers and by validating results. Field evidence to support the role of wave attenuation by vegetation is extremely limited for more extreme waves and surges. Therefore, implementation of formulations and results of field studies during non-storm conditions, scaled lab studies and results of large-scale salt marsh flume tests10 in (numerical) models are to date the best way to gain insights on the buffering role of coastal vegetation for more extreme conditions. Future large-scale mangrove flume experiments may provide a better evidence base and improved understanding of the role of mangroves on wave attenuation during extreme storms.
Our study focused on storm conditions, which have a low probability of occurrence and last for a short period. Sea-level rise (SLR), on the other hand, is a slower and long-term process with far-reaching consequences. The effect of SLR on intertidal areas depends on complex hydrodynamic interactions between tides, waves, wind, fresh water run-off (for estuarine wetlands) and sediment supply. These complex interactions hamper the long-term modelling and prediction of the effects of accelerated SLR52. On intertidal mudflats, wave energy is re-distributed and dissipated, which limits wave action on coastal vegetation present at the landward limits of the intertidal zone. The decay of intertidal flats will therefore put more stress on existing salt marshes and mangroves53. In addition, mechanisms such as coastal squeeze will determine the future state of coastal ecosystems and their coastal vegetation belts locally54. For vegetated areas the situation is even more complex, as also root growth, compaction and subsurface faunal processes play a role55. Despite these difficulties, multiple studies28,56,57,58 show the ability of vegetated foreshores to accumulate sediments and thereby keep up with SLR, but also underline the presence of tipping points. Studies indicate that once the rate of SLR exceeds a context dependent critical threshold, vegetated foreshores cannot keep up with SLR resulting in die-off of the lower vegetated part first57,59,60. Assessing the potential loss and degradation of coastal ecosystems due to the effects of SLR and coastal squeeze and the resulting effect of the partial loss of vegetationâs wave buffering capacity on coastal protection costs would be an interesting field for future study.
Global analyses involve inevitable assumptions and shortcuts. In this study a replacement costs method was used to express the benefits of applying hybrid coastal protection over purely âgreyâ coastal protection, hereby assuming that levees would be built along global coastlines with coastal protection needs. However, the answer to the question to protect or not to protect is not easily answered. For example, the protection level against flooding differs per country and can be supported by an analysis of investment costs versus avoided damages61. The local evaluation of different flood risk reduction strategies should be done case-by-case. Flood risk could depend on local deviations in topography that are not captured in this assessment that uses global data and 1âkm spaced coast-normal transects. In addition, we considered solely locations where vegetation is currently present. Rehabilitation attempts of vegetated foreshores should consider eco-morphodynamic requirements62 for a more sustainable use of ecosystems, as vegetated foreshores are part of the ecosystem and cannot be considered as isolated structures. Also, levee location is assumed directly behind the vegetated foreshores, but the exact levee position is best determined in a local assessment taking into account factors including the presence of infrastructure (e.g. houses and roads), subsidence63 and the resilience of the vegetated foreshores64 (including for example wave reflection from the levee that can hinder geomorphological development65). With SLR the latter becomes more important, as levee construction can limit landward migration of vegetation in the long run66.
Considering the current climate and biodiversity crisis, integration of coastal ecosystems in coastal protection schemes is essential. Currently, levees and seawalls are also regularly found along coastlines with little exposure or hazard (Fig. 4). In these cases, their negative impact by disturbing natural flows of water and sediment and thereby reducing coastal resilience may not outweigh their positive impacts in terms of protecting people or infrastructure. In addition, they are often constructed seaward of coastal vegetation, which means that (a) the survival of vegetation itself is threatened; and (b) the wave-reducing potential of the vegetation is not used. Moreover, to limit future coastal ecosystem loss and to work towards cost-effective and sustainable coastal protection schemes, hard infrastructure should be used with care and vegetated foreshores should be formally incorporated in coastal protection schemes.
Our study shows that integrating present coastal vegetation into infrastructural interventions could result in considerable cost savings. Ongoing removal of coastal vegetation will increase wave forces and is thereby likely to increase erosion, which may result in considerable land loss11,60,67. Allowing the removal of coastal vegetation for short-term individual or industrial profits might result in large community costs through resulting land loss and by raising the costs for future coastal protection infrastructure. Formal integration of vegetated foreshores with infrastructure design and sustainable management of coastal ecosystems is a cost-effective and sustainable route towards maintaining many coastlines around the world.
Methods
Coastline segments
For reasons of data availability and socioeconomic relevance, the analysis was limited to latitudes between 66° N and â60° S. In this area of interest, the world was divided in 1 arcmin (~2âkm) grid cells. To define a logical position for the establishment of an efficient levee, the coastline location was derived from the OpenStreetMap68, moved 100âm land inward and smoothed. For every cell containing a coastline segment, coastline length and a coast-normal transect were derived at the center of segments resulting in 495.361 transects that are on average 1.1âkm apart. Bootstrapping revealed that transect distances up to 2âkm give very similar results. All transects stretch 4âkm seaward and 4âkm inland to fully capture most foreshores.
Elevation data
A global intertidal bathymetry/elevation dataset from high-resolution EO data (USGS Landsat and Copernicus Sentinel-2), the Foreshore Assessment using Space Technology (FAST) intertidal elevation map69, was produced to compliment commonly used global data products with low resolution and higher inaccuracy in intertidal zones. Global coastlines were divided over 25000 tiles of each 40âÃâ40âkm2. For these tiles, all available images were collected for the period between 1997 and 2017. Surface water was identified, using normalized difference spectral indices (NDSI, here SWIR1 and Green band) for all images (median of 317 images per tile) covering various tidal conditions, and the per pixel mean calculated to derive time-ensemble average (TEA) NDSI images. We developed a new technique to transform TEA images to intertidal elevation independently of in situ calibration data. TEA-NDSI images were normalized by the spatially averaged NDSI values of regions identified (using global elevation datasets) as land and water, respectively. This resulted in a single image per tile that represented the inundation probability for each pixel in the intertidal zone. The inundation probability represents the long-term average tidal inundation, because it was derived from a collection of images that span a time period similar to the tidal epoch (period of 19 years). Pixels having a probability of 1 represent permanent water, and have elevations less than or equal to the lowest astronomical tide (LAT), whereas land (pâ=â0) represents elevations higher than or equal to the highest astronomical tide (HAT). By deduction, pâ=â0.5 is equivalent to local mean sea level (LMSL). Tidal statistics from the global tide model FES2012 were used to couple the derived inundation probability to an elevation. The main source of bed level data originates from this map and has a 20âm horizontal resolution and typically a 30â50âcm vertical accuracy (RMSEâ=â0.52âm, MAE 0.42âm, as assessed at a number of sites with high quality elevation data (Supplementary Fig. 7)). Bathymetry data (GEBCO35; 30 arc-second horizontally, tens of metres vertically) and topography data (MERIT36; 3 arc-seconds, 2âm vertically) were merged to create a continuous bathymetry-elevation map by changing the vertical datum of MERIT from EGM96 to MSL by assuming 0âm +MSL at the OSM coastline. Global bathymetry datasets (e.g. GEBCO) and elevation datasets (e.g. SRTM and MERIT) lack accuracy (especially nearshore), but are commonly used17,18,23,34. The final bed level was constructed using FAST intertidal data where sufficient valid data points were available, complemented by the merged GEBCO-MERIT data where these points were lacking.
Vegetation extent
The FAST coastal vegetation map69 was based on Landsat-8 and Sentinel-2 satellite images collected between 2013 and 2017. The map provides actual vegetation presence at 10âm resolution. Vegetation presence was obtained by applying an individual NDVI threshold per tile, with a total of 25,000 tiles, based on the yearly NDVI average and NDVI amplitude. The FAST coastal vegetation map is validated based on NDVI comparison with local measurements taken at Zuidgors, The Netherlands (R2â=â0.92) (Supplementary Fig. 8). If vegetation was present, the vegetation type was determined by global salt marsh32 and mangrove14 maps, complemented with Corine Land Cover30 (CLC, Europe only) and GlobCover v2.231 maps when there is no coverage. Determining global coastal vegetation extent is difficult and affected by eutrophication in coastal environments. This behaviour is observed on the coast along the Persian Gulf and the Red Sea. To improve accuracy only vegetated transects identified by the global salt marsh32 and mangrove14 map and confirmed by the FAST coastal vegetation map are included for these areas. Moreover, vegetated transects with a green belt width smaller than 250âm identified by GlobCover are excluded from the study for accuracy reasons (Supplementary Fig. 8). To avoid mixed vegetation types, the vegetation type was determined by the most dominant type. The vegetation width constituted of the sum of vegetated grid cells between the start and the end of the vegetated zone.
Water level and wave data
The design water levels were based on a combination of tide and storm surge for the selected probability of occurrence (return periods 2, 5, 10, 25, 50, 100 default, 250, 500, 1000 years) and came from the GTSR dataset34. SLR and subsidence were not taken into account because this study focuses on the present situation. Moreover, quantifying the future role of vegetated foreshores would not only require SLR scenarios but also an insight in the development of wetlands over time, which is strongly determined by local conditions such as sediment supply56,57,60. Offshore wave conditions were obtained from ERA-Interim33 re-analysis, based on data from 1979 till 2017 and reprojected to Dynamic Interactive Vulnerability Assessment (DIVA)70 points. Next, the Peak Over Threshold method was applied to construct representative values for the significant offshore wave height, Hs and the peak wave period Tp for all the return periods. The nearshore wave height was limited by the local water depth at the start of the (vegetated) foreshore using a breaker criterion (gammaâ=â0.55). This is a fairly low value considering the range of values cited in literature71 leading to conservative wave attenuation by vegetation results. Wave-bottom interactions in the sub-tidal zone and processes such as refraction and diffraction are not explicitly simulated. The conservative breaker criterion is chosen to implicitly account for these processes in a conservative manner. The wave period remained unchanged and the wave direction was assumed coast normal and wave growth along the transect due to wind effects was excluded. However, for the current study a more sophisticated approach to account for longshore wave variability based on topography was considered infeasible at the global scale and considered to yield limited outcome looking at the uncertainty in socioeconomic factors. The average Hs,offshoreâ=â4.6âm (stdâ=â2.0âm) and the average Hs,startforeshoreâ=â0.7âm (stdâ=â0.7âm).
Profile construction
The 8 kilometre coast-normal transects consisted of 321 gridpoints, thus a horizontal grid resolution of 25âm. We used four different methods: Foreshore method 1 (based on the FAST intertidal elevation map), Foreshore method 2â4 (based on MERIT-GEBCO). The properties of the FAST intertidal elevation map, MERIT and GEBCO are described under the header âElevation dataâ. Foreshore method 1 produced the most accurate profiles and foreshore method 4 the least accurate profiles. The profile construction steps are described hereafter. Validity checks were performed to identify false indications of intertidal area in the FAST intertidal elevation map. Individual data points were marked invalid and removed in case: (1) MERIT points were situated above the surge level with a return period of 2 years, while data from the intertidal map indicated a lower elevation. (2) Data from the FAST intertidal map was situated at open sea. (3) Data from the FAST intertidal map along the transect dropped below a minimum range threshold of 10âcm. A fourth check was performed based on the continuity of the data. Data from the FAST intertidal map contain discontinuities along the profile. These continuities exist on pixel level due to the use of the modified normalized difference water index and in some instances cloud coverage was preventing full coverage. Lastly, discontinuities arise due to the presence of (high elevated) tidal flats and banks in coastal areas. (4) Data length was defined as the length of continuous data points along the transect. If the data length of a patch decreased below a threshold of 100âm, the points were marked invalid. Gaps between valid data patches were filled using linear interpolation if the gap was smaller than 250âm. Eventually, one, none or multiple valid data patches were found along the transects. See Supplementary Fig. 2 for example transects.
Global coastline shapes range from straight sandy coastal stretches to complex coastlines often found in estuaries. With a transect length of 8âkm, the start and the end of the transects could both be situated on land, hampering an unambiguous identification of the foreshore of interest. We designed the algorithm such that the last foreshore was selected. For profiles using data from the FAST intertidal map (foreshore method 1, 50.9% of populated susceptible coastlines), the last valid patch corresponds to the last foreshore. The inclusion of tidal flats as part of the foreshore was determined based on the gap length. In case no (sufficient, thus not satisfying the minimum data length criterion of 100âm) valid data was available from the FAST intertidal map based on the four described checks, the profile was based on a merged GEBCO-MERIT set (methods 2, 3 and 4), respectively, 46.1%, 3.0% and 0.01%. For the second method, data points were selected between a minimum threshold of â2âm MSL and a maximum threshold equal to the surge level with a return period of 2 years. Next, for the selected points the direction of the slope was determined by comparing elevation between the data point concerned and the next data point. This resulted in patches of upward sloping sets of data points between the minimum and maximum threshold. Similar to foreshore method 1, the validity of the patches was checked using data length, gap length and the corresponding thresholds of 100âm and 250âm. The start and the end of the foreshore were determined by the first and last valid point of the last patch. Foreshore method 3 was used if not sufficient foreshore data were available to satisfy the minimum data length threshold (100âm). In these cases, the start of the foreshore was defined as the first upcrossing intersection with â2âm MSL along the transect. The end of the foreshore corresponded to the intersection between the elevation profile and the governing surge level with a return period of 2 years. Foreshore method 4 was used if no start and or end of the foreshore could be found. In this case the start and/or end point of the foreshore corresponded to the first and last data point, respectively.
In some cases, elevation for the end of the foreshore was missing due to several reasons. First, the upper part of the intertidal zone was sometimes missing from the FAST intertidal map, due to low frequency of inundation of the upper intertidal zone or cloud cover. Second, bed elevation in mangrove belts was hard to define based on satellite imagery, as the canopy is detected as the earth surface. These uncertainties were counteracted by consulting the mangrove and salt marsh maps. If vegetation was present in one of these maps, the derived foreshore was extended until the end of the vegetated zone. An elevation equal to the surge level with a return period of 2 years was chosen as elevation for extended foreshore points with an elevation exceeding this surge level.
Vegetation parameters
As deducting the type and size of mangrove trees and salt marshes from EO data at global scale is not possible (yet), the current modelling approach relies on field and literature observations. For the scope of this research the properties of the mangrove trees occurring at the seaward side of the mangrove belt are the most relevant. To avoid overestimation of wave attenuation in young mangrove forests, the mangrove dimensions are chosen such to be representative for young fringing pioneering mangroves up to a height of 3âm that are practically vertically uniform compared to mature trees. The modelling approach uses four parameters to represent vegetation: height, diameter, number of stems and drag coefficient. The exact characteristics are based on observations in literature8,9,72,73,74,75,76 (Nâ=â30âmâ2, dâ=â35âmm, hâ=â3.0âm).
High quality observations on wave attenuation by mangroves under storm conditions do not exist. For the drag coefficient the theoretical value, 1, of a rigid cylinder is chosen, because mangrove trunks can be considered rigid. For salt marshes a winter state representative as found in NW Europe is chosen. The values are defined based on FAST field tests (Romania, UK, Spain and the Netherlands) and literature10,24,77,78 (Nâ=â1225âmâ2, dâ=â1.25âmm, hâ=â0.30âm). A drag coefficient (CD) of 0.19 is chosen, which is the lower limit found during large-scale flume tests10. The drag coefficient depends on biophysical characters as well hydrodynamics. The drag coefficient represents drag due to skin friction and pressure differences, but also effects like swaying motion of stems24. The 1D modelling approach takes into account gaps in vegetation cover, e.g. due to the presence of channels. Zonation of vegetation types is not implemented, because this level of detail is insignificant in relation to the inaccuracies induced by the use of global datasets.
Wave attenuation model
To determine wave attenuation along the foreshore transects and the resulting significant wave heights relevant for the flood defence on a transect, we used a lookup-table approach. The lookup table was generated by combining 668,304 model output values for different combinations of foreshore slopes, vegetation covers and hydrodynamic conditions. The table contained wave heights modelled by XBeach79 in surfbeat mode (a nearshore numerical wave model that accounts for the presence of vegetation) at regular intervals along a steady slope, both with and without vegetation. XBeach uses for wave-vegetation interaction the rigid cylinder80 approach and includes an energy sink term to the wave energy balance to implement wave dampening81. We used conservative vegetation characteristics, winter state salt marshes and young pioneering mangroves. We characterized foreshores by their width and slope. The foreshore profile was the same for simulations with and without vegetation. The foreshore width was determined by calculating the distance between the start and the end of the foreshore. The slope was estimated using a linear regression. This approach has two advantages over detailed modelling of wave attenuation over all transects: it is much quicker, allowing for iterative improvements of the workflow and it does not suggest the precision one would expect from detailed models but cannot be delivered with global data. Average Hs,endforeshore,novegâ=â0.6âm (stdâ=â0.5âm) and Hs, endforeshore,vegâ=â0.3âm (stdâ=â0.4âm).
Coastline susceptible to flooding, urban and rural extents and population density
To assess the need for coastal flood defences, we made a distinction between areas susceptible to coastal flooding and higher, non-susceptible areas. We determined susceptible areas based on possible inundation using coastal flood maps of 1âkm resolution for a 1/1000 year surge level. These maps were created with a global geographic information system (GIS) based inundation model that is forced with a spatially varying sea level, accounting for attenuation of the water level due to land surface roughness82. A method that is more sophisticated compared to a simple âbathtubâ inundation method. Topographic features, as visible in MERIT, protecting the land from flooding are considered. To classify coastlines as urban or rural a distinction was made based on gridded population from the LandScan database83 using the 2UP model84. A transect is characterized âurbanâ if it intersects at least one cell with an urban population with a minimum of 1. Populated coasts have been identified by assigning the population density of the population susceptible to flooding in the proximity of the transects. We used WorldPop201785 population data and assigned population to the transects using a buffer of 15 kilometre radius. The population density is the division of the assigned population and the total area of the assigned cells. This procedure is repeated for buffer radius of 5, 10 and 20âkm, giving fairly comparable outcomes. Following this approach we found a ratio between rural and urban transects of 73/27.
Levee crest heights
The empirical EuroTop formulations47 gave the required levee heights with respect to water levels and wave heights, assuming the presence of a levee at the end of the vegetated foreshore. We hereby neglected the position and characteristics of levees present in the current situation, as no global dataset of coastal protection structures exists. The assumed levee had a standard 1:3 levee profile without berms and an allowed overtopping discharge of 1âlâsâ1âmâ1. These parameters are representative for simple, low-cost levees in developing countries but conservative for well-constructed and maintained levees. Consequently, savings on levee heights in countries with strict protection standards are overestimated, as reduction in required levee height due to vegetation presence is likely less than predicted here. However, this may be balanced out by the fact that we calculated with an average national construction cost per kilometre and levees applying to stricter protection standards may actually be more expensive (Supplementary Fig. 5).
Costs for levee construction and crest height reduction
The calculated levee crest height reductions were monetized using a levee unit price per kilometre length per metre heightening. We used an unit investment costs of levees (metre heightening per kilometre length) of USD 7.0 million42. This estimate represents an average of construction costs in the USA and the Netherlands stated in several studies86,87,88,89. It pertains to all investments costs, including ground work, construction, engineering costs, property or land acquisition, environmental compensation, and project management. Investment costs per metre heightening are well described by a linear function without intercept90. They concluded that for large-scale studies it is sufficient to assume linear costs for each metre of heightening, including the initial costs and the 95% confidence range is between 3x and x/3, where x is the unit cost value. Subsequently we applied three unit levee investment cost prices (low: USD 2.33 million, mid: USD 7.0 million, high: USD 21 million) in line with previous studies42,90. These cost estimates were then adjusted for all other countries by applying construction index multipliers (based on civil engineering construction costs91), to account for differences in construction costs across countries92. Costs were converted to USD2005 power purchasing parity (PPP), to be consistent with the SSPs, using GDP deflators from the World Bank (https://data.worldbank.org/), and annual average market exchange rates between Euros and USD taken from the European Central Bank (unit levee cost per countryâ=âunit levee cost x construction index per country / PPP MER rate 2005 index per country). Example: mid unit levee costsUSAâ=â7.0 Ã1 / 1â=â7.0 million USD2005 PPP km mâ1. If for a country data was not available in the database, we used the average of all countries in the same World Bank income group. For the reference year 2005, this applies to Western Sahara (ESH), North-Korea (PKR) and Somalia (SOM).
Reliability
A scoring table was used to get insight in the reliability of the results of the global analysis. Results were grouped into four reliability classes ranging from âpoorâ to âvery goodâ. Transects were placed in these classes based on data accuracy for three characteristics: hydrodynamics, vegetation and profile elevation. In Supplementary Fig. 6 the (sub) results of the analysis are presented. The first category, hydrodynamics, included known inaccuracies in the hydrodynamic data (GTSM and ERA-I). Data from the GTSM model was considered less reliable in areas with a low tidal range and/or with tropical storms, such as cyclones or hurricanes, as those were not included in our analyses. Also wave data from ERA-I are less reliable in these areas, because the effects of tropical storms are flattened due to the relatively coarse grid size. Hence, transects in these areas were pinpointed by linking them to NOAA data of historical hurricane tracks93. In Supplementary Fig. 6B, areas where tropical storms occur can clearly be recognized. In addition, the Mediterranean Sea, the Red Sea, the Black sea and the Caspian sea stand out in inaccuracy, because of limited tidal action.
Reliability of vegetation characteristics was determined by data source and vegetation width. For transects with extensive vegetation widths, crest height reduction was less sensitive for possible deviations of the vegetation width, due the non-linear relation between vegetation width and wave reduction. Vegetation cover proved most reliable in areas where data from the salt marsh32âand mangrove map14 were available. Hence, this resulted in a âgoodâ score (Supplementary Fig. 6C). Only in cases of extensive vegetation presence was a âvery goodâ score assigned. Transects were appointed as âvery goodâ if vegetation extended 500âm for mangroves, and 1000âm for salt marshes. These thresholds are chosen based on our model results, which show that after ~500âm (salt marshes) and 1000âm (mangroves) maximum reduced wave transmission by foreshore vegetation is reached. Vegetation cover reliability in Europe was classified as âgoodâ, due to reliable vegetation type classification based on CLC30 and the salt marsh map32 in combination accompanied by relatively small vegetation widths. The reliability of the derived vegetation characteristics is especially lacking at the east coast of Canada, at Latin Americaâs south coast, at Africaâs coasts facing the Mediterranean Sea, coasts along the Red Sea and the Persian Gulf, and along the coasts of China, Japan and Russia. For example, in the Persian Gulf states the vegetation presence map tends to falsely identify foreshores as vegetated.
The time-ensemble average (TEA) technique applied for the FAST intertidal elevation map relies on the availability of a reasonable number of images at different tidal stages where the differences in horizontal extent of water coverage can be identified, thus allowing a composite of inundation frequency to be derived. However, the technique is limited by the effective sensor resolution (~30âm, including uncertainty in georeferencing) relative to the horizontal extent of changes in inundation, a function of the tidal range and bed slope. Hence, changes in tidal water extent in microtidal or very high bed-slope regions tend to be too small for reliable discerning differences, leading to poor performance of the technique. However, the merged GEBCO-MERIT dataset was considered less reliable than the FAST intertidal map, based on the resolution and the merging of the two underlying datasets in the intertidal zone. In addition, MERIT tends to overestimate the elevation in mangrove areas, as it measures the canopies as the earthâs surface. Besides the elevation data, the foreshore definition method is used as a profile reliability indicator. The total score per transect is given by the sum of the sub-scores. The sub-scores are normalized to give equal weight to the scoring categories.
Validation
For validation of our method to assess vegetation presence, a comparison of 280 randomly located transects with aerial imagery was carried out. The area accessed in the global assessment was divided in tiles of 90 degrees longitude and 15 degrees latitude. From each tile 6 vegetated and 2 non-vegetated transects were selected. Next, a reference dataset was created by manually identifying vegetation presence using present imagery. Lastly, the vegetation width derived by the model and the manually derived set were compared (Supplementary Fig. 8). For this comparison we made three distinctions, based on (1) vegetation type, (2) foreshore derivation method and (3) vegetation cover source. Comparison showed that the used algorithm on global EO data performs satisfactorily (Supplementary Fig. 8), but in some cases tends to assign a vegetation cover of up to 250âm where there is none. Deviation between observation and the global assessment, is caused by methodological error in the global assessment and inaccuracy in the global datasets, e.g. different timestamps are inevitably compared. This would induce an exaggeration of the effect of vegetation. However, due to the limited dimension of the vegetation extent, the threshold for substantial crest height reduction is falsely exceeded in not more than 2.4% of the cases and the effect is largely balanced out by underestimation of the vegetation cover at larger lengths.
To validate wave reduction by vegetation calculated through our lookup table approach, we compared results with local modelling results for the South-Western part of the Netherlands for 38 vegetated transects. The numerical model SWAN94 in stationary mode was used to translate wave conditions from offshore to nearshore. The simulations were performed with a grid size of 0.01âdeg and bathymetry from EMODNET95. Extreme water levels were included by a water depth correction, using data from GTSR18. Both wind and wave boundary conditions were derived from the earlier described ERA-I re-analysis. The governing wave direction was based on the average of the fifteenth highest wave events in the available wave data. The wind direction was assumed to be aligned with the wave direction. A parametric JONSWAP spectrum shape was used, using a peak enhancement factor of 3.3 and directional spreading of 20 degrees. Foreshore profiles were constructed using an approach similar to foreshore method 2 in the global study but using local high-resolution bathymetry and topography data. Vegetation width was extracted from the salt marsh map32, which was confirmed locally using aerial imagery. Foreshore wave propagation was determined using XBeach in surfbeat mode79.
Our results showed an overestimation of the water depth at the start of the vegetated zone by 0.73âm on average. In addition, the global model derived milder slopes in comparison to the local analysis for narrow vegetated transects. The largest errors were found further away from the mouth of the estuary. Here, the deviation between the wave calculated by SWAN and the depth limited approach is largest. The wave height at the start of the vegetated zone was overestimated on average by 1.12âm, due to the complex geometry and the sheltered configuration of the estuary. The algorithm approximated the wave transmission reduction (RMSE 13%) and the levee crest height reduction relative to the required crest height without vegetation presence (RMSE 19%) with reasonable accuracy (Supplementary Fig. 9).
Sensitivity analysis
A sensitivity analysis has been performed to provide insight in the uncertainty in the presented potential global levee costs savings. The analysis focused specifically on single key parameters, such as the levee unit cost, the critical overtopping discharge and the wave breaker index. High, mid and low levee unit cost scenarios are taken from previous studies42,90. A high, mid, low for the critical overtopping discharge are respectively 10, 1 and 0.1âlâsâ1âmâ1 to incorporate the quality of the levee cover47. We chose RP10 and RP1000 for, respectively, the low and high storm return period scenario. The uncertainty spread of vegetation width is based on the 75% confidence intervals of the underestimated and overestimated vegetation widths of mangroves (+436âm, â136âm) and salt marshes (+597âm, â104âm) in the vegetation presence validation study. For the breaker index we solely chose a high scenario of 0.78, because the index of the global assessment (0.55) was already quite conservative71. For topography we applied a range corresponding to the typical vertical accuracy of the FAST intertidal elevation dataset (±50âcm). Two representative subsets of 500 transects for respectively mangroves and salt marshes have been derived using the clustering method k-means96, based on hydrodynamic conditions, vegetation cover, profile characteristics and geographical location. With these subsets, we repeated the analysis procedure of the global assessment for the sensitivity scenarios. The results point out that the largest spread is caused by the uncertainty in the unit levee cost with â66% and +200% for, respectively, the low and high scenario with respect to the global reference analysis. The other scenarios: topography (â39%, +47%), critical overtopping discharge (â40%, + 40%), storm return period (â28%, +34%), vegetation width (â28%, +39%), breaker index (+21%) (Supplementary Fig. 10). Larger water depths result in a decrease of depth-induced wave energy dissipation and more dissipation due to wave-vegetation interaction, which explains the outcomes of the topography sensitivity results. Similarly, an increase of the storm return period or the breaker index shifts the ratio of wave energy dissipation by wave-bottom interaction and wave-vegetation interaction. The coastal protection costs by vegetation are sensitive to critical overtopping discharge changes, because of the non-linear relation between the wave height in front of the levee and the overtopping discharge47.
Data availability
Data generated in this study (Data for Figs. 1, 2 and 3) have been deposited in the Zenodo database under accession code https://doi.org/10.5281/zenodo.5120632. The FAST intertidal elevation map and the FAST vegetation presence map data are available under restricted access due to the large size of the dataset, access can be obtained freely upon reasonable request to the corresponding author.
Code availability
All general-purpose software packages that we used are open source: Python 2.7.14 (https://www.python.org), NumPy (http://www.numpy.org/), GeoPandas (http://geopandas.org/), Xarray (http://xarray.pydata.org/en/stable/), SciPy (https://www.scipy.org/), Rasterio (https://rasterio.readthedocs.io/en/latest/), Shapely (https://pypi.org/project/Shapely/). Maps were created using Cartopy (v0.18.0. Met Office UK. https://pypi.python.org/pypi/Cartopy/0.18.0) and Matplotlib97 v3.3.4. The software written specifically for this project is available from the corresponding author on reasonable request. The numerical models used in this study are available at: (XBeach) https://oss.deltares.nl/web/xbeach/source-code-and-exe, (SWAN) http://swanmodel.sourceforge.net/.
References
McGranahan, G., Balk, D. & Anderson, B. The rising tide: assessing the risks of climate change and human settlements in low elevation coastal zones. Environ. Urban. 19, 17â37 (2007).
Hinkel, J. et al. Coastal flood damage and adaptation costs under 21st century sea-level rise. Proc. Natl Acad. Sci. USA 111, 3292â3297 (2014).
Güneralp, B., Güneralp, I. & Liu, Y. Changing global patterns of urban exposure to flood and drought hazards. Glob. Environ. Chang. 31, 217â225 (2015).
Borsje, B. W. et al. How ecological engineering can serve in coastal protection. Ecol. Eng. 37, 113â122 (2011).
Temmerman, S. et al. Ecosystem-based coastal defence in the face of global change. Nature 504, 79â83 (2013).
Zhu, Z. et al. Historic storms and the hidden value of coastal wetlands for nature-based flood defence. Nat. Sustain. 3, 853â862 (2020).
Anderson, M., Smith, J. & McKay, S. Wave dissipation by vegetation, CHETN-I-82. US Army Corps of Engineers Engineer Research and Development Center Coastal and Hydraulics Laboratory Vicksburg United States. http://repository.tudelft.nl/view/hydro/uuid:17f2915f-bb57-47f9-8e85-67ff0c82510c/ (2011).
Brinkman, R. M. Wave attenuation in mangrove forests: an investigation through field and theoretical studies (James Cook University, 2006).
Horstman, E. M. et al. Wave attenuation in mangroves: a quantitative approach to field observations. Coast. Eng. 94, 47â62 (2014).
Möller, I. et al. Wave attenuation over coastal salt marshes under storm surge conditions. Nat. Geosci. 7, 727â732 (2014).
Kirwan, M. L. et al. Limits on the adaptability of coastal marshes to rising sea level. Geophys. Res. Lett. 37, L23401 (2010).
Mitsch, W. J., Bernal, B. & Hernandez, M. E. Ecosystem services of wetlands. Int. J. Biodivers. Sci. 11, 1â4 (2015).
Duarte, C. M., Dennison, W. C., Orth, R. J. W. & Carruthers, T. J. B. The charisma of coastal ecosystems: addressing the imbalance. Estuaries Coasts 31, 233â238 (2008).
Giri, C. et al. Status and distribution of mangrove forests of the world using earth observation satellite data. Glob. Ecol. Biogeogr. 20, 154â159 (2011).
Luijendijk, A. et al. The State of the Worldâs Beaches. Sci. Rep. 8, 1â11 (2018).
Murray, N. J. et al. The global distribution and trajectory of tidal flats. Nature 565, 222â225 (2019).
Vousdoukas, M. I. et al. Economic motivation for raising coastal flood defenses in Europe. Nat. Commun. 11, 2119 (2020).
Vousdoukas, M. I. et al. Climatic and socioeconomic controls of future coastal flood risk in Europe. Nat. Clim. Chang. 8, 776â780 (2018).
Arkema, K. K. et al. Coastal habitats shield people and property from sea-level rise and storms. Nat. Clim. Chang. 3, 913â918 (2013).
Narayan, S. et al. The value of coastal wetlands for flood damage reduction in the Northeastern USA. Sci. Rep. 7, 9463 (2017).
Van Coppenolle, R. & Temmerman, S. A global exploration of tidal wetland creation for nature-based flood risk mitigation in coastal cities. Estuar. Coast. Shelf Sci. 226, 106262 (2019).
Chaplin-Kramer, R. et al. Global modeling of natureâs contributions to people. Science 366, 255â258 (2019).
Menéndez, P., Losada, I. J., Torres-Ortega, S., Narayan, S. & Beck, M. W. The global flood protection benefits of mangroves. Sci. Rep. 10, 1â11 (2020).
Vuik, V., Jonkman, S. N., Borsje, B. W. & Suzuki, T. Nature-based flood protection: the efficiency of vegetated foreshores for reducing wave loads on coastal dikes. Coast. Eng. 116, 42â56 (2016).
van Wesenbeeck, B. K., de Boer, W., Narayan, S., van der Star, W. R. L. & de Vries, M. B. Coastal and riverine ecosystems as adaptive flood defenses under a changing climate. Mitig. Adapt. Strateg. Glob. Chang. 22, 1087â1094 (2017).
Van Coppenolle, R., Schwarz, C. & Temmerman, S. Contribution of mangroves and salt marshes to nature-based mitigation of coastal flood risks in major deltas of the world. Estuaries Coasts 41, 1699â1711 (2018).
Barbier, E. B., Georgiou, I. Y., Enchelmeyer, B. & Reed, D. J. The value of wetlands in protecting Southeast Louisiana from Hurricane storm surges. PLoS ONE 8, 58715 (2013).
Schuerch, M. et al. Future response of global coastal wetlands to sea-level rise. Nature 561, 231â234 (2018).
Spencer, T., Brooks, S. M., Evans, B. R., Tempest, J. A. & Möller, I. Southern North Sea storm surge event of 5 December 2013: Water levels, waves and coastal impacts. Earth Sci. Rev. 146, 120â145 (2015).
European Union, Copernicus Land Monitoring Service 2018 & European Environment Agency (EEA). Corine Land Cover (CLC). https://land.copernicus.eu/pan-european/corine-land-cover/ (2018).
Arino, O. et al. Global Land Cover Map for 2009 (GlobCover 2009). European Space Agency (ESA) & Universit catholique de Louvain (UCL). https://doi.org/10.1594/PANGAEA.787668 (2012).
Mcowen, C. et al. A global map of saltmarshes. Biodivers. Data J. 5, e11764 (2017).
Dee, D. P. et al. The ERA-Interim reanalysis: configuration and performance of the data assimilation system. Q. J. R. Meteorol. Soc. 137, 553â597 (2011).
Muis, S., Verlaan, M., Winsemius, H. C., Aerts, J. C. J. H. & Ward, P. J. A global reanalysis of storm surges and extreme sea levels. Nat. Commun. 7, 11969 (2016).
Becker, J. J. et al. Global bathymetry and elevation data at 30 arc seconds resolution: SRTM30_PLUS. Mar. Geod. 32, 355â371 (2009).
Yamazaki, D. et al. A high-accuracy map of global terrain elevations. Geophys. Res. Lett. 44, 5844â5853 (2017).
Bao, T. Q. Effect of mangrove forest structures on wave attenuation in coastal Vietnam. Oceanologia 53, 807â818 (2011).
De Conto, R. M. & Pollard, D. Contribution of Antarctica to past and future sea-level rise. Nature 531, 591â597 (2016).
Silver, J. M. et al. Advancing coastal risk reduction science and implementation by accounting for climate, ecosystems, and people. Front. Mar. Sci. 6, 556 (2019).
Romañach, S. S. et al. Conservation and restoration of mangroves: global status, perspectives, and prognosis. Ocean Coast. Manag. 154, 72â82 (2018).
Lenk, S., Rybski, D., Heidrich, O., Dawson, R. J. & Kropp, J. P. Costs of sea dikes-regressions and uncertainty estimates. Nat. Hazards Earth Syst. Sci. 17, 765â779 (2017).
Ward, P. J. et al. A global framework for future costs and benefits of river-flood protection in urban areas. Nat. Clim. Chang. 7, 642â646 (2017).
Turner, R. K., Burgess, D., Hadley, D., Coombes, E. & Jackson, N. A cost-benefit appraisal of coastal managed realignment policy. Glob. Environ. Chang. 17, 397â407 (2007).
Broekx, S., Smets, S., Liekens, I., Bulckaen, D. & De Nocker, L. Designing a long-term flood risk management plan for the Scheldt estuary using a risk-based approach. Nat. Hazards 57, 245â266 (2011).
Spalding, M. & Parrett, C. L. Global patterns in mangrove recreation and tourism. Mar. Policy 110, 103540 (2019).
Ward, P. J. et al. Usefulness and limitations of global flood risk models. Nat. Clim. Chang. 5, 712â715 (2015).
Pullen, T., Allsop, N. W. H., Bruce, T., Kortenhaus, A. & Van der Meer, J. W. EurOtop wave overtopping of sea defences and related structures: assessment manual (Kuratorium für Forschung im Küsteningenieurwesen, 2007).
Hu, K., Chen, Q. & Wang, H. A numerical study of vegetation impact on reducing storm surge by wetlands in a semi-enclosed estuary. Coast. Eng. 95, 66â76 (2015).
Vuik, V., Heo, H. Y. S., Zhu, Z., Borsje, B. W. & Jonkman, S. N. Stem breakage of salt marsh vegetation under wave forcing: a field and model study. Estuar. Coast. Shelf Sci. 200, 41â58 (2018).
Primavera, J. H. et al. Preliminary assessment of post-Haiyan mangrove damage and short-term recovery in Eastern Samar, central Philippines. Mar. Pollut. Bull. 109, 744â750 (2016).
Simard, M. et al. Mangrove canopy height globally related to precipitation, temperature and cyclone frequency. Nat. Geosci. 12, 40â45 (2019).
Wang, Z. B., Elias, E. P. L., Van Der Spek, A. J. F. & Lodder, Q. J. Sediment budget and morphological development of the Dutch Wadden Sea: Impact of accelerated sea-level rise and subsidence until 2100. Geol. en. Mijnb./Neth. J. Geosci. 97, 183â214 (2018).
Cowell, P. J. & Thom, B. G. Morphodynamics Of Coastal Evolution. in Coastal Evolution: Late Quaternary Shoreline Morphodynamics 33â86 (Cambridge University Press, Cambridge, 1994).
Doody, J. P. Coastal squeeze and managed realignment in southeast England, does it tell us anything about the future? Ocean Coast. Manag 79, 34â41 (2013).
McIvor, A. L., Spencer, T., Möller, I. & Spalding, M. The response of mangrove soil surface elevation to sea level rise. Natural Coastal Protection Series: Report 3. Cambridge Coastal Research Unit Working Paper 42. Published by The Nature Conservancy and Wetlands International. 59 pages. ISSN 2050â7941. http://coastalresilience.org/science/mangroves/surface-elevation-and-sea-level-rise (2013).
Kirwan, M. L., Temmerman, S., Skeehan, E. E., Guntenspergen, G. R. & Fagherazzi, S. Overestimation of marsh vulnerability to sea level rise. Nat. Clim. Chang. 6, 253â260 (2016).
Lovelock, C. E. et al. The vulnerability of Indo-Pacific mangrove forests to sea-level rise. Nature 526, 559â563 (2015).
Willemsen, P. W. J. M., Horstman, E. M., Borsje, B. W., Friess, D. A. & Dohmen-Janssen, C. M. Sensitivity of the sediment trapping capacity of an estuarine mangrove forest. Geomorphology 273, 189â201 (2016).
Saintilan, N. et al. Thresholds of mangrove survival under rapid sea level rise. Science 368, 1118 LPâ1121 (2020).
Woodroffe, C. D. et al. Mangrove sedimentation and response to relative sea-level rise. Annu. Rev. Mar. Sci. 8, 243â266 (2016).
Barbier, E. B. Valuing the storm protection service of estuarine and coastal ecosystems. Ecosyst. Serv. 11, 32â38 (2015).
Winterwerp, J. C., Erftemeijer, P. L. A., Suryadiputra, N., Van Eijk, P. & Zhang, L. Defining eco-morphodynamic requirements for rehabilitating eroding mangrove-mud coasts. Wetlands 33, 515â526 (2013).
Nicholls, R. J. et al. A global analysis of subsidence, relative sea-level change and coastal flood exposure. Nat. Clim. Chang. 11, 338â342 (2021).
Best, S. N. et al. Do salt marshes survive sea level rise? Modelling wave action, morphodynamics and vegetation dynamics. Environ. Model. Softw. 109, 152â166 (2018).
Winterwerp, J. C. et al. Managing erosion of mangrove-mud coasts with permeable damsâlessons learned. Ecol. Eng. 158, 106078 (2020).
Borchert, S. M., Osland, M. J., Enwright, N. M., Kereen & Griffith, T. Coastal wetland adaptation to sea level rise: quantifing potential for landward migration and coastal squeeze. J. Appl. Ecol. 55, 2876â2887 (2018).
Gedan, K. B., Kirwan, M. L., Wolanski, E., Barbier, E. B. & Silliman, B. R. The present and future role of coastal wetland vegetation in protecting shorelines: answering recent challenges to the paradigm. Clim. Chang. 106, 7â29 (2011).
OpenStreetMap (OSM) contributors. Planet dump. https://planet.osm.org. https://osmdata.openstreetmap.de/data/coastlines.html (2015).
Calero, J. S. et al. FAST MI-SAFE platform: Foreshore Assessment using Space Technology. in Proceedings on the 2017 Conference on Big Data from Space (eds. Soille, P. & Marchetti, P. G.) 247â250 (Publications Office of the European Union, Luxembourg, 2017).
Vafeidis, A. T. et al. A new global coastal database for impact and vulnerability analysis to sea-level rise. J. Coast. Res. 244, 917â924 (2008).
Camenen, B. & Larson, M. Predictive formulas for breaker depth index and breaker type. J. Coast. Res. 23, 1028â1041 (2007).
Cole, T. G., Ewel, K. C. & Devoe, N. N. Structure of mangrove trees and forests in Micronesia. Forest Ecol. Manag. 117, 95â109 (1999).
Mazda, Y. et al. Drag force due to vegetation in mangrove swamps. Mangroves Salt Marshes 1, 193â199 (1997).
Krauss, K. W., Allen, J. A. & Cahoon, D. R. Differential rates of vertical accretion and elevation change among aerial root types in Micronesian mangrove forests. Estuar. Coast. Shelf Sci. 56, 251â259 (2003).
Narayan, S. et al. On the effectiveness of mangroves in attenuating cyclone-induced waves. in Proceedings of the 32nd International Conference on Coastal Engineering (eds. McKee-smith, J. & Lynett, P.) (American Society of Civil Engineers, 2010).
Zhang, X., Chua, V. P. & Cheong, H. F. Hydrodynamics in mangrove prop roots and their physical properties. J. Hydro Environ. Res. 9, 281â294 (2015).
Möller, I. Quantifying saltmarsh vegetation and its effect on wave height dissipation: results from a UK East coast saltmarsh. Estuar. Coast. Shelf Sci. 69, 337â351 (2006).
Jadhav, R. S., Chen, Q. & Smith, J. M. Spectral distribution of wave energy dissipation by salt marsh vegetation. Coast. Eng. 77, 99â107 (2013).
Van Rooijen, A. A. et al. Modeling the effect of wave-vegetation interaction on wave setup. J. Geophys. Res. Ocean. 121, 4341â4359 (2016).
Dalrymple, R. A., Kirby, J. T. & Hwang, P. A. Wave diffraction due to areas of energy-dissipation. J. Waterw. Port. Coast. Ocean Eng. 110, 67â79 (1984).
Mendez, F. J. & Losada, I. J. An empirical model to estimate the propagation of random breaking and nonbreaking waves over vegetation fields. Coast. Eng. 51, 103â118 (2004).
Haer, T. et al. Coastal and river flood risk analyses for guiding economically optimal flood adaptation policies: a country-scale study for Mexico. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 376, 20170329 (2018).
Bright, E. A., Coleman, P. R., Rose, A. N. & Urban, M. L. LandScan 2010. https://landscan.ornl.gov/ (2011).
van Huijstee, J., Bouwman, A. & van Rijn, F. Towards and Urban Preview: Modelling future urban growth with 2UP. https://www.pbl.nl/sites/default/files/cms/publicaties/pbl-2018-Towards-an-urban-preview_3255.pdf (2018).
Lloyd, C. T., Sorichetta, A. & Tatem, A. J. Data Descriptor: high resolution global gridded data for use in population studies. Sci. Data 4, 170001 (2017).
Bos, A. J. Optimal safety level for the New Orleans East polder; A preliminary risk analysis (University of Amsterdam, 2008).
de Grave, P. & Baarse, G. Kosten van maatregelen. Informatie ten behoeve van het projectWaterveiligheid 21e eeuw. https://puc.overheid.nl/rijkswaterstaat/doc/PUC_138415_31/ (2011).
Aerts, J. C. J. H., Botzen, W. J. W., de Moel, H. & Bowman, M. Cost estimates for flood resilience and protection strategies in New York City. Ann. N. Y. Acad. Sci. 1294, 1â104 (2013).
Jonkman, S. N., Hillen, M. M., Nicholls, R. J., Kanning, W. & van Ledden, M. Costs of adapting coastal defences to sea-level riseânew estimates and their implications. J. Coast. Res. 29, 1212â1226 (2013).
Lenk, S., Rybski, D., Heidrich, O., Dawson, R. J. & Kropp, J. P. Costs of sea dikesâregressions and uncertainty estimates. Nat. Hazards Earth Syst. Sci. 17, 765â779 (2017).
The 2009 Global Construction Cost and Reference Yearbook (Compass International Consultants, 2009).
Ward, P. J. et al. Partial costs of global climate change adaptation for the supply of raw industrial and municipal water: A methodology and application. Environ. Res. Lett. 5, 044011 (2010).
ESRI/National Oceanic and Atmospheric Association. NOAAHurricaneTracks [ESRI Shape Files]. https://maps4.arcgisonline.com/ArcGIS/rest/services/A-16/NOAA_Hurricane_Tracks/MapServer (2011).
Booij, N., Ris, R. C. & Holthuijsen, L. H. A third-generation wave model for coastal regions. J. Geophys. Res. 104, 7649â7666 (1999).
Schaap, D. & Moussat, E. in Geophysical Research Abstracts, vol. 15, 3620 (EGU General Assembly, Vienna, Austria, 2013).
Pedregosa, F. et al. Scikit-learn: machine learning in Python. J. Mach. Learn. Res. 12, 2825â2830 (2011).
Hunter, J. D. Matplotlib: a 2D graphics environment. Comput. Sci. Eng. 9, 90â95 (2007).
Acknowledgements
The research leading to these results received funding from the Aqueduct Global Flood Analyzer project, via subsidy 5000002722 from the Netherlands Ministry of Infrastructure and the Environment. The latter project is convened by the World Resources Institute. This research is part of Foreshore Assessment using Space Technology (FAST, 2014â2018). A project funded by the European Unionâs (EU) Seventh Framework Programme (FP7) for research, technological development and demonstration under grant agreement number 607131. FAST is developing downstream services for the European Earth Observation Programme Copernicus to support cost-effective, nature-based shoreline protection against flooding and erosion. P.J.W. received additional funding from the Dutch Research Council (NWO), in the form of a VIDI research grant 016.161.324.
Author information
Authors and Affiliations
Contributions
V.T.M.Z. wrote a substantial part of the manuscript, wrote and ran the (post) processing pipeline and produced all figures unless otherwise indicated. J.T.D. participated in conceiving the conceptual method, wrote a substantial part of the manuscript, setup the preliminary version of the final processing pipeline, wrote parts of the final processing pipeline. B.K.W. participated in conceiving the conceptual method, wrote a substantial part of the manuscript and provided photos for Fig. 4. D.E. wrote (major) parts of the processing pipeline and contributed to the manuscript. E.P.M. developed the FAST intertidal elevation map, assisted with the FAST coastal vegetation map, contributed to the manuscript and produced Supplementary Fig. 7. H.C.W. Conceived experiments and contributed to the manuscript. P.J.W. provided construction index multipliers and GDP deflators at country level for levee costs calculations and contributed to the manuscript. M.B.V. participated in conceiving the conceptual method and contributed to the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Communications thanks Dominik Paprotny and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
van Zelst, V.T.M., Dijkstra, J.T., van Wesenbeeck, B.K. et al. Cutting the costs of coastal protection by integrating vegetation in flood defences. Nat Commun 12, 6533 (2021). https://doi.org/10.1038/s41467-021-26887-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-021-26887-4
This article is cited by
-
The value of marsh restoration for flood risk reduction in an urban estuary
Scientific Reports (2024)
-
Urban flood risk management needs nature-based solutions: a coupled social-ecological system perspective
npj Urban Sustainability (2024)
-
Integrating mangrove growth and failure in coastal flood protection designs
Scientific Reports (2024)
-
Habitability of low-lying socio-ecological systems under a changing climate
Climatic Change (2024)
-
A paradigm shift in the quantification of wave energy attenuation due to saltmarshes based on their standing biomass
Scientific Reports (2022)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.