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.     

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.     

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.     

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.     

一种新的扇区拥挤告警指标及应用

针对我国管制扇区日益繁忙拥挤、空管保障系统的相对落后,扇区流量告警机制尚不健全的现状,在分析和研究美国、欧洲、日本扇区流量告警机制基础上,提出当量流量告警指标的概念;采用MMBB工作负荷评估方法确定不同特征航班对应的当量系数,计算出扇区的当量流量,建立一种新型的扇区流量告警方法;应用算例的结果验证了当量流量告警指标在扇区流量告警工作中的有效性和实用性。     

基于仿真的高速公路紧急事件交通流特征研究

研究高速公路紧急事件下的交通流特征,对于减少高速公路拥挤、降低二次事故、提高通行能力有重要的作用。由于高速公路紧急事件具有偶然性、不可预测性的特点,在实际调研中往往难以及时捕捉到交通事件的相应信息,导致源数据不足进而影响研究的进行。本文通过微观仿真手段,结合实地调研获得紧急状态下车辆运行参数,在AIMSUN仿真环境下加载紧急事件,建立了高速公路微观路段的仿真模型。并将路段平均车速、密度、最大排队长度等作为指标,定昔计算紧急事件下各种交通参数的变化规律。最后与调研结果进行了对比,表明仿真结果与实际情况较为一致。