Oceanic crust and island arc formation in Central Asia during late Neoproterozoic times - evidence from petrological and geochemical studies

Die vorliegende Arbeit behandelt die Entwicklung des 570 Ma alten, neoproterozoischen Agardagh - Tes-Chem Ophioliths (ATCO) in Zentralasien. Dieser Ophiolith liegt südwestlich des Baikalsees (50.5° N, 95° E) und wurde im frühen Stadium der Akkretion des Zentralasiatischen Mobilgürtels auf den nordwestlichen Rand des Tuvinisch-Mongolischen Mikrokontinentes aufgeschoben. Bei dem Zentralasiatische Mobilgürtel handelt es sich um einen riesigen Akkretions-Subduktionskomplex, der heute das größte zusammenhängende Orogen der Erde darstellt. Im Rahmen dieser Arbeit wurden eine Reihe plutonischer und vulkanischer Gesteine, sowie verschiedene Mantelgesteine des ATCO mittels mikroanalytischer und geochemischer Verfahren untersucht (Elektronenstrahlmikrosonde, Ionenstrahlmikrosonde, Spurenelement- und Isotopengeochemie). Die Auswertung dieser Daten ermöglichte die Entwicklung eines geodynamisch-petrologischen Modells zur Entstehung des ATCO. Die vulkanischen Gesteine lassen sich aufgrund ihrer Spurenelement- und Isotopenzusammensetzung in inselbogenbezogene und back-arc Becken bezogene Gesteine (IA-Gesteine und BAB-Gesteine) unterscheiden. Darüber hinaus gibt es eine weitere, nicht eindeutig zuzuordnende Gruppe, die hauptsächlich mafische Gänge umfasst. Der grösste Teil der untersuchen Vulkanite gehört zur Gruppe der IA-Gesteine. Es handelt sich um Al-reiche Basalte und basaltische Andesite, welche aus einem evolvierten Stammmagma mit Mg# 0.60, Cr ~ 180 µg/g und Ni ~ 95 µg/g hauptsächlich durch Klinopyroxenfraktionierung entstanden sind. Das Stammmagma selbst entstand durch Fraktionierung von ca. 12 % Olivin und geringen Anteilen von Cr-Spinell aus einer primären, aus dem Mantel abgeleiteten Schmelze. Die IA-Gesteine haben hohe Konzentrationen an inkompatiblen Spurenelementen (leichte-(L)- Seltenerdelement-(SEE)-Konzentrationen etwa 100-fach chondritisch, chondrit-normierte (La/Yb)c von 14.6 - 5.1), negative Nb-Anomalien (Nb/La = 0.37 - 0.62) und niedrige Zr/Nb Verhältnisse (7 - 14) relativ zu den BAB-Gesteinen. Initiale eNd Werte liegen bei etwa +5.5, initiale Bleiisotopenverhältnisse sind: 206Pb/204Pb = 17.39 - 18.45, 207Pb/204Pb = 15.49 - 15.61, 208Pb/204Pb = 37.06 - 38.05. Die Anreicherung lithophiler inkompatibler Spurenelemente (LILE) in dieser Gruppe ist signifikant (Ba/La = 11 - 130) und zeigt den Einfluss subduzierter Komponenten an. Die BAB-Gesteine repräsentieren Schmelzen, die sehr wahrscheinlich aus der gleichen Mantelquelle wie die IA-Gesteine stammen, aber durch höhere Aufschmelzgrade (8 - 15 %) und ohne den Einfluss subduzierter Komponenten entstanden sind. Sie haben niedrigere Konzentrationen an inkompatiblen Spurenelementen, flache SEE-Muster ((La/Yb)c = 0.6 - 2.4) und höhere initiale eNd Werte zwischen +7.8 und +8.5. Nb Anomalien existieren nicht und Zr/Nb Verhältnisse sind hoch (21 - 48). Um die geochemische Entwicklung der vulkanischen Gesteine des ATCO zu erklären, sind mindestens drei Komponenten erforderlich: (1) eine angereicherte, ozeaninselbasalt-ähnliche Komponente mit hoher Nb Konzentration über ~ 30 µg/g, einem niedrigen Zr/Nb Verhältnis (ca. 6.5), einem niedrigen initialen eNd Wert (um 0), aber mit radiogenen 206Pb/204Pb-, 207Pb/204Pb- und 208Pb/204Pb-Verhältnissen; (2) eine N-MORB ähnliche back-arc Becken Komponente mit flachem SEE-Muster und einem hohen initialen eNd Wert von mindestens +8.5, und (3) eine Inselbogen-Komponente aus einer verarmten Mantelquelle, welche durch die abtauchende Platte geochemisch modifiziert wurde. Die geochemische Entstehung der ATCO Vulkanite lässt sich dann am besten durch eine Kombination aus Quellenkontamination, fraktionierte Kristallisation und Magmenmischung erklären. Geodynamisch gesehen entstand der ATCO sehr wahrscheinlich in einem intraozeanischen Inselbogen - back-arc System. Bei den untersuchten Plutoniten handelt es sich um ultramafische Kumulate (Wehrlite und Pyroxenite) sowie um gabbroische Plutonite (Olivin-Gabbros bis Diorite). Die geochemischen Charakteristika der mafischen Plutonite sind deutlich unterschiedlich zu denen der vulkanischen Gesteine, weshalb sie sehr wahrscheinlich ein späteres Entwicklungsstadium des ATCO repräsentieren. Die Spurenelement-Konzentrationen in den Klinopyroxenen der ultramafischen Kumulate sind extrem niedrig, mit etwa 0.1- bis 1-fach chondritischen SEE-Konzentrationen und mit deutlich LSEE-verarmten Mustern ((La/Yb)c = 0.27 - 0.52). Berechnete Gleichgewichtsschmelzen der ultramafischen Kumulate zeigen grosse Ähnlichkeit zu primären boninitischen Schmelzen. Die primären Magmen waren daher boninitischer Zusammensetzung und entstanden in dem durch vorausgegangene Schmelzprozesse stark verarmten Mantelkeil über einer Subduktionszone. Niedrige Spurenelement-Konzentrationen zeigen einen geringen Einfluss der abtauchenden Platte an. Die Spurenelement-Konzentrationen der Gabbros sind ebenfalls niedrig, mit etwa 0.5 - 10-fach chondritischen SEE-Konzentrationen und mit variablen SEE-Mustern ((La/Yb)c = 0.25 - 2.6). Analog zu den Vulkaniten der IA-Gruppe haben alle Gabbros eine negative Nb-Anomalie mit Nb/La = 0.01 - 0.31. Die initialen eNd Werte der Gabbros variieren zwischen +4.8 und +7.1, mit einem Mittelwert von +5.9, und sind damit identisch mit denen der IA-Vulkanite. Bei den untersuchten Mantelgesteinen handelt es sich um teilweise serpentinisierte Dunite und Harzburgite, die alle durch hohe Mg/Si- und niedrige Al/Si-Verhältnisse gekennzeichnet sind. Dies zeigt einen refraktären Charakter an und steht in guter Übereinstimmung mit den hohen Cr-Zahlen (Cr#) der Spinelle (bis zu Cr# = 0.83), auf deren Basis der Aufschmelzgrad der residuellen Mantelgesteine berechnet wurde. Dieser beträgt etwa 25 %. Die geochemische Zusammensetzung und die petrologischen Daten der Ultramafite und Gabbros lassen sich am besten erklären, wenn man für die Entstehung dieser Gesteine einen zweistufigen Prozess annimmt. In einer ersten Stufe entstanden die ultramafischen Kumulate unter hohem Druck in einer Magmenkammer an der Krustenbasis, hauptsächlich durch Klinopyroxen-Fraktionierung. Bei dieser Magmenkammer handelte es sich um ein offenes System, dem von unten laufend neue Schmelze zugeführt wurde, und aus dem im oberen Bereich evolviertere Schmelzen geringerer Dichte entwichen. Diese evolvierten Schmelzen stiegen in flachere krustale Bereiche auf und bildeten dort meist isolierte Intrusionskörper. Diese Intrusionskörper erstarrten ohne Magmen-Nachschub, weshalb petrographisch sehr unterschiedliche Gesteine entstehen konnten. Eine geochemische Modifikation der abkühlenden Schmelzen erfolgte allerdings durch die Assimilation von Nebengestein. Da innerhalb der Gabbros keine signifikante Variation der initalen eNd Werte existiert, handelte es sich bei dem assimilierten Material hauptsächlich um vulkanische Gesteine des ATCO und nicht um ältere, möglicherweise kontinentale Kruste.

Isotope geochemical and geochronologic studies of basement units in the Central Western Carpathians (Slovakia)

In dieser Arbeit werden geochronologische und isotopen-geochemische Daten zur Entwicklung der Zentralen Westlichen Karpathen präsentiert. Die Karpathen bilden die östliche Fortsetzung der Alpen und können in drei Alpine Grundgebirgsdecken unterteilt werden, von denen zwei, die Veporische und die Gemerische, bearbeitet wurden. In der Veporischen Einheit wurden polymetamorphe Grundgebirgseinheiten untersucht, um deren genaue Altersstellung zu definieren und sie isotopengeochemisch zu klassifizieren. Dagegen wurde in der der Gemerischen Einheit, welche die Veporische Einheit überlagert, ein spezialisierter S-Typ Granit im Detail untersucht, um die petrogenetischen Prozesse, die zur magmatischen Entwicklung dieses Granits geführt haben, zu identifizieren. U-Pb Datierungen an Zirkonen der Veporischen Grundgebirgseinheiten zeigen für die gesamte Veporische Einheit ordovizische Entsehungsalter an (440-470 Ma). Diese Datierungen revidieren publizierte kambrische Entstehungsalter dieses Grundgebirges. Die Isotopensignatur (epsilon Nd und 87Sr/86Sr) der ordovizischen Grundgebirgseinheiten, bestehend aus stark überprägten Amphiboliten und Gneissen, ist von der Signatur der sich im Norden anschliessenden Tatrischen Einheit gut unterscheidbar. Die Bleiisotopenzusammensetzung dieser Gesteine ist stark krustal geprägt und überschneidet sich mit der der Tatrischen Einheit. Zusammen mit den T-DM Altern sind diese Einheiten vergleichbar mit prävariskischen Einheiten der Alpen. Somit kann das ordovizische Grundgebirge zu den peri-Gondwana Terranen gezählt werden, die an einem aktiven Kontinentalrand im Norden von Gondwana gebildet wurden. In den Gesteinen der Veporischen Einheit wurde im Weiteren eine starke metamorphe überprägung und intensiver felsischer Magmatismus karbonischen Alters erkannt (320-350 Ma). Dieses Ereignis ist zeitgleich mit dem Magmatismus, welcher hauptsächlich in der sich im Norden anschliessenden Tatrischen Einheit beobachtet wird. Dieser gehört der variskischen Orogenese an. Intensive alpine Deformation und Metamorphose konnte in der südlichen Veporischen Einheit anhand der Einzelzirkondatierungen und der Isotopendaten der ordovizischen Einheiten nachgewiesen werden. Am Dlha Dolina Granit in der Gemerischen Einheit können starke Fraktionierungs- und Auto-Metasomatose-Effekte beobachtet werden. Durch die magmatische Fraktionierung wird eine Anreicherung der SEE erzeugt, wogegen die Metasomatose die SEE stark verarmt. Es kommt sogar zur Ausbildung eines Tetraden Effektes im SEE Muster, welche den starken Einfluss von Fluiden während der spät-magmatischen Phase belegt. Gesamtgesteins Pb-Pb Daten beschränken das minimale Intrusionsalter dieses Granites auf 240 Ma. Dieses Alter ist in guter übereinstimmung mit den Sr-Isotopendaten der magmatisch dominierten Gesteine, wohingegen die stark metasomatisch geprägten Gesteine ein zu radiogenes 87Sr/86Sri aufweisen. Während dieser Arbeit wurde intensiv mit der Blei-Isotopenzusammensetzung von Gesamtgesteinsproben gearbeitet. Um die Auswertung dieser Daten optimieren zu können wurde ein Computerscript für das GPL Programm Octave erstellt. Die Hauptaufgabe dieses Scripts besteht darin, Regressionen für geochronologische Anwendungen gemäss York (1969) zu berechnen. Ausserdem können mu und kappa-Werte für diese Regressionen berechnet und eine Hauptkomponentenanalyse, welche hilfreich für den Vergleich von zwei Datensätzen ist, durchgeführt werden. Am Ende der vorliegenden Arbeit wird die analytische Methode für einen Mikrowellen beschleunigten Säureaufschluss von granitoidem Material zur Bestimmung der Sr- und Nd-Isotopenzusammensetzung und der Elementkonzentrationen vorgestellt. Diese kombinierte Methode nutzt ein TIMS für die Sr und Nd Isotopenmessungen und eine Einzelkollektor-ICPMS zur Bestimmung der SEE, Rb und Sr Konzentrationen, welche mithilfe von relativen Sensitivitätsfaktoren gegenüber einem internen Standard quantifiziert werden. Diese Methode wird durch Messungen von internationalen Referenzmaterialien bewertet. Die Ergebnisse zeigen eine Reproduzierbarkeit von <10% für die Elementkonzentrationen und von <5% für Elementverhältnisse.

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.

Continuous geochemical monitoring by mass-spectometer in the Campi Flegrei geothermal area. An application at Pisciarelli-Solfatara (diffuse and fumarolic gases) and at the mud gases during drilling of the CFDDP pilot hole

During this work has been developed an innovative methodology for continuous and in situ gas monitoring (24/24 h) of fumarolic and soil diffusive emissions applied to the geothermal and volcanic area of Pisciarelli near Agnano inside the Campi Flegrei caldera (CFc). In literature there are only scattered and in discrete data of the geochemical gas composition of fumarole at Campi Flegrei; it is only since the early ’80 that exist a systematic record of fumaroles with discrete sampling at Solfatara (Bocca Grande and Bocca Nuova fumaroles) and since 1999, even at the degassing areas of Pisciarelli. This type of sampling has resulted in a time series of geochemical analysis with discontinuous periods of time set (in average 2-3 measurements per month) completely inadequate for the purposes of Civil Defence in such high volcanic risk and densely populated areas. For this purpose, and to remedy this lack of data, during this study was introduced a new methodology of continuous and in situ sampling able to continuously detect data related and from its soil diffusive degassing. Due to its high sampling density (about one measurement per minute therefore producing 1440 data daily) and numerous species detected (CO2, Ar, 36Ar, CH4, He, H2S, N2, O2) allowing a good statistic record and the reconstruction of the gas composition evolution of the investigated area. This methodology is based on continuous sampling of fumaroles gases and soil degassing using an extraction line, which after undergoing a series of condensation processes of the water vapour content - better described hereinafter - is analyzed through using a quadrupole mass spectrometer

Neoproterozoic to middle Palaeozoic evolution of the Central Asian Orogenic Belt in south-central Mongolia: chronological and geochemical perspectives

Mongolia occupies a central position within the eastern branch of the large accretionary Central Asian Orogenic Belt (CAOB) or Altaids. The present work aims to outline the geodynamic environment and possible evolution of this part of the eastern CAOB, predominantly from the Cambrian to the middle Palaeozoic. The investigation primarily focussed on zircon geochronology as well as whole-rock geochemical and Sm–Nd isotopic analyses for a variety of metaigneous rocks from the southern Hangay and Gobi-Altai regions in south-central Mongolia. The southern slope of the Hangay Mountains in central Mongolia exposes a large NWSE-trending middle Neoproterozoic ophiolitic complex (c. 650 Ma), which is tectonically integrated within an accretionary complex developed between the Precambrian Baydrag and Hangay crustal blocks. Formation of the entire accretionary system along the north-eastern margin of the Baydrag block mainly occurred during the early Cambrian, but convergence within this orogenic zone continued until the early Ordovician, because of on-going southward subduction-accretion of the Baydrag block. An important discovery is the identification of a late Mesoproterozoic to early Neoproterozoic belt within the northern Gobi-Altai that was reworked during the late Cambrian and throughout the late Ordovician/Devonian. Early Silurian low-grade mafic and felsic metavolcanic rocks from the northern Gobi-Altai display subduction-related geochemical features and highly heterogeneous Nd isotopic compositions, which suggest an origin at a mature active continental margin. Early Devonian protoliths of granodioritic and mafic gneisses from the southern Gobi-Altai display geochemical and Nd isotopic compositions compatible with derivation and evolution from predominantly juvenile crustal and mantel sources and these rocks may have been emplaced within the outboard portion of the late Ordovician/early Silurian active continental margin. Moreover, middle Devonian low-grade metavolcanic rocks from the southwestern Gobi-Altai yielded geochemical and Nd isotopic data consistent with emplacement in a transitional arc-backarc setting. The combined U–Pb zircon ages and geochemical data obtained from the Gobi-Altai region suggest that magmatism across an active continental margin migrated oceanwards through time by way of subduction zone retreat throughout the Devonian. Progressive extension of the continental margin was associated with the opening of a backarc basin and culminated in the late Devonian with the formation of a Japan-type arc front facing a southward open oceanic realm (present-day coordinates).

Stream sediments analysis for geochemical mapping of Romagna Apennines (Northern Italy): monitoring and management tool of environmental resources at various scales

Geochemical mapping is a valuable tool for the control of territory that can be used not only in the identification of mineral resources and geological, agricultural and forestry studies but also in the monitoring of natural resources by giving solutions to environmental and economic problems. Stream sediments are widely used in the sampling campaigns carried out by the world's governments and research groups for their characteristics of broad representativeness of rocks and soils, for ease of sampling and for the possibility to conduct very detailed sampling In this context, the environmental role of stream sediments provides a good basis for the implementation of environmental management measures, in fact the composition of river sediments is an important factor in understanding the complex dynamics that develop within catchment basins therefore they represent a critical environmental compartment: they can persistently incorporate pollutants after a process of contamination and release into the biosphere if the environmental conditions change. It is essential to determine whether the concentrations of certain elements, in particular heavy metals, can be the result of natural erosion of rocks containing high concentrations of specific elements or are generated as residues of human activities related to a certain study area. This PhD thesis aims to extract from an extensive database on stream sediments of the Romagna rivers the widest spectrum of informations. The study involved low and high order stream in the mountain and hilly area, but also the sediments of the floodplain area, where intensive agriculture is active. The geochemical signals recorded by the stream sediments will be interpreted in order to reconstruct the natural variability related to bedrock and soil contribution, the effects of the river dynamics, the anomalous sites, and with the calculation of background values be able to evaluate their level of degradation and predict the environmental risk.

Geochemical and Strontium Isotope Characterization of Produced Waters from Marcellus Shale Natural Gas Extraction

Extraction of natural gas by hydraulic fracturing of the Middle Devonian Marcellus Shale, a major gas-bearing unit in the Appalachian Basin, results in significant quantities of produced water containing high total dissolved solids (TDS). We carried out a strontium (Sr) isotope investigation to determine the utility of Sr isotopes in identifying and quantifying the interaction of Marcellus Formation produced waters with other waters in the Appalachian Basin in the event of an accidental release, and to provide information about the source of the dissolved solids. Strontium isotopic ratios of Marcellus produced waters collected over a geographic range of ∼375 km from southwestern to northeastern Pennsylvania define a relatively narrow set of values (εSr SW = +13.8 to +41.6, where εSr SW is the deviation of the 87Sr/86Sr ratio from that of seawater in parts per 104); this isotopic range falls above that of Middle Devonian seawater, and is distinct from most western Pennsylvania acid mine drainage and Upper Devonian Venango Group oil and gas brines. The uniformity of the isotope ratios suggests a basin-wide source of dissolved solids with a component that is more radiogenic than seawater. Mixing models indicate that Sr isotope ratios can be used to sensitively differentiate between Marcellus Formation produced water and other potential sources of TDS into ground or surface waters.

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).