Global temperature is predicted to increase in the end of the century and one of the primary consequences of this warming is the sea level rise. Considering the vulnerabilities on coastal systems and water resources, it is important to evaluate the potential effects of this rising in coastal areas, since the saline intrusion on rivers would be intensified, leading to problems related to water quality. In this context, the present work aimed to verify saline intrusion changes along an important river, S ão Francisco Canal, located in Rio de Janeiro State, Brazil. For this purpose, a hydrodynamic modeling was performed using SisBaHiA, considering different sea levels and tide conditions. According to the results, it was verified the intensification on saline intrusion and higher salinity values due to a sea level rise of 0.5 m. These results show that new licenses for water withdrawals must be carefully analyzed as the fluvial flow plays an important role to contain the saltwater intrusion on the studied river. Accordingly, it is recommended the evaluation of climate change effects in order to choose best strategies to reduce coastal vulnerability, and the use of this theme on environmental licensing and territorial planning, integrating water planning with coastal management.
Climate change and its effects on the environment and human life have been increasingly discussed and considered in the present society. Despite the increasing on greenhouse gas (GHG) emissions after Industrial Revolution, the causes of the global warming are not straightforward. Nevertheless, even natural or hu- man-caused, the climate change effects should be considered and investigated in different scales and spheres both for global and regional planning.
This increasing on global surface temperature has already been observed in the present days. According to IPCC (Intergovernmental Panel on Climate Change), the global temperature increased 0.74˚C along the 20th century and its effects could also be observed on the global hydrological cycle [
These projections are an important issue, mainly considering the coastal vulnerability to sea level rise in addition to the effects of potential changes on winds, waves and on the frequency of climate extremes. Among others, the SL rise leads to problems associated to coastal flooding and contamination of coastal aquifers. However, these effects vary according to regional characteristics being important to study different environments and regions [
In Brazil, just few studies have been realized considering climate change effects on coastal zones. Despite of this, in some regions, comprising the coastal zone of the state of Rio de Janeiro, losses associated to extreme events are being experienced. These current evidences of SL rising effects on coastal environments demonstrate that further investigation is a necessary and strategic factor [
According to global climate models, the average prediction for global SL rise along the 21st century is about 0.47 m in relation to the mean SL for 1986-2005. For 2100, this change is projected in 0.53 m to the same reference period [
This projected rise would deteriorate water quality, since saline wedge reaches inshore regions of the coastal aquifers. Saline intrusion is a serious environmental condition as 80% of the world’s population live on coastal zones and use local aquifers for their water supply [
In some regions, saline intrusion on coastal aquifers may be intensified by lower rainfall rates, what decreases the water availability for human supply, as previously observed in Senegal [
São Francisco Canal is one of the main aquifers located in Rio de Janeiro State, Brazil (
To verify the saline intrusion on the studied region, hydrodynamic and dispersion modeling were performed using SisBaHia (Sistema Base de Hidrodinâmica Ambiental) system. This system corresponds to a finite element hydrodynamic numerical model governed by barotropic forcing and optimized for natural water bodies where the effects of variable density can be neglected. This model is based on principle of conservation of momentum what, in conjunction with the continuity equation, the equation of state and the respective constituent transport equations, constitute the fundamental mathematical model for any water body [
Eulerian transport models are widely used on environmental analysis and are suitable to simulate the transport of dissolved scalars in the water. They simulate advective and diffusive transports and represent decay or increasing rates of the target property. In the present study, this model was coupled to the hydrodynamic model in order to verify the salinity dispersion in the study area. This coupling is necessary since the salinity is an active tracer and it can interfere on the hydrodynamics [
The salinity dispersion model is depth integrated and defined as in:
where
The hydrodynamic and dispersion simulations were realized in an unstructured grid with a minimum resolution of 280 m on the open boundary. The local bathymetry was interpolated in the numerical grid using information from the Brazilian Navy nautical charts. The numerical grid and the local bathymetry are presented in
Twelve hydrodynamic scenarios were used to evaluate the saline intrusion in the studied river. These scenarios comprise different combination of astronomical tide conditions, storm surge amplitudes and the rise of 0.5 m on the SL, as shown in
Scenario | Astronomical tide | Storm surge (m) | SL (m) |
---|---|---|---|
1 | Spring | - | - |
2 | Neap | - | - |
3 | Spring | 0.40 | - |
4 | Neap | 0.40 | - |
5 | Spring | 0.80 | - |
6 | Neap | 0.80 | - |
7 | Spring | - | 0.5 |
8 | Neap | - | 0.5 |
9 | Spring | 0.40 | 0.5 |
10 | Neap | 0.40 | 0.5 |
11 | Spring | 0.80 | 0.5 |
12 | Neap | 0.80 | 0.5 |
Tides were included as boundary conditions based on the predictions for the tidal constituents from Itaguaí Port tidal station. Along the grid boundaries, this forcing presented a spatial variation by the use of a phase lag that was considered in function of the distance and the average depth.
In addition to the tidal forcing and the particular conditions for each scenario, the model was forced by hourly surface winds, measured in a meteorological station from INMET (Instituto Nacional de Meteorologia) located in Angra dos Reis (at 23.00˚S and 44.32˚W). The salinity values used as initial and boundary conditions were measured by the State Environmental Agency (INEA) during its monitoring campaigns along Sepetiba Bay.
The fluvial flow of some rivers including Saco do Engenho (0.5 m3/s), Mazomba River (1.4 m3/s), Itaguaí River (0.1 m3/s), Guarda River (5.5 m3/s), Santo Agostinho Canal (28.4 m3/s), Guandú Canal (2.8 m3/s) and Itá Canal (3.2 m3/s) were used in the modeling. The average river flows were considered as permanent and used as boundary conditions.
Water volume releases and uptakes by the main water users, mainly corresponding to energy and siderurgy sectors, were also included as sinks and sources of water volume. The sum of all uptakes in the modeled region of São Francisco Canal corresponds to 53.9 m3/s and the total release of water is realized to Guandú Canal corresponding to 20.1 m3/s.
Each experiment was performed using an integration time of 15 days. The experiment results were used to compare the differences between scenarios with SL rise and with the present SL.
In this section, the hydrodynamic and the salinity dispersion results are presented aiming to verify the main changes observed after the consideration of a SL rise due to a future climate scenario.
The Sepetiba bay is a semi-enclosed body and its hydrodynamics basically depends on the action of tides, winds and river flows. Because of this, they are the same factors influencing the water level on São Francisco Canal. The variations on the water level observed in the canal were verified with model results. In
Changes on the velocity fields in response to the increasing sea level due to climate change, tide condition and storm surges were observed along the canal. In
spring tide situations. Scenario 11 corresponds to the combined effect of mean SL rise and the presence of a storm surge event with 0.8 m of amplitude, and the second greatest intrusion was observed on the scenario that considers the mean SL rise in conjunction with the occurrence of storm surge event with 0.4 m of amplitude.
However, the third maximum intrusion was observed in the results obtained on Scenario 5, not on Scenario 7, showing, in this case, that the effect of the storm surge (0.80 m) during spring tide periods suppresses the effect of the mean SL change. Nevertheless, the impact of SL rise (Scenario 7) was greater than the impact with the scenario which just considers the effect of a storm surge with 0.4 m of amplitude (Scenario 3).
With the Scenario 1, that is the mildest of spring tide scenarios, the position of the upstream currents is near to the position verified on Scenario 12 results that represent a condition of neap tide with SL rise and 0.8 m of storm surge. Among the other neap tide scenarios, it was maintained the same relationships observed between spring tide scenarios.
In general, during neap tides, smaller impacts are observed on water levels and on upstream currents than during spring tide periods. However, considering worsening factors, as the SL rise and the presence of storm surge with higher amplitudes, the sea water from Sepetiba bay advances along the canal.
The differences on the salinity concentrations in São Francisco Canal are presented by the average and maximum values computed in each scenario. Considering the extreme values obtained with the model results, the displacement of the mean SL was not main responsible for the higher values, as can be observed in
Considering the average salinity obtained in each scenario, in general higher salinity values were computed on SL scenarios (
On scenarios with the effect of storm surges with 0.8 m of height and the SL rise (Scenarios 11 and 12), concentrations above the national normative limits (0.5) for this type of water body were observed for longer periods (
The average salinity values for these two scenarios also indicate this pattern (
Both average salinity values and the maximum inward position of the salt wedge indicate that the worse water quality condition is obtained in response to the advance of the surface currents along the river. Therefore, higher saline intrusions are observed with the SL rise and during the spring tide, being more intense in accordance to the amplitude of storm surges.
There are three strategies to reduce coastal vulnerability. These strategies are: protection, when the risk is minimized due to the reduction on its occurrence probability; retreat, when the risk is reduced by the limitation on its potential effects; and accommodation, when the capacity of the society to recover itself from the effects is increased [
In the case of SL rise, protection would mean the reduction on GHG emissions, but, even if the emissions are stabilized in the next decades, a significant increase on SL should occur due to the temperature differences between oceanic layers [
Retreat means the abandonment of present occupations and improvements, and the migration toward the continent. It is suitable for new occupations, when the new land uses would be projected according to the SL rising. However, considering intensely occupied waterfront areas, the alternative would be the combination of retreat and accommodation strategies, adapting the coastal zone to the new SL [
According to Titus (1991), in established coastal areas, engineering interventions should be used, and, when these protection measurements are expensive, they should be replaced by accommodation ones. Retreat strategies, according to Titus, are only appropriate in underdeveloped areas, with few investments in infrastructure, or where is possible to occur a natural adaptation [
In Brazil, the National Policy on Climate Change (PNMC, from Portuguese), recognized by the Federal Law n. 12,187/2009, establishes principles to avoid, predict and minimize climate change due to anthropogenic causes, aiming to adopt adaptation measures and integrated mitigation and adaptation strategies in local, regional and national spheres.
In Rio de Janeiro, the Law n. 5690/2010 establishes principles, aims, guidelines and instruments to adapt and prevent climate change negative effects, similar, in general, to those on National Policy. Among these instruments, economic-eco- logical zoning is important in light of the SL rise. However, environmental licensing, as contained in the Law, is just related to the control of GHG emissions by license applicants and it does not state about the local vulnerabilities to climate changes in installation site and the possible effect of this licenses on this vulnerability. Both in National and State policies, these issues are not addressed.
The state of Rio de Janeiro has been frequently affected by extreme rainfall events. Even though the frequency of these events is not related to the global climate change, they show the possible consequences of these changes [
Observing the changes in response to the future SL rise, the flood currents reached more than 19 km along the São Francisco Canal considering the most extreme scenario. Higher water levels were also observed along the canal, of about 2.72 m, while under normal conditions this value only reaches 0.96 m during spring tides.
Comparing the results from hydrodynamic and transport models, the highest salinity concentrations are not observed in scenarios with higher water elevations, but in scenarios where the stream inversion limits occurred upstream on the canal.
The most extreme scenario (storm surge, spring tide and SL rise) presented an average increase of 52% on salinity values in relation to the milder scenario and of 22% in relation to the corresponding scenario (storm surge and spring tide) with no SL rise. Thus, it is possible to verify that an increase of 0.5 m in the local SL would increase the salinity in the San Francisco Canal, what would deteriorate water quality. This increase is mainly evidenced by the further advancement of the salt wedge, affecting upstream areas.
According to the results, the SL rise would lead the salt wedge 450 m upstream along the river, comparing the results between Scenarios 1 and 7. Considering the effect of storm surges with amplitudes of 0.40 and 0.80 m, during the periods of spring tide, the differences in response to the SL rise are of 380 and 267 m, respectively.
The saline intrusion obtained in scenario 11 that corresponds to the most extreme situation, is positioned 4.4 km and 0.8 km upstream to the positions obtained with the milder scenarios, representing the period of neap and spring tides, respectively. The maximum intrusion was observed at 5.4 km from the river mouth.
This same type of analysis was performed by [
Despite the saline wedge intrusion in São Francisco Canal is not as extensive as that of the Gorai River, salinity concentrations ranged approximately in accord to the proportion cited by the authors. Near to the mouth of São Francisco Canal, the salinity increased on the average in a rate of 1.4 per meter of SL rise.
It is necessary to note that other changes besides SL rise were not considered in the present case study. However, it is easy to assume that in a future scenario, environmental conditions will be different from the current scenario, especially considering the effects from climate changes. Thus, future bathymetric and morphological conditions, as well as winds and rainfall rates, and therefore, river flows, could cause significant changes in the results presented in this work. These considerations may lead to the worsening or the attenuation of the results.
Thus, comprehensive studies on the local conditions are necessary to allow the strategic planning of water resources in the studied river basin, considering the imminent interference of climate change on water availability. Moreover, previous studies on the effects of climate change show that the water management in coastal areas should be made considering the potential saline wedge intrusions.
However, such planning is a complicated achievement as it is needed to consider the best water extraction locations and to control the water quality [
The sustainable exploitation of aquifers, according to some researchers, may be the solution to prevent the degradation by saltwater intrusion. However, it does not seem feasible for aquifers in broad expansion. São Francisco Canal, because of its location, becomes strategic for industrial use, especially because it is located nearby Itaguaí Port, which attracts the establishment of a range of economic sectors in the region.
According to the strategic plan for the region of Guandu river basin, some criteria must be considered in the implementation of new raw water users [
In this context, new water transfers must be carefully evaluated. As Guandu river flow depends on the present transfer from Paraíba do Sul River, even possible new transfers of water from this source could pose risks in maintaining the required flow to control the saltwater intrusion along the São Francisco Canal. Additionally, new water uptakes in Paraíba do Sul River could possibly exacerbate adverse conditions considering future climate change scenarios, such as drought extremes in river headwaters.
Thus, water resource and coastal managements must be integrated and in accordance with the principles stated in climate change policies to ensure that the environmental quality will be maintained. In addition, to turn possible the maintenance of the current water withdrawals and the establishment of new grants, it is necessary a thorough study of climate variables and considering them in the licensing process, in order to do not compromise the water supply and to choose strategies to reduce the vulnerability of the coastal zone.
In the present work, the effect of the sea level rise on water resources was evaluated in a river located in the region of Sepetiba Bay, Rio de Janeiro, Brazil. Using numerical modeling, changes on water levels and on canal hydrodynamics were observed. The upstream currents during rising tides would be more intrusive due to the intensification of pressure gradients in response to future SL rise.
This is an important result as climate change on the numerical modeling was just included as the SL displacement. Then, the modeled changes correspond just to the response to the modification on the SL and do not comprise the combined effects in consequence of this rise, as the increasing on wave heights near to the shore, due to the decreased bottom friction, and the morphological changes, due to the intensification on transport processes.
Resulting from these modeled hydrodynamical changes, the water quality along the river would be injured, in which more areas influenced by saline waters would be observed so as the increasing on salinity values at the river points previously affected by saline intrusion. The rise of 0.5 m on the mean SL would lead to a salt wedge intrusion up to 5.4 km away from the river mouth, affecting water uptakes of some industries established in the region.
The present results illustrate how this coastal environment would respond to sea level changes, however, further investigations are needed to better understand the water availability in a future climate scenario. Nevertheless, the potential future degradation on São Francisco Canal water quality indicates that future water grants should be evaluated by the local environmental agencies after considering climate scenarios. The decreasing on fluvial flow in the studied river would lead to more intense saline intrusion.
Finally, it is recommended the consideration of climate change effects along the environmental licensing, in order to establish suited strategies for each type of project and for each implementation area, and the integration between water resources and coastal management. Thus, it would be possible to determine the most vulnerable areas where major investments are needed and to manage the coastal environment according to charts of climate change vulnerability.
The authors would like to thank the CNPq (National Council of Technological and Scientific Development) for the financial support.
Toste, R., Rosman, P.C.C. and de Freitas, M.A.V. (2017) Saline Intrusion Response to Sea Level Rise and Its Implications on Water and Coastal Ma- nagement: A Case Study in Brazil. Journal of Water Resource and Protection, 9, 510- 522. https://doi.org/10.4236/jwarp.2017.95033