Meteorological Observations 2002 at the Arctic Station,Qeqertarsuaq (69°15'N), Central West Greenland

An automatic meteorological station has been operating at the Arctic Station (69°15'N, 53°31'W) in West Greenland since 1990. This paper summarises meteorological parameters during 2002, including snow cover, ground temperatures and active layer development, and air temperatures at the Station during the last 12 years are compared to large scale trends during the last century.

A compilation of 93 sedimentation rate determinations based on ²¹⁰Pb dating has been carried out for the North Sea-Baltic Sea transition area from a database containing 165 determinations carried out by Danish institutions. In the depositional parts of the area sedimentation rates generally range 25–6403 g m^?2 y^?1. An extreme rate of 13351 g m^?2 y^?1 is observed on a station in the Skagerrak. Sedimentation rates significantly increase with depth indicating that the Skagerrak and northern parts of the Kattegat as well as the deep basins in the Baltic Sea act as depocentres for fine-grained sediments. Apparently, sedimentation rates have increased in recent years.     

The future of the western Baltic Sea: two possible scenarios

Globally coupled climate models are generally capable of reproducing the observed trends in the globally averaged atmospheric temperature. However, the global models do not perform as well on regional scales. Here, we present results from four 100-year, high-resolution ocean model experiments (resolution less than 1 km) for the western Baltic Sea. The forcing is taken from a regional atmospheric model and a regional ocean model, imbedded into two global greenhouse gas emission scenarios, A1B and B1, for the period of 2000 to 2100 with each two realisations. Two control runs from 1960 to 2000 are used for validation. For both scenarios, the results show a warming with an increase of 0.5–2.5 K at the sea surface and 0.7–2.8 K below 40 m. The simulations further indicate a decrease in salinity by 1.5–2 practical salinity units. The increase in water temperature leads to a prolongation of heat waves based on present-day thresholds. This amounts to a doubling or even tripling of the heat wave duration. The simulations show a decrease in inflow events (barotropic/baroclinic), which will affect the deepwater generation and ventilation of the central Baltic Sea. The high spatial resolution allows us to diagnose the inflow events and the mechanism that will cause future changes. The reduction in barotropic inflow events correlates well with the increase in westerly winds. The changes in the baroclinic inflows can be consistently explained by the reduction of calm wind periods and thus a weakening of the necessary stratification in the western Baltic Sea and the Danish Straits.     

Magma plumbing beneath Anak Krakatau volcano, Indonesia: evidence for multiple magma storage regions

Understanding magma plumbing is essential for predicting the behaviour of explosive volcanoes. We investigate magma plumbing at the highly active Anak Krakatau volcano (Indonesia), situated on the rim of the 1883 Krakatau caldera by employing a suite of thermobarometric models. These include clinopyroxene-melt thermobarometry, plagioclase-melt thermobarometry, clinopyroxene composition barometry and olivine-melt thermometry. Petrological studies have previously identified shallow magma storage in the region of 2–8 km beneath Krakatau, while existing seismic evidence points towards mid- to deep-crustal storage zone(s), at 9 and 22 km, respectively. Our results show that clinopyroxene in Anak Krakatau lavas crystallized at a depth of 7–12 km, while plagioclase records both shallow crustal (3–7 km) and sub-Moho (23–28 km) levels of crystallization. These magma storage regions coincide with well-constrained major lithological boundaries in the crust, implying that magma ascent and storage at Anak Krakatau is strongly controlled by crustal properties. A tandem seismic tomography survey independently identified a separate upper crustal (<7 km) and a lower to mid-crustal magma storage region (>7 km). Both petrological and seismic methods are sensitive in detecting magma bodies in the crust, but suffer from various limitations. Combined geophysical and petrological surveys, in turn, offer increased potential for a comprehensive characterization of magma plumbing at active volcanic complexes.     

Implementation of high-order particle-tracking schemes in a water column model

Stochastic differential equations (SDEs) offer an attractively simple solution to complex transport-controlled problems, and have a wide range of physical, chemical, and biological applications, which are dominated by stochastic processes, such as diffusion. As for deterministic ordinary differential equations (ODEs), various numerical scheme exist for solving SDEs. In this paper various particle-tracking schemes are presented and tested for accuracy and efficiency (time vs. accuracy). To test the schemes, the particle tracking algorithms are implemented into a community wide used 1D water column model. Modelling individual particles allows a straightforward physical interpretation of the involved processes. Further, this approach is strictly mass conserving and does not suffer from the numerical diffusion that plagues grid-based methods. Moreover, the Lagrangian framework allows to assign properties to the individual particles, that can vary spatially and temporally. The movement of the particles is described by a stochastic differential equation, which is consistent with the advection–diffusion equation. Here, the concentration profile is represented by a set of independent moving particles, which are advected according to the velocity field, while the diffusive displacements of the particles are sampled from a random distribution, which is related to the eddy diffusivity field.The paper will show that especially the 2nd order schemes are accurate and highly efficient. At the same level of accuracy, the 2nd order scheme can be significantly faster than the 1st order counterpart. This gain in efficiency can be spent on a higher resolution for more accurate solutions at a lower cost.     

Provenance analysis of the Paleozoic sequences of the northern Gondwana margin in NW Iberia: Passive margin to Variscan collision and orocline development

The Cantabrian Zone of NW Iberia preserves a voluminous, almost continuous, sedimentary sequence that ranges from Neoproterozoic to Early Permian in age. Its tectonic setting is controversial and recent hypotheses include (i) passive margin deposition along the northern margin of Gondwana or (ii) an active continental margin or (iii) a drifting ribbon continent. In this paper we present detrital zircon U–Pb laser ablation age data from 13 samples taken in detrital rocks from the Cantabrian Zone sequence ranging from Early Silurian to Early Permian in depositional age. The obtained results, together with previously published detrital zircon ages from Ediacaran–Ordovician strata, allow a comprehensive analysis of changing provenance through time. Collectively, these data indicate that this portion of Iberia was part of the passive margin of Gondwana at least from Ordovician to Late Devonian times. Zircon populations in all samples show strong similarities with the Sahara Craton and with zircons found in Libya, suggesting that NW Iberia occupied a paleoposition close to those regions of present-day northern Africa during this time interval. Changes in provenance in the Late Devonian are attributed to the onset of the collision between Gondwana and Laurussia.Additionally, the Middle Carboniferous to Permian samples record populations consistent with the recycling of older sedimentary sequences and exhumation of the igneous rocks formed before and during the Variscan orogeny. Late-Devonian to Permian samples yield zircon populations that reflect topographic changes produced during the Variscan orogeny and development of the lithospheric scale oroclinal buckling.     

Mixing,hypersalinity and gradients in Hervey Bay,Australia

Hervey Bay, a large coastal embayment situated off the central eastern coast of Australia, is a shallow tidal area (average depth = 15 m), close to the continental shelf. It shows features of an inverse estuary, due to the high evaporation rate (approx. 2 m/year), low precipitation (less than 1 m/year) and on average almost no freshwater input from rivers that drain into the bay. The hydro- and thermodynamical structures of Hervey Bay and their variability are presented here for the first time, using a combination of four-dimensional modelling and observations from field studies. The numerical studies are performed with the Coupled Hydrodynamical Ecological Model for Regional Shelf Seas (COHERENS). Due to the high tidal range (>3.5 m), the bay is considered as a vertically well-mixed system, and therefore, only horizontal fronts are likely. Recent field measurements, but also the numerical simulations, indicate characteristic features of an inverse/hypersaline estuary with low salinity (35.5 psu) in the open ocean and peak values (>39.0 psu) in the head water of the bay. The model further predicts a nearly persistent mean salinity gradient of 0.5 psu across the bay (with higher salinities close to the shore). The investigation further shows that air temperature, wind direction and tidal regime are mainly responsible for the stability of the inverse circulation and the strength of the salinity gradient across the bay. Due to an ongoing drying trend, the occurrence of severe droughts at the central east coast of Australia and, therefore, a reduction in freshwater supply, the salinity flux out of the bay has increased, and the inverse circulation has also strengthened.     

Response of xenotime to prograde metamorphism

Xenotime is a widespread accessory mineral in lower greenschist to upper amphibolite facies metasedimentary rocks from the Palaeoproterozoic Mount Barren Group, southwestern Australia. Xenotime is closely associated with detrital zircon, commonly forming syntaxial outgrowths, in samples of sandstone, micaceous quartzite, slate, phyllite, garnet-bearing semi-pelites, and in kyanite-, garnet-, and staurolite-bearing mica schists. In situ geochronology of xenotime from lower greenschist sandstones has previously yielded multiple U–Pb ages with peaks at ~2.0, ~1.7, and ~1.65 Ga, interpreted to represent the age of detritus, early diagenesis, and a later thermal event, respectively. New U–Pb dating of xenotime in slate yields a major population at ~1.7 Ga with a minor population at ~1.2 Ga, reflecting diagenetic and metamorphic growth, respectively, whereas xenotime in phyllite forms a minor age population at ~1.7 Ga and a main peak at ~1.2 Ga. Mid-greenschist facies semi-pelitic schists (quartz-muscovite-garnet) contain xenotime that formed before 1.8 Ga and at 1.2 Ga, representing detrital and peak metamorphic ages, respectively. Xenotime in samples of amphibolite facies schist (650°C and ~8 kbars) yields U–Pb ages of ~1.2 Ga, coinciding with the time of peak metamorphism. A single analysis of a xenotime core from an amphibolite facies schist gave an age of ~1.8 Ga, consistent with the presence of detrital xenotime. Our results suggest that detrital xenotime may be preserved under greenschist facies conditions, but is largely replaced during upper amphibolite facies conditions. Detrital xenotime is replaced through dissolution–reprecipitation reactions forming compositionally distinct rims during greenschist and amphibolite facies metamorphism at 1.2 Ga. Diagenetic xenotime is present in lower greenschist facies samples, but was not observed in metasedimentary rocks that had experienced temperatures above mid-greenschist facies metamorphism (450°C). The apparent disappearance of detrital and diagenetic xenotime and appearance of metamorphic xenotime during prograde metamorphism indicates that some of the yttrium, heavy rare earth elements, and phosphorus needed for metamorphic xenotime growth are probably derived from the replacement of detrital and diagenetic xenotime.     

Linking gold mineralization to regional-scale drivers of mineral systems using in situ U–Pb geochronology and pyrite LA-ICP-MS element mapping

Proterozoic orogens commonly host a range of hydrothermal ores that form in diverse tectonic settings at different times. However, the link between mineralization and the regional-scale tectonothermal evolution of orogens is usually not well understood, especially in areas subject to multiple hydrothermal events.Regional-scale drivers for mineral systems vary between the different classes of hydrothermal ore, but all involve an energy source and a fluid pathway to focus mineralizing fluids into the upper crust. The Mount Olympus gold deposit in the Proterozoic Capricorn Orogen of Western Australia, was regarded as an orogenic gold deposit that formed at ca. 1738 Ma during the assembly of Proterozoic Australia. However,the trace element chemistry of the pyrite crystals closely resembles those of the Carlin deposits of Nevada,with rims that display solid solution gold accompanied by elevated As, Cu, Sb, Hg, and Tl, surrounding gold-poor cores. New SHRIMP UeP b dating of xenotime intergrown with auriferous pyrite and ore-stage alteration minerals provided a weighted mean~(207) Pb*/~(206) Pb* date of 1769 ± 5 Ma, interpreted as the age of gold mineralization. This was followed by two discrete episodes of hydrothermal alteration at 1727 ± 7 Ma and 1673 ± 8 Ma. The three ages are linked to multiple reactivation of the crustal-scale Nanjilgardy Fault during repeated episodes of intracratonic reworking. The regional-scale drivers for Carlin-like gold mineralization at Mount Olympus are related to a change in tectonic regime during the final stages of the intracratonic 1820 -1770 Ma Capricorn Orogeny. Our results suggest that substantial sized Carlin-like gold deposits can form in an intracratonic setting during regional-scale crustal reworking.     

Pb/Pb geochronology, petrography and chemistry of Zr-rich accessory minerals (zirconolite, tranquillityite and baddeleyite) in mare basalt 10047

In situ U-Pb isotopic measurements were carried out by ion microprobe on the Zr-rich accessory minerals zirconolite [CaZrTi₂O₇], tranquillityite [Fe₈²⁺(ZrY)₂Ti₃Si₃O₂₄] and baddeleyite [ZrO₂] in low-K, high-Ti mare basalt 10047 collected during the Apollo 11 mission. The analysed minerals are concentrated in pockets of late-stage mesostasis that comprises an intergrowth of silica, barian K-feldspar and Si-Al-K glass, from a phaneritic, subophitic, basalt comprising mainly pyroxene, plagioclase, ilmenite, cristobalite and troilite. Most Zr-rich minerals are unaltered, however, some tranquillityite is replaced by a complex intergrowth of zirconolite, baddeleyite, ilmenite and fayalite, suggesting that the mineral became unstable during crystallization. Several baddeleyite crystals have also undergone alteration to secondary zircon. Zirconolite was analysed in thin section 10047,11 and tranquillityite and baddeleyite in 10047,227, using a ∼6 μm primary ion beam. Both zirconolite and tranquillityite have significant U and low initial Pb contents, and are highly suitable for Pb/Pb dating. Fifteen analyses of zirconolite give a ²⁰⁷Pb/²⁰⁶Pb age of 3708 ± 7 Ma (²⁰⁷Pb/²⁰⁶Pb:²⁰⁴Pb/²⁰⁶Pb isochron; 95% confidence, including renormalisation of ratios) and twenty-five analyses of tranquillityite give 3710 ± 6 Ma. The ²⁰⁷Pb/²⁰⁶Pb dates are consistent with each other and refine results from an earlier study. Baddeleyite data were less precise, mainly due to lower secondary ionisation efficiency. Our results show that zirconolite and tranquillityite can provide precise isotopic dates and, given their presence in other samples, they represent important U-Pb chronometers for refining lunar geology.