Applying reactive models to column experiments to assess the hydrogeochemistry of seawater intrusion: optimising ACUAINTRUSION and selecting cation exchange coefficients with PHREEQC

Three sets of laboratory column experimental results concerning the hydrogeochemistry of seawater intrusion have been modelled using two codes: ACUAINTRUSION (Chemical Engineering Department, University of Alicante) and PHREEQC (U.S.G.S.). These reactive models utilise the hydrodynamic parameters determined using the ACUAINTRUSION TRANSPORT software and fit the chloride breakthrough curves perfectly. The ACUAINTRUSION code was improved, and the instabilities were studied relative to the discretisation. The relative square errors were obtained using different combinations of the spatial and temporal steps: the global error for the total experimental data and the partial error for each element. Good simulations for the three experiments were obtained using the ACUAINTRUSION software with slight variations in the selectivity coefficients for both sediments determined in batch experiments with fresh water. The cation exchange parameters included in ACUAINTRUSION are those reported by the Gapon convention with modified exponents for the Ca/Mg exchange. PHREEQC simulations performed using the Gains-Thomas convention were unsatisfactory, with the exchange coefficients from the database of PHREEQC (or range), but those determined with fresh water – natural sediment allowed only an approximation to be obtained. For the treated sediment, the adjusted exchange coefficients were determined to improve the simulation and are vastly different from those from the database of PHREEQC or batch experiment values; however, these values fall in an order similar to the others determined under dynamic conditions. Different cation concentrations were simulated using two different software packages; this disparity could be attributed to the defined selectivity coefficients that affect the gypsum equilibrium. Consequently, different calculated sulphate concentrations are obtained using each type of software; a smaller mismatch was predicted using ACUAINTRUSION. In general, the presented simulations by ACUAINTRUSION and PHREEQC produced similar results, making predictions consistent with the experimental data. However, the simulated results are not identical to the experimental data; sulphate (total S) is overpredicted by both models, most likely due to such factors as the kinetics of gypsum, the possible variations in the exchange coefficients due to salinity and the neglect of other processes.

Thermodynamic calculations for the salt crystallisation damage in porous built heritage using PHREEQC

This work considers the crystallisation mechanisms of the most common and aggressive salts that generate stress in porous building stones as a result of changing ambient conditions. These mechanisms include the salt crystallisation that result from decreasing relative humidity and changes in temperature and, in hydrated salts, the dissolution of the lower hydrated form and the subsequent precipitation of the hydrated salt. We propose a new methodology for thermodynamic calculations using PHREEQC that includes these crystallisation mechanisms. This approach permits the calculation of the equilibrium relative humidity and the parameterization of the critical relative humidity and crystallisation pressures for the dissolution–precipitation transitions. The influence of other salts on the effectives of salt crystallisation and chemical weathering is also assessed. We review the sodium and magnesium sulphate and sodium chloride systems, in both single and multicomponent solutions, and they are compared to the sodium carbonate and calcium carbonate systems. The variation of crystallisation pressure, the formation of new minerals and the chemical dissolution by the presence of other salts is also evaluated. Results for hydrated salt systems show that high crystallisation pressures are possible as lower hydrated salts dissolve and more hydrated salts precipitate. High stresses may be also produced by decreasing temperature, although it requires that porous materials are wet for long periods of time. The presence of other salts changes the temperature and relative humidity of salt transitions that generates stress rather than reducing the pressure of crystallisation, if any salt has previously precipitated. Several practical conclusions derive from proposed methodology and provide conservators and architects with information on the potential weathering activity of soluble salts. Furthermore, the model calculations might be coupled with projections of future climate to give as improved understanding of the likely changes in the frequency of phase transitions in salts within porous stone.

MINERALOGY AND GENESIS OF SMECTITES IN AN ALKALINE-SALINE ENVIRONMENT OF PANTANAL WETLAND, BRAZIL

Smectite formation in alkaline-saline environments has been attributed to direct precipitation from solution and/or transformation from precursor minerals, but these mechanisms are not universally agreed upon in the literature. The objective of this work was to investigate the mineralogy of smectites in the soils surrounding a representative alkaline-saline lake of Nhecolandia, a sub-region of the Pantanal wetland, Brazil, and then to identify the mechanisms of their formation. Soils were sampled along a toposequence and analyzed by X-ray diffraction, transmission electron microscopy-energy dispersive X-ray analysis, and inductively coupled plasma-mass spectrometry. Water was collected along a transect involving the studied toposequence and equilibrium diagrams were calculated using the databases PHREEQC and AQUA. The fine-clay fraction is dominated by smectite, mica, and kaolinite. Smectites are concentrated at two places in the toposequence: an upper zone, which includes the soil horizons rarely reached by the lake-level variation; and a lower zone, which includes the surface horizon within the area of seasonal lake-level variation. Within the upper zone, the smectite is dioctahedral, rich in Al and Fe, and is classified as ferribeidellite. This phase is interstratified with mica and vermiculite and has an Fe content similar to that of the mica identified. These characteristics suggest that the ferribeidellite originates from transformation of micas and that vermiculite is an intermediate phase in this transformation. Within the lower zone, smectites are dominantly trioctahedral, Mg-rich, and are saponitic and stevensitic minerals. In addition, samples enriched in these minerals have much smaller rare-earth element (REE) contents than other soil samples. The water chemistry shows a geochemical control of Mg and saturation with respect to Mg-smectites in the more saline waters. The REE contents, water chemistry, and the presence of Mg-smectite where maximum evaporation is expected, suggest that saponitic and stevensitic minerals originate by chemical precipitation from the water column of the alkaline-saline lake.

Three-dimensional model for multi-component reactive transport with variable density groundwater flow

PHWAT is a new model that couples a geochemical reaction model (PHREEQC-2) with a density-dependent groundwater flow and solute transport model (SEAWAT) using the split-operator approach. PHWAT was developed to simulate multi-component reactive transport in variable density groundwater flow. Fluid density in PHWAT depends not on only the concentration of a single species as in SEAWAT, but also the concentrations of other dissolved chemicals that can be subject to reactive processes. Simulation results of PHWAT and PHREEQC-2 were compared in their predictions of effluent concentration from a column experiment. Both models produced identical results, showing that PHWAT has correctly coupled the sub-packages. PHWAT was then applied to the simulation of a tank experiment in which seawater intrusion was accompanied by cation exchange. The density dependence of the intrusion and the snow-plough effect in the breakthrough curves were reflected in the model simulations, which were in good agreement with the measured breakthrough data. Comparison simulations that, in turn, excluded density effects and reactions allowed us to quantify the marked effect of ignoring these processes. Next, we explored numerical issues involved in the practical application of PHWAT using the example of a dense plume flowing into a tank containing fresh water. It was shown that PHWAT could model physically unstable flow and that numerical instabilities were suppressed. Physical instability developed in the model in accordance with the increase of the modified Rayleigh number for density-dependent flow, in agreement with previous research. (c) 2004 Elsevier Ltd. All rights reserved.

Heavy metals removal in dual media filters

Dissertação para obtenção do Grau de Mestre em Engenharia do Ambiente – Perfil de Engenharia Sanitária

Geochemistry and water quality of Lake Qarun, Egypt

Water geochemistry is a very important tool for studying the water quality in a given area. Geology and climate are the major natural factors controlling the chemistry of most natural waters. Anthropogenic impacts are the secondary sources of contamination in natural waters. This study presents the first integrative approach to the geochemistry and water quality of surface waters and Lake Qarun in the Fayoum catchment, Egypt. Moreover, geochemical modeling of Lake Qarun was firstly presented. The Nile River is the main source of water to the Fayoum watershed. To investigate the quality and geochemistry of this water, water samples from irrigation canals, drains and Lake Qarun were collected during the period 2010‒2013 from the whole Fayoum drainage basin to address the major processes and factors governing the evolution of water chemistry in the investigation area. About 34 physicochemical quality parameters, including major ions, oxygen isotopes, trace elements, nutrients and microbiological parameters were investigated in the water samples. Multivariable statistical analysis was used to interpret the interrelationship between the different studied parameters. Geochemical modeling of Lake Qarun was carried out using Hardie and Eugster’s evolutionary model and a model simulated by PHREEQC software. The crystallization sequence during evaporation of Lake Qarun brine was also studied using a Jänecke phase diagram involving the system Na‒K‒Mg‒ Cl‒SO4‒H2O. The results show that the chemistry of surface water in the Fayoum catchment evolves from Ca- Mg-HCO3 at the head waters to Ca‒Mg‒Cl‒SO4 and eventually to Na‒Cl downstream and at Lake Qarun. The main processes behind the high levels of Na, SO4 and Cl in downstream waters and in Lake Qarun are dissolution of evaporites from Fayoum soils followed by evapoconcentration. This was confirmed by binary plots between the different ions, Piper plot, Gibb’s plot and δ18O results. The modeled data proved that Lake Qarun brine evolves from drainage waters via an evaporation‒crystallization process. Through the precipitation of calcite and gypsum, the solution should reach the final composition "Na–Mg–SO4–Cl". As simulated by PHREEQC, further evaporation of lake brine can drive halite to precipitate in the final stages of evaporation. Significantly, the crystallization sequence during evaporation of the lake brine at the concentration ponds of the Egyptian Salts and Minerals Company (EMISAL) reflected the findings from both Hardie and Eugster’s evolutionary model and the PHREEQC simulated model. After crystallization of halite at the EMISAL ponds, the crystallization sequence during evaporation of the residual brine (bittern) was investigated using a Jänecke phase diagram at 35 °C. This diagram was more useful than PHREEQC for predicting the evaporation path especially in the case of this highly concentrated brine (bittern). The predicted crystallization path using a Jänecke phase diagram at 35 °C showed that halite, hexahydrite, kainite and kieserite should appear during bittern evaporation. Yet the actual crystallized mineral salts were only halite and hexahydrite. The absence of kainite was due to its metastability while the absence of kieserite was due to opposed relative humidity. The presence of a specific MgSO4.nH2O phase in ancient evaporite deposits can be used as a paleoclimatic indicator. Evaluation of surface water quality for agricultural purposes shows that some irrigation waters and all drainage waters have high salinities and therefore cannot be used for irrigation. Waters from irrigation canals used as a drinking water supply show higher concentrations of Al and suffer from high levels of total coliform (TC), fecal coliform (FC) and fecal streptococcus (FS). These waters cannot be used for drinking or agricultural purposes without treatment, because of their high health risk. Therefore it is crucial that environmental protection agencies and the media increase public awareness of this issue, especially in rural areas.

Migração de K+, Cu2+ e CI- em solos do subgrupo Itararé no Estado de São Paulo: ensaios sorcivos, percolação em colunas e simulações geoquímicas

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Effects of carbonation on the mineral composition of cement kiln dust

Due to their relatively high calcium oxide content, industrial mineral oxide wastes are potential candidates for mineral sequestration of carbon dioxide (CO2). Cement kiln dust (CKD), a byproduct of cement manufacturing contains 20-60% CaO making it a possible candidate for CO2 sequestration. In this study, three types of CKD are characterized, before and after carbonation, using environmental scanning electron microscopy and energy dispersive x-ray microanalysis to determine the mineralogical and morphological changes occurring due to carbonation. The reactants, products, and precipitation mechanisms were investigated to enhance understanding of the governing processes and allow better utilization of CKD for CO2 sequestration. The results of multiple independent analyses confirmed the formation of CaCO3 during carbonation. Examinations of the reaction pathways found that CaO and calcium hydroxide (Ca(OH)2) were the major reactants. Three types of CaCO3 precipitation mechanisms were observed: (1) diffusion of CO2 into Ca(OH)2 particles causing precipitation in the pores of the particle and the growth of a CaCO3 ring from the outside inward, (2) precipitation onto existing particles, and (3) precipitation from aqueous solution. The growth of a CaCO3 ring on the outside of a particle may slow further diffusion of CO2 into a particle slowing iv the overall sequestration rate. Additionally, changes caused by carbonation in the solubility of trace metals were studied by mixing pre- and post-carbonated CKD with water and analyzing the solution using inductively coupled plasma mass spectrometry. Decreases in the leaching of chromium, lead, and copper were observed, and is an incentive for use of CKD for CO2 sequestration. Equilibrium modeling using PHREEQC confirmed that CaO and Ca(OH)2 would carbonate readily and form CaCO3.

(Supplement) Groundwater quality in Lahore Pakistan: emphasis on Arsenic and Fluoride levels

The dataset provides detailed information on the study that was conducted in Lahore's 7 major towns. The sample was taken from 472 tubewells and analyzed for major cations and anions using APHA 2012 techniques as explained herein. Besides, E.coli determination was done to check for microbial contamination. The data includes results from PHREEQC modeling of As(III)/ As(V) species and saturation indices as well as Aquachem's computed hydrochemical water facies. The WHO (2011) and EPA standards included in Aquachem identified the parameters that where in violation. Bicarbonates dominated the groundwater types with 50.21% of the samples exceeding the EPA maximum permissible limit of 250 mg/L in drinking water. Similarly, 30.51% of the samples had TDS values greater than 500 mg/L while 85.38 % of the samples exceed 10 µg/L threshold limit value of arsenic. Also, instances of high magnesium hazard values were observed which requires constant assessment if the groundwater is used for irrigation. Higher than 50% MH values are detrimental to crops which may reduce the expected yields. The membrane filtration technique using m-Endo Agar indicated that 3.59% samples had TNC (too numerous to count) values for E.coli while 5.06% showed values higher than 0 cfu/ 100 ml acceptable value in drinking water. Any traces of E-coli in a groundwater sample indicate recent fecal contamination. Such outcomes signify presence of enteric pathogens. If the groundwater is not properly dosed with disinfectants it may cause harm to human health. It is concluded that more studies are needed and proper groundwater management implement to safeguard the lives of communities that depend solely on groundwater in the city.

Modelling the fate of chromated copper arsenate in a sandy soil

A pulse of chromated copper arsenate (CCA, a timber preservative) was applied in irrigation water to an undisturbed field soil in a laboratory column. Concentrations of various elements in the leachate from the column were measured during the experiment. Also, the remnants within the soil were measured at the end of the experiment. The geochemical modelling package, PHREEQC-2, was used to simulate the experimental data. Processes included in the CCA transport modelling were advection, dispersion, non-specific adsorption (cation exchange) and specific adsorption by clay minerals and organic matter, as well as other possible chemical reactions such as precipitation/dissolution. The modelling effort highlighted the possible complexities in CCA transport and reaction experiments. For example, the uneven dosing of CCA as well as incomplete knowledge of the soil properties resulted in simulations that gave only partial, although reasonable, agreement with the experimental data. Both the experimental data and simulations show that As and Cu are strongly adsorbed and therefore, will mostly remain at the top of the soil profile, with a small proportion appearing in leachate. On the other hand, Cr is more mobile and thus it is present in the soil column leachate. Further simulations show that both the quantity of CCA added to the soil and the pH of the irrigation water will influence CCA transport. Simulations suggest that application of larger doses of CCA to the soil will result in higher leachate concentrations, especially for Cu and As. Irrigation water with a lower pH will dramatically increase leaching of Cu. These results indicate that acidic rainfall or significant accidental spillage of CCA will increase the risk of groundwater pollution.

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.

Study of the effectiveness of crystal growth modifiers in the prevention of damage due to crystallization of sodium carbonates in stone artworks

The disintegration of stone materials used in sculpture and architecture due to the crystallization of salts is capable of irreparably damaging artistic objects and historic buildings. A number of phosphonates and carboxylates were tested here as potential crystallization modifiers for sodium carbonate crystallization. Precipitated phases during crystallization induced either by cooling or by evaporation tests were nahcolite (NaHCO3), natron (Na2CO3∙10H2O) and thermonatrite (Na2CO3∙H2O), identified using X-ray diffraction. By using the thermodynamic code PHREEQC and the calculation of the nucleation rate it was demonstrated that nahcolite had to be first phase formed during both tests. The formation of the other phases depended on the experimental conditions under which the two tests were conducted. Nahcolite nucleation is strongly inhibited in the presence of sodium citrate tribasic dihydrate (CA), polyacrylic acid 2100MW (PA) and etidronic acid (HEDP), when the additives are dosed at appropriate concentrations and the pH range of the resulting solution is about 8. Electrostatic attraction generated between the deprotonated organic additives and the cations present in solution appears to be the principal mechanism of additive-nahcolite interaction. Salt weathering tests, in addition to mercury intrusion porosimetry tests allowed to quantify the damage induced by such salts. FESEM observation of both salts grown on calcite single crystals and in limestone blocks subjected to salt crystallization tests allowed to identify the effect of these additives on crystal growth and development. The results show that PA seems to be the best inhibitor, while CA and HEDP, which show similar behaviors, are slightly less effective. The use of such effective crystallization inhibitors may lead to more efficient preventive conservation of ornamental stone affected by crystallization damage due to formation of sodium carbonate crystals.