Chalcophile element geochemistry of the Boggy Plain zoned pluton, southeastern Australia: a S-saturated barren compositionally diverse magmatic system

The behavior of the platinum group elements (PGE) and Re in felsic magmas is poorly understood due to scarcity of data. We report the concentrations of Ni, Cu, Re, and PGE in the compositionally diverse Boggy Plain zoned pluton (BPZP), which shows a variation of rock type from gabbro through granodiorite and granite to aplite with a SiO₂ range from 52 to 74 wt %. In addition, major silicate and oxide minerals were analyzed for Ni, Cu, and Re, and a systematic sulfide study was carried out to investigate the role of silicate, oxide, and sulfide minerals on chalcophile element geochemistry of the BPZP. Mass balance calculation shows that the whole rock Cu budget hosted by silicate and oxide minerals is <13 wt % and that Cu is dominantly located in sulfide phases, whereas most of the whole rock Ni budget (>70 wt %) is held in major silicate and oxide minerals. Rhenium is dominantly hosted by magnetite and ilmenite. Ovoid-shaped sulfide blebs occur at the boundary between pyroxene phenocrysts and neighboring interstitial phases or within interstitial minerals in the gabbro and the granodiorite. The blebs are composed of pyrrhotite, pyrite, chalcopyrite, and S-bearing Fe-oxide, which contain total trace metals (Co, Ni, Cu, Ag, Pb) up to ~16 wt %. The mineral assemblage, occurrence, shape, and composition of the sulfide blebs are a typical of magmatic sulfide. PGE concentrations in the BPZP vary by more than two orders of magnitude from gabbro (2.7–7.8 ppb Pd, 0.025–0.116 ppb Ir) to aplite (0.05 ppb Pd, 0.001 ppb Ir). Nickel, Cu, Re, and PGE concentrations are positively correlated with MgO in all the rock types although there is a clear discontinuity between the granodiorite and the granite in the trends for Ni, Rh, and Ir when plotted against MgO. Cu/Pd values gradually increase from 6,100 to 52,600 as the MgO content decreases. The sulfide petrology and chalcophile element geochemistry of the BPZP show that sulfide saturation occurred in the late gabbroic stage of magma differentiation. Segregation and distribution of these sulfide blebs controlled Cu and PGE variations within the BPZP rocks although the magma of each rock type may have experienced a different magma evolution history in terms of crustal assimilation and crystal fractionation. The sulfide melt locked in the cumulate rocks must have sequestered a significant portion of the chalcophile elements, which restricted the availability of these metals to magmatic-hydrothermal ore fluids. Therefore, we suggest that the roof rocks that overlay the BPZP were not prospective for magmatic-hydrothermal Cu, Au, or Cu–Au deposits.     

Diatom-inferred late Pleistocene and Holocene palaeolimnological changes in the Ioannina basin, northwest Greece

The character and impact of climate change since the last glacial maximum (LGM) in the eastern Mediterranean region remain poorly understood. Here, two new diatom records from the Ioannina basin in northwest Greece are presented alongside a pre-existing record and used to infer past changes in lake level, a proxy for the balance between precipitation and evaporation. Comparison of the three records indicates that lake-level fluctuations were the dominant driver of diatom assemblage composition change, whereas productivity variations had a secondary role. The reconstruction indicates low lake levels during the LGM. Late glacial lake deepening was underway by 15.0 cal kyr BP, implying that the climate was becoming wetter. During the Younger Dryas stadial, a lake-level decline is recorded, indicating arid climatic conditions. Lake Ioannina deepened rapidly in the early Holocene, but long-term lake-level decline commenced around 7.0 cal kyr BP. The pattern of lake-level change is broadly consistent with an existing lake-level reconstruction at Lake Xinias, central Greece. The timing of the apparent change, however, is different, with delayed early Holocene deepening at Xinias. This offset is attributed to uncertainties in the age models, and the position of Xinias in the rain shadow of the Pindus Mountains.     

Fluvial sedimentation in a salt‐controlled mini‐basin: stratal patterns and facies assemblages,Sivas Basin,Turkey

The Sivas Basin, located on the Central Anatolian Plateau in Turkey, is an elongate Oligo‐Miocene basin that contains numerous salt‐walled mini‐basins. Through field analysis, including stratigraphic section logging, facies analysis and geological mapping, a detailed tectono‐stratigraphic study of the Emirhan mini‐basin and its 2·6 km thick sediment fill has been undertaken. Three main palaeoenvironments are recognized – playa‐lake, braided stream and lacustrine – each corresponds to a relatively long‐lived depositional episode within a system that was dominated overall by the development of a distributive fluvial system. At local scale, this affects the geometry of the succession and influences facies distributions within preserved sequences. Sequences affected by wedge geometries are characterized by localized channelized sandstone bodies in the area of maximum subsidence and these pass laterally to floodplain mudstone towards the diaper; several internal unconformities are recognized. By contrast, sequences affected by hook geometries display narrow and steep drape‐fold geometries with no evidence of lateral facies change and apparent conformity in the preserved succession. The sediment fill of the Emirhan mini‐basin records the remobilization of diapir‐derived detritus and the presence of evaporitic bodies interbedded within the mini‐basin, implying the growth of salt walls expressed at the surface as palaeo‐topographic highs. The mini‐basin also records the signature of a regional change in stratigraphic assemblage, passing from playa‐lake facies to large‐scale highly amalgamated fluvial facies that represent progradation of the fluvial system. The initiation and evolution of this mini‐basin involves a variety of local and regional controls. Local factors include: (i) salt withdrawal, which influenced the rate and style of subsidence and consequently temporal and spatial variation in the stratigraphic assemblage and the stratal response related to halokinesis; and (ii) salt inflation, which influenced the topographic expression of the diapirs and consequently the occurrence of diapir‐derived detritus intercalated within the otherwise clastic‐dominated succession.     

Geochronological (U–Pb,U–Th–total Pb,Sm–Nd) and geochemical (REE, <Superscript>87</Superscript>Sr/<Superscript>86</Superscript>Sr, δ<Superscript>18</Superscript>O, δ<Superscript>13</Superscript>C) tracing of intraplate tectonism and associated fluid flow in the Warburton Basin,Australia

The Warburton Basin of central Australia has experienced a complex tectonic and fluid-flow history, resulting in the formation of various authigenic minerals. Geochemical and geochronological analyses were undertaken on vein carbonates from core samples of clastic sediments. Results were then integrated with zircon U–Pb dating and uraninite U–Th–total Pb dating from the underlying granite. Stable and radiogenic isotopes (δ¹⁸O, Sr and εNd), as well as trace element data of carbonate veins indicate that >200 °C basinal fluids of evolved meteoric origin circulated through the Warburton Basin. Almost coincidental ages of these carbonates (Sm–Nd; 432 ± 12 Ma) with primary zircon (421 ± 3.8 Ma) and uraninite (407 ± 16 Ma) ages from the granitic intrusion point towards a substantial period of active tectonism and an elevated thermal regime during the mid Silurian. We hypothesise that such a thermal regime may have resulted from extensional tectonism and concomitant magmatic activity following regional orogenesis. This study shows that the combined application of geochemical and geochronological analyses of both primary and secondary species may constrain the timing of tectonomagmatic events and associated fluid flow in intraplate sedimentary basins. Furthermore, this work suggests that the Sm–Nd-isotopic system is surprisingly robust and can record geologically meaningful age data from hydrothermal mineral species.     

Comparing mixing-length models of the diabatic wind profile over homogeneous terrain

Models of the diabatic wind profile over homogeneous terrain for the entire atmospheric boundary layer are developed using mixing-length theory and are compared to wind speed observations up to 300 m at the National Test Station for Wind Turbines at Høvsøre, Denmark. The measurements are performed within a wide range of atmospheric stability conditions, which allows a comparison of the models with the average wind profile computed in seven stability classes, showing a better agreement than compared to the traditional surface-layer wind profile. The wind profile is measured by combining cup anemometer and lidar observations, showing good agreement at the overlapping heights. The height of the boundary layer, a parameter required for the wind profile models, is estimated under neutral and stable conditions using surface-layer turbulence measurements, and under unstable conditions based on the aerosol backscatter profile from ceilometer observations.     

Current and future atmospheric circulation at 500 hPa over Greenland simulated by the CMIP3 and CMIP5 global models

The Greenland ice sheet is projected to be strongly affected by global warming. These projections are either issued from downscaling methods (such as Regional Climate Models) or they come directly from General Circulation Models (GCMs). In this context, it is necessary to evaluate the accuracy of the daily atmospheric circulation simulated by the GCMs, since it is used as forcing for downscaling methods. Thus, we use an automatic circulation type classification based on two indices (Euclidean distance and Spearman rank correlation using the daily 500 hPa geopotential height) to evaluate the ability of the GCMs from both CMIP3 and CMIP5 databases to simulate the main circulation types over Greenland during summer. For each circulation type, the GCMs are compared to three reanalysis datasets on the basis of their frequency and persistence differences. For the current climate (1961–1990), we show that most of the GCMs do not reproduce the expected frequency and the persistence of the circulation types and that they simulate poorly the observed daily variability of the general circulation. Only a few GCMs can be used as reliable forcings for downscaling methods over Greenland. Finally, when applying the same approach to the future projections of the GCMs, no significant change in the atmospheric circulation over Greenland is detected, besides a generalised increase of the geopotential height due to a uniform warming of the atmosphere.     

Evolution of a post-batholith dike swarm in central coastal Queensland, Australia: arc-front to backarc?

A swarm of felsic and mafic dikes cuts a Late Carboniferous–Permian batholith called the Urannah Suite in central coastal Queensland. Late Permian–Triassic westward thrusting (Hunter–Bowen Orogeny) exposed this mid-crustal Suite and the crosscutting dikes, thus enabling examination of dikes that range in age from syn- to post-batholithic. Although both mafic and felsic dikes have the same dominant northerly strike, field, geochronologic and geochemical examination reveal that the swarm is composite; felsic dikes are older (285 Ma) and geochemically and isotopically distinct from mafic dikes (273–229 Ma). Dike compositions are compared to those of the host plutonic rocks, and to volcanic rocks the same age as the dikes. Whereas the felsic (older) dikes are compositionally similar to their host granites (initial ⁸⁷Sr/⁸⁶Sr>0.7045), the mafic (younger) dikes are isotopically (Sr, Nd, Pb) less radiogenic. Moreover, several different types of mafic dikes are evident based on geochemistry; most of these have mixed characteristics in terms of tectonic classification. All but two dikes of basalt and basaltic andesite composition classify as ‘with-in plate' on Ti–Zr–Y tectonic classification plots. Many of the basalts have high TiO₂ contents (1.5–3.0 wt.%). Most of these have REE and spider diagram patterns like calc-alkaline tholeiites, the exceptions being a few alkali basalts recognized by their alkali content, and high Ti, Ce, Nb and Zr contents. When put into the context of all plutonic rocks in the area (late Paleozoic and Mesozoic), these dikes record a transition at 280 Ma, after which time, all magmatism in the region is less isotopically evolved (initial ⁸⁷Sr/⁸⁶Sr=0.7033–0.7044). A model of slab retreat and hinge movement to the east in the latest Permian explains the change of geochemical signature from arc-front to backarc at about 280 Ma.     

Regional Be production rate calibration for the past 12 ka deduced from the radiocarbon-dated Grøtlandsura and Russenes rock avalanches at 69° N, Norway

Two rock avalanches in Troms County – the Grøtlandsura and Russenes – were selected as CRONUS-EU natural cosmogenic ¹⁰Be production-rate calibration sites because they (a) preserve large boulders that have been continuously exposed to cosmic irradiation since their emplacement; (b) contain boulders with abundant quartz phenocrysts and veins with low concentrations of naturally-occurring ⁹Be (typically < 1.5 ppb); and (c) have reliable minimum radiocarbon ages of 11,424 ± 108 cal yr BP and 10,942 ± 77 cal yr BP (1σ), respectively. Quartz samples (n = 6) from these two sites contained between 4.28 × 10⁴ and 5.06 × 10⁴ at ¹⁰Be/g using the 1.387 Myr ¹⁰Be half-life. Determination of these concentrations accounts for topographic and self-shielding, and effects on nuclide production due to isostatic rebound are shown to be negligible. Persistent, constant snow and moss cover cannot be proven, but if taken into consideration they may have reduced ¹⁰Be concentrations by 10%. Using the ¹⁰Be half-life of 1.387 Myr and the Stone scaling scheme, and accounting for snow- and moss-cover, we calculate an error-weighted mean total ¹⁰Be production rate of 4.12 ± 0.19 at/g/yr (1σ). A corresponding error-weighted mean spallogenic ¹⁰Be production rate is 3.96 ± 0.16 at/g/yr (1σ), respectively. These are in agreement within uncertainty with other ¹⁰Be production rates in the literature, but are significantly, statistically lower than the global average ¹⁰Be production rate. This research indicates, like other recent studies, that the production of cosmogenic ¹⁰Be in quartz is lower than previously established by other production-rate calibration projects. Similarly, our findings indicate that regional cosmogenic production rates should be used for determining exposure ages of landforms in order to increase the accuracy of those ages. As such, using the total ¹⁰Be production rate from our study, we determine an error-weighted mean surface-exposure age of a third rock avalanche in Troms County (the Hølen avalanche) to be 7.5 ± 0.3 kyr (1σ). This age suggests that the rock avalanche occurred shortly after the 8.2 kyr cooling event, just as the radiocarbon ages of the Grøtlandsura and Russenes avalanches confirm field evidence that those rock-slope failures occurred shortly after deglaciation.     

Mechanisms controlling silicon isotope distribution in the Eastern Equatorial Pacific

The distribution of silicon isotopes along a meridional transect at 140°W longitude in the Eastern Equatorial Pacific was used to test the hypothesis that δ³⁰Si of silicic acid in surface waters should correlate with net silica production rates (gross silica production minus silica dissolution) rather than rates of gross silica production due to the opposing Si isotope fractionations associated with silica production and silica dissolution. Variations in δ³⁰Si appeared significantly correlated with net silica production rates in equatorial surface waters and not with gross production rates. Around the Equator, values of δ³⁰Si as low as deep water values occurred in the upper mesopelagic in a zone of net silica dissolution and high detrital biogenic silica content, where the release of low δ³⁰Si silicic acid from opal dissolution would be expected to decrease δ³⁰Si. The δ³⁰Si of the deep water at 140°W appears constant for depths >2000 m and is similar to the deep water at 110°W. This study brings to light the importance of considering Si fractionation during diatom silica dissolution, the biological fractionation during silica production and physical factors such as currents and mixing with adjacent water masses when interpreting silicon isotope distributions.