Visualizing the spatial water quality of Bentota,Sri Lanka in the presence of seawater intrusion

Seawater flows towards the inlands along with the rivers and canals, through the process of infiltration and leaks in the ground water characterized by high concentrations of soluble salts. High salinity concentrations can make groundwater unsuitable for public consumption and surface water unsuitable for irrigation and agricultural activities. This study envisages the fluctuations of ground and surface water quality of Bentota area in the presence of seawater intrusion. The temporal and spatial variations of eleven water parameters were monitored by collecting the water samples during one year period. Spatial distributions were assessed by applying the Inverse distance weighted (IDW) interpolation method in Arc GIS 10.5 software. Water quality is assessed on the integration of all parameters in terms of an index based on the World Health Organization (WHO) standards. The significant linear relationship between the considered parameters of surface water (SW) and groundwater (GW) were identified applying correlation analysis using SPSS software. All parameters of surface water were above the permissible limits of WHO standards. Surface water quality index values with respect to 60% of canals show very poor quality (>1 250) of surface water indicating their unsuitability for irrigation activities. Those surface water bodies indicated very highly saline conditions during dry months. The spatial distribution of ground water quality index with respect to the highest parameter values of each sampling location indicates that 52.2% of total land extent of Bentota Divisional Secretariat Division (DSD) has good quality of ground water which is suitable for drinking. Its 47.2% of total land extent has poor quality of ground water for drinking purpose and less than 0.5% of the area consists of excellent or very poor quality of ground water in each. This study helps to manage coastal aquifers by understanding the extreme water quality conditions and coastal salinity.     

Visualizing the spatial water quality of Bentota,Sri Lanka in the presence of seawater intrusion

Seawater flows towards the inlands along with the rivers and canals, through the process of infiltration and leaks in the ground water characterized by high concentrations of soluble salts. High salinity concentrations can make groundwater unsuitable for public consumption and surface water unsuitable for irrigation and agricultural activities. This study envisages the fluctuations of ground and surface water quality of Bentota area in the presence of seawater intrusion. The temporal and spatial variations of eleven water parameters were monitored by collecting the water samples during one year period. Spatial distributions were assessed by applying the Inverse distance weighted (IDW) interpolation method in Arc GIS 10.5 software. Water quality is assessed on the integration of all parameters in terms of an index based on the World Health Organization (WHO) standards. The significant linear relationship between the considered parameters of surface water (SW) and groundwater (GW) were identified applying correlation analysis using SPSS software. All parameters of surface water were above the permissible limits of WHO standards. Surface water quality index values with respect to 60% of canals show very poor quality (>1 250) of surface water indicating their unsuitability for irrigation activities. Those surface water bodies indicated very highly saline conditions during dry months. The spatial distribution of ground water quality index with respect to the highest parameter values of each sampling location indicates that 52.2% of total land extent of Bentota Divisional Secretariat Division (DSD) has good quality of ground water which is suitable for drinking. Its 47.2% of total land extent has poor quality of ground water for drinking purpose and less than 0.5% of the area consists of excellent or very poor quality of ground water in each. This study helps to manage coastal aquifers by understanding the extreme water quality conditions and coastal salinity.     

Application of multi-element relationships in stream sediments to mineral exploration: a case study of Walawe Ganga Basin,Sri Lanka

This study on the application of multi-element relationships in stream sediments to mineral exploration in the Walawe Ganga Basin presents one of the first of its kind in Sri Lanka. In order to determine the nature of these sediments, multi-element analysis of selected grain size (<63 μm, 63–125 μm, 125–177 μm and 177–250 μm) fractions was undertaken. The chemical composition of these fractions were compared with those of the upper continental crust.Some elements, notably Zr, Hf, Th, U, Ce and La, show very high enrichment factors (e.g. 54 for Zr) as compared to the upper crustal abundance. These are presumably associated with heavy minerals such as zircon, rutile and monazite found in abundance in the stream sediments of the Walawe Ganga Basin.Principal Component Analysis (PCA) of the geochemical data show that multi-element relationships could be effectively used to delineate target areas for mineral exploration. With the use of PCA, this study reveals that areas with associations of calc-silicate/marble and charnockitic rocks are probable source regions for mineral occurrences, particularly in the axial regions of anticlines.     

Metals in a lateritic peat deposit — A case study from Sri Lanka

Peat is known to show metal enrichment. The geochemical enrichment factor (GEF) of U on peat humic acid is, for example, 10,000. Apart from U, peat accumulates a large number of metals including Cu, V, Ni, Cr, Zn, Pb, Zr, Mo, etc.

A large laterite bound peat deposit in Sri Lanka was studied for its distribution of metals. Most metals, in particular Be, Zr, Li and the majority of the transition series metals, exhibit a strong positive correlation with K, Al, Fe and Mn. A noteworthy feature observed was the poor correlation of the metals with organic carbon, as contrasted with a very good correlation with elemental components of clay. The laterites found around the peat deposit appear to be the possible source materials for the metals in the peat deposit. The geochemical distribution of metals in the peat is governed by:

1. (1) the nature of the organic matter;

2. (2) clay component;

3. (3) pH and Eh changes;

4. (4) marine-continental sediment mixing;

5. (5) nature of source materials for the metals.

Most of the metals appeared to be bound to the clay matrix, coated perhaps with hydrous Fe- and Mn-oxides. The existence of brine-rich marine-based sediments in the peat deposit had an important bearing on the geochemical distribution of the metals, as evidenced by the fact that in saline conditions most metals form chloride complexes, particularly under the prevailing acidic conditions. The metals thus mobilized interact with humic acids. The peat deposit under investigation provides an ideal setting for the study of the distribution of metals in a laterite-clay-peat mixture.     

High resolution local Moho determination using gravity inversion: A case study in Sri Lanka

The seismic data incorporated in global Moho models are sparse and therefore the interpolation of global Moho depths on a local area may give unrealistic results, especially in regions without adequate seismic information. Gravity inversion is a useful tool that can be used to determine Moho depths in the mentioned regions. This paper describes an interactive way of local Moho depth determination using the gravity inversion method constrained with available seismic data. Before applying inversion algorithms, the Bouguer gravity data is filtered in various stages that reduce the potential bias usually expected in Moho depth determination using gravity methods with constant density contrast assumption. A test area with reliable seismic data is used to validate the results of Moho computation, and subsequently the same computation procedure is applied to the Sri Lankan region. The results of the test area are in better agreement with seismically determined Moho depths than those obtained by global Moho models. In the Sri Lankan region, Moho determination reveals a fairly uniform thin crust of average thickness around 20 km. The overall result suggests that our gravity inversion method is robust and may be suitable for local Moho determination in virgin regions, especially those without sufficient seismic data.     

Dynamic characteristics of saltwater intrusion in the Pearl River Estuary,China

River discharge, tide, wind, topography and other factors all have great impacts on the saltwater intrusion of Modaomen Waterway (MW), a major outlet of the Pearl River Estuary. A coupled 1D–3D numerical model was applied in this study to account for the dynamic characteristics of saltwater intrusion in the MW, and the impacts of tide and river discharge on the length of saltwater intrusion were uncovered. Results are as the followings: (1) River discharge from upstream induces an obvious dilution of salinity along the MW, whereas tide can exert a positive force that pushes salt water landward. The effects of river discharge and tide on the length of saltwater intrusion can be well described by a regression function; (2) the saltwater intrusion along the MW is generally aggravated by increases in tidal range from the South China Sea. The length of saltwater intrusion usually reaches a maximum 2 or 3 days before spring tide, and the hourly length of saltwater intrusion along the MW usually slows the tidal process for approximately 4 h, which can provide important information that the pumping operation along the MW to store freshwater in the backup storages needs to be at least 3 days ahead of the spring tide so as to avoid serious impact from saltwater intrusion; (3) the length of saltwater intrusion generally decreases with increasing river discharge. In 2005, 2009 and 2010, the average river discharge from upstream was 2680, 2630 and 3160 m³/s, respectively, with corresponding average lengths of saltwater intrusion of 32.7, 42.3 and 21.4 km. The inverse correlation between the water flow and the length of saltwater intrusion may provide some guidance for operations to maintain enough upstream flow to dilute the salinity and therefore satisfy the domestic water supply.     

The relationship between soil geochemistry and die back of montane forests in Sri Lanka: a case study

Tropical montane forests of Sri Lanka form a unique ecosystem with more than 50% of endemic plant species. It has been noted that trees, belonging to different size and age classes of these forest, have been dying due to a yet unknown factor. This phenomenon was first observed in the Horton Plains National Park, which is a high plateau, composed of tropical montane forests. Later dying of forests were observed at several areas including Hakgala montane forest. Physical parameters trace nutrients as well as toxic element concentrations in soils, were studied in order to identify the possible geochemical factors behind the forest die back. Systematic soil sampling was carried out covering the entire Horton Plains National Park and random samples were collected from Hakgala montane forest. Samples were analyzed for available Fe, Mn, Zn, Cu, Ni, Cd, Al and Pb using standard colorimetric and atomic absorption spectrometric procedures. Physical parameters such as pH, moisture content, and conductivity of the samples were also measured. Among extractable micro-nutrients Cu and Zn, and Ni show no deficiency or excess levels. However, the recorded available high concentrations of Fe, Mn and Al can be toxic to certain montane plant species. Acidic moist soil of the area may enhance the toxic effects of these elements. Possible source of these elements should be the underground lithology of the area. According to the results obtained, there is a relationship between forest die back and high Pb concentrations. The same phenomenon was also observed in the Hakgala forest. The distribution pattern of Pb in the Horton Plains coincides well with the die back distribution pattern. The observed Pb values at Horton Plains and Hakgala are almost similar to values observed at Pannipitiya and Dombagaskanda locations, which are located close to main roads carrying heavy traffic. It is quite possible therefore, that Pb toxicity may be a significant factor behind the forest die back even though other factors should not be completely ruled out at this juncture. Strong monsoon winds, bringing Pb from the polluted southwestern part of the country, can be the most possible source of soil Pb. Further studies however are indeed necessary to confirm the conclusion.     

Development planning and gender inequality: A case study in Mahaweli Development Project in Sri Lanka

This study points out that women's status in the development projects is lower than that in the traditional villages in terms of the women's household work load, their access to income generating opportunities, and their participation in development and community activities; and this can be attributed to the development planners' failure to pay due attention to gender differences, because development does not spontaneously benefit men and women equally unless specific women-focussed measures are adopted.     

Heavy metal abundances in the Kandy lake—An environmental case study from Sri Lanka

The Kandy lake, situated in the heart of Sri Lanka's second largest city with a population of nearly 120,000, has been monitored to probe the extent of heavy metal pollution. Although the lake is a source of drinking water to the city, a large number of effluent canals drain into the lake carrying a continuous flow of industrial and domestic waste matter. A total of 66 surface water samples were analyzed for their Fe²⁺, total Fe, total V, SO ₄ ²⁻ , Cd²⁺, and Pb²⁺ contents. Pb and Cd were found in high concentrations averaging 150 μg/l and 77 μg/l, respectively, and exhibit a marked positive correlation with each other (r=+0.94). Vehicular emissions and industrial waste matter contribute largely to the Pb and Cd contents of the lake, the anthropogenic influence outweighing the contributions made by geological materials. All field observations and laboratory experiments indicate a tendency of the Kandy Lake towards eutrophicity.     

Hydrogeochemistry of the groundwater flow system in a crystalline terrain: a study from the Kurunegala district,Sri Lanka

Water samples collected from dug wells and tube wells from the Kurunegala District of Sri Lanka have been studied for their major hydrogeochemical parameters to understand the chemical quality of water in the terrain. The region is composed of Precambrian metamorphic rocks where groundwater is only available in the regolith and along weak structural discontinuities. The study of the major chemical constituents of groundwater revealed several relationships with the aquifer lithology. Groundwater from mafic rocks have high dissolved solids, while quartzose metaclastic rocks yield water with low dissolved solids. The study area displays very low SO₄ ²⁻ contents of the groundwater. The chloride content is higher in the dry regions and in terrains underlain by pink granite and marble/calc gneiss while areas with marble, as expected, show high concentrations of Ca and Mg ions. The waters in the region can be classified into non-dominant cations to Na + K dominant and Cl⁻ and HCO₃ ⁻ dominant types. Water from charnockite-bearing areas tends to have non-dominant cations and more CO₃ ²⁻ + HCO₃ ⁻ types. Effects such as soluble salts in the regolith, fracture intensity and climatic variations play a significant role in the behavior of the hydrogeochemistry in the area.     

Management and environmental constraints to rice yield within a village irrigation system — A case study from Sri Lanka

Management and environmental constraints to rice yield in a paddy tract (lowland rice field) was studied within a village irrigation system in Sri Lanka. The relative importance of water coverage during the critical period and soil fertility variables on the yield variations were analysed. For the upper part of the paddy tract having adequate water supply, environmental and management factors seemed to explain the large variations in yield. For the middle and lower parts, with a shortage of water, water coverage was the only significant explanatory factor.For efficient use of a scarce water supply farmers allocated water first to the upper part, which resulted in water deficiency further down. An alternative strategy for water allocation, where the lower more fertile part of the paddy tract is favoured, would not avoid the risk of water scarcity for the individual farmer but there would be a potential for increase in the total production.     

Numerical simulation and field survey of inundation due to 2004 Indian Ocean tsunami in Trincomalee,Sri Lanka

This paper outlines the field measurements and numerical modelling carried out to develop a high-resolution tsunami inundation map, as a case study, for the city of Trincomalee on the east coast of Sri Lanka, which was devastated by the 2004 tsunami. We employ the deterministic approach together with numerical simulations based on the probable worst-case scenario to derive the inundation map. Linear and non-linear versions of shallow-water equations have been utilized to simulate tsunami propagation and onshore inundation, respectively. The field data considered in the present paper comprise the extent of inundation, the tsunami heights and the arrival times whilst the model results include the spatial distribution of the flow depth, the peak current speeds and the momentum flux. The computed extent of onshore inundation reproduces the observed overall pattern of inundation in most areas barring the south-eastern part of the city. Further, the model simulations suggest maximum flow depths up to about 2 m in most areas of the city whilst patches of flow depths exceeding 2 m can be seen in a narrow strip along the coastline. The computed current speeds also exceed 3 m/s at some locations adjacent to the shoreline.     

Petrogenesis of the Eppawala carbonatites, Sri Lanka: A cathodoluminescence and electron microprobe study

Field and petrographic investigations, cathodoluminescence (CL) studies as well as microprobe analyses of major rock-forming minerals were conducted to establish the crystallization processes in the Eppawala carbonatites, Sri Lanka. The well preserved magmatic textures and crystal morphologies combined with the chemistry of apatite, calcite and dolomite indicate two major stages of crystal growth, which were accompanied by dynamic crystallization conditions. Initially, nucleation of apatite, ilmenite and possibly olivine was associated with rapid crystal growth during slow cooling of the carbonatite melt at depth. The heat loss through the roof and crystallization processes induced the development of turbulent convective currents, which in turn prevented further nucleation and growth of crystals and led to the dispersion of these earlier formed crystals within the magma chamber. Then, rapid upward movement of magma along structural weaknesses led to (i) the transport of mineral clusters, (ii) deformation of ilmenite, (iii) fracturing of apatite and (iv) the emplacement of the carbonatite melt as dykes. Here, the conditions were favourable for the simultaneous crystallization of magnetite, calcite and dolomite in a non-turbulent environment. Subsequent subsolidus alteration caused the hydrothermal overprint of the documented mineral assemblages, particularly along grain boundaries. The study demonstrates that detailed textural examinations of carbonatites combined with mineral chemical analyses and CL investigations can reveal the crystallization processes within carbonatite melts.     

Induced polarization,resistivity, and self-potential: a case history of contamination evaluation due to landfill leakage

This article compares the efficiency of induced polarization (IP) and resistivity in characterizing a contamination plume due to landfill leakage in a typical tropical environment. The resistivity survey revealed denser electrical current flow that induced lower resistivity values due to the high ionic content. The increased ionic concentration diminished the distance of the ionic charges close to the membrane, causing a decrease in the IP phenomena. In addition, the self-potential (SP) method was used to characterize the preferential flow direction of the area. The SP method proved to be effective at determining the flow direction; it is also fast and economical. In this study, the resistivity results were better correlated with the presence of contamination (lower resistivity) than the IP (lower chargeability) data.     

Responses of water quality to land use in riparian buffers: a case study of Huangpu River,China

Along with the continuous urbanization process, massive riparian buffers are developed to satisfy the increasing demands for flood protection and economic growth. Such changes of land use types and landscape patterns pose a significant challenge to the protection of the eco-system of urban rivers. This research concentrates on revealing the potential relation between the water quality of the urban river and land characteristics including land use types and landscape patterns of the riparian area and discusses if the relationship varies with spatial scales. The study area of this research are the 800 m riparian buffer of the six river water quality monitoring sections along Huangpu River. Based on the long-term monthly collection data of six monitoring points in Huangpu River, Shanghai, China, and the land use data of 100 m downstream of 1 km upstream of the monitoring point. Using landscape indicators, GIS spatial analysis, and statistical methods, the responses of land use and landscape patterns to water quality in different scales are obtained. Results show that: (1) Urban Land (UL) and Agricultural Land (AL) have an enormous negative contribution to each buffer, which is the primary non-point source pollution and influencing factor. (2) The active contribution of Wetland (WL) and Forestry Land (FL) is the most significant in the small range buffer of 0–200 m. (3) The negative effect of Industrial Land (IL) on water quality is only significant in the small buffer. (4) The effect of landscape pattern on water quality is more significant on the large spatial scale than the small scale. At the same time, this paper extends the conclusion to discuss the land use of the 800 m buffer zone of Huang River by using point to the area, proving that the research result is helpful for the political decision of urban management and provides suggestions for the land optimization of riparian buffers.

     

Anthropogenic impacts recorded in the sediments of Lunawa, a small tropical estuary, Sri Lanka

The Lunawa Estuary has been selected as a model for studying the long-term diagnoss of tropical estuarine pollution and its impact on ecosystems. Surface and short-core sediments from the Lunawa Estuary were collected, and analyzed for total organic carbon (TOC), nitrogen, sulfur and hydrocarbon contents, hydrogen index and other related sediment properties. Water quality and surface sediment data indicate that the Lunawa Estuary is presently an anoxic water body. The upper part of the sediment cores, characterized by very high TOC values of more than 10%, and absence of bioturbation along with strong laminations indicate anoxic water since 1970 AD, determined by ¹⁴ C dating. However, before 1970 AD, lower TOC content and presence of bioturbation and homogenous sediments with animal traces suggest the presence of oxic bottom water and thus a much better physicochemical condition. The lower part of the cores may also suggest the influence of gradual global warming. Regular laminations in the upper portions of the cores probably resulted from monsoon pulses with material preserved under anoxic condition.     

Factors affecting the productivity of manufacturing sector in Sri Lanka: a spatial analysis

Two indexes of productivity — labour productivity and capital productivity — have been employed to measure and map productivity patterns of manufacturing industries in Sri Lanka. There exist extensive regional differences in the levels of productivity. Regression results show that the spatial variation of labour productivity is positively related to capital intensity and the level of urbanization and is negatively related to the nature of enterprises. Capital intensity and level of urbanisation have a negative influence on the capital productivity. The selected five variables in the regression analyses explain 65% of the spatial variation of labour productivity and 39% of capital productivity.     

Hydrological responses to rainfall variability and dam construction: a case study of the upper Senegal River basin

The understanding of the spatial and temporal dynamic of river systems is essential for developing sustainable water resource management plan. For the Senegal River, this subject is very complex according to the context of (1) transboundary basin, (2) several contrasted climatic zones (Guinea, South Sudanian, North Sudanian and Sahelian) with high rainfall variability and (3) high human pressures (dam construction and water uses). From 1954 to 2000, 80% (mean value) of the Senegal River flows recorded downstream part of the basin are provided by three majors tributaries (Bafing, Bakoye and Faléme) located in the upstream part. Then, in our study, this upper Senegal River basin was chosen in order to investigate the hydrological responses to rainfall variability and dam construction. Two nonparametric statistical methods, Mann–Kendall and Hubert test, were used to detect the long-term changes in the time series of precipitation and water discharge (1954–2000) at the annual and seasonal scales. The continuous wavelet transform (Morlet Wavelet) was employed to characterize the different mode in the water discharge variability. Flow duration curve and cumulative curve methods were used to assess the impact of dams on the hydrological regime of the Senegal River. Results showed that the Senegal River flows have been changing under the influence of both rainfall variation and dam construction. The long-term evolution of water discharge depend on long-term rainfall variability: The wet periods of the 1950s and 1960s correspond to periods of higher river flows, while the droughts of the 1970s and 1980s led to unprecedented river flows deficits. The new period, since 1994, show a high inter-annual variability of rainfall and discharge without clear trend. At seasonal scale, the results showed also a strong relationship between rainfall and runoff (R ² > 0.8) resulting from alternating wet and dry seasons and rapid hydrological responses according to annual rainfall. Nevertheless, the observed flows during dry seasons highlighted the influence of water storage and restitution of infiltrated waters in soils and surficial formations during wet seasons. In the dry seasons, the water budget of the three upstream tributaries showed a water deficit at the downstream gauging station. This deficit was characterized by water loss to underlying aquifers and highlighted the influence of geological setting on water balance. However, in this context, water restitution during the dry season remained dependent on climatic zone and on the total annual rainfall volume during the previous wet season. The results have highlighted an impact of the Manantali dam previously obscured: The dam has no effect on the regulation of high river flows. That is what explains that since its construction in 1988, flooding of coastal cities, like Saint-louis, by seasonal river floods has not ceased. The flooding risk in coastal cities is not avoided, and the dams caused hyper-salinization of the Senegal lower estuary. The breach created in the coastal barrier of the Langue of Barbary in October 2003 promotes direct export of excess floodwater to the sea and reduces this risk of flooding in the delta area. But, this solution led to considerable loss of potential water resources, and the authors recommend a new water management plan with a global focus. However, this study shows the positives impacts of the two dams. They allow the availability of freshwater in order to support agricultural irrigation in the valley and delta zone, in particular during low flows periods.     

Study of geothermal water intrusion due to groundwater exploitation in the Puebla Valley aquifer system, Mexico

Significant intrusion of geothermal water into fresh groundwater takes place in the Puebla Valley aquifer system, Mexico. The decline in the potentiometric surface due to the overexploitation of the groundwater induces this intrusion. This hydrological system comprises three aquifers located in Plio-Quaternary volcanic sediments and Mesozoic calcareous rocks. The hydraulic balance of the aquifer shows that the annual output exceeds the natural inputs by 12 million m³. Between 1973 and 2002, a drop in the potentiometric surface, with an 80 m cone of depression, was identified in a 5-km-wide area located southwest of the city of Puebla. Chemical analyses performed on water samples since 1990 have shown an increase in total dissolved solids (TDS) of more than 500 mg/L, coinciding with the region showing a cone of depression in the potentiometric surface. A three-dimensional flow and transport model, based on the hydrogeological and geophysical studies, was computed by using the MODFLOW and MT3D software. This model reproduces the evolution of the aquifer system during the last 30 years and predicts for 2010 an additional drawdown in the potentiometric surface of 15 m, and an increase in the geothermal water intrusion.     

An integrated approach to investigate saline water intrusion and to identify the salinity sources in the Central Godavari delta, Andhra Pradesh, India

The Central Godavari delta is located along the Bay of Bengal Coast, Andhra Pradesh, India, and is drained by Pikaleru, Kunavaram and Vasalatippa drains. There is no groundwater pumping for agriculture as wells as for domestic purpose due to the brackish nature of the groundwater at shallow depths. The groundwater table depths vary from 0.8 to 3.4 m and in the Ravva Onshore wells, 4.5 to 13.3 m. Electrical Resistivity Tomography (ERT) surveys were carried out at several locations in the delta to delineate the aquifer geometry and to identify saline water aquifer zones. Groundwater samples collected and analyzed for major ions for assessing the saline water intrusion and to identify the salinity origin in the delta region. The results derived from ERT indicated low resistivity values in the area, which can be attributed to the existence of thick marine clays from ground surface to 12–15 m below ground level near the coast and high resistivity values are due to the presence of coarse sand with freshwater away from the coast. The resistivity values similar to saline water <0.01 Ω m is attributed to the mixing of the saline water along surface water drains. In the Ravva Onshore Terminal low resistivity values indicated up coning of saline water and mixing of saline water from Pikaleru drain. The SO ₄ ^?2 /Cl^?and Na⁺²/Cl^?ratios did not indicate saline water intrusion and the salinity is due to marine palaeosalinity, dilution of marine clays and dissolution of evaporites.