The magmatic plumbing system beneath El Hierro (Canary Islands): constraints from phenocrysts and naturally quenched basaltic glasses in submarine rocks

A thermobarometric and petrologic study of basanites erupted from young volcanic cones along the submarine portions of the three El Hierro rift zones (NE-Rift, NW-Rift and S-Ridge) has been performed to reconstruct magma plumbing and storage beneath the island. Mineral-melt thermobarometry applied to naturally quenched glass and clinopyroxene rims yields pressures ranging from 350 to 1070 MPa with about 80% of the calculated pressures being in the range of 600–800 MPa. This corresponds to a depth range of 19–26 km, implying that the main level of final crystal fractionation is within the uppermost mantle. No systematic dependence between sample locality and fractionation pressures could be observed. Olivine and clinopyroxene crystals in the rocks are complexly zoned and have, on an inter-sample as well as on an intra-sample scale, highly variable core and rim compositions. This can best be explained by mixing of multiply saturated (olivine, magnetite, clinopyroxene, ilmenite), moderately evolved magmas with more mafic magmas being either only saturated with olivine + spinel or with olivine + spinel + clinopyroxene. The inter-sample differences indicate derivation from small, isolated magma chambers which have undergone distinct fractionation and mixing histories. This is in contrast to oceanic intraplate volcanoes situated on plumes with high melt supply rates, e.g. Kilauea Volcano (Hawaii), where magma is mainly transported through a central conduit system and stored in a shallow magma chamber prior to injection into the rift zones. The plumbing system beneath El Hierro rather resembles the magma storage systems beneath, e.g. Madeira or La Palma, indicating that small, intermittent magma chambers might be a common feature of oceanic islands fed by plumes with relatively low fluxes, which results in only limited and periodic magma supply.     

Soil properties controlling Zn speciation and fractionation in contaminated soils

We determined the speciation of Zn in 49 field soils differing widely in pH (4.1–7.7) and total Zn content (251–30,090 mg/kg) by using extended X-ray absorption fine structure (EXAFS) spectroscopy. All soils had been contaminated since several decades by inputs of aqueous Zn with runoff-water from galvanized power line towers. Pedogenic Zn species identified by EXAFS spectroscopy included Zn in hydroxy-interlayered minerals (Zn-HIM), Zn-rich phyllosilicates, Zn-layered double hydroxide (Zn-LDH), hydrozincite, and octahedrally and tetrahedrally coordinated sorbed or complexed Zn. Zn-HIM was only observed in (mostly acidic) soils containing less than 2000 mg/kg of Zn, reflecting the high affinity but limited sorption capacity of HIM. Zn-bearing precipitates, such as Zn-LDH and Zn-rich trioctahedral phyllosilicates, became more dominant with increasing pH and increasing total Zn content relative to available adsorption sites. Zn-LDH was the most abundant Zn-precipitate and was detected in soils with pH > 5.2. Zn-rich phyllosilicates were detected even at lower soil pH, but were generally less abundant than Zn-LDH. Hydrozincite was only identified in two calcareous soils with extremely high Zn contents. In addition to Zn-LDH, large amounts of Zn in highly contaminated soils were mainly accumulated as sorbed/complexed Zn in tetrahedral coordination. Soils grouped according to their Zn speciation inferred from EXAFS spectroscopy mainly differed with respect to soil pH and total Zn content. Clear differences were observed with respect to Zn fractionation by sequential extraction: From Zn-HIM containing soils, most of the total Zn was recovered in the exchangeable and the most recalcitrant fractions. In contrast, from soils containing the highest percentage of Zn-precipitates, Zn was mainly extracted in intermediate extraction steps. The results of this study demonstrate that soil pH and Zn contamination level relative to available adsorption sites are the most important factors controlling the formation of pedogenic Zn-species in aerobic soils and, consequently, Zn fractionation by sequential extraction.     

The importance of melt extraction for tracing mantle heterogeneity

Numerous isotope and trace element studies of mantle rocks and oceanic basalts show that the Earth’s mantle is heterogeneous. The isotopic variability in oceanic basalts indicates that most mantle sources consist of complex assemblages of two or more components with isolated long-term chemical evolution, on both global and local scales. The range in isotope and highly incompatible element ratios observed in oceanic basalts is commonly assumed to directly reflect that of their mantle sources. Accordingly, the end-points of isotope arrays are taken to represent the isotopic composition of the different components in the underlying mantle, which is then used to deduce the origin of mantle heterogeneity. Here, a melting model for heterogeneous mantle sources is presented that investigates how and to what extent isotope and trace element signatures are conveyed from source to melt. We model melting of a pyroxenite-bearing peridotite using recent experimental constrains for melting and partitioning of pyroxenite and peridotite. Identification of specific pyroxenite melting signatures allows finger-printing of pyroxenite melts and confirm the importance of lithological heterogeneity in the Earth’s mantle. The model results and the comparison of the calculated and observed trace element-isotope systematics in selected MORB and OIB suites (e.g. from the East Pacific Rise, Iceland, Tristan da Cunha, Gough and St.Helena) further show that factors such as the relative abundance of different source components, their difference in solidus temperature, and especially the extent, style and depth range of melt aggregation fundamentally influence the relationship between key trace element and isotope ratios (e.g. Ba/Th, La/Nb, Sr/Nd, La/Sm, Sm/Yb, ¹⁴³Nd/¹⁴⁴Nd). The reason for this is that any heterogeneity present in the mantle is averaged or, depending on the effectiveness of the melt mixing process, even homogenized during melting and melt extraction. Hence to what degree mantle heterogeneity is reflected in the erupted melts is not only a function of source and melting-induced variability. It also depends on the extent of mixing during melting and melt extraction and thus strongly on the relative incompatibility of the elements considered. The observed trace element variation in erupted melts can be greater or smaller than that of their mantle sources, depending on the incompatibility of the elements investigated. The isotopic variability in erupted melts, on the other hand, is generally smaller than that of their mantle source. Melt mixing during melt extraction consequently has an important influence on the relative extent of variation, and hence the degree of correlation between the isotope and trace element ratios. Overall fewer correlations between trace element and isotope ratios are expected whenever melts are extracted from a restricted depth range, as is the case for ocean island basalts, than for cases where melts are extracted over a larger depth interval (mid ocean ridges and especially ridge centered hotspots like Iceland). While the isotopic composition of the most enriched melts may correspond closely to those of the enriched source component, even the most depleted mid ocean ridge basalts are likely to underestimate the isotopic depletion of the depleted mantle component. These observations imply that using the chemical and isotopic range observed in oceanic basalts as directly representative of that in the corresponding mantle source can be misleading, since this assumption is strictly true only for homogeneous mantle sources. In addition to identifying source or partitioning-related differences in melts from different mantle sources, inferring the true composition, origin, and distribution of heterogeneous components in the Earth’s mantle therefore requires detailed knowledge about the mechanisms of melting and melt mixing during the melt extraction process. Only if these processes and their influence on the isotope-trace element relationship are understood, can the composition and origin of the different source components, and thus mantle heterogeneity, be accurately constrained.     

The Fernbank interglacial site near Ithaca, New York, USA

The Fernbank interglacial site, on the west side of Cayuga Lake, New York, has been recently subjected to more detailed study. To a lengthened mollusc list are added ostracodes, insects, fish, pollen, and plant macrofossils. Of these, plants are well preserved and diverse, whereas other groups are poorly preserved and incomplete. Nevertheless, all support the interglacial assignment (Sangamon), which is further supported by minimum age radiocarbon dates (> 50,000 ¹⁴C yr BP) and a TL date of 81 ± 11 ka. In the plant record near the top of the sequence, abundant tree charcoal indicates forest fires. Like the Toronto interglacial record, the plants show a declining July mean temperature from 24 to 18°C (according to transfer functions) through the sequence, from mixed deciduous forest to boreal forest.     

Sea ice induced changes in ocean circulation during the Eemian

We argue that Arctic sea ice played an important role during early stages of the last glacial inception. Two simulations of the Institut Pierre Simon Laplace coupled model 4 are analyzed, one for the time of maximum high latitude summer insolation during the last interglacial, the Eemian, and a second one for the subsequent summer insolation minimum, at the last glacial inception. During the inception, increased Arctic freshwater export by sea ice shuts down Labrador Sea convection and weakens overturning circulation and oceanic heat transport by 27 and 15%, respectively. A positive feedback of the Atlantic subpolar gyre enhances the initial freshening by sea ice. The reorganization of the subpolar surface circulation, however, makes the Atlantic inflow more saline and thereby maintains deep convection in the Nordic Seas. These results highlight the importance of an accurate representation of dynamic sea ice for the study of past and future climate changes.     

Evaluation of atmospheric thermal radiation algorithms for daylight hours

Existing simple but theoretically based clear-sky models for longwave down-welling radiation (LDR_c) and cloud impact algorithms transforming them to all-sky radiation (LDR) are checked against locally calibrated empirical algorithms. They are evaluated for daylight hours based on measurements in regionally differing climates of Germany. The Prata clear-sky scheme is additionally tested with adjusted coefficients so that LDR_c converges against a realistic emissivity for a completely dry atmosphere. This version is characterised by an improved modelled variance. Compared with locally calibrated schemes, bias and root mean square error (RMSE) of the more theoretical clear-sky schemes do not differ significantly and yield even better results at a mountain site. In contrast, the locally calibrated algorithms yield biases up to 9% and an increase in RMSE between 6% and 67%, if applied for other sites. For daylight hours, the cloud impact on LDR can be calculated via the ratio of observed to clear-sky global irradiation (CMF_sol). With CMF_sol, the Crawford and Duchon scheme reveals the lowest bias and a decrease in RMSE by 22% against the next best performing algorithms. Compared with synoptic cloud observations as input, the bias is reduced by 9 to 28 W m^?2 and the scattering of the residuals decreases by 20% to 30%. Based on published results for also non-European sites, it is inferred that the more theoretically based LDR_c schemes and cloud impact evaluated via CMF_sol are universally applicable and perform at least in the order of magnitude of locally calibrated empirical algorithms.     

A reconstructed dynamic Indian monsoon index extended back to 1880

The authors present a reconstruction of summer (June–July–August) mean dynamic Indian monsoon index (DIMI) back to 1880 based on a large number of historical surface observation data as well as information from the upper air data. The reconstruction shows a satisfying skill in terms of both the value of reduction of error and an evaluation against other independent monsoon indices. The skill of reconstruction increases over time with more predictor data (in particular upper-level data) becoming available. A comparison with the observed all Indian summer monsoon rainfall index (AIRI) shows a high consistence in both inter-decadal and inter-annual variability. The reconstruction shows stronger than normal monsoon during the 1880s, 1915–1925 (around 1920) and 1930–1945 (around 1940) as the AIRI. The El Nino/Southern Oscillation (ENSO)—monsoon relationship is reasonably captured in the reconstruction. Powers concentrating within quasi-biennial band stand out in the reconstruction as well as in the AIRI. A comparison of the reconstruction against an atmospheric general circulation model simulation with specified SST and external forcing agents spanning 1901–1999 indicates a slightly higher reproducibility of monsoon circulation than monsoon rainfall in terms of interannual variability. The relationship between the Asian continent warming and the ENSO–monsoon connection is also discussed by using the new dynamic index.     

Effects of land cover change on the tropical circulation in a GCM

Multivariate statistics are used to investigate sensitivity of the tropical atmospheric circulation to scenario-based global land cover change (LCC), with the largest changes occurring in the tropics. Three simulations performed with the fully coupled Parallel Climate Model (PCM) are compared: (1) a present day control run; (2) a simulation with present day land cover and Intergovernmental Panel on Climate Change (IPCC) Special Report on Emission Scenarios (SRES) A2 greenhouse gas (GHG) projections; and (3) a simulation with SRES A2 land cover and GHG projections. Dimensionality of PCM data is reduced by projection onto a priori specified eigenvectors, consisting of Rossby and Kelvin waves produced by a linearized, reduced gravity model of the tropical circulation. A Hotelling T ² test is performed on projection amplitudes. Effects of LCC evaluated by this method are limited to diabatic heating. A statistically significant and recurrent signal is detected for 33% of all tests performed for various combinations of parameters. Taking into account uncertainties and limitations of the present methodology, this signal can be interpreted as a Rossby wave response to prescribed LCC. The Rossby waves are shallow, large-scale motions, trapped at the equator and most pronounced in boreal summer. Differences in mass and flow fields indicate a shift of the tropical Walker circulation patterns with an anomalous subsidence over tropical South America.     

Thermal bioclimate in Strasbourg - the 2003 heat wave

This case study highlights the implications of the 2003 heat wave for the city of Strasbourg, France. The urban centers of France and other European countries were particularly affected by the heat wave. In some urban areas, the mortality rate was 60% above the expected value (Institute de Veille Sanitaire, 2003). The 2003 heat wave demonstrated once again that populations in urban centers are much more affected by extreme meteorological events than people living in rural areas. The aim of this analysis is to explore differences in thermal comfort conditions of (a) the city center of Strasbourg, and (b) its hinterland. The differences in thermal conditions existing between rural and urban areas are quantified by using a bio-climatological index termed physiologically equivalent temperature (PET). This index is based on the human energy balance and builds a relevant index for the quantification of the thermal environment of humans. We calculate the PET for the years 2003 and 2004 to highlight the temporal changes in the severity of climate extremes. The spatial scope of this study is improved compared to previous works in the field through the inclusion of PET calculations for five different sites on a central place in Strasbourg (Place Kléber). The calculations are characterized by different sky view factors and are compared to the reference site, which is located in a rural area. In the rural hinterland (Entzheim), the analysis of PET indicates a strong cold thermal stress during the winter months but no significant stress in summer. In 2003, summer temperatures were sensed as warmer compared to other years, but did not reach the extreme temperatures that may cause severe heat stress. For both the rural and the urban study sites PET was higher in the summer of 2003 than in 2004, which reflects the inferior thermal conditions in the urban area during the heat wave in 2003. For the entire study period, urban and rural day-time PET reached similar maximal values. Strong differences in PET, however, were observed between the rural and urban areas at night-time. The study of PET for several study sites on a central place in the city (Place Kléber) of Strasbourg for the years 2003 and 2004 showed that the sites with a higher sky view factor present higher values than sites with a lower sky view factor. The comparison of these PET values (Place Kléber) to the results for the rural area showed that during the day and the night the rural city of Entzheim has the lowest PET. During the day, the site at Place Kléber, which is located under a tree, has the lowest PET. The comparison of PET for the years 2003 and 2004 shows that PET in 2003 was about 5 to 7 K higher.     

The study of artificial intelligence for predicting land use changes in an arid ecosystem

Journal of Geographical Sciences - During the 21st century, artificial intelligence methods have been broadly applied in geosciences to simulate complex dynamic ecosystems, but the use of...