Turbulent flux observations and modelling over a shallow lake and a wet grassland in the Nam Co basin,Tibetan Plateau

The Tibetan Plateau plays an important role in the global water cycle and is strongly influenced by climate change. While energy and matter fluxes have been more intensely studied over land surfaces, a large proportion of lakes have either been neglected or parameterised with simple bulk approaches. Therefore, turbulent fluxes were measured over wet grassland and a shallow lake with a single eddy-covariance complex at the shoreline in the Nam Co basin in summer 2009. Footprint analysis was used to split observations according to the underlying surface, and two sophisticated surface models were utilised to derive gap-free time series. Results were then compared with observations and simulations from a nearby eddy-covariance station over dry grassland, yielding pronounced differences. Observations and footprint integrated simulations compared well, even for situations with flux contributions including grassland and lake. The accessibility problem for EC measurements on lakes can be overcome by combining standard meteorological measurements at the shoreline with model simulations, only requiring representative estimates of lake surface temperature.     

Seasonal prediction of global sea level anomalies using an ocean–atmosphere dynamical model

Advanced warning of extreme sea level events is an invaluable tool for coastal communities, allowing the implementation of management policies and strategies to minimise loss of life and infrastructure damage. This study is an initial attempt to apply a dynamical coupled ocean–atmosphere model to the prediction of seasonal sea level anomalies (SLA) globally for up to 7 months in advance. We assess the ability of the Australian Bureau of Meteorology’s operational seasonal dynamical forecast system, the Predictive Ocean Atmosphere Model for Australia (POAMA), to predict seasonal SLA, using gridded satellite altimeter observation-based analyses over the period 1993–2010 and model reanalysis over 1981–2010. Hindcasts from POAMA are based on a 33-member ensemble of seasonal forecasts that are initialised once per month for the period 1981–2010. Our results show POAMA demonstrates high skill in the equatorial Pacific basin and consistently exhibits more skill globally than a forecast based on persistence. Model predictability estimates indicate there is scope for improvement in the higher latitudes and in the Atlantic and Southern Oceans. Most characteristics of the asymmetric SLA fields generated by El-Nino/La Nina events are well represented by POAMA, although the forecast amplitude weakens with increasing lead-time.     

Responses of microorganisms and enzymes to soil contamination with metazachlor

A greenhouse pot experiment was conducted at the University of Warmia and Mazury in Olsztyn, with the aim of describing the influence of metazachlor on counts and biodiversity of soil microorganisms, soil enzymatic activity, physicochemical properties of soil and yield of spring oilseed rape. The first experimental factor was soil contamination with increasing rates of metazachlor: 0 (soil without the herbicide), 0.333 (recommended by the manufacturer), 6.666, 13.332, 26.665, 53.328, 106.656 and 213.312 mg kg^?1 dm of soil, while the second factor comprised two dates of determination: 30 and 60 days after starting the experiment. The tested herbicide had an adverse effect on reproduction of all analyzed microorganisms: oligotrophic bacteria and their endospore forms, Azotobacter spp. bacteria, organotrophic bacteria, actinomycetes and fungi. The values of the colony development and eco-physiological diversity indices decreased under the influence of excessive doses of the herbicide. Also, significant fluctuation in the enzymatic activity of soil was observed in response to the higher doses of metazachlor and depending on the date of determinations. The herbicide had an invariably negative influence on the activity of soil enzymes, causing the inhibition of dehydrogenases, catalase, urease, acid phosphatase, alkaline phosphatase, arylsulfatase and β-glucosidase. The physicochemical status of soil depended significantly on the degree of soil contamination with the herbicide, same as the yield of spring oilseed rape.     

Streamflow response to regional climate model output in the mountainous watershed: a case study from the Swiss Alps

Regional climate model (RCM) outputs are often used in hydrological modeling, in particular for streamflow forecasting. The heterogeneity of the meteorological variables such as precipitation, temperature, wind speed and solar radiation often limits the ability of the hydrological model performance. This paper assessed the sensitivity of RCM outputs from the PRUDENCE project and their performance in reproducing the streamflow. The soil and water assessment tool was used to simulate the streamflow of the Rhone River watershed located in the southwestern part of Switzerland, with the climate variables obtained from four RCMs. We analyzed the difference in magnitude of precipitation, maximum and minimum air temperature, and wind speed with respect to the observed values from the meteorological stations. In addition, we also focused on the impact of the grid resolution on model performance, by analyzing grids with resolutions of 50 × 50 and 25 × 25 km². The variability of the meteorological inputs from various RCMs is quite severe in the studied watershed. Among the four different RCMs, the Danish Meteorological Institute provided the best performance when simulating runoff. We found that temperature lapse rate is significantly important in the mountainous snow and glacier dominated watershed as compared to other variables like precipitation, and wind speed for hydrological performance. Therefore, emphasis should be given to minimum and maximum temperature in the bias correction studies for downscaling climatic data for impact modeling in the mountainous snow and glacier dominated complex watersheds.     

Behaviour of radiogenic Pb in zircon during ultrahigh-temperature metamorphism: an ion imaging and ion tomography case study from the Kerala Khondalite Belt,southern India

Zircon crystals from a locally charnockitized Paleoproterozoic high-K metagranite from the Kerala Khondalite Belt (KKB) of southern India have been investigated by high-spatial resolution secondary ion mass spectrometry analysis of U–Th–Pb and rare earth elements (REE), together with scanning ion imaging and scanning ion tomography (depth-profiled ion imaging). The spot analyses constrain the magmatic crystallization age of the metagranite to ca. 1,850 Ma, with ultrahigh-temperature (UHT) metamorphism occurring at ca. 570 Ma and superimposed charnockite formation at ca. 520–510 Ma, while the ion imaging reveals a patchy distribution of radiogenic Pb throughout the zircon cores. Middle- to heavy-REE depletion in ca. 570 Ma zircon rims suggests that these grew in equilibrium with garnet and therefore date the UHT metamorphism in the KKB. The maximum apparent ²⁰⁷Pb/²⁰⁶Pb age obtained from the unsupported radiogenic Pb concentrations is also consistent with formation of the Pb patches during this event. The superimposed charnockitization event appears to have caused additional Pb-loss in the cores and recrystallization of the rims. The results of depth-profiling of the scanning ion tomography image stack show that the Pb-rich domains range in size from <5 nm to several 10 nm (diameter if assumed to be spherical). The occurrence of such patchy Pb has previously been documented only from UHT metamorphic zircon, where it likely results from annealing of radiation-damaged zircon. The formation of a discrete, heterogeneously distributed and subsequently immobile Pb phase effectively arrests the normal Pb-loss process seen at lower grades of metamorphism.     

Petrogenesis of a voluminous Quaternary adakitic volcano: the case of Baru volcano

The origin of adakite magmas remains controversial because initially the term adakite had petrogenetic significance implying an origin by direct melting of the eclogitized subducting oceanic crust. Many models have been produced for their origin, and until now there has not been a straightforward method to discriminate between these models in a given adakite suite. Here, we use detailed chronological and geochemical studies of selected adakitic edifices that allows for the determination of the magmatic output rate parameter (Qe), which has been correlated with the rates of magma generation deep within subduction zones. By providing temporal and eruption rate estimates, we provide constraints on the possible petrogenetic processes involved in the generation of adakite-like signatures. Adakite magmas derived from the melting of the subducting slab should be volumetrically insignificant when compared to the adakite-like magmas produced by typical arc magma generation processes. In this study, we use this observation and the extraordinary stratigraphic exposure from Miocene to present in an adakitic volcano in Panama and to study the temporal and chemical variation in erupted magmas to estimate rates of magma generation. Detailed chemical and geochronological analyses of Baru volcano indicate that the volcanic edifice was constructed in its entirety during the Quaternary and magmas display adakite-like features such as steep rare earth elements patterns, pronounced depletions in the heavy rare earth elements, low Y, high Sr, and high Sr/Y. The magmatic output rates (Qe) that we have calculated show that compared to other typical adakitic volcanoes, most of the volcanic edifice of Baru volcano was constructed extremely rapidly (<~213 k.a.) and in time frames that are similar to typical arc volcanoes. The observed chemical and mineralogical variation, coupled with the high magma production rates, indicate that Baru volcano is more representative of a typical arc volcano than a small-volume melt of the subducting oceanic crust. The technique we outline may have broader application in determining the petrogenetic conditions of other adakite suites.     

Origin and significance of poikilitic and mosaic peridotite xenoliths in the western Pannonian Basin: geochemical and petrological evidences

Peridotite xenoliths erupted by alkali basaltic volcanoes in the western Pannonian Basin can be divided into two fundamentally contrasting groups. Geochemical characteristics of the abundant protogranular, porphyroclastic and equigranular nodules suggest that these samples originate from an old consolidated and moderately depleted lithospheric mantle domain. In contrast, the geochemical features of the worldwide rare, but in the Pannonian Basin relatively abundant, poikilitic xenoliths attest to a more complex evolution. It has been argued that the origin of the peculiar texture and chemistry may be intimately linked to melt/rock reactions at successively decreasing liquid volumes in a porous melt flow system. The most likely site where such reactions can take place is the asthenosphere–lithosphere boundary. In this context, poikilitic xenoliths may provide petrological and geochemical evidence for reactions between magmatic liquids issued from the uprising asthenosphere and the solid mantle rocks of the lithosphere. These reactions are important agents of the thermal erosion of the lithosphere; thus, they could have considerably contributed to the thinning of the lithosphere in the Pannonian region. We suggest that in the Pannonian Basin, there could be a strong relation between the unusual abundance of poikilitic mantle xenoliths and the strongly eroded lithosphere.     

Overstepping the garnet isograd: a comparison of QuiG barometry and thermodynamic modeling

The consequences of overstepping the garnet isograd reaction have been investigated by comparing the composition of garnet formed at overstepped P–T conditions (the overstep or “OS” model) with the P–T conditions that would be inferred by assuming garnet nucleated in equilibrium with the matrix assemblage at the isograd (the equilibrium or “EQ” model). The garnet nucleus composition formed at overstepped conditions is calculated as the composition that produces the maximum decrease in Gibbs free energy from the equilibrated, garnet-absent, matrix assemblage for the bulk composition under study. Isopleths were then calculated for this garnet nucleus composition assuming equilibrium with the matrix assemblage (the EQ model). Comparison of the actual P–T conditions of nucleation (the OS model) with those inferred from the EQ model reveals considerable discrepancy between the two. In general, the inferred garnet nucleation P–T conditions (the EQ model) are at a lower temperature and higher or lower pressure (depending on the coexisting calcic phase(s)) than the actual (OS model) nucleation conditions. Moreover, the degree of discrepancy increases with the degree of overstepping. Independent estimates of the pressure of nucleation of garnet were made using the Raman shift of quartz inclusions in garnet (quartz-in-garnet or QuiG barometry). To test the validity of this method, an experimental synthesis of garnet containing quartz inclusions was made at 800 °C, 20 kbar, and the measured Raman shift reproduced the synthesis conditions to within 120 bars. Raman band shifts from three natural samples were then used to calculate an isochore along which garnet was presumed to have nucleated. Model calculations were made at several temperatures along this isochore (the OS model), and these P–T conditions were compared to those computed assuming equilibrium nucleation (the EQ model) to estimate the degree of overstepping displayed by these samples. A sample from the garnet isograd in eastern Vermont is consistent with overstepping of around 10 degrees and 0.6 kbar (affinities of around 2 kJ/mole garnet). A sample from the staurolite–kyanite zone in the same terrane requires overstepping of around 50 °C and 2–5 kbar (affinities of around 10–18 kJ/mole garnet). A similar amount of overstepping was inferred for a blueschist sample from Sifnos, Greece. These results indicate that overstepping of garnet nucleation reactions may be common and pronounced in regionally metamorphosed terranes, and that the P–T conditions and paths inferred from garnet zoning studies may be egregiously in error.     

Thermal and rheological controls on the formation of mafic enclaves or banded pumice

Magma mixing can occur in a fluid manner to produce banded pumice or in a brittle manner to form enclaves. We propose that the critical control on mixing style is a competition between developing networks of crystals in the intruding magma that impart a strength to the magma and melting and disrupting those networks in the host. X-ray computed tomography analysis demonstrates that banded pumice from the 1915 Mt. Lassen eruption lacks crystal networks. In contrast, rhyodacite hosts with mafic enclaves from Chaos Crags contain well-developed networks of large crystals. We present a one-dimensional conductive cooling model that predicts mixing style, either ductile or brittle, as a function of magma compositions, temperatures, and the size of the intruding dike. Our model relies on three assumptions: (1) Mixing is initiated by the injection of a hot dike into a cooler magma body with a yield strength; (2) when magma crystallinity exceeds a critical value, 13 vol% plagioclase, the magma develops a yield strength; and (3) when total crystallinity exceeds 40 vol%, the magma has a penetrative crystal network and is effectively solid. Importantly, because the two magmas are of different compositions, their crystallinities and viscosities do not have the same variations with temperature. As the intruding magma cools, it crystallizes from the outside in, while simultaneously, host magma temperature near the intruder rises. Mixing of the two magmas begins when the host magma is heated sufficiently to (1) disrupt the crystal network and (2) initiate convection. If the shear stress exerted by the convecting host magma on the dike is greater than the yield strength of the dike margin (and dike crystallinity does not exceed 40 %), then fluid mixing occurs, otherwise enclaves form by brittle deformation of the dike. Application of the model to magma compositions representative of Lassen and Chaos Crags shows that emplacement of dikes <1 m thick should produce enclaves, whereas thicker dikes should generate fluid mixing and form banded pumice within days to weeks of emplacement. Similar relationships apply to other modeled magmatic systems, including Pinatubo, Unzen, and Ksudach/Shtuybel’ volcanoes. For all studied systems, the absolute size of the intruding dike, not just its proportion relative to the host, influences mixing style.     

Primary melt from Sannome-gata volcano,NE Japan arc: constraints on generation conditions of rear-arc magmas

The conditions under which rear-arc magmas are generated were estimated using primary basalts from the Sannome-gata volcano, located in the rear of the NE Japan arc. Scoriae from the volcano occur with abundant crustal and mantle xenoliths, suggesting that the magma ascended rapidly from the upper mantle. The scoriae show significant variations in their whole-rock compositions (7.9–11.1 wt% MgO). High-MgO scoriae (MgO > ~9.5 wt%) have mostly homogeneous ⁸⁷Sr/⁸⁶Sr ratios (0.70318–0.70320), whereas low-MgO scoriae (MgO < ~9 wt%) have higher ⁸⁷Sr/⁸⁶Sr ratios (>0.70327); ratios tend to increase with decreasing MgO content. The high-MgO scoriae are aphyric, containing ~5 vol% olivine microphenocrysts with Mg# [100 × Mg/(Mg + Fe²⁺)] of up to 90. In contrast, the low-MgO scoriae have crustal xenocrysts of plagioclase, alkali feldspar, and quartz, and the mineralogic modes correlate negatively with whole-rock MgO content. On the basis of these observations, it is inferred that the high-MgO scoriae represent primary or near-primary melts, while the low-MgO scoriae underwent considerable interaction with the crust. Using thermodynamic analysis of the observed petrological features of the high-MgO scoriae, the eruption temperature of the magmas was constrained to 1,160–1,220 °C. Given that the source mantle was depleted MORB-source mantle, the primary magma was plausibly generated by ~7 % melting of a garnet-bearing spinel peridotite; taking this into consideration, and considering the constraints of multi-component thermodynamics, we estimated that the primary Sannome-gata magma was generated in the source mantle with 0.5–0.6 wt% H₂O at 1,220–1,230 °C and at ~1.8 GPa, and that the H₂O content of the primary magma was 6–7 wt%. The rear-arc Sannome-gata magma was generated by a lower degree of melting of the mantle at greater depths and lower temperatures than the frontal-arc magma from the Iwate volcano, which was also estimated to be generated by ~15 % melting of the source mantle with 0.6–0.7 wt% H₂O at ~1,250 °C and at ~1.3 GPa.