Next Article in Journal
Effects of Sunflower Meal Supplementation as a Complementary Protein Source in the Laying Hen’s Diet on Productive Performance, Egg Quality, and Nutrient Digestibility
Next Article in Special Issue
Socio-Economic and Environmental Impacts of Biomass Valorisation: A Strategic Drive for Sustainable Bioeconomy
Previous Article in Journal
Model of the Consumer Switching Behavior Related to Healthy Food Products
Previous Article in Special Issue
Current Status and Potential of Tire Pyrolysis Oil Production as an Alternative Fuel in Developing Countries
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Multibranch Model for Metals in River Systems: Impacts and Control of Tannery Wastes in Bangladesh

1
School of Geography and the Environment, University of Oxford, South Parks Road, Oxford OX1 3QY, UK
2
Institute of Water and Flood Management (IWFM), Bangladesh University of Engineering and Technology (BUET), Dhaka 1000, Bangladesh
3
Department of Earth Sciences, University of Oxford, South Parks Road, Oxford OX1 3AN, UK
4
Oxford Molecular Biosensors, Centre for Innovation and Enterprise, Begbroke, Oxford OX5 1PF, UK
5
Geology Department, State University of New York College at Cortland, Cortland, NY 13045, USA
*
Author to whom correspondence should be addressed.
Sustainability 2021, 13(6), 3556; https://0-doi-org.brum.beds.ac.uk/10.3390/su13063556
Submission received: 19 February 2021 / Revised: 15 March 2021 / Accepted: 17 March 2021 / Published: 23 March 2021

Abstract

:
A new multibranch Integrated Catchment (INCA) model INCA-Metals has been developed to simulate the impact of tannery discharges on river systems. The model accounts for the key chemical reaction kinetic processes operating as well as sedimentation, resuspension, dilution, mixing and redistribution of pollutants in rivers downstream of tannery discharge points and for mine discharges or acid rock drainage sites. The model is dynamic and simulates the daily behaviour of hydrology and eight metals, including cadmium, mercury, copper, zinc, lead, arsenic, manganese and chromium, as well as cyanide and ammonia. The model is semi-distributed and can simulate catchments, tributaries and instream river behaviour. The model can also account for diffuse pollution from rural runoff as well as point sources from effluent and trade discharges. The model has been applied to the new Savar tannery complex on the Dhaleshwari River system in Bangladesh to assess the impacts on pollution levels in the river system and to evaluate a set of treatment scenarios for pollution control, particularly in the dry season. It is shown that the new effluent treatment plant at Savar needs to significantly improve its operation and treatment capability in order to alleviate metal pollution in the downstream Dhaleshwari River System and also protect the Meghna River System that falls in the Bay of Bengal.

1. Introduction

Globally, polluted surface and ground water causes harm to ecosystems and human health. Sustainable Development Goal (SDG) 6 [1] for water and sanitation builds on earlier frameworks established by the Millennium Development Goals to address the complexity of interconnected water quality, quantity and affordability concerns. SDG indicator 6.3 seeks to achieve improvement in ambient water quality by significantly reducing flows of pollution into water bodies. In many developing countries where progress toward SDG 6 is off-target, combinations of untreated industrial, domestic, urban and agricultural waste contribute to dangerous levels of toxicity. The reduction of heavy metal pollution from leather tanneries, the garment industry and mining activities is a key issue and is of major concern for the protection of river water quality and ecology, as well as for sustaining livelihoods for people living close to river systems.
Metal pollution from tanneries is a major issue in countries such as Bangladesh and India [2,3], where metals such as chromium, copper and cadmium are released from the tannery processes, often with minimal effluent treatment, causing severe pollution of local streams and rivers. Other metal pollution arises from different industries such as mine effluents and acid rock drainage (ARD) which can have significant impacts on river systems [4]. Treatment methods for both tanneries and the mining industry are often problematic, as it is expensive, technically demanding and insufficiently monitored by regulatory authorities. Many active and passive waste treatment systems have been investigated and applied to control tannery effluents and mine discharges [5,6]. The problem of pollution following mine closure is particularly acute in many European countries where historical mining has been undertaken with poor pollution control and where the high costs of treatment make restoration a major issue [7].
This paper focuses on pollution in the Dhaleshwari River system in Dhaka, Bangladesh, where a new complex of tanneries at Savar was established in 2016 (Figure 1). The Dhaleshwari river system is part of the Jamuna-Brahmaputra catchment and is the major source of freshwater to the Greater Dhaka Watershed. Flows in the river have decreased due to natural as well as human interventions, but augmentation of flow during the dry season in the Dhaleshwari and some other rivers is under consideration targeting pollution reduction and ecological restoration works planned under the Bangladesh River Masterplan and other initiatives.
A new novel multi-branch model, INCA-Metals, has been developed to model pollution impacts in rivers and their tributaries, and the model has been applied to assess metal pollution and a range of restoration strategies. The INCA-Metals model builds on the existing suite of INCA models that have been developed over the past 20 years to simulate both terrestrial and aquatic systems [8,9,10,11]. The structure of INCA has been tested on many catchments, including 20 UK catchments, 25 catchments across Europe, and globally on large river systems such as the Mekong, the Ganges and the Volta [12,13,14]. Early major applications of INCA have been published in two special volumes of international journals [15,16]. The advantage of INCA is that the model seeks to incorporate the dominant mechanisms and processes operating so that a realistic and rapid assessment of environmental change can be evaluated. In this paper, a new version of the model, INCA-Metals, is described for rivers and catchments to simulate eight metals (cadmium, lead, zinc, mercury, arsenic, copper, chromium, manganese), as well as cyanide and ammonia. Cyanide and ammonia are included so that the metals model can be applied to mining waste issues. Cyanide is a key component of gold metal extraction processes, and ammonia is often released from the mining effluent treatment process. Both can be extremely toxic to river systems and hence have been included in the model.

2. The INCA-Metals Model

The Integrated Catchment Framework Model (INCA) is used as the main platform for the new multibranch metals model. INCA is a catchment scale process-based model to calculate pollutant transfer from diffuse sources and point sources to the catchment outlet. To date, the INCA suite of models includes simulation of flow, nitrogen, phosphorus, sediments, chloride, carbon, pathogens, dissolved oxygen and mercury [8,9,10,11,17,18,19]. Many of the parameters and processes that have previously been successfully validated in other versions of INCA also apply to modelling of metals, such as the hydrology, landscape pathways, sediment settlement and entrainment, and flushing out components. Due to the dynamic nature of the model, variations in rainfall and temperature and changes of inputs, such as tannery treatment or mine discharges and runoff from diffuse pollution, can be investigated. This means the effects of climate change, population growth, land-use change and mitigation measures can be evaluated.
INCA-Metals is a mass-balance dynamic model that estimates the daily fluxes and concentrations of flow, ammonia, cyanide and eight metals in rivers. The eight metals are cadmium, lead, zinc, mercury, arsenic, copper, chromium and manganese. The flows and metals are simulated based on mass balance calculations accounting for contributions from various inputs and transformations. Key physico-chemical processes and stores are assumed in the landscape and within the river, and these are shown in Figure 1. The equations for INCA-Metals are based on those written for the Integrated Catchment model of Nitrogen (INCA-N [8,9,10]) but have been adapted to describe the metal adsorption to sediment, cyanide decay to ammonium and cyanide volatilization. The inputs and outputs are differentiated by landscape type and varied according to the environmental conditions, namely soil moisture and temperature. The model accounts for stocks of ammonium, cyanide and the eight metals in the soil and groundwater zones and in the stream reaches. The model also simulates the flow of water through the soil and groundwater zones from different soil and land use types to deliver mass to the river. This mass is then routed downstream after accounting for direct inputs from point sources, abstractions, within-river precipitation and sedimentation of metals, nitrification and cyanide decay and volatilization. Point sources such as discharges from tannery effluents or mine adits, waste dumps or tailings ponds, or effluent treatment facilities can be routed into any of the sub-catchments or reaches of the model.

2.1. The Simulation of Water Flow and Storage in the Landscape

The landscape mass balances of water, cyanide, ammonium and the eight metals are based on a 1 km2 cell (Figure 2). The inputs to the model can vary on a sub-catchment basis and according to soil or land-use type. These two factors allow the mass stored, process rates and hydrological pathways to vary spatially based on variations in soil moisture, temperature, adsorption potential and land management. The water volumes and the mass of cyanide, ammonium and the eight metals are summed based on the relative amounts of each land use or soil type within a sub-catchment and the output passed to the instream routing model. The hydrology component of the model is represented in Figure 3, where precipitation flows through the soil and groundwater zones, although there can be overland flow or direct runoff flow in high rainfall events.
The flow of water through the soil and groundwater zones is given by the following two equations:
Soil Zone
d q s z d t = p e f f q s z T s z
Groundwater Zone
d q g z d t = β q s z q z T g z
where qsz and qgz are the outflows from the soil and groundwater zones (m3 s−1 km−2), peff is the hydrologically effective rainfall (m3 s−1 km−2), β is the base flow index (Ø) and Tsz and Tgz are the response times associated soil and groundwater zones (days).
Within the soil zone, it is assumed the water can be partitioned into two volumes: drainage and retention. The drainage volume represents the water stored in the soil that responds rapidly to water inflow and drains under gravity; it may be thought of as macropore or drain flow: the flow that most strongly influences the rising hydrograph limb. The soil zone retention volume represents the water stored retained in the soil after gravity drainage; it responds more slowly than the drainage water and represents the majority of water in the soil.
The initial value of the soil zone drainage volume (VD, m3 km−2) is calculated from a user-supplied initial soil zone flow (qsz,initial, m3 km−2) and the soil water response time:
V D , i n i t i a l = q s z , i n i t i a l · T s z · 86 , 400
where VD,initial is the initial soil water drainage volume (m3 km−2), qsz,initial is the user-supplied initial soil water flow rate (m3 s−1) and Tsz is the user-supplied soil water time constant (days).
The initial value of the soil zone retention volume (VR, m3 km−2) is calculated from the supplied soil moisture deficit time-series (SMD, mm), an estimate of the maximum soil moisture deficit (SMDmax, mm) and a parameter that describes the linear relationship between the soil moisture deficit and the soil zone retention volume, C1. The value of this parameter represents the ratio of the total water in the retention volume relative to the easily available water, namely that about the wilting point. The value of the parameter is derived through calibration and typically takes a value of 1 to 3. The value of SMDinitial is estimated from the soil moisture deficit on day one of the simulation.
V R , i n i t i a l = C 1 ( S M D m a x S M D i n i t i a l ) · 1000
The groundwater initial volume (Vgw, m3 km−2) is estimated from the maximum size of the store and the proportion of pore space filled with water at the start of the model run:
V g w , i n i t i a l = d e f f , g w · C 2 · 10 6
where Vgw,initial is the initial groundwater volume (m3 km−2), deff,gw is the user-supplied maximum groundwater effective depth (m, active depth x effective porosity) and C2 is the user-supplied proportion of filled pore-space (Ø). Given the complexity of the geology in most model applications, there is no attempt to separate the groundwater into specific yield and retention components; the former is the water able to drain from the rock under gravity, and the latter is the water retained against the pull of gravity. The soil zone drainage and retention volumes and the groundwater zone volume are recalculated at each time step to account for the changes in input and output.

2.2. The Simulation of the Transport, Storage and Transformations of Cyanide, Ammonium and Metals in the Landscape

The change in cyanide mass in the soil, mcn,sz (kg CN km−2) and groundwater, mcn,gz (kg CN km−2) stores are given by Equations (6) and (7):
Soil Zone
d m c n , s z d t = m c n , s z q s z 86,400 V D + V R C 3 S S M D m c n , s z 10 6 V D + V R C 4 S S M D m c n , s z 10 6 V D + V R
Groundwater Zone
d m c n , g z d t = β m c n , g z q s z 86,400 V g w m c n , g z q g z 86,400 V g w C 5 m c n , g z 10 6 V g w
where C3, C4 and C5 are the rates of cyanide volatilisation and soil and groundwater zone conversion of cyanide to ammonium (Figure 4). All other terms have been defined previously except SSMD, which is the soil moisture factor and describes the linear dependency of the rate of the soil zone processes on the soil moisture. It is assumed that there are no diffuse inputs of cyanide to the catchment. The first term on the right-hand side of Equation (6) represents the lateral transport of cyanide with the soil water to the stream, the second term represents cyanide volatilisation and the third term represents the decay of cyanide to ammonium. The first and second terms in Equation (7) represent the flux of cyanide into the groundwater from the soil zone and the lateral flow of cyanide with the groundwater into the stream; the third term represents the decay of cyanide to ammonium in the groundwater. The soil moisture factor is calculated at each time step as
S S M D = S M D M A X S M D S M D M A X
where SMD is the user-supplied daily soil moisture deficit time series (mm). The factor is scaled between 0 and 1 and describes the situation in which, as the soil dries, the rate of the soil processes declines. In addition, each process rate parameter is a function of the soil temperature
C n = C n t Q 10 θ s t Q 10 b a z 10
where θs is the soil temperature (°C), Cn is the soil process parameter and tq10 (Ø) and tQ10bas (°C) are parameters determined by calibration. The parameter tq10 is the factor change in rate with a 10-degree change in temperature, and the parameter tQ10bas is the base temperature for the process at which the response is 1. The soil temperature is estimated from a seasonal relationship dependent on air temperature as follows
θ S = θ A C 6 sin ( 3 2 π d a y   n o . 365 )
where θs is the air temperature (°C) and C6 is the maximum temperature difference between summer and winter (°C). This relationship generates a seasonal pattern for each land use, which is controlled by the parameter C6. It is corrected for snow depth during winter months using Equation (4) from Rankinen et al. (2004) [20]. This temperature dependency is applied to all soil zone process parameters, not just those for cyanide.

2.3. Ammonium-N

The change in ammonium mass in the soil, mnh4,sz (kg N km−2) and groundwater, and mnh4,gz (kg N km−2) stores are given by Equations (11) and (12):
Soil Zone
d m n h 4 , s z d t = 100 m n h 4 , s z m n h 4 , s z q s z 86,400 V D + V R C 7 S S M D S P G I m n h 4 , s z 10 6 V D + V R C 8 S S M D m n h 4 , s z 10 6 V D + V R + C 4 S S M D m c n , s z 10 6 V D + V R
Groundwater Zone
d m n h 4 , g z d t = β m n h 4 , g z q s z 86,400 V g w m n h 4 , g z q g z 86,400 V g w C 9 m n h 4 , g z 10 6 V g w
where C7, C8 and C9 are the rates of ammonium uptake by plants from the soil and the soil and groundwater ammonium nitrification rates (Figure 4). The input mass of ammonium to the soil zone, mnh4,in (kg N ha day−1), represented by the first term on the right-hand side of Equation (11), includes livestock manure, fertiliser and wet and dry deposition. All other terms have been defined previously except SPGI, which is the plant growth index and describes the seasonal variations in the solar radiation and the subsequent control on plant growth during different seasons. The index is given by
S P G I = 0.66 + 0.34 sin ( 2 π [ d a y   n o . C 10 ] 365 )
where C10 is the day number associated with the start of the growing season. The second term on the right-hand side of Equation (11) represents the lateral movement of ammonium with flow from the soil zone to the stream, the third term represents the ammonium uptake by plants, the fourth term represents the conversion of ammonium to nitrate via nitrification and the fourth term is the increase in ammonium resultant from the decay of cyanide. The first and second terms on the right-hand side of Equation (12) represent the input of ammonium from the soil zone to groundwater and the lateral movement to the stream; the third term represents ammonium nitrification in the groundwater.

2.4. Metals

The following equations relate to each of the eight metals included in the model. Whilst the equations given here are generic to all the eight metals, within the model there are different parameters for each metal so that the transport, decay and storage of each are independent.
The change in mass of each metal in the soil, mmetal,sz (kg km−2) and groundwater, mmetal,gz (kg km−2) stores are given by Equations (14) and (15)
Soil Zone
d m m e t a l , s z d t = m m e t a l , s z q s z 86,400 V D + V R C 11 S S M D m m e t a l , s z 10 6 V D + V R
Groundwater Zone
d m m e t a l , g z d t = β m m e t a l , g z q s z 86,400 V g w m m e t a l , g z q g z 86,400 V g w C 12 m m e t a l , g z 10 6 V g w
where C11 and C12 are the rates of metal adsorption to the soil sediment and aquifer matrix. It is assumed that there are no diffuse inputs of metals to the catchment. The first term on the right-hand side of Equation (14) represents the lateral movement of metals transported with the flow to the stream. The second term represents the adsorption of a metal to the soil sediment. The first term on the right-hand side of Equation (15) represents the input of a metal from the soil to the groundwater zone by percolation, the second term represents the lateral flow of the metal from the groundwater to the stream and the third term represents the adsorption of metal to the aquifer matrix. The main flow paths for the metals are shown in Figure 5.

2.5. The Simulation of Water Flow and Storage in the River

The reach residence time constant, Treach (days) is calculated as
T r e a c h = L a q r e a c h , o u t b 86,400
where L is the reach length, qreach,out is the discharge from the reach and a and b are parameters relating the reach velocity to the discharge. The parameters a and b are determined through calibration, though b typically has a value of 0.67. The parameters can also be determined from tracer measurements. The change in the reach flow is calculated from the input–output mass balance of the form
d q r e a c h , o u t d t = q r e a c h , i n q r e a c h , o u t T r e a c h
where qreach,in is the sum of the input flows from upstream, if not the headwater reach, point source effluent, diffuse inputs from the soil and groundwater zones and a loss via abstraction.

2.6. Cyanide, Ammonium, and Metal Process Equation: River System

In the river, the key processes are cyanide volatilisation and degradation, ammonium nitrification and metal loss due to sedimentation or precipitation. The reach mass balance includes the upstream water quality together with diffuse inputs from the soil and groundwater zones, as well as direct effluent discharges and abstractions (Figure 6).
The cyanide mass, mcn,reach (kg) stored in a river reach is given by
d m c n , r e a c h d t = m c n , r e a c h _ i n m c n , r e a c h q r e a c h , o u t 86,400 V r e a c h C 13 c c n , r e a c h V r e a c h 1000 C 14 c c n , r e a c h V r e a c h 1000 m c n , a b s
where the mass into the reach, mcn,reach_in (kg day−1) is the sum of the upstream input, point source effluent and the diffuse inputs from the soil and groundwater. The second term on the right-hand side of Equation (18) represents the mass transfer downstream with the flow of water, the third term represents the volatilisation of cyanide, the fourth term represents the decay of cyanide to ammonium and the fifth term represents any mass that is removed from the reach by abstraction. The rate parameters describing cyanide volatilisation, C13, and decay, C14, are both temperatures dependent such that
C n = C n 1.047 θ W 20
where θw is the water temperature (°C), which is assumed equal to the input air temperature.
The ammonium mass, mcn,reach (kg N) stored in a river reach is given by
d m n h 4 , r e a c h d t = m n h 4 , r e a c h _ i n m n h 4 , r e a c h q r e a c h , o u t 86,400 V r e a c h C 15 c n h 4 , r e a c h V r e a c h 1000 C 14 c c n , r e a c h V r e a c h 1000 m n h 4 , a b s
where the mass into the reach, mnh4,reach_in (kg day−1) is the sum of the upstream input, point source effluent and the diffuse inputs from the soil and groundwater. The second term on the right-hand side of Equation (20) represents the mass transfer downstream with the flow of water, the third term represents the nitrification of ammonium, the fourth term represents the increase in ammonium from the decay of cyanide and the fifth term represents any mass that is removed from the reach by abstraction. The rate parameters describing nitrification C15 is temperature dependent.
The metal mass, mmetal,reach (kg N) stored in a river reach is given by
d m m e t a l , r e a c h d t = m m e t a l , r e a c h _ i n m m e t a l , r e a c h q r e a c h , o u t 86,400 V r e a c h C 16 c m e t a l , r e a c h V r e a c h 1000 m m e t a l , a b s
where the mass into the reach, mmetal,reach_in (kg day−1) is the sum of the upstream input, point source effluent and the diffuse inputs from the soil and groundwater. The second term on the right-hand side of Equation (21) represents the mass transfer downstream with the flow of water; the third term represents the sedimentation of the metal.; and the fourth term represents any mass that is removed from the reach by abstraction. The rate parameters describing sedimentation C16 is temperature-dependent.
Within the soil zone, it is assumed that the water can be partitioned into two volumes: drainage and retention. The drainage volume represents the water stored in the soil that responds rapidly to water inflow and drains under gravity; it may be thought of as macropore or drain flow (i.e., the flow that most strongly influences the rising hydrograph limb). The soil zone retention volume represents the water stored or retained in the soil after gravity drainage; it responds more slowly than the drainage water and represents the majority of water in the soil.
All the model constants and variables in the landscape equations and in the river equations are given in Table 1 and Table 2.

3. Pollution Impacts on the Dhaleshwari River (Bangladesh)

3.1. Tannery Pollution Globally and in Bangladesh

With about 300 million hides treated globally per year, leather tanning is one of the most environmentally problematic industries due to its high water consumption and toxic wastewater. Tannery effluents generally have a complex composition with very high levels of organic and inorganic pollutants as well as high concentrations of heavy metals, such as chromium, lead and cadmium. These can be extremely damaging to human health and aquatic ecosystems. Chromium tanning is the most common process due to its relatively short turnaround time and excellent quality yield. Around 90% of the world’s leather gets Cr-tanned [21] using over 15,000 tons of Cr-salts annually. Only 20% of the chemical used in the production process is fixed, and therefore large fluxes of chromium are released via tannery wastewater into river systems. In Bangladesh, the tannery industry has recently been relocated from Hazaribag to Savar in Dhaka (Figure 7), with a total of 155 factories moved. This move was designed to relieve serious pollution in the Buriganga River and provide a central effluent treatment plant (CETP) at Savar. However, ongoing delays in effective operation of the CETP have called into question whether the pollution problem has been just transferred from the Buriganga River to the Dhaleshwari River (Figure 7).

3.2. Impacts of Tannery Discharges on Pollution on Bangladesh Rivers

Tanneries generate significant polluting discharges into river systems unless effective effluent treatment is applied. Table 3 shows the mean chemistry from 16 published studies of tannery waste in Bangladesh. There are high levels of pollutants in the effluents, and Table 4 shows the average of these, with very low dissolved oxygen levels, corresponding with high BOD (Biological Oxygen Demand) and COD (Chemical Oxygen Demand) concentrations as well as dissolved solids discharges. Nutrients such as nitrate and phosphorus are also significant, as is chloride (or salt content). Most concerning are the metal concentrations, with extremely high lead, zinc, cadmium and chromium concentrations.
The high concentrations of metals, nutrients and organic compounds found in the tannery effluents have been shown to be deleterious to the aquatic eco-system. There are a wealth of studies on the harmful impacts of the tannery discharges on the Buriganga River water quality, sediments and fish populations [22]. These hazards can also harm human health through bio-magnification with metals such as chromium causing acute renal failure, gastrointestinal haemorrhage, lung cancer, asthma and nervous system damage, as shown in Table 5. Other metals such as lead, for example, can impair children’s development, cause memory loss and insomnia, cause vomiting, diarrhoea and kidney damage (Table 5).
Table 3. Characteristics of Tannery Wastewater in Bangladesh.
Table 3. Characteristics of Tannery Wastewater in Bangladesh.
References
[23][24][25][26][27][28][29][30][31][32][33][34][35][36][37][38]
Dissolved oxygen (mg/L) 0.5 2.73 BDL 2.71 3.085 1.751.3
Total dissolved solids (mg/L) 2711454032004119 719077,92321,300 7.54745213091513,00044.3
Total suspended solids (mg/L)12,3095010 850710 18159849125068880.86 23575500
Biological oxygen demand (mg/L)9291.110,400 92.1 <0.11364 4464920 90355048303621
Chemical oxygen demand (mg/L)15,30515,800 3380 12,8403980 220011,23837905213
pH7.785.47.2 6.327.892.78.37.57.57.0457.48.57.33.6
Phosphate (mg/L) 1.79 50.9 0.55 17.1
Sulphate (mg/L)1410.54648 960 4000 3850
Sulphide (mg/L)220.35040
Nitrate (mg/L)226.1120 0.7 0.14
Nitrite (mg/L) 0.36 1.17 <0.01
Chloride (mg/L)31,533960 83.6 2200 13.85000 104 19,015
Colour (Pt-Co)13,567792 275 1640
Sodium (mg/L) 1835 3400836 12,0063840
Lead (mg/L)0.680.1134.8 7.250.091 0.18 2.924.915 2.153.4
Copper (mg/L) 0.22BDL 0.058 0.41 102.05
Iron (mg/L) 8.49 0.27 7.6 14.7 54.968.4412.24
Chromium (mg/L)1445.118.9846,848 24.999943403810.35 79103391.92357.825.63190.1
Cadmium (mg/L) 0.0150.6 BDL<0.001 0.005 2.91.875<11.6
Zinc (mg/L)1.396 608.1 0.370.27 1.52 4.68 6.3
Nickel (mg/L) BDL BDL0.046 0.15 3.6BDL 3.11.35
Manganese (mg/L) 0.08 0.290.49 8.845
Arsenic (mg/L) <0.01 16
Table note: BDL = Below Detection Limit.
Table 4. The mean chemistry of Tannery Wastewater in Bangladesh.
Table 4. The mean chemistry of Tannery Wastewater in Bangladesh.
Dissolved oxygen (mg/L)1.4
Total dissolved solids (mg/L)17,283.3
Total suspended solids (mg/L)4835.1
Biological oxygen demand (mg/L)3790.2
Chemical oxygen demand (mg/L)10,261.0
pH6.9
Phosphate (mg/L)17.6
Sulphate (mg/L)2754.6
Sulphide (mg/L)2630.2
Nitrate (mg/L)61.7
Nitrite (mg/L)0.5
Chloride (mg/L)6631.7
Color (Pt-Co)7179.5
Sodium (mg/L)4383.4
Lead (mg/L)7.2
Copper (mg/L)0.2
Iron (mg/L)7.8
Chromium (mg/L)6678.4
Cadmium (mg/L)0.2
Zinc (mg/L)122.3
Nickel (mg/L)0.0
Manganese (mg/L)0.3

4. An Experiment on the Dhaleshwari River System and the Savar Tannery Complex

The pollution in Hazaribag has been of national concern since the 1980s, when the government first instructed tanneries to control their pollution footprint [40,41]. In 2001, The High Court of Bangladesh ordered the tanneries to apply the necessary measures to limit their toxic discharges within one year, following the submission of a written petition by the Bangladesh Environmental Lawyers Association (BELA). In 2002, the government approved a relocation plan to move the tanneries to Savar district (adjacent to the Dhaleshwari river), and in 2003 the Bangladesh Government signed an industry transfer agreement with the Small and Cottage Industries Corporation (BSCIC) and major leather tannery industry associations. Due to a series of political and economic issues, the relocation only formally took place in 2017. While currently over 90% of the tanneries have now moved to Savar, the new effluent treatment plant is still incomplete and not operating fully and efficiently. The introduction of tannery discharges in Savar has created additional pollution in the Dhaleshwari River system. Studies have revealed a gradual deterioration of its water quality over the past two decades [42,43,44,45], with several urban and industrial discharges along the river. Heavy metals have risen in water, sediments, and with alarming levels recorded in fish, which have been designated “unfit for human consumption” [43].
Few studies have been carried out to monitor and model pollution scenarios in the Dhaleshwari. Two studies assessed the relocation of the tanneries to Savar, but only focused on the impact on the Buriganga. Islam (2018) [40] used the Water Analysis Simulation Program (WASP) model to investigate the impact of the tanneries’ relocation on the physico-chemical parameters of the Buriganga. BOD, ammonia, nitrate and phosphate levels decreased significantly, but DO levels did not improve. Whitehead et al. [41] studied heavy metals (Cd, Cr, As, Cu) and treatment scenarios in the Buriganga River using the INCA-metals model. The simulations revealed that unless an 85% wastewater treatment level was attained, the heavy metals would still be above permissible limits in the river system.
In order to assess the pollution impacts of the new Savar tannery complex, a sampling and chemical analysis programme was conducted in August and September 2019. Figure 8 shows the sampling sites that are located upstream and downstream of the Savar tannery discharge point. The sampling strategy involved sampling at a depth of 0.6 m from water surface, using a Teflon-coated plastic bottle (autoclavable and acid-washed) with samples transfer in iceboxes during fieldwork, then in a refrigerator (at 4 °C) before their shipment to Oxford, UK. Chemical analysis was conducted using Inductively Coupled Plasma-Mass Spectrometry (ICP-MS). Measurements by ICP-MS were undertaken using a PerkinElmer NexION 350D ICP-MS at the Department of Earth Sciences, University of Oxford. Each metal was calibrated from a series of calibration standards, which were automatically produced by the coupled Elemental Scientific prepFAST M5 syringe-driven autosampler.
Samples were taken weekly for six weeks from all the sites (Figure 8), and chemical analysis conducted. Table 6 shows the average concentrations of metals in the upstream and downstream sites and the percent change between the two sites. The data show a consistent increase in pollution from upstream to downstream with chromium and cadmium, showing 254% and 65% increases, respectively. The metals concentrations at the upstream site are significant in themselves and this represents pollution from upstream discharges. So the new Savar tannary discharges are just exacerbating the existing pollution in the river system. The concentrations in the river are quite low because the river flows were high during the sampling programme and hence river dilution was significant. It should be noted that the sampling period coincided with the peak production season (post-Eid Al Adha), which follows a lunar calendar. As EID occurs at different times of the year it is important to consider the pollution levels in different seasons and especially during low flow periods when dilution is low. Therefore, in order to assess the impact of the river under low flow conditions, the river’s water quality was further investigated under low-flow conditions. Based on the historical flow dataset of the Dhaleshwari River system, Mimouni [39] has estimated the dilution factor to be of the order of 130. This implies that metal concentrations in the river would be approximately 130 times greater in the dry season compared to the wet season (Table 6). Of course, other environmental variables change on a seasonal basis such as temperature, and it is important to bear these in mind during any modelling study. As described above, the INCA equations are temperature-dependant and so will change with changing temperature. The model also considers the changing water velocities and hence residence times so that decay and sedimentation rates will vary from high-flow to low-flow conditions.
Since other polluting sources also release their effluents into the Dhaleshwari, attributing the pollution loads to the tanneries’ discharges remains challenging. In order to identify the heavy metals that were particularly impacted by the tanneries’ effluent during the sampling period, an amalgamation of different factors/conditions (spatial and temporal) that corroborate the potential contribution of the CETP effluent to the heavy metals concentrations recorded were taken into consideration. Any impact attributed to the CETP effluent was also based on their effect on the Dhaleshwari’s water quality vis-à-vis international/regional surface water quality standards (Table 7) and the link of each analysed heavy metal to the tanning process. Strong impacts of three heavy metals were recorded: Cr, Cd and As. Considering the environmental threats that high levels of these heavy metals entail, the improvement of CETP’s treatment performance is imperative to ensure the compliance of its Cr, Cd and As discharges. This must be prioritized during all seasons not to exhaust the river’s assimilative capacity and to minimize the toxicity threats to the wider ecosystem and local population.

5. Modelling the Savar Discharge and Mitigation Strategies

In order to determine the percentage level of upgrade needed in the CETP and highlight the impact of its improved performance on the water quality of the Dhaleshwari, different improved CETP treatment scenarios were simulated using INCA-Metals. This is the first modelling study to have targeted the impact of the relocated tanneries’ discharges at Savar on heavy metals in the Dhaleshwari River. This has involved setting up the process-based dynamic flow and quality model INCA Metals for the Dhaleshwari River catchment so that we could simulate a series of river reaches from its source down to the Buriganga. As shown in Figure 8, we model reaches of the river that correspond to the sampling locations and then add additional reaches down to the confluence point with the Buriganga River (Figure 7). Knowing the flow data for the river and the effluent discharge flow rate, then, by mass balance, we can estimate the effluent discharge concentrations. Using these concentrations, we can then evaluate the model performance, and the model simulates the concentrations of a range of metals with some accuracy, as shown in Figure 9. We can then use the model to simulate the impact of the discharge on the river system under a range of seasonal conditions and also a range of mitigation scenarios. Figure 10 shows a set of simulations along the river under low flow (i.e., dry season/worse case) conditions together with a set of treatment scenarios for the three heavy metals identified to have strongly impacted the Dhaleshwari during the sampling period. The concentrations along the river were simulated under the estimated low flow conditions to further highlight the gravity of the situation. It can be seen that levels of Cr, Cd and As are very high (black line) and well above international environmental standards if no treatment improvement is made (Table 7). These three key metals are of particular interest because of their impact on human health and fisheries. The results suggest that the CETP at Savar needs to operate at much higher levels of efficiency to control the river pollution, depicted by the other lines on the graphs representing improvements to the treatment plant under a range of scenarios from 65% to 95% reductions in discharge concentrations. In order to meet reasonable environmental conditions, an improvement of approximately 85% is required. Of course, a big issue with water quality modelling is that of uncertainty as decay parameters have not been measured explicitly but assumed from other studies [7,41]. Thus, the modelling results are subject to uncertainty and may vary along the river system under different assumptions. Another limitation is the lack of information regarding other discharges along the river, which were consequently not incorporated into the model used. Thus, the treatment improvement percentages suggested in our study are non-conservative and may vary, depending on many assumptions, but are nevertheless important in order to highlight the level of pollution control needed to safeguard the water quality of the Dhaleshwari.

6. Policy Recommendations

There are several clear messages for policy recommendations resulting from the modelling and monitoring analysis in this paper. The Savar plant may require a treatment upgrade from the current process solely based on a biological treatment to one that incorporates a chemical primary treatment to better handle the large concentrations of heavy metals released in the tanneries’ wastewater [39]. Such an upgrade could lead to significant improvements in water effluent quality. Operational management must be ensured to achieve a much higher treatment efficiency, and further attention should be reserved for the tanneries’ sludge management, which also impacts both the land and water quality of the area. New technology should be considered that can treat the particularly difficult waste from tanneries [46]. The Leather Working Group, a private multi-stakeholder group responsible for auditing and certifying leather producers, should consider adjusting their requirement, for tanneries to ensure wastewater treatment of their effluents, to the political, social and economic context of Bangladesh. It might be proactive to incorporate water-quality-oriented auditing of the Savar treatment plant to prompt action for improved control. For the government, it will be necessary to establish initiatives to enhance the uptake of clean production practices. The Department of Environment (DoE) as the main regulatory agency should liaise with other stakeholders to establish a realistic, time-bound approach to achieve SDG 6.3. Additionally, the DoE could establish ambient water quality standards based on the different pollutants released by local industries and the wastewater treatment capacity in Dhaka. It would be advisable to upgrade the capacity of the DoE to include real-time monitoring of effluent discharges and river water quality, as well as continuously ensuring public dissemination of water quality data. Long-term goals should target the analysis of riverine sediment and biota. The INCA-Metals model covers the whole of the Dhaka region, and this could form the basis of a real-time system to forecast pollution events and identify hotspots for decision-makers.

7. Conclusions

This paper provides the first modelling analysis of water quality changes in the Dhaleshwari since the relocation of tanneries from Hazaribag to Savar. It introduces the INCA-Metals model and its novel application in the Dhaka watershed. The paper then introduces the background of the study. It explains our field findings of increases in heavy metals concentrations from the upstream to downstream sites from the tanneries’ discharge point, expected to reach serious levels during the upcoming dry season, particularly for Cr, Cd and As. The results suggest that an 85% reduction in discharge concentrations is required to meet internationally approved standards for these particular metals. It should be emphasized that these findings need further investigation during peak and off-peak periods throughout the year in order to validate the results. The failure of the CETP to achieve sufficient pollution reduction calls into question the future viability of relying on central technical interventions, rather than considering other forms of recycling, reuse, or reduction of heavy metals in the process of leather tanning. For Bangladesh to achieve sufficient progress toward SDG 6, particularly target 6.3, addressing heavy metals pollution from leather tanneries will be a crucial step.
Effluent treatment facilities are considered a key step in ensuring reduced pollution of river systems from individual factories or from a cluster or sites, as is the case at Savar. At the moment, it appears that the treatment plant at Savar is underperforming. Sustainable Development Goal (SDG) 6 [1] for water and sanitation addresses the key issues of water quality and quantity. As in many countries, Bangladesh is striving to progress towards SDG 6.3 to ensure the adequate supply of good quality water for many uses including public supply, industrial supply and to restore rivers to good water quality and ecological status. Moving whole industrial sectors to new locations without putting in place working treatment plants is counterproductive as it merely pollutes yet another stretch of the river system. Urgent action needs to be taken to improve the treatment facilities at Savar.
The findings in this paper offer useful insights for other studies evaluating metals transport or the impact of particular industries upon the aquatic environment. Future work may usefully consider additional issues including speciation of metals, seasonality and other variables impacting water quality change over time and mechanisms for implementing improved monitoring systems.

Author Contributions

P.G.W. (Writing, conception, methodology), Z.M. (analytical work, writing), D.B. (Software development), G.B. (Data collation, modelling and writing), M.A.H. (conception, data and writing), R.P. (writing and policy), S.S. (Field sampling and lab analysis), P.H. (metals chemistry), C.P.N.R. (Lab analysis), L.J. (writing and data analysis), D.A. (Lab analysis). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by UK Aid from the UK Foreign, Commonwealth and Development Office (FCDO) for the benefit of developing countries (Programme Code 201880).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This document is an output from the REACH programme funded by UK Aid from the UK Foreign, Commonwealth and Development Office (FCDO) for the benefit of developing countries (Programme Code 201880). However, the views expressed, and information contained in it, are not necessarily those of or endorsed by FCDO, which can accept no responsibility for such views or information or for any reliance placed on them.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. United Nations. Transforming our World: The 2030 Agenda for Sustainable Development; United Nations: New York, NY, USA, 2015. [Google Scholar]
  2. Whitehead, P.G.; Bussi, G.; Hossain, M.A.; Dolk, M.; Das, P.; Comber, S.; Peters, R.; Charles, K.J.; Hope, R.; Hossain, S. Restoring water quality in the polluted Turag-Tongi-Balu river system, Dhaka: Modelling nutrient and total coliform intervention strategies. Sci. Total Environ. 2018, 631–632, 223–232. [Google Scholar] [CrossRef]
  3. Rampley, C.P.N.; Whitehead, P.G.; Softley, L.; Hossain, M.A.; Jin, L.; David, J.; Shawal, S.; Das, P.; Thompson, I.P.; Huang, W.E.; et al. River toxicity assessment using molecular biosensors: Heavy metal contamination in the Turag-Balu-Buriganga river systems, Dhaka, Bangladesh. Sci. Total Environ. 2020, 703, 134760. [Google Scholar] [CrossRef]
  4. Nordstrom, D.K. Advances in the Hydrogeochemistry and Microbiology of Acid Mine Waters. Int. Geol. Rev. 2000, 42, 499–515. [Google Scholar] [CrossRef]
  5. Johnson, D.B.; Hallberg, K.B. Acid mine drainage remediation options: A review. Sci. Total Environ. 2005, 338, 3–14. [Google Scholar] [CrossRef]
  6. Whitehead, P.G.; Prior, H. Bioremediation of acid mine drainage: An introduction to the Wheal Jane wetlands project. Sci. Total Environ. 2005, 338, 15–21. [Google Scholar] [CrossRef] [PubMed]
  7. Chapra, S.C.; Whitehead, P.G. Modelling impacts of pollution in river systems: A new dispersion model and a case study of mine discharges in the Abrud, Aries and Mures River System in Transylvania, Romania. Hydrol. Res. 2009, 40, 306–322. [Google Scholar] [CrossRef]
  8. Whitehead, P.G.; Wilson, E.; Butterfield, D. A semi-distributed integrated nitrogen model for multiple source assessment in catchments (INCA): Part I—Model structure and process equations. Sci. Total Environ. 1998, 210–211, 547–558. [Google Scholar] [CrossRef]
  9. Whitehead, P.G.; Wilson, E.; Butterfield, D.; Seed, K. A semi-distributed integrated flow and nitrogen model for multiple source assessment in catchments (INCA): Part II—Application to large river basins in south Wales and eastern England. Sci. Total Environ. 1998, 210–211, 559–583. [Google Scholar] [CrossRef]
  10. Wade, A.J.; Durand, P.; Beaujouan, V.; Wessel, W.; Raat, K.J.; Whitehead, P.G.; Butterfield, D.; Rankinen, K.; Lepisto, A. A nitrogen model for European catchments: INCA, new model structure and equations. Hydrol. Earth Syst. Sci. 2002, 6, 559–582. [Google Scholar] [CrossRef]
  11. Wade, A.J.; Whitehead, P.G.; Butterfield, D. The Integrated Catchments model of Phosphorus dynamics (INCA-P), a new approach for multiple source assessment in heterogeneous river systems: Model structure and equations. Hydrol. Earth Syst. Sci. 2002, 6, 583–606. [Google Scholar] [CrossRef] [Green Version]
  12. Bussi, G.; Darby, S.E.; Whitehead, P.G.; Jin, L.; Dadson, S.J.; Voepel, H.E.; Vasilopoulos, G.; Hackney, C.R.; Hutton, C.; Berchoux, T.; et al. Impact of dams and climate change on suspended sediment flux to the Mekong delta. Sci. Total Environ. 2021, 755, 142468. [Google Scholar] [CrossRef]
  13. Whitehead, P.G.; Sarkar, S.; Jin, L.; Futter, M.N.; Caesar, J.; Barbour, E.; Butterfield, D.; Sinha, R.; Nicholls, R.; Hutton, C.; et al. Dynamic modeling of the Ganga river system: Impacts of future climate and socio-economic change on flows and nitrogen fluxes in India and Bangladesh. Environ. Sci. Process. Impacts 2015, 17, 1082–1097. [Google Scholar] [CrossRef]
  14. Jin, L.; Whitehead, P.G.; Appeaning Addo, K.; Amisigo, B.; Macadam, I.; Janes, T.; Crossman, J.; Nicholls, R.J.; McCartney, M.; Rodda, H.J.E. Modeling future flows of the Volta River system: Impacts of climate change and socio-economic changes. Sci. Total Environ. 2018, 637–638, 1069–1080. [Google Scholar] [CrossRef]
  15. Wade, A.J.; Neal, C.; Butterfield, D.; Futter, M.N. Assessing nitrogen dynamics in European ecosystems, integrating measurement and modelling: Conclusions. Hydrol. Earth Syst. Sci. 2004, 8, 846–857. [Google Scholar] [CrossRef] [Green Version]
  16. Wade, A.J.; Butterfield, D.; Whitehead, P.G. Towards an improved understanding of the nitrate dynamics in lowland, permeable river-systems: Applications of INCA-N. J. Hydrol. 2006, 330, 185–203. [Google Scholar] [CrossRef]
  17. Lazar, A.N.; Butterfield, D.; Futter, M.N.; Rankinen, K.; Thouvenot-Korppoo, M.; Jarritt, N.; Lawrence, D.S.; Wade, A.J.; Whitehead, P.G. An assessment of the fine sediment dynamics in an upland river system: INCA-Sed modifications and implications for fisheries. Sci. Total Environ. 2010, 408, 2555–2566. [Google Scholar] [CrossRef]
  18. Futter, M.N.; Butterfield, D.; Cosby, B.J.; Dillon, P.J.; Wade, A.J.; Whitehead, P.G. Modeling the mechanisms that control in-stream dissolved organic carbon dynamics in upland and forested catchments. Water Resour. Res. 2007, 43. [Google Scholar] [CrossRef]
  19. Futter, M.N.; Erlandsson, M.A.; Butterfield, D.; Whitehead, P.G.; Oni, S.K.; Wade, A.J. PERSiST: A flexible rainfall-runoff modelling toolkit for use with the INCA family of models. Hydrol. Earth Syst. Sci. 2014, 18, 855–873. [Google Scholar] [CrossRef] [Green Version]
  20. Rankinen, K.; Kaste, Ø.; Butterfield, D. Adaptation of the Integrated Nitrogen Model for Catchments (INCA) to seasonally snow-covered catchments. Hydrol. Earth Syst. Sci. 2004, 8, 695–705. [Google Scholar] [CrossRef] [Green Version]
  21. Covington, A.D.; Wise, W.R. Tanning Chemistry: The Science of Leather—Anthony D Covington, William R Wise, 2nd ed.; The Royal Society of Chemistry: London, UK, 2020. [Google Scholar]
  22. Islam, S.; Ahmed, K.; Al-Mamun, H.; Masunaga, S. Assessment of trace metals in fish species of urban rivers in Bangladesh and health implications. Environ. Toxicol. Pharmacol. 2015, 39, 347–357. [Google Scholar] [CrossRef]
  23. Rahman, M.A. Characteristics of Major Industrial Liquid Pollutants in Bangladesh. Master’s Thesis, Bangladesh University of Engineering & Technology, Dhaka, Bangladesh, 1997. [Google Scholar]
  24. Saha, G.C. Assessment of Groundwater Contamination in Dhaka City from Tannery Waste. Master’s Thesis, Bangladesh University of Engineering and Technology, Dhaka, Bangladesh, 2001. [Google Scholar]
  25. Arias-Barreiro, C.R.; Nishizaki, H.; Okubo, K.; Aoyama, I.; Mori, I.C. Ecotoxicological characterization of tannery wastewater in Dhaka, Bangladesh. J. Environ. Biol. 2010, 31, 471–475. [Google Scholar]
  26. Ahmed, M.K.; Das, M.; Islam, M.M.; Akter, M.S.; Islam, S.; Al-Mansur, M.A. Physico-chemical properties of tannery and textile effluents and surface water of river Buriganga and Karnatoli, Bangladesh. World Appl. Sci. J. 2011, 12, 152–159. [Google Scholar]
  27. Siraj, S.; Islam, M.M.; Das, P.C.; Masum, S.M.; Jahan, I.A.; Ahsan, M.A.; Shajahan, M. Removal of Chromium from Tannery Effluent Using Chitosan-Charcoal Composite. J. Bangladesh Chem. Soc. 2012, 25, 53–61. [Google Scholar] [CrossRef]
  28. Karim, M.R.; Manshoven, S.; Islam, M.R.; Gascon, J.A.; Ibarra, M.; Diels, L.; Rahman, M.M. Assessment of an Urban Contaminated Site from Tannery Industries in Dhaka City, Bangladesh. J. Hazard. Toxic Radioact. Waste 2013, 17, 52–61. [Google Scholar] [CrossRef]
  29. Paul, H.L.; Antunes, A.P.M.; Covington, A.D.; Evans, P.; Phillips, P.S. Bangladeshi leather industry: An overview of recent sustainable developments. J. Soc. Leather Technol. Chem. 2013, 97, 25–32. [Google Scholar]
  30. Hashem, M.A.; Islam, A.; Mohsin, S.; Nur-A-Tomal, M.S. Green environment suffers by discharging of high-chromium-containing wastewater from the tanneries at Hazaribagh, Bangladesh. Sustain. Water Resour. Manag. 2015, 1, 343–347. [Google Scholar] [CrossRef]
  31. Jahan, M.; Akhtar, N.; Khan, N.; Roy, C.; Islam, R.; Nurunnabi, M. Characterization of tannery wastewater and its treatment by aquatic macrophytes and algae. Bangladesh J. Sci. Ind. Res. 2015, 49, 233–242. [Google Scholar] [CrossRef]
  32. Chowdhury, M.; Mostafa, M.G.; Biswas, T.K.; Mandal, A.; Saha, A.K. Characterization of the Effluents from Leather Processing Industries. Environ. Process. 2015, 2, 173–187. [Google Scholar] [CrossRef] [Green Version]
  33. Islam, S.; Islam, S.; Al-Mamun, H.; Islam, S.A.; Eaton, D.W. Total and dissolved metals in the industrial wastewater: A case study from Dhaka Metropolitan, Bangladesh. Environ. Nanotechnol. Monit. Manag. 2016, 5, 74–80. [Google Scholar] [CrossRef]
  34. Kabir, M.M.; Fakhruddin, A.N.M.; Chowdhury, M.A.Z.; Fardous, Z.; Islam, R. Characterization of tannery effluents of Hazaribagh area, Dhaka, Bangladesh. Pollution 2017, 3. [Google Scholar] [CrossRef]
  35. Abu Sayid Mia, M.A.S.; Nur-E-Alam, M.; Farid Ahmad, M.; Zahangir Alam, M.; Mafizur Rahman, M. Treatment of tannery wastewater by electrocoagulation technology. J. Sci. Innov. Res. 2017, 6, 129–134. [Google Scholar]
  36. Haque, M.A.; Chowdhury, R.A.; Chowdhury, W.A.; Baralaskar, A.H.; Bhowmik, S.; Islam, S. Immobilization possibility of tannery wastewater contaminants in the tiles fixing mortars for eco-friendly land disposal. J. Environ. Manag. 2019, 242, 298–308. [Google Scholar] [CrossRef]
  37. Abu Sayid Mia, M.; Uzzal Ali, M.; Nur-E-Alam, M.; Mafizur Rahman, M.; Zahangir Alam, M. Removal of Heavy Metals from Tannery Wastewater by using Sawdust and Spent Tea Leaves as an Adsorbent. Phys. Sci. Biophys. J. 2020, 4, 72–82. [Google Scholar]
  38. Hashem, M.A.; Hasan, M.; Momen, M.A.; Payel, S.; Nur-A-Tomal, M.S. Water hyacinth biochar for trivalent chromium adsorption from tannery wastewater. Environ. Sustain. Indic. 2020, 5, 100022. [Google Scholar] [CrossRef]
  39. Mimouni, Z. Assessment of Tanneries Wastewater Impact on Heavy Metals Post-Eid Al Adha in the Dhaleshwari River, Bangladesh. Master’s Thesis, University of Oxford, Oxford, UK, 2020. [Google Scholar]
  40. Islam, A.; Romić, D.; Akber, A.; Romić, M. Trace metals accumulation in soil irrigated with polluted water and assessment of human health risk from vegetable consumption in Bangladesh. Environ. Geochem. Health 2018, 40, 59–85. [Google Scholar] [CrossRef] [PubMed]
  41. Whitehead, P.G.; Bussi, G.; Peters, R.; Hossain, M.A.; Softley, L.; Shawal, S.; Jin, L.; Rampley, C.P.N.; Holdship, P.; Hope, R.; et al. Modelling heavy metals in the Buriganga River System, Dhaka, Bangladesh: Impacts of tannery pollution control. Sci. Total Environ. 2019, 697, 134090. [Google Scholar] [CrossRef] [PubMed]
  42. Kamal, M.M.; Malmgren-Hansen, A.; Badruzzaman, A.B.M. Assessment of pollution of the River Buriganga, Bangladesh, using a water quality model. Water Sci. Technol. 1999, 40, 129–136. [Google Scholar] [CrossRef]
  43. Ahmed, M.K.; Ahamed, S.; Rahman, S.; Haque, M.R.; Islam, M.M. Heavy Metals Concentration in Water, Sediments and their Bioaccumulations in Some Freshwater Fishes and Mussel in Dhaleshwari River, Banglades’. Terr. Aquat. Environ. Toxicol. 2009, 3, 33–41. [Google Scholar]
  44. Jolly, Y.N.; Choudhury, T.R.; Islam, A.; Suravi, N.I.; Sultana, M. Background Chemical Study of Relocated Hazaribagh Tannery Complex Environment, Savar. J. Bangladesh Acad. Sci. 2012, 36, 45–51. [Google Scholar] [CrossRef] [Green Version]
  45. Ahsan, M.A.; Siddique, M.A.B.; Munni, M.A.; Akbor, M.A.; Akter, S.; Mia, M.Y. Analysis of Physicochemical Parameters, Anions and Major Heavy Metals of the Dhaleshwari River Water, Tangail, Bangladesh. Am. J. Environ. Prot. 2018, 7, 29–39. [Google Scholar] [CrossRef] [Green Version]
  46. Ma, Y.; Liu, Y. Turning food waste to energy and resources towards a great environmental and economic sustainability: An innovative integrated biological approach. Biotechnol. Adv. 2019, 37, 107414. [Google Scholar] [CrossRef]
Figure 1. Maps showing Bangladesh above the Bay of Bengal and the main rivers in the Dhaka Region.
Figure 1. Maps showing Bangladesh above the Bay of Bengal and the main rivers in the Dhaka Region.
Sustainability 13 03556 g001
Figure 2. Overview of the INCA-Metals model structure.
Figure 2. Overview of the INCA-Metals model structure.
Sustainability 13 03556 g002
Figure 3. Hydrological stores and pathways in the landscape.
Figure 3. Hydrological stores and pathways in the landscape.
Sustainability 13 03556 g003
Figure 4. Cyanide and ammonium stores, transformations and pathways in the landscape.
Figure 4. Cyanide and ammonium stores, transformations and pathways in the landscape.
Sustainability 13 03556 g004
Figure 5. Metal stores, transformations and pathways in the landscape.
Figure 5. Metal stores, transformations and pathways in the landscape.
Sustainability 13 03556 g005
Figure 6. Stores, transformations and pathways in the river system.
Figure 6. Stores, transformations and pathways in the river system.
Sustainability 13 03556 g006
Figure 7. Hazaribag and Savar Tannery Locations on the Dhaleshwari and Buriganga Rivers in Dhaka.
Figure 7. Hazaribag and Savar Tannery Locations on the Dhaleshwari and Buriganga Rivers in Dhaka.
Sustainability 13 03556 g007
Figure 8. The sampling sites in the Dhaleshwari near Savar.
Figure 8. The sampling sites in the Dhaleshwari near Savar.
Sustainability 13 03556 g008
Figure 9. Simulated and observed metal concentrations (μg/L).
Figure 9. Simulated and observed metal concentrations (μg/L).
Sustainability 13 03556 g009
Figure 10. Heavy metals concentrations under low flow conditions with a set of savar treatment improvement strategies.
Figure 10. Heavy metals concentrations under low flow conditions with a set of savar treatment improvement strategies.
Sustainability 13 03556 g010
Table 1. Constants and variables in the landscape mass-balance equations.
Table 1. Constants and variables in the landscape mass-balance equations.
SymbolDefinitionUnits
User-supplied inputs as time series
peffHydrologically effective rainfallm3 s−1 km2
SMDSoil moisture deficitmm
θAAir temperatureoC
mnh4,inInput ammonium load (includes input from livestock, fertiliser and wet and dry deposition)kg N ha−1 day−1
User-supplied inputs as constants
βBase flow indexØ
TszSoil zone response timedays
TgzGroundwater zone response timedays
qsz,initialInitial outflow from soil zonem3 s−1 km2
SMDmaxThe maximum soil moisture deficit in the soilmm
deff,gwThe ‘effective’ depth of the groundwater zone equal to the depth of the active groundwater multiplied by the effective porositym
C1Ratio of soil retention volume to ‘easily available’ soil moistureØ
C2Fraction of total groundwater effective pore space occupied by waterØ
C3Cyanide volatilisation rate in the soil zonem day−1
C4Rate of decay of cyanide to ammonium in the soil zonem day−1
C5Rate of decay of cyanide to ammonium in the groundwater zonem day−1
C6Maximum temperature difference between summer and winteroC
C7Plant ammonium uptake ratem day−1
C8Soil zone ammonium nitrification ratem day−1
C9Groundwater zone ammonium nitrification ratem day−1
C10Day number associated with the start of the growing seasonØ
C11Rate of adsorption of metal to soil sedimentm day−1
C12Rate of adsorption of metal to aquifer matrixm day−1
tQ10Factor change in rate with a 10-degree change in temperatureØ
tQ10basBase temperature for the process rate at which the response is 1oC
Variables in the landscape equations
qszOutflow from soil zonem3 s−1 km2
qgwOutflow from groundwater zonem3 s−1 km2
VDSoil zone drainage volumem3 km2
VRSoil zone retention volumem3 km2
VgwGroundwater zone volumem3 km2
SSMDSoil moisture factorØ
SPGIPlant growth indexØ
mcn,szMass of cyanide in the soil zonekg km−2
mcn,gzMass of cyanide in the groundwater zonekg km−2
mnh4,szMass of ammonium in the soil zonekg N km−2
mnh4,gzMass of ammonium in the groundwater zonekg N km−2
mmetal,szMass of metal in the soil zonekg km−2
mmetal,gzMass of metal in the groundwater zonekg km−2
Table 2. Constants and variables in the river mass-balance equations.
Table 2. Constants and variables in the river mass-balance equations.
SymbolDefinitionUnits
User-supplied inputs as constants
TreachReach residence timedays
LReach lengthm
aDischarge-velocity parameterm−2
bDischarge-velocity parameterØ
C13Cyanide volatilisation ratem day−1
C14Rate of decay of cyanide to ammonium in the streamm day−1
C15Reach ammonium nitrification ratem day−1
C16Rate of metal sedimentation in the streamm day−1
Variables in the instream equations
qreach,outWater flow from the reachm−3 s−1
qreach,inWater flow to the reach from upstream, point source effluent, diffuse inputs from soil and groundwater zones and loss via abstractionm−3 s−1
mcn,reachMass of cyanide in the reachkg
mcn,reachMass of cyanide input to the reach from upstream, point source effluent, and diffuse inputs from soil and groundwater zoneskg
mcn,absMass of cyanide abstracted from the reachkg
mnh4,reachMass of ammonium in the reachkg
mnh4,reachMass of ammonium input to the reach from upstream, point source effluent, and diffuse inputs from soil and groundwater zoneskg
mnh4,absMass of ammonium abstracted from the reachkg
mmetal,reachMass of metal in the reachkg
mmetal,reachMass of metal input to the reach from upstream, point source effluent, and diffuse inputs from soil and groundwater zoneskg
mmetal,absMass of metal abstracted from the reachkg
ccn,reachConcentration of cyanide in the reachmg L−1
cnh4,reachConcentration of ammonium in the reachmg L−1
Cmetal,reachConcentration of metal in the reachmg L−1
VreachReach volumem3
Table 5. Impacts of various heavy metals on human health [39].
Table 5. Impacts of various heavy metals on human health [39].
Heavy MetalEffect
Arsenic (As)Nausea, diarrhoea, hyperpigmentation, hypopigmentation, skin cancer, bladder cancer, lung cancer, peripheral vascular disease, death
Cadmium (Cd)Pneumonitis, lung cancer, gastrointestinal tract disturbances,
liver and kidney dysfunction
Chromium (Cr)Dermatitis, acute renal failure, gastrointestinal haemorrhage,
lung cancer, asthma, nervous system damage
Copper (Cu)Gastrointestinal haemorrhage, liver and kidney damage,
mental disorders, anaemia, rheumatoid arthritis, nausea, hypertension, death
Mercury (Hg)Elemental fever, Vomiting, Diarrhea, Acute lung and kidney damage, Central nervous system damage, Nutritional disturbances,
Haemorrhagic gastritis
Nickel (Ni)Dermatitis, heart and liver damage, respiratory tract cancer
Lead (Pb)Impaired children development, memory loss, insomnia,
joints weakness, vomiting, diarrhea, abdominal pain, anaemia, kidney damage, miscarriages, coma, cardiovascular diseases, hypertension, death
Iron (Fe)Anaemia, heart disease, cancer, diabetes
Cobalt (Co)Sterility, hair loss, respiratory irritation, vomiting, diarrhoea, coma, death
Manganese (Mn)Sleepiness, weakness, leg cramps and paralysis, speech disturbances, emotional distress
Zinc (Zn)Gastrointestinal distress, diarrhoea
Table 6. Mean metal concentrations upstream and downstream of Savar discharge in wet and dry seasons with percentage change.
Table 6. Mean metal concentrations upstream and downstream of Savar discharge in wet and dry seasons with percentage change.
Mean Concentration Upstream μg/LMean Concentration Downstream μg/LEstimated Dry Season Mean Concentration Upstream mg/LEstimated Dry Season Mean Concentration Downstream mg/LPercentage Change %
Al5.545.560.720.720.36
Sr76.5079.379.9510.323.75
Ba23.1324.213.013.154.67
Cr0.180.640.020.08254.0
Co0.070.070.010.016.72
Zn6.688.170.871.0622.3
As2.813.010.370.397.27
Li2.322.450.300.325.50
Cd0.270.450.040.0665.0
Table 7. International Instream River Water Standards for metals (mg/L).
Table 7. International Instream River Water Standards for metals (mg/L).
USAEUINDIA
As0.340.10.2
Cd0.00180.0050.01
Cr 0.05
Cu 11.5
Fe 150
Pb0.0820.050.1
Ni0.47
Zn0.12515
Mn 1
Se0.0050.010.05
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Whitehead, P.G.; Mimouni, Z.; Butterfield, D.; Bussi, G.; Hossain, M.A.; Peters, R.; Shawal, S.; Holdship, P.; Rampley, C.P.N.; Jin, L.; et al. A New Multibranch Model for Metals in River Systems: Impacts and Control of Tannery Wastes in Bangladesh. Sustainability 2021, 13, 3556. https://0-doi-org.brum.beds.ac.uk/10.3390/su13063556

AMA Style

Whitehead PG, Mimouni Z, Butterfield D, Bussi G, Hossain MA, Peters R, Shawal S, Holdship P, Rampley CPN, Jin L, et al. A New Multibranch Model for Metals in River Systems: Impacts and Control of Tannery Wastes in Bangladesh. Sustainability. 2021; 13(6):3556. https://0-doi-org.brum.beds.ac.uk/10.3390/su13063556

Chicago/Turabian Style

Whitehead, Paul Geoffrey, Zineb Mimouni, Daniel Butterfield, Gianbattista Bussi, Mohammed Abed Hossain, Rebecca Peters, Shammi Shawal, Phillip Holdship, Cordelia Petra Nadine Rampley, Li Jin, and et al. 2021. "A New Multibranch Model for Metals in River Systems: Impacts and Control of Tannery Wastes in Bangladesh" Sustainability 13, no. 6: 3556. https://0-doi-org.brum.beds.ac.uk/10.3390/su13063556

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