4.4 billion years of crustal maturation: oxygen isotope ratios of magmatic zircon

Analysis of δ¹⁸O in igneous zircons of known age traces the evolution of intracrustal recycling and crust-mantle interaction through time. This record is especially sensitive because oxygen isotope ratios of igneous rocks are strongly affected by incorporation of supracrustal materials into melts, which commonly have δ¹⁸O values higher than in primitive mantle magmas. This study summarizes data for δ¹⁸O in zircons that have been analyzed from 1,200 dated rocks ranging over 96% of the age of Earth. Uniformly primitive to mildly evolved magmatic δ¹⁸O values are found from the first half of Earth history, but much more varied values are seen for younger magmas. The similarity of values throughout the Archean, and comparison to the composition of the “modern” mantle indicate that δ¹⁸O of primitive mantle melts have remained constant (±0.2‰) for the past 4.4 billion years. The range and variability of δ¹⁸O in all Archean zircon samples is subdued (δ¹⁸O(Zrc)=5–7.5‰) ranging from values in high temperature equilibrium with the mantle (5.3± 0.3‰) to slightly higher, more evolved compositions (6.5–7.5‰) including samples from: the Jack Hills (4.4–3.3 Ga), the Beartooth Mountains (4.0–2.9 Ga), Barberton (3.5–2.7 Ga), the Superior and Slave Provinces (3.0 to 2.7 Ga), and the Lewisian (2.7 Ga). No zircons from the Archean have been analyzed with magmatic δ¹⁸O above 7.5‰. The mildly evolved, higher Archean values (6.5–7.5‰) are interpreted to result from exchange of protoliths with surface waters at low temperature followed by melting or contamination to create mildly elevated magmas that host the zircons. During the Proterozoic, the range of δ¹⁸O(Zrc) and the highest values gradually increased in a secular change that documents maturation of the crust. After ∼1.5 Ga, high δ¹⁸O zircons (8 to >10‰) became common in many Proterozoic and Phanerozoic terranes reflecting δ¹⁸O(whole rock) values from 9 to over 12‰. The appearance of high δ¹⁸O magmas on Earth reflects nonuniformitarian changes in the composition of sediments, and rate and style of recycling of surface-derived material into magmas within the crust. Electronic Supplementary Material Supplementary material is available for this article at     

The Greer Lake leucogranite, Manitoba, and the origin of lepidolite-subtype granitic pegmatites

The Archean Greer Lake leucogranite intruded metabasalts of the Bird River Greenstone Belt in the southwestern part of the Superior Province of southeastern Manitoba. The considerably evolved, multiphase, peraluminous, B-, P-, and S-poor leucogranite (K/Rb 132 to 24) was probably generated by fault-friction-assisted anatexis of dominantly metatonalitic rocks and subsequent differentiation. The leucogranite produced interior, transitional, non-crosscutting pods of barren, beryl-columbite- and lepidolite-subtype pegmatites that solidified from local segregations of highly fractionated residual melt. Steep fractionation gradients characterize the granite-to-pegmatite transition, most conspicuously so in the case of the most evolved, Li, Rb, Cs, Be, Mn, Sn, Nb-Ta, F-rich, lepidolite-subtype pod AC #3 (with K/Rb ≥ 16 and Cs 330 ppmwt in accessory K-feldspar, ≥2.5 and ≤11,200 ppmwt, respectively, in lepidolite, Cs ≤28,000 ppmwt in beryl, and Ta/(Ta+Nb) at. ≤ 0.95 in manganotantalite). The Greer Lake example documents beyond any doubt the igneous derivation of lepidolite-subtype pegmatites from a plutonic parent. Most cases of generally very scarce lepidolite-subtype pegmatites obscure this relationship, as the volatile-rich, highly fluid melts stable to relatively low temperatures commonly migrate to great distances from their plutonic sources.     

Hot-fluid Driven Metasomatism of Samalpatti Carbonatites, South India: Evidence from Petrology, Mineral Chemistry, Trace Elements and Stable Isotope Compositions

The magmatic heritage of carbonatites can be identified on the basis of a combination of geological criteria such as, their mode of occurrence, the nature of associated igneous rocks, the presence of minerals of igneous origin, fenitization, characteristic trace element contents and isotopic composition. Late Proterozoic Samalpatti carbonatites were studied in view of these criteria, and were found to contain metamorphic minerals that normally form under thermal metamorphic conditions and which have unusual chemical compositions. A combination of criteria points clearly to a magmatic origin for these carbonatites. Field relations indicate that the dominant modes of intrusion of carbonatite into the encompassing pyroxenites and syenites include small dykes, veins, or lenses. The igneous nature of these carbonatites has been described elsewhere and chemically they are classified as calico-carbonatites. Currently, very little is known about the metamorphic textures and mineralogy observed in the Samalpatti carbonatites. In this study, several metamorphic minerals are reported including diopside, grossularite, vesuvianite, K-feldspar and wollastonite, and a hornfelsic texture is described. These mineral phases and texture characterize thermal metamorphism under low pressure and high temperature (LP-HT) metamorphic conditions (650°_750°C) or metasomatism aided by hot-fluid advection. The metamorphic nature of minerals reported is also confirmed by electron microprobe study. The Samalpatti carbonatite samples show much lower values of characteristic trace elements (P, Sr, Ba, Zr, Nb, Th, Y and REEs) than average concentrations for magmatic carbonatite. Stable isotopic (d¹³C and d¹⁸O) compositions of Samalpatti carbonatites do not fall in the primary igneous carbonatite (PIC) domain. The petrological and chemical signatures of these carbonatites suggest metasomatism in conjunction with fluid advection. Such a metasomatic process may drastically change the chemistry of the rocks in addition to enrichment of heavier stable isotopes. During this metasomatic process, characteristic elements would be dissolved in the high d¹⁸O fluid, and together with Rayleigh fractionation would contribute to enhanced concentrations of ¹³C and ¹⁸O in Samalpatti carbonatites.     

Tectonic reconstruction of the northern Andean blocks: Oblique convergence and rotations derived from the kinematics of the Piedras–Girardot area, Colombia

A detailed kinematic study in the Piedras–Girardot area reveals that approximately 32 km of ENE–WSW oblique convergence is accommodated within a northeast-trending transpressional shear zone with a shear strain of 0.8 and a convergence factor of 2. Early Campanian deformation is marked by the incipient propagation of northeast-trending faults that uplifted gentle domes where the accumulation of sandy units did not take place. Maastrichtian unroofing of a metamorphic terrane to the west is documented by a conglomerate that was deformed shortly after deposition developing a conspicuous intragranular fabric of microscopic veins that accommodates less than 5% extension. This extensional fabric, distortion of fossil molds, and a moderate cleavage accommodating less than 5% contraction, developed concurrently, but before large-scale faulting and folding. Paleogene folding and southwestward thrust sheet propagation are recorded by syntectonic strata. Neogene deformation took place only in the western flank of this foldbelt. The amount, direction, and timing of deformation documented here contradict current tectonic models for the Cordillera Oriental and demand a new tectonic framework to approach the study of the structure of the northern Andes. Thus, an alternative model was constructed by defining three continental blocks: the Maracaibo, Cordillera Central, and Cordillera Oriental blocks. Oblique deformation imposed by the relative eastward and northeastward motion of the Caribbean Plate was modeled as rigid-body rotation and translation for rigid blocks (derived from published paleomagnetic and kinematic data), and as internal distortion and dilation for weak blocks (derived from the Piedras–Girardot area). This model explains not only coeval dextral and sinistral transpression and transtension, but also large clockwise rotation documented by paleomagnetic studies in the Caribbean–northern Andean region.     

Geomorphology and geochronology of sackung features (uphill-facing scarps) in the Central Spanish Pyrenees

The paper analyses uphill-facing scarps and associated troughs developed in the oversteepened slopes of two neighbouring glacial valleys in the central Spanish Pyrenees. Previous studies of sackung landforms in the Pyrenees have argued for deglacial unloading as the genetic mechanism, but this causal and temporal relationship has not been proved due to the lack of chronological data. The antislope scarps in the two studied locations, Vallibierna and Estós, are developed in Palaeozoic metasedimentary rocks (parallel to the contour lines and the structural grain), occur in the intermediate sector of the hillslope, and are up to 0.5 km long and several meters high. A trench was excavated in a sackung trough fill in each of the valleys in order to gain information about their chronology and genesis. Charcoal from the lowermost unit in Vallibierna provided an age of 5.9 cal. ka for the sackung and extrapolation of the three dates obtained in Estós indicates that the trough formed ca. 7.6–7.8 cal. ka. Deglaciation of the studied sectors of the valleys occurred between 16 and 13 ka. The time lag of >5 ka suggests that glacial erosion and the subsequent debutressing of the oversteepened valley walls created slopes predisposed to sackung development, but did not initiate the movement. Seismic shaking is proposed as a probable triggering factor. This hypothesis, although supported by the sudden deformation event recorded by a failure plane exposed in Vallibierna trench, and by the seismic and neotectonic activity of the area, cannot be proved due to the lack of chronological information about paleoearthquakes.     

European VOC Emission Estimates Evaluated by Measurements and Model Calculations

Observations and model calculations of the concentration of hydrocarbonsat five Scandinavian rural sites during March–June 1993are reported.Decreasing concentrations from March to June are observedat all sites. The highest concentrations of hydrocarbons were found in air massescoming in from the southwest to southeast, indicating that long rangetransport fromcontinental Europe and the U.K. is important in pollution episodes. An episode of elevated concentrations of hydrocarbons observed at three of the sites in the middle of Marchis described and discussed in relation to the synoptic situation and thepresenceof other chemical compounds (NO₂, PAN, total nitrate andozone).A Lagrangian numerical model is used to calculate the concentrations of theindividual hydrocarbons at the fivesites and comparison with observations is made.The calculated concentrations for nonmethane hydrocarbons with quite longchemicallifetimes agree well with the observations.For the sum of observed and calculated hydrocarbons the correlationcoefficientsare in the range of 0.65–0.88 for the five sitesand the ratio between calculated and measured concentrations was0.72–0.97, indicating thatthe European VOC emission inventory is quite well estimated.     

The sensitivity of winter evaporation to the formulation of aerodynamic resistance in the ECMWF model

In atmospheric models, the roughness length for momentum, heat and moisture are often taken equal, and tuned to the momentum budget problem. In this paper, it is shown that the roughness lengths have considerable impact on the evaporation in winter. One-column simulations of the land-surface scheme are driven with a long time series of observations for Cabauw in The Netherlands. It is shown that with the operational roughness lengths for this location (as in use at ECMWF in May 1993), evaporation in January, February and March is overestimated by more than a factor 2. More realistic parameters, as documented for this site, virtually eliminate the error. This study shows the importance of the surface roughness lengths in determining evaporation from wet surfaces. It also illustrates the strength of long observational time series in identifying model deficiencies.     

Environmental and biological controls on elemental (Mg/Ca, Sr/Ca and Mn/Ca) ratios in shells of the king scallop Pecten maximus

The relationship between potential elemental proxies (Mg/Ca, Sr/Ca and Mn/Ca ratios) and environmental factors was investigated for the bivalve Pecten maximus in a detailed field study undertaken in the Menai Strait, Wales, U.K. An age model constructed for each shell by comparison of measured and predicted oxygen-isotope ratios allowed comparison on a calendar time scale of shell elemental data with environmental variables, as well as estimation of shell growth rates. The seasonal variation of shell Mn/Ca ratios followed a similar pattern to one previously described for dissolved Mn²⁺ in the Menai Strait, although further calibration work is needed to validate such a relationship. Shell Sr/Ca ratios unexpectedly were found to co-vary most significantly with calcification temperature, whilst shell Mg/Ca ratios were the next most significant control. The temporal variation in the factors that control shell Sr/Ca ratios strongly suggest the former observation most likely to be the result of a secondary influence on shell Sr/Ca ratios by kinetic effects, the latter driven by seasonal variation in shell growth rate that is in turn influenced in part by seawater temperature. P. maximus shell Mg/Ca ratio to calcification temperature relationships exhibit an inverse correlation during autumn to early spring (October to March-April) and a positive correlation from late spring through summer (May-June to September). No clear explanation is evident for the former trend, but the similarity of the records from the three shells analysed indicate that it is a real signal and not a spurious observation. These observations confirm that application of the Mg/Ca proxy in P. maximus shells remains problematic, even for seasonal or absolute temperature reconstructions. For the range of calcification temperatures of 5-19 °C, our shell Mg/Ca ratios in P. maximus are approximately one-fourth those in inorganic calcite, half those in the bivalve Pinna nobilis, twice those in the bivalve Mytilus trossulus, and four to five times higher than Mg/Ca ratios in planktonic and benthonic foraminifera. Our findings further support observations that Mg/Ca ratios in bivalve shell calcite are an unreliable temperature proxy, as well as substantial taxon- and species-specific variation in Mg incorporation into bivalves and other calcifying organisms, with profound implications for the application of this geochemical proxy to the bivalve fossil record.     

Strong Magnetic Levels in Lake Chapala Sediments (western Mexico): Their Mineralogy and Stratigraphic Significance

Lake Chapala, located 120 km northeast of Colima Volcano, lies at the north and northeast of the Citala rift in the Trans-Mexican Volcanic Belt. It belongs to the watershed of the Lerma River, which originates from the Mexico City area, 450 km to the east. Sediment cores, 0.5–2 m long, were collected from the lake. Magnetic susceptibility of the lake sediments generally ranges from 400 to 800 × 10⁻⁹ m³ kg⁻¹; but in some layers it exceeds 1000 or 1500 × 10⁻⁹ m³ kg⁻¹. The magnetic susceptibility vertical profiles display a thin peak (2–5 cm) or a double peak with magnetic susceptibility > 2000 × 10⁻⁹ m³ kg⁻¹. Scanning electron microscope analysis shows that the main mineral responsible for the magnetic susceptibility is titano-magnetite, which is relatively abundant in the magnetic layers. In most of the cores, the layer with magnetic susceptibility > 2000 × 10⁻⁹ m³ kg⁻¹ is coarser grained than the directly under and overlying sediments, which are composed of fine silt. But not all coarser levels are magnetic nor are all magnetic levels bound to coarser layers. The high titano-magnetic concentration probably originates from volcanic ash falls. Sedimentation rates, evaluated for several sites by means of the ²¹⁰Pb and ¹³⁷Cs methods, allow a date approximation (between 1535 and 1670) for the strong magnetic layer deposition. Since important eruptions of Colima Volcano, with ash fall, were reported from 1576 to 1623, the main peak of magnetic susceptibility in the sediment profiles is assumed to be related to these events. The main magnetic layer extends to greater depth in the profiles of the eastern part of the lake plain than at the west. It can, therefore, be suggested that a relative subsidence occurs in the eastern part of the lake; the axis of the eastern plain overlaps the area displaying the maximum subsidence rate and seems therefore to correspond to the prolongation of Citala rift.