Assessment of alternative adsorption models and global sensitivity analysis to characterize hexavalent chromium loss from soil to surface runoff

We investigate our ability to assess transfer of hexavalent chromium, Cr(VI), from the soil to surface runoff by considering the effect of coupling diverse adsorption models with a two‐layer solute transfer model. Our analyses are grounded on a set of two experiments associated with soils characterized by diverse particle size distributions. Our study is motivated by the observation that Cr(VI) is receiving much attention for the assessment of environmental risks due to its high solubility, mobility, and toxicological significance. Adsorption of Cr(VI) is considered to be at equilibrium in the mixing layer under our experimental conditions. Four adsorption models, that is, the Langmuir, Freundlich, Temkin, and linear models, constitute our set of alternative (competing) mathematical formulations. Experimental results reveal that the soil samples characterized by the finest grain sizes are associated with the highest release of Cr(VI) to runoff. We compare the relative abilities of the four models to interpret experimental results through maximum likelihood model calibration and four model identification criteria (i.e., the Akaike information criteria [AIC and AIC_C] and the Bayesian and Kashyap information criteria). Our study results enable us to rank the tested models on the basis of a set of posterior weights assigned to each of them. A classical variance‐based global sensitivity analysis is then performed to assess the relative importance of the uncertain parameters associated with each of the models considered, within subregions of the parameter space. In this context, the modelling strategy resulting from coupling the Langmuir isotherm with a two‐layer solute transfer model is then evaluated as the most skilful for the overall interpretation of both sets of experiments. Our results document that (a) the depth of the mixing layer is the most influential factor for all models tested, with the exception of the Freundlich isotherm, and (b) the total sensitivity of the adsorption parameters varies in time, with a trend to increase as time progresses for all of the models. These results suggest that adsorption has a significant effect on the uncertainty associated with the release of Cr(VI) from the soil to the surface runoff component.     

Variable stars in the field of open cluster NGC 6819

We report on the discovery of 25 variable stars plus 13 suspected variables found in the field of the open cluster NGC 6819. The stars were identified from time-series photometric data obtained on the Isaac Newton Telescope, La Palma, during two observing runs covering the 19 nights between 1999 June  22–30  and 1999 July  22–31  . The variables found include 12 eclipsing binaries with an additional three suspected, nine BY Draconis systems, plus four variables of other types, including one star believed to be a Cepheid. Three of the 15 eclipsing binaries are believed to be cluster members. Details of a further 10 suspected variable stars are also included.     

砂泥岩裂隙岩体埋深和岩性对渗透性影响分析

采用表征岩体渗透性的单位吸水量ω为参数,有效地分析了砂泥岩裂隙岩体中相对埋深和岩性对渗透性空间分布规律的影响.研究区内,在不同深度单位吸水量的最大值和最小值相差3-4个数量级,渗透性表现出显著的随机特征.通过将压水试验数据离散化并在不同深度上取logω的平均值,凸现了渗透性的结构特征.ω均值在一定范围内随相对深度表现出负指数分布规律;剔除logω的线性趋势项后,logω残差表现出与地层内部岩性相对应的规律.此外,研究还发现,渗透性参数随相对深度呈负指数分布规律是由含裂隙的砂岩介质引起的.这些规律为裂隙岩体的渗流计算和模拟奠定了基础.文章所提出的对压水试验数据离散化并在不同深度取均值以找出渗透性和深度关系,以及去除线性趋势项研究渗透性和岩性之间关系等方法具有广泛的应用价值.     

Using Sequential Self-Calibration Method to Identify Conductivity Distribution: Conditioning on Tracer Test data

An iterative inverse method, the sequential self-calibration method, is developed for mapping spatial distribution of a hydraulic conductivity field by conditioning on nonreactive tracer breakthrough curves. A streamline-based, semi-analytical simulator is adopted to simulate solute transport in a heterogeneous aquifer. The simulation is used as the forward modeling step. In this study, the hydraulic conductivity is assumed to be a deterministic or random variable. Within the framework of the streamline-based simulator, the efficient semi-analytical method is used to calculate sensitivity coefficients of the solute concentration with respect to the hydraulic conductivity variation. The calculated sensitivities account for spatial correlations between the solute concentration and parameters. The performance of the inverse method is assessed by two synthetic tracer tests conducted in an aquifer with a distinct spatial pattern of heterogeneity. The study results indicate that the developed iterative inverse method is able to identify and reproduce the large-scale heterogeneity pattern of the aquifer given appropriate observation wells in these synthetic cases.     

Polyoxymethylene co-polymers on grains

Conditions prevalent in dense molecular clouds are shown to favour the polymerization of H₂CO molecules and the deposition of formaldehyde co-polymer mantles, with typical radii 10^–5 cm, on smaller refractory grains. If a significant fraction of such co-polymer coated grains are expelled with systematic gas flows into the general interstellar medium, these moderately refractory grains may be responsible for the bulk of interstellar extinction and polarization at optical wavelengths. Mie calculations for a mixture consisting of iron, graphite and POM particles are presented as an example where POM grains of radii 0.15 dominate the extinction at optical wavelengths, providing a satisfactory overall fit to a range of extinction data. A size distribution of POM needles with a mean radius 0.15 also provides good agreement with data on interstellar linear as well as circular polarization. Suitably end-capped and stabilized co-polymer-coated grains, with either silicate or graphite cores, may survive at temperatures 450 K under interstellar ambient conditions and be responsible for the 10 emission feature in many sources. Theoretically computed band profiles of the 10 -feature in POM coated grains, in general, provide better agreement with observations than most types of silicate grains considered so far. We also note that an unexplained dip at 10 in the 8–12 feature of the infrared source OH 231.8+4.2 may be a signature of POM grains; likewise, a persistent 3.3 emission feature in many different types of infrared source could be attributed to the CH stretching mode in formaldehyde co-polymer grains.     

African monsoon variability during the previous interglacial maximum

Little is known about centennial- to millennial-scale climate variability during interglacial times, other than the Holocene. We here present high-resolution evidence from anoxic (unbioturbated) sediments in the eastern Mediterranean Sea that demonstrates a sustained ∼800-yr climate disturbance in the monsoonal latitudes during the Eemian interglacial maximum (∼125 ka BP). Results imply that before and after this event, the Intertropical Convergence Zone (ITCZ) penetrated sufficiently beyond the central Saharan watershed (∼21°N) during the summer monsoon to fuel flooding into the Mediterranean along the wider North African margin, through fossil river/wadi systems that to date have been considered only within a Holocene context. Relaxation in the ITCZ penetration during the intra-Eemian event curtailed this flux, but flow from the Nile - with its vast catchment area - was not affected. Previous work suggests a concomitant Eurasian cooling event, with intensified impact of the higher-latitude climate on the Mediterranean basin. The combined signals are very similar to those described for the Holocene cooling event around 8 ka BP. The apparent type of concurrent changes in the monsoon and higher-latitude climate may reflect a fundamental mechanism for variability in the transfer of energy (latent heat) between the tropics and higher latitudes.     

Can flat-ramp-flat fault geometry be inferred from fold shape?: A comparison of kinematic and mechanical folds

The inference of fault geometry from suprajacent fold shape relies on consistent and verified forward models of fault-cored folds, e.g. suites of models with differing fault boundary conditions demonstrate the range of possible folding. Results of kinematic (fault-parallel flow) and mechanical (boundary element method) models are compared to ascertain differences in the way the two methods simulate flexure associated with slip along flat-ramp-flat geometry. These differences are assessed by systematically altering fault parameters in each model and observing subsequent changes in the suprajacent fold shapes. Differences between the kinematic and mechanical fault-fold relationships highlight the differences between the methods. Additionally, a laboratory fold is simulated to determine which method might best predict fault parameters from fold shape. Although kinematic folds do not fully capture the three-dimensional nature of geologic folds, mechanical models have non-unique fold-fault relationships. Predicting fault geometry from fold shape is best accomplished by a combination of the two methods.