Development of an enhanced hydro-geochemical model to address mercury-speciation fate and transport in aquatic environments

An awareness of mercury (Hg) contamination in various aquatic environments around the world has increased over the past decade, mostly due to its ability to concentrate in the biota. Because the presence and distribution of Hg in aquatic systems depend on many factors (e.g., pe, pH, salinity, temperature, organic and inorganic ligands, sorbents, etc.), it is crucial to understand its fate and transport in the presence of complexing constituents and natural sorbents, under those different factors. An improved understanding of the subject will support the selection of monitoring, remediation, and restoration technologies. The coupling of equilibrium chemical reactions with transport processes in the model PHREEQC offers an advantage in simulating and predicting the fate and transport of aqueous chemical species of interest. Thus, a great variety of reactive transport problems could be addressed in aquatic systems with boundary conditions of specific interest. Nevertheless, PHREEQC lacks a comprehensive thermodynamic database for Hg. Therefore, in order to use PHREEQC to address the fate and transport of Hg in aquatic environments, it is necessary to expand its thermodynamic database, confirm it and then evaluate it in applications where potential exists for its calibration and continued validation. The objectives of this study were twofold: 1) to develop, expand, and confirm the Hg database of the hydrogeochemical PHREEQC to enhance its capability to simulate the fate of Hg species in the presence of complexing constituents and natural sorbents under different conditions of pH, redox, salinity and temperature; and 2) to apply and evaluate the new database in flow and transport scenarios, at two field test beds: Oak Ridge Reservation, Oak Ridge, TN and Everglades National Park, FL, where Hg is present and is of much concern. Overall, this research enhanced the capability of the PHREEQC model to simulate the coupling of the Hg reactions in transport conditions. It also demonstrated its usefulness when applied to field situations.

Chemical composition and stable sulfur isotope record of Black Sea sediments

The main terminal processes of organic matter mineralization in anoxic Black Sea sediments underlying the sulfidic water column are sulfate reduction in the upper 2-4 m and methanogenesis below the sulfate zone. The modern marine deposits comprise a ca. 1-m-deep layer of coccolith ooze and underlying sapropel, below which sea water ions penetrate deep down into the limnic Pleistocene deposits from >9000 years BP. Sulfate reduction rates have a subsurface maximum at the SO4[2-]-CH4 transition where H2S reaches maximum concentration. Because of an excess of reactive iron in the deep limnic deposits, most of the methane-derived H2S is drawn downward to a sulfidization front where it reacts with Fe(III) and with Fe2+ diffusing up from below. The H2S-Fe2+ transition is marked by a black band of amorphous iron sulfide above which distinct horizons of greigite and pyrite formation occur. The pore water gradients respond dynamically to environmental changes in the Black Sea with relatively short time constants of ca. 500 yr for SO4[2-] and 10 yr for H2S, whereas the FeS in the black band has taken ca. 3000 yr to accumulate. The dual diffusion interfaces of SO4[2-]-CH4 and H2S-Fe2+ cause the trapping of isotopically heavy iron sulfide with delta34S = +15 to +33 per mil at the sulfidization front. A diffusion model for sulfur isotopes shows that the SO4[2-] diffusing downward into the SO4[2-]-CH4 transition has an isotopic composition of +19 per mil, close to the +23 per mil of H2S diffusing upward. These isotopic compositions are, however, very different from the porewater SO4[2-] (+43 per mil) and H2S (-15 per mil) at the same depth. The model explains how methane-driven sulfate reduction combined with a deep H2S sink leads to isotopically heavy pyrite in a sediment open to diffusion. These results have general implications for the marine sulfur cycle and for the interpretation of sulfur isotopic data in modern sediments and in sedimentary rocks throughout earth's history.

Multi-impurity adsorption model for modeling crystal purity and shape evolution during crystallization processes in impure media

Peer reviewed

A new method to calculate the periodic steady state of sequencing batch reactors for biological wastewater treatment : model development and applications

Peer reviewed

Chlorine Contribution to Quantitative Structure and Activity Relationship Models of Disinfection By-Products' Quantum Chemical Descriptors and Toxicities

Quantitative Structure-Activity Relationship (QSAR) has been applied extensively in predicting toxicity of Disinfection By-Products (DBPs) in drinking water. Among many toxicological properties, acute and chronic toxicities of DBPs have been widely used in health risk assessment of DBPs. These toxicities are correlated with molecular properties, which are usually correlated with molecular descriptors. The primary goals of this thesis are: 1) to investigate the effects of molecular descriptors (e.g., chlorine number) on molecular properties such as energy of the lowest unoccupied molecular orbital (ELUMO) via QSAR modelling and analysis; 2) to validate the models by using internal and external cross-validation techniques; 3) to quantify the model uncertainties through Taylor and Monte Carlo Simulation. One of the very important ways to predict molecular properties such as ELUMO is using QSAR analysis. In this study, number of chlorine (NCl) and number of carbon (NC) as well as energy of the highest occupied molecular orbital (EHOMO) are used as molecular descriptors. There are typically three approaches used in QSAR model development: 1) Linear or Multi-linear Regression (MLR); 2) Partial Least Squares (PLS); and 3) Principle Component Regression (PCR). In QSAR analysis, a very critical step is model validation after QSAR models are established and before applying them to toxicity prediction. The DBPs to be studied include five chemical classes: chlorinated alkanes, alkenes, and aromatics. In addition, validated QSARs are developed to describe the toxicity of selected groups (i.e., chloro-alkane and aromatic compounds with a nitro- or cyano group) of DBP chemicals to three types of organisms (e.g., Fish, T. pyriformis, and P.pyosphoreum) based on experimental toxicity data from the literature. The results show that: 1) QSAR models to predict molecular property built by MLR, PLS or PCR can be used either to select valid data points or to eliminate outliers; 2) The Leave-One-Out Cross-Validation procedure by itself is not enough to give a reliable representation of the predictive ability of the QSAR models, however, Leave-Many-Out/K-fold cross-validation and external validation can be applied together to achieve more reliable results; 3) ELUMO are shown to correlate highly with the NCl for several classes of DBPs; and 4) According to uncertainty analysis using Taylor method, the uncertainty of QSAR models is contributed mostly from NCl for all DBP classes.

(Table S1) Tie-points used in age model construction for ODP Leg 172 sites

ix Ocean Drilling Program (ODP) sites, in the Northwest Atlantic have been used to investigate kinematic and chemical changes in the "Western Boundary Undercurrent" (WBUC) during the development of full glacial conditions across the Marine Isotope Stage 5a/4 boundary (~70,000 years ago). Sortable silt mean grain size(sort s) measurements are employed to examine changes in near bottom flow speeds, together with carbon isotopes measured in benthic foraminifera and % planktic foraminiferal fragmentation as proxies for changes in water-mass chemistry. A depth transect of cores, spanning 1.8-4.6 km depth, allows changes in both the strength and depth of the WBUC to be constrained across millennial scale events. Sort s measurements reveal that the flow speed structure of the WBUC during warm intervals ("interstadials") was comparable to modern (Holocene) conditions. However, significant differences are observed during cold intervals, with higher relative flow speeds inferred for the shallow component of the WBUC (~2 km depth) during all cold "stadial" intervals (including Heinrich Stadial 6), and a substantial weakening of the deep component (~3-4 km) during full glacial conditions. Our results therefore reveal that the onset of full glacial conditions was associated with a regime shift to a shallower mode of circulation (involving Glacial North Atlantic Intermediate Water) that was quantitatively distinct from preceding cold stadial events. Furthermore, our chemical proxy data show that the physical response of the WBUC during the last glacial inception was probably coupled to basin-wide changes in the water-mass composition of the deep Northwest Atlantic.

Age model and chemistry of sediment core STK7 from Steisslinger See

Stable isotope records for carbon and oxygen in bulk carbonates, carbon in bulk organic matter, and for total and chromium-reducible sulfur in a lacustrine sediment core from Lake Steisslingen (Southwest Germany) show several distinct and abrupt shifts during the last 15,000 years. Variations in the isotopic composition of authigenic carbonates indicate two major phases in the lake history. In the pre-Holocene, the hydrological budget of the lake was apparently stable. Variations of delta18O values of authigenic carbonates were, therefore, dominantly controlled by temperature changes. A decrease in the delta18Ocarb values of about 2 per mil at the Allerød/Younger Dryas transition is interpreted as a drop in mean annual air temperatures of approximately 5°C. An abrupt temperature increase of a similar magnitude is inferred at the Younger Dryas/Preboreal boundary. Throughout most of the Holocene, the isotopic composition of authigenic carbonates was influenced by marked changes in the hydrological budget of the lake. A major positive excursion in the delta13Ccarb and delta18Ocarb values at the beginning of the Atlantic and a smaller one in the Preboreal were related to evaporation effects, which indicate that dry climatic conditions must have prevailed at that time. A simultaneous increase in delta13C values of bulk organic matter at the beginning of the Atlantic suggests a high level of productivity in the lake. As a consequence, aqueous sulfate became limited as indicated by variations in the delta34S values of total and chromium-reducible sedimentary sulfur. Therefore, we conclude that the beginning of the Atlantic was characterized not only by dry but also by warm climatic conditions, which triggered a higher productivity in the lake. In the Subatlantic sediments, large variations in carbon, oxygen, and sulfur isotope ratios were observed as a result of human activities, causing considerable perturbations in the biogeochemical element cycling of Lake Steisslingen. Results obtained by the study of the continuous 15 ka record of Lake Steisslingen document clearly that isotopic proxy data from lacustrine sediments can provide useful information on environmental and climatic changes of local, regional, and in the case of the Younger Dryas event, of even hemispherical significance.

(Table 1) Chemical and isotopic composition of pore fluids in cores TGC-1 and TGC-6

The role of sediment diagenesis in the marine cycles of Li and B is poorly understood. Because Li and B are easily mobilized during burial and are consumed in authigenic clay mineral formation, their abundance in marine pore waters varies considerably. Exchange with the overlying ocean through diffusive fluxes should thus be common. Nevertheless, only a minor Li sink associated with the low-temperature alteration of volcanic ash has been observed. We describe a low-temperature diagenetic environment in the Black Sea dominated by the alteration of detrital plagioclase feldspars. Fluids expelled from the Odessa mud volcano in the Sorokin Trough originate from shallow (~100-400 m deep) sediments which are poor in volcanic materials but rich in anorthite. These fluids are depleted in Na+, K+, Li+, B, and 18O and enriched in Ca2+ and Sr2+, indicating that anorthite is dissolving and authigenic clays are forming. Using a simple chemical model, we calculate the pH and the partial pressure of CO2 (PCO2) in fluids associated with this alteration process. Our results show that the pH of these fluids is up to 1.5 pH units lower than in most deep marine sediments and that PCO2 levels are up to several hundred times higher than in the atmosphere. These conditions are similar to those which favor the weathering of silicate minerals in subaerial soil environments. We propose that in Black Sea sediments enhanced organic matter preservation favors CO2 production through methanogenesis and results in a low pore water pH, compared to most deep sea sediments. As a result, silicate mineral weathering, which is a sluggish process in most marine diagenetic environments, proceeds rapidly in Black Sea sediments. There is a potential for organic matter-rich continental shelf environments to host this type of diagenesis. Should such environments be widespread, this new Li and B sink could help balance the marine Li and Li isotope budgets but would imply an apparent imbalance in the B cycle.

New Method for the Estimation of Viscosity of Pure and Mixtures of Ionic Liquids Based on the UNIFAC–VISCO Model

A modified UNIFAC–VISCO group contribution method was developed for the correlation and prediction of viscosity of ionic liquids as a function of temperature at 0.1 MPa. In this original approach, cations and anions were regarded as peculiar molecular groups. The significance of this approach comes from the ability to calculate the viscosity of mixtures of ionic liquids as well as pure ionic liquids. Binary interaction parameters for selected cations and anions were determined by fitting the experimental viscosity data available in literature for selected ionic liquids. The temperature dependence on the viscosity of the cations and anions were fitted to a Vogel–Fulcher–Tamman behavior. Binary interaction parameters and VFT type fitting parameters were then used to determine the viscosity of pure and mixtures of ionic liquids with different combinations of cations and anions to ensure the validity of the prediction method. Consequently, the viscosities of binary ionic liquid mixtures were then calculated by using this prediction method. In this work, the viscosity data of pure ionic liquids and of binary mixtures of ionic liquids are successfully calculated from 293.15 K to 363.15 K at 0.1 MPa. All calculated viscosity data showed excellent agreement with experimental data with a relative absolute average deviation lower than 1.7%.

LES of flame acceleration and DDT in hydrogen-air mixture using artificially thickened flame approach and detailed chemical kinetics

A large eddy simulation is performed to study the deflagration to detonation transition phenomenon in an obstructed channel containing premixed stoichiometric hydrogen–air mixture. Two-dimensional filtered reactive Navier–Stokes equations are solved utilizing the artificially thickened flame approach (ATF) for modeling sub-grid scale combustion. To include the effect of induction time, a 27-step detailed mechanism is utilized along with an in situ adaptive tabulation (ISAT) method to reduce the computational cost due to the detailed chemistry. The results show that in the slow flame propagation regime, the flame–vortex interaction and the resulting flame folding and wrinkling are the main mechanisms for the increase of the flame surface and consequently acceleration of the flame. Furthermore, at high speed, the major mechanisms responsible for flame propagation are repeated reflected shock–flame interactions and the resulting baroclinic vorticity. These interactions intensify the rate of heat release and maintain the turbulence and flame speed at high level. During the flame acceleration, it is seen that the turbulent flame enters the ‘thickened reaction zones’ regime. Therefore, it is necessary to utilize the chemistry based combustion model with detailed chemical kinetics to properly capture the salient features of the fast deflagration propagation.

Simplified model for prediction of compressive strength of alkali-activated natural pozzolans

The assessment of pozzolanic activity is essential for estimating the reaction of a material as pozzolan. Natural pozzolans can be activated and condensed with sodium silicate in an alkaline environment to synthesize high performance cementitious construction materials with low environmental impact. In this paper, the pozzolanic activities of five natural pozzolans are studied. The correlation between type and chemical composition of natural pozzolan, which affects the formation of the geopolymer gel phase, both for the calcined and untreated natural pozzolans, have been reviewed. The improvement in pozzolanic properties was studied following heat treatment including calcinations and/or elevated curing temperature by using alkali solubility, and compressive strength tests. A model was developed to allow prediction of the alkali-activated pozzolan strength versus their chemical compositions, alkali solubility, and crystallinity.

New Method Based on the UNIFAC–VISCO Model for the Estimation of Ionic Liquids Viscosity Using the Experimental Data Recommended by Mathematical Gnostics

The viscosity of ionic liquids (ILs) has been modeled as a function of temperature and at atmospheric pressure using a new method based on the UNIFAC–VISCO method. This model extends the calculations previously reported by our group (see Zhao et al. J. Chem. Eng. Data 2016, 61, 2160–2169) which used 154 experimental viscosity data points of 25 ionic liquids for regression of a set of binary interaction parameters and ion Vogel–Fulcher–Tammann (VFT) parameters. Discrepancies in the experimental data of the same IL affect the quality of the correlation and thus the development of the predictive method. In this work, mathematical gnostics was used to analyze the experimental data from different sources and recommend one set of reliable data for each IL. These recommended data (totally 819 data points) for 70 ILs were correlated using this model to obtain an extended set of binary interaction parameters and ion VFT parameters, with a regression accuracy of 1.4%. In addition, 966 experimental viscosity data points for 11 binary mixtures of ILs were collected from literature to establish this model. All the binary data consist of 128 training data points used for the optimization of binary interaction parameters and 838 test data points used for the comparison of the pure evaluated values. The relative average absolute deviation (RAAD) for training and test is 2.9% and 3.9%, respectively.

Theoretical Performance Model and Initial Experimentation of a Baffled-Tube Ram Accelerator

Thesis (Master's)--University of Washington, 2016-08

The three dimensional model with the affection of stratified density (Caspian Sea)

The hydro dynamical actions in big Lakes directly influence dynamic, physical and chemical affairs. The circulation's models and temperature have something to do with the movements of fluids, and analysis for circulation in Caspian sea is because of the lack of observation through which the circulations and out comings are determined. Through the studies, three dimensional simulations (Large- Scale) are planned and performed, according to Smolakiewicz and Margolin works. This is a non- hydrostatic and Boussinesq approximation is used in its formulation is used in its formulation on the basis of Lipps (1990) theorem and curve lines, the fluid is constant adiabatic and stratified, and the wind power is considered zero. The profile of speed according to previous depth and before ridge can be drawn on the basis of density available between northern and southern ridges. The circulation field is drawn from 3 cm/s to 13 cm/s on the plate z= 5 cm , the vertical changes of speed on the plate is 0.02 m/s. Vertical profile , horizontal speed in previous on, and after the ridges on are drawn on different spaces. It changes from 0.5 cm/s to 1 cm/s before ridges.