No abstract

Abstract. The aim of the present study is to evaluate the potential ecological risk and trend of soil heavy-metal pollution around a coal gangue dump in Jilin Province (Northeast China). The concentrations of Cd, Pb, Cu, Cr and Zn were monitored by inductively coupled plasma mass spectrometry (ICP-MS). The potential ecological risk index method developed by Hakanson (1980) was employed to assess the potential risk of heavy-metal pollution. The potential ecological risk in the order of ER(Cd) > ER(Pb) > ER(Cu) > ER(Cr) > ER(Zn) have been obtained, which showed that Cd was the most important factor leading to risk. Based on the Cd pollution history, the cumulative acceleration and cumulative rate of Cd were estimated, then the fixed number of years exceeding the standard prediction model was established, which was used to predict the pollution trend of Cd under the accelerated accumulation mode and the uniform mode. Pearson correlation analysis and correspondence analysis are employed to identify the sources of heavy metals and the relationship between sampling points and variables. These findings provided some useful insights for making appropriate management strategies to prevent or decrease heavy-metal pollution around a coal gangue dump in the Yangcaogou coal mine and other similar areas elsewhere.

[1] Three conceptual models are evaluated for estimating transmissivity (T) fields using data from sequential pumping tests at a field site and data from similar tests simulated in a synthetic aquifer. The three approaches are (1) an equivalent homogeneous approach, (2) a heterogeneous approach based on a single pumping test, and (3) a heterogeneous approach based on joint interpretation of the sequential pumping tests (i.e., hydraulic tomography, HT). They are evaluated on the basis of their abilities to obtain representative estimates of the T field of the aquifer and, more importantly, on the ability of their estimates to predict drawdown distributions in the aquifer induced by independent validation pumping tests. Results show that the first approach yields scenario-dependent T estimates, which vary with the location of the pumping well. Independent validation tests show that the predicted drawdowns in both aquifers are biased and dispersed. While the second approach produces scenario-dependent T spatial distributions capturing the general pattern of the aquifer, the T fields consistently yield better drawdown predictions than those based on the first approach. Lastly, the joint interpretation approach reduces the scenario dependence of the T estimates and improves the quality of the T estimates as more data sets from sequential pumping tests are included. More importantly, the resultant T estimates lead to the best prediction of different flow events. The robustness of the joint interpretation is then elucidated.

[1] This paper investigates optimal sampling times of drawdowns for the analysis of hydraulic tomography (HT) survey. The investigation was carried out by analyzing the spatial and temporal evolution of cross-correlations between the head responses at an observation well and transmissivity (T) and storage coefficient (S) properties during a pumping test in homogeneous and heterogeneous aquifers. The analysis shows that the cross-correlation between the head and S values is limited to the region between the observation and the pumping well in the aquifers: It reaches the highest value near the early time (t m ), and decays to zero afterwards. The time t m is approximately equal to the time t 0 at which the extrapolated drawdown from the first straight line portion of an observed drawdown-log time plot becomes zero. At early times, the high cross-correlation between the head and T is confined to the region between the observation and the pumping well. This region then evolves into two humps: One on each side of the circular region encompassing the observation well and the pumping well. The size of the two humps expands and their values reach the maximum as flow reach steady-state. As a consequence, we hypothesize that pairs of head data at t 0 and those at either the steady-state or a late time during an HT survey could yield the best estimates of the heterogeneous T and S fields. Results from numerical experiments have verified this hypothesis and demonstrated that this sampling strategy is generally applicable even when the boundary condition is unknown. We, therefore, recommend in principle that (1) carrying out pumping tests of HT surveys for sufficiently long period of time such that drawdown reaches the entire area of interest and (2) using a constant head or zero drawdown for all boundaries during the inverse modeling analysis.

[1] In this study, the S-shaped log-log drawdown-time curve typical of pumping tests in unconfined aquifers is reinvestigated via numerical experiments. Like previous investigations, this study attributes the departure of the S shape from the drawdown-time behavior of the confined aquifer to the presence of an "additional" source of water. Unlike previous studies, this source of water is reinvestigated by examining the temporal and spatial evolution of the rate of change in storage in an unconfined aquifer during pumping. This evolution is then related to the transition of water release mechanisms from the expansion of water and compaction of the porous medium to the drainage of water from the unsaturated zone above the initial water table and initially saturated pores as the water table falls during the pumping of the aquifer. Afterward, the 1-D vertical drainage process in a soil column is simulated. Results of the simulation show that the transition of the water release mechanisms in the 1-D vertical flow without an initial unsaturated zone can also yield the S-shaped drawdown-time curve as in an unconfined aquifer. We therefore conclude that the transition of the water release mechanisms and vertical flow in the aquifer are the cause of the S-shaped drawdown-time curve observed during pumping in an unconfined aquifer. We also find that the moisture retention characteristics of the aquifer material have greater impact than its relative permeability characteristics on the drawdown-time curve. Furthermore, influences of the spatial variability of saturated hydraulic conductivity, specific storage, and saturated moisture content on the drawdown curve in the saturated zone are found to be more significant than those of other unsaturated properties. Finally, a cross-correlation analysis reveals that the drawdown at a location in a heterogeneous unconfined aquifer is mainly affected by local heterogeneity near the pumping and observation wells. Applications of a model assuming homogeneity to the estimation of aquifer parameters as such may require a large number of observation wells to obtain representative parameter values. In conclusion, we advocate that the governing equation for variably saturated flow through heterogeneous media is a more appropriate and realistic model that explains the S-shaped drawdown-time curves observed in the field.