Multicomponent ionic diffusion in porewaters: Coulombic effects revisited

The diffusion of an ion in porewaters cannot occur independently of the other ions in solution as a result of Coulombic coupling, as well as from other effects not considered here. Unfortunately, a longstanding disagreement exists about the correct form and meaning of the equations that describe Coulombic coupling in porewaters, i.e., Ben-Yaakov [Am. J. Sci. 281 (1981) 974] vs. Lasaga [Am. J. Sci. 281 (1981) 981]. This paper re-examines this controversy by reformulating the problem starting from fundamental concepts of mass and charge conservation. We show that these antagonistic formulations are both valid and, in fact, equivalent, when the different interpretations of charge balance are resolved. Most of the disagreements between Ben-Yaakov and Lasaga are then shown to result from differing methods of solution, not fundamental disparities in their models. We note, however, that the explanation for the concept of “stationary” gradients of nonreacting ions as given Ben-Yaakov is inaccurate, and such gradients do lead to diffusive fluxes that are counterbalanced by electrochemical migrational fluxes to produce no net flux (excluding advective flux). We further find that the bicarbonate diffusive flux will not balance the diffusional charge flux of sulfate during its reduction if advection is present. This latter imbalance generates compensating fluxes in the other nonreacting ions. We have applied our theory to a simplified case of sulfate reduction in a marine sediment. The results show that nonreacting ions do diffuse and that with normally expected values of porewater advection, the ratio of the bicarbonate to the sulfate flux can be far different than the ideal value of −2.     

Accuracy of numerical simulations of contaminant transport in heterogeneous aquifers: A comparative study

This work deals with a comparison of different numerical schemes for the simulation of contaminant transport in heterogeneous porous media. The numerical methods under consideration are Galerkin finite element (GFE), finite volume (FV), and mixed hybrid finite element (MHFE). Concerning the GFE we use linear and quadratic finite elements with and without upwind stabilization. Besides the classical MHFE a new and an upwind scheme are tested. We consider higher order finite volume schemes as well as two time discretization methods: backward Euler (BE) and the second order backward differentiation formula BDF (2). It is well known that numerical (or artificial) diffusion may cause large errors. Moreover, when the Péclet number is large, a numerical code without some stabilising techniques produces oscillating solutions. Upwind schemes increase the stability but show more numerical diffusion. In this paper we quantify the numerical diffusion for the different discretization schemes and its dependency on the Péclet number. We consider an academic example and a realistic simulation of solute transport in heterogeneous aquifer. In the latter case, the stochastic estimates used as reference were obtained with global random walk (GRW) simulations, free of numerical diffusion. The results presented can be used by researchers to test their numerical schemes and stabilization techniques for simulation of contaminant transport in groundwater.     

The Influence of Element Diffusion at the Sediment—Water Interface on the Basic Chemical Composition of Overlying Water Mass in Lake Lugu,China

Little work has been done on the influence of seiments on the basic chemical composition of overlying water mass.This paper deals with the vertical profile of the basic constituents such as Ca^ ,K^ ,Na^ ,and HCO3^-,as well as of pH in the overlying water mass and sediment porewater of Lake Lugu-a semi-closed,deep lake in Yunnan Province.The reand sediment porewater of Lake Lugu- a semi-closed,deep lake in Yunnan Province.The results revealed that those basic constituents may diffuse and transport from bottom sediments to overlying water mass through porewater.In the paper are also quantitatively evaluated the diffusive fluxes and the extent of their influence on overlying water mass,indicating that the lake sediment-water interface diffusion plays an important role in controlling the basic chemical composition of water in the whole lake.     

基于扩散滤波的多尺度三维变分研究

将扩散方程引入三维变分分析,揭示了传统3D-VAR不能有效提取多尺度观测信息的根本原因,即观测导致的目标函数梯度在空间的不连续分布.将扩散滤波融入基于梯度的最优化算法,发展了基于扩散滤波的多尺度3D-VAR.海表面温度数据同化试验结果表明,新方法可从长波至短波有效地提取多尺度的观测信息.     

Local And Global Diffusion Along Resonant Lines in Discrete Quasi-integrable Dynamical Systems

We detect and measure diffusion along resonances in a quasi-integrable symplectic map for different values of the perturbation parameter. As in a previously studied Hamiltonian case (Lega et al., 2003) results agree with the prediction of the Nekhoroshev theorem. Moreover, for values of the perturbation parameter slightly below the critical value of the transition between Nekhoroshev and Chirikov regime we have also found a diffusion of some orbits along macroscopic portions of the phase space. Such a diffusion follows in a spectacular way the peculiar structure of resonant lines.     

Comparison of parallel group velocity of ECH waves with electron resonant velocity: Implication for electron diffusion by ECH waves

We have investigated the role of group velocity in the calculation of pitch-angle diffusion coefficients by electron cyclotron harmonic (ECH) waves in planetary magnetospheres. The assumption which is generally made that the parallel group velocity can be neglected in comparison with particle parallel velocity is examined in detail. It is found that for lowest harmonic band this assumption is quite good. It is found that in general it is not possible to ignore the parallel group velocity. However, for lowest harmonic band this assumption is quite good at low electron temperatures.     

A discontinuous finite element suspended sediment transport model for water quality assessments in river networks

Suspended sediment plays an important role in the distribution and transport of many pollutants (such as radionuclides) in rivers. Pollutants may adsorb on fine suspended particles (e.g. clay) and spread according to the suspended sediment movement. Hence, the simulation of the suspended sediment mechanism is indispensable for realistic transport modelling. This paper presents and tests a simple mathematical model for predicting the suspended sediment transport in river networks. The model is based on the van Rijn suspended load formula and the advection–diffusion equation with a source or sink term that represents the erosion or deposition fluxes. The transport equation is solved numerically with the discontinuous finite element method. The model evaluation was performed in two steps, first by comparing model simulations with the measured suspended sediment concentrations in the Grote Nete–Molse Nete River in Belgium, and second by a model intercomparison with the sediment transport model NST MIKE 11. The simulations reflect the measurements with a Nash‐Sutcliffe model efficiency of 0.6, while the efficiency between the proposed model and the NST MIKE 11 simulations is 0.96. Both evaluations indicate that the proposed sediment transport model, that is sufficiently simple to be practical, is providing realistic results.     

井水温度微动态及其形成机制

介绍了我国地震地下水温度观测概况及已获得的微动态现象, 即水温固体潮、 水温同震效应、 水温前兆异常等, 并对其形成机理做了理论分析。 提出水温微动态形成机理研究应从二个地温梯度、 二个对流、 一个传导、 一个扩散等热系统与热过程去考虑。     

基于镭同位素分布的黄海和东海垂直混合速率计算

在黄海和东海采样测定了水体中的镭同位素分布,用平流扩散模型描述镭同位素分布,最小二乘方法计算了垂直涡动扩散系数和上升流或下降流流速.结果给出北黄海中部、南黄海中部、浙江沿岸和台湾北部海域存在上升流,流速分别为0.46×10-3cm·s-1、0.17×10-3~1.39×10-3cm·s-1、2.02×10-3~3.04×10-3 cm·s-1和1.06×10-3~2.51×10-3 cm·s-1.北黄海中部和东海东北部存在下降流.流速分别为-2.30×10-3 cm·s-1和-0.61×10-3~-2.10×10-3 cm·s-1.计算同时给出的垂直涡动扩散系数为5.84~48.2 cm2·s-1,平均值为22.3 cm2·s-1.北黄海和浙江沿岸上升流流速与文献的结果一致;北黄海中部存在下降流与文献的结论一致.本研究结果与文献结果一致是对所建立的方法的肯定,也是对文献研究结果的支持.     

各向异性扩散在电磁测深曲线平滑中的应用

电磁类地球物理方法由于极易受到各类噪声的干扰,使得估计的视电阻率曲线或相位曲线发生畸变,严重地影响了反演解译精度,如何对这类观测曲线进行合理有效地平滑,是目前数据处理中的重点,也是难点.本文把各向异性扩散（anisotropic diffusion AD）引入到曲线平滑中,提出了各向异性扩散的曲线平滑方法,以平滑点梯度值的降函数作为扩散速度,在梯度较大的位置予以较小的扩散速度以保护特征,在梯度较小的位置予以较大的扩散速度平滑噪声引起的扰动.同时为了减小“飞点”(outliers)对平滑的干扰,采用局部单调扩散进行预处理.实验结果显示,本文提出的平滑方法有效地平滑了噪声与“飞点”的干扰,恢复了曲线的基本形态,有效地保持了曲线基本特征.