Modelling the internal stress field in argillaceous rocks under humidification/desiccation

This paper deals with analytical and numerical modelling of the internal stress generated in argillaceous rocks during humidification/desiccation processes, which is an essential issue for damage study. This local stress field arises from two mechanisms: (i) complex interactions between free swelling/shrinking clay matrix and non‐strained inclusions of carbonate and quartz and (ii) a self‐restraint effect induced by the moisture gradient during the transient moisture exchange process. The inclusion–matrix interaction is investigated in different cases. Firstly, the analytical solution of the stress around a cylindrical inclusion embedded in an infinite swelling matrix is derived: The inclusion would suffer tension (compression) under humidification (desiccation), and the resulting cracking patterns are discussed. Then, the problem of two inclusions with different distances in an infinite swelling matrix is considered, and it is shown that the local stress around an inclusion will be perturbed and amplified by neighbouring inclusions. Finally, an inclusion outcropping at the free surface of a swelling matrix is modelled as to investigate the effect of free surface: The inclusion–matrix interface undergoes shear stresses of which the maximum is found at the free surface. In addition to the inclusion–matrix interaction, the self‐restraint effect is investigated: The induced stress is maximal at the beginning of humidification/desiccation processes and vanishes gradually with time. The quantity of the self‐restraint stress is strongly controlled by the hydric loading rate. Copyright © 2014 John Wiley & Sons, Ltd.     

Deep seismic reflection study over the Vindhyans of Rajasthan: Implications for geophysical setting of the basin

This paper presents results of high-resolution deep seismic reflection profiling of the Proterozoic Vindhyan basin of the Rajasthan area along the Chandli-Bundi-Kota-Kunjer profile. Seismic images have been used to estimate the thickness of Vindhyan strata as well as to understand the tectonic framework of the basin. The results are constrained by gravity, magnetic and magnetotelluric data. The study reveals gentle SE-dipping reflection bands representing the Vindhyan strata. The seismic sections depict gradual thickening of the Vindhyan succession towards southeast from Bundi. The velocities of the upper and lower Vindhyans are identified as 4.6-4.8 km/s and 5.1-5.3 km/s. The NW limit of the Vindhyan basin is demarcated by the Great Boundary Fault (GBF) that manifests as a 30 km wide NW dipping thrust fault extending to a depth of 30 km.     

A subduction wedge origin for Paleoarchean peridotitic diamonds and harzburgites from the Panda kimberlite, Slave craton: evidence from Re–Os isotope systematics

An extensive study of peridotitic sulfide inclusion bearing diamonds and their prospective harzburgitic host rocks from the 53 Ma Panda kimberlite pipe, Ekati Mine, NWT Canada, has been undertaken with the Re–Os system to establish their age and petrogenesis. Diamonds with peridotitic sulfide inclusions have poorly aggregated nitrogen (<30% N as B centers) at N contents of 200–800 ppm which differs from that of chromite and silicate bearing diamonds and indicates residence in the cooler portion of the Slave craton lithospheric mantle. For most of the sulfide inclusions, relatively low Re contents (average 0.457 ppm) and high Os contents (average 339 ppm) lead to extremely low ¹⁸⁷Re/¹⁸⁸Os, typically << 0.05. An age of 3.52 ± 0.17 Ga (MSWD = 0.46) and a precise initial ¹⁸⁷Os/¹⁸⁸Os of 0.1093 ± 0.0001 are given by a single regression of 11 inclusions from five diamonds that individually provide coincident internal isochrons. This initial Os isotopic composition is 6% enriched in ¹⁸⁷Os over 3.5 Ga chondritic or primitive mantle. Sulfide inclusions with less radiogenic initial Os isotopic compositions reflect isotopic heterogeneity in diamond forming fluids. The harzburgites have even lower initial ¹⁸⁷Os/¹⁸⁸Os than the sulfide inclusions, some approaching the isotopic composition of 3.5 Ga chondritic mantle. In several cases isotopically distinct sulfides occur in different growth zones of the same diamond. This supports a model where C–O–H–S fluids carrying a radiogenic Os signature were introduced into depleted harzburgite and produced diamonds containing sulfides conforming to the 3.5 Ga isochron. Reaction of this fluid with harzburgite led to diamonds with less radiogenic inclusions while elevating the Os isotope ratios of some harzburgites. Subduction is a viable way of introducing such fluids. This implies a role for subduction in creating early continental nuclei at 3.5 Ga and generating peridotitic diamonds.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Holistic assessment of groundwater resources and regional environmental problems in the North China Plain

Water balance components of the North China Plain (NCP) were analyzed, indicating the decrease both in precipitation and evaporation. The decreased precipitation and expansion of water use for agriculture, industrial and domestic purposes have caused a water crisis, which was managed until now by diverting water from the Yellow River and over exploitation of groundwater. The groundwater resource was assessed by estimating its recharge in both upper unconfined and lower confined layers, yielding a total value of 1.65 × 10¹⁰ m³/a. Total groundwater use was estimated and judged by the actual water table drawdown. Salt accumulation, water table decrease, fluoride and nitrate pollution were all found to be major regional environmental problems. Furthermore, heavy metals were found in high content in the soil and surface water in suburbs of large cities, posing a potential risk of pollution in the groundwater. It has been verified by isotropic data that dry conditions have occurred since 10 ka and are therefore part of the natural process. Possible solutions for water crises in the NCP are proposed.     

Three-dimensional supervirtual seismic refraction interferometry

Improved supervirtual interferometry (SVI) method that enhance the signal-to-noise ratio (SNR) of noisy seismic refraction data is presented. This 3D supervirtual algorithm is helpful in removal of side lobes caused by convolution in the conventional 3D SVI method. In this study, we extend the supervirtual seismic refraction interferometry method to 3D geometries commonly used in active seismic exploration. To achieve this objective, synthetic 3D seismic refraction data were created using single patch orthogonal geometry. Simple two-layers velocity model was used. A zero phase wavelet was convolved with the refraction travel time to create the seismic wiggle trace. The refraction arrival is computed using the equation for horizontal refractor case and is not adapted for dipping refractor. Noise of Gaussian distribution with zero mean and 0.25 standard deviation was added to simulate a case of moderate ambient noise. The supervirtual algorithm, consisting of crosscorrelation, alignment, summation, and first arrival calculation is performed. By aligning and summation of all the correlogram, the stationary position of source-receiver pairs as required in the 2D supervirtual method are eliminated in this case. Synthetic data presented in this study shows accurate first arrivals after the application of the 3D SVI and traces with much better SNR than the actual traces.     

A spatio-temporal three-dimensional conceptualization and simulation of Dera Ismail Khan alluvial aquifer in visual MODFLOW: a case study from Pakistan

Dera Ismail Khan (DIK) is situated in the Lower Indus Basin of Pakistan. The land use has been changed in the canal command area due to irrigation activities near the Indus River. To check the current status and predict the groundwater levels in the area, the unconfined aquifer has been simulated in Visual MODFLOW for a period of 35 years, i.e., from 1985 to 2020. The 2900-km² area has been modeled with a grid of 500 by 500 m and the depth set to 100 m. The aquifer in the study area has been divided vertically and laterally into three and ten zones, respectively, for the characterization. Water wells and streams were used as the sinks and hydrologic boundaries, respectively. The model was successfully calibrated in steady and the non-steady state. The simulation revealed that the whole simulation can be divided into two phases, i.e., before and after the construction of the Chashma Right Bank Canal (CRBC), whereas the results were summarized in the form of water table depth maps and groundwater budget calculations. To determine the groundwater sustainability, a conjunctive use scenario has been employed to simulate the aquifer dynamics till 2020. The simulation revealed incremental drawdowns till the end.     

Effects of biochar and bio-fertilizer on yield and qualitative properties of soybean and some chemical properties of soil

This study attempted to examine the effects of biochar amendment together with bio-fertilizer on soybean yield and its qualitative properties, as well as a few chemical properties of soil through a factorial randomized complete block design at three replications in east of Golestan Province (Iran) during 2014. The two factors under study included the following: (1) biochar amendment (in four levels of 0, 2.5, 8, and 16 tons per hectare), (2) bio-fertilizer containing phosphorus and sulfur growth-promoting rhizobacteria (in two levels of inoculation and non-inoculation) applied through foliar feeding. The results of analysis of variance indicated that interactions of biochar amendment and bio-fertilizer on harvest index and grain yield were significant (p ≤?0.01). According to the results of this study, the highest harvest index and oil content were 56.9, and 17.7%, respectively, in the treatment of 8 tons per hectare biochar and inoculation with bio-fertilizer. The lowest harvest index and the lowest oil content were in the control treatment. The interaction of biochar and bio-fertilizer on bulk density and cation exchange capacity was significant (p ≤?0.01). The results of this study demonstrated that biochar affected the amount of residual nitrogen in the soil after harvest, cation exchange capacity (CEC), acidity (pH), and electrical conductivity (EC). The highest grain yield (3440 kg/ha) was in the 8-ton biochar treatment with inoculated bio-fertilizer. Our study concludes that the biochar and bio-fertilizers can improve grain yield of soybean till 51% relative to the control.     

Evidence for the occurrence of hydrogeochemical processes in the groundwater of Khoy plain,northwestern Iran,using ionic ratios and geochemical modeling

Rapid population growth, industrialization, and agricultural expansion in the Khoy area (northwestern Iran) have led to its dependence on groundwater and degradation of groundwater quality. This study attempts to decipher the major processes and factors that degrade the groundwater quality of the Khoy plain. For this purpose, 54 groundwater samples from unconfined and confined aquifers of the plain were collected in July 2017 and analyzed for major cations and anions (Na, K, Ca, Mg, HCO₃, SO_4, and Cl), minor ions (NO₃ and F), and Al. Magnesium and bicarbonate were identified as the dominant cation and anion, respectively. Several ionic ratios and geochemical modeling using PHREEQC indicated that the most important hydrogeochemical processes to affect groundwater quality in the plain were weathering and dissolution of evaporitic and silicate minerals, mixing, and ion exchange. There were smaller effects from evaporation and anthropogenic factors (e.g., industries). Results showed that the high salinity of the groundwater in the northeast area of the plain was due to the high solubility of the evaporitic minerals, e.g., halite and gypsum. Reverse ion exchange and the contribution of mineral dissolution were more significant than ion exchange in the northeastern part of the plain. Elevated salinity of the groundwater in the southeast was attributed mostly to reverse ion exchange and somewhat to evaporation.     

Hydrogeochemical processes of thallium in a rural area of Southwest Guizhou, China

The water system in a rural area of Lanmuchang in Southwest Guizhou is facing a risk of thallium (Tl) contamination due to Tl mineralization around the area. The major trace elements and Tl in the water system are studied to understand the hydrogeochemical processes of Tl constrained by Tl mineralization. The results showed that the dispersion pattern of Tl follows a descending order in concentration from mine groundwater (deep groundwater) →stream water→shallow groundwater→background water, reflecting the impact of Tl mineralization on the hydrogeochemical composition. Tl concentrations in stream water in both regimes are remarkably higher (2-30 fold) downstream than up- and mid-streams, probably caused by the unidentified discharge of deep groundwater. Low Tl levels are detected in the current drinking water, however, the highly elevated Tl in stream water and ground water may pose a potential environmental risk through daily washing and agricultural irrigation. This study suggests that human activities, such as agricultural irrigation, could intensify the environmental risk of Tl.