Uncorrected Proof
© 2021 The Authors AQUA — Water Infrastructure, Ecosystems and Society Vol 00 No 0, 1 doi: 10.2166/aqua.2021.043
Influence of chemical fertilizers on arsenic mobilization in the alluvial Bengal delta plain: a critical review
Md. Shajedul Islam and M. G. Mostafa* Institute of Environmental Science, University of Rajshahi, Rajshahi 6205, Bangladesh *Corresponding author. E-mail: [email protected]
ABSTRACT
Arsenic contamination of alluvial aquifers of the Bengal delta plain causes a serious threat to human health for over 75 million people. The study aimed to explore the impacts of chemical fertilizer on arsenic mobilization in the sedimentary deposition of the alluvial Bengal delta plain. It selected ten comparatively higher affected Districts and the least affected two Divisions as a referral study site. The countrywide pooled concentration of arsenic in groundwater was 109.75 μg/L (52.59, 166.91) at a 95% confidence interval, which was double the national
guideline value (50 μg/L). The analysis results showed a strong positive correlation (r 0.5) of arsenic with NO3,NH4,PO4,SO4, Ca, and K,
where a portion of those species originated from fertilizer leaching into groundwater. The results showed that PO4 played a significant influ-
ence in arsenic mobilization, but the role of NO3,SO4, and NH4 was not clear at certain lithological conditions. It also showed that clay, peat, silt-clay, and rich microbial community with sufficiently organic carbon loaded soils could lead to an increase in arsenic mobilization. Finally, the study observed that the overall lithological conditions are the main reason for the high arsenic load in the study area.
Key words: aquifer sediment, arsenic mobilization, Bengal delta plain, fertilizer leaching, lithology, redox reaction
HIGHLIGHTS
• The groundwater of the Bengal delta basin is highly contaminated with arsenic.
• Arsenic significantly correlated with originated fertilizer species: NO3,NH4,PO4,SO4, Ca, and K. • Clayey, peaty, and silt-clay soil with heavy microbial and organic matters enhanced the arsenic mobilization by fertilizers leaching. • Lithological conditions are the major reasons for arsenic toxicity in groundwater.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY-NC-ND 4.0), which permits copying and redistribution for non-commercial purposes with no derivatives, provided the original work is properly cited (http://creativecommons.org/licenses/by-nc-nd/4.0/).
Downloaded from http://iwaponline.com/aqua/article-pdf/doi/10.2166/aqua.2021.043/929355/jws2021043.pdf?guestAccessKey=3101e41b-c033-4c5c-a3a1-55da9dd403ae by guest on 23 September 2021 Uncorrected Proof
AQUA — Water Infrastructure, Ecosystems and Society Vol 00 No 0, 2
GRAPHICAL ABSTRACT
INTRODUCTION Two Himalayan rivers, the Ganges and Brahmaputra, fall into the Bay of Bengal as a collective river transporting the largest sediment load. These rivers join in the central part of Bangladesh, with one more non-Himalayan River, the Meghna, which has created the largest delta in the world known as the Bengal delta. The Ganges–Brahmaputra River systems carry the largest sediment load in the world, about 80% of which is transported during the four rainy months (Goodbred & Kuehl 2000). More than 200 rivers and streams run through Bangladesh, with a mean annual discharge of water of about 38,000 m3/s (Haque et al. 2016), carrying over 2.4 109 MT of sediments every year before discharging into the Bay of Bengal (Talchabhadel et al. 2018). The huge alluvial sediment is carried away and deposited in basin areas and finally creates a large arable region in the world. Now this region has become the largest arsenic-affected area of the world. Bangladesh is a densely populated (2,890/mile2, positioned tenth in the world) and agrarian country in the world. Almost 90% of its lands are cultivable and 55% of the inhabitants are engaged directly in the agricultural sector (BBS 2019). Pre- sently, the application of chemical fertilizers has risen above 1,000% compared to the 1950s and has become a great threat to the environment and human health (Faruq 2018). The hydraulic conductivity of the deltaic alluvial and sandy land is sufficiently high and the residual portion of fertilizers can easily penetrate the topsoil by leaching. For this reason,
it can be assumed that through various chemical conversions, some water quality parameters such as pH, NO3,NO2,
PO4,NH4, K, Ca, etc. may increase in the sub-surface water. Several studies have shown that those parameters could influ- ence the release of arsenic from arsenic-rich sediment in aquifers (Anawar et al. 2003; Uddin & Kurosawa 2010; Kurosawa et al. 2013). The groundwater resource is the key factor for agricultural production in this fertile delta basin. There are 97% rural and 85% of urban people who directly use raw groundwater for drinking and other household purposes in Bangladesh (Mostafa et al. 2017; BBS 2018; MICS-B 2018). But this resource is under increasing threat from over-exploitation, population growth,
Downloaded from http://iwaponline.com/aqua/article-pdf/doi/10.2166/aqua.2021.043/929355/jws2021043.pdf?guestAccessKey=3101e41b-c033-4c5c-a3a1-55da9dd403ae by guest on 23 September 2021 Uncorrected Proof
AQUA — Water Infrastructure, Ecosystems and Society Vol 00 No 0, 3
rapid urbanization, and pollution from industries, domestic, and agricultural sources. Arsenic pollution is one of Bangladesh’s most severe environmental problems. According to the WHO, about 35–77 million people were regularly exposed to arsenic ( 50 μg/L) via drinking groundwater, and that was marked as the largest mass poisoning in history (Flanagan et al. 2012; Radfard et al. 2019; Saleh et al. 2019). The major source of arsenic in Bangladesh is geogenic since it is contained in the sedi- ments of the shallow Holocene aquifers of the Ganges delta basin (Ravenscroft et al. 2005). Millions of shallow tube-wells were drilled in this zone to supply drinking water, but the water was contaminated with a higher concentration of arsenic. Arsenic causes harm to the human body with various carcinogenic and non-carcinogenic problems. So far, approximately 40 thousand arsenicosis patients have been identified in Bangladesh (Johnston & Motaleb 2007). Several studies have been conducted on groundwater arsenic contamination in Bangladesh and neighboring West Bengal of India (Acharyya et al. 2000; Ravenscroft et al. 2001; UNICEF 2001; Chakraborti et al. 2010; CSISA-MI 2015; MICS-B 2018; World Bank Group 2019). Numerous studies, not only in Bangladesh but also in other Asian countries, have been con- ducted, including China (Sun 2004; Jiang et al. 2019; Sanjrani et al. 2019), India (Paul et al. 2015; Chandrashekhar et al. 2016; Shaji et al. 2020; Alsubih et al. 2021), Vietnam (Glodowska et al. 2020), and Japan (Hossain et al. 2016; Vongphuthone et al. 2017), where arsenic contamination in the environment reached an alarming position. However, the results have not produced the mitigation of arsenic poisoning. The mechanisms through which arsenic is released from soil or sediment into groundwater are still unknown and debatable. There are some hypotheses about the common mechanism of arsenic release in groundwater, i.e., the oxidation of arsenic-rich pyrite in the sediment, the reduction of iron oxo-hydroxide (FeOOH), and des- orbed arsenic from the sediment particles, and the ion exchange of adsorbed arsenic with phosphate from fertilizers. The bacterial activity in sediment was considered a driving factor to create reducing the environment in anoxic groundwater through oxidation of dissolved organic matter to the above mechanisms (Wang & Mulligan 2006). Further, soil conditions,
pH, Eh, NO3,HCO3,SO4, and NH4 of sediment and groundwater are considered to be influencing factors along with the above concept. Acharyya et al. (1999, 2000); Anawar et al. (2003), and Brömssen et al. (2014) reported that agricultural fer- tilizers may encourage arsenic mobilization by ion-exchange with P and N bearing ions resulting from fertilizers. According to
the lab experiments of Uddin & Kurosawa (2010) and Mahin et al. (2008), the NO3 and PO4 concentrations in sediment and groundwater were high in an arsenic-affected area in Bangladesh, and the source of these ions would be inorganic fertilizers. The mode of incidence and mobility of arsenic in sedimentary aquifers may be controlled by a complex interaction of micro- bially facilitated reactions and hydro-geochemical processes sensitive to site-specific hydrology and sediment type, as well as anthropogenic activities, such as application of chemical manure. Several studies were conducted in the last few decades but failed to explore a clear concept on the mechanisms of arsenic mobilization in the Ganges, Brahmaputra, and Meghna (GBM) delta aquifers and the role of chemical fertilizer in mobilization are not yet clearly understood. Hence, critical review studies are imperative to explore the geochemical mechanism controlling arsenic mobilization in the Bengal delta area. The study aimed to investigate the impact of chemical fertilizer on arsenic mobilization in the sedimentary deposition of the GBM delta system.
Hypothesis The application of chemical fertilizer is the major anthropogenic cause of arsenic contamination in groundwater of the Bengal delta plain (BDP).
Research questions 1. What is the country and District wise pooled concentration of arsenic in groundwater and what Districts are highly affected or relatively less affected? 2. How is the correlation of arsenic concentration with pH, nitrate, ammonium, phosphate, sulfate, calcium, and potassium in groundwater and how is their presence in the aquifer system due to the use of chemical fertilizer? 3. What is the role of those parameters in the mechanism of arsenic mobilization at certain lithological conditions? 4. Whether the chemical fertilizer is the main anthropogenic cause of high arsenic concentration in groundwater?
METHODS OF INVESTIGATION AND DATA SOURCE The mechanisms and influencing factors of arsenic mobilization from sediment to groundwater have been investigated world- wide since the 1980s, but a significant difference of those mechanisms between the papers before 2012 and after that time was
Downloaded from http://iwaponline.com/aqua/article-pdf/doi/10.2166/aqua.2021.043/929355/jws2021043.pdf?guestAccessKey=3101e41b-c033-4c5c-a3a1-55da9dd403ae by guest on 23 September 2021 Uncorrected Proof
AQUA — Water Infrastructure, Ecosystems and Society Vol 00 No 0, 4
observed. Thus, reviewed papers were divided into two groups: before and after 2012. Initially, the study identified over 150 published papers on arsenic contamination in aquifers but finally included 60 articles, 42 are before 2012 and 18 are after 2012. Here, the study follows the PRISMA statement (Moher et al. 2009) for every event of inclusion and exclusion criteria.
Bangladesh has 19 administrative Districts (old). The study pooled the district-wise value/concentration of pH, NO3,NH4,
PO4,SO4, Ca, and K of groundwater from secondary data and the pooled values of those parameters are included in Table 1. The study assumed that these ion concentrations were increased by inorganic fertilizer leaching and influenced the arsenic dissolution process in groundwater aquifers along with the geological processes. At the same time, the District-wise concen- tration of arsenic is shown in Table 1. It shows the build-up of correlation/association between those parameters with arsenic concentration. Also, the study selected the ten most arsenic-affected Districts out of a total of 64 (new), and their soil physio- graphy with other necessary data is shown later. Subsequently, Table 1 acts as referral material to find out the actual mechanism of arsenic mobilization in sediment and/or water body of the BDP. As well, the two least affected Divisions, Raj- shahi and Rangpur of Bangladesh are considered in the research hypothesis. For any statistical measurement in this study, XLSTAT and SPSS software was applied.
Study area: geomorphology, hydrology, and water type Tropical monsoon is a vital climatic feature of the BDP and has humic meteorological conditions. Average yearly rainfall is lowest in the north-west (1,430 mm) and increases both to the south-east (2,745 mm) and the north-east (4,178 mm). Despite the heavy rainfall, about 90% of river flows in Bangladesh originate in India, Nepal, and China, and those create a widespread floodplain. The Bengal delta, situated in Bangladesh, is the largest alluvial sedimentary basement of the world. It is typically made up of tertiary and quaternary deposits (Goodbred & Kuehl 2000). Sediment deposited on the Ganges–Brahmaputra river has created one of the largest deltas in the world (Mukherjee et al. 2001). The minerals of these areas are dominated
by detritus quartz and feldspar; and the main mineralization processes are Fe(OH)3 reduction, silicate weathering, and mixing with salty water (Mukherjee et al. 2009). The country has promising hydrogeologic conditions and tropical monsoon weather indicates potential storage of ground- water in aquifers. The unconsolidated near-surface Pleistocene to recent estuarine and fluvial sediments underlying most of the country usually form productive aquifers and deep semi-contracted to the uncontracted fluvial-deltaic deposit of Miocene age to the recent forms of many aquifers. The floodplains of the active/inactive delta plain and the major rivers of the GBM delta complex occupy 82% of the country (Figure 1)(BGS-DPHE 2001). The existing sub-surface geological studies indicate that the bulk of the good aquifers occur between 25 and 130 m depth but the upper portion of this layer is at high risk for arsenic contamination. The Pleistocene aquifers are typically free of arsenic, but this is not true for the Holocene aquifers, which are rich in arsenic concentration. The highest arsenic was detected in those Holocene aquifers, which are roughly 3,000 years old (Dowling et al. 2002). The water of the upper Holocene aquifers is nearly 100 years old and holds less arsenic concentration. However, the Holocene deposits are not homogenous and steady and are characterized by gaps and holes allowing the vertical extension of arsenic contamination (BGS-DPHE 2001). This explains the marked depth-dependence of the arsenic concentration. The maximum concentrations are found at 20–55 m, which correspond to young and shallow aquifers and the depth of aquifers is not enough to measure for waters being free of arsenic. In general, deeper aquifers (.120 m) produce groundwater with arsenic levels below the WHO standard of 10 μg/L (Bhattacharya et al. 2000). Arsenic is transported to the BDP by the phases adsorbed on suspended particulates because of the oxidation of primary arsenic-sul- fide minerals. It is typically adsorbed and co-precipitated with secondary Fe, Mn, and Al, solid phases in the form of arsenic (þ5) and finally deposited in the delta depending on the topography of the basin (Nickson et al. 1998, 2000). One of the pro- blems met in Bangladesh is the extreme inconsistency in the groundwater arsenic content between boreholes only 100 m away from each other (van Geen et al. 2011).
FERTILIZER CONSUMPTION IN BANGLADESH Based on chemical fertilizer consumption in the agricultural sector, Bangladesh is placed 18th out of 86 of the most world- wide consumer countries and 7th out of 32 Asian consumer countries (Index-mundi 2020). Fertilizer consumption (% of fertilizer production) in Bangladesh was 392.38 from 2016 to 2002. Its highest value over the past 14 years was 476.95 in 2015, while its lowest value was 121.21 in 2002. On the other hand, 175 kg/hectare chemical fertilizer was used in 2002 and this amount increased by 300 kg/hectare in 2016, but worldwide fertilizer consumption decreased by 2% during this same period (Figure 2)(WB 2017). The description of the most used fertilizers in Bangladesh is presented in Table 2.
Downloaded from http://iwaponline.com/aqua/article-pdf/doi/10.2166/aqua.2021.043/929355/jws2021043.pdf?guestAccessKey=3101e41b-c033-4c5c-a3a1-55da9dd403ae by guest on 23 September 2021 Table 1 | District-wise pooled concentration of arsenic and influencing species
Pooled concentration/value in the shallow aquifer Sl. Name of a No. District As (μg/L) pH NO3 (mg/L) NH4 (mg/L) PO4 (mg/L) SO4 (mg/L) Ca (mg/L) K (mg/L) Source
1 Dhaka 35.72 7.2 1.95 0.12 1.56 6.33 72.60 2.95 Alam et al. (2019); Mazeda et al. (2018) 2 Faridpur 141.34 7.53 3.98 – 1.92 5.27 103.8 5.01 Bodrud-doza et al. (2016); Islam et al. (2017a) 3 Mymensingh 31.76 8.19 0.06 0.96 0.16 0.16 41.2 2.36 Shahidullah et al. (2000); Nizam et al. (2016) 4 Narayanganj 178.12 6.9 4.71 1.05 3.55 8.9 68.0 3.58 Sarker & Zaman (2003); UNICEF-BBS (2009) norce Proof Uncorrected 5 Tangail 57.88 7.29 5.34 – 1.90 5.94 31.04 1.49 Proshad et al. (2017); Ahmed et al. (2021) 6 Pabna 102.11 6.94 7.77 – 5.86 2.95 75.0 7.30 Sarkar & Hassan (2006); Hossain et al. (2010) 7 Bogra 19.30 7.3 1.12 – 0.22 1.50 48.0 3.04 Shamsad (2010); Islam et al. (2017b) 8 Rangpur 18.30 7.5 1.10 0.09 0.38 1.90 12.6 2.10 Saha et al. (2019); UNICEF-BBS (2009) AQUA 9 Rajshahi 55.62 6.91 2.80 0.33 1.83 1.17 78.54 2.22 Mostafa et al. (2017); Rahman et al. (2017) –
10 Dinujpur 12.43 7.4 0.92 0.35 1.55 24.22 1.11 Islam et al. (2016); Bhuiyan et al. (2015) —
11 Chittagong 12.30 7.51 0.39 – 4.34 5.56 72.33 2.89 Mojumder & Islam (1995); Ahmed et al. 5 0, No 00 Vol Society and Ecosystems Infrastructure, Water (2010) 12 Comilla 359.76 7.1 4.50 9.29 8.14 5.34 111.6 9.52 Prodip et al. (2016); Brömssen et al. (2014) 13 Noakhali 362.38 7.2 3.91 2.80 3.50 5.39 72.00 10.90 Bhuiyan et al. (2010); Ahmed et al. (2010) 14 Sylhet 38.49 7.1 3.10 – 0.30 2.90 7.10 2.30 Islam et al. (2017c); UNICEF-BBS (2009) 15 Jessor 73.21 7.4 1.80 4.25 0.60 2.80 93.0 3.00 Uddin & Kurosawa (2010); Shaibur et al. (2012) 16 Khulna 156.40 7.5 6.20 – 4.98 10.25 101.0 17.0 Islam et al. (2017d) 17 Kushtia 83.55 7.5 3.10 5.71 5.05 3.19 88.50 5.50 Hossain et al. (2013); Islam & Mostafa (2021a, 2021b) 18 Barishal 333.10 7.58 4.23 – 5.3 46.86 44 9.46 Sukhen et al. (2017); UNICEF-BBS (2009) 19 Patuakhali 13.54 7.42 4.12 – 0.9 1.20 14.0 3.01 Islam et al. (2017e); UNICEF-BBS (2009)
aRef. DPHE/BGS/DFID (2000); Chakrabarti et al. (2010); CSISA-MI (2015); MICS-B (2018); and World Bank Group (2019).
Downloaded from http://iwaponline.com/aqua/article-pdf/doi/10.2166/aqua.2021.043/929355/jws2021043.pdf?guestAccessKey=3101e41b-c033-4c5c-a3a1-55da9dd403ae by guest on 23 September 2021 Uncorrected Proof
AQUA — Water Infrastructure, Ecosystems and Society Vol 00 No 0, 6
Figure 1 | Map of the Ganges–Brahmaputra–Meghna (GBM) plain.
Some raw materials of fertilizer production contain arsenic minerals, and those types of fertilizer contribute to the arsenic in the soil, which then leaches to groundwater. For example, phosphate fertilizers are possible sources of arsenic. The concen- tration of arsenic in the fertilizer will vary with the source rock for phosphate to produce the fertilizer. The arsenic content of phosphate fertilizer, widely used by farmers in Bangladesh, ranges from 2.42 to 6.79 mg/kg (Mahin et al. 2008; Uddin & Kurosawa 2010). Campos (2002) reported that NPK fertilizer carries 6.12 + 0.04 mg/kg arsenic. Thus, chemical methods contribute to arsenic in the soil system as well as in the groundwater body. Chemical fertilizers that contain toxic pollutants or species leach through the soils and enter into groundwater and can increase the arsenic mobilization in aquifers (Table 2).
Downloaded from http://iwaponline.com/aqua/article-pdf/doi/10.2166/aqua.2021.043/929355/jws2021043.pdf?guestAccessKey=3101e41b-c033-4c5c-a3a1-55da9dd403ae by guest on 23 September 2021 Uncorrected Proof
AQUA — Water Infrastructure, Ecosystems and Society Vol 00 No 0, 7
Figure 2 | Rate of the world and nationwide chemical fertilizer consumption.
Table 2 | Descriptions of mostly used chemical fertilizer in Bangladesh (Gowariker et al. 2009; Jones 2012; Rajani 2019)
Fertilizer class Type Example Chemical formula Probably released ions
þ Single nutrient Nitrogenous Urea NH2CONH2 H ,NH4, NO3
Ammonium nitrate NH4NO3 NH4,NO3
Calcium nitrate Ca(NO3)2 Ca, NO3 þ Phosphorus Single super phosphate, SSP Ca(H2PO4)2;14–18% P2O5 Ca, PO4,H
Gypsum CaSO4.2H2O Ca, SO4 þ Triple superphosphate Ca(H2PO4)2.H2O; 43–48% P2O5 Ca, PO4,H Potassium Muriate of potash (MoP) 95–99% KCl K, Cl Lime Lime, milk of magnesia, soda lime Cao, Ca(OH)2, soda Ca, Na, OH þ Multi-nutrient NP Mono-ammonium phosphate, MAP NH4H2PO4 NH4,PO4,H þ Di-ammonium phosphate, DAP (NH4)2HPO4 NH4,PO4,H þ NPK Ratio: 10-10-10 or 16-4-8 N þ K2O þ P2O5 NO3,NO2,K ,PO4,NH4
DISTRIBUTION OF ARSENIC IN BANGLADESH’S GROUNDWATER Arsenic, a category-1 carcinogenic element (Driscoll et al. 2000), occurs naturally in groundwater supplies through all parts of south-east Asia. About 25% of water wells of shallow aquifers in the Ganges basin of Bangladesh and West Bengal of India are severely contaminated by naturally occurring arsenic (Ravenscroft 2007a, 2007b). It was first recognized as a thought-provoking problem in Bangladesh in the 1990s. The arsenic concentration in most of the areas exceeded the maximum level of 10 μg/L suggested by the WHO. Over 75 million people, from 59 out of 64 Districts, were thought to be at risk of drinking water con- taminated by arsenic in Bangladesh (Safiuddin 2011). Every year, an assessed 43,000 people are affected by arsenic poisoning in the country (Jahan 2016). The authority has taken several steps and made strategies to try to address the calamity. Nevertheless, a countrywide movement and social mobilization activities by the government and NGOs to raise consciousness and knowledge levels among the public remain far below expectations. Several studies have been conducted in Southeast Asia and other countries of the world, including Bangladesh, but arsenic mobilization towards groundwater aquifers is not well understood. The complete scenario of the countrywide arsenic level of groundwater is summarized in Table 3.
MOBILITY OF ARSENIC IN BENGAL BASIN AQUIFER The mobility of arsenic in soils to groundwater depends on several factors including soil mineralogy, redox potential, pH, Eh, and the presence of other components that compete with arsenic for soil retention sites. Arsenic contamination in aquifers
Downloaded from http://iwaponline.com/aqua/article-pdf/doi/10.2166/aqua.2021.043/929355/jws2021043.pdf?guestAccessKey=3101e41b-c033-4c5c-a3a1-55da9dd403ae by guest on 23 September 2021 Uncorrected Proof
AQUA — Water Infrastructure, Ecosystems and Society Vol 00 No 0, 8
Table 3 | Statistics of the arsenic calamity in Bangladesh (Ref: Aggarwal et al. 2000; Fazal et al. 2001; Faruque 2006; MICS-B/BBS/UNICEF 2009; Chakraborti et al. 2010)
Basic data Value Percentage (%)
Area of Bangladesh 147,570 km2 – The population of Bangladesh (June 2019) 165.6 million – Number of Districts in Bangladesh (new) 64 – WHO arsenic potable water standard 10 μg/L – Bangladesh arsenic potable water standard 50 μg/L – Number of Districts surveyed 64 100 Number of Districts having arsenic above 50 μg/L in groundwater 59 92.2 Total area of affected 59 Districts 126,134 km2 85.5 Total population at risk .80 million .48 Potentially exposed population .35 million .21 Total number of patients suffering from arsenicosis .38,000 – The total number of patients who died 10 – Heavy contaminated districts: Chandpur (90%), Monshiganj (83%), Gopalganj (79%), Madaripur (69%), Noakhali (69%), Satkhira (67%), Comilla (65%), Faridpur (65%), Shariatpur (65%), Bagerhat (60%), and Meherpur (60%) Least affected districts: Thakurgaon, Dinajpur, Panchagarh, Nilphamari, Natore, Lalmonirhat, Patuakhali, and Barguna
can arise both naturally and anthropogenically and the method of the mechanism of arsenic mobilization be completely regional. This contamination in Bangladesh arises because of an unfortunate combination of four natural aspects: a source of arsenic (arsenic present in the aquifer sediment), adsorption of arsenic (adsorbed with metal in sediment), mobil- ization (arsenic release from the sediment to neighboring water), and transport (arsenic circulating in the groundwater) (Ravenscroft et al. 2005; Mostafa et al. 2010). In Bangladesh, two prevailing hypotheses expressing the mobilization of arsenic into groundwater are the oxidation of pyrite and the reduction of metal oxy-hydroxide (Mostafa et al. 2011). A study performed in Bangladesh showed a decline in groundwater level at 0.1–0.5 m year 1, indicating a reduction of the sto- rage in aquifer resulting from the unsustainable withdrawal of groundwater for both irrigation and domestic use (Zahid 2015). The over-extraction of groundwater and extensive application of inorganic fertilizer and arsenic-rich pesticides were the main suspects for the anthropogenic cause of arsenic contamination in Bangladesh, but recent investigations strongly oppose this thought by respecting the facts that are discussed below. 1. Oxidation of pyrite: Several studies conducted in the last two decades illustrated that iron-pyrite and arsenopyrite in the þ þ aquifer sediments oxidized to As3 and As5 , and were released into the groundwater aquifer (Chowdhury et al. 1999; Mazumdar 2000; Nickson et al. 2000; Mukherjee & Bhattacharya 2001; Ravenscroft et al. 2001; Hasan et al. 2007; Das et al. 2015). Furthermore, the oxidation might have occurred due to the entry of atmospheric oxygen into the aquifers during groundwater extraction through deep and shallow tube-wells.
þ þ ! þ 2 þ þ FeS O2 H2O Fe(OH)3 SO4 H
þ 3þ þ ! 2þ þ 2 þ þ þ : FeAsS Fe H2O Fe SO4 H H3AsO4(aq )
The rationality of this hypothesis has been keenly criticized by researchers since the groundwater in this region is: (i) anoxic
type, with negligible dissolved oxygen, (ii) near-neutral pH or moderate alkaline, (iii) categorized by low SO4 concentrations
with no relationship to arsenic, and (iv) characterized by very low concentration of NO3 (Ravenscroft et al. 2001; Pal et al. 2009). This mechanism was postponed at that time, and the following processes demonstrated later. 2. FeOOH reduction: Bacterial metabolism of dissolved organic matter in the sub-surface leading to anoxic conditions, thus leading to reduction of lepidocrocite mineral, FeO(OH), and then resulting in the release of adsorbed arsenic into the groundwater (McArthur et al. 2004; Swartz et al. 2004; Harvey et al. 2005).
Downloaded from http://iwaponline.com/aqua/article-pdf/doi/10.2166/aqua.2021.043/929355/jws2021043.pdf?guestAccessKey=3101e41b-c033-4c5c-a3a1-55da9dd403ae by guest on 23 September 2021 Uncorrected Proof
AQUA — Water Infrastructure, Ecosystems and Society Vol 00 No 0, 9
3. Redox process: Arsenic is released by the reductive dissolution of arsenic-rich FeO(OH), and this is a redox process that
occurs after microbial oxidation of organic matter has consumed DO, NO3, and NH4 (as the nutrient of microorganisms) (Zheng et al. 2004; Harvey et al. 2006; Mukherjee et al. 2006; Pal et al. 2009). 4. Ion exchange: Arsenate anion adsorbed to aquifer minerals is being displaced into the aqueous medium by a competitive exchange of phosphate anion resulting in arsenic contamination in groundwater (Mahin et al. 2008; Anawar et al. 2013; Shankar et al. 2014). The above mechanisms of 2, 3, and 4 were established after the 2000s. According to recently published papers, the avail- ability of dissolved organic matter, sufficient microbial community, an abundance of highly positive charge baring metal- minerals, high phosphate loads, and low Eh with moderately high pH in sediment or aquifer are the major criteria to boost arsenic concentration in a shallow groundwater system. Rather than anthropogenic causes, it is assumed that the geo- genic condition is a major cause for arsenic mobilization in the study regions. Leached chemical species such as nitrate, phosphate, and ammonium, which originate from inorganic fertilizer, can affect the above mechanisms 2, 3, and 4. Also, due to the application of fertilizer, pH and Eh may change indirectly in sediment. How the arsenic adsorption-desorption mechanisms in the Bengal delta aquifer are affected by chemical fertilizers are critically discussed in the next sections.
DISCUSSION About 85% of irrigation and drinking water is being extracted from shallow aquifers in Bangladesh (DPHE 2017). These aqui- fers mainly consist of clayey, silt-clay, or sandy sediment. Several studies have confirmed that aquifers enriched in Fe-, Mn-,
Al-oxides, PO4,NH4, and NO3 with organic matter and deposited in the middle part of shallow aquifers contain moderately high concentrations of arsenic, whereas the sediments of deep aquifers contain a low level of arsenic (Nath et al. 2008; Anawar et al. 2011; Uddin et al. 2011; Mohammadi et al. 2020). Among the total 64 administrative Districts (new) of Ban- gladesh, over 60% of groundwater samples of ten Districts exceed the national arsenic standard concentration of 50 μg/L (Table 3), and those Districts are considered as reference sites of the present study. These study sites represent three of five major geological units of the GBM delta: young Holocene terrace, fluvial flood plain, and delta plain (BBS 2018). The soil physiography and land types of most of the vulnerable areas are of relatively young sediment formation (∼3,000 years old) and these follow the order clayey . silt-clay . silt-loam, and most of the areas are river basins with a low-lying area having clay-muddy land (Table 4). Thus, it can be seen that high arsenic levels in groundwater or sediment are often associ- ated with geologically young sediments (Holocene aquifer) and flat low-lying zones where groundwater flow is slow. Such conditions favor the reductive dissolution of iron oxides and/or desorption of arsenic from metal oxides associated with sedi- ments. Several components are leached into the soil through the decomposition of chemical fertilizers. For this reason, the
concentration of NO3,NH4,PO4,SO4, Ca, and K in shallow groundwater is found in elevated levels because of fertilizers used in agriculture to increase crop yields. In Bangladesh, silt-clay, clayey, peaty-clay, and peat sediments are mainly located within 10–15 m of the surface, followed by sandy sediment going deeper. This part of the soil layer is typically important to desorption of arsenic and around this layer the impacting factors are applying their actions (Uddin & Kurosawa 2010). Stollenwerk et al. (2007), Mahin et al. (2008), Uddin & Kurosawa (2010), and Anawar et al. (2011) have shown by their exper- iments conducted in a separate region of the country that the N and P containing ions generated from nitrogen and phosphate fertilizer accumulate mainly in peat or clayey zone of the soil. They also found the concentration of arsenic in this part of the soil profile is significantly higher, and the arsenic loads are strongly correlated with those ions. The present study and previous reports (e.g., Sutton et al. 2009) show that the groundwater in the upper aquifers of the
Holocene delta is characterized by high levels of As, Fe, PO4,NH4,NO3, K, total organic carbon with a high pH, EC, and low value of Eh, whereas the deeper aquifers in the Pleistocene terrace deposits are characterized by a low level of
arsenic, dissolved carbon, NH4, Fe, and PO4, along with a relatively low pH and EC. Like the present study, many prior inves- tigations have noticed similar occurrences of high and low arsenic concentrations in groundwater depending on the depth of the aquifer and the redox conditions in the different areas of Bangladesh (Harvey et al. 2006; Van Geen et al. 2011). Arsenic was inversely correlated with Eh values in the upper aquifer, whereas no relationship was found in the deeper aquifer (Zheng et al. 2004; McArthur 2018). The results suggested that arsenic mobilization was linked to the improvement of reducing conditions.
Here, the study developed a linear correlation matrix of arsenic with such parameters as pH, NO3,NH4,PO4,SO4, Ca, and K in shallow aquifer water of Bangladesh, and this is shown in Table 5. The results indicated that the values of the correlation
Downloaded from http://iwaponline.com/aqua/article-pdf/doi/10.2166/aqua.2021.043/929355/jws2021043.pdf?guestAccessKey=3101e41b-c033-4c5c-a3a1-55da9dd403ae by guest on 23 September 2021 Uncorrected Proof
AQUA — Water Infrastructure, Ecosystems and Society Vol 00 No 0, 10
Table 4 | Ten most arsenic-affected Districts (new) of Bangladesh with soil properties and other necessary data (Table constructed from the report of BBS 2018)
% Samples exceed Net Most vulnerable the national standard Sub-soil nature Cultivable irrigated Sl.No Districts of arsenica Physiographic unit Soil physiography/Land type (acidic/alkaline) land area (%) land (%)
1 Chandpur 90 Lower Meghna Clayey – 20%; Silt-loam – Neutral, pH up 61.95 61.30 river floodplain 73%; Other – 7% to 7.1 2 Monshiganj 83 Meghna estuarine Clayey – 21%; Silt-clay – Neutral, pH up 57.66 57.47 floodplain 33%; Silt-loam – 24%; to 7.2 Other – 22% 3 Gopalganj 79 Gopalganj- Clayey – 82%; Other – Alkaline, pH 57.58 64.46 Khulna Bill 18% up to 7.7 4 Madaripur 69 Meghna estuarine Clayey – 21%; Loamy – Alkaline, pH 61.15 45.42 floodplain 33%; Silt-loam – 24%; up to 7.7 Other – 22% 5 Noakhali 69 Meghna estuarine Silt-loam – 45%; Other – Alkaline, pH 37.10 32.84 floodplain 55% up to 7.6 6 Satkhira 67 Gangues tidal Clayey – 18%; Silt-loam – Alkaline, pH 41.68 61.78 floodplain 78%; Other – 4% up to 7.9 7 Comilla 65 Meghna estuarine Clayey – 21%; Loamy – Neutral, pH up 65.88 76.54 floodplain 33%; Silt-loam – 24%; to 7.2 Other – 22% 8 Faridpur 65 Low Gangues Clayey – 31%; Silt-clay – Alkaline, pH 65.68 53.29 River 29%; Silt-loam – 13%; up to 7.7 floodplain Other – 27% 9 Shariatpur 65 Low Gangues Clayey – 31%; Silt-clay – Alkaline, pH 51.38 32.34 River 29%; Silt-loam – 13%; up to 7.7 floodplain Other – 27% 10 Bagerhat 60 Gangues tidal Clayey – 18%; Silt-loam – Neutral/ 29.35 21.99 floodplain 78%; Other – 4% alkaline, pH up to 7.5
aRef. DPHE/BGS/DFID (2000); Chakrabarti et al. (2010); CSISA-MI (2015); MICS-B (2018); and World Bank Group (2019).
Table 5 | Correlation matrix of variables with countrywide arsenic concentration
Parameters As pH NO3 PO4 SO4 Ca K
As 1 pH 0.1801 1
NO3 0.4798 0.4756 1
PO4 0.6899 0.2528 0.5937 1
SO4 0.5528 0.1177 0.2455 0.3921 1 Ca 0.4448 0.1603 0.2925 0.6117 0.0078 1 K 0.7089 0.0147 0.5859 0.6904 0.4051 0.528 1
Bold numbers indicate a strong correlation with each other.
coefficient (r) are about 0.5 or higher, and the concentration of these ions is strongly correlated with arsenic at p ¼ ,0.05, except for the pH. As well, an inter-component correlation was measured among the eight different variables through the principal component analysis (PCA) method (Table 6). The result showed a total variance of 65.35% and 17.159% for PC1 and PC2, respectively, with an Eigenvalue .1, as determined by two PCs of R-mode. The positive and negative values in PCA clarified that the water samples were affected or unaffected by the presence of extracted loads on a specific
Downloaded from http://iwaponline.com/aqua/article-pdf/doi/10.2166/aqua.2021.043/929355/jws2021043.pdf?guestAccessKey=3101e41b-c033-4c5c-a3a1-55da9dd403ae by guest on 23 September 2021 Uncorrected Proof
AQUA — Water Infrastructure, Ecosystems and Society Vol 00 No 0, 11
Table 6 | Principal component (two components extracted) loadings of the variables (sorted by size)
Variables PCA1 PCA2
PO4 0.934 0.132
NO3 0.919 0.338 As 0.913 0.234 K 0.827 0.350 Ca 0.759 0.088
NH4 0.743 0.559
SO4 0.666 0.559 pH 0.649 0.656 %Variance 65.350 17.159 %Cumulative 65.350 82.509 Eigenvalue 5.228 1.373
Bold numbers indicate strong loading.
constituent. Except for the pH, PO4,NO3, As, K, Ca, and NH4 showed a strong association (bold font) with PC1. However,
only NH4 and pH are strongly loaded for PC2. Again, this result confirmed the strong association among these parameters that revealed the inter-dependency with each other. Hence, the study assumed that the above-mentioned ions were respon- sible for the desorption of excess arsenic from its minerals into groundwater. However, various factors such as Fe, Al, and Mn load along with Eh, organic carbon, and microbial activities in the sediment can affect this process and are discussed earlier. The mechanism of influencing the above parameters (generated through fertilizer leaching) on arsenic dissolution is dis- cussed below.
Arsenic vs phosphorus
Arsenic mobility in groundwater is a function of P concentration. The pooled concentration of PO4 in a countrywide shallow
aquifer is shown in Table 1. A strong positive correlation (r ¼ 0.69; p ¼ ,0.05 at 95% CI) between arsenic and PO4 was observed according to Figure 3. Nonetheless, the actual r value is greater than this observed value because some portion
of free PO4 in groundwater can be settled out by metal-PO4 precipitation, but did not occur for arsenic. Campos (2002) and Anawar et al. (2011) conducted studies in a region of Brazil and Bangladesh, respectively, and they also found the
same correlation between arsenic and PO4 in the groundwater. Under certain conditions (high pH and/or low iron and
Eh), it may increase arsenic solubility in soils contaminated with P-containing fertilizers because of competitive PO4-AsO4
exchange (Roberts et al. 2004; Tyrovola et al. 2006). An experimental study has shown that a 0.1M KH2PO4 solution can leach as much as 150 mg/kg of arsenic from Fe-coated sediments (Loeppert et al. 2002). Most likely, arsenic (þ5) was des- 2 orbed from the surfaces of the Fe coatings and replaced by the sorption of HPO4 or H2PO4 . The elemental P and As are placed in the same group of the periodic table, and these have similar electronic configuration and properties. Besides, phos- phates and arsenates showed almost a similar physicochemical behavior in soils. These two ions compete for sorption sites on the sediment or soil granules. Phosphate may significantly reduce the adsorption of arsenic because it typically presents
greater concentrations than arsenates. Also, the charge density of PO4 is higher than arsenate and for this reason, PO4 has a greater attraction to highly positive Fe, Al, and/or Mn ions. Thus, it can easily replace the adsorbed arsenate from these metal mineral/ore surfaces. Manning & Goldberg (1997) stated that phosphorus and arsenic adsorption envelopes
exhibited a similar adsorption at pH 5.0 on kaolinite minerals and failed suddenly at pH . 6.5. The transport of PO4 in the sub-surface probably depends on the physical and chemical conditions and depth of the aquifer (Stollenwerk et al.
2007). For instance, higher pH in groundwater causes decreased sorption of PO4 on aquifer sediments. Reza et al. (2010)
and Uddin & Kurosawa (2010) confirmed that if PO4 was coming from fertilizers, its concentration decreased as the soil
depth increased. The capacity of the sub-surface soils to adsorb PO4 applied as fertilizers depends, in part, on the PO4 load applied (Uddin & Kurosawa 2010). From the above discussion, it may be concluded that in the ten most arsenic-affected Districts of Bangladesh, phosphates fertilizer may be one of the major causes of heavy arsenic loads in groundwater. It is not the only cause but also depends on
Downloaded from http://iwaponline.com/aqua/article-pdf/doi/10.2166/aqua.2021.043/929355/jws2021043.pdf?guestAccessKey=3101e41b-c033-4c5c-a3a1-55da9dd403ae by guest on 23 September 2021 Uncorrected Proof
AQUA — Water Infrastructure, Ecosystems and Society Vol 00 No 0, 12
Figure 3 | Scattered plot of pH, NO3,NH4,PO4,SO4, and K with arsenic concentration.
local soil characteristics. If the soil is peaty, muddy, or clayey, then it can retain the phosphates, otherwise not. Sandy soil cannot adsorb this ion. This argument was supported by several reports (Reza et al. 2010; Uddin & Kurosawa 2010). The soils of the other two Divisions, Rajshahi and Rangpur of Bangladesh, are sandy/sandy-loam (Table 7). Despite relatively more application of chemical fertilizer in both these divisions, they are much less affected by arsenic. In these areas, phos- phate-fertilizer did not accumulate to the upper layer of the soil and cannot participate in the process of arsenic mobilization. Reza et al. (2010) also created a plot with r ¼ 0.7, p , 0.05 for the same GBM floodplain areas, which means high arsenic concentrations are associated with the higher range of phosphate concentration. Phosphorus is able to replace arsenic in soils and sediments, and thus the groundwater of the areas can have a higher possibility of contamination with arsenic. Another report showed that phosphate-fertilizer uses lower amounts in the ten highly arsenic-affected Districts of Bangladesh (BBS 2018). Therefore, it can be said that the geogenic conditions are the primary and vital cause of arsenic con- tamination in the groundwater of the Bengal delta plain rather than phosphate concentration.
Arsenic vs nitrogen It is well known that nitrate may be engaged in numerous microbe-mediated anoxic redox reactions in arsenic-rich soils and aquifers. Nitrate can reduce into nitrite and then ammonium, which coupled with the oxidations of several electron donors such as Fe(þ2), As(þ3), sulfide, lactate, pyruvate, acetate, and other organic/inorganic materials (Jiang et al. 2013; Zhang
Downloaded from http://iwaponline.com/aqua/article-pdf/doi/10.2166/aqua.2021.043/929355/jws2021043.pdf?guestAccessKey=3101e41b-c033-4c5c-a3a1-55da9dd403ae by guest on 23 September 2021 Uncorrected Proof
AQUA — Water Infrastructure, Ecosystems and Society Vol 00 No 0, 13
Table 7 | Comparable dataset of the ten most prevalent Districts and Rajshahi and Rangpur Divisions
Pooled value
Samples exceed the Rate of % Area of Soil national Cultivable Net inorganic Mean Mean the physiography/ standard of Sub-soil land area irrigated fertilizer uses soil DOC soil Eh Location country Land type arsenic nature (%) land (%) (%) (mg/L) (V)
Ten most prevalent 16.62 Clayey/silt- Over 70% Neutral- 53.46 52.33 87 4.63 þ0.16 Districts (see clay/silt- alkaline, Table 5) loam pH: 7.1– 7.9 Rajshahi and 23.14 Loamy/sandy- Below 10% Acidic- 72.12 58.21 113 2.19 þ0.34 Rangpur loam/silt- neutral, Divisions loam pH: 4.5– (consisting of 16 7.2 Districts)
et al. 2017). The desorption of arsenic from FeO(OH) mineral in soil or sediment, nitrate, nitrite, and ammonium has a com- plex, complicated, and debatable functioning capacity. Some influencing factors such as organic matter, microbial community, soil profile, pH, and Eh also affect the efficiency of the functioning capacity of these ions. Nitrate and ammonium may accumulate in subsoil by the leaching of nitrogen-rich fertilizer through various chemical decomposition. Like phos-
phate, the pooled concentration of both these ions is significantly correlated (NO3: r ¼ 0.48, and NH4: r ¼ 0.71 at p , 0.05; 95% CI) with arsenic in groundwater of the Bengal delta plain (Figure 3). Nickson et al. (2000), Kiyoshi et al. (2008), Kurosawa et al. (2008, 2013), and Uddin & Kurosawa (2010) carried out their investigations in separate regions of Bangla-
desh and found a positive correlation of NO3,NH4, and pH with arsenic in groundwater, but the relationships among them
were not easy to understand. They assumed that a rise in the NO3 and NH4 concentration enhances microbial activity through the consumption of dissolved oxygen and then arsenic was released from FeO(OH)-rich sediments to the surround-
ing water body in a reducing condition. In this procedure, the NO3 and NH4 acted as a nutrient for microbes and the þ þ groundwater became a reducing environment. This reducing condition of the Fe (Fe3 to Fe2 ) was boosted by microbial metabolism of sedimentary organic matter and makes arsenic mobilization easier. This assumption may support the findings
of the present study that found the significantly positive correlation of both NO3 and NH4 with arsenic in the countrywide groundwater samples. They have also identified the source of these ions in groundwater mainly as N-fertilizer, based on the δ15N analysis. Therefore, the application of N-fertilizer appears to have a positive effect on the arsenic concentration in
groundwater. Uddin & Kurosawa (2010) stated that the deposition of NH4 in peat/peat-clay sediment of shallow aquifer
(10–15 m) was higher than deep aquifer like PO4. They measured up to 60% organic matter in the shallow layer and assumed that this large amount of organic carbon may be able to raise the total nitrogen. Also, in the peaty-clay layer, organic matter can oxidize in the presence of micro-organisms and serve electrons which are received by nitrate to produce nitrite ion and
then NH4 with N2 by the complete redox process. This combined amount of NH4 may enhance the arsenic release in sediment.
In the recent past, several studies established the different rules of NO3 and NH4 to the mobilization of arsenic in soil or water bodies by their experiments (Smith et al. 2017; Park et al. 2018; Shakya & Ghosh 2019; Yao et al. 2019; Zhu et al. 2020). The result of these investigations differs from the above mechanism. Mayorga et al. (2013) stated that the arsenic in water samples of Spain decreased along with the increase in nitrate concentrations and assumed that the increased use of nitrogenous fertilizer and pig manure in agricultural activities increased the nitrate content in groundwater, which would favor the precipitation of FeO(OH) and arsenic adsorption. As well, Smith et al. (2017) reported that the accumulation of nitrate in anoxic groundwater improved bacterial mediated processes in the aquifer that ultimately result in arsenic immobil-
ization. Another study illustrated that the NO3 concentration increase in groundwater results in a decrease in the total arsenic concentration through adsorption of arsenic and iron sulfide precipitation (Shakya & Ghosh 2019), although, there was a fixed amount of metal-oxyhydroxides (MOOH) yielded following the nitrate additions, which resulted in an increased arsenic level like those measured in the groundwater before nitrate addition. The nitrate itself was removed in the process, being
reduced first to N2O gas and then to N2 gas.
Downloaded from http://iwaponline.com/aqua/article-pdf/doi/10.2166/aqua.2021.043/929355/jws2021043.pdf?guestAccessKey=3101e41b-c033-4c5c-a3a1-55da9dd403ae by guest on 23 September 2021 Uncorrected Proof
AQUA — Water Infrastructure, Ecosystems and Society Vol 00 No 0, 14