Modeling of ore-forming and geoenvironmental systems: Roles of fluid flow and chemical reaction processes

Ore-forming and geoenviromental systems commonly involve coupled fluid flowand chemical reaction processes. The advanced numerical methods and computational modeling have become indispensable tools for simulating such processes in recent years. This enables many hitherto unsolvable geoscience problems to be addressed using numerical methods and computational modeling approaches. For example, computational modeling has been successfully used to solve ore-forming and mine site contamination/remediation problems, in which fluid flow and geochemical processes play important roles in the controlling dynamic mechanisms. The main purpose of this paper is to present a generalized overview of: (1) the various classes and models associated with fluid flow/chemically reacting systems in order to highlight possible opportunities and developments for the future; (2) some more general issues that need attention in the development of computational models and codes for simulating ore-forming and geoenviromental systems; (3) the related progresses achieved on the geochemical modeling over the past 50 years or so; (4) the general methodology for modeling of oreforming and geoenvironmental systems; and (5) the future development directions associated with modeling of ore-forming and geoenviromental systems.

Reactive mass transport modelling of a three-dimensional vertical fault zone with a finger-like convective flow regime

For a three-dimensional vertically-oriented fault zone, we consider the coupled effects of fluid flow, heat transfer and reactive mass transport, to investigate the patterns of fluid flow, temperature distribution, mineral alteration and chemically induced porosity changes. We show, analytically and numerically, that finger-like convection patterns can arise in a vertically-oriented fault zone. The onset and patterns of convective fluid flow are controlled by the Rayleigh number which is a function of the thermal properties of the fluid and the rock, the vertical temperature gradient, and the height and the permeability of the fault zone. Vigorous fluid flow causes low temperature gradients over a large region of the fault zone. In such a case, flow across lithological interfaces becomes the most important mechanism for the formation of sharp chemical reaction fronts. The degree of rock buffering, the extent and intensity of alteration, the alteration mineralogy and in some cases the formation of ore deposits are controlled by the magnitude of the flow velocity across these compositional interfaces in the rock. This indicates that alteration patterns along compositional boundaries in the rock may provide some insights into the convection pattern. The advective mass and heat exchanges between the fault zone and the wallrock depend on the permeability contrast between the fault zone and the wallrock. A high permeability contrast promotes focussed convective flow within the fault zone and diffusive exchange of heat and chemical reactants between the fault zone and the wallrock. However, a more gradual permeability change may lead to a regional-scale convective flow system where the flow pattern in the fault affects large-scale fluid flow, mass transport and chemical alteration in the wallrocks

ROSINA/DFMS and IES observations of 67P: Ion-neutral chemistry in the coma of a weakly outgassing comet

Context. The Rosetta encounter with comet 67P/Churyumov-Gerasimenko provides a unique opportunity for an in situ, up-close investigation of ion-neutral chemistry in the coma of a weakly outgassing comet far from the Sun. Aims. Observations of primary and secondary ions and modeling are used to investigate the role of ion-neutral chemistry within the thin coma. Methods. Observations from late October through mid-December 2014 show the continuous presence of the solar wind 30 km from the comet nucleus. These and other observations indicate that there is no contact surface and the solar wind has direct access to the nucleus. On several occasions during this time period, the Rosetta/ROSINA/Double Focusing Mass Spectrometer measured the low-energy ion composition in the coma. Organic volatiles and water group ions and their breakup products (masses 14 through 19), CO2+ (masses 28 and 44) and other mass peaks (at masses 26, 27, and possibly 30) were observed. Secondary ions include H3O+ and HCO+ (masses 19 and 29). These secondary ions indicate ion-neutral chemistry in the thin coma of the comet. A relatively simple model is constructed to account for the low H3O+/H2O+ and HCO+/CO+ ratios observed in a water dominated coma. Results from this simple model are compared with results from models that include a more detailed chemical reaction network. Results. At low outgassing rates, predictions from the simple model agree with observations and with results from more complex models that include much more chemistry. At higher outgassing rates, the ion-neutral chemistry is still limited and high HCO+/CO+ ratios are predicted and observed. However, at higher outgassing rates, the model predicts high H3O+/H2O+ ratios and the observed ratios are often low. These low ratios may be the result of the highly heterogeneous nature of the coma, where CO and CO2 number densities can exceed that of water.

Proximity-accelerated chemical coupling reaction in the benzodiazepine-binding site of gamma-aminobutyric acid type A receptors: superposition of different allosteric modulators

Benzodiazepines are widely used drugs. They exert sedative/hypnotic, anxiolytic, muscle relaxant, and anticonvulsant effects and act through a specific high affinity binding site on the major inhibitory neurotransmitter receptor, the gamma-aminobutyric acid type A (GABA(A)) receptor. Ligands of the benzodiazepine-binding site are classified into three groups depending on their mode of action: positive and negative allosteric modulators and antagonists. To rationally design ligands of the benzodiazepine site in different isoforms of the GABA(A) receptor, we need to understand the relative positioning and overlap of modulators of different allosteric properties. To solve these questions, we used a proximity-accelerated irreversible chemical coupling reaction. GABA(A) receptor residues thought to reside in the benzodiazepine-binding site were individually mutated to cysteine and combined with a cysteine-reactive benzodiazepine site ligand. Direct apposition of reaction partners is expected to lead to a covalent reaction. We describe here such a reaction of predominantly alpha(1)H101C and also three other mutants (alpha(1)G157C, alpha(1)V202C, and alpha(1)V211C) with an Imid-NCS derivative in which a reactive isothiocyanate group (-NCS) replaces the azide group (-N(3)) in the partial negative allosteric modulator Ro15-4513. Our results show four contact points of imidazobenzodiazepines with the receptor, alpha(1)H101C being shared by classical benzodiazepines. Taken together with previous data, a similar orientation of these ligands within the benzodiazepine-binding pocket may be proposed.

Molecular crosstalk between a chemical and a biological stressor and consequences on disease manifestation in rainbow trout.

The aim of the present study was to examine the molecular and organism reaction of rainbow trout, Oncorhynchus mykiss, to the combined impact of two environmental stressors. The two stressors were the myxozoan parasite, Tetracapsuloides bryosalmonae, which is the etiological agent of proliferative kidney disease (PKD) and a natural stressor to salmonid populations, and 17β-estradiol (E2) as prototype of estrogen-active chemical stressors in the aquatic environment. Both stressors, the parasite and estrogenic contaminants, co-exist in Swiss rivers and are discussed as factors contributing to the decline of Swiss brown trout populations over the last decades. Using a microarray approach contrasting parasite-infected and non-infected rainbow trout at low or high estrogen levels, it was observed that molecular response patterns under joint exposure differed from those to the single stressors. More specifically, three major response patterns were present: (i) expression responses of gene transcripts to one stressor are weakened by the presence of the second stressor; (ii) expression responses of gene transcripts to one stressor are enhanced by the presence of the second stressor; (iii) expression responses of gene transcripts at joint treatment are dominated by one of the two stressors. Organism-level responses to concurrent E2 and parasite treatment - assessed through measuring parasite loads in the fish host and cumulative mortalities of trout - were dominated by the pathogen, with no modulating influence of E2. The findings reveal function- and level-specific responses of rainbow trout to stressor combinations, which are only partly predictable from the response to the single stressors.

A review of air–ice chemical and physical interactions (AICI): liquids, quasi-liquids, and solids in snow

Snow in the environment acts as a host to rich chemistry and provides a matrix for physical exchange of contaminants within the ecosystem. The goal of this review is to summarise the current state of knowledge of physical processes and chemical reactivity in surface snow with relevance to polar regions. It focuses on a description of impurities in distinct compartments present in surface snow, such as snow crystals, grain boundaries, crystal surfaces, and liquid parts. It emphasises the microscopic description of the ice surface and its link with the environment. Distinct differences between the disordered air–ice interface, often termed quasi-liquid layer, and a liquid phase are highlighted. The reactivity in these different compartments of surface snow is discussed using many experimental studies, simulations, and selected snow models from the molecular to the macro-scale. Although new experimental techniques have extended our knowledge of the surface properties of ice and their impact on some single reactions and processes, others occurring on, at or within snow grains remain unquantified. The presence of liquid or liquid-like compartments either due to the formation of brine or disorder at surfaces of snow crystals below the freezing point may strongly modify reaction rates. Therefore, future experiments should include a detailed characterisation of the surface properties of the ice matrices. A further point that remains largely unresolved is the distribution of impurities between the different domains of the condensed phase inside the snowpack, i.e. in the bulk solid, in liquid at the surface or trapped in confined pockets within or between grains, or at the surface. While surface-sensitive laboratory techniques may in the future help to resolve this point for equilibrium conditions, additional uncertainty for the environmental snowpack may be caused by the highly dynamic nature of the snowpack due to the fast metamorphism occurring under certain environmental conditions. Due to these gaps in knowledge the first snow chemistry models have attempted to reproduce certain processes like the long-term incorporation of volatile compounds in snow and firn or the release of reactive species from the snowpack. Although so far none of the models offers a coupled approach of physical and chemical processes or a detailed representation of the different compartments, they have successfully been used to reproduce some field experiments. A fully coupled snow chemistry and physics model remains to be developed.

Probing the Electrocatalytic Oxygen Reduction Reaction Reactivity of Immobilized Multicopper Oxidase CueO

The bioelectrocatalytic (oxygen reduction reaction, ORR) properties of the multicopper oxidase CueO immobilized on gold electrodes were investigated. Macroscopic electrochemical techniques were combined with in situ scanning tunneling microscopy (STM) and surface-enhanced Raman spectroscopy at the ensemble and at the single-molecule level. Self-assembled monolayer of mercaptopropionic acid, cysteamine, and p-aminothiophenol were chosen as redox mediators. The highest ORR activity was observed for the protein attached to amino-terminated adlayers. In situ STM experiments revealed that the presence of oxygen causes distinct structure and electronic changes in the metallic centers of the enzyme, which determine the rate of intramolecular electron transfer and, consequently, affect the rate of electron tunneling through the protein. Complementary Raman spectroscopy experiments provided access for monitoring structural changes in the redox state of the type 1 copper center of the immobilized enzyme during the CueO-catalyzed oxygen reduction cycle. These results unequivocally demonstrate the existence of a direct electronic communication between the electrode substrate and the type 1 copper center.

Synthetic studies towards a novel, chemical stable, abasic site analogue of DNA

We synthesized the phosphinate 7 via photoaddition of methanol to the alpha, beta unsaturated deoxyribono lactone as the key step, followed by an Arbusov reaction for the introduction of phosphorous. Precursor 7 serves for the synthesis and incorporation into DNA of a novel chemically stable abasic site analogue that might act as an inhibitor for DNA glycosylases

Synthesis of pyrrolidine C-nucleosides via Heck reaction

A novel method for the synthesis of pyrrolidine C-nucleosides has been developed. The key step of the synthesis is the palladium(0)-mediated coupling of a disubstituted N-protected 2-pyrroline and 5-iodouracil. C-Nucleoside 14 and its N-methyl derivative 15 can easily be converted to the corresponding phosphoramidite building blocks for DNA synthesis

Chemistry of pyrrocorphins: methylative opening of the macrocycle between rings A and D, a side reaction in the peripheral C-methylation of a 20-methyl-pyrrocorphinate

One of the minor products from the previously described peripheral -methylation of a magnesium()-20-methyl--pyrrocorphinate is a C-19-methylated 19,20-seco-corphinoid derivative which, on complexation with nickel() acetate, recyclizes to a nickel()-tetradehydro-corrinate.