Predicting DNAPL entrapment and recovery: the influence of hydraulic property correlation

The influence of aquifer property correlation on multiphase fluid migration, entrapment and recovery was explored by incorporating correlated and uncorrelated porosity, permeability, and capillary pressure-saturation (P_c-Sat) parameter fields in a cross-sectional numerical multiphase flow model. Comparison of two-dimensional entrapped organic saturation distributions for a simulated tetrachloroethylene (PCE) spill in ensembles of aquifer realizations suggests that the degree of spatial correlation in P_c-Sat parameters exerts a controlling influence on dense nonaqueous phase liquid (DNAPL) spreading and redistribution in saturated aquifers. The predicted evolution of DNAPL source zones and resultant remediation efficiency under surfactant enhanced aquifer remediation (SEAR) also appear to be strongly influenced by the spatial correlation of aquifer parameters and multiphase flow constitutive relationships. Results for a limited number of realizations selected from each ensemble showed that removal of 60% to 99% of entrapped PCE could reduce dissolved contaminant concentration and mass flux by approximately two orders of magnitude under natural gradient conditions. Aqueous phase contaminant mass flux did not vary uniformly as a function of % DNAPL removed, however, and notable differences in behavior were observed for models incorporating correlated versus uncorrelated P_c-Sat and permeability fields. Although these results must be confirmed through analysis of additional realizations, it is likely that similar or larger differences between correlated and uncorrelated system behavior will be observed in aquifers with greater spatially variability than that of the nonuniform, homogeneous sand aquifer studied here. Funding for this research was provided by the United States Environmental Protection Agency, Great Lakes and Mid-Atlantic Center for Hazardous Substance Research under Grant No. R-825540, the Michigan Department of Environmental Quality under Contract No. Y80011, and the Strategic Environmental Research and Development Program under Project No. CU-1293. The content of this publication does not necessarily represent the views of these agencies and has not been subject to agency review.     

Bald kelp: natural and harvesting-induced frond loss in the South African kelp Ecklonia maxima

Ecklonia maxima dominates kelp beds on the west coast of southern Africa, and is commercially and ecologically valuable. Some plants lose all their secondary blades (fronds), leaving only the gas-filled bulb and short primary blade at the top of the stipe. These ‘bald’ kelps may persist for months or longer, occupying substratum space but contributing little to productivity. We investigated natural causes of balding in five kelp beds in False Bay, South Africa, and the effect of simulated commercial frond-harvesting on loss of fronds. Densities of sporophytes with and without fronds were measured at the water surface during low tides. Generalised linear modelling analysis showed a significant relationship between balding and site (whether sheltered from or exposed to swell/wind), position in the kelp bed (shoreward or seaward) and total E. maxima density. We conclude that balding is caused mainly by drying of meristematic basal regions of fronds during emersion of kelp heads at low tide. After partial harvesting of fronds, kelp heads floated higher out of the water, but after 95 days frond loss was significantly higher only when fronds were cut too close to the primary blade, possibly because of damage to meristematic tissue. Nevertheless, increased emersion that results from frond-harvesting may increase desiccation and blade loss and we suggest caution in setting harvest limits for this kelp.     

桐柏碰撞造山带及其邻区变形特征与构造格局

桐柏碰撞造山带及其邻区可以划分为九个大地构造单元,自北向南分别是:华北克拉通南缘岩石构造单元——宽坪岩群、具弧后盆地性质的二郎坪岩石构造单元、具岛弧性质的秦岭杂岩单元、龟山岩组和南湾岩组构成的俯冲前缘楔构造带、构造混杂岩带、桐柏北部高压岩片单元、桐柏核部杂岩单元、桐柏南部高压岩片单元以及随州构造变形带。根据详细的构造解析以及新的地质年代学资料,本文将中生代以来的构造变形划分为五幕,前两幕变形主要发育在构造混杂岩带以南的各个岩石构造单元中,之后的三幕变形则波及整个研究区。第一幕变形的时间约为255～238Ma,以发育区域上透入性的片理及北西西向的拉伸线理为主,并导致了高压岩片早期自西向东的挤出。第二幕变形的时间约为230～215Ma,以自北向南的逆冲推覆构造为主,使得高压岩片进一步垂向抬升。第三幕变形应早于下侏罗统,以近北西西向的宽缓褶皱为主要特征,该幕变形期间桐柏核部杂岩及其两侧高压岩片单元发生同步的抬升。第四幕变形大致发生在140～130Ma之间,主要表现为桐柏核部杂岩两侧走滑型韧性剪切带的活动,桐柏核部杂岩表现出向东的挤出。第五幕变形发生在120～80Ma,表现为北西向及北东向的脆性断裂活动,并切割以上所有构造形迹。桐柏高压岩片的抬升剥露受多幕变形控制,呈阶段性的抬升。     

Matching solute breakthrough with deterministic and stochastic aquifer models

Two different deterministic and two alternative stochastic (i.e., geostatistical) approaches to modeling the distribution of hydraulic conductivity (K) in a nonuniform (sigma2ln(K)) = 0.29) glacial sand aquifer were used to explore the influence of conceptual model selection on simulations of three-dimensional tracer movement. The deterministic K models employed included a homogeneous effective K and a perfectly stratified 14 layer model. Stochastic K models were constructed using sequential Gaussian simulation and sequential i ndicator simulation conditioned to available K values estimated from measured grain size distributions. Standard simulation software packages MODFLOW, MT3DMS, and MODPATH were used to model three-dimensional ground water flow and transport in a field tracer test, where a pulse of bromide was injected through an array of three fully screened wells and extracted through a single fully screened well approximately 8 m away. Agreement between observed and simulated transport behavior was assessed through direct comparison of breakthrough curves (BTCs) and selected breakthrough metrics at the extraction well and at 26 individual multilevel sample ports distributed irregularly between the injection and extraction wells. Results indicate that conceptual models incorporating formation variability are better able to capture observed breakthrough behavior. Root mean square (RMS) error of the deterministic models bracketed the ensemble mean RMS error of stochastic models for simulated concentration vs. time series, but not for individual BTC characteristic metrics. The spatial variability models evaluated here may be better suited to simulating breakthrough behavior measured in wells screened over large intervals than at arbitrarily distributed observation points within a nonuniform aquifer domain.     

The influence of dimensionality on simulations of mass recovery from nonuniform dense non-aqueous phase liquid (DNAPL) source zones

The influence of model dimensionality on predictions of mass recovery from dense non-aqueous phase liquid (DNAPL) source zones in nonuniform permeability fields was investigated using a modified version of the modular three-dimensional transport simulator (MT3DMS). Thirty-two initial two- (2D) and three-dimensional (3D) tetrachloroethene–DNAPL source zone architectures, taken from a recent modeling study, were used as initial conditions for this analysis. Commonly employed source zone metrics were analyzed to determine differences between 2D and 3D predictions: (i) down-gradient flux-averaged contaminant concentration, (ii) reductions in contaminant mass flux through a down-gradient boundary, (iii) source zone ganglia-to-pool (GTP) ratio, and (iv) time required to achieve a remediation objective. 3D flux-averaged contaminant concentrations were approximately 3.5 times lower than concentrations simulated in 2D. This difference was attributed to dilution of the contaminant concentrations down gradient of the source zone. Contaminant flux reduction predictions for a given mass recovery were generally 5% higher in 3D simulations than in 2D simulations. The GTP ratio declined over time as mass was recovered in both 2D and 3D simulations. Although the source longevity (i.e., time required to achieve 99.99% mass recovery) differed between individual 2D and 3D realizations, the mean source longevity for the 2D and 3D simulation ensembles was within 2%. 2D simulations tended to over-predict the time required to achieve lower mass recovery levels (e.g. 50% mass recovery) due to a smaller contaminated area exposed to uncontaminated water. These findings suggest that ensemble averages of 2D numerical simulations of DNAPL migration, entrapment, dissolution, and mass recovery in statistically homogenous, nonuniform media may provide reasonable approximations to average behavior obtained using simulations conducted in fully three-dimensional domains.     

陕西省略阳县杨家坝多金属矿区成矿作用地球化学示踪

本文以陕西省略阳县杨家坝多金属矿区为例,通过对该矿区岩(矿)石和各类脉体的岩相学、稀土和微量元素及流体包裹体地球化学示踪发现,多金属中重晶石的形成及硅化、碳酸盐化、滑石化、纤闪石化、绢云母化和绿泥石化和透闪石岩的含氟浅闪石及隐晶硅质玻璃和微晶硅质,这表明与成矿蚀变有关的热液流体是一种不同于一般地壳流体的富硅、钛、铁、碱质和挥发份,并具备熔体性质的成矿流体。各类样品的稀土元素配分模式明显富集LREE,在Ce、Eu、Yb异常的组合上也各具特征。LREE富集是地幔流体作用的显著特征之一;矿石和脉体的负Ce、正Eu和负Yb异常以及流体包裹体中H2O-C6H6成分,是高温还原地幔流体的重要标志;正Ce和负Eu异常的出现,以及Ce、Eu、Yb异常的减弱和消失则是壳幔混染叠加改造的显示。综合研究表明:矿区成矿过程可能统一受制于秦岭地区碰撞造山背景下具高温还原性质的地幔流体作用,且由此引发壳幔强烈混染的叠加改造作用在成矿过程中发挥了重要作用。     

Bio-treatment and water reuse as feasible treatment approaches for improving wastewater management during flotation of copper ores

Established and emerging technologies for treatment of flotation wastewaters are compared and discussed in the context of their applicability during the management of wastewaters from Cu–Co ores processing through flotation in the Katanga province of the Democratic Republic of Congo. The effects from water recycling on water quality and flotation performance are briefly presented in parallel. The ore processing schemes and the wastewater management practices at two operational concentrators are illustrated as study cases and their peculiarities outlined. A reference to a copper concentrator in a nearby Zambia is addressed for comparative purpose. Based on analysis of the findings, the clarification, bio-treatment inside the tailing ponds or the use of constructed wetlands as polishing stage prior to water reuse are suggested as feasible treatment approaches in view improving the management practice of flotation wastewaters during the dressing of copper ores in the Katanga province.