Identification of hydrogeochemical processes and pollution sources of groundwater resources in the Marand plain,northwest of Iran

The main aims of the present study are to identify the major factors affecting hydrogeochemistry of groundwater resources in the Marand plain, NW Iran and to evaluate the potential sources of major and trace elements using multivariate statistical analysis such as hierarchical clustering analysis (HCA) and factor analysis (FA). To achieve these goals, groundwater samples were collected in three sampling periods in September 2013, May 2014 and September 2014 and analyzed with regard to ions (e.g., Ca²⁺, Mg²⁺, Na⁺ and K⁺, HCO₃ ^?, SO₄ ^2?, Cl^?, F^? and NO₃ ^?) and trace metals (e.g., Cr, Pb, Cd, Mn, Fe, Al and As). The piper diagrams show that the majority of samples belong to Na–Cl water type and are followed by Ca–HCO₃ and mixed Ca–Na–HCO₃. Cross-plots show that weathering and dissolution of different rocks and minerals, ion exchange, reverse ion exchange and anthropogenic activities, especially agricultural activities, influence the hydrogeochemistry of the study area. The results of the FA demonstrate that 6 factors with 81.7% of total variance are effective in the overall hydrogeochemistry, which are attributed to geogenic and anthropogenic impacts. The HCA categorizes the samples into two clusters. Samples of cluster C1, which appear to have higher values of some trace metals like Pb and As, are spatially located at the eastern and central parts of the plain, while samples of cluster C2, which express the salinization of the groundwater, are situated mainly westward with few local exceptions.     

Analysing the influence of Farakka Barrage Project on channel dynamics and meander geometry of Bhagirathi river of West Bengal,India

Farakka Barrage Project (FBP) constructed in 1975 on the river Ganga in West Bengal, India has largely altered hydro-geomorphic processes inducing a different kind of impetus in channel dynamics especially channel oscillation of the Bhagirathi river. Controlled hydro-geomorphological regime drastically altering natural regime has significantly influenced channel process with resultant change in channel geometry. Geospatial analysis using topographical maps of 1927, 1954 and 1974 along with satellite images of 1990, 2004, and 2014 clearly portrays rapidity in channel evolution in the post-Farakka period which is explicitly proved from four meander cut-offs after 1975. The meander geometry analysis of 66 loops of lower Bhagirathi river during 1927–2014 in terms of four basic parameters (radius of curvature, wavelength, amplitude, and channel length) and four indices (sinuosity index, radius–wavelength ratio, meander shape index and meander form index) depicts telltale pattern of channel straightening with increased number of wide and open loops after the construction of barrage. The average sinuosity of the loops has increased from 2.27 in 1927 to 2.35 in 1974 and dwindles to 1.80 in 2014. This distinct outcome induced by FBP is also found for other indices explained with special reference new hydraulic regime in post-Farakka period.     

Heterogeneity analysis of reservoir porosity and permeability in the Late Ordovician glacio-fluvial Sarah Formation paleovalleys,central Saudi Arabia

The Late Ordovician glacio-fluvial Sarah Formation is an important tight gas reservoir target in Saudi Arabia. This study uses statistical methods to characterize the petrophysical heterogeneity of the paleovalleys of the Sarah Formation that crop out in central Saudi Arabia. Four paleovalleys were studied: Bukayriyah, Hanadir, Sarah, and Khanasir Sarah. Several lithofacies were identified in each that vary in texture, porosity, permeability, and facies abundance that reflect periods of ice advance and retreat. The heterogeneity analysis is based on three statistical measures, namely, the coefficient of variation, the Dykstra-Parsons coefficient, and the Lorenz coefficient. The coefficient of variation values is in the 0.62–1.94 range, indicating an extremely heterogeneous distribution. The Dykstra-Parsons coefficient values are in the 0.56–0.88 range, suggesting very high to extremely high heterogeneity in the reservoirs. The Lorenz coefficient correlates well with the Dykstra-Parsons coefficient for paleovalleys of the Sarah Formation. The heterogeneity parameters studied here indicate that the outcrops of Sarah Formation paleovalleys represent heterogeneous to very heterogeneous reservoirs, which may be attributed to complex depositional and diagenetic variations that have affected the porosity and permeability distribution.     

Climatology of landfalling tropical cyclones in Bangladesh 1877–2003

Bangladesh is highly susceptible to tropical cyclones. Unfortunately, there is a dearth of climatological studies on the tropical cyclones of Bangladesh. The Global Tropical Cyclone Climatic Atlas (GTCCA) lists historical storm track information for all the seven tropical cyclone ocean basins including the North Indian Ocean. Using GIS, tropical cyclones that made landfall in Bangladesh during 1877–2003 are identified and examined from the climatological perspective. For the convenience of study, the coast of Bangladesh is divided into five segments and comparisons are made among the coastal segments in terms of cyclone landfall and vulnerability. There is a large variability in the year-to-year occurrence of landfalling tropical cyclones in Bangladesh. Most of the tropical cyclones (70%) hit in the months of May–June and October–November generally show the well-known pattern of pre- and post-monsoon cyclone seasons in that region.     

Groundwater contamination with arsenic in Sherajdikhan,Bangladesh: geochemical and hydrological implications

An integrated study has been carried out to elucidate the distribution and occurrence of arsenic in selected groundwater samples in the area of Sherajdikhan, Bangladesh. Arsenic and other parameters (T, pH, EC, Na⁺, K⁺, Ca²⁺, Mg²⁺, Cl⁻, NO₃ ⁻, SO₄ ²⁻, HCO₃ ⁻, PO₄ ³⁻, Fe, Mn and DOC) have been measured in groundwater samples collected from shallow/deep tube wells at different depths. Hydrogeochemical data suggest that the groundwaters are generally Ca–Mg–HCO₃ and Mg–Ca–HCO₃ types with bicarbonate (HCO₃ ⁻) as the dominant anion, though the other type of water has also been observed. Dissolved arsenic in groundwater ranged from 0.006 to 0.461 mg/l, with 69% groundwater samples exceeded the Bangladesh limit for safe drinking water (0.05 mg/l). Correlation and principal component analysis have been performed to find out possible relationships among the examined parameters in groundwater. Low concentrations of NO₃ ⁻ and SO₄ ²⁻, and high concentrations of DOC, HCO₃ ⁻ and PO₄ ³⁻ indicate the reducing condition of subsurface aquifer where sediments are deposited with abundant organic matter. Distinct relationship of As with Fe and Mn, and strong correlation with DOC suggests that the biodegradation of organic matter along with reductive dissolution of Fe–Mn oxyhydroxides has being considered the dominant process to release As in the aquifers studied herein.     

Numerical simulation of stress distributions and displacements around an entry roadway with igneous intrusion and potential sources of seam gas emission of the Barapukuria coal mine, NW Bangladesh

This paper uses two-dimensional boundary element method (BEM) numerical modeling to analyze the deformation and failure behavior of a coal seam and to understand the nature of gas flow into a roadway entering the Barapukuria coal mine in Bangladesh. The Barapukuria basin contains Permian-aged Gondwana coals with high volatile B bituminous rank. Three models (A, B, and C) are presented here. Model A assumes horseshoe-shaped geometry, model B assumes trapezoid-shaped geometry, and model C assumes horseshoe-shaped geometry coupled with a roof fall-induced cave generated by the break-up of rock materials along the vertical dimension of an igneous dyke. The simulation results show that there is little difference in strata deformation between models A and B. In model A, there is no horizontal tensional stress and the overall horizontal stress patterns are compressive, while the distribution and magnitude of vertical stress show higher tensional stresses on the immediate rib sides and floor. In model B, both horizontal and vertical stress distributions indicate low to medium tensional stresses on the immediate roof, floor, and rib sides, but compressive stresses are prominent toward the interior of the coal seam. Deformation vectors indicate that failure extends laterally to about 7.5 m around the excavation geometry.On the contrary, for model C, the distributions and magnitudes of horizontal and vertical stress show higher tensional stresses in both rib sides of the roof fall zone. The deformation around the dyke-induced perturbation zone affects a large volume of coal. The deformation vectors with high magnitudes are nearly horizontal and propagate laterally up to 30 m; whereas, low-magnitude deformation vectors extend about 25 m toward the roof and 20 m toward the floor. The vertical tensional displacement, which is concentrated in the floor and the left and right hand sides of the roof, propagates about 30 m on both sides and about 22 m in the floor. From these simulation results, it is thought that the extension of the dyke-induced perturbation zone toward the roof, floor, and rib sides of the entry roadway initially creates small tensional cracks that gradually grow into large-scale tensional features. These features could also be responsible for high concentrations of gas, which are emitted into the mine from fractured coals due to insufficient mine ventilation and low atmospheric pressure.     

The influence of wettability on NAPL dissolution fingering

Laboratory experiments and numerical simulations in homogeneous porous media were used to investigate the influence of porous medium wettability on the formation and growth of preferential dissolution pathways, dissolution fingers, during nonaqueous phase liquid (NAPL) dissolution. As the porous medium became increasingly NAPL-wet, dissolution fingers grew wider and slower. This result was observed in physical experiments with 0% and 100% NAPL-wet conditions and confirmed with numerical simulations at these and intermediate wettabilities. A previously derived expression for an upscaled mass transfer rate coefficient that accounts for the growth of dissolution fingers was used to quantify the effect of fingering on overall NAPL removal rates. For the test cases evaluated, NAPL dissolution fingering controlled the overall rate of NAPL dissolution after the dissolution front moved 4 cm in 0% NAPL-wet conditions and 18 cm in 100% NAPL-wet conditions. Thus, even in completely NAPL-wet media dissolution fingering may control the overall rate of NAPL dissolution after relatively short travel distances. The importance of NAPL dissolution fingering in heterogeneous systems with spatially varying NAPL saturations, though, remains an important question for future work.     

Flood susceptibility modelling using advanced ensemble machine learning models

Floods are one of nature's most destructive disasters because of the immense damage to land, buildings, and human fatalities.It is difficult to forecast the areas that are vulnerable to flash flooding due to the dynamic and complex nature of the flash floods.Therefore, earlier identification of flash flood susceptible sites can be performed using advanced machine learning models for managing flood disasters.In this study, we applied and assessed two new hybrid ensemble models, namely Dagging and Random Subspace(RS) coupled with Artificial Neural Network(ANN), Random Forest(RF), and Support Vector Machine(SVM) which are the other three state-of-the-art machine learning models for modelling flood susceptibility maps at the Teesta River basin, the northern region of Bangladesh.The application of these models includes twelve flood influencing factors with 413 current and former flooding points, which were transferred in a GIS environment.The information gain ratio, the multicollinearity diagnostics tests were employed to determine the association between the occurrences and flood influential factors.For the validation and the comparison of these models, for the ability to predict the statistical appraisal measures such as Freidman, Wilcoxon signed-rank, and t-paired tests and Receiver Operating Characteristic Curve(ROC) were employed.The value of the Area Under the Curve(AUC) of ROC was above 0.80 for all models.For flood susceptibility modelling, the Dagging model performs superior, followed by RF,the ANN, the SVM, and the RS, then the several benchmark models.The approach and solution-oriented outcomes outlined in this paper will assist state and local authorities as well as policy makers in reducing flood-related threats and will also assist in the implementation of effective mitigation strategies to mitigate future damage.     

Pan-African Magmatism, and Sedimentation in the NW Himalaya

Correlation of early Palaeozoic, Pan-African (500 ± 50 Ma) granites that intruded the Chail, Salkhala, Haimanta Formations in the Lesser Himalaya, Zanskar crystallines, and Lower Taglang La of Tso-Morari crystallines in the northwestern Himalaya, is based on the field relationship, tectonic setting, mineralogical, and geochemical characteristics, and isotope dating of the granites. These granite plutons exhibit identical petrographical, and geochemical character. The mineralogical composition of the granites is quite similar, consisting of quartz, K-feldspar, plagioclase feldspar, biotite, muscovite, garnet, tourmaline, ± cordierite, andalusite, and sillimanite fibrolite. The granite which are massive, and inequigranular in the core of the plutons, show strongly foliated character indicating development of ductile shear zone at the margins. These are peraluminous S-type granites having high A/CNK value (> 1). Presence of normative corundum, rounded shape of zircon, and high initial Sr ratio suggest crustal source of the granites. Mantle normalized spider-diagram exhibits similar characters for all these granitoids. The intrusion of the Pan-African granites mark an abrupt end of the sedimentation that continued virtually uninterrupted from Palaeoproterozoic. The sudden break in sedimentation towards the terminal phases of the Lower Cambrian has been observed in almost all parts in Lesser as well as the Tethys Himalaya. Occurrences of large number of plutons along different tectonic belts of northwestern Himalaya are indicative of widespread tectono-thermal event during early Palaeozoic (500 ± 50 Ma). The bracketing of the two features like, the break in sedimentation during post-Late Cambrian, and the intrusion of granites around 500 ± 50 Ma, is considered to be the result of a strong diastrophic orogenic event correlatable to the late phases of the Pan-African Orogeny in Africa.     

Mineralogical changes during intense chemical weathering of sedimentary rocks in Bangladesh

Bangladesh is situated in a subtropical to tropical climatic zone. A recently weathered crust has developed on sedimentary bedrock (sandstone, siltstone, shale and claystones) of Tertiary–Quaternary age. Weathered samples were collected from 16 sections totaling 68 samples and were analyzed mineralogically. The main primary minerals identified in the weathered crust of sedimentary rocks are quartz, plagioclase, K-feldspar, biotite, muscovite, sparse carbonate and epidote. The secondary minerals are kaolinite, illite, chlorite, gibbsite and goethite. Weathering initiated along the grain boundaries and cleavage planes of the minerals, forming small cloudy materials which were very difficult to identify. In the advanced stage of weathering, these cloudy materials have turned into secondary minerals. In region 1, high rain fall (7100 mm/yr) and monsoonic climate resulted in a kaolinite–gibbsite–goethite suite through the weathering of feldspars and biotite. The occurrence of gibbsite in the relatively elevated lands of Sylhet and Fe-kaolinite throughout the study areas is indicative of a humid–tropical climate during formation of the weathered crust.