Vertical geochemical variations and arsenic mobilization in the shallow alluvial aquifers of the Chapai-Nawabganj District, northwestern Bangladesh: implication of siderite precipitation

Core sediments from three disturbed boreholes (JOR, GHAT, and RAJ) and two undisturbed boreholes (DW1 and DW2) were collected in the study area of the Chapai-Nawabganj district of northwestern Bangladesh for geochemical analyses. In the study area, groundwater samples from fourteen As-contained private wells and five nested piezometers at both the DW1 and DW2 boreholes were also collected and analyzed. The groundwater arsenic concentrations in the uppermost aquifer (10–40 m of depth) range from 3 to 315 μg/L (mean 47.73 ± 73.41 μg/L), while the arsenic content in sediments range from 2 to 14 mg/kg (mean 4.36 ± 3.34 mg/kg). An environmental scanning electron microscope (ESEM) with an energy dispersive X-ray spectrometer was used to investigate the presence of major and trace elements in the sediments. Groundwaters in the study area are generally the Ca–HCO₃ type with high concentrations of As, but low levels of Fe, Mn, NO₃ ^? and SO ₄ ^?2 . The concentrations of As, Fe, Mn decrease with depth in the groundwater, showing vertical geochemical variations in the study area. Statistical analysis clearly shows that As is closely associated with Fe and Mn in the sediments of the JOR core (r = 0.87, p < 0.05 for Fe and r = 0.78, p < 0.05 for Mn) and GHAT core (r = 0.95, p < 0.05 for Fe and r = 0.93, p < 0.05 for Mn), while As is not correlated with Fe and Mn in groundwater. The comparatively low Fe and Mn concentrations in some groundwater and the ESEM image revealed that siderite precipitated as a secondary mineral on the surface of the sediment particles. The correlations along with results of sequential extraction experiments indicated that reductive dissolution of FeOOH and MnOOH represents a mechanism for releasing arsenic into the groundwater.     

Generating non-stationary random fields of auto-correlated,normally distributed CPT profile by matrix decomposition method

This paper emphasises the true realisation of Cone Penetration Test (CPT) profiles considering non-stationary nature of the data. Formulation of stationary random field theory has been modified and adapted to non-stationary state in order to take into account the mean and variance variability for soil properties. Multi-variance correlation matrix along with the Cholesky decomposition technique was employed to produce realisations of non-homogenous and non-stationary random fields of CPT profiles. A piecewise and segmental data realisation according to the lithology and site class specifications acquired directly from CPT data is adopted in this study so as to render an accurate data simulation. For validation of proposed method 8 CPT test profiles collected from Urmia Lake site have been introduced and simulated by the stationary and non-stationary algorithms. The mean correlation coefficient between the actual CPT data profiles and related realisations along with some other important statistical parameters and their coefficients of variation strongly demonstrate that non-stationary random field generation technique gives quite better accuracy, by comparison to the conventional stationary random field generation scheme.     

Probabilistic rainfall thresholds for flood forecasting: evaluating different methodologies for modelling rainfall spatial correlation (or dependence)

Rainfall threshold (RT) method is one of the evolving flood forecasting approaches. When the cumulative rainfall depth for a given initial soil moisture condition intersects the threshold rainfall curve, the peak discharge is expected to be equal or greater than the threshold discharge for flooding at the target site. Besides the total rainfall depth, spatial and temporal distribution of rainfall impacts the flood peak discharge and the time to peak. To revisit a previous study conducted by the authors, in which spatially independent rainfall pattern was assumed, the spatial distribution of rainfall was simulated following a Monte Carlo approach. The structure of the spatial dependence among sub‐watersheds' rainfalls was taken into account under three different scenarios, namely independent, bivariate copula (2copula) and multivariate Gaussian copula (MGC). For each set of generated random dimensionless rainfalls, the probabilistic RT curves were derived for dry moisture condition. Results were evaluated with both historical and simulated events. For the simulated events, threshold curves were assessed by means of categorical statistics, such as hit rate, false rate and critical success index (CSI). Results revealed that the best performance based on the CSI criterion corresponded to 50% curve in 2copula and MGC scenarios as well as 90% curve in the independent scenario. The recognition of 50% curve in 2copula and MGC scenarios is in agreement with our expectations that the mean probable curve should have the best performance. Moreover, the proposed inclusion of spatially dependent rainfall scenario improved the performance of RT curves by about 25% in comparison with the presumed spatially uniform rainfall scenario. Copyright © 2011 John Wiley & Sons, Ltd.     

Steady critical discharge rates from vertical and horizontal wells in fresh–saline aquifers with sharp interfaces

In this study, the critical (or maximum) discharge rates before saline water enters a well were determined for vertical and horizontal wells in a freshwater aquifer which is separated from a static saline aquifer by a sharp interface. Flow around the well was solved by integration of a point sink solution along the well axis, and both the critical discharge rate and critical interface rise were determined through a comparison of the heads and vertical gradients at the saline–fresh water interface. The rates were determined for vertical and horizontal wells with various lengths and depths for different aquifer salinities. Results were generalized by drawing dimensionless type curves. The results showed that the dimensionless total critical discharge rates are higher for the longer horizontal wells and longer vertical wells with a certain bottom depth, and they almost linearly decrease with well depth at rates of 0.7–0.9. For the dimensionless well length of 0.2, the dimensionless total discharge rate of a horizontal well is about 0.1 more than that of a vertical well with the same length and well-top depth. Also, the critical discharge rates per unit length of well are inversely proportional to well length and remarkably higher for shallower wells. Additionally, the critical pumping rate is proportional to the salinity difference of the aquifers. These results were confirmed by comparison to existing solutions for vertical wells with dimensionless lengths of 0.2, 0.5 and 0.6, and for critical interface rises in the range of 0.75–0.9.     

Synoptic analysis of dust storms in the Middle East

Dust storm in the Middle East and south-west Asia is a natural hazard and the Tigris-Euphrates alluvial plain has been recognized as the main dust source in this area. In this study, more than 60 dust storms that occurred during the period 2003–2011 are investigated on the basis of MODIS satellite images, and 12 of the dust storms are selected for synoptic analysis using the NCEP-NCAR Reanalysis Data. The potential dust sources in the Middle East and south-west Asian region (20°E to 80°E, 5°N to 50°N) are analyzed and used in the synoptic analysis. Dust storms in the region can be grouped into two main categories, i.e., the Shamal dust storms and the frontal dust storms. Synoptic systems, associated with the two categories, are distinguished and the frequency of the patterns is identified. For 68% of the Shamal dust storms, a high pressure system is situated between 0°E to 30°E and 27°N to 45°N, and a low pressure system between 50°E to 70°E and 23°N to 43°N. For 86% of the frontal dust storms, a high is located between 51°E to 67°E and 18°N to 33°N and a low between 28°E to 48°E and 32°N to 43°N. Three main patterns for Shamal dust storms are identified, which represent about 60% of the Shamal dust storms. This analysis confirms that the Shamal is related to the anticyclones located over northern Africa to Eastern Europe and the monsoon trough over Iraq, southern Iran, Pakistan and the Indian Subcontinent. The analysis also shows that the main dust sink for the frontal dust storms in Tigris and Euphrates alluvial plain extends from center of Iraq to west and center of Iran and, in most severe cases, to northern Iran and the southern coast of the Caspian Sea.     

The effect of specimen diameter size on uniaxial compressive strength,P-wave velocity and the correlation between them

Estimation of uniaxial compressive strength (UCS) by P-wave velocity (V_P) is of great interest to geotechnical engineers in various design projects. The specimen diameter size is one of the main factors that influence rock parameters such as UCS and V_P. In this study, the diameter size of specimens that effect UCS and V_P is investigated. Moreover, the correlation between UCS and V_P are examined via empirical analysis. For this purpose, 15 travertine samples were collected and core specimens with a diameters size of 38, 44, 54, 64 and 74 mm were prepared. Then, uniaxial compressive strength and P-wave velocity tests were conducted according to the procedure suggested by ISRM (1981). It is concluded that the diameter size of the specimen has a significant effect on UCS and V_P. Moreover, it was found that the best correlation between relevant parameters obtained for the specimen diameter of 38 mm.     

An investigation of trends in precipitation volume for the last three decades in different regions of Fars province, Iran

Under condition of climate changes as global warming, monitoring and detecting trend of precipitation volume is essential and will be useful for agricultural sections. Considering the fact that there were not enough research related to precipitation volume, this study aimed to determine trends in precipitation volume, monthly and annually in different regions of Fars province for the last three decades (33?years period; 1978–2010). Fars province is located in arid and semi-arid regions of Iran, and it plays an important role in agricultural production. Inverse distance weighting interpolation method was used to provide precipitation data for all regions. To analyze the trends of precipitation volume, Mann–Kendall test, Sen’s slope estimator, and 10-year moving average low-pass filter (within time series) were used. The negative trends were identified by the Sen’s slope estimator as well as Mann–Kendall test. However, all the trends were insignificant at the surveyed confidence level (95%). With regards to the application of 10-year moving average low-pass filter, a considerable decreasing trend was observed after around year 1994. Since one of the most important restrictions in agricultural development of the Fars province is lack of sufficient water resources, any changes onward to lack of sufficient precipitation impose impressive pressure and stress on valuable resources and subsequently agricultural production.