Next Article in Journal
Dengue Seroprevalence and Seroconversion in Urban and Rural Populations in Northeastern Thailand and Southern Laos
Next Article in Special Issue
Screening of Emerging Pollutants (EPs) in Estuarine Water and Phytoremediation Capacity of Tripolium pannonicum under Controlled Conditions
Previous Article in Journal
Prevalence of Sarcopenia in Community-Dwelling Older Adults in Valencia, Spain
Previous Article in Special Issue
Mitigation of Nickel Toxicity and Growth Promotion in Sesame through the Application of a Bacterial Endophyte and Zeolite in Nickel Contaminated Soil
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Long-Term Heavy Metal Retention by Mangroves and Effect on Its Growth: A Field Inventory and Scenario Simulation

1
Institute for Environment and Resources, Vietnam National University Ho Chi Minh City, 142 To Hien Thanh, District 10, Ho Chi Minh City 72506, Vietnam
2
Institute of Geoecology, Technical University of Braunschweig, 38106 Braunschweig, Germany
3
National Institute of Applied Mechanics and Informatics, Vietnam Academy of Science and Technology, 291 Dien Bien Phu, District 3, Ho Chi Minh City 722007, Vietnam
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2020, 17(23), 9131; https://doi.org/10.3390/ijerph17239131
Submission received: 20 October 2020 / Revised: 30 November 2020 / Accepted: 5 December 2020 / Published: 7 December 2020

Abstract

:
The ability of mangroves in taking up and storing heavy metal (HM) helps in reducing HM pollution. However, HMs likewise adversely affect the growth of mangroves. We assess the effects of the long-term soil HMs enrichment on the growth of Rhizophora apiculata forest in the Can Gio Mangrove Forest (Southern Vietnam) in different environmental conditions of soil salinity, ground elevation, and tree density based on a novel set of measured data. These data were analyzed and were used to calibrate and validate for a tree growth model with influencing factors salinity, elevation, tree density, and heavy metals content. Three scenario simulations were performed to predict the mangrove dynamics under different levels of heavy metal pollution in combined environmental conditions of salinity and elevation. Simulation results show the decline of total forest biomass from 1,750,000 tons (baseline scenario with no HM pollution) down to 850,000 tons and 350,000 tons for the current HM pollution and double HM pollution scenarios, respectively. Both data analysis and simulations have shown that although mangroves can assist in reducing HM pollution, the quality and health of this ecosystem will be severely affected if the environment is excessively polluted. In addition, a data-and-model driven management tool is devised for the sustainable management of the mangrove environmental resources.

Graphical Abstract

1. Introduction

Mangrove ecosystems are adversely affected by global and local stresses. Environmental stress, such as pollution, and especially heavy metal (HM) pollution, has been extensively addressed worldwide because of its adverse effects on ecosystems and human health. Excess of sedimentary HM concentrations constitute a potential risk of ecosystem health. It is a challenge to devise techniques for the reduction and cleanup of pollutants in the highly dynamic environment of mangrove forests.
Naturally, mangrove trees have the ability for storing metals, transferring these elements from the sediment and concentrating them in their tissues [1]. They act as a biological uptake medium and are a convenient means for phytoremediation [2]. Mangroves remove metals through absorption, cation exchange, filtration, and chemical changes through root accumulation processes [3]. There are many cases in which wetland plants (mangroves) are utilized for removal of pollutants, including metals. This functionality of mangroves has been studied and proven to be a good practical method to treat pollutants [4,5,6,7,8,9,10].
However, excessive heavy metal pollution causes a series of physiological and biochemical alterations affecting the growth, metabolism, and cell structure of plants [11] with the implication of potential ecological system risks [12,13]. Plants may alter the speciation of metals and may also suffer toxic effects as a result of accumulating them. Plant’s behaviors when being stressed by heavy metals were examined by greenhouse experimental studies, which applied for young mangroves, including oxidative stress (malondialdehyde (MDA) [14]), osmotic stress [15], inhibition of photosynthesis, inhibition of growth, and increased mortality [16].
To assess plants stress by heavy metals in the field, extensive empirical observations are required. However, in field observations, spatial heterogeneity concerning ecotope composition, initial contaminant levels and other factors make a consistent data analysis difficult. This holds especially true for the situation in large ecosystems with spatially varying environmental and pollutant conditions such as in mangrove areas. The structure of mangrove forests, their health and their ecosystem functions are strongly influenced by the factors resources (nutrients), regulators (soil salinity, pollutants) and hydro-regimes (duration, frequency and depth of inundation) and their mutual interactions [17,18,19,20,21].
This problem can be overcome by employing mathematical modelling tools which incorporate different factors into an integrated system capable of predicting ecosystem behavior under multiple stressors. While field inventories are resource-demanding, advances in mathematical modeling are now providing globally coterminous datasets at high spatial resolutions. There are quite a few different modelling approaches. Individual-based models (IBMs) for mangroves were developed for the assessment of the effects of natural threats and applied to the management of mangrove forest stability for coastal protection using a combination of windthrow models [22,23,24]. Quantitative statistical and a qualitative categorical ranking approaches were also used to identify risk factors for mangrove lost in Kenya (Rideout et al. [25]). Box models were applied for establishing mass balances for heavy metals. In the paper of Mandal et al. [26] the up-taken as budget for the Sundarbans mangroves was calculated. A phytoremediation model for the mangrove species Rhizophora apiculata was developed [27] coupling growth dynamics of plant organs with Cr uptake taking into account adverse effects of Cr in the plant. Also, a phytoremediation model at catchment scale combining hydrodynamics, Cr transport, and fate of chromium in the soil-plant domain [28] was realized by integration of an environmental fate model into the hydrodynamic model Delft3D and applied to the Thi Vai catchment in Vietnam.
The emerging problems in research on the mangrove retention of HM in relation to its growth concern the tradeoff between the phytoremediation function and the adverse effects of HM. In this context, we focus on R. apiculata forest in the Can Gio area, southern Vietnam. Previous studies in this area mainly focused on the potential of biological effects and ecological risk of heavy metal accumulation in the sediments. Due to the important role of ecosystem services provided by the Can Gio Mangrove [29,30], it was said that more attention should be paid on long-term trends of HMs and the complex interplay between hydrodynamics, anthropogenic stressors and HM concentrations, as well as ecosystem health [29,31,32]. Thus, our study aims to:
(i)
Review and re-evaluate HMs contamination degrees in mangrove soil;
(ii)
assess the effects of the long-term HM pollution on the growth of R. apiculata in the interactions with the effects of the natural factors soil salinity, ground elevation and tree density; and
(iii)
perform scenario simulations of mangrove dynamics under different levels of HM pollution and other natural stressors.
Establishment of reasonable model parameters is a major issue because of the lack of appropriately designed long term field experiments. In our study area we are in a favorable situation since the mangroves were planted and their age were recorded and the data on elevation and salinity are available. So this situation constitutes, though not planned in the first place, an “experimental design” suitable for the purpose of parameter estimation. To this end, the set of parameters for growth response functions of a previously developed mangrove dynamic model CGMM [33] were estimated based on the measured data in 208 sites in the Can Gio mangrove forest.

2. Materials and Methods

2.1. The Can Gio Mangrove Forest

The Can Gio Mangrove Forest (in the South of Vietnam) is located ca. 40 km off the biggest industrial city in Vietnam—Ho Chi Minh City (Figure 1). This forest is one of the most extensive rehabilitated mangrove forests in the world. The planted mangrove species Rhizophora apiculata Blume covers most of the forest. Rhizophora apiculata in this forest has been growing over 40 years since 1978. The structure of this mangrove forest is remarkably diverse and provides an excellent case study for research. Regeneration is assisted by planting of seedlings/saplings, but the trees also commence natural generative reproduction and contribute to the propagule input. The silvicultural practice used to consist of thinning events which removed competing dominant trees. Since this forest was nominated as “Biosphere Reserve” by UNESCO in 2000, the forest was strictly protected, and the activity of forest thinning was inhibited. The mangrove ecosystem in this area was said providing not only various goods and services to local people such as timber, seedlings, medicines, but also foods because variety of marine organisms live within and around this mangrove, e.g., shrimp, fish, snail, or crab [30].
This mangrove area is subject to asymmetric semidiurnal tidal regimes and typical tropical monsoon climate with two distinct seasons, dry season (extends from November to next April) and rainy season (from May to October). Mean water salinity varies from 12 to 30 part per thousand (ppt). Salinity changes by season, highest salinity occurs in dry season and reaches the value of more than 30 ppt in March and April. Flooding frequency and elevation of the study area affect primary production and species distribution in Can Gio [34,35]. van Loon et al. [36] classified Can Gio mangroves into five classes based on tidal regimes: High tidal areas (no mangroves observed), medium high tidal areas (Avicennia spp., Sonneratia), normal high tidal areas (Rhizophora spp., Ceriops, Brugueira), spring high tidal areas (Lumnitzera, Brugueira, Acrostichum aureum), and equinoctial tidal areas (Ceriops spp., Phoenix paludosa).
The area has been found being polluted by HMs. Our previous publication [37] reported high concentrations of Cr, Cu, and Ni in the soils and tree organs among the 18 elements extracted from samples in the 11 locations as shown in Figure 1 (the triangle blue points). Costa-Böddeker et al. [31,32,38] also gave evidences of the long-term pollution of these HMs in this area. Mean sedimentation rates calculated from the sediment cores are remarkably high (9.2 cm year−1 upstream and 4.7 cm year−1 downstream) due to the rapid urbanization and industrialization in the recent two decades as was reported by Costa-Böddeker et al. [32,38]. However, natural processes such as tidal currents and heavy rainfall also play a key role on sediment dynamics in this area [39]. These high sediment fluxes were also associated with increasing HMs contaminants in the area.

2.2. Sampling Positions, Samples Preparation, and Analysis

A comprehensive field investigation was conducted in 2018 for 208 sampling plots of 10 × 10 square meters (the points in Figure 1) located in the area. The key growth driving factors inundation frequency (as expressed via ground elevation) and soil water salinity were determined for each plot. In addition, measurements of the mangrove tree heights (H), diameters at breast height (dbh), and tree density were conducted. The age of the forest stands is well known because it is a planted monospecific forest mostly covered by R. apiculata. The monospecific structure and the even age favor the insight into species–environment interactions. Thus, this monospecific planted forest plays an important role as a long-term experiment allowing the study of influences of environmental drivers on tree development.
The salinity (in part per thousand, ppt) of pore water in the soils was determined using an Atago refractometer. After extracting a drop of interstitial water from sediments, pore-water from the sediments was collected using a syringe. In the 11 out of the 208 sampling plots (marked by triangle blue points in Figure 1) HMs concentrations in soils, roots and leaves were determined. These plots are located at the same positions as the 2013-sampling plots from our previous work (Nguyen et al. [37]).
Samples of soils (0–30 cm depth) were collected and covered carefully with plastic scrapers. Soil samples were homogenized in an agate mortar. Every 5 g of the oven dried homogenized soil samples were directly digested with aqua regia. The concentrations of HMs (Cr, Cu, Fe, Ni) in the aqua regia extracts were determined with an ICP OES (Optima 2100 DV, PerkinElmer, USA). The aqua regia extraction for Fe is according to the ISO 12914. Up to 0.5 g of sample was placed in a Teflon vessel with 12 mL of aqua regia (37% HCl:70% HNO3 (3:1) mixture). The vessels were heated in a microwave apparatus (MARS 6-CEM, USA) up to 180 °C within 6 min and remained at 180 °C for 10 min. For the other HMs, the extraction was based on the ISO 11466.3. Up to 5.0 g of sample was placed in a 250 mL Pyrex beaker with 40 mL of aqua regia (37% HCl:70% HNO3 (3:1) mixture). The beaker was heated at 130 °C for about 30 min and evaporated almost to dryness. The residue was re-suspended using HNO3 1.6 M. After centrifuging the supernatant was diluted to 25 mL [40,41].
Tree roots and leaves were collected on one adult tree in each of the 11 sampling plots. The samples were analyzed at the Hochiminh city Institute of Resource Geography, Vietnam Academy of Science and Technology, Vietnam. All plant samples were rinsed with pure water to eliminate residues of sea water and soil attached. Representative parts of the plant samples and aliquots of the soil samples were air dried at room temperature. The cut root samples and leaves were milled under addition of liquid nitrogen. Aliquots of 1 g were dried in an oven at 105 °C according to DIN EN 212 (2003) until constant mass was reached. The dry mass was used as base for the calculations of the concentrations of elements. Aliquots of 3 g to 10 g of the air dried and milled plant samples were reduced to ashes according to DIN 19684-3 (2000) at 550 °C and the residues were digested with aqua regia according to DIN EN 13346 (2000).

2.3. Tree Growth Model

The structure and conceptual basis of the growth model are based on the classical Botkin model [42]. The model equations are given by the differential equation for the tree diameter (dbh) and the relationship between diameter and tree height (H).

2.3.1. Tree Growth Equation

The growth equation is expressed in the Equations (1)–(4) [42].
d d b h d t = G o p t × d b h · 1 d b h · H D m a x · H m a x 2 b 1 + 3 b 2 · d b h 4 b 3 · d b h 2 × M U L ,
H = b 1 + b 2 · d b h b 3 · d b h 2
b 2 = 2 H m a x b 1 D m a x
b 3 = H m a x b 1 D m a x 2
B i o m = a 1 · d b h c 1
where: dbh is the diameter at breast height in cm; H is the tree height in cm; H m a x and D m a x are the species specific maximum values of tree height and tree diameter; b 1 , b 2 , and b 3 are the parameters characterizing the relationship between H and dbh; G o p t is the species specific growth rate under optimal environmental condition; a 1 and c 1 are scaling parameters for calculating tree biomass B i o m from dbh value; MUL is the growth multiplier depending on environmental conditions and tree competition. It takes values between 0 and 1. In this study, MUL comprises of the influencing factors salinity, ground elevation, tree density and heavy metal pollutants.
Here we omit the inner physiological differentiation among trees of the same species and assume that mangrove structure and productivity are governed by environmental factors and by their own interactions (competition). Thus, at a given age, trees of the same species will attain the same size under optimal conditions for tree growth.
The form of multiplier function for salinity and density factors is
g i x = 1 a 0 i 1 + e d i t r i x i + a 0 i .
This multiplier function is determined by the threshold value t r i , the slope d i and the asymptotic value a 0 i ensuring that the growth rate remains at a small level even at high values of salinity or tree density. The index i refers to the salinity and tree density, respectively. This approach is motivated by the salinity dependence of the growth of mangroves as published by Chen and Twilley [22] or by Berger and Hildenbrandt [24].
Elevation can be considered as the factor that reflects the frequency and duration of tidal flooding, thus, instead of considering two factors that influence tree development (tidal frequency and duration) we integrated both into only one factor (elevation). The multiplier function for elevation is (from Anh [33]).
f e l = a m a x a 1 e 1 e e l e l m i n e l 1 α + a 1 e e e l e l m i n e l 2 β + a 2 e 1 e e l e l m i n e l 2 β .
It is conceived such that it takes on the minimal value a 1 e at low elevations and the minimal value a 2 e at high elevations, exhibiting an optimal value or a plateau of optimal values in between. It is set equal to zero below the minimum elevation e l m i n .
The form of multiplier function for pollutant is
P x = 1 e x t h 1 α 1 · e x t h 2 α 2 .
This function for pollutants takes into account hormetic growth stimulation. Hormetic growth stimulation stands for stimulatory effects at low doses and was reported for Cd, Cr, Al, Se, and Pb [43,44]. The form of this response curve is determined by the threshold value t h i ,   i = 1 , 2 and the form factors α 1 ,   and   α 2 . Below threshold value t h 1 of HM concentration tree growth rate increases, above threshold value t h 2 it decreases.

2.3.2. Calculation of Tree Growth Rate G

Under influences of environmental factors, the actual value of the growth rate G of a tree is always smaller than the optimum G o p t value and is reduced by multipliers MUL reflecting the influence of environmental conditions: G = G o p t · M U L .
In order to analyze relationships between influencing factors and tree development, values of tree growth rate G must be known. Analysis of tree development based on G values is more precise than based on size of trees (diameter and height) since G averages the growth of a tree taking into account also the factor tree age and influencing factors ( M U L ).
Integration of Equation (1) yields an expression of G as a function of the diameter at the time of planting, d b h 0 , and at the time of measuring (age at the time of measurement), d b h a g e [33]
G = 2 · H m a x b 1 + a 3 e x 2 e 2 x a + b 1 e x b 1 + a 2 e x e 2 x d x C t ,
where a = H m a x b 1 , x 0 = l n d b h 0 / D m a x , x 1 = l n d b h a g e / D m a x , C is an integration constant.

2.4. Data Analysis and Model Parameters Estimation

2.4.1. Data Analysis

HM levels were evaluated by reference to the Enrichment Factor (EF) [45,46]. This factor normalizes the element concentration with respect to a reference element which is not expected to be enriched from anthropogenic sources. Therefore, our assessment was based on the EF index. We used iron (Fe) as reference element. The EF is calculated as follow (Equation (10)):
E F = M e F e S a m p l e M e F e B a c k g r o u n d ,
where M e F e S a m p l e is the ratio of the concentrations of the element to Fe in the sample and M e F e B a c k g r o u n d is the ratio of the concentrations of the reference material (background values). EF values lower than 1.5 suggest that the element is derived mainly from natural sources, EF values higher than 1.5 suggest anthropogenic sources [47] whereas EF > 5 shows a significant HM enrichment [46].
Continental shale (CS) averages [48] were used as background values for the calculations of this index. The use of CS instead of regional levels as background values could lead to a conflicting pollution classification which may hamper the assessment of environmental impacts in this region [38,49]. Unfortunately, it was not possible to get data before the period of anthropogenic pressures which considerably modified this area.

2.4.2. Model Parameter Estimation

The growth model as explained in the Section 2.3 needs to be parameterized for the mangrove species R. apiculata. This parameterization is based on the 208 measured data of tree stem heights and stem diameters at breast height with age ranging from 19 to 40 years.
Among the analyzed 208 data plots, 65 plots were chosen for model parameter estimation. The locations were selected according to ranging salinity, elevation and tree density. Other 11 plots which were sampled for the determination of HM contents were used to estimate parameters for the multiplier equation of pollutant (Equation (8)). The remaining 132 plots were used for model validation. Biomass data were taken from the work of Nam et al. [50].
Parameters were estimated by nonlinear regression techniques. This regression problem was solved by embedding an ODE (Ordinary differential equation) solver in the optimization toolbox of Matlab R2018b. For the estimation of parameters of biomass function and multiplier functions as shown in Equations (5)–(8), the regression problem was solved by standard nonlinear regression techniques (using Matlab R2018b).

2.5. Scenarios Setting

Scenarios modeling is a useful tool to understand the responses of a system to various factors and their combinations. In this study, we explored the trunk biomass of R. apiculata forest in the Can Gio area, which was given by Equation (5). Tree trunks represent the largest contribution to the total dry weight (DW) above-ground biomass. Three scenarios were set up:
(i)
The baseline scenario (SC0) which assumed the unpolluted environmental condition. This scenario projects how the system may have developed within a four-decade period (1978–2020) under ideal conditions of no HM pollution.
(ii)
Scenario based on actual chromium pollution (SC0P1) in the area in accordance with our observational data.
(iii)
Worst case scenario assuming a twofold chromium load (SC0P2).
Tree densities in the forest stands were randomly generated by a uniform distribution in the range of 10–115 trees per plot (10 × 10 sq. meters) in accordance with our observational data.
The optimal number of scenarios is generally considered to be three or four [51]. The appropriate scenarios design should provide consistent outcomes that enable both the assessment of current conditions and assist future planning. We determined that comparing a baseline (no pollution) scenario with two clearly distinct scenarios of different HM levels would provide clear, useful, and visible alternatives for management planning.
In the simulation, the tree growth model was embedded into a life cycle model CGMM including the phases “Establishment–Growth–Competition–Reproduction–Mortality” (Anh [33]. Detailed explanation and equations for the life cycle are shown in Anh [33] Appendix C). Input data for the simulations include a topographic/bathymetric map of the study area, a map of salinity distribution and a map of root Cr content distribution. Note that we choose root concentrations of Cr as prediction variable rather than soil concentrations. The reason is that soil data reflect only the current environmental pollution state whereas root concentrations are the result of the HMs accumulation into the trees. Besides, a map of R. apiculata forest was also prepared. This map is available because this mangrove species was planted, and its data was well recorded. All these maps are shown in Figure 2 with a spatial resolution of 100 m × 100 m.
Chromium was chosen as the HM pollutant for our simulation because Cr is the main pollutant in our study area (in soil and in plant tissues). We also conducted experimental studies on the growth of young R. apiculata in dependence on Cr concentrations in combination with other factors [27,52]. Thus, data of Cr effect on plant growth are available to support the evaluation of simulation results. In our study, Cr oxidation forms were not specified; however, in general, both speciation forms Cr(III) and Cr(VI) occur in this area [31].
Other assumptions were: (1) All locations/plots used for simulations were fully forested with R. apiculata; (2) current mangroves at all sites were rehabilitated at the beginning and were naturally regenerated from seedlings of adult trees; (3) there were no disturbances following stand establishment. The simulation covers the time period from the year 1978 when R. apiculata forest was planted until the year 2020. The runs were initialized by a random uniform distribution of seedlings in each stand.

3. Results and Discussion

3.1. HMs Pollution in the Soil and Potential of HMs Retention by R. apiculata

3.1.1. Soil HMs Pollution in the Area

Table 1 shows the HM data in soil, root, and leaf at the 11 sampling positions. The mean values of these HMs contents in the soil were 77.47 mg/kg for Cu, 9.82 mg/kg for Cr, and 5.87 mg/kg for Ni.
For conveniently describing the emerging polluted situation in the area, we divide the study area into two regions, the Thi Vai area and the Can Gio area (see Figure 1). Thi Vai catchment (Figure 1) is affected by the industrial activities on the East bank and the Can Gio area is directly influenced by activities from the center of Ho Chi Minh city (see Figure 1). The distribution of HMs contents among the different locations in the Can Gio (positions NP1, DK4, DK6, DK7, DK8, NP2) is more homogeneous than in the Thi Vai catchment (positions DN1, DN2, M5, M3, M1 in Figure 1). Standard deviation values of HM contents in the soil in Can Gio locations are smaller than those in the Thi Vai ones, they range from 0.28 to 6.29 for Can Gio and 2.48 to 29.79 for Thi Vai catchment (cf. Table 1). The HMs contents in the soils in the Thi Vai catchment are higher than those in the Can Gio, mean values for Cu, Cr and Ni respectively in the Thi Vai are 83.76, 16.35, and 9.30 mg/kg and in the Can Gio are 72.23, 4.38, and 3.01 mg/kg. Cu content is the highest among the three HMs both in Can Gio and in Thi Vai (the highest content is 134.03 mg/kg at the location DN2). Compared to our measurement campaign in 2013 [37], the pollution situation has generally improved. Chromium and Ni contents are lower than those in 2013 and their EF values are smaller than 1.5 (Table 2) falling in the EF class “not enriched”. But Cu (in 2018) in average is 3.2 times higher than this in 2013 (Cu content ranges from 60.4–134.03 mg/kg in 2018 and from 19.7–32.7 mg/kg in 2013 respectively, cf. Table 3). Its EF values are larger than 1.5 at all locations except DN1 (cf. Table 2) showing that Cu is moderately to highly enriched in the whole area. However, although concentrations of Cr and Ni below the limits for “enriched” (EF value < 1.5), they may not be related to low contamination levels as dilution effects may mask high HM rates of deposition (as was discussed by Costa-Böddeker et al. [38]).
The HM pollution in this region has happened for about two decades since ~2000. The Thi Vai river and its surroundings have been severely polluted because of the direct discharge of wastewater from the industrial companies. Costa-Boddeker et al. [32] analyzed the sediment cores in this region and reported that the violations of HM pollution especially for Cu, Cr, and Ni were observed particularly after the establishment of the industries in the area in the late-1990. Concentrations of these HMs from 2003 to 2012 ranged from 107.70–208.80 mg/kg for Cr, 26.74–82.32 mg/kg for Cu, and 56.25–82.99 for Ni respectively in the upstream section. These concentrations were lower in the downstream section with values range from 51.6–82.5 mg/kg for Cr, 11.5–38.3 mg/kg for Cu, and 24.7–46.5 mg/kg for Ni respectively (Costa-Boddeker et al. [32], see Table 3). These HMs contents are in consist with the increment number of industrial companies in this catchment. The number of industrial companies in the Thi Vai catchment increased from 10 companies in 2002 to 135 companies in 2017 (data from the Ho Chi Minh City Department of Natural Resources and Environment) with main activities in tannery, mechanical engineering, electrical engineering, packaging, textile and dyes industry, oil activities and cement production. Until 2006, over 60% of the industrial zones had not implemented an adequate wastewater treatment system and the untreated wastewater was discharged directly into the river [31].
For the Can Gio area, Costa-Böddeker et al. [38] reported the increasing trend of HMs pollution was from ~2002. Chromium, Cu and Ni were also the highest concentrations in sediments ranging from 27.1 μg/g to 71.5 μg/g for Cr, 7.1 μg/g to 27.0 μg/g for Cu, and 11.7 μg/g to 56. 3 μg/g for Ni (Table 3). Beside those HMs, the area has also been polluted by Pb and Co with the order Cr > Ni > Cu > Pb > Co [38]. Following Costa-Böddeker et al. [38], enhanced HM concentrations in the sediment of the Can Gio area may have originated from both natural and anthropogenic sources, with the latter being the main contributor in recent years. The anthropogenic sources are mostly associated with the industrialization process in Southern Vietnam, the impact of shrimp farms as well as oil spills and untreated wastewater. Factors such as shoreline erosion caused by tides and deforestation may also have contributed to these HMs distribution [38].
HM pollution occurring in the soils of mangrove areas worldwide are mostly attributed to industrial activities. Compared to other mangrove regions, the contents of HMs associated mostly with anthropogenic sources in our region are lower than in some other industrialized mangrove regions, see for instance the Maipo mangrove area in Hong Kong or Surabaya in Indonesia and Yanbu in Saudi Arabia (Table 3). However, Cu contents are higher than in other contaminated areas (with values from 1.2 to 4-fold higher, Table 3), particularly within the sites upstream in Can Gio (NP1, DK4) and downstream in Thi Vai catchment (DN2, M3, M1). The intensified industrialization over the past two decades has brought about the moderate to high enrichment of the HMs Cu, Cr and Ni. It was said that the complex interplay of the dynamics of sediments such as increasing erosion, shoreline retreat and transportation of HMs may be hiding the real state of HM pollution and ecological disturbances in this region [32].

3.1.2. Potential of the HMs Retention by the Mangrove Species R. apiculata

In the tissues of R. apiculata, the concentrations of Cr, Cu, and Ni in the root and the leaf are shown in Table 1. In the root, the mean values of Cr, Cu, and Ni are 2.63, 2.52, and 2.01 mg/kg respectively, ranging from 0.26 to 11.69 mg/kg for Cr, 0.21 to 8.17 mg/kg for Cu and 0.11 to 8.54 mg/kg for Ni. In the leaf, the average contents of Cr, Cu, and Ni are 0.83, 2.1 and 0.39 mg/kg respectively, ranging from 0.23 to 2.05 mg/kg for Cr, 0.18 to 7.83 mg/kg for Cu and 0.005 to 1.81 mg/kg for Ni. Among the positions, Cu contents in tree tissues (root and leaf) at locations along the Thi Vai river (locations DN1, DN2, M5, M3, M1) are higher than those in the Can Gio (NP1, DK4, DK6, DK7, DK8, NP2). Even though the concentration of Cu in the soil is high compared to this measured in 2013 (Table 3), this content is low in plant tissues both in leaf and root compared to the tissue values measured in 2013. This could be explained by the incident time of HM. Copper might have just recently entered the environment. Similar to what was found in the soil, the contents of HMs are nearly homogeneous for all locations in the Can Gio but not for the Thi Vai region (Table 1, standard deviation ranges from 0.15 to 0.41 for HMs in the root and 0.19 to 0.56 for HMs in the leaf in the Can Gio area, respectively; while this value ranges from 2.01 to 4.03 for HMs in the root and from 0.66 to 2.16 for HMs in the leaf respectively in the Thi Vai catchment).
Average Cr concentration in the root of R. apiculata is the highest concentration among the three HMs in the root, but this value is lower than in the roots of other mangrove species of other regions (see Table 3) except in some mangrove roots in Maipo, Hong Kong (~2 times higher). While in the leaf, Cu concentration is the highest among the three HMs. However, these three HM concentrations in our leaves are still lower than in the leaves in other mangrove regions as shown in Table 3.
The HMs accumulated in root and leaf tissues clearly demonstrate the potential of R. apiculata in the retention of HM. Because these HM contaminants in the soil occurred since ~2003, the trees accumulated high concentrations of these soil enriched HMs in their tissues. It was found in our previous work Nguyen et al. [37] that the plant tissues at 11 sampling positions (the same positions as in this work and is shown in Figure 1) accumulated high amounts of Cr, Ni, and Cu especially in the tree cores [37]. Tree stems include the sapwood which engages in transport of water and minerals to the crown of the tree. The sapwood is periodically converted to heartwood (in the tree cores) which constitutes the layers of death sapwood cells. Therefore, by time, plants accumulate elements and immobilize them in their dead cells in the heartwood (in tree cores). The high contents of these soil contaminants in the tree cores corroborate the assumption that these elements have been accumulated since a long time.
The HMs retention by mangroves has been reported in many countries. Mangrove wetlands are often considered as sinks for contaminants. There are many cases in which wetland plants (mangroves) are utilized for removal of pollutants, including metals. This functionality of mangroves has been studied and proven to be a good practical method to treat pollutants. The study of Kamaruzzaman et al. [4] with R. apiculata in Malaysia corroborated the potential of metal accumulation especially in their root. MacFarlane et al. [5] conducted experiments and reported that Cu was accumulated in root tissues of Avicennia marina under a certain concentration limit (200 µg/g). Nazli and Hashim [6] found the accumulation of Cu and Pb in roots and leaves of Sonneratia caseolaris and concluded that the roots of this species have high capacity in taking up heavy metals and could be a potential phytoremediation species for heavy metal treatment in Malaysian mangrove ecosystems. Chowdhury et al. [7] concluded that the species Sonneratia apetala and Avicennia officinalis are potential hyperaccumulators which accumulated high amounts of Cu (25.89 mg/kg), Fe (1376.7 mg/kg), and Cr (2.85 mg/kg) in the pneumatophore. Fengzhong et al. [8] found that Avicennia marina has a high removal ability for lead (Pb, 83.8%) and Cd (74.2%) and that K. candel has ability to absorb copper (Cu, 70.5%) and nickel (Ni, 50.5%). Wen-jiao et al. [9] found that Rhizophora stylosa in Yinglou bay (China) absorbed high amounts of cadmium (Cd). Titah and Pratikno [10] discussed that Avicennia alba can be considered for use in phyto-monitoring and phytoremediation of Cr in coastal areas. In general, these authors stressed the potential role of mangroves in HM retention. Thus, there is no doubt that mangroves in general and R. apiculata in our study area possess a remarkable capacity to retain HMs and tolerate relatively high levels of metals.
Mangrove plants develop a series of enzymatic and nonenzymatic antioxidants to scavenge reactive oxygen species (ROS) and reduce the toxic effects of abiotic stresses, including HMs [56]. Field inventories and laboratory experiments showed that the mangrove trapping of HMs was a very efficient and fast phenomenon [56,57]. The tolerance of mangrove plants to HM stress is normally a mixture of metal avoidance and scavenging of reactive oxygen species (ROS) [56].
Mangrove soil with the anoxic nature, usually with negative redox potential and high amounts of sulfide, iron, and organic matter, decrease the metal solubility and bioavailability that makes sediment be a long-term sink [56,58,59]. Mangrove sediments generally have a high capacity for adsorbing and holding trace metals, however, excess loading of HMs may exceed the binding capacity of the sediment [56,60]. Salinity is an important variable for the distribution and absorption of HM in the soil [31,32]. It has a key role to control the bioavailability and toxicity of HM in sediments. Changes in salinity may cause the mangrove sediments to lose their metal-binding capacity, resulting in mobilization of metals [56]. The competition between major cations and HM for sorption sites is enhanced with increasing salinity [61]. The mangrove sediments then shift from a HM sink to a HM source, and such disruption is often associated with human activities [56,58].
Among the tree organs, root of mangrove tree has an important role in depurating water and sediment by retaining metals. Metals are generally highly bioavailable and accumulated in mangrove roots. The common studied HMs, such as Cu, Zn, Cd, Cr, and Hg, generally showed high bioconcentration factors in roots. Thus, under increasing HM pollution pressure from the rapid industrialization and urbanization, it is said that mangroves can serve as a natural wastewater treatment system, which could be considered as a phytoremediation process.

3.2. Growth of R. apiculata in Different Environmental Conditions

3.2.1. Tree Growth Rate (G) and Influencing Factors

The structure of mangrove forests, their health and their ecosystem functions are strongly influenced by the factors resources (nutrients), regulators (soil salinity, pollutants), and hydro-regimes (duration, frequency, and depth of inundation) and their mutual interactions [17,18,19,20,21]. The long-term influence of these driving factors on tree development is in a way summarized by the growth rate G. The data of G values can be extracted by inserting the growth parameters (208 data plots) into Equation (9).
The growth rate (G) extracted from 76 plots among the 208 data plots were used to estimate parameters for the multiplier Equations (6)–(8) which expressed the influences of above driving factors. Table 4 shows the estimated parameter values.
The extracted tree growth rate (G) from the data shows a negative correlation (r = −0.78) with the pore water salinity. Figure 3 shows growth rates of trees in the plots located on large salinity gradients and have similar tree density. Along the salinity gradient, growth rates (G) of trees decrease. This effect is enhanced at high tree density. This relationship agrees with other previous studies on the impacts of salinity on the growth of mangroves [62,63]. Reef et al. [64] supposed that nutrient availability of mangroves is controlled by the redox state of the soil, which may be influenced by salinity. Dai et al. [65] reported that salt stress influences mangrove productivity. A common response to salt stress is that the leaf surface expansion is reduced, thereby reducing photosynthesis, which in turn reduces carbohydrate production [66,67,68]. Experiments from our previous study [69] on young R. apiculata in this area showed that plants growing on salt-leached soil have the highest growth rates and on natural saline soil have the lowest growth rates. Salinity in the natural saline soil can be considered as a toxic agent influencing the uptake of nutrients and thus plant growth. This might be due to competitive inhibition of ammonium uptake by sodium. Thus, in high salinity, growth rate (G) of mangrove trees decreases.
Beside the physical environment, tree density also affects tree growth since the trees compete to each other for nutrient and light. Tree densities decline with forest age, notably because of the silvicultural practice of intermediary thinning for removing competing dominant trees [70]. However, tree density also declined in natural forests due to physiological tolerances and competitive interactions [71].
Another factor which directly influences mangrove growth is the inundation. Inundation stress and increasing sea level are influencing the growth of mangroves. These factors are reflected by the ground elevation. The higher water level provokes changes of elevation gradients. The main consequences are the decline of biomass, tree density, basal area, and C pool. Figure 4 shows the fitted curve of model growth rate G of R. apiculata and the G extracted from data. The large gradient of elevation influences the growth of mangroves such as those observed in our samples. This result shows that R. apiculata develops well within elevation range from −0.5 to 0.5 m. Gilman et al. [72] noted that mangrove species maintain their preferred hydroperiod, they may also expand laterally into areas of suitable elevation which are the environmental conditions for recruitment and establishment of mangroves in new areas. Schwarzer et al. [39] argued that any loss in elevation makes mangroves more susceptible to submergence from sea-level rise. Ellison and Farnsworth [73] reported that mangroves need to grow in a stable condition.

3.2.2. Tree Growth Rate (G) in Polluted Condition

The data of growth parameters (mean dbh, mean height H, age, mean growth rate G) and the environmental factors (soil salinity and EF values for Cr, Cu, and Ni in the soil) at the 11 measured forest stands are shown in Table 2. These growth parameters reflect the prevailing environmental conditions during the lifetime of a tree. The highly contaminated stands (Table 2, positions DK6, DK7, DK8, NP2, and M1 with EF (Cu) value > 2.5) have G values lower than the other stands. The lowest value of G appears at stand M1, where trees grow on interactive negative conditions of high salinity (17.1 ppt), high density (17 trees) and highly contaminated soil (EF (Cu) = 3.9, Table 2).
We also found that the HMs Cu, Cr, and Ni accumulated in plant roots negatively affect the plant growth at high concentration levels (Figure 5). These results suggest that mangrove growth is limited by the excessive presence of HMs which are taken up by the roots. Our data indicate that in average, Cr, Cu, and Ni concentrations accumulated in root over a concentration of 20 mg/kg are negatively correlated with growth of R. apiculata as expressed in Figure 5a–c. However, the lower G values at low concentrations may hint to a hormesis effect. Homesis effects of the heavy metals Pb, Cu and Cd were found in the experiments of Zhang et al. [74] with Sonneratia apetala. Unfortunately, the duration of their experiment was not long enough to observe fully plant behavior. The interactive effect of these HMs accumulated in tree root on tree growth rate G in relation to the soil salinity is also shown in Figure 5d.
For Cu, MacFarlane and Burchett [75] reported that Avicennia marina was found to be highly tolerant to Cu. It was accumulated in root tissue but at concentrations of 200 μg/g and higher, no further increases in root Cu levels occurred. Significant reductions in seedling height, leaf number and area were found with significant increases in Cu concentrations at 100 μg/g sediment Cu. At Cu sediment levels of 400 μg/g, a decrease in total biomass and root growth inhibition was observed. Chawla and Patel [76] observed changes in chlorophyll contents (total chlorophyll and ratio chl a/b were examined) by the application of different concentrations of Cu in seedlings of the mangrove species Avicennia marina (Forsk.)Vierh. The results indicated that the content of total chlorophyll increased at low levels of Cu but started declining at higher concentrations.
For Cr, our previous experiments in Richter et al. [27] and Nguyen et al. [69] showed that the phytoremediation performance of young R. apiculata in treating Cr was at the maximum Cr concentration at 600 mg/L in industrial wastewater. Oliveira [77] discussed that is Cr structurally similar to other essential elements and may therefore displace the nutrients from physiological binding sites.
For Ni, to our knowledge, there have been limiting studies on effects of excess Ni on mangrove growth. Nickel and Cu as micronutrients for plants play a role as brokers of redox transformations [78]. However, Ahmad and Ashraf [79] discussed that although Ni is essential to several plant metabolic phenomena, it is extremely toxic to plants when present at excessive levels in the soil or in nutrient solutions to which plants are exposed. The general signs associated with Ni toxicity in plants include reduced shoot and root growth, poor development of the branching system, and deformation of various plant parts.

3.3. Simulation Results

3.3.1. Model Validation

As already mentioned in Section 2.4.2, the data set was split into one data set for parameter estimation (76 plots) and one for validation (132 plots). In Figure 6 the dbhs of R. apiculata predicted from the model are plotted against the measured dbhs. There are no apparent systematic deviations between data and model predictions and the correlation coefficient is R2 = 0.87. The predictions of G values in dependence of the variables’ density and salinity together with the 95% confidence surface are shown in Figure 7. The coefficient of variation is R2 = 0.8902. One clearly recognizes the nonlinearity of the response surface both with respect to density and salinity (Figure 7).
It should be emphasized that the role of mono specific forest plantation here is especially important as a long-term experiment allowing the assessment of tree development of a single species in different environments. The validation results justify the application of the model for predicting the production of mangrove trees at each age in dependence on different environmental gradients and HM pollution loads. The model can thus be employed as a decision support tool for the management of mangrove forest.

3.3.2. Scenarios Simulation Results

Figure 8 shows snapshots of the spatial distribution of biomasses at four time points for the three scenarios described above. The histograms of biomasses at the final stage of simulations obtained for the three scenarios are given in Figure 9. The histograms are shifting to the left towards smaller biomass values with growing pollution levels and have lower variances.
Under the baseline scenario a rather inhomogeneous biomass distribution with large variation is obtained reflecting the unpolluted natural environmental conditions. Biomasses range from 5000 (~1800 stands) to 175,000 (~100 stands) kg/ha in 2020. Under the presence of HM, the range of biomasses becomes smaller and the distribution becomes more homogeneous (3000 stands and 6000 stands with trunk biomass 5000 kg/ha for SC0P1 and SC0P2, respectively). The occurrence of stands with high biomass are very limited (90 stands have the maximal biomass 125,000 kg/ha in SC0P1 and ~100 stands have the biomass >25,000 kg/ha for SC0P2, Figure 9).
The biomasses obtained from the simulations compare well with empirical data of other studies. In the work of Vinh et al. [80] the average measured trunk dry biomass of R. apiculata in Can Gio was reported 76,250 kg/ha (76.25 ± 5.13 tons/ha). It resulted from dry weighing 36 felled trees with ranging diameters from 7–36.2 cm at one specific site. In the study of Hoan [81] with 45 felled trees in different locations, the average trunk dry biomass of the R. apiculata was in the range of 47,600 to 139,800 kg/ha with diameters ranging from 9 to 21 cm. Maximum trunk biomasses in 2020 obtained from our simulations are 175,000 kg/ha in the baseline scenario (ideal environmental condition), 140,000 kg/ha and 100,000 kg/ha in the SCOP1 and SC0P2 scenarios respectively (cf. Figure 9). The slight difference in our predicted biomass compared to the data from the above authors are probably due to the large variation in tree size. In our study, biomass was obtained from the summation over all trees in one-hectare plots regardless of tree size (there exist young trees to adult trees), while in the referenced work above biomass was weighed based on the adult trees with starting diameter of 7 cm.
By our analysis of field data and simulation studies we have assessed the influence of HMs and other factors tree density, soil salinity and ground elevation on the trunk biomass of R. apiculata. Variation of trunk biomass of R. apiculata forest is much higher among the plots in SC0 with no HM pollution. This suggests that the interactions of soil condition and tree density with HM pollutants were responsible for the slower growth and for the smaller biomass in the two HM polluted scenarios (SC0P1 and SC0P2).
The simulations show the order of magnitudes involved in the reduction of biomass due to HM pollution. Total trunk biomasses in 2020 of trees growing under HM stress since 1978 amount to 850,000 tons and 350,000 tons for the SC0P1 and SC0P2 respectively, in contrast to 1,750,000 tons achieved under no HM stress conditions (for the SC0). The tree growth in the SC0 is influenced just by the natural factors of salinity, elevation, and tree density. These natural environmental factors have changed only gradually by the complex dynamics of the estuarine environment, so we can consider these factors quasi as stable environmental conditions (without considering the effects of climate change). Thus, to achieve mangrove stands with high biomass, we should control the pollution by reducing the amount of HMs release.

3.4. Toward a Balancing Management Approach on Using Mangroves to Clean Up Polluted Environment and to Protect the Mangrove Ecosystem

The Ministry of Agriculture and Rural Development (MARD) in Vietnam reported five main causes of the Vietnamese mangrove loss: (1) Non-sustainable use of mangrove ecosystems for aquaculture, (2) storms and natural disasters, (3) deforestation of mangrove forests for timber and natural resources, (4) multi-source pollutants from agriculture, industrial and urban areas, and (5) lack of sufficient regulatory mechanisms for the protection and sustainable development of mangrove ecosystems. Mangrove plants show high metal retention abilities. Because of this ability, mangroves have been proposed for utilizing for the removal of pollutants, including metals. The problem is, on the one hand, mangroves have an important role and function as a natural phytoremediation system; however, on the other hand mangrove ecosystems themselves are highly sensitive to all kinds of disturbances.
The high rate of heavy metal contamination in mangrove areas has become a global problem owing to its high toxic ability, abundance, and persistence [2,31,45]. However, approaches set up for monitoring and evaluating components are often absent, or methods adopted are mostly inconsistent between projects [82,83,84]. This reveals the significant lack of current agreement on holistic ecosystem health monitoring. There is therefore a pressing need for innovative approaches to improve the management of the remaining mangroves by developing a reliable and thorough assessment method for determining their ecosystem status and their likely future trajectory.
The simulation model presented in this paper is a first step towards the support of a balancing management approach involving both the phytoremediation role of mangroves by absorption of HM and its vulnerability by HM. The model is capable of analyzing quantitatively the tradeoff between pollutant removal and adverse effects on mangrove growth.
Unhealthy vegetation leads to the breakdown of important functions of the mangrove wetlands. The next step in the development of a decision support system is the use of remotely sensed data as the input for mangrove dynamic models (such as this of our model). Such consideration is not now a routine part of mangrove management planning and rehabilitation processes. While field work is more challenging in mangrove habitat, remotely sensed data is capable of mapping mangrove distribution, estimating biophysical parameters and characterizing ecosystem process. The combination of remotely sensed data and growth models can open the way to a data-and-model driven management of mangrove forest stability for coastal protection.

4. Conclusions

We have assessed the effects and interactions of HM pollutants and other environmental stressors on trunk biomass of R. apiculata forest in the Can Gio Mangrove forest and have shown that the tradeoff between using mangroves as a natural phytoremediation system and reserving this highly sensitive ecosystem is a challenge. It should be emphasized that the role of mono specific forest plantation here is particularly important as a long-term experiment allowing the assessment of tree development of a single species in different environments with well documented age of the stands. Data acquisition in natural forests with unknown tree ages demands much more efforts concerning the design of appropriate sampling schemes.
Our result may support mangrove ecosystem managers in their difficult task to protect this vulnerable ecosystem. In order to avoid further negative effects for the mangrove ecosystem of this protected area, the pollution load must be decreased in the future. Further studies are required to advance the understanding of the long-term trends and the complex interplay between contaminant concentrations, hydrodynamics, anthropogenic stressors and the mangrove ecosystem. Also, adequate solutions for management policies in cooperation with stakeholders from industry, administration, environmental protection and forestry based on a data-model driven network of monitoring systems need to be developed.

Author Contributions

A.N.: Conceptualization, methodology, software, original draft preparation, writing—reviewing and editing. O.R.: Supervision, data analysis, validation, writing—reviewing and editing. B.V.Q.L.: visualization, investigation, data curation, writing. N.T.K.P.: Laboratory sample analysis. K.C.D.: Data analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant number 105.99-2017.12.

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. MacFarlane, G.R.; Koller, C.E.; Blomberg, S.P. Accumulation and partitioning of heavy metals in mangroves: A synthesis of field-based studies. Chemosphere 2007, 69, 1454–1464. [Google Scholar] [CrossRef] [PubMed]
  2. Rahman, M.S.; Hossain, M.B.; Babu, S.M.O.F.; Rahman, M.; Ahmed, A.S.S.; Jolly, Y.N.; Choudhury, T.R.; Begum, B.A.; Kabir, J.; Akter, S. Source of metal contamination in sediment, their ecological risk, and phytoremediation ability of the studied mangrove plants in ship breaking area, Bangladesh. Mar. Pollut. Bull. 2019, 141, 137–146. [Google Scholar] [CrossRef] [PubMed]
  3. Wright, D.J.; Otte, M.L. Wetland plant effects on the biogeochemistry of metals beyond the rhizosphere. Biol. Environ. Proc. R. Ir. Acad. 1999, 99B, 3–10. [Google Scholar]
  4. Kamaruzzaman, B.Y.; Ong, M.C.; Jalal, K.C.A.; Shahbudin, S.; Nor, O.M. Accumulation of lead and copper in Rhizophora apiculata from Setiu mangrove forest, Terengganu, Malaysia. J. Environ. Biol. 2009, 30, 821. [Google Scholar] [PubMed]
  5. MacFarlane, G.R.; Pulkownik, A.; Burchett, M.D. Accumulation and distribution of heavy metals in the grey mangrove, Avicennia marina (Forsk.)Vierh.: Biological indication potential. Environ. Pollut. 2003, 123, 139–151. [Google Scholar] [CrossRef]
  6. Nazli, M.F.; Hashim, N.R. Heavy Metal Concentrations in an Important Mangrove Species, Sonneratia caseolaris, in Peninsular Malaysia. Environ. Asia 2010, 3, 50–55. [Google Scholar] [CrossRef]
  7. Chowdhury, R.; Favas, P.J.C.; Jonathan, M.P.; Venkatachalam, P.; Raja, P.; Sarkar, S.K. Bioremoval of trace metals from rhizosediment by mangrove plants in Indian Sundarban Wetland. Mar. Pollut. Bull. 2017, 124, 1078–1088. [Google Scholar] [CrossRef]
  8. Fengzhong, Z.; Liyu, H.; Wenjiao, Z. A primary study on adsorption of certain heavy metals on the litter leaf detritus of some mangrove species. J. Xiamen Univ. (Nat. Sci.) 1998, 37, 137–141. [Google Scholar]
  9. Wen-jiao, Z.; Xiao-yong, C.; Peng, L. Accumulation and biological cycling of heavy metal elements in Rhizophora stylosa mangroves in Yingluo Bay, China. Mar. Ecol. Prog. Ser. 1997, 159, 293–301. [Google Scholar]
  10. Titah, H.; Pratikno, H. Chromium Accumulation by Avicennia alba Growing at Ecotourism Mangrove Forest in Surabaya, Indonesia. J. Ecol. Eng. 2020, 21, 222–227. [Google Scholar] [CrossRef]
  11. Wang, S.; Yang, Z.; Xu, L. Mechanisms of copper toxicity and resistance of plants. Ecol. Environ. 2003, 3, 336–341. [Google Scholar]
  12. Sandilyan, S.; Kathiresan, K. Decline of mangroves—A threat of heavy metal poisoning in Asia. Ocean Coast. Manag. 2014, 102, 161–168. [Google Scholar] [CrossRef]
  13. Dubey, S.; Shri, M.; Gupta, A.; Rani, V.; Chakrabarty, D. Toxicity and detoxification of heavy metals during plant growth and metabolism. Environ. Chem. Lett. 2018, 16, 1169–1192. [Google Scholar] [CrossRef]
  14. Yan, Z.; Tam, N.F.Y. Effects of lead stress on anti-oxidative enzymes and stress-related hormones in seedlings of Excoecaria agallocha Linn. Plant Soil 2013, 367, 327–338. [Google Scholar] [CrossRef]
  15. Guo, X.Y.; Yan, C.L.; Ye, B.B. Effect of Cd–Zn Combined stress on contents of osmotic substances in Kandelia candel (L.) druce seedlings. Chin. J. Appl. Environ. Biol. 2009, 15, 308–312. [Google Scholar] [CrossRef]
  16. Li, P.; Zhao, C.; Zhang, Y.; Wang, X.; Wang, X.; Wang, J.; Wang, F.; Bi, Y. Calcium alleviates cadmium-induced inhibition on root growth by maintaining auxin homeostasis in Arabidopsis seedlings. Protoplasma 2016, 253, 185–200. [Google Scholar] [CrossRef]
  17. Lugo, A.E.; Snedaker, S.C. The Ecology of Mangroves. Annu. Rev. Ecol. Syst. 1974, 5, 39–64. [Google Scholar] [CrossRef]
  18. Thom, B.G. Mangrove ecology—A geomorphological perspective. In Mangrove Ecosystems in Australia; Clough, B.F., Ed.; Australian National University Press: Canberra, Australia, 1982; pp. 3–17. [Google Scholar]
  19. Feller, I.C.; Lovelock, C.E.; Berger, U.; McKee, K.L.; Joye, S.B.; Ball, M.C. Biocomplexity in Mangrove Ecosystems. Annu. Rev. Mar. Sci. 2009, 2, 395–417. [Google Scholar] [CrossRef] [Green Version]
  20. Twilley, R.R.; Rivera-Monroy, V.H. Developing performance measures of mangrove wetlands using simulation models of hydrology, nutrient biogeochemistry, and community dynamics. J. Coast. Res. 2005, 40, 79–93. [Google Scholar]
  21. Ewel, K.; Bourgeois, J.; Cole, T.; Zheng, S. Variation in environmental characteristics and vegetation in high-rainfall mangrove forests, Kosrae, Micronesia. Glob. Ecol. Biogeogr. Lett. 1998, 7, 49–56. [Google Scholar] [CrossRef]
  22. Chen, R.; Twilley, R.R. A gap dynamic model of mangrove forest development along gradients of soil salinity and nutrient resources. J. Ecol. 1998, 86, 37–51. [Google Scholar] [CrossRef]
  23. Grueters, U.; Ibrahim, M.R.; Satyanarayana, B.; Dahdouh-Guebas, F. Individual-based modeling of mangrove forest growth: MesoFON—Recent calibration and future direction. Estuar. Coast. Shelf Sci. 2019, 227, 106302. [Google Scholar] [CrossRef] [Green Version]
  24. Berger, U.; Hildenbrandt, H. A new approach to spatially explicit modelling of forest dynamics: Spacing, ageing and neighbourhood competition of mangrove trees. Ecol. Model. 2000, 132, 287–302. [Google Scholar] [CrossRef]
  25. Rideout, A.J.R.; Joshi, N.P.; Viergever, K.M.; Huxham, M.; Briers, R.A. Making predictions of mangrove deforestation: A comparison of two methods in Kenya. Glob. Chang. Biol. 2013, 19, 3493–3501. [Google Scholar] [CrossRef] [PubMed]
  26. Mandal, S.K.; Ray, R.; González, A.G.; Pokrovsky, O.S.; Mavromatis, V.; Jana, T.K. Accumulation, transport and toxicity of arsenic in the Sundarbans mangrove, India. Geoderma 2019, 354, 113891. [Google Scholar] [CrossRef]
  27. Richter, O.; Nguyen, H.A.; Nguyen, K.L.; Nguyen, V.P.; Biester, H.; Schmidt, P. Phytoremediation by mangrove trees: Experimental studies and model development. Chem. Eng. J. 2016, 294, 389–399. [Google Scholar] [CrossRef]
  28. Nguyen, A.; Le, B.V.Q.; Richter, O. The Role of Mangroves in the Retention of Heavy Metal (Chromium): A Simulation Study in the Thi Vai River Catchment, Vietnam. Int. J. Environ. Res. Public Health 2020, 17, 5823. [Google Scholar] [CrossRef]
  29. Cormier Salem, M.; Nguyen, T.; Ariadna, B.; Durand, J.-D.; Bettarel, Y.; Klein, J.; Hoang, H.; Panfili, J. The mangrove’s contribution to people: Interdisciplinary pilot study of the Can Gio Mangrove Biosphere Reserve in Viet Nam. Comptes Rendus Geosci. 2017, 349, 341–350. [Google Scholar] [CrossRef]
  30. Kuenzer, C.; Tuan, V.Q. Assessing the ecosystem services value of Can Gio Mangrove Biosphere Reserve: Combining earth-observation- and household-survey-based analyses. Appl. Geogr. 2013, 45, 167–184. [Google Scholar] [CrossRef]
  31. Costa-Boddeker, S.; Hoelzmann, P.; Thuyen, L.X.; Huy, H.D.; Nguyen, H.A.; Richter, O.; Schwalb, A. Ecological risk assessment of a coastal zone in Southern Vietnam: Spatial distribution and content of heavy metals in water and surface sediments of the Thi Vai Estuary and Can Gio Mangrove Forest. Mar. Pollut. Bull. 2017, 114, 1141–1151. [Google Scholar] [CrossRef]
  32. Costa-Boddeker, S.; Thuyen, L.X.; Hoelzmann, P.; de Stigter, H.C.; van Gaever, P.; Huy, H.D.; Schwalb, A. The hidden threat of heavy metal pollution in high sedimentation and highly dynamic environment: Assessment of metal accumulation rates in the Thi Vai Estuary, Southern Vietnam. Environ. Pollut. 2018, 242, 348–356. [Google Scholar] [CrossRef] [PubMed]
  33. Anh, N.H. A Model for Predicting Mangrove Forest Dynamics under Variable Environmental Conditions—A Case Study of the Estuary of Dongnai—Saigon River System, Vietnam; Technischen Universität Carolo-Wilhelmina: Braunschweig, Germany, 2011. [Google Scholar]
  34. Hang, H.T.; Anh, N.H. Study on using vegetation indicator in monitoring the change of soil features in Can Gio area with the help of remote sensing. J. Sci. Technol. Dev. Vietnam 2002, 5, 10–17. (In Vietnamese) [Google Scholar]
  35. Phan Nguyen, H. Effects of mangrove restoration and conservation on the biodiversity and environment in Can Gio District. In Proceedings of the Mangrove Management and Conservation Workshop, Okinawa, Japan, 26–30 March 2000; pp. 111–137. [Google Scholar]
  36. van Loon, A.F.; Dijksma, R.; van Mensvoort, M.E.F. Hydrological classification in mangrove areas: A case study in Can Gio, Vietnam. Aquat. Bot. 2007, 87, 80–82. [Google Scholar] [CrossRef]
  37. Nguyen, H.A.; Richter, O.; Huynh, D.H.; Nguyen, K.L.; Kolb, M.; Nguyen, V.P.; Bao, T.T. Accumulation of contaminants in mangrove species Rhizophora apiculate along Thi Vai River in the South of Vietnam. In EWATEC-COAST: Technologies for Environmental and Water Protection of Coastal Regions in Vietnam: Contributions to 4th International Conference for Environment and Natural Resources—ICENR 2014; Meon, G., Pätsch, M., Phuoc, N.V., Quan, N.H., Eds.; Cuvillier Verlag: Göttingen, Germany, 2014. [Google Scholar]
  38. Costa-Böddeker, S.; Thuyên, L.X.; Hoelzmann, P.; de Stigter, H.C.; van Gaever, P.; Huy, H.Đ.; Smol, J.P.; Schwalb, A. Heavy metal pollution in a reforested mangrove ecosystem (Can Gio Biosphere Reserve, Southern Vietnam): Effects of natural and anthropogenic stressors over a thirty-year history. Sci. Total Environ. 2020, 716, 137035. [Google Scholar] [CrossRef] [PubMed]
  39. Schwarzer, K.; Nguyen, T.; Ricklefs, K. Sediment re-deposition in the mangrove environment of Can Gio, Saigon River estuary (Vietnam). J. Coast. Res. 2016, SPI 75, 138–142. [Google Scholar] [CrossRef] [Green Version]
  40. ISO. 12914:2012 Soil Quality—Microwave Assisted Extraction of the Aqua Regia Soluble Fraction for the Determination of Elements; International Organization for Standardization: Geneva, Switzerland, 2012. [Google Scholar]
  41. ISO. 11466:1995 Soil Quality—Extraction of Trace Elements Soluble in Aqua Regia; International Organization for Standardization: Geneva, Switzerland, 1995. [Google Scholar]
  42. Botkin, D.B. Forest Dynamics: An Ecological Model; Oxford University Press: New York, NY, USA, 1993. [Google Scholar]
  43. Poschenrieder, C.; Cabot, C.; Martos, S.; Gallego Páramo, B.; Barceló, J. Do toxic ions induce hormesis in plants? Plant Sci. Int. J. Exp. Plant Biol. 2013, 212, 15–25. [Google Scholar] [CrossRef]
  44. Calabrese, E.J.; Blain, R.B. Hormesis and plant biology. Environ. Pollut. 2009, 157, 42–48. [Google Scholar] [CrossRef]
  45. Sinex, S.A.; Wright, D.A. Distribution of trace metals in the sediments and biota of Chesapeake Bay. Mar. Pollut. Bull. 1988, 19, 425–431. [Google Scholar] [CrossRef]
  46. Yongming, H.; Peixuan, D.; Junji, C.; Posmentier, E.S. Multivariate analysis of heavy metal contamination in urban dusts of Xi’an, Central China. Sci. Total Environ. 2006, 355, 176–186. [Google Scholar] [CrossRef]
  47. Zhang, J.; Liu, C.L. Riverine Composition and Estuarine Geochemistry of Particulate Metals in China—Weathering Features, Anthropogenic Impact and Chemical Fluxes. Estuar. Coast. Shelf Sci. 2002, 54, 1051–1070. [Google Scholar] [CrossRef]
  48. Taylor, S.R.; McLennan, S.M. The geochemical evolution of the continental crust. Rev. Geophys. 1995, 33, 241–265. [Google Scholar] [CrossRef]
  49. Wang, Z.-H.; Feng, J.; Nie, X.-P. Recent environmental changes reflected by metals and biogenic elements in sediments from the Guishan Island, the Pearl River Estuary, China. Estuar. Coast. Shelf Sci. 2015, 164, 493–505. [Google Scholar] [CrossRef]
  50. Nam, V.N.; Thuy, N.S.; Thanh, C.V.; Ha, H.S.; Tra, V.L.; Hoan, H.D. Development of Allometric Relationship between Biomass and Growth Params to Build Up a Reference Table for Mangrove Silviculture Practice and Management; Technical Report; Department of Agriculture and Rural Development of Ho Chi Minh City: Ho Chi Minh City, Vietnam, December 2004; p. 74. (In Vietnamese)
  51. Peterson, G.; Cumming, G.; Carpenter, S. Scenario Planning: A Tool for Conservation in an Uncertain World. Conserv. Biol. 2003, 17, 358–366. [Google Scholar] [CrossRef] [Green Version]
  52. Nguyen, H.A.; Richter, O.; Nguyen, V.P. Multiphysics Modeling of Pollutant Uptake by Mangroves. In Modeling, Simulation and Optimization of Complex Processes HPSC 2015; Springer: Cham, Switzerland, 2017; pp. 127–138. [Google Scholar]
  53. Ong Che, R.G. Concentration of 7 Heavy Metals in Sediments and Mangrove Root Samples from Mai Po, Hong Kong. Mar. Pollut. Bull. 1999, 39, 269–279. [Google Scholar] [CrossRef]
  54. Pratikno, H.; Titah, H.S.; Handayanu, H.; Harnani, B.R.D. The ability of Avicennia marina (Api-api putih) to Uptake Heavy Metal of Chromium at Wonorejo Coastal in Surabaya. In Proceedings of the International Conference on Coastal and Delta Areas, Semarang, Indonesia, 26 September 2017; pp. 164–169. [Google Scholar]
  55. Alharbi, O.M.L.; Khattab, R.A.; Ali, I.; Binnaser, Y.S.; Aqeel, A. Assessment of heavy metals contamination in the sediments and mangroves (Avicennia marina) at Yanbu coast, Red Sea, Saudi Arabia. Mar. Pollut. Bull. 2019, 149, 110669. [Google Scholar] [CrossRef]
  56. Yan, Z.; Sun, X.; Xu, Y.; Zhang, Q.; Li, X. Accumulation and Tolerance of Mangroves to Heavy Metals: A Review. Curr. Pollut. Rep. 2017, 3, 302–317. [Google Scholar] [CrossRef]
  57. Tam, N.F.Y.; Wong, Y.S. Retention of nutrients and heavy metals in mangrove sediment receiving wastewater of different strengths. Environ. Technol. 1993, 14, 719–729. [Google Scholar] [CrossRef]
  58. Lacerda, L.D.; Carvalho, C.E.V.; Tanizaki, K.F.; Ovalle, A.R.C.; Rezende, C.E. The Biogeochemistry and Trace Metals Distribution of Mangrove Rhizospheres. Biotropica 1993, 25, 252–257. [Google Scholar] [CrossRef]
  59. Carlos Augusto, R.S.; Lacerda, L.D.; Rezende, C.E. Metals Reservoir in a Red Mangrove Forest. Biotropica 1990, 22, 339–345. [Google Scholar] [CrossRef]
  60. Stigliani, W.M. Global Perspectives and Risk Assessment. In Biogeodynamics of Pollutants in Soils and Sediments: Risk Assessment of Delayed and Non-Linear Responses; Salomons, W., Stigliani, W.M., Eds.; Springer: Berlin/Heidelberg, Germany, 1995; pp. 331–343. [Google Scholar]
  61. Acosta, J.A.; Jansen, B.; Kalbitz, K.; Faz, A.; Martínez-Martínez, S. Salinity increases mobility of heavy metals in soils. Chemosphere 2011, 85, 1318–1324. [Google Scholar] [CrossRef]
  62. Pezeshki, S.R.; Delaune, R.D.; Patrick, W.H. Flooding and saltwater intrusion: Potential effects on survival and productivity of wetland forests along the U.S. Gulf Coast. For. Ecol. Manag. 1990, 33–34, 287–301. [Google Scholar] [CrossRef]
  63. Ye, Y.; Tam, N.F.-Y.; Lu, C.-Y.; Wong, Y.-S. Effects of salinity on germination, seedling growth and physiology of three salt-secreting mangrove species. Aquat. Bot. 2005, 83, 193–205. [Google Scholar] [CrossRef]
  64. Reef, R.; Feller, I.; Lovelock, C. Nutrition of mangroves. Tree Physiol. 2010, 30, 1148–1160. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Dai, Z.; Trettin, C.C.; Frolking, S.; Birdsey, R.A. Mangrove carbon assessment tool: Model development and sensitivity analysis. Estuar. Coast. Shelf Sci. 2018, 208, 23–35. [Google Scholar] [CrossRef]
  66. Ball, M.C.; Pidsley, S.M. Growth Responses to Salinity in Relation to Distribution of Two Mangrove Species, Sonneratia alba and S. lanceolata, in Northern Australia. Funct. Ecol. 1995, 9, 77–85. [Google Scholar] [CrossRef]
  67. Takemura, T.; Hanagata, N.; Sugihara, K.; Baba, S.; Karube, I.; Dubinsky, Z. Physiological and biochemical responses to salt stress in the mangrove, Bruguiera gymnorrhiza. Aquat. Bot. 2000, 68, 15–28. [Google Scholar] [CrossRef]
  68. Suárez, N.; Medina, E. Salinity effect on plant growth and leaf demography of the mangrove, Avicennia germinans L. Trees 2005, 19, 722. [Google Scholar] [CrossRef]
  69. Nguyen, K.L.; Nguyen, H.A.; Richter, O.; Pham, M.T.; Nguyen, V.P. Ecophysiological responses of young mangrove species Rhizophora apiculata (Blume) to different chromium contaminated environments. Sci. Total Environ. 2017, 574, 369–380. [Google Scholar] [CrossRef]
  70. Tri, N.H.; Hong, P.N.; Cuc, L.T. Can Gio Mangrove Biosphere Reserve; Cartographic Publishing House: Hanoi, Vietnam, 2000. [Google Scholar]
  71. Alongi, D.M. Present state and future of the world’s mangrove forests. Environ. Conserv. 2002, 29, 331–349. [Google Scholar] [CrossRef] [Green Version]
  72. Gilman, E.L.; Ellison, J.; Duke, N.C.; Field, C. Threats to mangroves from climate change and adaptation options: A review. Aquat. Bot. 2008, 89, 237–250. [Google Scholar] [CrossRef]
  73. Ellison, A.M.; Farnsworth, E.J. Simulated sea level change alters anatomy, physiology, growth, and reproduction of red mangrove (Rhizophora mangle L.). Oecologia 1997, 112, 435–446. [Google Scholar] [CrossRef] [PubMed]
  74. Zhang, J.-E.; Liu, J.-L.; Ouyang, Y.; Liao, B.-W.; Zhao, B.-L. Physiological responses of mangrove Sonneratia apetala Buch-Ham plant to wastewater nutrients and heavy metals. Int. J. Phytoremediat. 2011, 13, 456–464. [Google Scholar] [CrossRef]
  75. MacFarlane, G.R.; Burchett, M.D. Toxicity, growth and accumulation relationships of copper, lead and zinc in the grey mangrove Avicennia marina (Forsk.) Vierh. Mar. Environ. Res. 2002, 54, 65–84. [Google Scholar] [CrossRef]
  76. Chawla, A.; Patel, K. Effect of Manganese and Copper stress on a leaf biochemical parameter of a mangrove Avicennia marina (Forsk.) Vierh. Eur. J. Biomed. Pharm. Sci. 2018, 5, 250–253. [Google Scholar]
  77. Oliveira, H. Chromium as an Environmental Pollutant: Insights on Induced Plant Toxicity. J. Bot. 2012, 2012, 375843. [Google Scholar] [CrossRef]
  78. Merchant, S.S. The Elements of Plant Micronutrients. Plant Physiol. 2010, 154, 512. [Google Scholar] [CrossRef] [Green Version]
  79. Ahmad, M.S.A.; Ashraf, M. Essential Roles and Hazardous Effects of Nickel in Plants. In Reviews of Environmental Contamination and Toxicology; Whitacre, D.M., Ed.; Springer: New York, NY, USA, 2011; pp. 125–167. [Google Scholar]
  80. Vinh, T.V.; Marchand, C.; Linh, T.V.K.; Vinh, D.D.; Allenbach, M. Allometric models to estimate above-ground biomass and carbon stocks in Rhizophora apiculata tropical managed mangrove forests (Southern Viet Nam). For. Ecol. Manag. 2019, 434, 131–141. [Google Scholar] [CrossRef]
  81. Hoan, H.D. Determine on the capacity of carbon accumulation of Rhizophora apiculata Blume plantation forests in Can Gio Mangrove Biosphere Reserve, Ho Chi Minh City. Ph.D. Thesis, Vietnamese Academy of Forest Sciences, Hanoi, Vietnam, January 2019. [Google Scholar]
  82. Bakhtiyari, M.; Lee, S.Y.; Warnken, J. Seeing the forest as well as the trees: An expert opinion approach to identifying holistic condition indicators for mangrove ecosystems. Estuar. Coast. Shelf Sci. 2019, 222, 183–194. [Google Scholar] [CrossRef]
  83. Sutula, M.A.; Stein, E.D.; Collins, J.N.; Fetscher, A.E.; Clark, R. A Practical Guide for the Development of a Wetland Assessment Method: The California Experience1. JAWRA J. Am. Water Resour. Assoc. 2006, 42, 157–175. [Google Scholar] [CrossRef]
  84. Sheaves, M.; Brookes, J.; Coles, R.; Freckelton, M.; Groves, P.; Johnston, R.; Winberg, P. Repair and revitalisation of Australia’s tropical estuaries and coastal wetlands: Opportunities and constraints for the reinstatement of lost function and productivity. Mar. Policy 2014, 47, 23–38. [Google Scholar] [CrossRef]
Figure 1. The Can Gio mangrove area and the mensuration plots. Thi Vai catchment (inside the red frame) is affected by the industrial activities on the East bank and the Can Gio area is directly influenced by activities from the center of Ho Chi Minh city. The round yellow points are the measurement plots without taking samples. The triangle blue points are the measuring plus sampling plots for heavy metals (HMs) content determination.
Figure 1. The Can Gio mangrove area and the mensuration plots. Thi Vai catchment (inside the red frame) is affected by the industrial activities on the East bank and the Can Gio area is directly influenced by activities from the center of Ho Chi Minh city. The round yellow points are the measurement plots without taking samples. The triangle blue points are the measuring plus sampling plots for heavy metals (HMs) content determination.
Ijerph 17 09131 g001
Figure 2. Spatial distribution of input data for the growth model, including an elevation map (A), a salinity map (B), a Cr content distribution map (C), and a map of Rhizophora apiculata forest (D).
Figure 2. Spatial distribution of input data for the growth model, including an elevation map (A), a salinity map (B), a Cr content distribution map (C), and a map of Rhizophora apiculata forest (D).
Ijerph 17 09131 g002
Figure 3. Growth rate of the trees in the plots which are located on large salinity gradients and have the same range of tree density. The red points and curves are the data plots and the section and 95% confidence bands constructed at a density of 80 trees per plot, respectively. The blue points and curves are the data plots and the section and 95% confidence bands constructed at a density of 15 trees per plot, respectively.
Figure 3. Growth rate of the trees in the plots which are located on large salinity gradients and have the same range of tree density. The red points and curves are the data plots and the section and 95% confidence bands constructed at a density of 80 trees per plot, respectively. The blue points and curves are the data plots and the section and 95% confidence bands constructed at a density of 15 trees per plot, respectively.
Ijerph 17 09131 g003
Figure 4. Fitted curve of model growth rate G (the red curve) of R. apiculata and the G values extracted from data (the round points) along elevation gradient.
Figure 4. Fitted curve of model growth rate G (the red curve) of R. apiculata and the G values extracted from data (the round points) along elevation gradient.
Ijerph 17 09131 g004
Figure 5. Tree growth rate (G) is dependent on concentration (mg/kg) of the contaminants Ni (a), Cu (b), and Cr (c). This is a fit to the multiplier function for metal contaminants (Equation (8)). (d) Response surface of G in dependence on the total root concentration of Ni, Cu, and Cr and on the soil salinity.
Figure 5. Tree growth rate (G) is dependent on concentration (mg/kg) of the contaminants Ni (a), Cu (b), and Cr (c). This is a fit to the multiplier function for metal contaminants (Equation (8)). (d) Response surface of G in dependence on the total root concentration of Ni, Cu, and Cr and on the soil salinity.
Ijerph 17 09131 g005
Figure 6. Measured dbh values vs. model predictions for the validation data. The coefficient of variation value is R2 = 0.87.
Figure 6. Measured dbh values vs. model predictions for the validation data. The coefficient of variation value is R2 = 0.87.
Ijerph 17 09131 g006
Figure 7. 95% confidence surfaces together with the G values of the validation data set.
Figure 7. 95% confidence surfaces together with the G values of the validation data set.
Ijerph 17 09131 g007
Figure 8. Distribution patterns of trunk biomass over time in different scenarios. The columns show the biomasses at different times, the rows show the results for the three scenarios.
Figure 8. Distribution patterns of trunk biomass over time in different scenarios. The columns show the biomasses at different times, the rows show the results for the three scenarios.
Ijerph 17 09131 g008
Figure 9. Biomass (kg/ha) histogram of three scenarios in the last phase of simulations (in 2020).
Figure 9. Biomass (kg/ha) histogram of three scenarios in the last phase of simulations (in 2020).
Ijerph 17 09131 g009
Table 1. Contents (value, max value, mean value, min value and standard deviation std) of the heavy metals Cu, Cr, and Ni (mg/kg) in soil and in tissues of the mangrove species R. apiculata at the 11 sampling positions.
Table 1. Contents (value, max value, mean value, min value and standard deviation std) of the heavy metals Cu, Cr, and Ni (mg/kg) in soil and in tissues of the mangrove species R. apiculata at the 11 sampling positions.
SOILROOTLEAF
CuCrNiCuCrNiCuCrNi
CGUpstreamNP180.044.963.310.211.060.570.311.510.71
DK470.823.762.720.600.700.450.720.370.12
DK664.514.352.630.310.320.110.250.290.01
DK765.654.443.080.540.260.170.390.400.09
DK875.764.283.250.371.200.660.180.230.01
DownstreamNP276.624.493.060.470.320.220.361.270.46
TVUpstreamDN162.4112.398.493.074.464.903.530.390.14
DN2134.0319.8411.358.1711.698.542.120.490.01
M560.4021.1412.015.552.192.194.031.040.27
M379.4613.628.824.881.521.303.402.051.82
DownstreamM182.5114.755.823.565.243.087.831.080.76
CG max80.044.963.310.601.200.660.721.510.71
mean72.234.383.010.420.640.360.370.680.23
min64.513.762.630.210.260.110.180.230.01
std6.290.390.280.150.410.230.190.560.29
TV max134.0321.1412.018.1711.698.547.832.051.82
mean83.7616.359.305.055.024.004.181.010.60
min82.5114.755.823.565.243.087.831.080.76
std29.793.902.482.014.032.872.160.660.74
Note: bold: to press the ID of the positions.
Table 2. Values of tree growth parameters and environmental factors at the 11 measurement plots.
Table 2. Values of tree growth parameters and environmental factors at the 11 measurement plots.
PositionMean dbh (cm)Mean Height (cm)Density/Plot (10 × 10 sq·m)Age (Year)Porosity Salinity (ppt)Mean Tree Growth Rate (G)EF (Cu)EF (Cr)EF (Ni)
Can GioNP1Upstream18.621989.9726268.0390.533.410.110.09
DK4 10.461232.91423212.4277.084.130.110.10
DK6 22.292293.7883812.9275.462.530.090.07
DK7 15.061673.72133715.5201.612.420.080.08
DK8 19.712082.55114015.5197.742.620.070.07
NP2Downstream13.141501.13163918.2166.913.290.100.09
Thi VaiDN1Upstream16.131770.9213258.0332.131.180.120.11
DN2 16.161757.5120259.6331.823.490.260.20
M5 13.991574.53113111.3317.361.830.320.24
M3 11.971381.85202313.5300.982.530.220.19
M1Downstream12.451428.23173217.1152.973.900.350.18
Note: bold italic numbers are the EF values larger than 1.5.
Table 3. Range of heavy metal concentrations in the soil, mangroves leaf and root of the Can Gio Mangrove Forest and of the other mangroves worldwide.
Table 3. Range of heavy metal concentrations in the soil, mangroves leaf and root of the Can Gio Mangrove Forest and of the other mangroves worldwide.
LocationSpeciesPartCr (mg/kg)Cu (mg/kg)Ni (mg/kg)Reference
This study, Can Gio, Vietnam (2018)SoilSoil3.76–21.1460.4–134.032.63–12.01This study
Rhizophora apiculataLeaves0.23–2.050.18–7.830.005–1.82
Roots0.26–11.70.21–8.20.11–8.54
Can Gio, Vietnam (2013)SoilSoil86–24119.7–32.756.6–247[37]
Rhizophora apiculataLeaves1.48–23.62.77–6.411.4–25.7
Roots2.14–90.60.746–18.31.79–62.9
Can Gio, Vietnam (2002–2012)Sediment coreSediment107.70–208.8026.74–82.3256.25–82.99[32]
Note:
Thi Vai side, upstream
Can Gio, Vietnam (2002–2012)Sediment coreSediment51.6–82.511.5–38.324.7–46.5[32]
Note: Thi Vai side, downstream
Can Gio, Vietnam (1977–2011)Sediment coreSediment27.1–71.57.1–27.011.7–56.3[38]
Note: Can Gio side
Mai po, Hong KongSoilSoil20–74.651.1–87.443.9–86.9[53]
Acanthus ilicifoliusRoots1.6–6.225.2–65.48–20.1
Aegicerus corniculatumRoots1.8–5.822.4–46.26–32.01
Kandelia candelRoots1.8–6.419.6–29.44–28.001
Surabaya, IndonesiaSoilSoil47–79.3 [10,54]
Avicennia albaRoots25.4–55.3
Avicennia marinaRoots28–92.25
Yanbu, Red Sea, Saudi ArabiaSoilSoil14.9–28917.2–217.227.3–241.8[55]
Avicennia marinaLeaves14.2–50.118.1–40.216.1–56.3
Roots16.3–40.516.8–37.317.2–38.2
Note: bold and italic text and numbers are values of this study.
Table 4. Result of the parameter estimation.
Table 4. Result of the parameter estimation.
ParameterValue
Density multiplier a 0 1 0.402741
d 1 0.156347
t r 1 (trees/100 m2) 66.9241
Salinity multiplier a 0 2 0.21276
d 2 −0.4
t r 2 (ppt) 11.348
Elevation multiplier a 1 e 0.027
e l 1 −0.96
α 13.1
e l 2 1.48
β 6.02
a 2 e 0.11
Pollutant multiplier CuCrNi
t h 1 3.281.890.53
t h 2 54.47170.25105.43
α 1 0.680.120.06
α 2 0.811.692.53
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nguyen, A.; Richter, O.; Le, B.V.Q.; Phuong, N.T.K.; Dinh, K.C. Long-Term Heavy Metal Retention by Mangroves and Effect on Its Growth: A Field Inventory and Scenario Simulation. Int. J. Environ. Res. Public Health 2020, 17, 9131. https://doi.org/10.3390/ijerph17239131

AMA Style

Nguyen A, Richter O, Le BVQ, Phuong NTK, Dinh KC. Long-Term Heavy Metal Retention by Mangroves and Effect on Its Growth: A Field Inventory and Scenario Simulation. International Journal of Environmental Research and Public Health. 2020; 17(23):9131. https://doi.org/10.3390/ijerph17239131

Chicago/Turabian Style

Nguyen, Anh, Otto Richter, Bao V.Q. Le, Nguyen Thi Kim Phuong, and Kim Chi Dinh. 2020. "Long-Term Heavy Metal Retention by Mangroves and Effect on Its Growth: A Field Inventory and Scenario Simulation" International Journal of Environmental Research and Public Health 17, no. 23: 9131. https://doi.org/10.3390/ijerph17239131

APA Style

Nguyen, A., Richter, O., Le, B. V. Q., Phuong, N. T. K., & Dinh, K. C. (2020). Long-Term Heavy Metal Retention by Mangroves and Effect on Its Growth: A Field Inventory and Scenario Simulation. International Journal of Environmental Research and Public Health, 17(23), 9131. https://doi.org/10.3390/ijerph17239131

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop