Element mobility during pyrite weathering: implications for acid and heavy metal pollution at mining-impacted sites

Based on back scattered electron images and electron micro-probe analysis results, four alteration layers, including a transition layer, a reticulated ferric oxide layer, a nubby ferric oxide layer and a cellular ferric oxide layer, were identified in the naturally weathering products of pyrite. These layers represent a progressive alteration sequence of pyrite under weathering conditions. The cellular ferric oxide layer correlates with the strongest weathering phase and results from the dissolution of nubby ferric oxide by acidic porewater. Leaching coefficient was introduced to better express the response of element mobility to the degree of pyrite weathering. Its variation shows that the mobility of S, Co and Bi is stronger than As, Cu and Zn. Sulfur in pyrite is oxidized to sulfuric acid and sulfate that are basically released into to porewater, and heavy metals Co and Bi are evidently released by acid dissolution. As, Cu and Zn are enriched in ferric oxide by adsorption and by co-precipitation, but they would re-release to the environment via desorption or dissolution when porewater pH becomes low enough. Consequently, Co, Bi, As, Cu and Zn may pose a substantial impact on water quality. Considering that metal mobility and its concentration in mine waste are two important factors influencing heavy metal pollution at mining-impacted sites, Bi and Co are more important pollutants in this case.     

Hydrochemical characteristics and salinity of groundwater in the Ejina Basin, Northwestern China

A hydrochemical investigation was conducted in the Ejina Basin to identify the hydrochemical characteristics and the salinity of groundwater. The results indicate that groundwater in the area is brackish and are significantly zonation in salinity and water types from the recharge area to the discharge area. The ionic ration plot and saturation index (SI) calculation suggest that the silicate rock weathering and evaporation deposition are the dominant processes that determine the major ionic composition in the study area. Most of the stable isotope δ¹⁸O and δD compositions in the groundwater is a meteoric water feature, indicating that the groundwater mainly sources from meteoric water and most groundwater undergoes a long history of evaporation. Based on radioactive isotope tritium (³H) analysis, the groundwater ages were approximately estimated in different aquifers. The groundwater age ranges from less than 5 years, between 5 years and 50 years, and more than 50 years. Within 1 km of the river water influence zone, the groundwater recharges from recent Heihe river water and the groundwater age is about less than 5 years in shallow aquifer. From 1 km to 10 km of the river water influence zone, the groundwater sources from the mixture waters and the groundwater age is between 5 years and 50 years in shallow aquifer. The groundwater age is more than 50 years in deep confined aquifer.     

青藏高原腹地植物碳同位素组成对环境条件的响应

现代植物碳同位素组成是特定环境影响的结果,通过对植物碳同位素组成的研究可以揭示植物生长期环境信息。针对青藏高原腹地高寒草甸~高寒草原过渡区植被碳同位素组成进行研究;该区高山嵩草样δ13C值在-25.63‰~-27.95‰间,平均值-26.63‰;高寒草原区混合样δ13C值于-26.29‰~-27.73‰间,平均值-27.04‰。高山嵩草样δ13C值总体呈现由南东往北西方向正偏趋势,研究区北部高寒草原区混合植物样也呈现出由南向北富重碳同位素趋势。这些变化规律被认为是主要受降水环境影响的结果,而区域内降水条件的展布规律则是受高原夏季风运移方式的控制。对植物δ13C值与地理位置的回归分析表明,该区植被碳同位素组成与地理位置相关,高山嵩草样(r=0.44603,n=29,p<0.05)和混合样(r=0.8112,n=5,p<0.1)均表现出对区域降水环境条件的良好响应。据此,以该区植物δ13C值为背景,进行合理推算,拟定了研究区内干旱区和湿润区界限的位置。     

Mixing of event and pre‐event water in a shallow Entisol in sloping farmland based on isotopic and hydrometric measurements,SW China

Water percolation and flow processes in subsurface geologic media play an important role in determining the water source for plants and the transport of contaminants or nutrients, which is essential for water resource management and the development of measures for pollution mitigation. During June 2013, the dynamics of the rainwater, soil water, subsurface flows and groundwater in a shallow Entisol on sloping farmland were monitored using a hydrometric and isotopic approach. The results showed that effective mixing of rainwater and soil water occurred in hours. The rebound phenomenon of δD profiles in soils showed that most isotope‐depleted rainwater largely bypassed the soil matrix when the water saturation in the soil was high. Preferential‐flow, which was the dominant water movement pattern in the vadose zone, occurred through the whole soil profile, and infrequent piston‐flow was mainly found at 20–40 cm in depth. The interflow in the soil layer, composed of 75.2% rainwater, was only generated when the soil profile had been saturated. Underflow in the fractured mudrock was the dominant flow type in this hillslope, and outflow was dominated by base flow (groundwater flow) with a mean contribution of 76.7%. The generation mechanism of underflow was groundwater ridging, which was superimposed upon preferential‐flow composed mainly of rainwater. The quick mixing process of rainwater and soil water and the rapid movement of the mixture through preferential channels in the study soil, which shows a typical bimodal pore size distribution, can explain the prompt release of pre‐event water in subsurface flow. Water sources of subsurface flows at peak discharge could be affected by the antecedent soil water content, rain characteristics and antecedent groundwater levels. Copyright © 2016 John Wiley & Sons, Ltd.     

Assessing storm surge impacts on coastal inundation due to climate change: case studies of Baltimore and Dorchester County in Maryland

Natural Hazards - Hurricane Isabel (2003) generated record flooding around Chesapeake Bay and caused extensive damage in rural Eastern Shore of Maryland and metropolitan cities like Baltimore....     

Implicit integration under mixed controls of a breakage model for unsaturated crushable soils

This paper discusses a series of stress point algorithms for a breakage model for unsaturated granular soils. Such model is characterized by highly nonlinear coupling terms introduced by breakage‐dependent hydro‐mechanical energy potentials. To integrate accurately and efficiently its constitutive equations, specific algorithms have been formulated using a backward Euler scheme. In particular, because implementation and verification of unsaturated soil models often require the use of mixed controls, the incorporation of various hydro‐mechanical conditions has been tackled. First, it is shown that the degree of saturation can be replaced with suction in the constitutive equations through a partial Legendre transformation of the energy potentials, thus changing the thermomechanical state variables and enabling a straightforward implementation of a different control mode. Then, to accommodate more complex control scenarios without redefining the energy potentials, a hybrid strategy has been used, combining the return mapping scheme with linearized constraints. It is shown that this linearization strategy guarantees similar levels of accuracy compared with a conventional strain–suction‐controlled implicit integration. In addition, it is shown that the use of linearized constraints offers the possibility to use the same framework to integrate a variety of control conditions (e.g., net stress and/or water‐content control). The convergence profiles indicate that both schemes preserve the advantages of implicit integration, that is, asymptotic quadratic convergence and unconditional stability. Finally, the performance of the two implicit schemes has been compared with that of an explicit algorithm with automatic sub‐stepping and error control, showing that for the selected breakage model, implicit integration leads to a significant reduction of the computational cost. Such features support the use of the proposed hybrid scheme also in other modeling contexts, especially when strongly nonlinear models have to be implemented and/or validated by using non‐standard hydro‐mechanical control conditions. Copyright © 2015 John Wiley & Sons, Ltd.     

Development of a soil moisture‐based distributed hydrologic model for determining hydrologically based critical source areas

A simple grid cell‐based distributed hydrologic model was developed to provide spatial information on hydrologic components for determining hydrologically based critical source areas. The model represents the critical process (soil moisture variation) to run‐off generation accounting for both local and global water balance. In this way, it simulates both infiltration excess run‐off and saturation excess run‐off. The model was tested by multisite and multivariable evaluation on the 50‐km² Little River Experimental Watershed I in Georgia, U.S. and 2 smaller nested subwatersheds. Water balance, hydrograph, and soil moisture were simulated and compared to observed data. For streamflow calibration, the daily Nash‐Sutcliffe coefficient was 0.78 at the watershed outlet and 0.56 and 0.75 at the 2 nested subwatersheds. For the validation period, the Nash‐Sutcliffe coefficients were 0.79 at the watershed outlet and 0.85 and 0.83 at the 2 subwatersheds. The per cent bias was less than 15% for all sites. For soil moisture, the model also predicted the rising and declining trends at 4 of the 5 measurement sites. The spatial distribution of surface run‐off simulated by the model was mainly controlled by local characteristics (precipitation, soil properties, and land cover) on dry days and by global watershed characteristics (relative position within the watershed and hydrologic connectivity) on wet days when saturation excess run‐off was simulated. The spatial details of run‐off generation and travel time along flow paths provided by the model are helpful for watershed managers to further identify critical source areas of non‐point source pollution and develop best management practices.