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.     

Evaluation of a Watershed Model for Estimating Daily Flow Using Limited Flow Measurements1

Abstract: The Soil and Water Assessment Tool (SWAT) model was evaluated for estimation of continuous daily flow based on limited flow measurements in the Upper Oyster Creek (UOC) watershed. SWAT was calibrated against limited measured flow data and then validated. The Nash‐Sutcliffe model Efficiency (NSE) and mean relative error values of daily flow estimations were 0.66 and 15% for calibration, and 0.56 and 4% for validation, respectively. Also, further evaluation of the model’s estimation of flow at multiple locations was conducted with parametric paired t‐test and nonparametric sign test at a 95% confidence level. Among the five main stem stations, four stations were statistically shown to have good agreement between predicted and measured flows. SWAT underestimated the flow of the fifth main stem station possibly because of the existence of complex flood control measures near to the station. SWAT estimated the daily flow at one tributary station well, but with relatively large errors for the other two tributaries. The spatial pattern of predicted flows matched the measured ones well. Overall, it was concluded from the graphical comparisons and statistical analyses of the model results that SWAT was capable of reproducing continuous daily flows based on limited flow data as is the case in the UOC watershed.     

Model prediction on the rise of pCO2 in uniform flows by leakage of CO2 purposefully stored under the seabed

A numerical study was conducted to predict pCO₂ change in the ocean on a continental shelf by the leakage of CO₂, which is originally stored in the aquifer under the seabed, in the case that a large fault connects the CO₂ reservoir and the seabed by an earthquake or other diastrophism. The leakage rate was set to be 6.025 × 10⁻⁴ kg/m²/sec from 2 m × 100 m fault band, which corresponds to 3800 t-CO₂/year, referring to the monitored seepage rate from an existing EOR field. The target space in this study was limited to the ocean above the seabed, the depth of which was 200 or 500 m. The computational domain was idealistically rectangular with the seabed fault-band perpendicular to the uniform flow. The CO₂ takes a form of bubbles or droplets, depending on the depth of water, and their behaviour and dissolution were numerically simulated during their rise in seawater flow. The advection–diffusion of dissolved CO₂ was also simulated. As a result, it was suggested that the leaked CO₂ droplets/bubbles all dissolve in the seawater before spouting up to the atmosphere, and that the increase in pCO₂ in the seawater was smaller than 500 μ atm.     

The contemporary cement cycle of the United States

A country-level stock and flow model for cement, an important construction material, was developed based on a material flow analysis framework. Using this model, the contemporary cement cycle of the United States was constructed by analyzing production, import, and export data for different stages of the cement cycle. The United States currently supplies approximately 80% of its cement consumption through domestic production and the rest is imported. The average annual net addition of in-use new cement stock over the period 2000–2004 was approximately 83 million metric tons and amounts to 2.3 tons per capita of concrete. Nonfuel carbon dioxide emissions (42 million metric tons per year) from the calcination phase of cement manufacture account for 62% of the total 68 million tons per year of cement production residues. The end-of-life cement discards are estimated to be 33 million metric tons per year, of which between 30% and 80% is recycled. A significant portion of the infrastructure in the United States is reaching the end of its useful life and will need to be replaced or rehabilitated; this could require far more cement than might be expected from economic forecasts of demand for cement.