Front. Mar. Sci. Frontiers in Marine Science Front. Mar. Sci. 2296-7745 Frontiers Media S.A. 10.3389/fmars.2025.1526082 Marine Science Original Research Evaluating mangroves as nature-based solutions for coastal protection under current and future sea level rise scenarios Jayson-Quashigah Philip-Neri 1 2 * Staneva Joanna 1 Chen Wei 1 Djath Bughsin’ 1 Mahu Edem 3 Appeaning Addo Kwasi 3 1 Institute of Coastal Systems-Analysis and Modeling, Helmholtz-Zentrum, Hereon, Geesthacht, Germany 2 Institute for Environment and Sanitation Studies, University of Ghana, Accra, Ghana 3 Department of Marine and Fisheries Sciences, University of Ghana, Accra, Ghana

Edited by: Sergio Maldonado, Tecnologico de Monterrey, Mexico

Reviewed by: Songdong Shao, Dongguan University of Technology, China

Denny Nugroho Sugianto, Diponegoro University, Indonesia

*Correspondence: Philip-Neri Jayson-Quashigah, pnjquashigah@ug.edu.gh

19 05 2025 2025 12 1526082 11 11 2024 28 04 2025 Copyright © 2025 Jayson-Quashigah, Staneva, Chen, Djath, Mahu and Appeaning Addo 2025 Jayson-Quashigah, Staneva, Chen, Djath, Mahu and Appeaning Addo

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Coastal hazards, such as erosion and flooding, pose significant threats to many coastal areas, and in extreme cases, have led to the decimation of some coastal communities. Historically, management has mostly focused on the use of grey infrastructure such as seawalls and groynes. However, these interventions are costly and can cause unintended consequences, promoting a shift toward Nature-based Solutions (NbS), such as the use of mangroves. Mangroves, particularly, have been proven to protect shorelines due to their ability to attenuate waves and trap sediment. However, there is limited research, particularly in data-deficient regions such as Ghana, West Africa, to support such NbS initiatives. This study explored the potential of mangroves as NbS to mitigate coastal erosion, using the 1D morphodynamic model XBeach. The baseline model was validated against measured coastal profiles, and the results show accurate predictions of sediment volume changes with an overall RMSE of 0.75 m. Based on the calibrations, we explored the effects of mature mangroves at varying densities on the berm and within the intertidal zone under current and projected sea level conditions, on coastal erosion. The results show a significant reduction in sediment volume erosion from 28 m3 to 0.9 m3 in the current situation, representing 97% protection; and from 468 m3 to 2.6 m3 under future sea level rise of 0.233 m by 2040, indicating 99% protection. Notably, high densities of mangroves, introduced on the berm, which is more practical for the area, provided up to 53% reduction in erosion for the current situation and 97% for the future. These scenario-based simulations demonstrate the potential of mangroves as a dynamic coastal defense strategy, with the approach providing a valuable tool for testing and optimizing NbS interventions.

nature-based solutions XBeach mangroves coastal protection sea level rise Volta Delta section-in-acceptance Coastal Ocean Processes

香京julia种子在线播放

    1. <form id=HxFbUHhlv><nobr id=HxFbUHhlv></nobr></form>
      <address id=HxFbUHhlv><nobr id=HxFbUHhlv><nobr id=HxFbUHhlv></nobr></nobr></address>

      Introduction

      Recently, coastlines have come under intense pressure from coastal hazards, with changing climate and anthropogenic factors playing a key role. Climate change is driving global sea level rise, leading to increased flooding, inundation, and more frequent storm surges, all of which contribute to severe erosion in many vulnerable places. Consequently, coastal erosion, flooding, and inundation are on the rise, particularly along sandy coasts, which make up about one-third of the world’s coastlines. Historical assessments have revealed a dramatic loss of approximately 28,000 km² of land globally between 1984 and 2015, with some beaches experiencing losses of over 600 m cross-shore (Mentaschi et al., 2018). The study identified a combination of sea level rise, subsidence, and human activity as the underlying cause of this global land loss. Further, projections show that by 2100, sea level rise could be responsible for 73 to 85% of coastline retreat, under a business-as-usual scenario (Vousdoukas et al., 2020). Under the same scenario, the global average 100-year extreme sea levels of 58–172 mm are projected to increase between 2000 and 2100 (Vousdoukas et al., 2018).

      Within these vulnerable sandy coastlines, deltaic regions have been recognized as more prone to coastal erosion and flooding due to their low-lying nature and geomorphology (Ericson et al., 2006; Scown et al., 2023). The Volta Delta in Ghana, for example, is considered at high risk of relative sea level rise, leading to increased exposure to flooding and erosion (Tessler et al., 2015). Historically, coastal erosion has been a major problem in the Volta Delta (Appeaning Addo et al., 2018; Aagaard et al., 2021; Mann et al., 2023), which has led to the destruction of coastal infrastructure, the displacement of people, and in some instances, the decimation of communities (Jayson-Quashigah et al., 2021). Traditional interventions, predominantly “grey” or “hard” engineered structures, have been implemented to manage this situation. Examples include the Keta Sea Defense Project, which comprised groynes, revetments, and beach nourishment (Nairn and Dibajnia, 2004), and the Ada Sea Defense Project, which similarly included nourishment and groynes (Bolle et al., 2015). These projects have been able to hold the shoreline to some extent at the sites, but a common problem identified is the transfer of erosion downdrift of the structures (Angnuureng et al., 2013; Jayson-Quashigah et al., 2013, 2021; Mann et al., 2023). Additionally, the cost of building and maintaining such defense structures is exorbitant.

      With the projected increase in the threats of coastal erosion, flooding, and inundation, and the negative impacts of grey interventions, the focus in recent decades has shifted to the adoption of nature-based solutions (NbS) and other non-traditional, integrated interventions to tackle the problem (Charoenlerkthawin et al., 2022; Singhvi et al., 2022). Over the years, several studies have explored the capability of NbS, such as the use of seagrass and other coastal vegetation, including mangroves, to protect the coast against storms, erosion, and flooding (James et al., 2021; van Zelst et al., 2021; Amos and Akib, 2023; Thao et al., 2023). Chen et al. (2022), for example, demonstrated the concept of green nourishment where seagrass is planted on nearshore beach nourishment to assess the impact on waves and hence erosion (see also Vuik et al., 2016; Chen et al., 2024). Other studies such as van Zelst et al. (2021), have demonstrated that vegetated foreshores can be used to supplement conventional engineering approaches. Furthermore, mangroves have demonstrated the capability of attenuating waves and therefore protecting the shoreline from storms (Khanh Phan, 2019; Thao et al., 2023; van Hespen et al., 2023).

      Along the West African coast, including the coast of Ghana, there have been calls by scientists for a shift towards more nature-based solutions to manage coastal erosion, promoting a more environmentally friendly and cost-effective approach (Hagedoorn et al., 2021; Wegman et al., 2023). However, within this region, there is little data and research currently supporting this call. There remains a substantial gap in ocean observation and model assessments in the region. Existing model assessments have mainly focused on the design and implementation of grey infrastructure (Nairn and Dibajnia, 2004; Bolle et al., 2015). Assessing the NbS options is, however, critical to guide policy direction and implementation. Addressing this gap, the UN Ocean Decade endorsed as a decade action the project “Mangroves as a Nature-based Solution for Coastal Hazards (MANCOGA)”, for which this study has been conducted.

      This pioneering study aims to assess the effectiveness of using mangroves that grow naturally within the region as a means of protecting the eastern coast of Ghana from increasing erosion. Utilizing the open-access XBeach model with the vegetation component, the approach simulates What-if Scenarios (WiS) for both current and future sea levels. By introducing mature mangroves along the beach at varying densities, the study estimates their protective ability. Sediment dynamics and shoreline stability are assessed under each scenario, providing relevant data for decision-making regarding coastal management strategies in the Volta Delta and similar coastal regions worldwide. This study is novel in its application of process-based modeling to simulated NbS in a region where empirical and modeling research remains scarce. This provides the first region-specific modeling evidence to support mangrove restoration and integration into coastal management policies.

      The structure of this paper is as follows: The materials and methods section explains the study site and the morphodynamic modelling used to simulate the protective effects of mangroves. The results section presents findings from the What-if Scenarios that demonstrate how different mangrove densities influence erosion. The discussion explores the broader implications of the results, including the benefits of using mangroves as an alternative to grey infrastructure with rising sea levels. Finally, the conclusion summarizes the findings and their relevance to coastal management and future research directions.

      Materials and methods Site description

      This case study was carried out along a sandy coastline of the Volta Delta of Ghana, with a selected profile east of the Volta estuary ( Figure 1 ). The area is considered a high-energy zone with nearshore wave heights in the range of 1–2 m 75% of the time (Verheyen et al., 2014). The waves approach the coast in a predominantly SSW direction (between 170°-210°) with periods exceeding 10 s 80% of the time (Verheyen et al., 2014). These waves generate strong longshore currents in the eastward direction and are responsible for the high rates of longshore sediment transport reported in the area (Nairn and Dibajnia, 2004; Anthony et al., 2016). The tide is semi-diurnal with an average tidal range of about 1 m (Wiafe et al., 2013; Appeaning Addo et al., 2018). Generally, the beaches are considered relatively steep with reported slopes ranging between 1:3 and 1:15 (Bollen et al., 2011; Roest, 2018). Studies have shown that the area is historically eroding (Mann et al., 2023), with recent short-term rates for some areas reaching as high as 30 m/year (Jayson-Quashigah et al., 2019), leading to the destruction of communities such as Fuveme and Agavedzi along this coast. The location falls within the ongoing project looking at the use of mangroves as NbS for coastal hazards in eastern Ghana (MANCOGA) project (MANCOGA, 2024), as well as the ongoing Harmony Coast project, which has begun regular mapping of bathymetry and beach profiles along this coast (Angnuureng, 2023; Angnuureng et al., 2024).

      Location of the study area. The red line indicates the model domain.

      The Volta Delta area is known to host some of the most extensive and dense mangroves along the coast of Ghana (Awuku-Sowah et al., 2023; Ofori et al., 2023). Among the over five mangrove species identified in Ghana (Ofori et al., 2023), three dominant species are present in the study area, namely Avicennia germinans (black mangrove), Rhizophora racemosa (red mangrove), and Laguncularia racemosa (white mangrove) (Nunoo and Agyekumhene, 2022). Though the red mangroves are more dominant in the area, the white mangrove is noted to be more prevalent closer to the coastline, and this is confirmed by our observations ( Figure 2 ).

      Parches of Laguncularia racemosa (white mangrove) along the beach within the Volta Delta. Sediment-trapping effects can be observed (elevating the dune).

      Model framework

      A simple framework was adapted for this study, based on XBeach, similar to that of Chen et al. (2022) and Chen et al. (2024) who employed the approach to assess the use of seagrass as a nature-based solution for coastal erosion. The model has also been successfully applied at the global level to assess the role of vegetation in coastal defense (van Zelst et al., 2021). Here, a near-shore morphodynamic model is set up and forced with global ocean conditions from ERA5 reanalysis data ( Figure 3 ). This choice ensures that the model mimics ecological realism as closely as possible, enhancing the transferability of the results for practical application. The XBeach model simultaneously solves the time-dependent short-wave action balance, the roller energy equations, the nonlinear shallow water equations of mass and momentum, sediment transport formulations, and bed updates on the scale of wave groups (Roelvink et al., 2009). The study adopts the 1D surfbeat mode (instationary) to simulate coastal erosion under various scenarios (Chen et al., 2024).

      Model Framework.

      Input data (current situation)

      Data for the modelling was compiled from various sources. Nearshore bathymetry was based on work that was carried out as part of the Harmony Coast project (Angnuureng, 2023; Angnuureng et al., 2024). The Harmony Coast project has initiated regular mapping of the nearshore bathymetry within the Volta delta using a Valeport Midas Echo-sounder with depths ranging between 3–10 m. The average spacing between the transects used for the mapping is 500 m. Beyond the 10 m depth, the GEBCO global data was utilized (IHO-IOC, 2019) to complement the nearshore data. Beach profile data was also acquired from the Harmony Coast project. The profiles were measured using a Real-Time Kinematic Differential Global Positioning System (RTK-DGPS) coordinated with existing ground control points established along the coast (Angnuureng et al., 2024). The profiles were measured in June and July 2023 for the study area. Due to the unavailability of measurements across the surf zone, an equilibrium beach profile was estimated using the Equation 1:

      h ( y ) = A y 2 / 3

      where h is the depth at seaward distance y, and A is a scale parameter that depends on sediment characteristics (Dean, 1991).

      The equilibrium beach profile was used to test possible surf zone behavior, and the most probable initial surf zone profile (based on shoreline response) was selected and used for the final model setup. For boundary conditions, the fifth-generation European Centre for Medium-Range Weather Forecasts (ECMWF) atmospheric reanalysis (ERA5) produced by the Copernicus Climate Change Service (C3S) is used (Hersbach et al., 2023; Soci et al., 2024). The hourly data was downloaded for the closest point to the coast (1° E, 5.5° N). The parameters extracted include the significant height of combined wind waves and swell, peak wave period, and mean wave direction. The data was adjusted based on earlier studies by Giardino et al. (2018) (see also Jayson-Quashigah et al., 2021), where the calibration of the ERA data using altimeter data led to an increase of wave heights by 10% and wave periods by 4.8%. For this study, the offshore boundary wave heights were therefore increased by 10% to compensate for the underestimation reported by these studies. Water level from the tidal gauge at the Tema Port (Ghana Ports, 2025) was used ( Figure 4 ).

      Hydrodynamic variables used in the modeling include sea level, significant wave height (Hs), mean wave period (Tp) and mean wave direction (Dir).

      Input data (future scenario)

      Using the 2023 dataset as the current situation, a future scenario was also simulated. For the future scenario, the Intergovernmental Panel on Climate Change’s (IPCC’s) regional sea level rise (SLR) projection based on the Sixth Assessment Report (AR6) was used (IPCC, 2021). The near-term (2021–2040) timeframe was selected. This time frame provides a more immediate and actionable timeframe for policymakers and stakeholders to implement adaptive measures. Additionally, only the Shared Socioeconomic Pathway (SSP) 5–8.5, which is marked by extreme greenhouse gas emissions, was chosen, accounting for the worst-case scenario. Further, only the 83rd percentile of sea level rise was chosen to capture the higher probability range and account for significant uncertainties. For this scenario, a sea level rise of 0.233 m is projected (Avornyo et al., 2024). The mean sea level was therefore increased by 0.233 m based on the 2023 baseline level.

      For the future scenario, the wave conditions were assumed to remain the same. Hence, the same wave boundary conditions were used. Scott et al. (2021) discuss how atmospheric processes can modulate wave height and frequency, but emphasize that the overarching influence of SLR is expected to remain dominant in driving long-term coastal erosion. The attribution of erosion to sea-level dynamics is further supported by the findings of Mentaschi et al. (2018) which show long-term observations revealing an unequivocal linkage between SLR and shoreline erosion.

      Model setup, calibration, and validation

      A cross-shore profile was generated from the bathymetry and topography data extending approximately 5 km offshore ( Figure 5 ) with a varying grid size ranging from 20 m offshore to 1 m nearshore. The model simulation spans two (2) months between June and July 2023, where the highest erosion was recorded based on the historical profiles that were mapped for the location (Angnuureng et al., 2024). The higher eroding period was chosen to test the ability of mangroves to protect the shoreline during such extreme events.

      Model domain indicating initial bathymetry with equilibrium profile interpolation for the surf zone. The Mean Sea level is at 0.8 m.

      The boundary conditions were forced with the ERA5 reanalysis data as described using the time-varying JONSWAP spectra with the tidal data from the Tema port.

      For the XBeach model, the Chezy bed roughness coefficient is implemented, with a value of 55.0 (XBeach Team, 2023; version 1.23). A grain size of 0.8 mm (D50) and 0.12 mm (D90) was used based on the average of the measurements reported for the area from earlier studies (Verheyen et al., 2014; Jayson-Quashigah et al., 2019). Calibration was conducted using several parameters to determine which combinations produced the closest result of erosion based on the measured June and July profiles. Among all the parameters, three were critical, as also indicated by other studies (Roelvink and Costas, 2019; Kombiadou et al., 2021); the bermslope, which allows the slope of the profile near the waterline to be nudged towards the given value; the facAs, which is the time-averaged flows due to wave asymmetry and delta, which represents the fraction of wave height added to the water depth to adjust maximum wave height in wave breaking formulations (Kombiadou et al., 2020). For this model, the longshore transport component(lsgrad) was activated and calibrated with reference to other studies (Pender and Karunarathna, 2013; Roelvink and Costas, 2019). This was necessary to take care of the longshore transport, which is strong in the area (Nairn and Dibajnia, 2004; Anthony et al., 2019) and their sensitivity was tested with several values to determine the optimal value, ensuring the robustness and reliability of the results. The key parameter settings that worked best for the study area are summarized in Table 1 . All other parameters were left at the recommended default settings (XBeach Team, 2023).

      Key model parameters.

      Parameter Best value Source
      facAs 0.35 Gruwez et al., 2014; Verheyen et al., 2014; Vousdoukas et al., 2012
      morfac 5 Kombiadou et al., 2021
      turb wave_averaged Vousdoukas et al., 2012
      D50 & 90 0.8 mm & 1.2 mm Jayson-Quashigah et al., 2019; Gruwez et al., 2014
      delta 0.4 Based on calibration
      gammax 1.5 Kombiadou et al., 2021
      lsgrad -0.0002 m-1 Based on calibration (ref., Roelvink and Costas, 2019)
      kx 0.01 Gruwez et al., 2014; Verheyen et al., 2014
      bermslope 0.18 Based on calibration

      The calibration of the model was mainly guided by the morphological changes observed due to the absence of adequate in situ historical data on the ocean state. The modelled changes were validated against the measured profile changes (Angnuureng et al., 2024). Shoreline erosion based on the location of the berm crest and volume erosion across the profile was estimated between June and July 2023. The profiles, however, do not extend beyond the low water mark, hence, changes in the surf zone were not considered for the calibration and validation.

      NbS integration (what-if scenarios)

      The vegetation module of the XBeach was activated to introduce mangroves into the model. It has been used elsewhere to assess the role of vegetation in wave attenuation (Chen et al., 2022, 2024; van Hespen et al., 2023). The mangrove vegetation is represented as rigid cylinders with parameters including stem height, diameter, and density (Burger, 2005; Chen et al., 2022). The parameters were set based on in situ data collected from the study area on mangroves (with specific reference to the white mangrove). A conservative drag coefficient (CD ) of 1 was adopted (Adytia et al., 2019; Yoshikai et al., 2022; Lopez-Arias et al., 2024) and increased by a factor of 0.1 for every higher level of mangrove density. The height of the mangroves measured in the area varies, ranging between 2 m and 8 m. However, those observed close to the beach hardly exceed 5 m in height.

      Table 2 summarizes the mangrove characteristics for the various density levels used for the model based on field measurements and Thao et al. (2023).

      Mangrove characteristics.

      Density Root Trunk Canopy
      N (m2) bv (m) ah(m) N (m2) bv (m) ah(m) N (m2) bv (m) ah(m)
      Sparse 15 0.02 0.04 0.3 0.04 1.5 30 0.01 1.5
      Medium 45 0.05 0.10 0.6 0.08 2 90 0.01 2
      Dense 70 0.07 0.16 1.2 0.16 2.5 140 0.01 2.5

      N, number per square meter; bv, diameter; ah, height.

      Three (3) broad scenarios were considered, (1) no presence of mangroves on the beach (Baseline), for both the current situation and projected sea level rise for 2040 (2) the introduction of matured mangroves on the berm (behind the crest) at varying densities (Scn_Man-I), for the current situation and projected sea level rise for 2040 and (3) the introduction of matured mangroves in the intertidal zone at varying densities for both current and projected sea level rise for 2040 (Scn_Man-II).

      Predominantly, the mangroves occur at the back of the beach (normally where there are creeks or lagoons behind the dune system), but there are some instances of them occurring on the dune along pockets of the beach (see Figure 2 ). However, their ability to protect the coastline has not been tested. The third scenario, though, does not occur naturally in the area, but was considered a possible option ( Figure 6 ).

      Scenarios (Scn_Man-I; mangroves on the berm and Scn_Man-II; mangroves in the intertidal zone).

      Results Model calibration and sensitivity

      The model’s sensitivity to coastal erosion was tested using three main parameters, namely the bermslope, facAs, and delta (Roelvink and Costas, 2019; Kombiadou et al., 2020). A range of these variables was simulated using the June 2023 profile as a baseline, and the resultant profiles were compared to the July 2023 profile. For facAs, the values tested were 0.25, 0.30, 0.35, 0.40, and 0.45, for delta; 0.30, 0.35, 0.40, 0.45, and 0.50 were tested, however, results for 0.45 and 0.50 were inconsistent and therefore not reported. Bermslope values tested were 0.14, 0.15, 0.16, 0.17, and 0.18, respectively. From the results, the model was sensitive to all three variables tested with the most significant changes observed with variations in facAs ( Figure 7 ).

      Calibration results for facAs 0.1-0.5, delta 0.3- 0.7, and bermslope 0.12-0.2. The area represented by dashed lines was not used for validation since there was no measured data. The best values are represented by the green color.

      The Root Mean Square Error (RMSE) for facAs calibration ranged between 0.71 m to 0.90 m, with the highest recorded for facAs 0.45. For variations in delta, the RSME ranged between 0.79 m to 0.89 m, with delta 0.4 producing the lowest value. RSME values for the bermslope evaluated ranged between 0.73 m to 0.85 m, with the lowest value captured for the bermslope value of 18. It should be noted that these errors were only calculated considering the backshore and foreshore zones where there was measured data.

      Scenario results for the current situation (2023) Baseline (no mangroves on the beach)

      For the baseline, the focus was to simulate the coastal erosion without any intervention, mimicking the existing situation as closely as possible. The profile measurements conducted in June and July 2023 were used to initiate and validate the morphological simulations. Using the berm as a proxy, the measured profiles show erosion of approximately 11 m from June to July ( Figure 8 ) with an average depth of erosion of 1.4 m. This translates to approximately 28 m3 of sediment volume erosion, with 46% of that (13 m3) being deposited at the back beach through overwash (this does not include the surf zone due to lack of measured data). The beach face slope was approximately 4.7 degrees.

      Erosion of the bed and dune without mangrove.

      For the baseline simulations, the best results were obtained with the calibration factors: facAs 0.35, delta 0.4, and bermslope 0.18. This setup was able to accurately predict the berm erosion of 11 m observed from June to July, 2023 ( Figure 8 ) and a maximum depth of erosion of 1.43 m. Sediment loss was also accurately predicted with a slight overestimate of 3% (29 m3), with only 24% of the sediment (7 m3) deposited at the back beach. Consequently, the model underestimated sediment overwash by approximately 47%. The beach face slope was also accurately predicted to be approximately 4.7 degrees. Overall, the model accurately predicted the July profile with an RMSE of 0.75 m.

      Effects of mangroves on coastal erosion (current situation)

      Two mangrove scenarios were evaluated by introducing mangroves on the berm (Scn_Man_I) and within the intertidal zone (Scn_Man_II) at three density levels: sparse, medium, and dense.

      Mangroves on the berm (current situation)

      For Scn_Man_I (mangroves on the berm), the results show a reduction in the maximum depth and volume of erosion for all three density levels ( Figure 9 ; Table 3 ). With the sparse mangrove, there was approximately a 22% reduction in the volume of erosion from 28 m3 to approximately 22 m3, with the maximum depth of erosion reducing by only 0.1 m. At a medium density of mangroves, the volume of sediment eroded reduced further to 14 m3, representing a 50% reduction in erosion and a lower maximum depth of erosion of 1 m. With dense mangroves on the beach, erosion was reduced by 53% (13.3 m3), and the maximum depth of erosion recorded was 0.9 m. This represents only a 3% reduction from the medium-density scenario.

      The response of the beach profile to mangroves on the berm (Scn_Man_I) at varying densities. The green bar represents the location where mangroves were introduced in the model.

      Summary of the effects of Mangroves on the shoreline dynamics for the current situation.

      Density Levels Scenario 2 Scenario 3
      Maximum depth of erosion (m) Volume erosion (m3) Volume accretion (m3) Maximum depth of erosion (m) Volume erosion (m3) Volume accretion (m3) Wave dissipation (W/m2)
      Sparse 1.3 21.8 5.9 1.4 27.1 6.8 6
      Medium 1.0 14.0 2.7 0.9 12.9 4.4 153
      Dense 0.9 13.3 1.7 0.2 0.9 5.9 450
      Mangroves in the intertidal zone (current situation)

      For the Scn_Man_II, mangroves were introduced in the intertidal zone at the same density levels (sparse, medium, and dense). The result again shows a reduction in erosion from sparse to dense. In this scenario, sparse mangroves offered little protection, with sediment volume erosion reducing by only 3.2% from 28 m3 to 27.1 m3 and no significant change in maximum depth of erosion. With a medium density of mangroves, the protection level already exceeds what is recorded for Scn_Man_I, with approximately 54% protection. The maximum depth of erosion was also lower at 0.9 m, comparable to the dense mangroves of Scn_Man_I. At high densities, approximately 97% protection was achieved with only 0.9 m3 of sediment volume erosion. The maximum depth of erosion also reduced significantly to 0.2 m.

      With the introduction of dense mangroves in the intertidal zone, the berm was almost completely preserved (with no significant erosion observed), and a sediment gain was recorded, indicating the trapping of sediment within the intertidal zone by the mangroves ( Figure 10 ). The wave dissipation values increase from approximately 6 W/m2 to 153 W/m2 and then 450 W/m2 for sparse, medium, and dense mangroves, respectively. Notably, the dense mangrove within the intertidal zone promoted high dissipation, leading to the deposition of sediments and the formation of a nearshore sandbar protecting the coast from high-energy waves, hence the holding of the shoreline. The results of the Scenario are summarized in Table 3 .

      The response of the beach profile to mangroves in the intertidal zone (Scn_Man_II) at varying densities. The green bar represents the location where mangroves were introduced in the model.

      Scenario results under sea level rise (2040) Baseline (without mangroves)

      For the future scenario, with a sea level rise of 0.233 m, the entire profile was inundated and lowered by up to 5 m ( Figure 9 ) without mangroves introduced. This resulted in a total sediment erosion of approximately 469 m3 with no deposition observed ( Figure 11 ).

      Baseline erosion by 2040 based on sea level rise of 0.233 m without mangroves.

      Effects of mangrove on coastal erosion (future scenario) Mangroves on the berm (future scenario)

      With mangroves on the berm, relatively higher protection levels were observed following a similar pattern as the current situation. For sparse mangroves, sediment volume erosion reduced drastically from 469 m3 for the baseline to approximately 34 m3, representing a 93% reduction in erosion. The maximum depth of erosion recorded was 1.5 m. With medium to dense mangroves, sediment erosion was reduced further to 12 m3 (for medium) and 13 m3 (for dense), representing approximately 97% reduction in erosion in both scenarios. Also, the maximum depth of erosion was reduced to 1 m in both cases ( Figure 12 ; Table 4 ).

      The effects of mangroves on the berm at various densities with sea level rise of 0.233 m by the year 2040. The green bar represents the location where mangroves were introduced in the model.

      Effects of mangroves on coastal dynamics under sea level rise.

      Density Levels Scenario 2 Scenario 3
      Maximum depth of erosion (m) Volume erosion (m3) Volume accretion (m3) Maximum depth of erosion (m) Volume erosion (m3) Volume accretion (m3) Max wave dissipation (W/m2)
      Low 1.5 34.3 7.9 4.5 406.9 0 263
      Medium 1.0 13.4 4.4 3.6 293.5 0 391
      High 1.0 12.2 3.6 0.4 2.6 5.5 506
      Mangroves in the intertidal zone (future scenario)

      In this future scenario, sparse to medium-density mangroves within the intertidal zone will offer little protection, between 13% and 37% respectively. This can be seen in the relatively high volumes of erosion, 407 m3 and 294 m3 of erosion, and a maximum depth of erosion of 4.5 m and 3.6 m, respectively ( Table 4 ). In both cases, the mangroves are completely eroded, and the whole area is inundated ( Figure 13 ).

      The effects of mangroves on beach erosion under a future sea level rise of 0.233 m by the year 2040. The green bar indicates where mangroves were introduced.

      However, at higher densities, the mangroves were able to stabilize the shoreline, reducing volume erosion from 469 m3 to 3 m3, representing a 99% reduction in erosion. The high density means higher heights, which can attenuate waves better (506 W/m2) and trap sediment, leading to a higher protective ability ( Table 4 ).

      Discussion Model calibration and performance evaluation

      The model calibration carried out shows the model was sensitive to all three parameters, namely facAs, delta, and bermslope. Similar results have been reported by previous studies, both in the area and elsewhere (Verheyen et al., 2014; Roelvink and Costas, 2019; Kombiadou et al., 2021). The best values (facAs 0.35, delta 0.4, and bermslope 0.18) were able to capture the erosion trend more accurately, with a 3% overestimation of sediment loss. However, sediment overwash was significantly underestimated (47%), which can be attributed to the underestimation of wave-runup. Other studies, such as De Beer et al. (2020); Kombiadou et al. (2021), and Roelvink et al. (2017), have established that with the surfbeat mode of XBeach, there is a general underestimation of runup, especially for steep beaches, which plays a role in berm dynamics. Overall, the model performed well, with an overall RMSE of 0.75 m. The results of this model assessment present improvements compared to previous attempts reported by Verheyen et al. (2014). The main limitation of this model (also for previous attempts) was the absence of data for the surf zone profile; hence, idealized surf zone profiles were used based on existing profiles elsewhere and initial model behavior using an equilibrium interpolation. Also, there is limited observational data on the ocean state for calibration and validation of the wave dynamics. With a future sea level rise of 0.233 m by 2040 under the worst-case scenario (SSP5-8.5), the model shows a complete inundation and erosion of the beach. This result is also consistent with other studies that indicate that even with a modest sea level rise of 0.1 m, most of the delta area will be inundated, especially along the coastline (Wiafe et al., 2013; Brempong et al., 2023; Avornyo et al., 2024).

      Impact of mangroves on coastal erosion mitigation

      With a reliable prediction of erosion from the model, What-if Scenarios (WiS) of mangroves were evaluated to assess their ability to protect the coast against coastal erosion both under current and future sea levels. The results demonstrate the ability of mangroves to protect the shoreline with a significant reduction in erosion as mangroves were introduced on the beach and in the intertidal zone. The introduction of mangroves behind the berm is considered a viable option, as there is already evidence of mangroves on the berm along pockets of beaches on the eastern coast of Ghana. For the current situation, this option led to a 53% reduction in sediment volume erosion with dense mangroves introduced on the beach. With a rise in sea level, a similar trend of protection is observed, with high-density mangroves protecting the coast by 97%. The mangroves on the berm serve as a barrier, diminishing overland flow velocities and facilitating the deposition of sediments being transported by the water. This gradually increases the berm height, eventually reducing wave overtopping and erosion. This unique ability to adapt to rising sea levels through building up the elevation has been noted by other studies, such as Krauss et al., 2014 and Mitra, 2020. This effect can also observed along pockets of beaches within the Delta where mangroves are on the berm (see Figure 2 ).

      The second scenario, which does not exist naturally along this coast, is the mangroves within the intertidal zone. However, this scenario was able to halt erosion completely with the introduction of dense mangroves for both current and future sea level rise. Generally, the higher densities of mangroves were very effective at dissipating wave energy (up to 506 W/m2), which is consistent with other studies (Bao, 2011; Spalding et al., 2014; Kamil et al., 2021). However, though sparse to medium-density mangroves can offer some level of protection for the current situation, higher densities of mangroves are required to protect the coast under rising sea levels. Overall, introducing mangroves into the intertidal zone, while not typical for the region, could serve as an enhanced protective measure when facing extreme wave water level conditions. In particular, the formation of natural sand barriers facilitated by wave dissipation is a promising outcome, creating an added layer of protection for vulnerable coastlines. However, the practicality of this scenario would need to be examined further, as the elevated salinity and high-energy wave environment could limit mangrove growth and effectiveness in this zone. Long-term studies could explore strategies to mitigate these challenges, such as pairing the mangrove growth with salt-tolerant vegetation or temporary protection measures.

      In other studies, mangroves have been demonstrated to offer shoreline protection, for example, in the case of Hurricane Harvey, where areas with mangrove cover experienced little erosion compared to other areas without mangroves (Pennings et al., 2021). Sánchez-Núñez et al. (2019) also demonstrated through field experiments, how mangroves can reduce erosion 3 to 15 times with higher wave energies. Several laboratory studies have also demonstrated how mangroves can attenuate waves, thereby reducing coastal erosion (see Amos and Akib (2023). The dense roots of mangroves help stabilize sediments, consequently diminishing sediment suspension and erosion. The above-ground roots and stems are also able to slow down flow and trap sediment, hence building the beach where they are present. This ability has been demonstrated through this assessment and can be leveraged to offer nature-based coastal protection for the eastern coast of Ghana and similar coasts.

      Integrating nature-based solutions in coastal management

      From a coastal management perspective, the results highlight the value of integrating nature-based solutions, such as mangroves, into broader coastal protection strategies. Grey infrastructure solutions, such as seawalls and groynes, though effective in localized settings, tend to disrupt sediment transport and create long-term environmental issues. By incorporating mangroves into existing management plans, especially in areas such as the Volta Delta, decision-makers can reduce the reliance on costly and environmentally disruptive “hard” engineering structures. The success of the medium to high-density mangrove scenarios in reducing erosion, for instance, suggests that coastal managers could implement mangrove plantations in areas where beach loss is currently mitigated by engineered defenses, creating a hybrid solution that maximizes both ecological and structural benefits. Furthermore, the use of WiS in this study offers significant insight into future coastal resilience planning. The ability to test different densities and placements of mangroves allows decision-makers to make informed choices that consider both ecological and hydrodynamic conditions. As coastal threats increase due to sea-level rise and storm surges, the ability to simulate these scenarios offers a proactive approach to managing risks rather than reacting to disasters. This approach can be extended to other regions facing similar challenges, using adaptive management practices that integrate scenario testing into long-term coastal planning.

      This study acknowledges the challenges imminent with implementing mangroves for coastal protection (some of which have been discussed) and requires further investigation. However, there have been attempts elsewhere to investigate these challenges. For example, elsewhere, attempts have been made to introduce temporal structures to protect mangroves until they reach a growth stage that can offer protection to the coast (Yuanita et al., 2019; Amos and Akib, 2023). They have also been effectively combined with other engineering approaches to reduce the costs and challenges (Tusinski and Jan Verhagen, 2014). This hybrid approach, combining grey infrastructure with natural buffers like mangroves, could be particularly useful in highly dynamic and at-risk coastal regions like the Volta Delta. The strategic placement of mangroves alongside groynes or revetments, for instance, could help mitigate downdrift erosion caused by engineered structures, creating a more sustainable and balanced solution. Such methods not only reduce maintenance costs but also enhance biodiversity and long-term coastal resilience. This approach can be investigated and adopted to help promote the use of mangroves as a nature-based solution for coastal erosion along the eastern coast of Ghana.

      Conclusions

      This study evaluated the effectiveness of mangroves as a Nature-based Solution (NBS) for mitigating coastal erosion in the Volta Delta, using a modelling approach to simulate erosion under different scenarios. The results demonstrate the significant protective capacity of mangroves, with a reduction in sediment volume erosion of up to 53% for the current situation and 97% for the future scenario when high-density mangroves were introduced on the berm. In the more novel scenario, where dense mangroves were placed in the intertidal zone, erosion was almost entirely mitigated, with a 97% reduction in sediment volume erosion for the current situation and a 99% reduction for the future scenario. These findings align with broader studies on the role of mangroves in coastal protection, particularly their ability to attenuate wave energy and promote sediment deposition.

      The scenario-based approach enabled a nuanced assessment of the different densities and placements of mangroves. In the scenario where mangroves were introduced behind the berm, erosion reduction was comparatively lower. Meanwhile, the intertidal placement of dense mangroves proved highly effective, not only dissipating wave energy but also facilitating the development of offshore sandbars, which further protected the coastline from wave-induced erosion. Quantitatively, the wave dissipation values increased significantly across the scenarios, from 6 W/m² with sparse mangroves to 506 W/m² in the dense mangrove scenario, showcasing the enhanced protective capability of higher-density mangrove stands. These results highlight the ability of mangroves to serve as both a physical barrier and a dynamic ecosystem that stabilizes coastal zones.

      From a coastal management perspective, the findings underscore the viability of integrating NbS into traditional coastal defense strategies. Given the escalating threats from sea-level rise and increasing storm surges, the use of mangroves offers a more sustainable and cost-effective alternative or complement to conventional “grey” infrastructure such as seawalls and groynes. Mangroves provide additional ecosystem services, such as carbon sequestration and habitat provision, which makes them a multifaceted solution for coastal protection. Coastal managers can consider implementing mangrove plantations in combination with engineered defenses to create a hybrid solution that maximizes both ecological and structural benefits. Furthermore, the use of What-if Scenarios (WiS) proved instrumental in demonstrating the potential of NBS under different environmental conditions. These scenarios allowed for a flexible and adaptive decision-making process, which is crucial in dynamic coastal environments like the Volta Delta, where erosion is exacerbated by anthropogenic activities and climate change. Scenario-based planning, coupled with quantitative model outputs, offers coastal managers a robust tool for designing interventions that are both ecologically sustainable and economically viable.

      While this study confirms the efficacy of mangroves as an NBS for coastal erosion management, several areas require further investigation to address critical needs. These include continuous monitoring of mangrove growth and coastal dynamics to understand the long-term effectiveness of mangroves as NbS; research on hybrid approaches that combine mangroves with engineered structures; and evaluating the cost-effectiveness of mangrove-based solutions compared to the traditional methods used in the area.

      Data availability statement

      The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

      Author contributions

      P-NJ-Q: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing. JS: Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Supervision, Writing – review & editing. WC: Conceptualization, Methodology, Writing – review & editing. BD: Conceptualization, Methodology, Writing – review & editing. EM: Writing – review & editing, Supervision, Project administration, Resources. KA: Writing – review & editing, Supervision, Project administration, Resources.

      Funding

      The author(s) declare that financial support was received for the research and/or publication of this article. This research is made possible through funding from GIZ MeerWissen through the MANCOGA project (Grant Agreement No. 81298219 awarded to the University of Ghana and Helmholtz-Zentrum Hereon). WC was funded by REST-COAST. This study was conducted under the MANCOGA (Ocean Decade Endorsed) project, supported by the support of MEERWISSEN and Research visit funding from the Helmholtz Information and Data Science Academy (HIDA) Trainee Network. We are also grateful to the Harmony Coast Project, sponsored by the National Geographic Society, grant number NGS-97885R-23, for providing data. WC and JS acknowledge the EU Green Deal project REST-COAST: Large-scale restoration of coastal ecosystems through rivers to sea connectivity (grant agreement 101037097 and the EU Project HORIZON-MISS-2021-OCEAN-05-01, EDITO-Model Lab: Lab Underlying models for the European Digital Twin Ocean (grant agreement 101093293).

      Conflict of interest

      The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

      The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

      Generative AI statement

      The author(s) declare that no Generative AI was used in the creation of this manuscript.

      Publisher’s note

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

      References Aagaard T. Anthony E. J. Gillies B. Laursen S. N. Sukstorf F. N. Breuning-madsen H. (2021). Geomorphology Holocene development and coastal dynamics at the Keta Sand Spit, Volta River delta, Ghana. Geomorphology 387, 107766. doi: 10.1016/j.geomorph.2021.107766 Adytia D. Husrin S. Latifah A. L. (2019). Dissipation of solitary wave due to mangrove forest: A numerical study by using non-dispersive wave model. Ilmu Kelaut 24, 4150. doi: 10.14710/ik.ijms.24.1.41-50 Amos D. Akib S. (2023). A review of coastal protection using artificial and natural countermeasures—Mangrove vegetation and polymers. Eng 4, 941953. doi: 10.3390/eng4010055 Angnuureng D. B. (2023). Bight of Benin coastal monitoring programme. Available online at: https://oceandata.ucc.edu.gh/ (Accessed December 22, 2023). Angnuureng D. B. Appeaning Addo K. Wiafe G. (2013). Impact of sea defense structures on downdrift coasts: The case of Keta in Ghana. Academia J. Environ. Sci. 1, 104121. doi: 10.15413/ajes.2013.0102 Angnuureng B. D. Brempong E. K. Almar R. Tutu G. O. Attipoe E. Klubi E. . (2024). The state of the lower Volta Delta Beaches in Ghana from field observations (Zenodo). doi: 10.5281/zenodo.10554597 Anthony E. J. Almar R. Aagaard T. (2016). Recent shoreline changes in the Volta River delta, West Africa: the roles of natural processes and human impacts. Afr J. Aquat Sci. 41, 8187. doi: 10.2989/16085914.2015.1115751 Anthony E. J. Almar R. Besset M. Reyns J. Laibi R. Ranasinghe R. . (2019). Response of the Bight of Benin (Gulf of Guinea, West Africa) coastline to anthropogenic and natural forcing, Part 2: Sources and patterns of sediment supply, sediment cells, and recent shoreline change. Cont Shelf Res. 173, 93103. doi: 10.1016/j.csr.2018.12.006 Appeaning Addo K. Nicholls R. J. Codjoe S. N. A. Abu M. (2018). A biophysical and socioeconomic review of the Volta delta, Ghana. J. Coast Res. 345, 12161226. doi: 10.2112/jcoastres-d-17-00129.1 Avornyo S. Y. Minderhoud P. S. J. Teatini P. Seeger K. Hauser L. T. Woillez M. N. . (2024). The contribution of coastal land subsidence to potential sea-level rise impact in data-sparse settings: The case of Ghana’s Volta delta. Quaternary Sci. Adv. 14, 100175. doi: 10.1016/J.QSA.2024.100175 Awuku-Sowah E. M. Graham N. A. J. Watson N. M. (2023). The contributions of mangroves to physiological health in Ghana: Insights from a qualitative study of key informants. Space Soc. 4, 100137. doi: 10.1016/j.wss.2023.100137 Bao T. Q. (2011). Effect of mangrove forest structures on wave attenuation in coastal Vietnam. Oceanologia 53, 807818. doi: 10.5697/oc.53-3.807 Bolle A. das Neves L. Rooseleer J. (2015). Coastal protection for Ada, Ghana: A case study. Proc. Institution Civil Engineers: Maritime Eng. 168, 19. doi: 10.1680/jmaen.15.00013 Bollen M. Trouw K. Lerouge F. Gruwez V. Bolle A. Hoffman B. . (2011). Design of a coastal protection scheme for Ada at the Volta River mouth (Ghana). Coastal Eng. Proc. 1, 36. doi: 10.9753/icce.v32.management.36 Brempong E. K. Almar R. Angnuureng D. B. Mattah P. A. D. Avornyo S. Y. Jayson-Quashigah P. N. . (2023). Future flooding of the Volta Delta caused by sea level rise and land subsidence. J. Coast Conserv 27, 24. doi: 10.1007/s11852-023-00952-0 Burger B. (2005). Wave Attenuation in Mangrove Forests: Numerical modelling of wave attenuation by implementation of a physical description of vegetation in SWAN (Delft: Delft University of Technology). Charoenlerkthawin W. Bidorn K. Burnett W. C. Sasaki J. Panneerselvam B. Bidorn B. (2022). Effectiveness of grey and green engineered solutions for protecting the low-lying muddy coast of the Chao Phraya Delta, Thailand. Sci. Rep. 12 (1), 20448. doi: 10.1038/s41598-022-24842-x Chen W. L. Muller P. Grabowski R. C. Dodd N. (2022). Green nourishment: an innovative nature-based solution for coastal erosion. Front. Mar Sci. 8. doi: 10.3389/fmars.2021.814589 Chen W. Staneva J. Jacob B. Sánchez-Artús X. Wurpts A. (2024). What-if nature-based storm buffers on mitigating coastal erosion. Sci. Total Environ. 928, 172247. doi: 10.1016/J.SCITOTENV.2024.172247 Dean R. G. (1991). Equilibrium beach profiles: characteristics and applications. J. Coast. Res. 7 (1), 5384. Available online at: http://www.jstor.org/stable/4297805 De Beer A. F. Mccall R. T. Long J. W. Tissier M. F. S. Reniers A. J. H. M. (2020). Simulating wave runup on an intermediate-reflective beach using a wave-resolving and a wave-averaged version of XBeach. Coast. Eng. 163, 103788. doi: 10.1016/j.coastaleng.2020.103788 Ericson J. P. Vörösmarty C. J. Dingman S. L. Ward L. G. Meybeck M. (2006). Effective sea-level rise and deltas: Causes of change and human dimension implications. Glob Planet Change 50, 6382. doi: 10.1016/j.gloplacha.2005.07.004 Ghana Ports (2025). Port-log - tides. Available online at: https://Ghanaports.port-log.net/live/Display.php?Dataset=1 (Accessed February 28, 2025). Giardino A. Schrijvershof R. Nederhoff C. M. de Vroeg H. Brière C. Tonnon P. K. . (2018). A quantitative assessment of human interventions and climate change on the West African sediment budget. Ocean Coast Manag 156, 249265. doi: 10.1016/j.ocecoaman.2017.11.008 Gruwez V. Verheyen B. Wauters P. Bolle A. (2014). 2DH Morphodynamic Time-Dependent Hindcast Modelling of a Groyne System in Ghana. In Proceedings of the 11th International Conference on Hydroscience & Engineering, eds. Lehfeldt R. Kopmann R. (Hamburg: Karlsruhe: Bundesanstalt für Wasserbau). Hagedoorn L. C. Appeaning Addo K. Koetse M. J. Kinney K. van Beukering P. J. H. (2021). Angry waves that eat the coast: An economic analysis of nature-based and engineering solutions to coastal erosion. Ocean Coast. Manag. 214, 105945. doi: 10.1016/j.ocecoaman.2021.105945 Hersbach H. Bell B. Berrisford P. Biavati G. Horányi A. Muñoz Sabater J. . (2023). ERA5 hourly data on single levels from 1940 to present (Copernicus Climate Change Service (C3S) Climate Data Store (CDS). doi: 10.24381/cds.adbb2d47 IHO-IOC (2019). The IHO-IOC GEBCO cook book (Monaco: IHO Publication B-11). Available online at: http://www.gebco.net/. IPCC . (2021). Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change. Eds. Masson-Delmotte V. Zhai P. Pirani A. Connors S. L. Péan C. Berger S. . (Cambridge, United Kingdom and New York, NY, USA: Cambridge University Press). doi: 10.1017/9781009157896 James R. K. Lynch A. Herman P. M. J. van Katwijk M. M. van Tussenbroek B. I. Dijkstra H. A. . (2021). Tropical biogeomorphic seagrass landscapes for coastal protection: persistence and wave attenuation during major storm events. Ecosystems 24, 301318. doi: 10.1007/s10021-020-00519-2 Jayson-Quashigah P.-N. Addo K. A. Kodzo K. S. (2013). Medium resolution satellite imagery as a tool for monitoring shoreline change. Case study of the Eastern coast of Ghana. J. Coast Res. (65), 511516. doi: 10.2112/SI65-087.1 Jayson-Quashigah P.-N. Appeaning Addo K. Amisigo B. Wiafe G. (2019). Assessment of short-term beach sediment change in the Volta Delta coast in Ghana using data from Unmanned Aerial Vehicles (Drone). Ocean Coast Manag 182, 104952. doi: 10.1016/j.ocecoaman.2019.104952 Jayson-Quashigah P.-N. Appeaning Addo K. Wiafe G. Amisigo B. A. Brempong E. K. Kay S. . (2021). Wave dynamics and shoreline evolution in deltas: A case study of sandy coasts in the Volta delta of Ghana. Interpretation 9, SH99SH113. doi: 10.1190/int-2021-0028.1 Kamil E. A. Takaijudin H. Hashim A. M. (2021). Mangroves as coastal bio-shield: A review of mangroves performance in wave attenuation. Civil Eng. J. (Iran) 7, 19641981. doi: 10.28991/cej-2021-03091772 Khanh Phan L. (2019). Wave attenuation in coastal mangroves: Mangrove squeeze in the Mekong Delta (Vietnam: Delft University of Technology). doi: 10.4233/uuid:9397d964-1674-4838-a13a-504742dba55e Kombiadou K. Costas S. Roelvink D. (2021). Simulating destructive and constructive morphodynamic processes in steep beaches. J. Mar Sci. Eng 9, 119. doi: 10.3390/jmse9010086 Kombiadou K. Costas S. Roelvink D. McCall R. (2020). Post-storm recuperation as a stepping-stone towards long-term integrated modelling in steep beaches (EGU General Assembly). doi: 10.5194/egusphere-egu2020-17470 Krauss K. W. McKee K. L. Lovelock C. E. Cahoon D. R. Saintilan N. Reef R. . (2014). How mangrove forests adjust to rising sea level. New Phytol. 202, 1934. doi: 10.1111/nph.12605 Lopez-Arias F. Maza M. Calleja F. Govaere G. Lara J. L. (2024). Integrated drag coefficient formula for estimating the wave attenuation capacity of Rhizophora sp. mangrove forests. Front. Mar Sci. 11. doi: 10.3389/fmars.2024.1383368 MANCOGA (2024). Mangroves as nature-based solution for coastal hazards in eastern Ghana (MANCOGA). Available online at: https://mancoga.com/. Mann T. Serwa A. Rovere A. Casella E. Appeaning-Addo K. Jayson-Quashigah P.-N. . (2023). Multi-decadal shoreline changes in Eastern Ghana—natural dynamics versus human interventions. Geo-Marine Lett. 43, 17. doi: 10.1007/s00367-023-00758-x Mentaschi L. Vousdoukas M. I. Pekel J. F. Voukouvalas E. Feyen L. (2018). Global long-term observations of coastal erosion and accretion. Sci. Rep. 8, 111. doi: 10.1038/s41598-018-30904-w Mitra A. (2020). “Mangroves: A potential vegetation against sea level rise,” in Mangrove forests in India: exploring ecosystem services (Springer International Publishing, Cham), 157187. doi: 10.1007/978-3-030-20595-9_6 Nairn R. B. Dibajnia M. (2004). Design and construction of a large headland system, Keta Sea Defence Project, West Africa. J. Coast Res. 33 (S1), 294314. Available online at: http://www.jstor.org/stable/25736261. Nunoo F. K. E. Agyekumhene A. (2022). Mangrove degradation and management practices along the coast of Ghana. Agric. Sci. 13, 10571079. doi: 10.4236/as.2022.1310065 Ofori S. A. Asante F. Ama T. Boateng B. Dahdouh-Guebas F. (2023). The composition, distribution, and socio-economic dimensions of Ghana’s mangrove ecosystems. J. Environ. Manage 345, 3014797. doi: 10.1016/j.jenvman.2023.118622 Pender D. Karunarathna H. (2013). A statistical-process based approach for modelling beach profile variability. Coastal Eng. 81, 1929. doi: 10.1016/j.coastaleng.2013.06.006 Pennings S. C. Glazner R. M. Hughes Z. J. Kominoski J. S. Armitage A. R. (2021). Effects of mangrove cover on coastal erosion during a hurricane in Texas, USA. Ecology 102 (4), e03309. doi: 10.1002/ecy.3309 Roelvink D. Costas S. (2019). Coupling nearshore and aeolian processes: XBeach and duna process-based models. Environ. Modelling Software 115, 98112. doi: 10.1016/j.envsoft.2019.02.010 Roelvink D. Mccall R. Mehvar S. Nederhoff K. Dastgheib A. (2018). Improving predictions of swash dynamics in XBeach: The role of groupiness and incident-band runup. Coast. Eng. 134, 103123. doi: 10.1016/j.coastaleng.2017.07.004 Roelvink D. Reniers A. van Dongeren A. van Thiel de Vries J. McCall R. Lescinski J. (2009). Modelling storm impacts on beaches, dunes and barrier islands. Coastal Eng. 56, 11331152. doi: 10.1016/J.COASTALENG.2009.08.006 Roest L. W. M. (2018). The coastal system of the Volta delta, Ghana Opportunities and strategies for development. TU Delft Delta Infrastructures and Mobility Initiative (DIMI). Available online at: https://pure.tudelft.nl/ws/files/37464456/Roest_2018_The_coastal_system_of_the_Volta_delta.pdf. Sánchez-Núñez D. A. Bernal G. Mancera Pineda J. E. (2019). The relative role of mangroves on wave erosion mitigation and sediment properties. Estuaries Coasts 42, 21242138. doi: 10.1007/s12237-019-00628-9 Scott T. McCarroll R. J. Masselink G. Castelle B. Dodet G. Saulter A. . (2021). Role of atmospheric indices in describing inshore directional wave climate in the United Kingdom and Ireland. Earths Future 9 (5). doi: 10.1029/2020EF001625 Scown M. W. Dunn F. E. Dekker S. C. van Vuuren D. P. Karabil S. Sutanudjaja E. H. . (2023). Global change scenarios in coastal river deltas and their sustainable development implications. Global Environ. Change 82, 102736. doi: 10.1016/j.gloenvcha.2023.102736 Singhvi A. Luijendijk A. P. van Oudenhoven A. P. E. (2022). The grey–green spectrum: A review of coastal protection interventions. J. Environ. Manage 311, 114824. doi: 10.1016/j.jenvman.2022.114824 Soci C. Hersbach H. Simmons A. Poli P. Bell B. Berrisford P. . (2024). The ERA5 global reanalysis from 1940 to 2022. Q. J. R. Meteorological Soc. 150 (764), 40144048. doi: 10.1002/qj.4803 Spalding M. McIvor A. Tonneijck F. H. Tol S. van Eijk P. (2014). Mangroves for coastal defence: Guidelines for coastal managers & policy makers. The Nature Conservancy. Available online at: www.nature.org (Accessed March 3, 2025). Tessler Z. D. Vörösmarty C. J. Grossberg M. Gladkova I. Aizenman H. Syvitski J. P. M. . (2015). Profiling risk and sustainability in coastal deltas of the world. Sci. (1979) 349, 638643. doi: 10.1126/science.aab3574 Thao N. D. Kiet N. Anh L. T. Huynh P. T. D. Hiep L. D. Trang D. D. T. (2023). “Wave attenuation in coastal mangroves: the case of mekong delta, Vietnam,” in IOP Conference Series: Earth and Environmental Science (IOP Publishing) 1226, 012026. doi: 10.1088/1755-1315/1226/1/012026 Tusinski A. Jan Verhagen H. (2014). Coastal Engineering Proceedings. 1 (34), management.45. doi: 10.9753/icce.v34.management.45 van Hespen R. Hu Z. Borsje B. De Dominicis M. Friess D. A. Jevrejeva S. . (2023). Mangrove forests as a nature-based solution for coastal flood protection: Biophysical and ecological considerations. Water Sci. Eng. 16, 113. doi: 10.1016/j.wse.2022.10.004 van Zelst V. T. M. Dijkstra J. T. van Wesenbeeck B. K. Eilander D. Morris E. P. Winsemius H. C. . (2021). Cutting the costs of coastal protection by integrating vegetation in flood defences. Nat. Commun. 12, 6533. doi: 10.1038/s41467-021-26887-4 Verheyen B. Gruwez V. Zimmermann N. Bolle A. Wauters P. (2014). “Medium term time-dependent morphodynamic modelling of beach profile evolution in ada, Ghana,” in 11th International Conference on Hydroscience & Engineering. 701708 (Hamburg). Available online at: http://vzb.baw.de/e-medien/iche-2014/PDF/07SedimentTransportandMorphodynamics/07_23.pdf (Accessed September 28-October 2, 2014). Vousdoukas M. I. Mentaschi L. Voukouvalas E. Verlaan M. Jevrejeva S. Jackson L. P. . (2018). Global probabilistic projections of extreme sea levels show intensification of coastal flood hazard. Nat. Commun. 9, 112. doi: 10.1038/s41467-018-04692-w Vousdoukas M. I. Ranasinghe R. Mentaschi L. Plomaritis T. A. Athanasiou P. Luijendijk A. . (2020). Sandy coastlines under threat of erosion. Nat. Clim Chang 10, 260263. doi: 10.1038/s41558-020-0697-0 Vousdoukas M. I. Almeida L. P. M. Ferreira Ó. (2012). Beach erosion and recovery during consecutive storms at a steep-sloping, meso-tidal beach. Earth Surf Process Landf 37, 583593. doi: 10.1002/esp.2264 Vuik V. Jonkman S. N. Borsje B. W. Suzuki T. (2016). Nature-based flood protection: The efficiency of vegetated foreshores for reducing wave loads on coastal dikes. Coastal Eng. 116, 4256. doi: 10.1016/j.coastaleng.2016.06.001 Wegman C. Wilms T. Angnuureng D. Lambregts P. Ritsema J. (2023). Baseline study West African Case Studies for Coastal Nature Based Solutions. Netherlands Enterprise Agency. Available online at: https://www.rvo.nl/sites/default/files/2024-02/Report_BaselinestudyWestAfricanCaseStudiesforCoastalNBS_final.pdf. Wiafe G. Boateng I. Appeaning-Addo K. Jayson-Quashigah P.-N. Ababio S. D. Sowah L. (2013). Handbook of coastal processes and management in Ghana. 1st Edn (Gloucester, UK: The Choir Press). XBeach Team (2023). XBeach manual: Release XBeach 1.24.6057 Halloween. Available online at: https://svn.oss.deltares.nl/repos/xbeach/trunk (Accessed June 4, 2024). Yoshikai M. Nakamura T. Bautista D. M. Herrera E. C. Baloloy A. Suwa R. . (2022). Field measurement and prediction of drag in a planted rhizophora mangrove forest. J. Geophys Res. Oceans 127 (11). doi: 10.1029/2021JC018320 Yuanita N. Kurniawan A. Setiawan H. Hasan F. Khasanah M. (2019). Physical model of natural coastal protection system: Wave transmission over mangrove seedling trees. J. Coast Res. 91, 176180. doi: 10.2112/SI91-036.1
      ‘Oh, my dear Thomas, you haven’t heard the terrible news then?’ she said. ‘I thought you would be sure to have seen it placarded somewhere. Alice went straight to her room, and I haven’t seen her since, though I repeatedly knocked at the door, which she has locked on the inside, and I’m sure it’s most unnatural of her not to let her own mother comfort her. It all happened in a moment: I have always said those great motor-cars shouldn’t be allowed to career about the streets, especially when they are all paved with cobbles as they are at Easton Haven, which are{331} so slippery when it’s wet. He slipped, and it went over him in a moment.’ My thanks were few and awkward, for there still hung to the missive a basting thread, and it was as warm as a nestling bird. I bent low--everybody was emotional in those days--kissed the fragrant thing, thrust it into my bosom, and blushed worse than Camille. "What, the Corner House victim? Is that really a fact?" "My dear child, I don't look upon it in that light at all. The child gave our picturesque friend a certain distinction--'My husband is dead, and this is my only child,' and all that sort of thing. It pays in society." leave them on the steps of a foundling asylum in order to insure [See larger version] Interoffice guff says you're planning definite moves on your own, J. O., and against some opposition. Is the Colonel so poor or so grasping—or what? Albert could not speak, for he felt as if his brains and teeth were rattling about inside his head. The rest of[Pg 188] the family hunched together by the door, the boys gaping idiotically, the girls in tears. "Now you're married." The host was called in, and unlocked a drawer in which they were deposited. The galleyman, with visible reluctance, arrayed himself in the garments, and he was observed to shudder more than once during the investiture of the dead man's apparel. HoME香京julia种子在线播放 ENTER NUMBET 0016www.gcqjwc.com.cn
      gcxpfn.com.cn
      www.goseven.net.cn
      www.protestant.com.cn
      tgergz.com.cn
      www.shenyanyi.com.cn
      tflv7.net.cn
      www.qhll.com.cn
      www.weconz.com.cn
      qmesub.com.cn
      处女被大鸡巴操 强奸乱伦小说图片 俄罗斯美女爱爱图 调教强奸学生 亚洲女的穴 夜来香图片大全 美女性强奸电影 手机版色中阁 男性人体艺术素描图 16p成人 欧美性爱360 电影区 亚洲电影 欧美电影 经典三级 偷拍自拍 动漫电影 乱伦电影 变态另类 全部电 类似狠狠鲁的网站 黑吊操白逼图片 韩国黄片种子下载 操逼逼逼逼逼 人妻 小说 p 偷拍10幼女自慰 极品淫水很多 黄色做i爱 日本女人人体电影快播看 大福国小 我爱肏屄美女 mmcrwcom 欧美多人性交图片 肥臀乱伦老头舔阴帝 d09a4343000019c5 西欧人体艺术b xxoo激情短片 未成年人的 插泰国人夭图片 第770弾み1 24p 日本美女性 交动态 eee色播 yantasythunder 操无毛少女屄 亚洲图片你懂的女人 鸡巴插姨娘 特级黄 色大片播 左耳影音先锋 冢本友希全集 日本人体艺术绿色 我爱被舔逼 内射 幼 美阴图 喷水妹子高潮迭起 和后妈 操逼 美女吞鸡巴 鸭个自慰 中国女裸名单 操逼肥臀出水换妻 色站裸体义术 中国行上的漏毛美女叫什么 亚洲妹性交图 欧美美女人裸体人艺照 成人色妹妹直播 WWW_JXCT_COM r日本女人性淫乱 大胆人艺体艺图片 女同接吻av 碰碰哥免费自拍打炮 艳舞写真duppid1 88电影街拍视频 日本自拍做爱qvod 实拍美女性爱组图 少女高清av 浙江真实乱伦迅雷 台湾luanlunxiaoshuo 洛克王国宠物排行榜 皇瑟电影yy频道大全 红孩儿连连看 阴毛摄影 大胆美女写真人体艺术摄影 和风骚三个媳妇在家做爱 性爱办公室高清 18p2p木耳 大波撸影音 大鸡巴插嫩穴小说 一剧不超两个黑人 阿姨诱惑我快播 幼香阁千叶县小学生 少女妇女被狗强奸 曰人体妹妹 十二岁性感幼女 超级乱伦qvod 97爱蜜桃ccc336 日本淫妇阴液 av海量资源999 凤凰影视成仁 辰溪四中艳照门照片 先锋模特裸体展示影片 成人片免费看 自拍百度云 肥白老妇女 女爱人体图片 妈妈一女穴 星野美夏 日本少女dachidu 妹子私处人体图片 yinmindahuitang 舔无毛逼影片快播 田莹疑的裸体照片 三级电影影音先锋02222 妻子被外国老头操 观月雏乃泥鳅 韩国成人偷拍自拍图片 强奸5一9岁幼女小说 汤姆影院av图片 妹妹人艺体图 美女大驱 和女友做爱图片自拍p 绫川まどか在线先锋 那么嫩的逼很少见了 小女孩做爱 处女好逼连连看图图 性感美女在家做爱 近距离抽插骚逼逼 黑屌肏金毛屄 日韩av美少女 看喝尿尿小姐日逼色色色网图片 欧美肛交新视频 美女吃逼逼 av30线上免费 伊人在线三级经典 新视觉影院t6090影院 最新淫色电影网址 天龙影院远古手机版 搞老太影院 插进美女的大屁股里 私人影院加盟费用 www258dd 求一部电影里面有一个二猛哥 深肛交 日本萌妹子人体艺术写真图片 插入屄眼 美女的木奶 中文字幕黄色网址影视先锋 九号女神裸 和骚人妻偷情 和潘晓婷做爱 国模大尺度蜜桃 欧美大逼50p 西西人体成人 李宗瑞继母做爱原图物处理 nianhuawang 男鸡巴的视屏 � 97免费色伦电影 好色网成人 大姨子先锋 淫荡巨乳美女教师妈妈 性nuexiaoshuo WWW36YYYCOM 长春继续给力进屋就操小女儿套干破内射对白淫荡 农夫激情社区 日韩无码bt 欧美美女手掰嫩穴图片 日本援交偷拍自拍 入侵者日本在线播放 亚洲白虎偷拍自拍 常州高见泽日屄 寂寞少妇自卫视频 人体露逼图片 多毛外国老太 变态乱轮手机在线 淫荡妈妈和儿子操逼 伦理片大奶少女 看片神器最新登入地址sqvheqi345com账号群 麻美学姐无头 圣诞老人射小妞和强奸小妞动话片 亚洲AV女老师 先锋影音欧美成人资源 33344iucoom zV天堂电影网 宾馆美女打炮视频 色五月丁香五月magnet 嫂子淫乱小说 张歆艺的老公 吃奶男人视频在线播放 欧美色图男女乱伦 avtt2014ccvom 性插色欲香影院 青青草撸死你青青草 99热久久第一时间 激情套图卡通动漫 幼女裸聊做爱口交 日本女人被强奸乱伦 草榴社区快播 2kkk正在播放兽骑 啊不要人家小穴都湿了 www猎奇影视 A片www245vvcomwwwchnrwhmhzcn 搜索宜春院av wwwsee78co 逼奶鸡巴插 好吊日AV在线视频19gancom 熟女伦乱图片小说 日本免费av无码片在线开苞 鲁大妈撸到爆 裸聊官网 德国熟女xxx 新不夜城论坛首页手机 女虐男网址 男女做爱视频华为网盘 激情午夜天亚洲色图 内裤哥mangent 吉沢明歩制服丝袜WWWHHH710COM 屌逼在线试看 人体艺体阿娇艳照 推荐一个可以免费看片的网站如果被QQ拦截请复制链接在其它浏览器打开xxxyyy5comintr2a2cb551573a2b2e 欧美360精品粉红鲍鱼 教师调教第一页 聚美屋精品图 中韩淫乱群交 俄罗斯撸撸片 把鸡巴插进小姨子的阴道 干干AV成人网 aolasoohpnbcn www84ytom 高清大量潮喷www27dyycom 宝贝开心成人 freefronvideos人母 嫩穴成人网gggg29com 逼着舅妈给我口交肛交彩漫画 欧美色色aV88wwwgangguanscom 老太太操逼自拍视频 777亚洲手机在线播放 有没有夫妻3p小说 色列漫画淫女 午间色站导航 欧美成人处女色大图 童颜巨乳亚洲综合 桃色性欲草 色眯眯射逼 无码中文字幕塞外青楼这是一个 狂日美女老师人妻 爱碰网官网 亚洲图片雅蠛蝶 快播35怎么搜片 2000XXXX电影 新谷露性家庭影院 深深候dvd播放 幼齿用英语怎么说 不雅伦理无需播放器 国外淫荡图片 国外网站幼幼嫩网址 成年人就去色色视频快播 我鲁日日鲁老老老我爱 caoshaonvbi 人体艺术avav 性感性色导航 韩国黄色哥来嫖网站 成人网站美逼 淫荡熟妇自拍 欧美色惰图片 北京空姐透明照 狼堡免费av视频 www776eom 亚洲无码av欧美天堂网男人天堂 欧美激情爆操 a片kk266co 色尼姑成人极速在线视频 国语家庭系列 蒋雯雯 越南伦理 色CC伦理影院手机版 99jbbcom 大鸡巴舅妈 国产偷拍自拍淫荡对话视频 少妇春梦射精 开心激动网 自拍偷牌成人 色桃隐 撸狗网性交视频 淫荡的三位老师 伦理电影wwwqiuxia6commqiuxia6com 怡春院分站 丝袜超短裙露脸迅雷下载 色制服电影院 97超碰好吊色男人 yy6080理论在线宅男日韩福利大全 大嫂丝袜 500人群交手机在线 5sav 偷拍熟女吧 口述我和妹妹的欲望 50p电脑版 wwwavtttcon 3p3com 伦理无码片在线看 欧美成人电影图片岛国性爱伦理电影 先锋影音AV成人欧美 我爱好色 淫电影网 WWW19MMCOM 玛丽罗斯3d同人动画h在线看 动漫女孩裸体 超级丝袜美腿乱伦 1919gogo欣赏 大色逼淫色 www就是撸 激情文学网好骚 A级黄片免费 xedd5com 国内的b是黑的 快播美国成年人片黄 av高跟丝袜视频 上原保奈美巨乳女教师在线观看 校园春色都市激情fefegancom 偷窥自拍XXOO 搜索看马操美女 人本女优视频 日日吧淫淫 人妻巨乳影院 美国女子性爱学校 大肥屁股重口味 啪啪啪啊啊啊不要 操碰 japanfreevideoshome国产 亚州淫荡老熟女人体 伦奸毛片免费在线看 天天影视se 樱桃做爱视频 亚卅av在线视频 x奸小说下载 亚洲色图图片在线 217av天堂网 东方在线撸撸-百度 幼幼丝袜集 灰姑娘的姐姐 青青草在线视频观看对华 86papa路con 亚洲1AV 综合图片2区亚洲 美国美女大逼电影 010插插av成人网站 www色comwww821kxwcom 播乐子成人网免费视频在线观看 大炮撸在线影院 ,www4KkKcom 野花鲁最近30部 wwwCC213wapwww2233ww2download 三客优最新地址 母亲让儿子爽的无码视频 全国黄色片子 欧美色图美国十次 超碰在线直播 性感妖娆操 亚洲肉感熟女色图 a片A毛片管看视频 8vaa褋芯屑 333kk 川岛和津实视频 在线母子乱伦对白 妹妹肥逼五月 亚洲美女自拍 老婆在我面前小说 韩国空姐堪比情趣内衣 干小姐综合 淫妻色五月 添骚穴 WM62COM 23456影视播放器 成人午夜剧场 尼姑福利网 AV区亚洲AV欧美AV512qucomwwwc5508com 经典欧美骚妇 震动棒露出 日韩丝袜美臀巨乳在线 av无限吧看 就去干少妇 色艺无间正面是哪集 校园春色我和老师做爱 漫画夜色 天海丽白色吊带 黄色淫荡性虐小说 午夜高清播放器 文20岁女性荫道口图片 热国产热无码热有码 2015小明发布看看算你色 百度云播影视 美女肏屄屄乱轮小说 家族舔阴AV影片 邪恶在线av有码 父女之交 关于处女破处的三级片 极品护士91在线 欧美虐待女人视频的网站 享受老太太的丝袜 aaazhibuo 8dfvodcom成人 真实自拍足交 群交男女猛插逼 妓女爱爱动态 lin35com是什么网站 abp159 亚洲色图偷拍自拍乱伦熟女抠逼自慰 朝国三级篇 淫三国幻想 免费的av小电影网站 日本阿v视频免费按摩师 av750c0m 黄色片操一下 巨乳少女车震在线观看 操逼 免费 囗述情感一乱伦岳母和女婿 WWW_FAMITSU_COM 偷拍中国少妇在公车被操视频 花也真衣论理电影 大鸡鸡插p洞 新片欧美十八岁美少 进击的巨人神thunderftp 西方美女15p 深圳哪里易找到老女人玩视频 在线成人有声小说 365rrr 女尿图片 我和淫荡的小姨做爱 � 做爱技术体照 淫妇性爱 大学生私拍b 第四射狠狠射小说 色中色成人av社区 和小姨子乱伦肛交 wwwppp62com 俄罗斯巨乳人体艺术 骚逼阿娇 汤芳人体图片大胆 大胆人体艺术bb私处 性感大胸骚货 哪个网站幼女的片多 日本美女本子把 色 五月天 婷婷 快播 美女 美穴艺术 色百合电影导航 大鸡巴用力 孙悟空操美少女战士 狠狠撸美女手掰穴图片 古代女子与兽类交 沙耶香套图 激情成人网区 暴风影音av播放 动漫女孩怎么插第3个 mmmpp44 黑木麻衣无码ed2k 淫荡学姐少妇 乱伦操少女屄 高中性爱故事 骚妹妹爱爱图网 韩国模特剪长发 大鸡巴把我逼日了 中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片 大胆女人下体艺术图片 789sss 影音先锋在线国内情侣野外性事自拍普通话对白 群撸图库 闪现君打阿乐 ady 小说 插入表妹嫩穴小说 推荐成人资源 网络播放器 成人台 149大胆人体艺术 大屌图片 骚美女成人av 春暖花开春色性吧 女亭婷五月 我上了同桌的姐姐 恋夜秀场主播自慰视频 yzppp 屄茎 操屄女图 美女鲍鱼大特写 淫乱的日本人妻山口玲子 偷拍射精图 性感美女人体艺木图片 种马小说完本 免费电影院 骑士福利导航导航网站 骚老婆足交 国产性爱一级电影 欧美免费成人花花性都 欧美大肥妞性爱视频 家庭乱伦网站快播 偷拍自拍国产毛片 金发美女也用大吊来开包 缔D杏那 yentiyishu人体艺术ytys WWWUUKKMCOM 女人露奶 � 苍井空露逼 老荡妇高跟丝袜足交 偷偷和女友的朋友做爱迅雷 做爱七十二尺 朱丹人体合成 麻腾由纪妃 帅哥撸播种子图 鸡巴插逼动态图片 羙国十次啦中文 WWW137AVCOM 神斗片欧美版华语 有气质女人人休艺术 由美老师放屁电影 欧美女人肉肏图片 白虎种子快播 国产自拍90后女孩 美女在床上疯狂嫩b 饭岛爱最后之作 幼幼强奸摸奶 色97成人动漫 两性性爱打鸡巴插逼 新视觉影院4080青苹果影院 嗯好爽插死我了 阴口艺术照 李宗瑞电影qvod38 爆操舅母 亚洲色图七七影院 被大鸡巴操菊花 怡红院肿么了 成人极品影院删除 欧美性爱大图色图强奸乱 欧美女子与狗随便性交 苍井空的bt种子无码 熟女乱伦长篇小说 大色虫 兽交幼女影音先锋播放 44aad be0ca93900121f9b 先锋天耗ばさ无码 欧毛毛女三级黄色片图 干女人黑木耳照 日本美女少妇嫩逼人体艺术 sesechangchang 色屄屄网 久久撸app下载 色图色噜 美女鸡巴大奶 好吊日在线视频在线观看 透明丝袜脚偷拍自拍 中山怡红院菜单 wcwwwcom下载 骑嫂子 亚洲大色妣 成人故事365ahnet 丝袜家庭教mp4 幼交肛交 妹妹撸撸大妈 日本毛爽 caoprom超碰在email 关于中国古代偷窥的黄片 第一会所老熟女下载 wwwhuangsecome 狼人干综合新地址HD播放 变态儿子强奸乱伦图 强奸电影名字 2wwwer37com 日本毛片基地一亚洲AVmzddcxcn 暗黑圣经仙桃影院 37tpcocn 持月真由xfplay 好吊日在线视频三级网 我爱背入李丽珍 电影师傅床戏在线观看 96插妹妹sexsex88com 豪放家庭在线播放 桃花宝典极夜著豆瓜网 安卓系统播放神器 美美网丝袜诱惑 人人干全免费视频xulawyercn av无插件一本道 全国色五月 操逼电影小说网 good在线wwwyuyuelvcom www18avmmd 撸波波影视无插件 伊人幼女成人电影 会看射的图片 小明插看看 全裸美女扒开粉嫩b 国人自拍性交网站 萝莉白丝足交本子 七草ちとせ巨乳视频 摇摇晃晃的成人电影 兰桂坊成社人区小说www68kqcom 舔阴论坛 久撸客一撸客色国内外成人激情在线 明星门 欧美大胆嫩肉穴爽大片 www牛逼插 性吧星云 少妇性奴的屁眼 人体艺术大胆mscbaidu1imgcn 最新久久色色成人版 l女同在线 小泽玛利亚高潮图片搜索 女性裸b图 肛交bt种子 最热门有声小说 人间添春色 春色猜谜字 樱井莉亚钢管舞视频 小泽玛利亚直美6p 能用的h网 还能看的h网 bl动漫h网 开心五月激 东京热401 男色女色第四色酒色网 怎么下载黄色小说 黄色小说小栽 和谐图城 乐乐影院 色哥导航 特色导航 依依社区 爱窝窝在线 色狼谷成人 91porn 包要你射电影 色色3A丝袜 丝袜妹妹淫网 爱色导航(荐) 好男人激情影院 坏哥哥 第七色 色久久 人格分裂 急先锋 撸撸射中文网 第一会所综合社区 91影院老师机 东方成人激情 怼莪影院吹潮 老鸭窝伊人无码不卡无码一本道 av女柳晶电影 91天生爱风流作品 深爱激情小说私房婷婷网 擼奶av 567pao 里番3d一家人野外 上原在线电影 水岛津实透明丝袜 1314酒色 网旧网俺也去 0855影院 在线无码私人影院 搜索 国产自拍 神马dy888午夜伦理达达兔 农民工黄晓婷 日韩裸体黑丝御姐 屈臣氏的燕窝面膜怎么样つぼみ晶エリーの早漏チ○ポ强化合宿 老熟女人性视频 影音先锋 三上悠亚ol 妹妹影院福利片 hhhhhhhhsxo 午夜天堂热的国产 强奸剧场 全裸香蕉视频无码 亚欧伦理视频 秋霞为什么给封了 日本在线视频空天使 日韩成人aⅴ在线 日本日屌日屄导航视频 在线福利视频 日本推油无码av magnet 在线免费视频 樱井梨吮东 日本一本道在线无码DVD 日本性感诱惑美女做爱阴道流水视频 日本一级av 汤姆avtom在线视频 台湾佬中文娱乐线20 阿v播播下载 橙色影院 奴隶少女护士cg视频 汤姆在线影院无码 偷拍宾馆 业面紧急生级访问 色和尚有线 厕所偷拍一族 av女l 公交色狼优酷视频 裸体视频AV 人与兽肉肉网 董美香ol 花井美纱链接 magnet 西瓜影音 亚洲 自拍 日韩女优欧美激情偷拍自拍 亚洲成年人免费视频 荷兰免费成人电影 深喉呕吐XXⅩX 操石榴在线视频 天天色成人免费视频 314hu四虎 涩久免费视频在线观看 成人电影迅雷下载 能看见整个奶子的香蕉影院 水菜丽百度影音 gwaz079百度云 噜死你们资源站 主播走光视频合集迅雷下载 thumbzilla jappen 精品Av 古川伊织star598在线 假面女皇vip在线视频播放 国产自拍迷情校园 啪啪啪公寓漫画 日本阿AV 黄色手机电影 欧美在线Av影院 华裔电击女神91在线 亚洲欧美专区 1日本1000部免费视频 开放90后 波多野结衣 东方 影院av 页面升级紧急访问每天正常更新 4438Xchengeren 老炮色 a k福利电影 色欲影视色天天视频 高老庄aV 259LUXU-683 magnet 手机在线电影 国产区 欧美激情人人操网 国产 偷拍 直播 日韩 国内外激情在线视频网给 站长统计一本道人妻 光棍影院被封 紫竹铃取汁 ftp 狂插空姐嫩 xfplay 丈夫面前 穿靴子伪街 XXOO视频在线免费 大香蕉道久在线播放 电棒漏电嗨过头 充气娃能看下毛和洞吗 夫妻牲交 福利云点墦 yukun瑟妃 疯狂交换女友 国产自拍26页 腐女资源 百度云 日本DVD高清无码视频 偷拍,自拍AV伦理电影 A片小视频福利站。 大奶肥婆自拍偷拍图片 交配伊甸园 超碰在线视频自拍偷拍国产 小热巴91大神 rctd 045 类似于A片 超美大奶大学生美女直播被男友操 男友问 你的衣服怎么脱掉的 亚洲女与黑人群交视频一 在线黄涩 木内美保步兵番号 鸡巴插入欧美美女的b舒服 激情在线国产自拍日韩欧美 国语福利小视频在线观看 作爱小视颍 潮喷合集丝袜无码mp4 做爱的无码高清视频 牛牛精品 伊aⅤ在线观看 savk12 哥哥搞在线播放 在线电一本道影 一级谍片 250pp亚洲情艺中心,88 欧美一本道九色在线一 wwwseavbacom色av吧 cos美女在线 欧美17,18ⅹⅹⅹ视频 自拍嫩逼 小电影在线观看网站 筱田优 贼 水电工 5358x视频 日本69式视频有码 b雪福利导航 韩国女主播19tvclub在线 操逼清晰视频 丝袜美女国产视频网址导航 水菜丽颜射房间 台湾妹中文娱乐网 风吟岛视频 口交 伦理 日本熟妇色五十路免费视频 A级片互舔 川村真矢Av在线观看 亚洲日韩av 色和尚国产自拍 sea8 mp4 aV天堂2018手机在线 免费版国产偷拍a在线播放 狠狠 婷婷 丁香 小视频福利在线观看平台 思妍白衣小仙女被邻居强上 萝莉自拍有水 4484新视觉 永久发布页 977成人影视在线观看 小清新影院在线观 小鸟酱后丝后入百度云 旋风魅影四级 香蕉影院小黄片免费看 性爱直播磁力链接 小骚逼第一色影院 性交流的视频 小雪小视频bd 小视频TV禁看视频 迷奸AV在线看 nba直播 任你在干线 汤姆影院在线视频国产 624u在线播放 成人 一级a做爰片就在线看狐狸视频 小香蕉AV视频 www182、com 腿模简小育 学生做爱视频 秘密搜查官 快播 成人福利网午夜 一级黄色夫妻录像片 直接看的gav久久播放器 国产自拍400首页 sm老爹影院 谁知道隔壁老王网址在线 综合网 123西瓜影音 米奇丁香 人人澡人人漠大学生 色久悠 夜色视频你今天寂寞了吗? 菲菲影视城美国 被抄的影院 变态另类 欧美 成人 国产偷拍自拍在线小说 不用下载安装就能看的吃男人鸡巴视频 插屄视频 大贯杏里播放 wwwhhh50 233若菜奈央 伦理片天海翼秘密搜查官 大香蕉在线万色屋视频 那种漫画小说你懂的 祥仔电影合集一区 那里可以看澳门皇冠酒店a片 色自啪 亚洲aV电影天堂 谷露影院ar toupaizaixian sexbj。com 毕业生 zaixian mianfei 朝桐光视频 成人短视频在线直接观看 陈美霖 沈阳音乐学院 导航女 www26yjjcom 1大尺度视频 开平虐女视频 菅野雪松协和影视在线视频 华人play在线视频bbb 鸡吧操屄视频 多啪啪免费视频 悠草影院 金兰策划网 (969) 橘佑金短视频 国内一极刺激自拍片 日本制服番号大全magnet 成人动漫母系 电脑怎么清理内存 黄色福利1000 dy88午夜 偷拍中学生洗澡磁力链接 花椒相机福利美女视频 站长推荐磁力下载 mp4 三洞轮流插视频 玉兔miki热舞视频 夜生活小视频 爆乳人妖小视频 国内网红主播自拍福利迅雷下载 不用app的裸裸体美女操逼视频 变态SM影片在线观看 草溜影院元气吧 - 百度 - 百度 波推全套视频 国产双飞集合ftp 日本在线AV网 笔国毛片 神马影院女主播是我的邻居 影音资源 激情乱伦电影 799pao 亚洲第一色第一影院 av视频大香蕉 老梁故事汇希斯莱杰 水中人体磁力链接 下载 大香蕉黄片免费看 济南谭崔 避开屏蔽的岛a片 草破福利 要看大鸡巴操小骚逼的人的视频 黑丝少妇影音先锋 欧美巨乳熟女磁力链接 美国黄网站色大全 伦蕉在线久播 极品女厕沟 激情五月bd韩国电影 混血美女自摸和男友激情啪啪自拍诱人呻吟福利视频 人人摸人人妻做人人看 44kknn 娸娸原网 伊人欧美 恋夜影院视频列表安卓青青 57k影院 如果电话亭 avi 插爆骚女精品自拍 青青草在线免费视频1769TV 令人惹火的邻家美眉 影音先锋 真人妹子被捅动态图 男人女人做完爱视频15 表姐合租两人共处一室晚上她竟爬上了我的床 性爱教学视频 北条麻妃bd在线播放版 国产老师和师生 magnet wwwcctv1024 女神自慰 ftp 女同性恋做激情视频 欧美大胆露阴视频 欧美无码影视 好女色在线观看 后入肥臀18p 百度影视屏福利 厕所超碰视频 强奸mp magnet 欧美妹aⅴ免费线上看 2016年妞干网视频 5手机在线福利 超在线最视频 800av:cOm magnet 欧美性爱免播放器在线播放 91大款肥汤的性感美乳90后邻家美眉趴着窗台后入啪啪 秋霞日本毛片网站 cheng ren 在线视频 上原亚衣肛门无码解禁影音先锋 美脚家庭教师在线播放 尤酷伦理片 熟女性生活视频在线观看 欧美av在线播放喷潮 194avav 凤凰AV成人 - 百度 kbb9999 AV片AV在线AV无码 爱爱视频高清免费观看 黄色男女操b视频 观看 18AV清纯视频在线播放平台 成人性爱视频久久操 女性真人生殖系统双性人视频 下身插入b射精视频 明星潜规测视频 mp4 免賛a片直播绪 国内 自己 偷拍 在线 国内真实偷拍 手机在线 国产主播户外勾在线 三桥杏奈高清无码迅雷下载 2五福电影院凸凹频频 男主拿鱼打女主,高宝宝 色哥午夜影院 川村まや痴汉 草溜影院费全过程免费 淫小弟影院在线视频 laohantuiche 啪啪啪喷潮XXOO视频 青娱乐成人国产 蓝沢润 一本道 亚洲青涩中文欧美 神马影院线理论 米娅卡莉法的av 在线福利65535 欧美粉色在线 欧美性受群交视频1在线播放 极品喷奶熟妇在线播放 变态另类无码福利影院92 天津小姐被偷拍 磁力下载 台湾三级电髟全部 丝袜美腿偷拍自拍 偷拍女生性行为图 妻子的乱伦 白虎少妇 肏婶骚屄 外国大妈会阴照片 美少女操屄图片 妹妹自慰11p 操老熟女的b 361美女人体 360电影院樱桃 爱色妹妹亚洲色图 性交卖淫姿势高清图片一级 欧美一黑对二白 大色网无毛一线天 射小妹网站 寂寞穴 西西人体模特苍井空 操的大白逼吧 骚穴让我操 拉好友干女朋友3p