Geochemical Modeling and Analysis of the Frac Water Used in the Hydraulic Fracturing of the Marcellus Formation, Pennsylvania

The hydraulic fracturing of the Marcellus Formation creates a byproduct known as frac water. Five frac water samples were collected in Bradford County, PA. Inorganic chemical analysis, field parameters analysis, alkalinity titrations, total dissolved solids(TDS), total suspended solids (TSS), biological oxygen demand (BOD), and chemical oxygen demand (COD) were conducted on each sample to characterize frac water. A database of frac water chemistry results from across the state of Pennsylvania from multiple sources was compiled in order to provide the public and research communitywith an accurate characterization of frac water. Four geochemical models were created to model the reactions between frac water and the Marcellus Formation, Purcell Limestone, and the oil field brines presumed present in the formations. The average concentrations of chloride and TDS in the five frac water samples were 1.1 �± 0.5 x 105 mg/L (5.5X average seawater) and 140,000 mg/L (4X average seawater). BOD values for frac water immediately upon flow back were over 10X greater than the BOD of typical wastewater, but decreased into the range of typical wastewater after a short period of time. The COD of frac water decreases dramatically with an increase in elapsed time from flow back, but remain considerably higher than typicalwastewater. Different alkalinity calculation methods produced a range of alkalinity values for frac water: this result is most likely due to high concentrations of aliphatic acid anions present in the samples. Laboratory analyses indicate that the frac watercomposition is quite variable depending on the companies from which the water was collected, the geology of the local area, and number of fracturing jobs in which the frac water was used, but will require more treatment than typical wastewater regardless of theprecise composition of each sample. The geochemical models created suggest that the presence of organic complexes in an oil field brine and Marcellus Formation aid in the dissolution of ions such as bariumand strontium into the solution. Although equilibration reactions between the Marcellus Formation and the slickwater account for some of the final frac water composition, the predominant control of frac water composition appears to be the ratio of the mixture between the oil field brine and slickwater. The high concentration of barium in the frac water is likely due to the abundance of barite nodules in the Purcell Limestone, and the lack of sulfate in the frac water samples is due to the reducing, anoxic conditions in the earth's subsurface that allow for the degassing of H2S(g).

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.

Geochemical determination of the fate and transport of injected fresh wastewater to a deep saline aquifer

Deep well injection into non-potable saline aquifers of treated domestic wastewater has been used in Florida for decades as a safe and effective alternative to ocean outfall disposal. The objectives of this study were to determine the fate and transport of injected wastewater at two deep well injection sites in Miami Dade County, Florida, USA. Detection of ammonium in the Middle Confining units of the Floridan aquifer above the injection zone at both sites has been interpreted as evidence of upward migration of injected wastewater, posing a risk to underground sources of drinking water. Historical water quality data, including ammonia, chloride, temperature, and pH from existing monitoring wells at both sites from 1983 to 2008, major ions collected monthly from 2006 and 2008, and a synoptic sampling event for stable isotopes, tritium, and dissolved gases in 2008, were used to determine the source of ammonium in groundwater and possible migration pathways. Geochemical modeling was used to determine possible effects of injected wastewater on native water and aquifer matrix geochemistry. Injected wastewater was determined to be the source of elevated ammonium concentrations above ambient water levels, based on the results of major ion concentrations, tritium, dissolved noble gases and 15N isotopes analyses. Various possible fluid migration pathways were identified at the sites. Data for the south site suggest buoyancy-driven vertical pathways to overlying aquifers bypassing the confining units, with little mixing of injected wastewater with native water as it migrated upward. Once it is introduced into an aquifer, the injectate appeared to migrate advectively with the regional groundwater flow. Geochemical modeling indicated that CO 2-enriched injected wastewater allowed for carbonate dissolution along the vertical pathways, enhancing permeability along these flowpaths. At the north site, diffusive upward flow through the confining units or offsite vertical pathways were determined to be possible, however no evidence was detected for any on-site confining unit bypass pathway. No evidence was observed at either site of injected wastewater migration to the Upper Floridan aquifer, which is used as a municipal water supply and for aquifer storage and recovery.

Role of fluid mixing in deep dissolution of carbonates

The presence of cavities filled with new minerals in carbonate rocks is a common feature in oil reservoirs and lead-zinc deposits. Since groundwater equilibrates rapidly with carbonates, the presence of dissolution cavities in deep carbonate host rocks is a paradox. Two alternative geochemical processes have been proposed to dissolve carbonates at depth: hydrogen sulfide oxidation to sulfuric acid, and metal sulfide precipitation. With the aid of geochemical modeling we show that mixing two warm solutions saturated with carbonate results in a new solution that dissolves limestone. Variations in the proportion of the end-member fluids can also form a supersaturated mixture and fill the cavity with a new generation of carbonate. Mixing is in general more effective in dissolving carbonates than the aforementioned processes. Moreover, mixing is consistent with the wide set of textures and mineral proportions observed in cavity infillings.

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.

Geologia, Geocronologia, Geoquímica e Petrogênese das Rochas Ígneas Cretácicas da Província Magmática do Cabo e suas Relações com as Unidades Sedimentares da Bacia de Pernambuco (NE do Brasil)

The area studied forms a thin NNE-directed belt situated south of Recife town (Pernambuco state), northeastern Brazil. Geologically, it comprises the Pernambuco Basin (PB), which is limited by the Pernambuco Lineament to the north, the Maragogi high to the south and the Pernambuco Alagoas massif to the west, all of them with Precambrian age. This thesis reports the results obtained for the Cabo Magmatic Province (CMP), aiming the characterization of the geology, stratigraphy, geochronology, geochemistry and petrogenesis of the Cretaceous igneous rocks presented in the PB. The PB is composed of the Cabo Formation (rift phase) at the base (polymictic conglomerates, sandstones, shales), an intermediate unit, the Estiva Formation (marbles and argillites), and, at the top, the Algodoais Formation (monomictic conglomerates, sandstones, shales). The CMP is represented by trachytes, rhyolites, pyroclastics (ignimbrites), basalts / trachy-andesites, monzonites and alkali-feldspar granite, which occur as dykes, flows, sills, laccoliths and plugs. Field observations and well descriptions show that the majority of the magmatic rocks have intrusive contacts with the Cabo Formation, although some occurrences are also suggestive of synchronism between volcanism and siliciclastic sedimentation. 40Ar/39Ar and zircon fission tracks for the magmatic rocks indicate an average age of 102 r 1 Ma for the CMP. This age represents an expressive event in the province and is detected in all igneous dated materials. It is considered as a minimum age (Albian) for the magmatic episode and the peak of the rift phase in the PB. The 40Ar/39Ar dates are about 10-14 Ma younger than published palynologic ages for this basin. Geochemically, the CMP may be divided in two major groups; i) a transitional to alkaline suite, constituted by basalts to trachy-andesites (types with fine-grained textures and phenocrysts of sanidine and plagioclase), trachytes (porphyrytic texture, with phenocrysts of sanidine and plagioclase) and monzonites; ii) a alkaline suite, highly fractionated, acidic volcano-plutonic association, formed by four subtypes (pyroclastic flows ignimbrites, fine-to medium-grained rhyolites, a high level granite, and later rhyolites). These four types are distinguished essentially by field aspects and petrographic and textural features. Compatible versus incompatible trace element concentrations and geochemical modeling based on both major and trace elements suggest the evolution through low pressure fractional crystallization for trachytes and other acidic rocks, whereas basalts / trachy-andesites and monzonites evolved by partial melting from a mantle source. Sr and Nd isotopes reveal two distinct sources for the rocks of the CMP. Concerning the acidic ones, the high initial Sr ratios (ISr = 0.7064-1.2295) and the negative HNd (-0.43 to -3.67) indicate a crustal source with mesoproterozoic model ages (TDM from 0.92 to 1.04 Ga). On the other hand, the basic to intermediate rocks have low ISr (0.7031-0.7042) and positive HNd (+1.28 to +1.98), which requires the depleted mantle as the most probable source; their model ages are in the range 0.61-0.66 Ga. However, the light rare earth enrichment of these rocks and partial melting modeling point to an incompatible-enriched lherzolitic mantle with very low quantity of garnet (1-3%). This apparent difference between geochemical and Nd isotopes may be resolved by assuming that the metasomatizing agent did not obliterate the original isotopic characteristics of the magmas. A 2 to 5% partial melting of this mantle at approximately 14 kbar and 1269oC account very well the basalts and trachy-andesites studied. By using these pressure and temperatures estimates for the generation of the basaltic to trachy-andesitic magma, it is determined a lithospheric stretching (E) of 2.5. This E value is an appropriated estimate for the sub-crustal stretching (astenospheric or the base of the lithosphere?) region under the Pernambuco Basin, the crustal stretching probably being lower. The integration of all data obtained in this thesis permits to interpret the magmatic evolution of the PB as follows; 1st) the partial melting of a garnet-bearing lherzolite generates incompatible-enriched basaltic, trachy-andesitic and monzonitic magmas; 2nd) the underplating of these basaltic magmas at the base of the continental crust triggers the partial melting of this crust, and thus originating the acidic magmas; 3rd) concomitantly with the previous stage, trachytic magmas were produced by fractionation from a monzonitic to trachy-andesitic liquid; 4th) the emplacement of the several magmas in superficial (e.g. flows) or sub-superficial (e.g. dykes, sills, domes, laccoliths) depths was almost synchronically, at about 102 r 1 Ma, and usually crosscutting the sedimentary rocks of the Cabo Formation. The presence of garnet in the lherzolitic mantle does not agree with pressures of about 14 kbar for the generation of the basaltic magma, as calculated based on chemical parameters. This can be resolved by admitting the astenospheric uplifting under the rift, which would place deep and hot material (mantle plume?) at sub-crustal depths. The generation of the magmas and their subsequent emplacement would be coupled with the crustal rifting of the PB, the border (NNE-SSW directed) and transfer (NW-SE directed) faults serving as conduits for the magma emplacement. Based on the E parameter and the integration of 40Ar/39Ar and palynologic data it is interpreted a maximum duration of 10-14 Ma for the rift phase (Cabo Formation clastic sedimentation and basic to acidic magmatism) of the PB

Strontium concentrations and isotopic compositions of anhydrites from the TAG active mound

Sr isotope analyses have been conducted on anhydrite samples from the TAG (Trans-Atlantic Geotraverse) active hydrothermal mound (26°08?N, Mid-Atlantic Ridge) that have previously been shown to exhibit two distinct patterns of REE behavior when normalized to TAG end-member hydrothermal fluid. Despite differences in REE patterns, the Sr isotope data indicate that all the anhydrites precipitated from fluids with a similar range of hydrothermal fluid and seawater components, and all but one were seawater-dominated (52%-75%). Speciation calculations using the EQ3/6 software package for geochemical modeling of aqueous systems suggest that the REE complexation behavior in different fluid mixing scenarios can explain the variations in the REE patterns. Anhydrites that exhibit relatively flat REE patterns [(La_bs)/(Yb_bs) = 0.8-2.0; subscript bs indicates normalization to end-member black smoker hydrothermal fluid] and a small or no Eu anomaly [(Eu_bs)/(Eu*_bs) = 0.8-2.0] are inferred to have precipitated from mixes of end-member hydrothermal fluid and cold seawater. REE complexes with hard ligands (e.g., fluoride and chloride) are less stable at low temperatures and trivalent Eu has an ionic radius similar to that of Ca2+ and the other REE, and so they behave coherently. In contrast, anhydrites that exhibit slight LREE-depletion [(La_bs)/(Yb_bs) = 0.4-1.4] and a distinct negative anomaly [(Eu_bs)/(Eu*_bs) = 0.2-0.8] are inferred to have precipitated from mixes of end-member hydrothermal fluid and conductively heated seawater. The LREE depletion results from the presence of very stable LREE chloro-complexes that effectively limit the availability of the LREE for partitioning into anhydrite. Above 250°C, Eu is present only in divalent form as chloride complexes, and discrimination against Eu2+ is likely due to both the mismatch in ionic radii between Eu2+ and Ca2+, and the strong chloro-complexation of divalent Eu which promotes stability in the fluid and inhibits partitioning of Eu2+ into precipitating anhydrite. These variations in REE behavior attest to rapid fluctuations in thermal regime, fluid flow and mixing in the subsurface of the TAG mound that give rise to heterogeneity in the formation conditions of individual anhydrite crystals.

Negative pH, efflorescent mineralogy, and consequences for environmental restoration at the Iron Mountain Superfund site, California

The Richmond Mine of the Iron Mountain copper deposit contains some of the most acid mine waters ever reported. Values of pH have been measured as low as −3.6, combined metal concentrations as high as 200 g/liter, and sulfate concentrations as high as 760 g/liter. Copious quantities of soluble metal sulfate salts such as melanterite, chalcanthite, coquimbite, rhomboclase, voltaite, copiapite, and halotrichite have been identified, and some of these are forming from negative-pH mine waters. Geochemical calculations show that, under a mine-plugging remediation scenario, these salts would dissolve and the resultant 600,000-m3 mine pool would have a pH of 1 or less and contain several grams of dissolved metals per liter, much like the current portal effluent water. In the absence of plugging or other at-source control, current weathering rates indicate that the portal effluent will continue for approximately 3,000 years. Other remedial actions have greatly reduced metal loads into downstream drainages and the Sacramento River, primarily by capturing the major acidic discharges and routing them to a lime neutralization plant. Incorporation of geochemical modeling and mineralogical expertise into the decision-making process for remediation can save time, save money, and reduce the likelihood of deleterious consequences.

Life cycle of Deccan trap magma chambers: A crystal scale elemental and strontium isotopic investigation

The Deccan Trap basalts are the remnants of a massive series of lava flows that erupted at the K/T boundary and covered 1-2 million km2 of west-central India. This eruptive event is of global interest because of its possible link to the major mass extinction event, and there is much debate about the duration of this massive volcanic event. In contrast to isotopic or paleomagnetic dating methods, I explore an alternative approach to determine the lifecycle of the magma chambers that supplied the lavas, and extend the concept to obtain a tighter constraint on Deccan’s duration. My method relies on extracting time information from elemental and isotopic diffusion across zone boundaries in individual crystals. I determined elemental and Sr-isotopic variations across abnormally large (2-5 cm) plagioclase crystals from the Thalghat and Kashele “Giant Plagioclase Basalts” from the lowermost Jawhar and Igatpuri Formations respectively in the thickest Western Ghats section near Mumbai. I also obtained bulk rock major, trace and rare earth element chemistry of each lava flow from the two formations. Thalghat flows contain only 12% zoned crystals, with 87 Sr/86Sr ratios of 0.7096 in the core and 0.7106 in the rim, separated by a sharp boundary. In contrast, all Kashele crystals have a wider range of 87Sr/86Sr values, with multiple zones. Geochemical modeling of the data suggests that the two types of crystals grew in distinct magmatic environments. Modeling intracrystalline diffusive equilibration between the core and rim of Thalghat crystals led me to obtain a crystal growth rate of 2.03x10-10 cm/s and a residence time of 780 years for the crystals in the magma chamber(s). Employing some assumptions based on field and geochronologic evidence, I extrapolated this residence time to the entire Western Ghats and obtained an estimate of 25,000–35,000 years for the duration of Western Ghats volcanism. This gave an eruptive rate of 30–40 km3/yr, which is much higher than any presently erupting volcano. This result will remain speculative until a similarly detailed analytical-modeling study is performed for the rest of the Western Ghats formations.

Groundwater-surface water interactions on tree islands in the Everglades, south Florida

The marked decline in tree island cover across the Everglades over the last century, has been attributed to landscape-scale hydrologic degradation. To preserve and restore Everglades tree islands, a clear understanding of tree island groundwater-surface water interactions is needed, as these interactions strongly influence the chemistry of shallow groundwater and the location and patterns of vegetation in many wetlands. The goal of this work was to define the relationship between groundwater-surface water interactions, plant-water uptake, and the groundwater geochemical condition of tree islands. Groundwater and surface water levels, temperature, and chemistry were monitored on eight constructed and one natural tree island in the Everglades from 2007–2010. Sap flow, diurnal water table fluctuations and stable oxygen isotopes of stem, ground and soil water were used to determine the effect of plant-water uptake on groundwater-surface water interactions. Hydrologic and geochemical modeling was used to further explore the effect of plant-groundwater-surface water interactions on ion concentrations and potential mineral formation.^

Arsenic Geochemitry in Warm Spring Ponds: New Field and Experimental Results

Silver Bow Creek (SBC) flows into the Warm Springs Ponds Operable Unit (WSPOU), where various containment cells are used to precipitate copper and other metals (e.g., Cd, Cu, Mn, Pb, Zn). Lime is added seasonally to increase the pH and assist in removal of metals from the water column. Although the WSPOU is effective at removing copper and other cationic trace metals, concentrations of dissolved arsenic exiting the facility are often above the site specific standard, 20 20 ug/L, during low-flow periods each summer and fall. This thesis is a continuation of arsenic geochemistry studies by Montana Tech in the WSPOU. Field work focused on Pond 3, the largest and first in the series of treatment ponds. Shallow groundwater was sampled from 8 PVC piezometers located near the south end of Pond 3. Three sediment pore-water diffusion samplers (“peepers”) were also deployed at the south end of Pond 3 to examine vertical gradients in chemistry in the top 25 cm of the pond sediment. In general, the pH and Eh values of the shallow groundwater and sediment pore-water were less than in the pond water. Concentrations of arsenic were generally higher in subsurface water, and tended to pass through a maximum (up to 530 g/L) about 10 cm below the sediment-water interface. In the peeper cells, there was a strong positive correlation between dissolved As and dissolved Fe, and an inverse correlation with sulfate. Therefore, the zone of arsenic release corresponds to a zone of bacterial Fe and sulfate reduction in the shallow, organic-rich sediment. Redox speciation of arsenic shows that arsenate (As(V)) is dominant in the pond, and arsenite (As(III)) is dominant in the subsurface water. A series of laboratory experiments with pH adjustment were completed using SBC water collected near the inlet to the WSPOU as well as water and shallow sediment collected from Pond 3. Water ± sediment mesocosms were set up in 1-L Nalgene bottles (closed system) or a 20-L aquarium (open system), both with continuous stirring. The pH of the mesocosm was adjusted by addition of NaOH or HNO3 acid. The closed system provided better pH control since the water was not in contact with the atmosphere, which prevented exchange of carbon dioxide. In both the closed and open systems, dissolved arsenic concentrations either decreased or stayed roughly the same with increase in pH to values > 11. Therefore, the release of dissolved As into the treatment ponds in low-flow periods is not due to changes in pH alone. All of these results support the hypothesis that the arsenic release in WSPOU is linked to microbial reduction of ferric oxide minerals in the organic-rich sediment. Upwards diffusion of dissolved As from the sediment pore-water into the pond water is the most likely explanation for the increase in As concentration of the WSPOU in low-flow periods.

Life Cycle of Deccan Trap Magma Chambers: A Crystal Scale Elemental and Strontium Isotopic Investigation

The Deccan Trap basalts are the remnants of a massive series of lava flows that erupted at the K/T boundary and covered 1-2 million km2 of west-central India. This eruptive event is of global interest because of its possible link to the major mass extinction event, and there is much debate about the duration of this massive volcanic event. In contrast to isotopic or paleomagnetic dating methods, I explore an alternative approach to determine the lifecycle of the magma chambers that supplied the lavas, and extend the concept to obtain a tighter constraint on Deccan’s duration. My method relies on extracting time information from elemental and isotopic diffusion across zone boundary in an individual crystal. I determined elemental and Sr-isotopic variations across abnormally large (2-5 cm) plagioclase crystals from the Thalghat and Kashele “Giant Plagioclase Basalts” from the lowermost Jawhar and Igatpuri Formations respectively in the thickest Western Ghats section near Mumbai. I also obtained bulk rock major, trace and rare earth element chemistry of each lava flow from the two formations. Thalghat flows contain only 12% zoned crystals, with 87Sr/86Sr ratios of 0.7096 in the core and 0.7106 in the rim, separated by a sharp boundary. In contrast, all Kashele crystals have a wider range of 87Sr/86Sr values, with multiple zones. Geochemical modeling of the data suggests that the two types of crystals grew in distinct magmatic environments. Modeling intracrystalline diffusive equilibration between the core and rim of Thalghat crystals led me to obtain a crystal growth rate of 2.03x10-10 cm/s and a residence time of 780 years for the crystals in the magma chamber(s). Employing some assumptions based on field and geochronologic evidence, I extrapolated this residence time to the entire Western Ghats and obtained an estimate of 25,000 – 35,000 years for the duration of Western Ghats volcanism. This gave an eruptive rate of 30 – 40 km3/yr, which is much higher than any presently erupting volcano. This result will remain speculative until a similarly detailed analytical-modeling study is performed for the rest of the Western Ghats formations.

Geochemical characterization of heavy metal contaminated area using multivariate factorial kriging

This paper describes a geostatistical method, known as factorial kriging analysis, which is well suited for analyzing multivariate spatial information. The method involves multivariate variogram modeling, principal component analysis, and cokriging. It uses several separate correlation structures, each corresponding to a specific spatial scale, and yields a set of regionalized factors summarizing the main features of the data for each spatial scale. This method is applied to an area of high manganese-ore mining activity in Amapa State, North Brazil. Two scales of spatial variation (0.33 and 2.0 km) are identified and interpreted. The results indicate that, for the short-range structure, manganese, arsenic, iron, and cadmium are associated with human activities due to the mining work, while for the long-range structure, the high aluminum, selenium, copper, and lead concentrations, seem to be related to the natural environment. At each scale, the correlation structure is analyzed, and regionalized factors are estimated by cokriging and then mapped.

Modeling hydrothermal system: deriving observables and hydrothermal instability in volcanic and non-volcanic setting

Hydrothermal fluids are a fundamental resource for understanding and monitoring volcanic and non-volcanic systems. This thesis is focused on the study of hydrothermal system through numerical modeling with the geothermal simulator TOUGH2. Several simulations are presented, and geophysical and geochemical observables, arising from fluids circulation, are analyzed in detail throughout the thesis. In a volcanic setting, fluids feeding fumaroles and hot spring may play a key role in the hazard evaluation. The evolution of the fluids circulation is caused by a strong interaction between magmatic and hydrothermal systems. A simultaneous analysis of different geophysical and geochemical observables is a sound approach for interpreting monitored data and to infer a consistent conceptual model. Analyzed observables are ground displacement, gravity changes, electrical conductivity, amount, composition and temperature of the emitted gases at surface, and extent of degassing area. Results highlight the different temporal response of the considered observables, as well as the different radial pattern of variation. However, magnitude, temporal response and radial pattern of these signals depend not only on the evolution of fluid circulation, but a main role is played by the considered rock properties. Numerical simulations highlight differences that arise from the assumption of different permeabilities, for both homogeneous and heterogeneous systems. Rock properties affect hydrothermal fluid circulation, controlling both the range of variation and the temporal evolution of the observable signals. Low temperature fumaroles and low discharge rate may be affected by atmospheric conditions. Detailed parametric simulations were performed, aimed to understand the effects of system properties, such as permeability and gas reservoir overpressure, on diffuse degassing when air temperature and barometric pressure changes are applied to the ground surface. Hydrothermal circulation, however, is not only a characteristic of volcanic system. Hot fluids may be involved in several mankind problems, such as studies on geothermal engineering, nuclear waste propagation in porous medium, and Geological Carbon Sequestration (GCS). The current concept for large-scale GCS is the direct injection of supercritical carbon dioxide into deep geological formations which typically contain brine. Upward displacement of such brine from deep reservoirs driven by pressure increases resulting from carbon dioxide injection may occur through abandoned wells, permeable faults or permeable channels. Brine intrusion into aquifers may degrade groundwater resources. Numerical results show that pressure rise drives dense water up to the conduits, and does not necessarily result in continuous flow. Rather, overpressure leads to new hydrostatic equilibrium if fluids are initially density stratified. If warm and salty fluid does not cool passing through the conduit, an oscillatory solution is then possible. Parameter studies delineate steady-state (static) and oscillatory solutions.

Partial melting within the lower crust – Constraints from bulk rock geochemical data and trace element mineral analyses of granulites and migmatites from Central Finland

This PhD thesis concerns geochemical constraints on recycling and partial melting of Archean continental crust. A natural example of such processes was found in the Iisalmi area of Central Finland. The rocks from this area are Middle to Late Archean in age and experienced metamorphism and partial melting between 2.7-2.63 Ga. The work is based on extensive field work. It is furthermore founded on bulk rock geochemical data as well as in-situ analyses of minerals. All geochemical data were obtained at the Institute of Geosciences, University of Mainz using X-ray fluorescence, solution ICP-MS and laser ablation-ICP-MS for bulk rock geochemical analyses. Mineral analyses were accomplished by electron microprobe and laser ablation ICP-MS. Fluid inclusions were studied by microscope on a heating-freezing-stage at the Geoscience Center, University Göttingen. Part I focuses on the development of a new analytical method for bulk rock trace element determination by laser ablation-ICP-MS using homogeneous glasses fused from rock powder on an Iridium strip heater. This method is applicable for mafic rock samples whose melts have low viscosities and homogenize quickly at temperatures of ~1200°C. Highly viscous melts of felsic samples prevent melting and homogenization at comparable temperatures. Fusion of felsic samples can be enabled by addition of MgO to the rock powder and adjustment of melting temperature and melting duration to the rock composition. Advantages of the fusion method are low detection limits compared to XRF analyses and avoidance of wet-chemical processing and use of strong acids as in solution ICP-MS as well as smaller sample volumes compared to the other methods. Part II of the thesis uses bulk rock geochemical data and results from fluid inclusion studies for discrimination of melting processes observed in different rock types. Fluid inclusion studies demonstrate a major change in fluid composition from CO2-dominated fluids in granulites to aqueous fluids in TTG gneisses and amphibolites. Partial melts were generated in the dry, CO2-rich environment by dehydration melting reactions of amphibole which in addition to tonalitic melts produced the anhydrous mineral assemblages of granulites (grt + cpx + pl ± amph or opx + cpx + pl + amph). Trace element modeling showed that mafic granulites are residues of 10-30 % melt extraction from amphibolitic precursor rocks. The maximum degree of melting in intermediate granulites was ~10 % as inferred from modal abundances of amphibole, clinopyroxene and orthopyroxene. Carbonic inclusions are absent in upper-amphibolite facies migmatites whereas aqueous inclusion with up to 20 wt% NaCl are abundant. This suggests that melting within TTG gneisses and amphibolites took place in the presence of an aqueous fluid phase that enabled melting at the wet solidus at temperatures of 700-750°C. The strong disruption of pre-metamorphic structures in some outcrops suggests that the maximum amount of melt in TTG gneisses was ~25 vol%. The presence of leucosomes in all rock types is taken as the principle evidence for melt formation. However, mineralogical appearance as well as major and trace element composition of many leucosomes imply that leucosomes seldom represent frozen in-situ melts. They are better considered as remnants of the melt channel network, e.g. ways on which melts escaped from the system. Part III of the thesis describes how analyses of minerals from a specific rock type (granulite) can be used to determine partition coefficients between different minerals and between minerals and melt suitable for lower crustal conditions. The trace element analyses by laser ablation-ICP-MS show coherent distribution among the principal mineral phases independent of rock composition. REE contents in amphibole are about 3 times higher than REE contents in clinopyroxene from the same sample. This consistency has to be taken into consideration in models of lower crustal melting where amphibole is replaced by clinopyroxene in the course of melting. A lack of equilibrium is observed between matrix clinopyroxene / amphibole and garnet porphyroblasts which suggests a late stage growth of garnet and slow diffusion and equilibration of the REE during metamorphism. The data provide a first set of distribution coefficients of the transition metals (Sc, V, Cr, Ni) in the lower crust. In addition, analyses of ilmenite and apatite demonstrate the strong influence of accessory phases on trace element distribution. Apatite contains high amounts of REE and Sr while ilmenite incorporates about 20-30 times higher amounts of Nb and Ta than amphibole. Furthermore, trace element mineral analyses provide evidence for magmatic processes such as melt depletion, melt segregation, accumulation and fractionation as well as metasomatism having operated in this high-grade anatectic area.