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.

Thermodynamics of Forming Water Clusters at Various Temperatures and Pressures by Gaussian-2, Gaussian-3, Complete Basis Set-QB3, and Complete Basis Set-APNO Model Chemistries; Implications for Atmospheric Chemistry

The Gaussian-2, Gaussian-3, complete basis set- (CBS-) QB3, and CBS-APNO methods have been used to calculate ΔH° and ΔG° values for neutral clusters of water, (H2O)n, where n = 2−6. The structures are similar to those determined from experiment and from previous high-level calculations. The thermodynamic calculations by the G2, G3, and CBS-APNO methods compare well against the estimated MP2(CBS) limit. The cyclic pentamer and hexamer structures release the most heat per hydrogen bond formed of any of the clusters. While the cage and prism forms of the hexamer are the lowest energy structures at very low temperatures, as temperature is increased the cyclic structure is favored. The free energies of cluster formation at different temperatures reveal interesting insights, the most striking being that the cyclic trimer, cyclic tetramer, and cyclic pentamer, like the dimer, should be detectable in the lower troposphere. We predict water dimer concentrations of 9 × 1014 molecules/cm3, water trimer concentrations of 2.6 × 1012 molecules/cm3, tetramer concentrations of approximately 5.8 × 1011 molecules/cm3, and pentamer concentrations of approximately 3.5 × 1010 molecules/cm3 in saturated air at 298 K. These results have important implications for understanding the gas-phase chemistry of the lower troposphere.

Determining gypsum growth temperatures using monophase fluid inclusions-Application to the giant gypsum crystals of Naica, Mexico

Determining the formation temperature of minerals using fluid inclusions is a crucial step in understanding rock-forming scenarios. Unfortunately, fluid inclusions in minerals formed at low temperature, such as gypsum, are commonly in a metastable monophase liquid state. To overcome this problem, ultra-short laser pulses can be used to induce vapor bubble nucleation, thus creating a stable two-phase fluid inclusion appropriate for subsequent measurements of the liquid-vapor homogenization temperature, T-h. In this study we evaluate the applicability of T-h data to accurately determine gypsum formation temperatures. We used fluid inclusions in synthetic gypsum crystals grown in the laboratory at different temperatures between 40 degrees C and 80 degrees C under atmospheric pressure conditions. We found an asymmetric distribution of the T-h values, which are systematically lower than the actual crystal growth temperatures, T-g; this is due to (1) the effect of surface tension on liquid-vapor homogenization, and (2) plastic deformation of the inclusion walls due to internal tensile stress occurring in the metastable state of the inclusions. Based on this understanding, we have determined growth temperatures of natural giant gypsum crystals from Naica (Mexico), yielding 47 +/- 1.5 degrees C for crystals grown in the Cave of Swords (120 m below surface) and 54.5 +/- 2 degrees C for giant crystals grown in the Cave of Crystals (290 m below surface). These results support the earlier hypothesis that the population and the size of the Naica crystals were controlled by temperature. In addition, this experimental method opens a door to determining the growth temperature of minerals forming in low-temperature environments.

Lake surface water temperatures of European Alpine lakes (1989–2013) based on the Advanced Very High Resolution Radiometer (AVHRR) 1 km data set

Lake water temperature (LWT) is an important driver of lake ecosystems and it has been identified as an indicator of climate change. Consequently, the Global Climate Observing System (GCOS) lists LWT as an essential climate variable. Although for some European lakes long in situ time series of LWT do exist, many lakes are not observed or only on a non-regular basis making these observations insufficient for climate monitoring. Satellite data can provide the information needed. However, only few satellite sensors offer the possibility to analyse time series which cover 25 years or more. The Advanced Very High Resolution Radiometer (AVHRR) is among these and has been flown as a heritage instrument for almost 35 years. It will be carried on for at least ten more years, offering a unique opportunity for satellite-based climate studies. Herein we present a satellite-based lake surface water temperature (LSWT) data set for European water bodies in or near the Alps based on the extensive AVHRR 1 km data record (1989–2013) of the Remote Sensing Research Group at the University of Bern. It has been compiled out of AVHRR/2 (NOAA-07, -09, -11, -14) and AVHRR/3 (NOAA-16, -17, -18, -19 and MetOp-A) data. The high accuracy needed for climate related studies requires careful pre-processing and consideration of the atmospheric state. The LSWT retrieval is based on a simulation-based scheme making use of the Radiative Transfer for TOVS (RTTOV) Version 10 together with ERA-interim reanalysis data from the European Centre for Medium-range Weather Forecasts. The resulting LSWTs were extensively compared with in situ measurements from lakes with various sizes between 14 and 580 km2 and the resulting biases and RMSEs were found to be within the range of −0.5 to 0.6 K and 1.0 to 1.6 K, respectively. The upper limits of the reported errors could be rather attributed to uncertainties in the data comparison between in situ and satellite observations than inaccuracies of the satellite retrieval. An inter-comparison with the standard Moderate-resolution Imaging Spectroradiometer (MODIS) Land Surface Temperature product exhibits RMSEs and biases in the range of 0.6 to 0.9 and −0.5 to 0.2 K, respectively. The cross-platform consistency of the retrieval was found to be within ~ 0.3 K. For one lake, the satellite-derived trend was compared with the trend of in situ measurements and both were found to be similar. Thus, orbital drift is not causing artificial temperature trends in the data set. A comparison with LSWT derived through global sea surface temperature (SST) algorithms shows lower RMSEs and biases for the simulation-based approach. A running project will apply the developed method to retrieve LSWT for all of Europe to derive the climate signal of the last 30 years. The data are available at doi:10.1594/PANGAEA.831007.

Cool Glacial Temperatures and Changes in Moisture Source Recorded in Oman Groundwaters

Concentrations of atmospheric noble gases (neon, argon, krypton, and xenon) dissolved in groundwaters from northern Oman indicate that the average ground temperature during the Late Pleistocene (15,000 to 24,000 years before present) was 6.5° ± 0.6°C lower than that of today. Stable oxygen and hydrogen isotopic groundwater data show that the origin of atmospheric water vapor changed from a primarily southern, Indian Ocean source during the Late Pleistocene to a dominantly northern, Mediterranean source today. The reduced northern water vapor source is consistent with a drier Last Glacial Maximum through much of northern Africa and Arabia.

(Table 1) Oxygen isotope temperatures from hydrothermally altered breccias at DSDP Hole 60-453

Gabbro-metabasalt polymict breccias cored in Deep Sea Drilling Project Hole 453 are cemented in part by hydrothermal alteration to lower greenschist facies (chlorite-epidote-actinolite) mineral assemblages. Temperature estimates for this alteration, based on oxygen isotope determinations of secondary minerals, are nearly 100°C at the top of the breccias and over 200°C in a zone of intense alteration near the base.

Cenomanian/Turonian sea surface temperatures of the equatorial Atlantic reconstructed from geochemical data

Oceanic anoxic event 2 (OAE-2) occurring during the Cenomanian/Turonian (C/T) transition is evident from a globally recognized positive stable carbon isotopic excursion and is thought to represent one of the most extreme carbon cycle perturbations of the last 100 Myr. However, the impact of this major perturbation on and interaction with global climate remains unclear. Here we report new high-resolution records of sea surface temperature (SST) based on TEX86 and d 18O of excellently preserved planktic foraminifera and stable organic carbon isotopes across the C/T transition from black shales located offshore Suriname/French Guiana (Demerara Rise, Ocean Drilling Program Leg 207 Site 1260) and offshore Senegal (Cape Verde Basin, Deep Sea Drilling Project Leg 41 Site 367). At Site 1260, where both SST proxy records can be determined, a good match between conservative SST estimates from TEX86 and d 18O is observed. We find that late Cenomanian SSTs in the equatorial Atlantic Ocean (33°C) were substantially warmer than today (27°-29°C) and that the onset of OAE-2 coincided with a rapid shift to an even warmer (35°-36°C) regime. Within the early stages of the OAE a marked (4°C) cooling to temperatures lower than pre-OAE conditions is observed. However, well before the termination of OAE-2 the warm regime was reestablished and persisted into the Turonian. Our findings corroborate the view that the C/T transition represents the onset of the interval of peak Cretaceous warmth. More importantly, they are consistent with the hypotheses that mid-Cretaceous warmth can be attributed to high levels of atmospheric carbon dioxide (CO2) and that major OAEs were capable of triggering global cooling through the negative feedback effect of organic carbon-burial-led CO2 sequestration. Evidently, however, the factors that gave rise to the observed shift to a warmer climate regime at the onset of OAE-2 were sufficiently powerful that they were only briefly counterbalanced by the high rates of carbon burial attained during even the most extreme interval of organic carbon burial in the last 100 Myr.

Sea surface temperatures during Marine Isotope Stage (MIS) 11

The Marine Isotope Stage (MIS) 11 (424-374 ka) was characterized by a protracted deglaciation and an unusually long climatic optimum. It remains unclear to what degree the climate development during this interglacial reflects the unusually weak orbital forcing or greenhouse gas trends. Previously, arguments about the duration and timing of the MIS11 climatic optimum and about the pace of the deglacial warming were based on a small number of key records, which appear to show regional differences. In order to obtain a global signal of climate evolution during MIS11, we compiled a database of 78 sea surface temperature (SST) records from 57 sites spanning MIS11, aligned these individually on the basis of benthic (N = 28) or planktonic (N = 31) stable oxygen isotope curves to a common time frame and subjected 48 of them to an empirical orthogonal function (EOF) analysis. The analysis revealed a high commonality among all records, with the principal SST trend explaining almost 49% of the variability. This trend indicates that on the global scale, the surface ocean underwent rapid deglacial warming during Termination V, in pace with carbon dioxide rise, followed by a broad SST optimum centered at ~410 kyr. The second EOF, which explained ~18% of the variability, revealed the existence of a different SST trend, characterized by a delayed onset of the temperature optimum during MIS11 at ~398 kyr, followed by a prolonged warm period lasting beyond 380 kyr. This trend is most consistently manifested in the mid-latitude North Atlantic and Mediterranean Sea and is here attributed to the strength of the Atlantic meridional overturning circulation. A sensitivity analysis indicates that these results are robust to record selection and to age-model uncertainties of up to 3-6 kyr, but more sensitive to SST seasonal attribution and SST uncertainties >1 °C. In order to validate the CCSM3 (Community Climate System Model, version 3) predictive potential, the annual and seasonal SST anomalies recorded in a total of 74 proxy records were compared with runs for three time slices representing orbital configuration extremes during the peak interglacial of MIS11. The modeled SST anomalies are characterized by a significantly lower variance compared to the reconstructions. Nevertheless, significant correlations between proxy and model data are found in comparisons on the seasonal basis, indicating that the model captures part of the long-term variability induced by astronomical forcing, which appears to have left a detectable signature in SST trends.

Mean air temperatures of Greenland

Melt rate and surface temperature on the Greenland ice sheet are parameterized in terms of snow accumulation, mean annual air temperatur and mean July air temperature. Melt rates are calculated using positive degree-days, and firn warming (i.e. the positive deviation of the temperature at 10-15 m depth from the mean annual air temperature) is estimated from the calculated amount of refrozen melt water in the firn. A comparison between observed and calculated melt rates shows that the parameterization provides a reasonable estimate of the present ablation rates in West Greenland between 61°N and 76°N. The average equilibrium line elevation is estimated to be about 1150 m and 1000 m for West and East Greenland respectively, which is several hundred meter lower than previous estimates. However, the total annual ablation from the ice sheet is found to be about 280 km**3 of water per year which agrees well with most other estimates. The melt-rate model predicts significant melting and consequently significant firn warming even at the highest elevations of the South Greenland ice sheet, whereas a large region of central Greenland north of 70° N experiences little or no summer melting. This agrees with the distribution of the dry snow facics as given by BENSON (1962).

Tab. 1: Januar temperatures on location of the high northern hemisphere (Lustrenwerte)

Während sich Mitteleuropa bereits seit 1920 in einer Phase sinkender Januartemperatur befindet, trat dieser Klima-Umschwung in höheren Breiten immer später ein, so waren in Königsberg die beiden Januar Lustren 1916-1920 und 1921-1925 gleich warm, und in Upsala und Haparanda waren die Jahrfünfte mildester Januartemperatur in beiden Fällen 1931-35 - eine Erscheinung von übrigens weltweitem Charakter, denn auch weiteste Teile der U.S.A. hatten in dieser Zeit die mildesten Januarmonate.

Alkenone-derived UK'37 data and sea surface temperatures (SST) of sediment cores from the Gulf of Taranto (Italy)

Auxiliary data include one file with alkenone-derived UK'37 data and sea surface temperatures (SST). On these data Figs. 7 and 8 of the manuscript are based. The SST are derived from UK'37 by using the transfer function: SST = 29.876 UK'37 - 1.334 of Conte et al. (2006). The data are against the ages (in A.D.) of samples derived from cores GT91-1 (39[deg]59'23"N, 17[deg]45'25"E), GT89-3 and GT90-3 (both 39[deg]45'43"N, 17[deg]53'55"E ). Also included are composite records for UK'37 and SST. For creating the composite records, GT-89-3 was taken as reference core. In the overlapping period the GT89-3 data seem in general lower than the GT91-1 data. To accommodate for this in the composite record, the average difference (0.0343 UK'37 units; equivalent to 1.023 [deg]C) was subtracted from the GT91-1 record. Hereafter, for each depth in the overlapping interval the respective values (UK'37 or SST) of GT89-3 and GT91-1 were averaged. We have also averaged with 16 additional alkenone measurements, from 1793 to 1851, performed in the GT90-3 core.

(Table 1) Viscosity coefficients at various temperatures of samples from ODP Hole 176-735B

Ocean Drilling Program (ODP) Hole 735B was drilled to a depth of 1.5 km in a tectonic window of gabbroic lower oceanic crust created at the Southwest Indian Ridge. The gabbros have a very stable natural remanent magnetization (NRM) of reversed polarity with most unblocking temperatures slightly below the Curie temperature of magnetite. The NRM includes a drilling-induced overprint but its intensity decays strongly towards the interior of the drill core. The demagnetization data yield no or only a very small secondary magnetization component acquired during the present Brunhes chron or an earlier normal chron, suggesting cooling through most of the blocking temperature range during chron C5r and a strong resistance against the acquisition of thermoviscous magnetization. A novel furnace has been designed to measure magnetizations and their time dependences at high temperatures (up to 580 deg C) inside a commercial SQUID magnetometer. Magnetic viscosity experiments have been conducted on the gabbros at temperatures up to 550 deg C to determine the time and temperature stability of remanent magnetization. Viscosities are generally small and increase little with temperature below the main blocking temperature, where the increase becomes almost an order of magnitude. Extrapolations to geological times infer viscous acquisitions that would be 5-25% of a thermoremanence in 100 kyr and at temperatures of 200-500 deg C. At ocean bottom temperature the predicted magnetization of one sample acquired in the present Brunhes chron should be 10% of the NRM. However, this is not recognized during NRM demagnetization and partial thermoremanent magnetization (pTRM) acquisitions at 250 deg C are also much smaller than predicted. It thus appears that the NRMs are generally magnetically harder than magnetizations acquired after heating to 570 deg C in the laboratory. Susceptibility changes during heating are small (<5%) indicating a seemingly stable magneto-mineralogy, but conspicuous minima occur after heating to 520 deg C. Also, quasi paleointensity experiments reveal characteristic patterns in the NRM/pTRM ratios and also large increases in pTRM capacity after heating to 570 deg C. Moreover, anhysteretic remanent magnetization acquisition in the low field range (<=10 mT) is strongly enhanced after heating by factors up to three. The alteration of the magneto-mineralogy is interpreted to result from the annealing of defects in magnetite that originate from tectonically induced strain. The oceanic gabbros of Hole 735B are thus ideal source layer material for marine magnetic anomalies, and secondary thermoviscous acquisition, as a possible cause for anomalous skewness, is essentially absent.

Calculated sea level records and deep ocean temperatures from DSDP and ODP sites

During the mid-Pleistocene transition the dominant 41 ka periodicity of glacial cycles transitioned to a quasi-100 ka periodicity for reasons not yet known. This study investigates the potential role of deep ocean hydrography by examining oxygen isotope ratios in benthic foraminifera. Oxygen isotope records from the Atlantic, Pacific and Indian Ocean basins are separated into their ice volume and local temperature/hydrography components using a piece-wise linear transfer function and a temperature calibration. Although our method has certain limitations, the deep ocean hydrography reconstructions show that glacial deep ocean temperatures approached freezing point as the mid-Pleistocene transition progressed. Further analysis suggests that water mass reorganisation could have been responsible for these temperature changes, leading to such stable conditions in the deep ocean that some obliquity cycles were skipped until precessional forcing triggered deglaciation, creating the apparent quasi-100 ka pattern. This study supports previous work that suggests multiples of obliquity cycles dominate the quasi-100 ka glacial cycles with precession components driving deglaciations.