Hydrogeochemical Indicators of Groundwater Flow Systems in the Yangwu River Alluvial Fan,Xinzhou Basin,Shanxi, China

Based on analysis of groundwater hydrochemical and isotopic indicators, this article aims to identify the groundwater flow systems in the Yangwu River alluvial fan, in the Xinzhou Basin, China. Groundwater δ²H and δ¹⁸O values indicate that the origin of groundwater is mainly from precipitation, with local evaporative influence. d-excess values lower than 10% in most groundwaters suggest a cold climate during recharge in the area. Major ion chemistry, including rCa/rMg and rNa/rCl ratios, show that groundwater salinization is probably dominated by water–rock interaction (e.g., silicate mineral weathering, dissolution of calcite and dolomite and cation exchange) in the Yangwu River alluvial fan, and locally by intensive evapotranspiration in the Hutuo River valley. Cl and Sr concentrations follow an increasing trend in shallow groundwater affected by evaporation, and a decreasing trend in deep groundwater. ⁸⁷Sr/⁸⁶Sr ratios reflect the variety of lithologies encountered during throughflow. The groundwater flow systems (GFS) of the Yangwu River alluvial fan include local and intermediate flow systems. Hydrogeochemical modeling results, simulated using PHREEQC, reveal water–rock interaction processes along different flow paths. This modeling method is more effective for characterizing flow paths in the intermediate system than in the local system. Artificial exploitation on groundwater in the alluvial fan enhances mixing between different groundwater flow systems.     

VIC模型与SWAT模型在中小流域径流模拟中的对比研究

结合GIS与RS技术的分布式水文模型已成为当今水文界研究的重点。从气象与水文水资源学科交叉的角度对分布式水文模型VIC模型与SWAT模型进行研究,并将其应用于白莲河流域,以此探讨该模型在中小流域的适用性。模拟结果表明,VIC模型与SWAT模型在白莲河流域率定期与检验期的模拟效果相差很小。SWAT模型的效率系数与相关系数略高一些,SWAT模型的模拟效果比较平均,每年相差不大;但VIC模型在1995年和1999年模拟效果明显好于其它年份,尤其在2002年,VIC模型模拟的洪峰与实测的相差较大,从而影响总体的效率系数偏低,而SWAT模型模拟的更接近实测值。研究结果表明两种模型对于我国中小流域的径流模拟具有一定的适用性.     

On causes of gigantism in herbaceous plants

Landscape-geochemical features providing for manifestation of gigantism in herbaceous plants have been revealed in natural habitats in the south of Sakhalin and Kunashir islands. Tall herb assemblages have proved to be associated with geochemical landscapes characterized by reducing (gley or hydrogen sulfide) conditions and increased contents of petroleum hydrocarbons and some trace elements (total and movable forms). A hypothesis is put forward that gigantism in herbaceous plants is manifested in zones of active faults, which serve as a kind of conduits supplying endogenous heat, matter, and water to the root systems.     

Numerical analysis of water and solute transport in variably-saturated fractured clayey till

This study numerically investigates the influence of initial water content and rain intensities on the preferential migration of two fluorescent tracers, Acid Yellow 7 (AY7) and Sulforhodamine B (SB), through variably-saturated fractured clayey till. The simulations are based on the numerical model HydroGeoSphere, which solves 3D variably-saturated flow and solute transport in discretely-fractured porous media. Using detailed knowledge of the matrix, fracture, and biopore properties, the numerical model is calibrated and validated against experimental high-resolution tracer images/data collected under dry and wet soil conditions and for three different rain events. The model could reproduce reasonably well the observed preferential migration of AY7 and SB through the fractured till, although it did not capture the exact depth of migration and the negligible impact of the dead-end biopores in a near-saturated matrix. A sensitivity analysis suggests fast flow mechanisms and dynamic surface coating in the biopores, and the presence of a plough pan in the till.     

Converging hydrostatic and hydromechanic concepts of preferential flow definitions

The boundary between preferential flow and Richards-type flow is a priori set at a volumetric soil water content θ at which soil water diffusivity D (θ) = η (= 10^− 6 m² s^− 1), where η is the kinematic viscosity. First we estimated with a hydrostatic approach from soil water retention curves the boundary, θ^K, between the structural pore domain, in which preferential flow occurs, and the matrix pore domain, in which Richards-type flow occurs. We then compared θ^K with θ that was derived from the respective soil hydrological property functions of same soil sample. Second, from in situ investigations we determined 96 values of θ^G as the terminal soil water contents that established themselves when the corresponding water-content waves of preferential flow have practically ceased. We compared the frequency distribution of θ^G with the one of θ that was calculated from the respective soil hydrological property functions of 32 soil samples that were determined with pressure plate apparatuses in the laboratory. There is support of the notion that θ^K ≈ θ^G ≈ θ, thus indicating the potential of θ to explain more generally what constitutes preferential flow. However, the support is assessed as working hypothesis on which to base further research rather than a procedure to a clear-cut identification of preferential flow and associated flow paths.     

Impact of initial and boundary conditions on preferential flow

Preferential flow in soil is approached by a water-content wave, WCW, that proceeds downward from the ground surface. WCWs were obtained from sprinkler experiments with infiltration rates varying from 5 to 40 mm h^− 1. TDR-probes and tensiometers measured volumetric water contents θ(z,t) at seven depths, and capillary heads, h(z,t) at six depths in a column of an undisturbed soil. The wave is characterized by the velocity of the wetting front, c_W, the amplitude, w_S, and the final water content, θ. We tested with uni-variate and bi-variate linear regressions the impacts of initial volumetric water contents, θ_ini, and input rates, q_S, on c_W, w_S and θ.The test showed that θ_ini influenced θ and w_S and q_S effected c_W. The expected proportionality of w_S ≈ qs^1/3 was weak and c_W ≈ qs^2/3 was strong.     

Incorporating layer- and local-scale heterogeneities in numerical simulation of unsaturated flow and tracer transport

This study characterizes layer- and local-scale heterogeneities in hydraulic parameters (i.e., matrix permeability and porosity) and investigates the relative effect of layer- and local-scale heterogeneities on the uncertainty assessment of unsaturated flow and tracer transport in the unsaturated zone of Yucca Mountain, USA. The layer-scale heterogeneity is specific to hydrogeologic layers with layerwise properties, while the local-scale heterogeneity refers to the spatial variation of hydraulic properties within a layer. A Monte Carlo method is used to estimate mean, variance, and 5th, and 95th percentiles for the quantities of interest (e.g., matrix saturation and normalized cumulative mass arrival). Model simulations of unsaturated flow are evaluated by comparing the simulated and observed matrix saturations. Local-scale heterogeneity is examined by comparing the results of this study with those of the previous study that only considers layer-scale heterogeneity. We find that local-scale heterogeneity significantly increases predictive uncertainty in the percolation fluxes and tracer plumes, whereas the mean predictions are only slightly affected by the local-scale heterogeneity. The mean travel time of the conservative and reactive tracers to the water table in the early stage increases significantly due to the local-scale heterogeneity, while the influence of local-scale heterogeneity on travel time gradually decreases over time. Layer-scale heterogeneity is more important than local-scale heterogeneity for simulating overall tracer travel time, suggesting that it would be more cost-effective to reduce the layer-scale parameter uncertainty in order to reduce predictive uncertainty in tracer transport.     

Influence of temporally variable groundwater flow conditions on point measurements and contaminant mass flux estimations

Monitoring of contaminant concentrations, e.g., for the estimation of mass discharge or contaminant degradation rates, often is based on point measurements at observation wells. In addition to the problem, that point measurements may not be spatially representative, a further complication may arise due to the temporal dynamics of groundwater flow, which may cause a concentration measurement to be not temporally representative. This paper presents results from a numerical modeling study focusing on temporal variations of the groundwater flow direction. “Measurements” are obtained from point information representing observation wells installed along control planes using different well frequencies and configurations. Results of the scenario simulations show that temporally variable flow conditions can lead to significant temporal fluctuations of the concentration and thus are a substantial source of uncertainty for point measurements. Temporal variation of point concentration measurements may be as high as the average concentration determined, especially near the plume fringe, even when assuming a homogeneous distribution of the hydraulic conductivity. If a heterogeneous hydraulic conductivity field is present, the concentration variability due to a fluctuating groundwater flow direction varies significantly within the control plane and between the different realizations. Determination of contaminant mass fluxes is also influenced by the temporal variability of the concentration measurement, especially for large spacings of the observation wells. Passive dosimeter sampling is found to be appropriate for evaluating the stationarity of contaminant plumes as well as for estimating average concentrations over time when the plume has fully developed. Representative sampling has to be performed over several periods of groundwater flow fluctuation. For the determination of mass fluxes at heterogeneous sites, however, local fluxes, which may vary considerably along a control plane, have to be accounted for. Here, dosimeter sampling in combination with time integrated local water flux measurements can improve mass flux estimates under dynamic flow conditions.     

基于超临界水氧化过程的能源环境系统设计

阐述了超临界水氧化过程的工艺路线,介绍了超临界水氧化过程的特点,提出了基于超临界水氧化过程的能源环境系统,设计了几种热量与能量回收系统程的耦合工艺,为提高超临界水氧化过程的经济性奠定了基础。     

Modeling potential herbicide loss to surface waters on the Swiss plateau

Lack of sufficiently detailed data often limits the applicability of complex transport-reaction models for estimating potential herbicide loss to surface waters. Therefore, there is also a need for simple models that are easy to apply but still capture the main features of the underlying processes.In this study, a simple regression model was developed to assess the vulnerability of catchments in the Swiss Plateau to diffuse herbicide loss to surface waters. The model is designed as a screening tool to rank the catchments in a relative sense and not to calculate Predicted Environmental Concentrations (PEC) of pesticides. The main goal is to capture two dominating factors controlling diffuse herbicide transport into streams and rivers. These factors are herbicide application and fast flow processes that are mainly responsible for herbicide transport. In a first step vulnerability of sites to herbicide loss is estimated based on site-specific conditions irrespective of actual herbicide application. In the second step, this vulnerability assessment is combined with actual herbicide application data to estimate the potential herbicide loss.The fast flow index (FFI), derived from discharge data using a base flow separation method, was applied as a proxy for the amount of fast flow occurring. The influence of catchment attributes (including topographic, climatic and soil data) on the FFI was analyzed using a multiple regression approach based on data from 57 catchments of the Swiss Plateau. By combining regression analysis with mechanistic knowledge, a two factor non-linear model based on river density and soil permeability as dominant input factors was selected as the best model for FFI prediction given the available data. Higher dimensional models had to be excluded because the strong correlation between the potential input factors led to unrealistic dependences while only minimally improving the quality of the fit.The spatial pattern of the predicted FFI as a measure for the vulnerability to diffuse herbicide losses shows a clearly increasing trend from the western to the eastern part of the Swiss Plateau and towards the pre-alpine/alpine regions in the south.In general the pattern of herbicide use corresponds to site conditions typical of a low FFI. However, the spatial analysis revealed exceptions, namely areas in which high actual herbicide use coincides with a high FFI.Despite the uncertainties in the model, this simple approach seems to be useful for supporting site-adapted agricultural practice whenever the higher accuracy of more detailed models is not required or too expensive to achieve. In addition, in combination with data on actual herbicide application, it can support the design of monitoring strategies by identifying critical areas of actual herbicide loss.